kBJpsiDiMuon, kBJpsiDiElectron,
kBPsiPrimeDiMuon, kBPsiPrimeDiElectron, kPiToMu, kKaToMu,
kNoDecay, kHadronicD, kHadronicDWithout4Bodies, kOmega, kLambda, kPhiKK,
- kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
+ kAll, kNoDecayHeavy, kHardMuons, kBJpsi, kBJpsiUndecayed,
kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kZDiElectron, kNeutralPion, kAllMuonic,
kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty, kPsiPrimeJpsiDiElectron,
kElectronEM, kGammaEM, kDiElectronEM, kBeautyUpgrade
fChildSelect[0]=kKPlus;
break;
case kBJpsi:
+ case kBJpsiUndecayed:
fChildSelect[0]= 443;
break;
case kChiToJpsiGammaToMuonMuon:
ForceParticleDecay( 531,443,1); // B_s
ForceParticleDecay( 5122,443,1); // Lambda_b
break;
+ case kBJpsiUndecayed:
+ ForceParticleDecay( 511,443,1); // B0
+ ForceParticleDecay( 521,443,1); // B+/-
+ ForceParticleDecay( 531,443,1); // B_s
+ ForceParticleDecay( 5122,443,1); // Lambda_b
+ fPythia->SetMDCY(fPythia->Pycomp(443),1,0); // switch-off J/psi
+ break;
case kBPsiPrimeDiElectron:
ForceParticleDecay( 511,100443,1); // B0
ForceParticleDecay( 521,100443,1); // B+/-
case kNeutralPion:
case kBPsiPrimeDiElectron:
case kBeautyUpgrade:
- AliWarning(Form("Warning: case %d not implemented for this class!",(int)decay));
+ case kBJpsiUndecayed:
+ AliWarning(Form("Warning: case %d not implemented for this class!",(int)decay));
break;
}
ReadDecayTable();
//check if particle is already decayed by Pythia
if(part->GetStatusCode() != 1 || part->GetNDaughters()>0)
{
- Info("AliGenEvtGen","Attention: particle %d is already decayed by Pythia!",pdg);
+ AliDebug(1,Form("Attention: particle %d is already decayed by Pythia!",pdg));
continue;
}