Access to POEVT1 common block.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 9 Mar 2010 08:40:13 +0000 (08:40 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 9 Mar 2010 08:40:13 +0000 (08:40 +0000)
TDPMjet/DPMcommon.h
TDPMjet/TDPMjet.cxx
TDPMjet/TDPMjet.h

index 3abee3a..d5bd816 100644 (file)
@@ -237,5 +237,26 @@ COMMON_BLOCK_DEF(PoprcsCommon,POPRCS);
 /**********************************************************/
 /*          D E S C R I P T I O N :                      */
 /*--------------------------------------------------------*/
+/*
+
+    COMMON /POEVT1/ NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
+     &                JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),
+     &                VHEP(4,NMXHEP)
+*/
+const Int_t nmxhep = 4000;
+
+typedef struct {
+  Int_t    nevhep;
+  Int_t    nhep;
+  Int_t    isthep[nmxhep];
+  Int_t    idhep [nmxhep];
+  Int_t    jmohep[nmxhep][2];
+  Int_t    jdahep[nmxhep][2];
+  Double_t phep  [nmxhep][5];
+  Double_t vhep  [nmxhep][4];
+} Poevt1Common;
+
+#define POEVT1 COMMON_BLOCK(POEVT1,poevt1)
+COMMON_BLOCK_DEF(Poevt1Common,POEVT1);
 
 #endif
index 2fef2bb..4667239 100644 (file)
@@ -258,6 +258,7 @@ Int_t TDPMjet::ImportParticles(TClonesArray *particles, Option_t *option)
       } // Particle loop  
   }
   return nump;
+  
 }
 
 
@@ -558,3 +559,9 @@ void TDPMjet::Dt_Rndmou(int u, int c, int cd, int cm, int i, int j)
        dt_rndmou(u, c, cd, cm, i, j);
 }
 
+
+Int_t TDPMjet::NHEP()                   const {return POEVT1.nhep;}
+Int_t TDPMjet::ISTHEP(Int_t i)          const {return POEVT1.isthep[i];}
+Int_t TDPMjet::IDHEP(Int_t i)           const {return POEVT1.idhep[i];}
+Int_t TDPMjet::PHEP(Int_t i, Int_t j)   const {return POEVT1.phep[i][j];}
+   
index d50c1fe..1454672 100644 (file)
@@ -99,6 +99,15 @@ public:
    virtual Int_t GetTargSpectators() const;
    virtual Int_t    GetProcessCode() const;
 
+   // Phojet access (POEVT1)
+   virtual Int_t NHEP()                   const;
+   virtual Int_t ISTHEP(Int_t i)          const;
+   virtual Int_t IDHEP(Int_t i)           const;
+   virtual Int_t PHEP(Int_t i, Int_t j)   const;
+   
+      
+          
+   
    // Access to DPMJET routines:
    virtual void Dt_Dtuini(int nevts, double epn, int npmass, int npchar, 
                          int ntmass, int ntchar, int idp, int iemu);
@@ -115,6 +124,7 @@ public:
    virtual void Dt_Rndmin(int u, int c, int cd, int cm, int i, int j);   
    virtual void Dt_Rndmou(int u, int c, int cd, int cm, int i, int j);   
 
+   
 protected:
 
    Int_t        fNEvent;  // Event number to be generated