kNoDecay, kHadronicD, kHadronicDWithout4Bodies, kOmega, kLambda, kPhiKK,
kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kZDiElectron, kNeutralPion, kAllMuonic,
- kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty
+ kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty, kPsiPrimeJpsiDiElectron
} Decay_t;
#endif
fChildSelect[0]= kProton;
fChildSelect[1]= 211;
break;
+ case kPsiPrimeJpsiDiElectron:
+ fChildSelect[0]= 211;
+ fChildSelect[1]= 11;
+ break;
case kOmega:
case kAll:
case kAllMuonic:
ForceParticleDecay(100553,11,2);// Upsilon'
ForceParticleDecay(200553,11,2);// Upsilon''
break;
+ case kPsiPrimeJpsiDiElectron:
+ products[0] = 443;
+ products[1] = 211;
+ mult[0] = 1;
+ mult[1] = 2;
+ ForceParticleDecay( 100443, products, mult, 2, 1);
+ ForceParticleDecay( 443,11,2);
+ break;
case kBJpsiDiMuon:
products[0] = 443;
if (norm > 0.) fBraPart[kc] /= norm;
}
-void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart)
+void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart, Bool_t flag)
{
//
// Force decay of particle into products with multiplicity mult
for (Int_t i = 0; i < npart; i++) {
nprod += (CountProducts(channel, products[i]) >= mult[i]);
}
- if (nprod) {
+ if ((nprod && !flag) || ((nprod == npart) && flag)) {
fPythia->SetMDME(channel,1,1);
- } else {
+ } else { //
fPythia->SetMDME(channel,1,0);
fBraPart[kc] -= fPythia->GetBRAT(channel);
}
Double_t mom, Double_t theta, Double_t phi);
Int_t CountProducts(Int_t channel, Int_t particle);
void ForceParticleDecay(Int_t particle, Int_t product, Int_t mult);
- void ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart);
+ void ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart, Bool_t flag = 0);
void ForceHadronicD(Int_t optUser4Bodies=1);
void ForceOmega();
void ForceLambda();
case kNoDecayHeavy:
case kNoDecayBeauty:
case kNeutralPion:
+ case kPsiPrimeJpsiDiElectron:
break;
}
}