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