| 354 |
|
if(sigmas[0]>-0.5) { // negative sigmas are the error signature of do_limit_wrapper, so we want to exclude them. |
| 355 |
|
limitmap->Fill(mglu,mlsp,sigmas[0]); |
| 356 |
|
dout << "A limit has been added at " << sigmas[0] << " for m_{glu}="<<mglu << " and m_{lsp}="<<mlsp<<endl; |
| 357 |
< |
|
| 358 |
< |
} else { //end of if sigma is positive |
| 357 |
> |
} //end of if sigma is positive |
| 358 |
> |
//end of not systematics only condition |
| 359 |
> |
} else { |
| 360 |
|
sysjesmap->Fill(mglu,mlsp,sys_jes); |
| 361 |
|
sysjsumap->Fill(mglu,mlsp,sys_jsu); |
| 362 |
|
sysresmap->Fill(mglu,mlsp,sys_res); |
| 363 |
< |
} |
| 363 |
< |
}//end of not systematics only condition |
| 363 |
> |
} |
| 364 |
|
} |
| 365 |
|
} |
| 366 |
|
} |
| 378 |
|
if(runninglocally) { |
| 379 |
|
CompleteSave(limcanvas,"SUSYScan/Limits_JZB_geq"+any2string(jzb_cut[ibin])); |
| 380 |
|
} else { |
| 381 |
< |
TFile *outputfile=new TFile(("output/DistributedLimits_job"+string(any2string(jobnumber))+"_of_"+string(any2string(njobs))+".root").c_str(),"UPDATE"); |
| 381 |
> |
TFile *outputfile=new TFile(("output/DistributedLimits_job"+string(any2string(jobnumber))+"_of_"+string(any2string(njobs))+".root").c_str(),"UPDATE");//needs to be "UPDATE" as we can get to this point for different JZB cuts and don't want to erase the previous data :-) |
| 382 |
|
limitmap->Write(); |
| 383 |
|
outputfile->Close(); |
| 384 |
|
} |
