SusyScan/PlotScript/GeneratorMasses.cc


#include "GeneratorMasses.h"

#include <fstream>
#include <iostream>
#include <cmath>

std::vector<GeneratorMasses> FillGeneratorMasses(std::string file)
{
   std::ifstream masses_file;
   masses_file.open(file.c_str());
   std::vector<GeneratorMasses> points;
   while (1) {
      GeneratorMasses p;
      masses_file >> p.MZERO      
                  >> p.MHALF      
                  >> p.TANB       
                  >> p.SGNMU      
                  >> p.AZERO      
                  >> p.MTOP       
                  >> p.muQ        
                  >> p.Q          
                  >> p.M1         
                  >> p.M2
                  >> p.M3         
                  >> p.MGL        
                  >> p.MUL        
                  >> p.MB1        
                  >> p.MSN        
                  >> p.MNTAU      
                  >> p.MZ1        
                  >> p.MW1        
                  >> p.MHL        
                  >> p.MUR        
                  >> p.MB2        
                  >> p.MEL        
                  >> p.MTAU1      
                  >> p.MZ2        
                  >> p.MW2        
                  >> p.MHH        
                  >> p.MDL        
                  >> p.MT1        
                  >> p.MER        
                  >> p.MTAU2      
                  >> p.MZ3        
                  >> p.MHA        
                  >> p.MDR        
                  >> p.MT2        
                  >> p.MZ4        
                  >> p.MHp;

      if (!masses_file.good()) break;
      if (fabs(p.SGNMU)!=1.) {
         std::cerr << "check lines near m0=" << p.MZERO << ", m1/2=" << p.MHALF  << std::endl;
         break;
      } 
      points.push_back( p );
   }
   return points;
}
Revision:	1.1
Committed:	Sat Nov 27 10:41:09 2010 UTC (14 years, 5 months ago) by auterman
Content type:	text/plain
Branch:	MAIN
Log Message:	code cleanup
#	User	Rev	Content
1	auterman	1.1
2			#include "GeneratorMasses.h"
3
4			#include <fstream>
5			#include <iostream>
6			#include <cmath>
7
8			std::vector<GeneratorMasses> FillGeneratorMasses(std::string file)
9			{
10			std::ifstream masses_file;
11			masses_file.open(file.c_str());
12			std::vector<GeneratorMasses> points;
13			while (1) {
14			GeneratorMasses p;
15			masses_file >> p.MZERO
16			>> p.MHALF
17			>> p.TANB
18			>> p.SGNMU
19			>> p.AZERO
20			>> p.MTOP
21			>> p.muQ
22			>> p.Q
23			>> p.M1
24			>> p.M2
25			>> p.M3
26			>> p.MGL
27			>> p.MUL
28			>> p.MB1
29			>> p.MSN
30			>> p.MNTAU
31			>> p.MZ1
32			>> p.MW1
33			>> p.MHL
34			>> p.MUR
35			>> p.MB2
36			>> p.MEL
37			>> p.MTAU1
38			>> p.MZ2
39			>> p.MW2
40			>> p.MHH
41			>> p.MDL
42			>> p.MT1
43			>> p.MER
44			>> p.MTAU2
45			>> p.MZ3
46			>> p.MHA
47			>> p.MDR
48			>> p.MT2
49			>> p.MZ4
50			>> p.MHp;
51
52			if (!masses_file.good()) break;
53			if (fabs(p.SGNMU)!=1.) {
54			std::cerr << "check lines near m0=" << p.MZERO << ", m1/2=" << p.MHALF << std::endl;
55			break;
56			}
57			points.push_back( p );
58			}
59			return points;
60			}