712 |
|
return gStyle->GetColorPalette(theColor); |
713 |
|
} |
714 |
|
|
715 |
+ |
float findmaxentry(vector<TH2F*> histos, int i, int j) { |
716 |
+ |
float max = histos[0]->GetBinContent(i,j); |
717 |
+ |
for(int k=0;k<histos.size();k++) { |
718 |
+ |
float entry = histos[k]->GetBinContent(i,j); |
719 |
+ |
if(entry>=max) max=entry; |
720 |
+ |
} |
721 |
+ |
|
722 |
+ |
return max; |
723 |
+ |
} |
724 |
|
|
725 |
|
TH2F *make_best_limits(vector<TH2F*> explimits,vector<TH2F*> obslimits, int scantype, vector<TH2F*> exp1mlimits, vector<TH2F*> exp1plimits, vector<TH2F*> exp2mlimits, vector<TH2F*> exp2plimits, vector<TH2F*> &allbestexplimits) { |
726 |
|
if(obslimits.size()==0) { |
742 |
|
|
743 |
|
for(int i=1;i<=explimits[0]->GetNbinsX();i++) { |
744 |
|
for(int j=1;j<=explimits[0]->GetNbinsY();j++) { |
745 |
< |
float min=explimits[0]->GetBinContent(i,j); |
745 |
> |
float min=findmaxentry(explimits,i,j); |
746 |
|
float omin=obslimits[0]->GetBinContent(i,j); |
747 |
< |
int source=1; |
747 |
> |
int source=0; |
748 |
|
for(int k=0;k<explimits.size();k++) { |
749 |
|
float currlim=explimits[k]->GetBinContent(i,j); |
750 |
|
if(currlim<min&&currlim>0) { |