Skip to content
Snippets Groups Projects
Commit 787ae1ca authored by Jihee Kim's avatar Jihee Kim
Browse files

Histogram setting

parent 5a76c249
No related branches found
No related tags found
No related merge requests found
......@@ -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);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment