Decayer options for W added. (Zaida Conesa)

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

diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx
index 1a7b16c34be35947e046a43832fb9b7ed7aef513..6a7d662f5498c4837518e84fbd5a9663eebe8550 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.cxx
+++ b/PYTHIA6/AliDecayerPythia.cxx
@@ -223,10 +223,10 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(    443,13,2);             // J/Psi -> mu+ mu-   
        break;
     case kBPsiPrimeDiMuon:
-       ForceParticleDecay(  511,30443,1); // B0     
-       ForceParticleDecay(  521,30443,1); // B+/-     
-       ForceParticleDecay(  531,30443,1); // B_s     
-       ForceParticleDecay( 5122,30443,1); // Lambda_b 
+       ForceParticleDecay(  511,100443,1); // B0     
+       ForceParticleDecay(  521,100443,1); // B+/-     
+       ForceParticleDecay(  531,100443,1); // B_s     
+       ForceParticleDecay( 5122,100443,1); // Lambda_b 
        ForceParticleDecay(100443,13,2);    // Psi'   
        break;
     case kBJpsiDiElectron:
@@ -236,11 +236,17 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay( 5122,443,1); // Lambda_b
        ForceParticleDecay(  443,11,2);  // J/Psi    
        break;
+    case kBJpsi:
+       ForceParticleDecay(  511,443,1); // B0     
+       ForceParticleDecay(  521,443,1); // B+/-     
+       ForceParticleDecay(  531,443,1); // B_s     
+       ForceParticleDecay( 5122,443,1); // Lambda_b
+       break;
     case kBPsiPrimeDiElectron:
-       ForceParticleDecay(  511,30443,1); // B0     
-       ForceParticleDecay(  521,30443,1); // B+/-     
-       ForceParticleDecay(  531,30443,1); // B_s     
-       ForceParticleDecay( 5122,30443,1); // Lambda_b 
+       ForceParticleDecay(  511,100443,1); // B0     
+       ForceParticleDecay(  521,100443,1); // B+/-     
+       ForceParticleDecay(  531,100443,1); // B_s     
+       ForceParticleDecay( 5122,100443,1); // Lambda_b 
        ForceParticleDecay(100443,11,2);   // Psi'   
        break;
     case kPiToMu:
@@ -248,7 +254,23 @@ void AliDecayerPythia::ForceDecay()
        break;
     case kKaToMu:
        ForceParticleDecay(321,13,1); // K->mu     
+       break;     
+    case kWToMuon:
+        ForceParticleDecay(  24, 13,1); // W -> mu
+       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;
@@ -514,7 +536,7 @@ void AliDecayerPythia::Streamer(TBuffer &R__b)
 }
 #endif
 
-void AliDecayerPythia::Copy(AliDecayerPythia &) const
+void AliDecayerPythia::Copy(TObject &) const
 {
   //
   // Copy *this onto AliDecayerPythia -- not implemented
@@ -1494,6 +1516,7 @@ void AliDecayerPythia::Copy(AliDecayerPythia &) const
            797    0   42    0.001000    e+              nu_e            pi0                                             
            798    0   42    0.001000    e+              nu_e            eta                                             
            799    0   42    0.001000    e+              nu_e            eta'                                            
+
            800    0   42    0.001000    e+              nu_e            rho0                                            
            801    0   42    0.001000    e+              nu_e            omega                                           
            802    1   42    0.070000    mu+             nu_mu           Kbar0                                           
@@ -3674,3 +3697,4 @@ void AliDecayerPythia::Copy(AliDecayerPythia &) const
           2553    1    0    0.610139    W+              e-                                                              
 
 */
+