kZDiMuon added. (Zaida Conesa)

diff --git a/EVGEN/AliDecayer.h b/EVGEN/AliDecayer.h
index 5cb7cf3..cab3b1c 100644 (file)
--- a/EVGEN/AliDecayer.h
+++ b/EVGEN/AliDecayer.h
@@ -18,7 +18,7 @@ typedef enum
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu, 
     kNoDecay, kHadronicD, kOmega, kPhiKK, 
     kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
-    kWToMuon,kWToCharm, kWToCharmToMuon
+    kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon
 } Decay_t;
 #endif
 

diff --git a/EVGEN/AliGenMC.cxx b/EVGEN/AliGenMC.cxx
index bdbb1d2..a8e4e36 100644 (file)
--- a/EVGEN/AliGenMC.cxx
+++ b/EVGEN/AliGenMC.cxx
@@ -107,6 +107,7 @@ void AliGenMC::Init()
     case kKaToMu:
     case kWToMuon:
     case kWToCharmToMuon:
+    case kZDiMuon:
        fChildSelect[0]=kMuonMinus;
        break;
     case kWToCharm:

diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx
index c31924f..36a24b1 100644 (file)
--- a/PYTHIA6/AliDecayerPythia.cxx
+++ b/PYTHIA6/AliDecayerPythia.cxx
@@ -256,7 +256,7 @@ void AliDecayerPythia::ForceDecay()
        break;
     case kKaToMu:
        ForceParticleDecay(321,13,1); // K->mu     
-       break;     
+       break;
     case kWToMuon:
         ForceParticleDecay(  24, 13,1); // W -> mu
        break;
@@ -264,15 +264,18 @@ void AliDecayerPythia::ForceDecay()
         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;
+        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 kZDiMuon:
+        ForceParticleDecay(  23, 13,2); // Z -> mu+ mu-
+       break;
     case kHadronicD:
        ForceHadronicD();
        break;