diff --git a/ecal/scripts/makeplot.C b/ecal/scripts/makeplot.C index c94fe49b8ce91fd639a72cb0d0a1115166829fa4..11a66bfa91e1eeeb47abb9b444ce4987de356fb1 100644 --- a/ecal/scripts/makeplot.C +++ b/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);