PH_HEPEVT_Interface.h

Go to the documentation of this file.

#ifndef _PH_HEPEVT_Interface_included_
#define _PH_HEPEVT_Interface_included_

#include <vector>
#include "PhotosBranch.h"
#include "PhotosParticle.h"

namespace Photospp {

const static int NMXHEP = 10000;
const static double NO_BOOST_THRESHOLD = 1.0e-8;

extern "C" {
extern struct {
  int nevhep;
  int nhep;
  int isthep[NMXHEP];
  int idhep[NMXHEP];
  int jmohep[NMXHEP][2];
  int jdahep[NMXHEP][2];
  double phep[NMXHEP][5];
  double vhep[NMXHEP][4];
  //      COMMON/PH_HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
  //  &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
} ph_hepevt_;

extern struct {
  int qedrad[NMXHEP];
} ph_phoqed_;

extern "C" void phodmp_();
}

class PH_HEPEVT_Interface {
 public:
  static int set(PhotosBranch *branch);

  static void get();

  static void prepare();

  static void check_ME_channel();

  static void complete();

  static void clear();

 public:
  static int decay_idx;
  static int ME_channel;

 private:
  static void add_particle(int i, PhotosParticle *particle, int first_mother, int last_mother, int first_daughter, int last_daughter);

  static std::vector<PhotosParticle *> m_particle_list;
};

} // namespace Photospp
#endif