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.10 by auterman, Wed Jun 22 15:03:36 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 +   std::cout << "Read " << _scan.size() << " points from argument list." <<std::endl;
16 + }
17 +
18 + void TheLimits::Fill(const std::string filelist)
19 + {
20 +   std::ifstream masses_file;
21 +   masses_file.open(filelist.c_str());
22 +   std::string file;
23 +   while (1) {
24 +      GeneratorMasses * p = new GeneratorMasses;
25 +      masses_file >> file;
26 +      if (!masses_file.good()) break;
27 +      add( new SusyScan(file));
28 +   }
29 +   std::cout << "Read " << _scan.size() << " points from file " << filelist <<std::endl;
30 +   masses_file.close();
31 + }
32  
33   void TheLimits::match()
34   {
# Line 46 | Line 69 | void TheLimits::match()
69          (*it)->MZ4 = (*gt)->MZ4;
70          (*it)->MHp = (*gt)->MHp;
71          match = true;
72 +        //std::cout <<"m0="<<(*it)->Mzero<<", m1/2="<<(*it)->Mhalf<< ", sq="<<(*it)->MUL << ", gl="<<(*it)->MGL <<std::endl;
73        }  
74      }
75      //if (!match) std::cout << "No match for M0="<<(*it)->Mzero
# Line 97 | Line 121 | void TheLimits::FillGeneratorMasses(std:
121                    >> p->MHp;
122  
123        if (!masses_file.good()) break;
124 +      //std::cout << p->Mzero<<", "<<p->Mhalf<<", "<<p->TanBeta<<std::endl;
125        if (fabs(p->Mu)!=1.) {
126           std::cerr << "check lines near m0=" << p->Mzero << ", m1/2=" << p->Mhalf << std::endl;
127           break;
128        }
129 +      //std::cout <<"m0="<< p->Mzero<<", m1/2="<< p->Mhalf<< ", sq="<<p->MUL << ", gl="<<p->MGL <<std::endl;
130        _masses.push_back( p );
131     }
132 +   std::cout << "Read " << _masses.size() << " points from file " << file <<std::endl;
133 +
134 + }
135 +
136 +
137 + void TheLimits::OverwriteLimits(std::string flag)
138 + {
139 +  std::cout<<"WARNING: OVERWRITING LIMITS!!!!"<<std::endl;
140 +  for (std::vector<SusyScan*>::iterator it=_scan.begin();it!=_scan.end();++it){
141 +    if (flag=="ABCD_MHT") {
142 +      (*it)->ExpNsigLimit = 17.85;
143 +      (*it)->PLExpNsigLimit = 12.558;
144 +      (*it)->FCExpNsigLimit = 15.75;
145 +      (*it)->MCMCExpNsigLimit = 15.8085;
146 +      (*it)->ObsNsigLimit = 10.0637;
147 +      (*it)->PLObsNsigLimit = 4.22124;
148 +      (*it)->FCObsNsigLimit  = 0.75;
149 +      (*it)->MCMCObsNsigLimit  = 8.75403;
150 +    }
151 +    else if (flag=="ABCD_HT") {
152 +      (*it)->ExpNsigLimit = 24.0438;
153 +      (*it)->PLExpNsigLimit = 17.8537;
154 +      (*it)->FCExpNsigLimit = 22.05;
155 +      (*it)->MCMCExpNsigLimit = 20.781;
156 +      (*it)->ObsNsigLimit    = 18.24;
157 +      (*it)->PLObsNsigLimit = 9.70375;
158 +      (*it)->FCObsNsigLimit  = 11.85;
159 +      (*it)->MCMCObsNsigLimit  = 14.34;
160 +    }
161 +    
162 +    
163 +      (*it)->ExpXsecLimit    = (*it)->ExpNsigLimit * (*it)->Xsection / (*it)->signal;
164 +      (*it)->PLExpXsecLimit  = (*it)->PLExpNsigLimit  * (*it)->Xsection / (*it)->signal;
165 +      (*it)->FCExpXsecLimit  = (*it)->FCExpNsigLimit  * (*it)->Xsection / (*it)->signal;
166 +      (*it)->MCMCExpXsecLimit= (*it)->MCMCExpNsigLimit  * (*it)->Xsection / (*it)->signal;
167 +      (*it)->ObsXsecLimit    = (*it)->ObsNsigLimit  * (*it)->Xsection / (*it)->signal;
168 +      (*it)->PLObsXsecLimit  = (*it)->PLObsNsigLimit  * (*it)->Xsection / (*it)->signal;
169 +      (*it)->FCObsXsecLimit  = (*it)->FCObsNsigLimit   * (*it)->Xsection / (*it)->signal;
170 +      (*it)->MCMCObsXsecLimit= (*it)->MCMCObsNsigLimit  * (*it)->Xsection / (*it)->signal;
171 +    
172 +  
173 +  }
174 +
175 +
176 + }
177 +
178 + void TheLimits::ExpandGrid(int s)
179 + {
180 +  std::vector<SusyScan*> new_grid;
181 +  for (std::vector<SusyScan*>::iterator it=_scan.begin(); it!=_scan.end(); ++it){
182 +    double x=9999, y=9999;
183 +    std::vector<SusyScan*>::iterator next_x=_scan.end(), next_y=_scan.end();
184 +    for (std::vector<SusyScan*>::iterator nx=_scan.begin(); nx!=_scan.end(); ++nx)
185 +      if ((*it)->Mzero<(*nx)->Mzero && fabs((*it)->Mzero-(*nx)->Mzero)<x) {
186 +        x=fabs((*it)->Mzero-(*nx)->Mzero);
187 +        next_x=nx;
188 +      }
189 +    for (std::vector<SusyScan*>::iterator ny=_scan.begin(); ny!=_scan.end(); ++ny)
190 +      if ( (*it)->Mzero==(*ny)->Mzero && (*it)->Mhalf<(*ny)->Mhalf && fabs((*it)->Mhalf-(*ny)->Mhalf)<y) {
191 +        y=fabs((*it)->Mhalf-(*ny)->Mhalf);
192 +        next_y=ny;
193 +      }
194 +    if (next_y!=_scan.end()) {  
195 +        SusyScan*  n = new SusyScan((**it + **next_y) * 0.5);
196 +        new_grid.push_back( n );
197 +        //std::cout
198 +        //<< "l:"<<(*it)->Mzero<<","<<(*it)->Mhalf
199 +        //<< " <>"<<    n->Mzero<<","<<n->Mhalf
200 +        //<< " <> r:"<<(*next_y)->Mzero<<","<<(*next_y)->Mhalf
201 +        //<<std::endl;
202 +    }  
203 +      
204 +  }
205 +  _scan.insert(_scan.end(), new_grid.begin(), new_grid.end());
206   }

Diff Legend

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