ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/Modules/ExclusionPlot.C
(Generate patch)

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/ExclusionPlot.C (file contents):
Revision 1.16 by buchmann, Tue Dec 13 09:54:38 2011 UTC vs.
Revision 1.17 by buchmann, Tue Dec 13 10:44:44 2011 UTC

# Line 712 | Line 712 | Int_t get_exclusion_region_color(double
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) {
# Line 733 | Line 742 | TH2F *make_best_limits(vector<TH2F*> exp
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) {

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines