| 503 |
|
float mceff,mcefferr,jzbeff,jzbefferr; |
| 504 |
|
if(!automatized) dout << "MC efficiencies:" << endl; |
| 505 |
|
Value mceff_nosigcont = MCefficiency(events,mceff,mcefferr,mcjzb,requireZ,Neventsinfile,addcut,-1); |
| 506 |
< |
cout << "no signal contamionation result : " << mceff_nosigcont << " which means value = " << mceff_nosigcont.getValue() << endl; |
| 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 |
< |
results.push_back(res); |
| 579 |
< |
cout << mceff_nosigcont.getValue() << endl; |
| 579 |
< |
// results.push_back(mceff_nosigcont.getValue()); |
| 580 |
< |
// results.push_back(mceff_nosigcont.getError()); |
| 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 |
|
} |
| 583 |
|
|
| 584 |
|
vector<vector<float> > compute_systematics(string mcjzb, float mcpeakerror, string datajzb, samplecollection &signalsamples, vector<float> bins, bool requireZ=false) { |
