Electron data set

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

ecal/scripts/makeplot.C

@@ -72,7 +72,7 @@ int makeplot(double e_start = 1.0, double e_end = 1.0, const char* input_fname =
   TCanvas *c7  = new TCanvas("c7", "c7",  500, 500);
   TCanvas *c8  = new TCanvas("c8", "c8",  500, 500);
   TCanvas *c9  = new TCanvas("c9", "c9",  500, 500);
-
+  cout << "Before Histogram" << endl;
   // 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);
@@ -106,7 +106,7 @@ int makeplot(double e_start = 1.0, double e_end = 1.0, const char* input_fname =
   for (int ievent = 0; ievent < nentries; ievent++)
   {
 	t->GetEntry(ievent);
-
+	cout << "Inside loop event by event" << endl;
 	Int_t ncluster   = EcalClusters_;
 	Int_t nreconhits = RecoEcalHits_;
 
@@ -149,7 +149,7 @@ int makeplot(double e_start = 1.0, double e_end = 1.0, const char* input_fname =
         for(int ireconhit=0; ireconhit < nreconhits; ireconhit++)
 		h9->Fill(rec_x_pos[ireconhit],rec_y_pos[ireconhit], 1.0);
   }
-
+  cout << "Before Drawing" << endl;
   // Drawing and Saving figures
   c1->cd();
   h1->SetLineColor(kBlue);