ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/TheLimits.cc
(Generate patch)

Comparing UserCode/auterman/SusyScan/PlotScript/TheLimits.cc (file contents):
Revision 1.4 by auterman, Sat Nov 27 12:22:28 2010 UTC vs.
Revision 1.6 by auterman, Sat Jan 29 11:00:19 2011 UTC

# Line 6 | Line 6
6   #include <iostream>
7   #include <cmath>
8  
9 + void TheLimits::Fill(int argc, char** argv)
10 + {
11 +   for (int i = 1; i<argc; ++i)
12 +   {
13 +     add( new SusyScan(argv[i]) );
14 +   }
15 + }
16 +
17 + void TheLimits::Fill(const std::string filelist)
18 + {
19 +   std::ifstream masses_file;
20 +   masses_file.open(filelist.c_str());
21 +   std::string file;
22 +   while (1) {
23 +      GeneratorMasses * p = new GeneratorMasses;
24 +      masses_file >> file;
25 +      if (!masses_file.good()) break;
26 +      add( new SusyScan(file));
27 +   }
28 +   masses_file.close();
29 + }
30  
31   void TheLimits::match()
32   {
# Line 104 | Line 125 | void TheLimits::FillGeneratorMasses(std:
125        _masses.push_back( p );
126     }
127   }
128 +
129 +
130 + void TheLimits::OverwriteLimits(std::string flag)
131 + {
132 +  for (std::vector<SusyScan*>::iterator it=_scan.begin();it!=_scan.end();++it){
133 +    if (flag=="ABCD_MHT") {
134 +      (*it)->ExpNsigLimit = 17.85;
135 +      (*it)->PLExpNsigLimit = 12.558;
136 +      (*it)->FCExpNsigLimit = 15.75;
137 +      (*it)->MCMCExpNsigLimit = 15.8085;
138 +      (*it)->ObsNsigLimit = 10.0637;
139 +      (*it)->PLObsNsigLimit = 4.22124;
140 +      (*it)->FCObsNsigLimit  = 0.75;
141 +      (*it)->MCMCObsNsigLimit  = 8.75403;
142 +    }
143 +    else if (flag=="ABCD_HT") {
144 +      (*it)->ExpNsigLimit = 24.0438;
145 +      (*it)->PLExpNsigLimit = 17.8537;
146 +      (*it)->FCExpNsigLimit = 22.05;
147 +      (*it)->MCMCExpNsigLimit = 20.781;
148 +      (*it)->ObsNsigLimit    = 18.24;
149 +      (*it)->PLObsNsigLimit = 9.70375;
150 +      (*it)->FCObsNsigLimit  = 11.85;
151 +      (*it)->MCMCObsNsigLimit  = 14.34;
152 +    }
153 +    
154 +    
155 +      (*it)->ExpXsecLimit    = (*it)->ExpNsigLimit * (*it)->Xsection / (*it)->signal;
156 +      (*it)->PLExpXsecLimit  = (*it)->PLExpNsigLimit  * (*it)->Xsection / (*it)->signal;
157 +      (*it)->FCExpXsecLimit  = (*it)->FCExpNsigLimit  * (*it)->Xsection / (*it)->signal;
158 +      (*it)->MCMCExpXsecLimit= (*it)->MCMCExpNsigLimit  * (*it)->Xsection / (*it)->signal;
159 +      (*it)->ObsXsecLimit    = (*it)->ObsNsigLimit  * (*it)->Xsection / (*it)->signal;
160 +      (*it)->PLObsXsecLimit  = (*it)->PLObsNsigLimit  * (*it)->Xsection / (*it)->signal;
161 +      (*it)->FCObsXsecLimit  = (*it)->FCObsNsigLimit   * (*it)->Xsection / (*it)->signal;
162 +      (*it)->MCMCObsXsecLimit= (*it)->MCMCObsNsigLimit  * (*it)->Xsection / (*it)->signal;
163 +    
164 +  
165 +  }
166 +
167 +
168 + }
169 +

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines