| 52 |
|
nt->SetAlias("w","pt1> -999"); |
| 53 |
|
ntMix->SetAlias("w","weight*(pthat > 50)"); |
| 54 |
|
|
| 55 |
– |
TCut lead(Form("pt1>%d && abs(eta1) < 2",leadCut)); |
| 56 |
– |
TCut dijet(Form("pt1>%d && pt2>%d && abs(eta1) < 2 && abs(eta2) < 2",leadCut,subleadCut)); |
| 57 |
– |
|
| 58 |
– |
TCut deltaPhi("abs(dphi)>2.0944"); |
| 59 |
– |
TCut jetID("trkMax1 > 4 || trkMax2 > 4"); |
| 60 |
– |
TCut noise("noise < 0"); |
| 55 |
|
TCut weight("w"); |
| 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); |
| 70 |
|
h[i]->GetXaxis()->CenterTitle(); |
| 71 |
|
h[i]->GetYaxis()->CenterTitle(); |
| 72 |
|
|
| 73 |
< |
t[i]->Draw(Form("ntrk:cent>>%s",h[i]->GetName()),weight*dijet&&deltaPhi&&jetID&&noise,"colz"); |
| 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 |
|
|
