334 |
|
if(!automatized && k>0 ) dout << "PDF assessment: "; |
335 |
|
if(!automatized) dout << " MC efficiency: " << result << "+-" << resulterr << " ( JZB>" << jzbSel << " : " << sel << " , JZB<-" << jzbSel << " : " << nsel << " and nevents=" << tot << ") with normFact=" << normFactor << std::endl; |
336 |
|
delete effh; |
337 |
+ |
return result_wo_signalcont; |
338 |
|
} |
339 |
|
|
340 |
|
|
502 |
|
|
503 |
|
float mceff,mcefferr,jzbeff,jzbefferr; |
504 |
|
if(!automatized) dout << "MC efficiencies:" << endl; |
505 |
< |
MCefficiency(events,mceff,mcefferr,mcjzb,requireZ,Neventsinfile,addcut,-1); |
505 |
> |
Value mceff_nosigcont = MCefficiency(events,mceff,mcefferr,mcjzb,requireZ,Neventsinfile,addcut,-1); |
506 |
> |
if(!automatized) cout << " Without signal contamination, we find an efficiency of " << mceff_nosigcont << endl; |
507 |
> |
|
508 |
|
if(PlottingSetup::computeJZBefficiency) JZBefficiency(events,informalname,jzbeff,jzbefferr,requireZ,addcut); |
509 |
|
if(!automatized) dout << "JZB efficiency: " << jzbeff << "+/-" << jzbefferr << endl; |
510 |
|
|
575 |
|
if(fabs(jesup)>fabs(jesdown)) res.push_back(fabs(jesup)); else res.push_back(fabs(jesdown)); |
576 |
|
if(fabs(scaleup)>fabs(scaledown)) res.push_back(fabs(scaleup)); else res.push_back(fabs(scaledown)); |
577 |
|
res.push_back(fabs(resolution)); |
578 |
+ |
res.push_back(mceff_nosigcont.getValue()); |
579 |
+ |
res.push_back(mceff_nosigcont.getError()); |
580 |
|
if(ismSUGRA) res.push_back(PDFuncert); |
581 |
|
results.push_back(res); |
582 |
|
} |