Process type added to header. (A. Bogdanov)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 3 Nov 2006 17:04:03 +0000 (17:04 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 3 Nov 2006 17:04:03 +0000 (17:04 +0000)
TDPMjet/AliGenDPMjet.cxx
TDPMjet/AliGenDPMjet.h
TDPMjet/DPMcommon.h
TDPMjet/TDPMjet.cxx
TDPMjet/TDPMjet.h

index d528ab2..9f6af3a 100644 (file)
 #include <TParticleClassPDG.h>
 #include <TPDGCode.h>
 #include <TLorentzVector.h>
-
+#include "AliRunLoader.h"
 #include "AliGenDPMjet.h"
 #include "AliGenDPMjetEventHeader.h"
-#include "AliPythia.h"
 #include "AliRun.h"
 #include "AliDpmJetRndm.h"
+#include "AliHeader.h"
+#include "AliStack.h"
+#include "AliMC.h"
 
 ClassImp(AliGenDPMjet)
 
@@ -93,8 +95,6 @@ AliGenDPMjet::AliGenDPMjet(Int_t npart)
     SetTarget();
     SetProjectile();
     fVertex.Set(3);
-    // Instance AliPythia
-    AliPythia::Instance(); 
     AliDpmJetRndm::SetDpmJetRandom(GetRandom());
 }
 
@@ -402,14 +402,24 @@ void AliGenDPMjet::MakeHeader()
     ((AliGenDPMjetEventHeader*) header)->SetTotalEnergy(fDPMjet->GetTotEnergy());
     ((AliGenDPMjetEventHeader*) header)->SetParticipants(fDPMjet->GetfIp(), 
                                                         fDPMjet->GetfIt());
-
+ ((AliGenDPMjetEventHeader*) header)->SetProcessType(fDPMjet->GetProcessCode());
 // Bookkeeping for kinematic bias
     ((AliGenDPMjetEventHeader*) header)->SetTrials(fTrials);
 // Event Vertex
     header->SetPrimaryVertex(fVertex);
     gAlice->SetGenEventHeader(header);    
+ AddHeader(header);
 }
 
+void AliGenDPMjet::AddHeader(AliGenEventHeader* header)
+{
+    // Add header to container or runloader
+    if (fContainer) {
+        fContainer->AddHeader(header);
+    } else {
+        AliRunLoader::GetRunLoader()->GetHeader()->SetGenEventHeader(header);
+    }
+}
 
 
 //______________________________________________________________________________
index 74482a0..fcbad1b 100644 (file)
@@ -16,6 +16,12 @@ class TArrayI;
 class TParticle;
 class TClonesArray;
 class TGraph;
+class AliGenEventHeader;
+class AliStack;
+class AliRunLoader;
+class AliGenDPMjetEventHeader;
+
+
 
 class AliGenDPMjet : public AliGenMC
 {
@@ -28,7 +34,7 @@ class AliGenDPMjet : public AliGenMC
     virtual void    Generate();
     virtual void    Init();
     virtual void    SetEnergyCMS(Float_t energy = 14000.) {fEnergyCMS = energy; fBeamEn = energy / 2.;}
-    virtual void    SetImpactParameterRange(Float_t bmin=0., Float_t bmax=16.)
+    virtual void    SetImpactParameterRange(Float_t bmin=0., Float_t bmax=1.)
                        {fMinImpactParam=bmin; fMaxImpactParam=bmax;}
     virtual void    SetProcess(DpmProcess_t iproc) {fProcess = iproc;}
     virtual void    SetCentral(Int_t icentr=-2) {fICentr = icentr;}
@@ -52,6 +58,7 @@ class AliGenDPMjet : public AliGenMC
     virtual Float_t GetGenImpPar() {return fGenImpPar;}
     
     AliGenDPMjet &  operator=(const AliGenDPMjet & rhs);
+    void     AddHeader(AliGenEventHeader* header);
 
  protected:
     Bool_t SelectFlavor(Int_t pid);
index effd213..3abee3a 100644 (file)
@@ -217,5 +217,25 @@ COMMON_BLOCK_DEF(DtglcpCommon,DTGLCP);
 /* NWTAAC = total number of interacting nucleons in proj. */
 /* NWTBAC = total number of interacting nucleons in target*/
 /*========================================================*/
+/*========================================================*/
+/* COMMON /POPRCS/IPROCE,IDNODF,IDIFR1,IDIFR2,IDDPOM,     */
+/*               IPRON(15,4)                                     */
+/*--------------------------------------------------------*/
+typedef struct {
+   Int_t       iproce;  
+   Int_t        idnodf;
+   Int_t        idifr1;
+   Int_t        idifr2;
+   Int_t        iddpom;
+   Int_t        ipron[15][4];
+} PoprcsCommon;
+
+#define POPRCS COMMON_BLOCK(POPRCS,poprcs)
+COMMON_BLOCK_DEF(PoprcsCommon,POPRCS);
+
+
+/**********************************************************/
+/*          D E S C R I P T I O N :                      */
+/*--------------------------------------------------------*/
 
 #endif
index 5fc4942..c456288 100644 (file)
@@ -465,7 +465,11 @@ Int_t TDPMjet::GetTargSpectators() const
 {
        return DTGLCP.nwtbac;
 }
-
+//______________________________________________________________________________
+Int_t TDPMjet::GetProcessCode() const
+{
+               return POPRCS.iproce;
+}
 //______________________________________________________________________________
 void TDPMjet::Dt_Rndm(int idummy)
 {
index fbb8d43..f5f8f1d 100644 (file)
@@ -96,8 +96,8 @@ public:
    virtual Int_t GetTargWounded() const;
    virtual Int_t GetProjSpectators() const;
    virtual Int_t GetTargSpectators() const;
-   
-  
+   virtual Int_t    GetProcessCode() 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);