Channel kHardMuons added.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 8 Oct 2003 09:28:15 +0000 (09:28 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 8 Oct 2003 09:28:15 +0000 (09:28 +0000)
EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx

index 9fc1deb..0e80c5b 100644 (file)
@@ -17,7 +17,7 @@ typedef enum
     kBJpsiDiMuon, kBJpsiDiElectron, 
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu,
     kKaToMu, kNoDecay,
-    kHadronicD, kOmega, kPhiKK, kAll, kNoDecayHeavy
+    kHadronicD, kOmega, kPhiKK, kAll, kNoDecayHeavy, kHardMuons
 } Decay_t;
 #endif
 
index 04ebde2..662a0e4 100644 (file)
@@ -92,6 +92,7 @@ void AliGenMC::Init()
     case kBPsiPrimeDiElectron:
        fChildSelect[0] = kElectron;    
        break;
+    case kHardMuons:   
     case kSemiMuonic:
     case kDiMuon:
     case kBJpsiDiMuon:
index 2372c66..4c4f32a 100644 (file)
@@ -121,9 +121,36 @@ void AliDecayerPythia::ForceDecay()
 // select mode    
     Int_t products[2];
     Int_t mult[2];
-
+    Int_t products1[3];
+    Int_t mult1[3];
+    
     switch (decay) 
     {
+    case kHardMuons:
+       products1[0] =     13;
+       products1[1] =    443;
+       products1[2] = 100443;
+       mult1[0] = 1;
+       mult1[1] = 1;
+       mult1[2] = 1;
+       ForceParticleDecay(  511, products1, mult1, 3); 
+       ForceParticleDecay(  521, products1, mult1, 3); 
+       ForceParticleDecay(  531, products1, mult1, 3); 
+       ForceParticleDecay( 5122, products1, mult1, 3); 
+       ForceParticleDecay( 5132, products1, mult1, 3); 
+       ForceParticleDecay( 5232, products1, mult1, 3); 
+       ForceParticleDecay( 5332, products1, mult1, 3); 
+       ForceParticleDecay( 100443, 443, 1);  // Psi'  -> J/Psi X    
+       ForceParticleDecay(    443,  13, 2);  // J/Psi -> mu+ mu-   
+
+       ForceParticleDecay(  411,13,1); // D+/-     
+       ForceParticleDecay(  421,13,1); // D0     
+       ForceParticleDecay(  431,13,1); // D_s     
+       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 kSemiMuonic:
        ForceParticleDecay(  411,13,1); // D+/-     
        ForceParticleDecay(  421,13,1); // D0     
@@ -209,7 +236,7 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(  521,30443,1); // B+/-     
        ForceParticleDecay(  531,30443,1); // B_s     
        ForceParticleDecay( 5122,30443,1); // Lambda_b 
-       ForceParticleDecay(100443,11,2);    // Psi'   
+       ForceParticleDecay(100443,11,2);   // Psi'   
        break;
     case kPiToMu:
        ForceParticleDecay(211,13,1); // pi->mu