Electron data set

• Energy range as an input
• A script to read input ROOT file
• A script for energy resolution

ecal/scripts/makeplot.C

@@ -41,7 +41,7 @@ int makeplot(const char* input_fname = "sim_output/sim_emcal_electrons_output.ro
   t->SetBranchStatus("RecoEcalHits", 1);
   t->SetBranchAddress("RecoEcalHits", &RecoEcalHits_);
 
-  const Int_t kMaxEcalClusters = 1000000;
+  const Int_t kMaxEcalClusters = 100000;
   Double_t cluster_x_pos[kMaxEcalClusters];
   Double_t cluster_y_pos[kMaxEcalClusters];
   Double_t cluster_z_pos[kMaxEcalClusters];
@@ -55,7 +55,7 @@ int makeplot(const char* input_fname = "sim_output/sim_emcal_electrons_output.ro
   t->SetBranchAddress("EcalClusters.position.z",cluster_z_pos);
   t->SetBranchAddress("EcalClusters.energy",cluster_energy);
 
-  const Int_t kMaxRecoEcalHits = 1000000;
+  const Int_t kMaxRecoEcalHits = 100000;
   Double_t rec_x_pos[kMaxRecoEcalHits];
   Double_t rec_y_pos[kMaxRecoEcalHits];
   Double_t rec_energy[kMaxRecoEcalHits];