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.7 by auterman, Wed Feb 9 08:51:32 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 +   //std::cout << filelist << ": "<< _scan.size() <<std::endl;
29 +   masses_file.close();
30 + }
31  
32   void TheLimits::match()
33   {
# Line 104 | Line 126 | void TheLimits::FillGeneratorMasses(std:
126        _masses.push_back( p );
127     }
128   }
129 +
130 +
131 + void TheLimits::OverwriteLimits(std::string flag)
132 + {
133 +  for (std::vector<SusyScan*>::iterator it=_scan.begin();it!=_scan.end();++it){
134 +    if (flag=="ABCD_MHT") {
135 +      (*it)->ExpNsigLimit = 17.85;
136 +      (*it)->PLExpNsigLimit = 12.558;
137 +      (*it)->FCExpNsigLimit = 15.75;
138 +      (*it)->MCMCExpNsigLimit = 15.8085;
139 +      (*it)->ObsNsigLimit = 10.0637;
140 +      (*it)->PLObsNsigLimit = 4.22124;
141 +      (*it)->FCObsNsigLimit  = 0.75;
142 +      (*it)->MCMCObsNsigLimit  = 8.75403;
143 +    }
144 +    else if (flag=="ABCD_HT") {
145 +      (*it)->ExpNsigLimit = 24.0438;
146 +      (*it)->PLExpNsigLimit = 17.8537;
147 +      (*it)->FCExpNsigLimit = 22.05;
148 +      (*it)->MCMCExpNsigLimit = 20.781;
149 +      (*it)->ObsNsigLimit    = 18.24;
150 +      (*it)->PLObsNsigLimit = 9.70375;
151 +      (*it)->FCObsNsigLimit  = 11.85;
152 +      (*it)->MCMCObsNsigLimit  = 14.34;
153 +    }
154 +    
155 +    
156 +      (*it)->ExpXsecLimit    = (*it)->ExpNsigLimit * (*it)->Xsection / (*it)->signal;
157 +      (*it)->PLExpXsecLimit  = (*it)->PLExpNsigLimit  * (*it)->Xsection / (*it)->signal;
158 +      (*it)->FCExpXsecLimit  = (*it)->FCExpNsigLimit  * (*it)->Xsection / (*it)->signal;
159 +      (*it)->MCMCExpXsecLimit= (*it)->MCMCExpNsigLimit  * (*it)->Xsection / (*it)->signal;
160 +      (*it)->ObsXsecLimit    = (*it)->ObsNsigLimit  * (*it)->Xsection / (*it)->signal;
161 +      (*it)->PLObsXsecLimit  = (*it)->PLObsNsigLimit  * (*it)->Xsection / (*it)->signal;
162 +      (*it)->FCObsXsecLimit  = (*it)->FCObsNsigLimit   * (*it)->Xsection / (*it)->signal;
163 +      (*it)->MCMCObsXsecLimit= (*it)->MCMCObsNsigLimit  * (*it)->Xsection / (*it)->signal;
164 +    
165 +  
166 +  }
167 +
168 +
169 + }
170 +

Diff Legend

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