137 |
|
return histo; |
138 |
|
} |
139 |
|
|
140 |
+ |
TH2F *hardlimit(TH2F *limit) { |
141 |
+ |
TH2F *co = (TH2F*)limit->Clone("hcopy"); |
142 |
+ |
for(int i=1;i<limit->GetNbinsX();i++) { |
143 |
+ |
for(int j=1;j<limit->GetNbinsY();j++) { |
144 |
+ |
if(limit->GetBinContent(i,j)<1&&limit->GetBinContent(i,j)>0) co->SetBinContent(i,j,1); |
145 |
+ |
else co->SetBinContent(i,j,0); |
146 |
+ |
} |
147 |
+ |
} |
148 |
+ |
return co; |
149 |
+ |
} |
150 |
+ |
|
151 |
|
TH2F* make_exclusion_shape(TH2F *excl, int isprimary) { |
152 |
|
TH2F *exclusion = (TH2F*)excl->Clone("exclusion"); |
153 |
|
for(int i=1; i<(excl->GetNbinsX()+1); i++) { |
163 |
|
} |
164 |
|
|
165 |
|
void produce_extensive_plots(TH2F *xsec,TH2F *limits,TH2F *rellimits, TH2F *rellimitsd3, TH2F *rellimitst3, int scantype) { |
166 |
< |
TCanvas *ca = new TCanvas("ca","ca",1200,1200); |
167 |
< |
ca->Divide(2,2); |
166 |
> |
TCanvas *ca = new TCanvas("ca","ca",2400,1200); |
167 |
> |
ca->Divide(4,2); |
168 |
|
ca->cd(1); |
169 |
|
ca->cd(1)->SetLogz(1); |
170 |
|
xsec->GetZaxis()->SetRangeUser(0.001,1000); |
217 |
|
partial << "Limits/exclusion__" << limits->GetName(); |
218 |
|
fill_with_text(exclline,excllined3,excllinet3,ca->cd(4),scantype); |
219 |
|
CompleteSave(ca,partial.str()); |
220 |
+ |
|
221 |
+ |
ca->cd(5); |
222 |
+ |
(hardlimit(rellimits))->Draw("COL"); |
223 |
+ |
exclline->Draw("same"); |
224 |
+ |
ca->cd(6); |
225 |
+ |
(hardlimit(rellimitst3))->Draw("COL"); |
226 |
+ |
excllinet3->Draw("same"); |
227 |
+ |
ca->cd(7); |
228 |
+ |
(hardlimit(rellimitsd3))->Draw("COLZ"); |
229 |
+ |
excllined3->Draw("same"); |
230 |
+ |
|
231 |
+ |
CompleteSave(ca,partial.str()+"__PlusInfo"); |
232 |
|
delete ca; |
233 |
|
} |
234 |
|
|