| 266 |
|
stringstream partial; |
| 267 |
|
partial << "Limits/exclusion__" << limits->GetName(); |
| 268 |
|
fill_with_text(exclline,excllined3,excllinet3,ca->cd(4),scantype); |
| 269 |
< |
CompleteSave(ca,partial.str()); |
| 269 |
> |
// CompleteSave(ca,partial.str()); |
| 270 |
|
|
| 271 |
|
ca->cd(5); |
| 272 |
|
(hardlimit(rellimits))->Draw("COL"); |
| 273 |
|
exclline->Draw("same"); |
| 274 |
+ |
TText *a = write_title("Exclusion shape for 3x#sigma_{ref}"); |
| 275 |
+ |
a->Draw(); |
| 276 |
|
ca->cd(6); |
| 277 |
|
(hardlimit(rellimitst3))->Draw("COL"); |
| 278 |
|
excllinet3->Draw("same"); |
| 279 |
+ |
TText *b = write_title("Exclusion shape for #sigma_{ref}"); |
| 280 |
+ |
b->Draw(); |
| 281 |
|
ca->cd(7); |
| 282 |
|
(hardlimit(rellimitsd3))->Draw("COLZ"); |
| 283 |
+ |
TText *c = write_title("Exclusion shape for #sigma_{ref}/3"); |
| 284 |
+ |
c->Draw(); |
| 285 |
|
excllined3->Draw("same"); |
| 286 |
|
|
| 287 |
|
CompleteSave(ca,partial.str()+"__PlusInfo"); |
| 871 |
|
int mappoint=name.find("map"); |
| 872 |
|
if(mappoint<0||mappoint>500) return name; // this mean that something weird is happening |
| 873 |
|
stringstream nice_label; |
| 874 |
< |
nice_label << "JZB > " << name.substr(mappoint+3,name.size()) << " GeV"; |
| 874 |
> |
string identifier=name.substr(mappoint+3,name.size()); |
| 875 |
> |
if(identifier!="0") nice_label << "JZB > " << name.substr(mappoint+3,name.size()) << " GeV"; |
| 876 |
> |
else nice_label << "JZB shape"; |
| 877 |
|
return nice_label.str(); |
| 878 |
|
} |
| 879 |
|
|
| 967 |
|
TLegend *sourceleg = new TLegend(0.2,0.5,0.55,0.75); |
| 968 |
|
for(int i=0;i<(int)explimits.size();i++) { |
| 969 |
|
stringstream legendentry; |
| 970 |
< |
legendentry << i+1 << " = " << give_nice_source_label(explimits[i]->GetName()); |
| 970 |
> |
legendentry << give_nice_source_label(explimits[i]->GetName()); |
| 971 |
|
Int_t currcol=get_exclusion_region_color(i+1,bestlimitsource); |
| 972 |
|
explimits[i]->SetFillColor(currcol); |
| 973 |
|
explimits[i]->SetLineColor(currcol); |
| 1222 |
|
} |
| 1223 |
|
} |
| 1224 |
|
|
| 1225 |
+ |
string IdentifyScan(TString histoname) { |
| 1226 |
+ |
cout << "We are dealing with a histo called " << histoname << endl; |
| 1227 |
+ |
if (histoname.Contains("T5zzl")) return "0.75"; |
| 1228 |
+ |
if(histoname.Contains("T5zzh")) return "0.25"; |
| 1229 |
+ |
if(histoname.Contains("T5zz")) return "0.5"; |
| 1230 |
+ |
|
| 1231 |
+ |
return "IdentifyScanError"+string(histoname.Data()); |
| 1232 |
+ |
} |
| 1233 |
+ |
|
| 1234 |
|
void process_file(TFile* file, float stdmargin) { |
| 1235 |
|
standardmargin=stdmargin; |
| 1236 |
|
xsecfilename=PlottingSetup::cbafbasedir+"/"+PlottingSetup::SMSReferenceXSFile; |
| 1241 |
|
|
| 1242 |
|
int scantype=PlottingSetup::SMS; |
| 1243 |
|
std::string scanx = "0"; // Just for the legend |
| 1227 |
– |
if (TString(file->GetName()).Contains("T5zzl")) scanx = "0.75"; |
| 1228 |
– |
else if(TString(file->GetName()).Contains("T5zzh")) scanx = "0.25"; |
| 1229 |
– |
else if(TString(file->GetName()).Contains("T5zz")) scanx = "0.5"; |
| 1244 |
|
|
| 1245 |
|
vector<TH2F*> systematics_histos; |
| 1246 |
|
vector<TH2F*> limits_histos; |
| 1285 |
|
else if(is_systematic) systematics_histos.push_back((TH2F*) obj); |
| 1286 |
|
if(name.Contains("mSUGRA")) scantype=PlottingSetup::mSUGRA; |
| 1287 |
|
if(name.Contains("GMSB")) scantype=PlottingSetup::GMSB; |
| 1288 |
+ |
if(name.Contains("ScanIdentifier")) scanx=IdentifyScan(name); |
| 1289 |
|
} |
| 1290 |
+ |
// if (TString(file->GetName()).Contains("T5zzl")) scanx = "0.75"; |
| 1291 |
+ |
// else if(TString(file->GetName()).Contains("T5zzh")) scanx = "0.25"; |
| 1292 |
+ |
// else if(TString(file->GetName()).Contains("T5zz")) scanx = "0.5"; |
| 1293 |
|
write_warning(__FUNCTION__,"Deactivated systematics plots"); |
| 1294 |
|
// if(systematics_histos.size()>0) make_all_syst_plots(systematics_histos,scantype,scanx); |
| 1295 |
|
if(limits_histos.size()>0) create_exclusion_plots(limits_histos,scantype,scanx); |
