| 201 |
|
float staterr=uncertainties[isample*jzbcuts.size()+ibin][2]; |
| 202 |
|
float systerr=uncertainties[isample*jzbcuts.size()+ibin][3]; |
| 203 |
|
float toterr =uncertainties[isample*jzbcuts.size()+ibin][4]; |
| 204 |
< |
float observed,null,result; |
| 205 |
< |
fill_result_histos(observed, null,null,null,null,null,null,null,mcjzb,JZBcutat,(int)5,result,(signalsamples.FindSample(signalsamples.collection[isample].filename)),signalsamples); |
| 204 |
> |
float observed,observederr,null,result; |
| 205 |
> |
fill_result_histos(observed,observederr, null,null,null,null,null,null,null,mcjzb,JZBcutat,14000,(int)5,result,(signalsamples.FindSample(signalsamples.collection[isample].filename)),signalsamples); |
| 206 |
|
observed-=result;//this is the actual excess we see! |
| 207 |
|
float expected=observed/luminosity; |
| 208 |
|
|
| 442 |
|
limfile->Close(); |
| 443 |
|
write_info("prepare_limits","limitfile.root and datacard.txt have been generated. You can now use them to calculate limits!"); |
| 444 |
|
|
| 445 |
< |
} |
| 445 |
> |
} |
