kBJpsiDiMuon, kBJpsiDiElectron,
kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu,
kKaToMu, kNoDecay,
- kHadronicD, kOmega, kPhiKK, kAll, kNoDecayHeavy, kHardMuons
+ kHadronicD, kOmega, kPhiKK, kAll, kNoDecayHeavy, kHardMuons, kBJpsi
} Decay_t;
#endif
Int_t mult[2];
Int_t products1[3];
Int_t mult1[3];
+ printf("Forcedecay %d %d \n", decay, kAll);
switch (decay)
{
ForceParticleDecay( 5122,443,1); // Lambda_b
ForceParticleDecay( 443,11,2); // J/Psi
break;
+ case kBJpsi:
+ ForceParticleDecay( 511,443,1); // B0
+ ForceParticleDecay( 521,443,1); // B+/-
+ ForceParticleDecay( 531,443,1); // B_s
+ ForceParticleDecay( 5122,443,1); // Lambda_b
+ break;
case kBPsiPrimeDiElectron:
ForceParticleDecay( 511,30443,1); // B0
ForceParticleDecay( 521,30443,1); // B+/-