CS Update

ClusterShapeAnalysis/ClusterHists.h

@@ -1,5 +1,6 @@
 #pragma once
 
+#include <TH3.h>
 #include <TH2.h>
 #include <TH1.h>
 
@@ -8,6 +9,7 @@
 namespace EVENT
 {
   class TrackerHit;
+  class SimTrackerHit;
 }
 
 //! Histograms for cluster analysis
@@ -22,6 +24,7 @@ class ClusterHists
 
   // Fill histograms with a single track hit
   void fill(const EVENT::TrackerHit* trkhit);
+  void fill(const EVENT::SimTrackerHit* simtrkhit);
 
   TH1* h_cluster_edep_BX;
 
@@ -115,4 +118,88 @@ class ClusterHists
   TH2* h_z_r_vx;
   TH2* h_x_y_vx;
 
+  //time and edep cut for all layers -- juliet
+  TH1* h_cluster_edep_Tcut;
+  TH1* h_cluster_edep_Tcut_layer0;
+  TH1* h_cluster_edep_Tcut_layer1;
+  TH1* h_cluster_edep_Tcut_layer2;   
+  TH1* h_cluster_edep_Tcut_layer3;
+  TH1* h_cluster_edep_Tcut_layer4;
+  TH1* h_cluster_edep_Tcut_layer5;
+  TH1* h_cluster_edep_Tcut_layer6;
+  TH1* h_cluster_edep_Tcut_layer7;
+  TH1* h_cluster_edep_Tcut_layer8;
+
+  TH1* h_trackerhit_time_Tcut;
+  TH1* h_trackerhit_time_Tcut_layer0;
+  TH1* h_trackerhit_time_Tcut_layer1;
+  TH1* h_trackerhit_time_Tcut_layer2;
+  TH1* h_trackerhit_time_Tcut_layer3;
+  TH1* h_trackerhit_time_Tcut_layer4;
+  TH1* h_trackerhit_time_Tcut_layer5;
+  TH1* h_trackerhit_time_Tcut_layer6;
+  TH1* h_trackerhit_time_Tcut_layer7;
+
+  //time and edep for 1st layer -- juliet
+  TH1* h_cluster_edep_layer0;
+  TH1* h_cluster_edep_layer1;
+  TH1* h_cluster_edep_layer2;   
+  TH1* h_cluster_edep_layer3;
+  TH1* h_cluster_edep_layer4;
+  TH1* h_cluster_edep_layer5;
+  TH1* h_cluster_edep_layer6;
+  TH1* h_cluster_edep_layer7;
+  TH1* h_cluster_edep_layer8;
+
+//hit energy depo per layer
+  TH1* h_hit_edep_layer0;
+  TH1* h_hit_edep_layer1;
+  TH1* h_hit_edep_layer2;   
+  TH1* h_hit_edep_layer3;
+  TH1* h_hit_edep_layer4;
+  TH1* h_hit_edep_layer5;
+  TH1* h_hit_edep_layer6;
+  TH1* h_hit_edep_layer7;
+  TH1* h_hit_edep_layer8;
+//time of arriver per layer
+  TH1* h_trackerhit_time;
+  TH1* h_trackerhit_time_layer0;
+  TH1* h_trackerhit_time_layer1;
+  TH1* h_trackerhit_time_layer2;
+  TH1* h_trackerhit_time_layer3;
+  TH1* h_trackerhit_time_layer4;
+  TH1* h_trackerhit_time_layer5;
+  TH1* h_trackerhit_time_layer6;
+  TH1* h_trackerhit_time_layer7;
+  TH1* h_trackerhit_time_layer8;
+  //hits per cluster per layer: 
+  TH1* h_thclen;
+  TH1* h_thclen_layer0;
+  TH1* h_thclen_layer1;
+  TH1* h_thclen_layer2;
+  TH1* h_thclen_layer3;
+  TH1* h_thclen_layer4;
+  TH1* h_thclen_layer5;
+  TH1* h_thclen_layer6;
+  TH1* h_thclen_layer7;
+  TH1* h_thclen_layer8;
+  //hits per cluster per layer cut: 
+  TH1* h_thclen_cut;
+  TH1* h_thclen_layer0_cut;
+  TH1* h_thclen_layer1_cut;
+  TH1* h_thclen_layer2_cut;
+  TH1* h_thclen_layer3_cut;
+  TH1* h_thclen_layer4_cut;
+  TH1* h_thclen_layer5_cut;
+  TH1* h_thclen_layer6_cut;
+  TH1* h_thclen_layer7_cut;
+  TH1* h_thclen_layer8_cut;
+  //2D
+  TH2* h_toa_edepCluster;
+  TH2* h_cluster_edep_thlen;
+  TH2* h_cluster_edep_thlen_cut;
+  //3D Histos: 
+  TH3* h_3DPosition_digi;
+  TH3* h_3DPosition_cdigi;  
+
 };