track_EEC/macros/draw_TrackEEC.C

0001 const float jet_pTs[6] = {10,20,30,40,60,1000};
0002 const int marks[5] = {20, 21, 33, 34, 47};
0003 const int colors[5] = {1, 632, 600, 418, 617};
0004
0005
0006 void draw_TrackEEC()
0007 {
0008     TFile *f = new TFile("/sphenix/tg/tg01/jets/bkimelman/track_EEC/trackDSTs/trackEECs_July15_2026_cuts.root","READ");
0009
0010     TH1D *inJet[5];
0011
0012     double maxC = 0.0;
0013
0014     for(int i=0; i<5; i++)
0015     {
0016         inJet[i] = (TH1D*)f->Get(std::format("hEEC_inJetAll_Q2_{:.0f}_{:.0f}",jet_pTs[i], jet_pTs[i+1]).c_str());
0017         //inJet[i] = (TH1D*)f->Get(std::format("hEEC_inJet_{:.0f}_{:.0f}",jet_pTs[i], jet_pTs[i+1]).c_str());
0018
0019         inJet[i]->Scale(1.0/inJet[i]->Integral());
0020         inJet[i]->Scale(1.0,"width");
0021
0022         for(int j=1; j<=inJet[i]->GetNbinsX(); j++)
0023         {
0024             double bc = inJet[i]->GetBinContent(j);
0025             if(bc > maxC) maxC = bc;
0026         }
0027
0028         inJet[i]->SetMarkerStyle(marks[i]);
0029         inJet[i]->SetMarkerColor(colors[i]);
0030         inJet[i]->SetLineColor(colors[i]);
0031     }
0032
0033     TH2D *h = new TH2D("h",";#Delta R;EEC",1,inJet[0]->GetXaxis()->GetXmin(),inJet[0]->GetXaxis()->GetXmax(),1,1e-3,1.5*maxC);
0034
0035     TCanvas *c1 = new TCanvas();
0036     gStyle->SetOptStat(0);
0037     c1->SetLogx();
0038     c1->SetLogy();
0039
0040     TLegend *leg = new TLegend(0.15,0.15,0.5,0.45);
0041
0042     h->Draw();
0043
0044     for(int i=0; i<5; i++)
0045     {
0046         inJet[i]->Draw("PSAME");
0047         leg->AddEntry(inJet[i],std::format("{:.0f} < p_{{T,jet}} < {:.0f} GeV",jet_pTs[i], jet_pTs[i+1]).c_str(),"p");
0048     }
0049
0050     leg->Draw("same");
0051
0052     c1->SaveAs("/sphenix/tg/tg01/jets/bkimelman/track_EEC/trackDSTs/inJetEEC_pT_July15_2026_cuts.pdf");
0053
0054
0055 }