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

Comparing UserCode/auterman/SusyScan/PlotScript/scan.C (file contents):
Revision 1.1 by auterman, Fri Nov 30 13:10:50 2007 UTC vs.
Revision 1.2 by auterman, Tue Dec 4 12:57:34 2007 UTC

# Line 176 | Line 176 | void TScan::ReadGeneratedMasses(std::str
176     }
177   }
178  
179 + void TScan::ReadGeneratedXsects(std::string file)
180 + {
181 +   ifstream xsects_file;
182 +   xsects_file.open(file.c_str());
183 +   SUSY_XSECS x;
184 +   while (1) {
185 +      xsects_file >> x.MZERO      
186 +                  >> x.MHALF      
187 +                  >> x.TANB      
188 +                  >> x.SGNMU      
189 +                  >> x.AZERO      
190 +                  >> x.MTOP      
191 +                  >> x.XS0
192 +                  >> x.XS201
193 +                  >> x.XS202
194 +                  >> x.XS204
195 +                  >> x.XS205
196 +                  >> x.XS207
197 +                  >> x.XS208
198 +                  >> x.XS209
199 +                  >> x.XS210
200 +                  >> x.XS211
201 +                  >> x.XS212
202 +                  >> x.XS213
203 +                  >> x.XS214
204 +                  >> x.XS216
205 +                  >> x.XS217
206 +                  >> x.XS218
207 +                  >> x.XS219
208 +                  >> x.XS220
209 +                  >> x.XS221
210 +                  >> x.XS222
211 +                  >> x.XS223
212 +                  >> x.XS224
213 +                  >> x.XS225
214 +                  >> x.XS226
215 +                  >> x.XS227
216 +                  >> x.XS228
217 +                  >> x.XS229
218 +                  >> x.XS230
219 +                  >> x.XS231
220 +                  >> x.XS232
221 +                  >> x.XS233
222 +                  >> x.XS234
223 +                  >> x.XS235
224 +                  >> x.XS236
225 +                  >> x.XS237
226 +                  >> x.XS238
227 +                  >> x.XS239
228 +                  >> x.XS240
229 +                  >> x.XS241
230 +                  >> x.XS242
231 +                  >> x.XS243
232 +                  >> x.XS244
233 +                  >> x.XS246
234 +                  >> x.XS247
235 +                  >> x.XS248
236 +                  >> x.XS249
237 +                  >> x.XS250
238 +                  >> x.XS251
239 +                  >> x.XS252
240 +                  >> x.XS253
241 +                  >> x.XS254
242 +                  >> x.XS256
243 +                  >> x.XS258
244 +                  >> x.XS259
245 +                  >> x.XS261
246 +                  >> x.XS262
247 +                  >> x.XS263
248 +                  >> x.XS264
249 +                  >> x.XS265
250 +                  >> x.XS271
251 +                  >> x.XS272
252 +                  >> x.XS273
253 +                  >> x.XS274
254 +                  >> x.XS275
255 +                  >> x.XS276
256 +                  >> x.XS277
257 +                  >> x.XS278
258 +                  >> x.XS279
259 +                  >> x.XS280
260 +                  >> x.XS281
261 +                  >> x.XS282
262 +                  >> x.XS283
263 +                  >> x.XS284
264 +                  >> x.XS285
265 +                  >> x.XS286
266 +                  >> x.XS287
267 +                  >> x.XS288
268 +                  >> x.XS289
269 +                  >> x.XS290
270 +                  >> x.XS291
271 +                  >> x.XS292
272 +                  >> x.XS293
273 +                  >> x.XS294
274 +                  >> x.XS295
275 +                  >> x.XS296
276 +                  >> x.XS297
277 +                  >> x.XS298
278 +                  >> x.XS299
279 +                  >> x.XS300
280 +                  >> x.XS301;
281 +      if (!xsects_file.good()) break;
282 +      if (fabs(x.SGNMU)!=1.) {
283 +         cerr << "check lines near m0=" << x.MZERO << ", m1/2=" << x.MHALF << endl;
284 +         break;
285 +      }
286 +      xsects.push_back(x);
287 +   }
288 + }
289  
290   int TScan::DoStuff()
291   {
# Line 183 | Line 293 | int TScan::DoStuff()
293     ReadGeneratedMasses( "mass_scan.dat" );
294     cout << "Read in " << points.size() << " mass points." <<endl;
295  
296 +   ///Read in the generated x-sections
297 +   ReadGeneratedXsects( "xsec_scan.dat" );
298 +   cout << "Read in " << xsects.size() << " x-section points." <<endl;
299 +
300     ///Where are no mass-points generated?
301     TH2F   * no_solution         = Area51(  );
302     MyTGraph * contour_no_solution = GetContour( no_solution );
# Line 249 | Line 363 | int TScan::DoStuff()
363  
364  
365  
252 –   //c2->Draw();             //interactive
253 –   //gApplication->Run();    //interactive
254 –
366     c1->SaveAs("dirt.eps");
367 +   gApplication->Run();    //interactive
368 +
369     //ps.NewPage();
370  
371   }

Diff Legend

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