Typo corrected.

[u/mrichter/AliRoot.git] / EVGEN / AliDecayerPythia.cxx

diff --git a/EVGEN/AliDecayerPythia.cxx b/EVGEN/AliDecayerPythia.cxx
index e2c060d00be8431a625866ee33bfb65f494559dd..c91f0493e3e77eed3e323fa53ef2f4d2ad8f320b 100644 (file)
--- a/EVGEN/AliDecayerPythia.cxx
+++ b/EVGEN/AliDecayerPythia.cxx
@@ -15,6 +15,12 @@
 
 /*
 $Log$
+Revision 1.11  2002/04/26 10:32:59  morsch
+Option kNoDecayHeavy added.
+
+Revision 1.10  2002/02/22 17:28:05  morsch
+ReadDecayTable() and WriteDecayTable() methods added.
+
 Revision 1.9  2001/12/20 10:37:13  morsch
 - Add omega forced decay.
 - Semileptonic decays for some more B and D baryons.
@@ -104,11 +110,15 @@ void AliDecayerPythia::Init()
     fPythia->ResetDecayTable();
     for (j=0; j < 14; j++) {
        kc=fPythia->Pycomp(heavy[j]);
-       fPythia->SetMDCY(kc,1,1);
-       for (i=fPythia->GetMDCY(kc,2); 
-            i<fPythia->GetMDCY(kc,2)+fPythia->GetMDCY(kc,3); 
-            i++) {
+       if (fDecay == kNoDecayHeavy) {
+         fPythia->SetMDCY(kc,1,0);
+       } else {
+         fPythia->SetMDCY(kc,1,1);
+         for (i=fPythia->GetMDCY(kc,2); 
+              i<fPythia->GetMDCY(kc,2)+fPythia->GetMDCY(kc,3); 
+              i++) {
            fPythia->SetMDME(i,1,1);
+         }
        }
     }
 
@@ -141,7 +151,9 @@ void AliDecayerPythia::ForceDecay()
 {
 // Force a particle decay mode
     Decay_t decay=fDecay;
-    
+    fPythia->SetMSTJ(21,2);
+    if (decay == kNoDecayHeavy) return;
+
 //
 // select mode    
 
@@ -235,11 +247,17 @@ void AliDecayerPythia::ForceDecay()
     case kHadronicD:
        ForceHadronicD();
        break;
+    case kPhiKK:
+       ForceParticleDecay(333,321,2); // Phi->K+K-
+       break;
     case kOmega:
        ForceOmega();
     case kAll:
        break;
     case kNoDecay:
+       fPythia->SetMSTJ(21,0);
+       break;
+    case kNoDecayHeavy:
        break;
     }
 }