| 1 |
|
| 2 |
void analyze(const char* infile = "/d101/yetkin/work/zdcScint/r150476.root", const char* outfile = "plotsRun150476.root"){
|
| 3 |
|
| 4 |
const char* trig = "HLT_HIMinBiasHF";
|
| 5 |
double towerCut = 3;
|
| 6 |
double ebCut = 3;
|
| 7 |
|
| 8 |
bool MC = false;
|
| 9 |
|
| 10 |
TFile * inf = new TFile(infile);
|
| 11 |
TFile* outf = new TFile(outfile,"recreate");
|
| 12 |
|
| 13 |
TH1D* ha[100];
|
| 14 |
TH1D* hb[100];
|
| 15 |
TH1D* hc[100];
|
| 16 |
TH1D* he[100];
|
| 17 |
TH1D* hn[100];
|
| 18 |
|
| 19 |
TH1D* hneb[100];
|
| 20 |
TH1D* hnee[100];
|
| 21 |
|
| 22 |
TH2D* h2o[100];
|
| 23 |
TH2D* h2a[100];
|
| 24 |
TH2D* h2b[100];
|
| 25 |
TH2D* h2c[100];
|
| 26 |
TH2D* h2d[100];
|
| 27 |
TH2D* h2pm[100];
|
| 28 |
TH2D* h2ebhf[100];
|
| 29 |
TH2D* h2z[100];
|
| 30 |
TH2D* h2zl[100];
|
| 31 |
TH2D* h2zpm[100];
|
| 32 |
TH2D* h2zz[100];
|
| 33 |
TH2D* h2zzp[100];
|
| 34 |
TH2D* h2zzm[100];
|
| 35 |
TH2D* h2zlp[100];
|
| 36 |
TH2D* h2zlm[100];
|
| 37 |
|
| 38 |
for(int i = 0; i < 20; ++i){
|
| 39 |
ha[i] = new TH1D(Form("ha%d",i),";N HF towers above throshold;Events",100,0,1000);
|
| 40 |
hb[i] = new TH1D(Form("hb%d",i),";Sum Tower E_{T} [GeV];Events",100,0,4000);
|
| 41 |
hc[i] = new TH1D(Form("hc%d",i),";Sum HF E_{T} [GeV];Events",100,0,5000);
|
| 42 |
|
| 43 |
|
| 44 |
he[i] = new TH1D(Form("he%d",i),"",1000,0,4000);
|
| 45 |
|
| 46 |
h2z[i] = new TH2D(Form("h2z%d",i),";Sum HF E_{T} [GeV]; Sum ZDC Energy [TeV]",100,0,5000,100,0,350);
|
| 47 |
h2zpm[i] = new TH2D(Form("h2zpm%d",i),";Sum ZDC- Energy [TeV]; Sum ZDC+ Energy [TeV]",100,-5,200,100,-5,200);
|
| 48 |
h2zz[i] = new TH2D(Form("h2zz%d",i),";Sum HF E_{T} [GeV]; Sum ZDC Energy [TeV]",100,0,500,100,0,350);
|
| 49 |
h2zzp[i] = new TH2D(Form("h2zzp%d",i),";Sum HF- E_{T} [GeV]; Sum ZDC+ Energy [TeV]",100,0,500,100,0,200);
|
| 50 |
h2zzm[i] = new TH2D(Form("h2zzm%d",i),";Sum HF+ E_{T} [GeV]; Sum ZDC- Energy [TeV]",100,0,500,100,0,200);
|
| 51 |
h2zl[i] = new TH2D(Form("h2zl%d",i),";log_{e}(Sum HF E_{T} [GeV]); log_{e}(Sum ZDC Energy [GeV])",100,0,9,100,0,14);
|
| 52 |
|
| 53 |
h2zlp[i] = new TH2D(Form("h2zlp%d",i),";log_{e}(Sum HF E_{T} [GeV]); log_{e}(Sum ZDC+ Energy [GeV])",100,0,9,100,0,14);
|
| 54 |
|
| 55 |
h2zlm[i] = new TH2D(Form("h2zlm%d",i),";log_{e}(Sum HF E_{T} [GeV]); log_{e}(Sum ZDC- Energy [GeV])",100,0,9,100,0,14);
|
| 56 |
|
| 57 |
|
| 58 |
hn[i] = new TH1D(Form("hn%d",i),"",1000,0,10000);
|
| 59 |
|
| 60 |
hneb[i] = new TH1D(Form("hneb%d",i),";N_{hits} EE > 0.3 GeV; Events",100,0,10000);
|
| 61 |
hnee[i] = new TH1D(Form("hnee%d",i),";N_{hits} EE > 0.3 GeV; Events",100,0,10000);
|
| 62 |
|
| 63 |
h2o[i] = new TH2D(Form("h2o%d",i),"",420,0,420,600,0,6000);
|
| 64 |
h2a[i] = new TH2D(Form("h2a%d",i),"",420,0,420,900,0,9000);
|
| 65 |
h2b[i] = new TH2D(Form("h2b%d",i),"",900,0,9000,600,0,6000);
|
| 66 |
h2c[i] = new TH2D(Form("h2c%d",i),"",420,0,420,370,0,3700);
|
| 67 |
h2d[i] = new TH2D(Form("h2d%d",i),"",370,0,3700,600,0,6000);
|
| 68 |
|
| 69 |
h2pm[i] = new TH2D(Form("h2pm%d",i),"; E_{T} Sum HF- [GeV]; E_{T} Sum HF+ [GeV]",100,0,2500,100,0,2500);
|
| 70 |
h2ebhf[i] = new TH2D(Form("h2ebhf%d",i),"; E_{T} Sum HF [GeV]; E_{T} Sum EB [GeV]",100,0,5000,100,0,500);
|
| 71 |
|
| 72 |
}
|
| 73 |
|
| 74 |
// ha[1]->SetLineColor(2);
|
| 75 |
// ha[1]->SetMarkerColor(2);
|
| 76 |
|
| 77 |
TTree* t1 = (TTree*)inf->Get("hltanalysis/HltTree");
|
| 78 |
|
| 79 |
TCanvas* c1 = new TCanvas("c1","",1200,800);
|
| 80 |
c1->Divide(5,3);
|
| 81 |
|
| 82 |
c1->cd(5);
|
| 83 |
t1->Draw(Form("hiHFplus:hiHFminus>>h2pm0"),trig,"colz");
|
| 84 |
|
| 85 |
c1->cd(6);
|
| 86 |
t1->Draw(Form("hiEB:hiHF>>h2ebhf0"),trig,"colz");
|
| 87 |
|
| 88 |
c1->cd(7);
|
| 89 |
t1->Draw(Form("hiZDCplus/1000.:hiZDCminus/1000.>>h2zpm0"),trig,"colz");
|
| 90 |
|
| 91 |
|
| 92 |
c1->cd(8);
|
| 93 |
t1->Draw(Form("hiZDC/1000.:hiHF>>h2z0"),trig,"colz");
|
| 94 |
|
| 95 |
t1->Draw(Form("log(hiZDC):log(hiHF)>>h2zl0"),trig,"colz");
|
| 96 |
t1->Draw(Form("log(hiZDCplus):log(hiHF)>>h2zlp0"),trig,"colz");
|
| 97 |
t1->Draw(Form("log(hiZDCminus):log(hiHF)>>h2zlm0"),trig,"colz");
|
| 98 |
|
| 99 |
c1->cd(9);
|
| 100 |
t1->Draw(Form("hiZDC/1000.:hiHF>>h2zz0"),trig,"colz");
|
| 101 |
t1->Draw(Form("hiZDCplus/1000.:hiHFminus>>h2zzp0"),trig,"colz");
|
| 102 |
t1->Draw(Form("hiZDCminus/1000.:hiHFplus>>h2zzm0"),trig,"colz");
|
| 103 |
|
| 104 |
TCanvas* cc = new TCanvas("cc","",400,400);
|
| 105 |
h2pm0->Draw("colz");
|
| 106 |
cc->SetLogz();
|
| 107 |
cc->Print("HFplusMinus.gif");
|
| 108 |
|
| 109 |
TCanvas* cc2 = new TCanvas("cc2","",400,400);
|
| 110 |
cc2->SetLogz();
|
| 111 |
|
| 112 |
h2ebhf0->Draw("colz");
|
| 113 |
cc2->Print("HFvsEB.gif");
|
| 114 |
|
| 115 |
TCanvas* cc3 = new TCanvas("cc3","",400,400);
|
| 116 |
cc3->SetLogz();
|
| 117 |
|
| 118 |
h2zpm0->Draw("colz");
|
| 119 |
cc3->Print("ZDCpm.gif");
|
| 120 |
|
| 121 |
TCanvas* cc4 = new TCanvas("cc4","",400,400);
|
| 122 |
cc4->SetLogz();
|
| 123 |
h2z0->SetNdivisions(5);
|
| 124 |
h2z0->Draw("colz");
|
| 125 |
cc4->Print("ZDChf.gif");
|
| 126 |
|
| 127 |
TCanvas* cc5 = new TCanvas("cc5","",400,400);
|
| 128 |
cc5->SetLogz();
|
| 129 |
h2zz0->Draw("colz");
|
| 130 |
cc5->Print("ZDChfZoom.gif");
|
| 131 |
|
| 132 |
TCanvas* ccpz5 = new TCanvas("ccpz5","",400,400);
|
| 133 |
ccpz5->SetLogz();
|
| 134 |
h2zzp0->Draw("colz");
|
| 135 |
ccpz5->Print("ZDCplusHFZoom.gif");
|
| 136 |
|
| 137 |
TCanvas* ccmz5 = new TCanvas("ccmz5","",400,400);
|
| 138 |
ccmz5->SetLogz();
|
| 139 |
h2zzm0->Draw("colz");
|
| 140 |
ccmz5->Print("ZDCminusHFZoom.gif");
|
| 141 |
|
| 142 |
|
| 143 |
TCanvas* ccl5 = new TCanvas("ccl5","",400,400);
|
| 144 |
ccl5->SetLogz();
|
| 145 |
h2zl0->Draw("colz");
|
| 146 |
ccl5->Print("ZDChfLog.gif");
|
| 147 |
|
| 148 |
TCanvas* ccl5p = new TCanvas("cclp5","",400,400);
|
| 149 |
ccl5p->SetLogz();
|
| 150 |
h2zlp0->Draw("colz");
|
| 151 |
ccl5p->Print("ZDChfLogPlus.gif");
|
| 152 |
|
| 153 |
TCanvas* ccl5m = new TCanvas("cclm5","",400,400);
|
| 154 |
ccl5m->SetLogz();
|
| 155 |
h2zlm0->Draw("colz");
|
| 156 |
ccl5m->Print("ZDChfLogMinus.gif");
|
| 157 |
|
| 158 |
outf->Write();
|
| 159 |
|
| 160 |
}
|
| 161 |
|
