kZDiMuon added. (Zaida Conesa)

diff --git a/EVGEN/AliDecayer.h b/EVGEN/AliDecayer.h
index 5cb7cf313b5afdfa3f00851c74987d6269fc687e..cab3b1c3ac490a68f8d1872cb636e1b6183f7fdf 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 bdbb1d224d5d444cf2b0d176af84993c0beb4f19..a8e4e36e61e0f9083bd3bb14c53292dd3804e8b6 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 c31924f9f00d64af3c414bf84c999c90f8971b0d..36a24b1b08f4c5a1b718982197a7d5247d31756d 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;