Decayer options for W added. (Zaida Conesa)

diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx
index 94ba15c33a6b96ccf17747d08d4b665f77bcdabd..6a7d662f5498c4837518e84fbd5a9663eebe8550 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.cxx
+++ b/PYTHIA6/AliDecayerPythia.cxx
@@ -170,8 +170,7 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay( 5122,13,1); // Lambda_b    
        ForceParticleDecay( 5132,13,1); // Xsi_b    
        ForceParticleDecay( 5232,13,1); // Sigma_b    
        ForceParticleDecay( 5122,13,1); // Lambda_b    
        ForceParticleDecay( 5132,13,1); // Xsi_b    
        ForceParticleDecay( 5232,13,1); // Sigma_b    

-       ForceParticleDecay( 5332,13,1); // Omega_b  
-       ForceParticleDecay(   24,13,1); // W+/- -> mu+/-      
+       ForceParticleDecay( 5332,13,1); // Omega_b    

        break;
     case kDiMuon:
        ForceParticleDecay(  221,13,2); // eta
        break;
     case kDiMuon:
        ForceParticleDecay(  221,13,2); // eta


@@ -255,7 +254,23 @@ void AliDecayerPythia::ForceDecay()
        break;
     case kKaToMu:
        ForceParticleDecay(321,13,1); // K->mu     
        break;
     case kKaToMu:
        ForceParticleDecay(321,13,1); // K->mu     

+       break;     
+    case kWToMuon:
+        ForceParticleDecay(  24, 13,1); // W -> mu

        break;
        break;

+    case kWToCharm:
+        ForceParticleDecay(   24, 4,1); // W -> c
+       break;
+    case kWToCharmToMuon:
+         ForceParticleDecay(   24, 4,1); // W -> c
+        ForceParticleDecay(  411,13,1); // D+/- -> mu
+        ForceParticleDecay(  421,13,1); // D0  -> mu
+        ForceParticleDecay(  431,13,1); // D_s  -> mu
+        ForceParticleDecay( 4122,13,1); // Lambda_c
+        ForceParticleDecay( 4132,13,1); // Xsi_c
+        ForceParticleDecay( 4232,13,1); // Sigma_c
+        ForceParticleDecay( 4332,13,1); // Omega_c
+        break;

     case kHadronicD:
        ForceHadronicD();
        break;
     case kHadronicD:
        ForceHadronicD();
        break;