Counters for diffractive collisions
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 29 Sep 2012 07:16:40 +0000 (07:16 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sat, 29 Sep 2012 07:16:40 +0000 (07:16 +0000)
TDPMjet/AliGenDPMjet.cxx
TDPMjet/DPMcommon.h

index 7e8c954..3fd2c40 100644 (file)
@@ -19,6 +19,7 @@
 // Uses the TDPMjet implementation of TGenerator.
 
 #include <TDPMjet.h>
+#include "DPMcommon.h"
 #include <TRandom.h>
 #include <TArrayI.h>
 #include <TParticle.h>
@@ -481,6 +482,7 @@ void AliGenDPMjet::MakeHeader()
     // Event Vertex
     header->SetPrimaryVertex(fVertex);
     header->SetInteractionTime(fTime);
+    ((AliGenDPMjetEventHeader*) header)->SetNDiffractive(POEVT1.nsd1, POEVT1.nsd2, POEVT1.ndd);
     gAlice->SetGenEventHeader(header);    
     AddHeader(header);
 }
index d5bd816..aaff398 100644 (file)
@@ -241,7 +241,7 @@ COMMON_BLOCK_DEF(PoprcsCommon,POPRCS);
 
     COMMON /POEVT1/ NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
      &                JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),
-     &                VHEP(4,NMXHEP)
+     &                VHEP(4,NMXHEP), NSD1, NSD2, NDD
 */
 const Int_t nmxhep = 4000;
 
@@ -254,6 +254,9 @@ typedef struct {
   Int_t    jdahep[nmxhep][2];
   Double_t phep  [nmxhep][5];
   Double_t vhep  [nmxhep][4];
+  Int_t    nsd1;
+  Int_t    nsd2;
+  Int_t    ndd;
 } Poevt1Common;
 
 #define POEVT1 COMMON_BLOCK(POEVT1,poevt1)