| 309 |
|
} |
| 310 |
|
|
| 311 |
|
char jzbSelStr[256]; sprintf(jzbSelStr,"%f",jzbSel); |
| 312 |
+ |
char metSel[256]; sprintf(metSel, "met[4] > %f", jzbSel); |
| 313 |
+ |
string metSelection(metSel); |
| 314 |
|
// All acceptance cuts at gen. level |
| 315 |
|
//TCut kbase("abs(genMll-91.2)<20&&genNjets>2&&genZPt>0&&genJZB"+geq_or_leq()+TString(jzbSelStr)+"&&genId1==-genId2"); |
| 316 |
|
TCut kbase(""); |
| 323 |
|
TCut njets(cutnJets); |
| 324 |
|
TCut jzbp; |
| 325 |
|
TCut jzbn; |
| 326 |
+ |
TCut met(("pfJetGoodNum > 1 && abs(mll-91.2) < 10.0 && id1 == id2 &&" + metSelection).c_str()); |
| 327 |
|
if(flipped==0) { |
| 328 |
|
jzbp=TCut((TString(mcjzb)+geq_or_leq()+TString(jzbSelStr))); |
| 329 |
|
jzbn=TCut((TString(mcjzb)+ngeq_or_leq()+TString("-")+TString(jzbSelStr))); |
| 331 |
|
jzbp=TCut(TString(mcjzb)+ngeq_or_leq()+TString("-")+TString(jzbSelStr)); |
| 332 |
|
jzbn=TCut(TString(mcjzb)+geq_or_leq()+TString(jzbSelStr)); |
| 333 |
|
} |
| 334 |
< |
float ntotal = events->Draw("pt1", kbase, "goff"); |
| 334 |
> |
float ntotal = events->Draw("pt1", addcut.c_str(), "goff"); |
| 335 |
|
TCut theCut; |
| 336 |
|
switch(type) { |
| 337 |
|
case 1: |
| 346 |
|
case 4: |
| 347 |
|
theCut = kbase+massId+njets+jzbn; |
| 348 |
|
break; |
| 349 |
+ |
case 5: |
| 350 |
+ |
theCut = kbase + met; |
| 351 |
+ |
break; |
| 352 |
|
default: |
| 353 |
|
theCut = kbase+massId+njets+jzbn; |
| 354 |
|
break; |
| 355 |
|
} |
| 356 |
< |
|
| 356 |
> |
|
| 357 |
|
string stheCut(theCut); |
| 358 |
|
char var[20]; |
| 359 |
|
sprintf(var, "pdfW[%d]", k); |
