kBSemiElectronic and kBSemiMuonic added.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 6 Feb 2007 08:37:49 +0000 (08:37 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 6 Feb 2007 08:37:49 +0000 (08:37 +0000)
EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx

index cab3b1c..a5126ce 100644 (file)
@@ -13,7 +13,7 @@ typedef TVirtualMCDecayer AliDecayer;
 #if ROOT_VERSION_CODE >= 197633  //Corresponds to Root v3-04-01
 typedef enum 
 {
-    kSemiElectronic, kDiElectron, kSemiMuonic, kDiMuon,
+    kBSemiElectronic, kSemiElectronic, kDiElectron, kBSemiMuonic, kSemiMuonic, kDiMuon,
     kBJpsiDiMuon, kBJpsiDiElectron, 
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu, 
     kNoDecay, kHadronicD, kOmega, kPhiKK, 
index 677b8aa..84f642f 100644 (file)
@@ -111,6 +111,7 @@ void AliGenMC::Init()
 //
 //  Initialization
     switch (fForceDecay) {
+    case kBSemiElectronic:
     case kSemiElectronic:
     case kDiElectron:
     case kBJpsiDiElectron:
@@ -118,6 +119,7 @@ void AliGenMC::Init()
        fChildSelect[0] = kElectron;    
        break;
     case kHardMuons:   
+    case kBSemiMuonic:
     case kSemiMuonic:
     case kDiMuon:
     case kBJpsiDiMuon:
index 2093bd4..75bde8a 100644 (file)
@@ -176,6 +176,15 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay( 4232,13,1); // Sigma_c 
        ForceParticleDecay( 4332,13,1); // Omega_c     
        break;
+    case kBSemiMuonic:
+       ForceParticleDecay(  511,13,1); // B0     
+       ForceParticleDecay(  521,13,1); // B+/-     
+       ForceParticleDecay(  531,13,1); // B_s     
+       ForceParticleDecay( 5122,13,1); // Lambda_b    
+       ForceParticleDecay( 5132,13,1); // Xsi_b    
+       ForceParticleDecay( 5232,13,1); // Sigma_b    
+       ForceParticleDecay( 5332,13,1); // Omega_b    
+       break;
     case kSemiMuonic:
        ForceParticleDecay(  411,13,1); // D+/-     
        ForceParticleDecay(  421,13,1); // D0     
@@ -203,6 +212,15 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay(100553,13,2);// Upsilon'
        ForceParticleDecay(200553,13,2);// Upsilon''
        break;
+    case kBSemiElectronic:
+       ForceParticleDecay(  511,11,1); // B0     
+       ForceParticleDecay(  521,11,1); // B+/-     
+       ForceParticleDecay(  531,11,1); // B_s     
+       ForceParticleDecay( 5122,11,1); // Lambda_b     
+       ForceParticleDecay( 5132,11,1); // Xsi_b    
+       ForceParticleDecay( 5232,11,1); // Sigma_b    
+       ForceParticleDecay( 5332,11,1); // Omega_b    
+       break;
     case kSemiElectronic:
        ForceParticleDecay(  411,11,1); // D+/-     
        ForceParticleDecay(  421,11,1); // D0