| 255 |
|
|
| 256 |
|
if(!doobserved) { |
| 257 |
|
dout << endl << endl << "LIMITS: (Tex)" << endl; |
| 258 |
+ |
tout << "\\begin{table}[hbtp]" << endl; |
| 259 |
+ |
tout << "\\renewcommand{\arraystretch}{1.3}" << endl; |
| 260 |
+ |
tout << "\\begin{center}" << endl; |
| 261 |
+ |
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; |
| 262 |
+ |
tout << "" << endl; |
| 263 |
|
tout << "\\begin{tabular}{ | l | "; |
| 264 |
|
for (int irow=0;irow<jzbcuts.size();irow++) tout << " l |"; |
| 265 |
|
tout << "} " << endl << " \\hline " << endl << "& \t "; |
| 277 |
|
tout << " \\\\ \\hline " << endl; |
| 278 |
|
} |
| 279 |
|
tout << "\\end{tabular}" << endl; |
| 280 |
+ |
tout << " \\end{tabular}"<< endl; |
| 281 |
+ |
tout << "\\end{center}"<< endl; |
| 282 |
+ |
tout << "\\end{table} "<< endl; |
| 283 |
+ |
|
| 284 |
|
}//do observed |
| 285 |
|
|
| 286 |
|
dout << endl << endl << "Final selection efficiencies with total statistical and systematic errors, and corresponding observed and expected upper limits (UL) on ($\\sigma\\times$ BR $\\times$ acceptance) for the LM4 and LM8 scenarios, in the different regions. The last column contains the predicted ($\\sigma \\times $BR$\\times$ acceptance) at NLO obtained from Monte Carlo simulation." << endl; |
| 287 |
|
dout << "Scenario \t Efficiency [%] \t Upper limits [pb] \t \\sigma [pb]" << endl; |
| 288 |
|
for(int icut=0;icut<jzbcuts.size();icut++) { |
| 289 |
< |
dout << "Region with JZB>" << jzbcuts[icut] << endl; |
| 289 |
> |
dout << "Region with JZB>" << jzbcuts[icut] << (ConsiderSignalContaminationForLimits?" (accounting for signal contamination)":" (not accounting for signal contamination)") << endl; |
| 290 |
|
for(int isample=0;isample<signalsamples.collection.size();isample++) { |
| 291 |
|
dout << limits[isample][0] << "\t" << Round(100*uncertainties[isample*jzbcuts.size()+icut][1],3) << "+/-" << Round(100*uncertainties[isample*jzbcuts.size()+icut][2],3) << " (stat) +/- " << Round(100*uncertainties[isample*jzbcuts.size()+icut][3],3) << " (syst) \t" << Round((vlimits[isample][2*icut]),3) << "\t" << Round(vlimits[isample][2*icut+1],3) << endl; |
| 292 |
|
} |
