kAll, kNoDecayHeavy, kHardMuons, kBJpsi,
kWToMuon,kWToCharm, kWToCharmToMuon, kZDiMuon, kZDiElectron, kNeutralPion, kAllMuonic,
kChiToJpsiGammaToMuonMuon, kChiToJpsiGammaToElectronElectron, kNoDecayBeauty, kPsiPrimeJpsiDiElectron,
- kElectronEM, kGammaEM, kDiElectronEM, kLambdaBLambda
+ kElectronEM, kGammaEM, kDiElectronEM, kBeautyUpgrade
} Decay_t;
#endif
case kLb: return id=IpLb;
case -kLb: return id=IpLbBar;
case kXi_c: return id=IpXic;
+ case kB: return id=IpBPlus;
+ case -kB: return id=IpBMinus;
default : AliFatal(Form("Unknown particle type: %i",iPID)); id=0;
}
static Int_t IpLb(TRandom * /*ran*/) {return (int)kLb;}
static Int_t IpLbBar(TRandom * /*ran*/) {return -(int)kLb;}
static Int_t IpXic(TRandom * /*ran*/) {return (int)kXi_c;}
+ static Int_t IpBPlus(TRandom * /*ran*/) {return (int)kB;}
+ static Int_t IpBMinus(TRandom * /*ran*/) {return -(int)kB;}
static Double_t PtFlat(const Double_t * /*px*/, const Double_t * /*dummy*/) {return 1;}
static Double_t YFlat (const Double_t * /*py*/, const Double_t * /*dummy*/) {return 1;}
case kNoDecayHeavy:
case kNoDecayBeauty:
case kNeutralPion:
- case kLambdaBLambda:
+ case kBeautyUpgrade:
break;
}
ForceParticleDecay( 331,22,1); // etaprime
ForceParticleDecay( 333,22,1); // phi
break;
- case kLambdaBLambda:
- ForceParticleDecay( 5122, 4122, 1);
- ForceParticleDecay( 4122, 3312, 1);
+ case kBeautyUpgrade:
+ ForceBeautyUpgrade();
break;
}
}
for (Int_t i = 156; i <= 160; i++) fPythia->SetMDME(i, 1, 0);
}
+void AliDecayerPythia::ForceBeautyUpgrade()
+{
+ //
+ // Force dedicated decay channels of signals ineresting
+ // for the ITS upgrade (Lb, Lc, Xi_c, B)
+ //
+
+ ForceParticleDecay( 5122, 4122, 1);
+ ForceParticleDecay( 4232, 3312, 1);
+ const Int_t prod[2]={421,211};
+ Int_t mult[2]={1,1};
+ ForceParticleDecay(521,prod,mult,2,1);
+ ForceHadronicD(1);
+}
+
void AliDecayerPythia::Lu1Ent(Int_t flag, Int_t idpart,
Double_t mom, Double_t theta, Double_t phi)
{
void ForceOmega();
void ForceLambda();
void SwitchOffHeavyFlavour();
+ void ForceBeautyUpgrade();
Float_t GetBraPart(Int_t kf);
void Copy(TObject &decayer) const;
AliTPythia8::Instance()->ReadString("3122:onMode = off");
AliTPythia8::Instance()->ReadString("3122:onIfAll = 2212 211 ");
break;
- case kLambdaBLambda:
+ case kBeautyUpgrade:
AliTPythia8::Instance()->ReadString("5122:onMode = off");
AliTPythia8::Instance()->ReadString("4122:onMode = off");
AliTPythia8::Instance()->ReadString("5122:onIfAll = 4122");