Option to switch pi0 decay off
authormorsch <andreas.morsch@cern.ch>
Fri, 28 Feb 2014 09:13:39 +0000 (10:13 +0100)
committermorsch <andreas.morsch@cern.ch>
Fri, 28 Feb 2014 09:13:39 +0000 (10:13 +0100)
PYTHIA6/AliDecayerPythia.cxx
PYTHIA6/AliDecayerPythia.h

index c9dae46..6b9a023 100644 (file)
@@ -61,7 +61,8 @@ AliDecayerPythia::AliDecayerPythia():
     fDecay(kAll),
     fHeavyFlavour(kTRUE),
     fLongLived(kFALSE),
-    fPatchOmegaDalitz(0)
+    fPatchOmegaDalitz(0),
+    fPi0(1)
 {
 // Default Constructor
     for (Int_t i=0; i< 501; i++) fBraPart[i]= 1.;
@@ -74,7 +75,8 @@ AliDecayerPythia::AliDecayerPythia(const AliDecayerPythia &decayer):
     fDecay(kAll),
     fHeavyFlavour(kTRUE),
     fLongLived(kFALSE),
-    fPatchOmegaDalitz(0)
+    fPatchOmegaDalitz(0),
+    fPi0(1)
 {
     // Copy Constructor
     decayer.Copy(*this);
@@ -118,8 +120,8 @@ void AliDecayerPythia::Init()
        fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
     }
 */
+    if (fPi0) fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
 
-    fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
     Int_t isw = 0;
     if (fLongLived) isw = 1;
     
index 1026c92..5545a54 100644 (file)
@@ -34,6 +34,7 @@ public AliDecayer
     virtual Float_t GetPartialBranchingRatio(Int_t ipart);
     virtual Float_t GetLifetime(Int_t kf);
     virtual void    SwitchOffBDecay();
+    virtual void    SwitchOffPi0() {fPi0 = 0;}
     virtual void    WriteDecayTable();
     virtual void    ReadDecayTable();
     
@@ -63,9 +64,10 @@ public AliDecayer
     Bool_t      fHeavyFlavour;    //! Flag for heavy flavors
     Bool_t      fLongLived;       //! Flag for long lived particle decay
     Bool_t      fPatchOmegaDalitz;//! Flag to patch the omega Dalitz decays 
+    Bool_t      fPi0;             //! Flag for pi0 decay 
     static Bool_t fgInit;         //! initialization flag 
     
-    ClassDef(AliDecayerPythia, 3) // AliDecayer implementation using Pythia  
+    ClassDef(AliDecayerPythia, 4) // AliDecayer implementation using Pythia  
 };
 #endif