66 |
|
lepcut.push_back("two10to15"); |
67 |
|
lepcut.push_back("two10to20"); |
68 |
|
lepcut.push_back("two20to7000"); |
69 |
+ |
lepcut.push_back("conf1010"); |
70 |
+ |
lepcut.push_back("conf1510"); |
71 |
+ |
lepcut.push_back("conf2010"); |
72 |
+ |
lepcut.push_back("conf1515"); |
73 |
+ |
lepcut.push_back("conf2015"); |
74 |
|
|
75 |
|
int Mass0min = 0; |
76 |
|
int Mass0max = 3000; |
134 |
|
float rjzb=(jzb[1]+0.06*pt-3.58612); |
135 |
|
int rpredobs; |
136 |
|
int rpredobssign; |
137 |
< |
int categories[6]; |
137 |
> |
int categories[lepcut.size()]; |
138 |
|
int mllcategories[mllcut.size()]; |
139 |
|
int jzbregion[jzbcut.size()]; |
140 |
|
int predobscat[predobs.size()]; |
151 |
|
else categories[4]=0; |
152 |
|
if((pt1>20)&&(pt2>20)) categories[5]=1; |
153 |
|
else categories[5]=0; |
154 |
+ |
if(pt1>10&&pt2>10) categories[6]=1; //conf1010 |
155 |
+ |
else categories[6]=0; |
156 |
+ |
if(pt1>15&&pt2>10) categories[7]=1; //conf1510 |
157 |
+ |
else categories[7]=0; |
158 |
+ |
if(pt1>20&&pt2>10) categories[8]=1; //conf2010 |
159 |
+ |
else categories[8]=0; |
160 |
+ |
if(pt1>15&&pt2>15) categories[9]=1; //conf1515 |
161 |
+ |
else categories[9]=0; |
162 |
+ |
if(pt1>20&&pt2>15) categories[10]=1; //conf2015 |
163 |
+ |
else categories[10]=0; |
164 |
|
|
165 |
|
if(!(rjzb>0&&id1==id2)) { |
166 |
|
predobscat[0]=1; |
241 |
|
} |
242 |
|
|
243 |
|
cout << "Have processed " << Nevents->GetEntries() << " events." << endl; |
244 |
< |
TFile *out = new TFile("NEWalloutputWithMll2.root","RECREATE"); |
244 |
> |
TFile *out = new TFile("NEWalloutputWithMll2complete.root","RECREATE"); |
245 |
|
for(int ip=0;ip<predobs.size();ip++) { |
246 |
|
for(int il=0;il<lepcut.size();il++) { |
247 |
|
for(int im=0;im<mllcut.size();im++) { |