Possibility to leave pi0 decays enabled during primary event generation.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Feb 2007 10:37:12 +0000 (10:37 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Feb 2007 10:37:12 +0000 (10:37 +0000)
EVGEN/AliDecayer.h
EVGEN/AliGenHIJINGpara.cxx
PYTHIA6/AliDecayerPythia.cxx

index a5126ce..b8da294 100644 (file)
@@ -18,7 +18,7 @@ typedef enum
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu, 
     kNoDecay, kHadronicD, kOmega, kPhiKK, 
     kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
-    kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon
+    kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kNeutralPion
 } Decay_t;
 #endif
 
index 061b2a6..20cabde 100644 (file)
@@ -281,6 +281,13 @@ void AliGenHIJINGpara::Init()
 //
     if (fPi0Decays && gMC)
        fDecayer = gMC->GetDecayer();
+
+    if (fPi0Decays)
+    {
+       fDecayer->SetForceDecay(kNeutralPion);
+       fDecayer->Init();
+    }
+    
 }
 
 
index 75bde8a..6712c61 100644 (file)
@@ -100,7 +100,14 @@ void AliDecayerPythia::Init()
        }
     }
 //...Switch off decay of pi0, K0S, Lambda, Sigma+-, Xi0-, Omega-.
-    fPythia->SetMDCY(fPythia->Pycomp(111) ,1,0);
+
+    if (fDecay != kNeutralPion) {
+       fPythia->SetMDCY(fPythia->Pycomp(111) ,1,0);
+    } else {
+       fPythia->SetMDCY(fPythia->Pycomp(111) ,1,1);
+    }
+    
+     
     fPythia->SetMDCY(fPythia->Pycomp(310) ,1,0);
     fPythia->SetMDCY(fPythia->Pycomp(3122),1,0);
     fPythia->SetMDCY(fPythia->Pycomp(3112),1,0);
@@ -328,7 +335,9 @@ void AliDecayerPythia::ForceDecay()
        fPythia->SetMSTJ(21,0);
        break;
     case kNoDecayHeavy:
+    case kNeutralPion:
        break;
+
     }
 }