TauolaHEPEVTEvent.h

#ifndef _TauolaHEPEVTEvent_h_included_
#define _TauolaHEPEVTEvent_h_included_

#include <iostream>
#include "TauolaEvent.h"
#include "TauolaParticle.h"
#include "TauolaHEPEVTParticle.h"

namespace Tauolapp
{

// Uncomment this line to use interface to common block HEPEVT
// But first be sure about suitable for you value of NMXHEP
// and whether phep, vhep should be declared float or double
//#define USE_HEPEVT_INTERFACE

#ifdef USE_HEPEVT_INTERFACE

// Change this value to match HEPEVT size
const int NMXHEP = 10000;

extern "C" struct {
  int    nevhep;            // serial number
  int    nhep;              // number of particles
  int    isthep[NMXHEP];    // status code
  int    idhep [NMXHEP];    // particle PDG ID
  int    jmohep[NMXHEP][2]; // parent particles
  int    jdahep[NMXHEP][2]; // childreen particles
  double phep  [NMXHEP][5]; // four-momentum, mass [GeV]
  double vhep  [NMXHEP][4]; // vertex [mm]
} hepevt_;

#endif

class TauolaHEPEVTParticle;

class TauolaHEPEVTEvent : public TauolaEvent {

 public:

  ~TauolaHEPEVTEvent();

  TauolaHEPEVTEvent();

  void addParticle(TauolaHEPEVTParticle *p);

  TauolaHEPEVTParticle *getParticle(int i);

  int getParticleCount();

  std::vector<TauolaParticle*> findParticles(int pdgID);

  std::vector<TauolaParticle*> findStableParticles(int pdgID);

  void print();
  
  void clear();

#ifdef USE_HEPEVT_INTERFACE

  static void read_event_from_HEPEVT(TauolaHEPEVTEvent *evt);
  
  static void write_event_to_HEPEVT(TauolaHEPEVTEvent *evt);
#endif

 private:

  std::vector<TauolaHEPEVTParticle*> particle_list;
};

} // namespace Tauolapp
#endif