Electron data set

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

ecal/scripts/makeplot.C

@@ -22,8 +22,8 @@ int makeplot(const char* input_fname = "sim_output/sim_emcal_electrons_output.ro
   //if(!fp)
   //      fprintf(stderr,"error: can't open %s\n", fout);
 
-  Double_t e_start = 1.0;
-  Double_t e_end   = 1.0;
+  //Double_t e_start = 1.0;
+  //Double_t e_end   = 1.0;
 
   // Input ROOT file
   TFile *f = new TFile(input_fname,"READ");
@@ -80,8 +80,10 @@ int makeplot(const char* input_fname = "sim_output/sim_emcal_electrons_output.ro
   // Declare histograms
   TH1D *h1 = new TH1D("h1","Scattering Angle(#theta)",          100,130.0,180.0);
   TH1D *h2 = new TH1D("h2","Pseudo-rapidity(#eta)",             100,-5.0,0.0);
-  TH2D *h3 = new TH2D("h3","Cluster E vs Pseudo-rapidity",      100,e_start-0.5,e_end+0.5,100,-5.0,0.0);
-  TH1D *h4 = new TH1D("h4","Reconstructed energy per event",    100,e_start-0.5,e_end+0.5);
+  TH2D *h3 = new TH2D("h3","Cluster E vs Pseudo-rapidity",      100,-0.5,30.5,100,-5.0,0.0);
+  //TH2D *h3 = new TH2D("h3","Cluster E vs Pseudo-rapidity",      100,e_start-0.5,e_end+0.5,100,-5.0,0.0);
+  TH1D *h4 = new TH1D("h4","Reconstructed energy per event",    100,-0.5,30.5);
+  //TH1D *h4 = new TH1D("h4","Reconstructed energy per event",    100,e_start-0.5,e_end+0.5);
   TH1D *h5 = new TH1D("h5","Number of Clusters per event",      5,-0.5,4.5);
   TH1D *h6 = new TH1D("h6","Scattering Angle(#theta) with CUT", 100,130.0,180.0);
   TH1D *h7 = new TH1D("h7","Pseudo-rapidity(#eta) with CUT",    100,-5.0,0.0);