45 |
|
ntMix->SetAlias("phi1","jtphi1"); |
46 |
|
ntMix->SetAlias("phi2","jtphi2"); |
47 |
|
|
48 |
< |
nt->SetAlias("cent","hf"); |
49 |
< |
ntMix->SetAlias("cent","hf"); |
48 |
> |
nt->SetAlias("cent","hfp"); |
49 |
> |
ntMix->SetAlias("cent","hfp"); |
50 |
|
|
51 |
|
|
52 |
|
nt->SetAlias("w","pt1> -999"); |
56 |
|
|
57 |
|
TTree* t[2] = {nt,ntMix}; |
58 |
|
|
59 |
+ |
TH2D* h[2]; |
60 |
+ |
TH1D* hMean[3]; |
61 |
+ |
TH1D* hSigma[3]; |
62 |
+ |
|
63 |
|
|
64 |
|
TCanvas* c1 = new TCanvas("c1","",1200,600); |
65 |
|
c1->Divide(2,1); |
72 |
|
|
73 |
|
t[i]->Draw(Form("ntrk:cent>>%s",h[i]->GetName()),weight*dijet&&deltaPhi&&jetID&&vtx&&noise,"colz"); |
74 |
|
|
75 |
+ |
fitSlices(h[i]); |
76 |
+ |
hMean[i] = (TH1D*)gDirectory->Get(Form("%s_1",h[i]->GetName())); |
77 |
+ |
hSigma[i] = (TH1D*)gDirectory->Get(Form("%s_2",h[i]->GetName())); |
78 |
+ |
|
79 |
|
} |
80 |
|
|
81 |
|
|