| 213 |
|
- JZB cut [9] |
| 214 |
|
- plot name [10]*/ |
| 215 |
|
|
| 216 |
< |
cout << "Calling limit capsule instead of calling : CL95(" << luminosity << "," << lumiuncert*luminosity << "," << mceff << "," << mcefferr << "," << Npred[ibin] << "," << Nprederr[ibin] << "," << Nobs[ibin] << "," << false << "," << nuisancemodel<< ") " << endl; |
| 216 |
> |
dout << "Calling limit capsule instead of calling : CL95(" << luminosity << "," << lumiuncert*luminosity << "," << mceff << "," << mcefferr << "," << Npred[ibin] << "," << Nprederr[ibin] << "," << Nobs[ibin] << "," << false << "," << nuisancemodel<< ") " << endl; |
| 217 |
|
|
| 218 |
|
stringstream command; |
| 219 |
|
command << PlottingSetup::cbafbasedir << "/DistributedModelCalculations/Limits/LimitCapsule.exec " << repname.str() << " " << luminosity << " " << luminosity*lumiuncert << " " << mceff << " " << mcefferr << " " << Npred[ibin] << " " << Nprederr[ibin] << " " << Nobs[ibin] << " " << -1 << " " << plotfilename; |
| 220 |
< |
std::cout << command.str() << std::endl; |
| 220 |
> |
dout << command.str() << endl; |
| 221 |
|
|
| 222 |
|
int retval = 256; |
| 223 |
|
int attempts=0; |
| 224 |
|
while(!(retval==0||attempts>=5)) {//try up to 5 times |
| 225 |
|
attempts++; |
| 226 |
< |
cout << "Starting limit calculation (LimitCapsule) now : Attempt " << attempts << endl; |
| 226 |
> |
dout << "Starting limit calculation (LimitCapsule) now : Attempt " << attempts << endl; |
| 227 |
|
retval=gSystem->Exec(command.str().c_str()); |
| 228 |
|
} |
| 229 |
|
|
| 230 |
|
LimitDroplet limres; |
| 231 |
|
limres.readDroplet(repname.str()); |
| 232 |
< |
cout << limres << endl; |
| 232 |
> |
dout << limres << endl; |
| 233 |
|
remove(repname.str().c_str()); |
| 234 |
|
sigma95=limres.expected; |
| 235 |
|
|
| 336 |
|
if(!doobserved) { |
| 337 |
|
dout << endl << endl << "LIMITS: (Tex)" << endl; |
| 338 |
|
tout << "\\begin{table}[hbtp]" << endl; |
| 339 |
< |
tout << "\\renewcommand{\arraystretch}{1.3}" << endl; |
| 339 |
> |
tout << "\\renewcommand{\\arraystretch}{1.3}" << endl; |
| 340 |
|
tout << "\\begin{center}" << endl; |
| 341 |
|
tout << "\\caption{Observed upper limits on the cross section of different LM benchmark points " << (ConsiderSignalContaminationForLimits?" (accounting for signal contamination)":" (not accounting for signal contamination)") << "}\\label{tab:lmresults}" << endl; |
| 342 |
|
tout << "" << endl; |
