PhotosHEPEVTEvent.h

Go to the documentation of this file.

#ifndef _PhotosHEPEVTEvent_h_included_
#define _PhotosHEPEVTEvent_h_included_

#include <iostream>
#include "PhotosEvent.h"
#include "PhotosParticle.h"
#include "PhotosHEPEVTParticle.h"

namespace Photospp {

// 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 PhotosHEPEVTParticle;

class PhotosHEPEVTEvent : public PhotosEvent {

 public:
  ~PhotosHEPEVTEvent();

  PhotosHEPEVTEvent();

  void addParticle(PhotosHEPEVTParticle *p);

  PhotosHEPEVTParticle *getParticle(int i);

  void setParticle(int i, PhotosHEPEVTParticle *p);

  int getParticleCount();

  std::vector<PhotosParticle *> getParticleList();

  void print();

  void clear();

#ifdef USE_HEPEVT_INTERFACE

  static void read_event_from_HEPEVT(PhotosHEPEVTEvent *evt);

  static void write_event_to_HEPEVT(PhotosHEPEVTEvent *evt);
#endif

 private:
  std::vector<PhotosHEPEVTParticle *> particle_list;
};

} // namespace Photospp
#endif