31 {
32
33
37 {
38 std::cout << "No track data file is available! " << std::endl;
39 abort();
40 }
42
43 m_inputFile.open( m_inputFileName.c_str() );
44 if ( !m_inputFile )
45 {
46 cout << "EvtTrackGen: PROBLEMS OPENING FILE " << m_inputFileName << endl;
47 exit( 0 );
48 }
49
50 Evt.clear();
51 while ( !m_inputFile.eof() )
52 {
53 m_inputFile >> nParticles;
54 std::vector<EvtVector4R> vp4;
55 for ( int i = 0; i < nParticles; i++ )
56 {
57 EvtVector4R p4;
58 m_inputFile >> idParticles[i];
59 m_inputFile >> pxParticle;
60 p4.
set( 1, pxParticle );
61 m_inputFile >> pyParticle;
62 p4.
set( 2, pyParticle );
63 m_inputFile >> pzParticle;
64 p4.
set( 3, pzParticle );
65 m_inputFile >> eParticle;
66 p4.
set( 0, eParticle );
67 vp4.push_back( p4 );
68 }
69 Evt.push_back( vp4 );
70 }
71 Evt.pop_back();
72
74 {
75 std::cout << "The number of daughters are not cosistent with that the data file"
76 << std::endl;
77 abort();
78 }
79 for ( int i = 0; i < nParticles; i++ )
80 {
83 {
84 std::cout << "The daughter particle pdg in your data file is not consistent with you "
85 "decay card."
86 << std::endl;
87 abort();
88 }
89 }
90
91
92
93
94
95
96
97
98
99
100
101}
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
static std::vector< std::string > SV
static EvtId evtIdFromStdHep(int stdhep)
void set(int i, double d)