kDSemiMuonic decay option added
authormorsch <andreas.morsch@cern.ch>
Wed, 1 Oct 2014 05:38:08 +0000 (07:38 +0200)
committermorsch <andreas.morsch@cern.ch>
Wed, 1 Oct 2014 05:38:08 +0000 (07:38 +0200)
EVGEN/AliDecayer.h
EVGEN/AliGenMC.cxx
PYTHIA6/AliDecayerPythia.cxx
PYTHIA8/AliDecayerPythia8.cxx

index 0144d82..b9eaefb 100644 (file)
@@ -13,7 +13,7 @@ typedef TVirtualMCDecayer AliDecayer;
 #if ROOT_VERSION_CODE >= 197633  //Corresponds to Root v3-04-01
 typedef enum 
 {
-    kBSemiElectronic, kSemiElectronic, kDiElectron, kBSemiMuonic, kSemiMuonic, kDiMuon, kJpsiDiMuon,
+  kBSemiElectronic, kSemiElectronic, kDiElectron, kBSemiMuonic, kDSemiMuonic, kSemiMuonic, kDiMuon, kJpsiDiMuon,
     kBJpsiDiMuon, kBJpsiDiElectron, 
     kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu, 
     kNoDecay, kHadronicD, kHadronicDWithout4Bodies, kOmega, kLambda, kPhiKK, 
index 1157369..c7f1b95 100644 (file)
@@ -129,6 +129,7 @@ void AliGenMC::Init()
        break;
     case kHardMuons:   
     case kBSemiMuonic:
+    case kDSemiMuonic:
     case kSemiMuonic:
     case kDiMuon:
     case kJpsiDiMuon:
index 7119714..9a8a528 100644 (file)
@@ -291,6 +291,15 @@ void AliDecayerPythia::ForceDecay()
        ForceParticleDecay( 5232,13,1); // Sigma_b    
        ForceParticleDecay( 5332,13,1); // Omega_b    
        break;
+    case kDSemiMuonic:
+       ForceParticleDecay(  411,13,1); // D0     
+       ForceParticleDecay(  421,13,1); // D+/-     
+       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     
index 7fb2265..097c7d3 100644 (file)
@@ -226,6 +226,29 @@ void AliDecayerPythia8::ForceDecay()
        fPythia8->ReadString("5332:onMode = off");
        fPythia8->ReadString("5332:onIfAny = 13");
        break;
+    case kDSemiMuonic:
+//      D+- -> mu X
+       fPythia8->ReadString("411:onMode = off");
+       fPythia8->ReadString("411:onIfAll = 13");
+//      D0  -> mu X
+       fPythia8->ReadString("421:onMode = off");
+       fPythia8->ReadString("421:onIfAll = 13");
+//      D_s  -> mu X
+       fPythia8->ReadString("431:onMode = off");
+       fPythia8->ReadString("431:onIfAll = 13");
+//      Lambda_c -> mu X
+       fPythia8->ReadString("4122:onMode = off");
+       fPythia8->ReadString("4122:onIfAll = 13");
+//      Sigma_c  -> mu X
+       fPythia8->ReadString("4132:onMode = off");
+       fPythia8->ReadString("4132:onIfAll = 13");
+//      Sigma  -> mu X
+       fPythia8->ReadString("4232:onMode = off");
+       fPythia8->ReadString("4232:onIfAll = 13");
+//      Omega_c  -> mu X
+       fPythia8->ReadString("4332:onMode = off");
+       fPythia8->ReadString("4332:onIfAll = 13");
+       break;
     case kSemiMuonic:
 //      D+- -> mu X
        fPythia8->ReadString("411:onMode = off");