#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"
{
	/** Definition of the PH_HEPEVT common block */
	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_;

	/** Fortran debug routine */
	extern "C" void phodmp_();
}

class PH_HEPEVT_Interface
{
public:
	/** Convert PhotosBranch to HEPEVT */
	static int  set(PhotosBranch* branch);

	/** Update event record with data from HEPEVT */
	static void get();

	/** Prepare particles for processing */
	static void prepare();

	/** Check channel for complete matrix element calculation */
	static void check_ME_channel();

	/** Finalize processing */
	static void complete();

	/** Clear HEPEVT */
	static void clear();
public:
	/** Index of decaying particle*/
	static int decay_idx;
	/** Number of channel to be used - flag for fortran routines */
	static int ME_channel;
private:
	/** Add single particle to HEPEVT */
	static void add_particle(int i, PhotosParticle * particle,
	                         int first_mother, int last_mother,
	                         int first_daughter, int last_daughter);

	/** List of particles added to HEPEVT */
	static std::vector<PhotosParticle*> m_particle_list;
};

} // namespace Photospp
#endif