6 #ifndef HEPMC3_HEPEVT_WRAPPER_H 7 #define HEPMC3_HEPEVT_WRAPPER_H 24 #ifndef HEPMC3_HEPEVT_NMXHEP 26 #define HEPMC3_HEPEVT_NMXHEP 10000 29 #ifndef HEPMC3_HEPEVT_PRECISION 31 #define HEPMC3_HEPEVT_PRECISION double 71 #ifndef HEPEVT_WRAPPER_HEADER_ONLY 78 extern struct HEPEVT* hepevtptr;
88 static void print_hepevt( std::ostream& ostr = std::cout );
91 static void print_hepevt_particle(
int index, std::ostream& ostr = std::cout );
94 static bool check_hepevt_consistency( std::ostream& ostr = std::cout );
97 static void zero_everything();
98 #ifndef HEPEVT_WRAPPER_HEADER_ONLY 100 static bool GenEvent_to_HEPEVT(
const GenEvent* evt );
102 static bool HEPEVT_to_GenEvent(
GenEvent* evt );
105 static bool fix_daughters();
114 static int status(
const int& index ) {
return hepevtptr->
isthep[index-1]; }
115 static int id(
const int& index ) {
return hepevtptr->
idhep[index-1]; }
120 static double px(
const int& index ) {
return hepevtptr->
phep[index-1][0]; }
121 static double py(
const int& index ) {
return hepevtptr->
phep[index-1][1]; }
122 static double pz(
const int& index ) {
return hepevtptr->
phep[index-1][2]; }
123 static double e(
const int& index ) {
return hepevtptr->
phep[index-1][3]; }
124 static double m(
const int& index ) {
return hepevtptr->
phep[index-1][4]; }
125 static double x(
const int& index ) {
return hepevtptr->
vhep[index-1][0]; }
126 static double y(
const int& index ) {
return hepevtptr->
vhep[index-1][1]; }
127 static double z(
const int& index ) {
return hepevtptr->
vhep[index-1][2]; }
128 static double t(
const int& index ) {
return hepevtptr->
vhep[index-1][3]; }
129 static int number_parents(
const int& index );
130 static int number_children(
const int& index );
131 static int number_children_exact(
const int& index );
135 static void set_status(
const int& index,
const int& status ) { hepevtptr->
isthep[index-1] = status; }
136 static void set_id(
const int& index,
const int&
id ) { hepevtptr->
idhep[index-1] = id; }
137 static void set_parents(
const int& index,
const int& firstparent,
const int& lastparent );
138 static void set_children(
const int& index,
const int& firstchild,
const int& lastchild );
139 static void set_momentum(
const int& index,
const double& px,
const double& py,
const double& pz,
const double& e );
140 static void set_mass(
const int& index,
double mass );
141 static void set_position(
const int& index,
const double& x,
const double& y,
const double& z,
const double& t );
147 inline void HEPEVT_Wrapper::print_hepevt( std::ostream& ostr )
149 ostr <<
" Event No.: " << hepevtptr->
nevhep << endl;
150 ostr<<
" Nr Type Parent(s) Daughter(s) Px Py Pz E Inv. M." << endl;
151 for(
int i=1; i<=hepevtptr->
nhep; ++i )
153 HEPEVT_Wrapper::print_hepevt_particle( i, ostr );
157 inline void HEPEVT_Wrapper::print_hepevt_particle(
int index, std::ostream& ostr )
161 sprintf(buf,
"%5i %6i",index,hepevtptr->
idhep[index-1]);
163 sprintf(buf,
"%4i - %4i ",hepevtptr->
jmohep[index-1][0],hepevtptr->
jmohep[index-1][1]);
165 sprintf(buf,
"%4i - %4i ",hepevtptr->
jdahep[index-1][0],hepevtptr->
jdahep[index-1][1]);
168 sprintf(buf,
"%8.2f %8.2f %8.2f %8.2f %8.2f",hepevtptr->
phep[index-1][0],hepevtptr->
phep[index-1][1],hepevtptr->
phep[index-1][2],hepevtptr->
phep[index-1][3],hepevtptr->
phep[index-1][4]);
173 inline bool HEPEVT_Wrapper::check_hepevt_consistency( std::ostream& )
176 printf(
"HEPEVT_Wrapper::check_hepevt_consistency is not implemented!\n");
180 inline void HEPEVT_Wrapper::zero_everything()
182 memset(hepevtptr,0,
sizeof(
struct HEPEVT));
185 inline int HEPEVT_Wrapper::number_parents(
const int& index )
187 return (hepevtptr->
jmohep[index-1][0]) ? (hepevtptr->
jmohep[index-1][1]) ? hepevtptr->
jmohep[index-1][1]-hepevtptr->
jmohep[index-1][0] : 1 : 0;
190 inline int HEPEVT_Wrapper::number_children(
const int& index )
192 return (hepevtptr->
jdahep[index-1][0]) ? (hepevtptr->
jdahep[index-1][1]) ? hepevtptr->
jdahep[index-1][1]-hepevtptr->
jdahep[index-1][0] : 1 : 0;
195 inline int HEPEVT_Wrapper::number_children_exact(
const int& index )
198 for(
int i=1; i<=hepevtptr->
nhep; ++i )
199 if (((hepevtptr->
jmohep[i-1][0]<=index&&hepevtptr->
jmohep[i-1][1]>=index))||(hepevtptr->
jmohep[i-1][0]==index)||(hepevtptr->
jmohep[i-1][1]==index)) nc++;
205 inline void HEPEVT_Wrapper::set_parents(
const int& index,
const int& firstparent,
const int&lastparent )
207 hepevtptr->
jmohep[index-1][0] = firstparent;
208 hepevtptr->
jmohep[index-1][1] = lastparent;
211 inline void HEPEVT_Wrapper::set_children(
const int& index,
const int& firstchild,
const int& lastchild )
213 hepevtptr->
jdahep[index-1][0] = firstchild;
214 hepevtptr->
jdahep[index-1][1] = lastchild;
217 inline void HEPEVT_Wrapper::set_momentum(
const int& index,
const double& px,
const double& py,
const double& pz,
const double& e )
219 hepevtptr->
phep[index-1][0] = px;
220 hepevtptr->
phep[index-1][1] = py;
221 hepevtptr->
phep[index-1][2] = pz;
222 hepevtptr->
phep[index-1][3] = e;
225 inline void HEPEVT_Wrapper::set_mass(
const int& index,
double mass )
227 hepevtptr->
phep[index-1][4] = mass;
230 inline void HEPEVT_Wrapper::set_position(
const int& index,
const double& x,
const double& y,
const double& z,
const double& t )
232 hepevtptr->
vhep[index-1][0] = x;
233 hepevtptr->
vhep[index-1][1] = y;
234 hepevtptr->
vhep[index-1][2] = z;
235 hepevtptr->
vhep[index-1][3] = t;
240 inline bool HEPEVT_Wrapper::fix_daughters()
248 for(
int i=1; i<=HEPEVT_Wrapper::number_entries(); i++ )
249 for(
int k=1; k<=HEPEVT_Wrapper::number_entries(); k++ )
if (i!=k)
250 if ((HEPEVT_Wrapper::first_parent(k)<=i)&&(i<=HEPEVT_Wrapper::last_parent(k)))
251 HEPEVT_Wrapper::set_children(i,(HEPEVT_Wrapper::first_child(i)==0?k:std::min(HEPEVT_Wrapper::first_child(i),k)),(HEPEVT_Wrapper::last_child(i)==0?k:std::max(HEPEVT_Wrapper::last_child(i),k)));
253 for(
int i=1; i<=HEPEVT_Wrapper::number_entries(); i++ )
254 is_fixed=(is_fixed&&(HEPEVT_Wrapper::number_children_exact(i)==HEPEVT_Wrapper::number_children(i)));
static double m(const int &index)
Get generated mass.
static int last_child(const int &index)
Get index of last daughter.
static double z(const int &index)
Get Z Production vertex.
static double y(const int &index)
Get Y Production vertex.
static double t(const int &index)
Get production time.
static int status(const int &index)
Get status code.
static int event_number()
Get event number.
int nhep
Number of entries in the event.
static double x(const int &index)
Get X Production vertex.
#define HEPMC3_HEPEVT_NMXHEP
static const int NMXHEP
Number of particles in the HEPEVT structure.
static double pz(const int &index)
Get Z momentum.
momentum_t phep[NMXHEP][5]
Momentum: px, py, pz, e, m.
static double e(const int &index)
Get Energy.
momentum_t vhep[NMXHEP][4]
Time-space position: x, y, z, t.
static void set_hepevt_address(char *c)
Set Fortran block address.
static void set_id(const int &index, const int &id)
Set PDG particle id.
HEPMC3_HEPEVT_PRECISION momentum_t
Precision of the 4-momentum, time-space position and mass.
Stores event-related information.
static double py(const int &index)
Get Y momentum.
Fortran common block HEPEVT.
static int first_child(const int &index)
Get index of 1st daughter.
static void set_number_entries(const int &noentries)
Set number of entries.
An interface to HEPEVT common block.
static int id(const int &index)
Get PDG particle id.
static void set_event_number(const int &evtno)
Set event number.
int jmohep[NMXHEP][2]
Pointer to position of 1st and 2nd (or last!) mother.
static int first_parent(const int &index)
Get index of 1st mother.
static double px(const int &index)
Get X momentum.
int isthep[NMXHEP]
Status code.
static void set_status(const int &index, const int &status)
Set status code.
Definition of class GenEvent.
static int last_parent(const int &index)
Get index of last mother.
int jdahep[NMXHEP][2]
Pointer to position of 1nd and 2nd (or last!) daughter.
#define HEPMC3_HEPEVT_PRECISION
static int number_entries()
Get number of entries.
static int max_number_entries()
Block size.