More user settings for decays.

diff --git a/TDPMjet/AliGenDPMjet.cxx b/TDPMjet/AliGenDPMjet.cxx
index 32624e68615d9b2c7996123fc376eda906eb275d..a315e3d3f39be8366506a47937c6c8f3b3d926f7 100644 (file)
--- a/TDPMjet/AliGenDPMjet.cxx
+++ b/TDPMjet/AliGenDPMjet.cxx
@@ -57,6 +57,7 @@ AliGenDPMjet::AliGenDPMjet()
      fNoGammas(0),
      fLHC(0),
      fPi0Decay(0),
+     fDecayAll(0),
      fGenImpPar(0.),
      fProcess(kDpmMb)
 {
@@ -83,6 +84,7 @@ AliGenDPMjet::AliGenDPMjet(Int_t npart)
      fNoGammas(0),
      fLHC(0),
      fPi0Decay(0),
+     fDecayAll(0),
      fGenImpPar(0.),
      fProcess(kDpmMb)
 {
@@ -113,6 +115,7 @@ AliGenDPMjet::AliGenDPMjet(const AliGenDPMjet &/*Dpmjet*/)
      fNoGammas(0),
      fLHC(0),
      fPi0Decay(0),
+     fDecayAll(0),
      fGenImpPar(0.),
      fProcess(kDpmMb)
 {
@@ -146,6 +149,7 @@ void AliGenDPMjet::Init()
     fDPMjet->SetfFCentr(fICentr);  
     fDPMjet->SetbRange(fMinImpactParam, fMaxImpactParam); 
     fDPMjet->SetPi0Decay(fPi0Decay);
+    fDPMjet->SetDecayAll(fDecayAll);
 //
 //  Initialize DPMjet  
 //    

diff --git a/TDPMjet/AliGenDPMjet.h b/TDPMjet/AliGenDPMjet.h
index 588da78bf6ff41eb5db27a1f90e12448e7c86506..4dc0307c7229d959c55550ff29ff391950fe8e3a 100644 (file)
--- a/TDPMjet/AliGenDPMjet.h
+++ b/TDPMjet/AliGenDPMjet.h
@@ -44,6 +44,7 @@ class AliGenDPMjet : public AliGenMC
     virtual void    SetSpectators(Int_t spects=1)     {fSpectators = spects;}
     virtual void    SetBoostLHC(Int_t flag=0)         {fLHC        = flag;}
     virtual void    SetPi0Decay(Int_t iPi0)  {fPi0Decay = iPi0;}
+    virtual void    SetDecayAll(Int_t iDec)  {fDecayAll = iDec;}
     virtual void    GetImpactParameterRange(Float_t& bmin, Float_t& bmax)
                        {bmin = fMinImpactParam; bmax = fMaxImpactParam;}
     virtual Int_t   GetSpectators()        {return fSpectators;}
@@ -74,6 +75,7 @@ class AliGenDPMjet : public AliGenMC
     Int_t         fNoGammas;       // Don't write gammas if flag "on"
     Int_t         fLHC;            // Assume LHC as lab frame
     Int_t         fPi0Decay;       // Flag for pi0 decays
+    Int_t         fDecayAll;       // Flag to switch on long-lived particle decays
     Float_t      fGenImpPar;      // GeneratedImpactParameter
     DpmProcess_t  fProcess;        // Process type
     

diff --git a/TDPMjet/TDPMjet.cxx b/TDPMjet/TDPMjet.cxx
index aa3524960c9b63f39ddc534f0e1651b8ef929587..55ebdd1fea6acaf1fc9df0576e653e16bc85aa18 100644 (file)
--- a/TDPMjet/TDPMjet.cxx
+++ b/TDPMjet/TDPMjet.cxx
@@ -111,6 +111,7 @@ ClassImp(TDPMjet)
        fBmax(0.),
        fFCentr(0),
        fPi0Decay(0),
+       fDecayAll(0),
        fProcess(kDpmMb)
 {
 // Default Constructor
@@ -133,6 +134,7 @@ TDPMjet::TDPMjet(DpmProcess_t  iproc, Int_t Ip=208, Int_t Ipz=82, Int_t It=208,
       fBmax(0.),
       fFCentr(0),
       fPi0Decay(0),
+      fDecayAll(0),
       fProcess(iproc)
 {  
     printf("TDPMJet Constructor %d %d %d %d \n", Ip, Ipz, It, Itz);
@@ -272,10 +274,13 @@ void TDPMjet::Initialize()
 //  Beam energy and crossing-angle
     fprintf(out, "BEAM      %10.1f%10.1f%10.1f%10.1f%10.1f%10.1f\n",fEpn, fEpn, 0., 0., 0., 0.);
 //  Centrality
-    fprintf(out, "CENTRAL   %10.1f%10.1f%10.1f%10.1f%10.1f%10.1f\n",-1., fBmin, fBmax, 0., 0., 0.);
+    if (fIp > 1. && fIt > 1) 
+       fprintf(out, "CENTRAL   %10.1f%10.1f%10.1f%10.1f%10.1f%10.1f\n",-1., fBmin, fBmax, 0., 0., 0.);
 //  Particle decays
     if (fPi0Decay) 
-    fprintf(out, "PARDECAY  %10.1f%10.1f%10.1f%10.1f%10.1f%10.1f\n", 2., 0., 0., 0., 0., 0.);    
+       fprintf(out, "PARDECAY  %10.1f%10.1f%10.1f%10.1f%10.1f%10.1f\n", 2., 0., 0., 0., 0., 0.);    
+
+    
 //
 //  PHOJET specific
     fprintf(out, "PHOINPUT\n");
@@ -293,6 +298,26 @@ void TDPMjet::Initialize()
         fprintf(out, "PROCESS           0 0 0 0 0 0 1 0\n");
     }
     
+    Int_t iPDG[19] = 
+       {
+//          K0s   pi0  lam   sig+  sig-  tet0
+       310,  111, 3122, 3222, 3112, 3322,
+//          tet- om-    D+      D0     Ds+
+       3312, 3334,  411,  421,  431,
+//          etac lamc+ sigc++ sigc+ sigc0 Ksic+
+        441, 4122, 4222, 4212, 4112, 4232,
+//         Ksic0 sig0 
+       4132, 3212
+       };
+    
+    
+    Int_t iON = (fDecayAll) ? 1:0;
+    for (Int_t i = 0; i < 19; i++) {
+       fprintf(out, "LUND-DECAY%5d %5d\n",  iPDG[i], iON);    
+    }
+       
+    fprintf(out, "LUND-MSTJ %5d %5d\n",   22, 1);    
+
     fprintf(out, "ENDINPUT\n");
 //
 //  START card

diff --git a/TDPMjet/TDPMjet.h b/TDPMjet/TDPMjet.h
index 867280e745e06eb0a99000ef3a973a4768a9d785..209af0ac25cbed64ce4e6cc5c577ccd8133b5f67 100644 (file)
--- a/TDPMjet/TDPMjet.h
+++ b/TDPMjet/TDPMjet.h
@@ -70,8 +70,8 @@ public:
    virtual void  SetfFCentr(Int_t icentr)  {fFCentr = icentr;}
    virtual Int_t GetfFCentr() const        {return fFCentr;}
 
-   virtual void  SetPi0Decay(Int_t iPi0)  {fPi0Decay = iPi0;}
-
+   virtual void  SetPi0Decay(Int_t iPi0)    {fPi0Decay = iPi0;}
+   virtual void  SetDecayAll(Int_t iDecAll) {fDecayAll = iDecAll;}
    
    // Access to DPMJET common blocks:
    virtual Int_t    GetEvNum() const;            
@@ -130,6 +130,7 @@ protected:
    Float_t      fBmax;   // Maximum impact parameter
    Int_t        fFCentr;  // Flag to force central collisions
    Int_t        fPi0Decay;// Flag for pi0 decays
+   Int_t        fDecayAll;// Flag to decay also long-lived particles
    DpmProcess_t fProcess; // Process type
    
    ClassDef(TDPMjet,2)  //Interface to DPMJET Event Generator