-// Make clean
-// AllowAllDecays();
-//
-// select mode
-
- switch (decay)
- {
- case semimuonic:
- if (fProcess==charm || fProcess == charm_unforced) {
- ForceParticleDecay( 411,13,1); // D+/-
- ForceParticleDecay( 421,13,1); // D0
- ForceParticleDecay( 431,13,1); // D_s
- ForceParticleDecay( 4122,13,1); // Lambda_c
- }
- if (fProcess==beauty || fProcess == beauty_unforced) {
- ForceParticleDecay( 511,13,1); // B0
- ForceParticleDecay( 521,13,1); // B+/-
- ForceParticleDecay( 531,13,1); // B_s
- ForceParticleDecay( 5122,13,1); // Lambda_b
- }
- break;
- case dimuon:
- ForceParticleDecay( 41,13,2); // phi
- ForceParticleDecay( 443,13,2); // J/Psi
- ForceParticleDecay(30443,13,2); // Psi'
- ForceParticleDecay( 553,13,2); // Upsilon
- ForceParticleDecay(30553,13,2); // Upsilon'
- break;
- case semielectronic:
-
- ForceParticleDecay( 411,11,1); // D+/-
- ForceParticleDecay( 421,11,1); // D0
- ForceParticleDecay( 431,11,1); // D_s
- ForceParticleDecay( 4122,11,1); // Lambda_c
-
- ForceParticleDecay( 511,11,1); // B0
- ForceParticleDecay( 521,11,1); // B+/-
- ForceParticleDecay( 531,11,1); // B_s
- ForceParticleDecay( 5122,11,1); // Lambda_b
- break;
- case dielectron:
-
- ForceParticleDecay( 41,11,2); // phi
- ForceParticleDecay( 443,11,2); // J/Psi
- ForceParticleDecay(30443,11,2); // Psi'
- ForceParticleDecay( 553,11,2); // Upsilon
- ForceParticleDecay(30553,11,2); // Upsilon'
- break;
- case b_jpsi_dimuon:
- ForceParticleDecay( 511,443,1); // B0
- ForceParticleDecay( 521,443,1); // B+/-
- ForceParticleDecay( 531,443,1); // B_s
- ForceParticleDecay( 5122,443,1); // Lambda_b
- ForceParticleDecay( 443,13,2); // J/Psi
- break;
- case b_psip_dimuon:
- ForceParticleDecay( 511,30443,1); // B0
- ForceParticleDecay( 521,30443,1); // B+/-
- ForceParticleDecay( 531,30443,1); // B_s
- ForceParticleDecay( 5122,30443,1); // Lambda_b
- ForceParticleDecay(30443,13,2); // Psi'
- break;
- case b_jpsi_dielectron:
- ForceParticleDecay( 511,443,1); // B0
- ForceParticleDecay( 521,443,1); // B+/-
- ForceParticleDecay( 531,443,1); // B_s
- ForceParticleDecay( 5122,443,1); // Lambda_b
- ForceParticleDecay( 443,11,2); // J/Psi
- break;
- case b_psip_dielectron:
- ForceParticleDecay( 511,30443,1); // B0
- ForceParticleDecay( 521,30443,1); // B+/-
- ForceParticleDecay( 531,30443,1); // B_s
- ForceParticleDecay( 5122,30443,1); // Lambda_b
- ForceParticleDecay(30443,11,2); // Psi'
- break;
- case pitomu:
- ForceParticleDecay(211,13,1); // pi->mu
- break;
- case katomu:
- ForceParticleDecay(321,13,1); // K->mu
- break;
- }