#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,
break;
case kHardMuons:
case kBSemiMuonic:
+ case kDSemiMuonic:
case kSemiMuonic:
case kDiMuon:
case kJpsiDiMuon:
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
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");