| 265 |
|
} else { |
| 266 |
|
dout << "Configuration (mglu="<<mglu<<" , mlsp="<<mlsp << ") passed the mcefficiency test." << endl; |
| 267 |
|
} |
| 268 |
< |
if(efficiencyonly) { |
| 269 |
< |
efficiencymap->Fill(mglu,mlsp,result); |
| 270 |
< |
Neventsmap->Fill(mglu,mlsp,nevents); |
| 271 |
< |
continue; |
| 272 |
< |
} |
| 268 |
> |
efficiencymap->Fill(mglu,mlsp,result); |
| 269 |
> |
Neventsmap->Fill(mglu,mlsp,nevents); |
| 270 |
> |
if(efficiencyonly) continue; |
| 271 |
> |
|
| 272 |
|
do_systematics_for_one_file((scansample.collection)[0].events,nevents,"SUSY SCAN", systematics,mcjzb,datajzb,peakerror,requireZ, addcut.str()); |
| 273 |
|
float JZBcutat = systematics[0][0]; |
| 274 |
|
float mceff = systematics[0][1]; |
| 290 |
|
cout << "back in " << __FUNCTION__ << endl; |
| 291 |
|
if(sigmas[0]>-0.5) { // negative sigmas are the error signature of do_limit_wrapper, so we want to exclude them. |
| 292 |
|
limitmap->Fill(mglu,mlsp,sigmas[0]); |
| 293 |
+ |
sysjesmap->Fill(mglu,mlsp,sys_jes); |
| 294 |
+ |
sysjsumap->Fill(mglu,mlsp,sys_jsu); |
| 295 |
+ |
sysresmap->Fill(mglu,mlsp,sys_res); |
| 296 |
+ |
efficiencymap->Fill(mglu,mlsp,result); |
| 297 |
+ |
Neventsmap->Fill(mglu,mlsp,nevents); |
| 298 |
|
dout << "A limit has been added at " << sigmas[0] << " for m_{glu}="<<mglu << " and m_{lsp}="<<mlsp<<endl; |
| 299 |
|
} //end of if sigma is positive |
| 300 |
|
//end of not systematics only condition |
| 303 |
|
sysjesmap->Fill(mglu,mlsp,sys_jes); |
| 304 |
|
sysjsumap->Fill(mglu,mlsp,sys_jsu); |
| 305 |
|
sysresmap->Fill(mglu,mlsp,sys_res); |
| 306 |
+ |
efficiencymap->Fill(mglu,mlsp,result); |
| 307 |
+ |
Neventsmap->Fill(mglu,mlsp,nevents); |
| 308 |
|
} |
| 309 |
|
}//efficiency is valid |
| 310 |
|
} |
