40 {
41
42
46 std::cout<<"No track data file is available! "<<std::endl;abort();
48
49 m_inputFile.open(m_inputFileName.c_str());
50 if (!m_inputFile){
51 cout << "EvtTrackGen: PROBLEMS OPENING FILE "
52 << m_inputFileName << endl;
53 exit(0);
54 }
55
56 Evt.clear();
57 while(!m_inputFile.eof()){
58 m_inputFile >> nParticles;
59 std::vector<EvtVector4R> vp4;
60 for (int i = 0; i < nParticles; i++){
62 m_inputFile >> idParticles[i];
63 m_inputFile >> pxParticle;
p4.set(1,pxParticle);
64 m_inputFile >> pyParticle;
p4.set(2,pyParticle);
65 m_inputFile >> pzParticle;
p4.set(3,pzParticle);
66 m_inputFile >> eParticle;
p4.set(0,eParticle);
68 }
69 Evt.push_back(vp4);
70 }
71 Evt.pop_back();
72
73 if(nParticles!=
getNDaug()){std::cout<<
"The number of daughters are not cosistent with that the data file"<<std::endl;abort();}
74 for(int i=0;i<nParticles;i++){
76 if(pid!=
getDaug(i) ){std::cout<<
"The daughter particle pdg in your data file is not consistent with you decay card."<<std::endl;abort();}
77 }
78
79
80
81
82
83
84
85
86
87
88}
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
static std::vector< std::string > SV
static EvtId evtIdFromStdHep(int stdhep)
double double double * p4