---

HepEvt.h

Go to the documentation of this file.

// -*- C++ -*-  $Id: HepEvt.h,v 1.1 2000/01/17 23:37:29 burnett Exp $
//

// define the HEPEVT commom
#ifndef __HEPEVT_H
#define __HEPEVT_H

#include "gismo/Generator.h"
#include "CLHEP/Utilities/fortran.h"
typedef long integer;
typedef float real;

unsigned const nmxhep = 50;
extern "C"{
//BEGIN_FORTRAN_COMMON(hepevt)
extern struct HEPEVT{
 integer nevhep;                // event number
 integer        nhep;                  // number of entries
 integer isthep[nmxhep];        // status code
 integer idhep[nmxhep];         // PDG particle id
 integer jmohep[nmxhep][2];     // position of first, second mother
 integer jdahep[nmxhep][2];     // position of first, last daughter
 real      phep[nmxhep][5];       // 4-momemtum, mass (GeV)
 real    vhep[nmxhep][4];       // vertex, production time in mm
END_FORTRAN_COMMON(hepevt)
}

class HepEvt  : public Generator
{
public:
        static long& nevhep(){return FORTRAN_COMMON(hepevt).nevhep;}
        static long& nhep(){ return FORTRAN_COMMON(hepevt).nhep;}
        static long& isthep(int i){return FORTRAN_COMMON(hepevt).isthep[i-1];}
        static long& idhep(int i){return FORTRAN_COMMON(hepevt).idhep[i-1];}
        static long* jmohep(int i){return FORTRAN_COMMON(hepevt).jmohep[i-1];}
        static long* jdahep(int i){return FORTRAN_COMMON(hepevt).jdahep[i-1];}

        static float* phep(int i){return FORTRAN_COMMON(hepevt).phep[i-1];}
        static float* vhep(int i){return FORTRAN_COMMON(hepevt).vhep[i-1];}
//      static float& mass(int i){return FORTRAN_COMMON(hepevt).phep[i-1][4];}

        HepEvt(std::istream& _in);
        // constructor that reads from an open stream

        MCParticle * getEvent();
        // creates structure with next event

        int readEvent();
        // fill the common from input file. failed if n()==0

        void propagate(const Medium*);
        // special propagate that should find the unpropagated daughters

private:
        std::istream& inputFile;
        bool _hasVertices;
};



#endif

---