]>
Commit | Line | Data |
---|---|---|
37017262 | 1 | void PlotMatchedTracksFromData(TString datafilename="rootFiles/LHC11a4_bis/Et.ESD.simPbPb.EMCAL.LHC11a4_bis.root",TString simfilename="rootFiles/LHC11a4_bis/Et.ESD.simPbPb.EMCAL.LHC11a4_bis.root", int bin = 10, int binLast = 10, TString det = "EMCal"){ |
2 | gStyle->SetOptTitle(0); | |
3 | gStyle->SetOptStat(0); | |
4 | gStyle->SetOptFit(0); | |
5 | TString outname = ""; | |
6 | TString outnamebin = Form("%iTo%i",bin,binLast); | |
7 | ||
8 | TFile *fsim = TFile::Open(simfilename, "READ"); | |
9 | TList *lsim = (TList*)fsim->Get("out1"); | |
10 | TH3F *fHistMatchedTracksEvspTBkgdvsMult = lsim->FindObject("fHistMatchedTracksEvspTBkgdMult"); | |
11 | fHistMatchedTracksEvspTBkgdvsMult->GetZaxis()->SetRange(bin,binLast); | |
12 | TH2D *hBkgd2D = (TProfile2D*) fHistMatchedTracksEvspTBkgdvsMult->Project3D("yx"); | |
13 | TProfile * profBkgd2D = hBkgd2D->ProfileX(); | |
14 | profBkgd2D->SetLineColor(2); | |
15 | profBkgd2D->SetLineWidth(2); | |
16 | ||
17 | TFile *f = TFile::Open(datafilename, "READ"); | |
18 | ||
19 | f->cd(); | |
20 | ||
21 | TList *l = (TList*)f->Get("out1"); | |
22 | TH3F *fHistMatchedTracksEvspTvsMult = l->FindObject("fHistMatchedTracksEvspTvsMult"); | |
23 | fHistMatchedTracksEvspTvsMult->GetZaxis()->SetRange(bin,binLast); | |
24 | TH2D *hMatchedTracks2D = (TProfile2D*) fHistMatchedTracksEvspTvsMult->Project3D("yx"); | |
25 | TProfile * profMatchedTracks2D = hMatchedTracks2D->ProfileX(); | |
26 | hMatchedTracks2D->GetXaxis()->SetTitle("p"); | |
27 | hMatchedTracks2D->GetYaxis()->SetTitle("E^{cluster}"); | |
28 | profMatchedTracks2D->SetLineWidth(2); | |
29 | TCanvas *c1 = new TCanvas("c1","c1",600,400); | |
30 | c1->SetTopMargin(0.02); | |
31 | c1->SetRightMargin(0.02); | |
32 | c1->SetBorderSize(0); | |
33 | c1->SetFillColor(0); | |
34 | c1->SetFillColor(0); | |
35 | c1->SetBorderMode(0); | |
36 | c1->SetFrameFillColor(0); | |
37 | c1->SetFrameBorderMode(0); | |
38 | hMatchedTracks2D->Draw("colz"); | |
39 | profMatchedTracks2D->Draw("same"); | |
40 | profBkgd2D->Draw("same"); | |
41 | outname = "/tmp/TrackMatchingData2D"+det+outnamebin+".png"; | |
42 | c1->SaveAs(outname.Data()); | |
43 | } |