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}; |
66 |
|
h[i]->GetXaxis()->CenterTitle(); |
67 |
|
h[i]->GetYaxis()->CenterTitle(); |
68 |
|
|
69 |
< |
t[i]->Draw(Form("ntrk:cent>>%s",h[i]->GetName()),weight*dijet&&deltaPhi&&jetID&&noise,"colz"); |
69 |
> |
t[i]->Draw(Form("ntrk:cent>>%s",h[i]->GetName()),weight*dijet&&deltaPhi&&jetID&&vtx&&noise,"colz"); |
70 |
|
|
71 |
|
} |
72 |
|
|