From 77a1cc123cd1754989df8d6d62126878c56342a6 Mon Sep 17 00:00:00 2001 From: Bong-Hwi Lim Date: Thu, 29 Jun 2023 18:23:57 +0200 Subject: [PATCH] Use decayer --- .../PWGLF/Pythia8_Monash2013_Rsn_f0f1.C | 41 +++++-------------- 1 file changed, 10 insertions(+), 31 deletions(-) diff --git a/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C b/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C index 693c3a01..2804ebda 100644 --- a/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C +++ b/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C @@ -13,7 +13,6 @@ AliGenerator *GeneratorCustom() // pythia8 AliGenerator *py8 = GeneratorPythia8(kPythia8Tune_Monash2013); ctl->AddGenerator(py8, "Pythia8 (Monash2013)", 1.); - // randomly injected particles const int nParticles = 7; particle_inj particleList[nParticles] = { @@ -26,37 +25,11 @@ AliGenerator *GeneratorCustom() {1, "f1(1285)", 20223, 20., 0.9}, {1, "f1(1420)", 20333, 20., 0.9}}; - // Force decay mode of f0(1500) to K0s K0s - disabled since it's configured in the decay table directly - // (AliPythia8::Instance())->ReadString("9030221:onMode = off"); - // (AliPythia8::Instance())->ReadString("9030221:oneChannel = 1 1.00 0 310 -310"); - // (AliPythia8::Instance())->ReadString("9030221:onIfMatch = 310 -310"); - // Force decay mode of f0(1710) to K0s K0s - (AliPythia8::Instance())->ReadString("10331:onMode = off"); - (AliPythia8::Instance())->ReadString("10331:oneChannel = 1 1.00 0 310 -310"); - (AliPythia8::Instance())->ReadString("10331:onIfMatch = 310 -310"); - // force decay of f2(1525) to K0s K0s - (AliPythia8::Instance())->ReadString("335:onMode = off"); - (AliPythia8::Instance())->ReadString("335:oneChannel = 1 1.00 0 310 -310"); - (AliPythia8::Instance())->ReadString("335:onIfMatch = 310 -310"); - // force decay of f2(1270) to K0s K0s - (AliPythia8::Instance())->ReadString("225:onMode = off"); - (AliPythia8::Instance())->ReadString("225:oneChannel = 1 1.00 0 310 -310"); - (AliPythia8::Instance())->ReadString("225:onIfMatch = 310 -310"); - // force decay of f0(1370) to K0s K0s - (AliPythia8::Instance())->ReadString("10221:onMode = off"); - (AliPythia8::Instance())->ReadString("10221:oneChannel = 1 1.00 0 310 -310"); - (AliPythia8::Instance())->ReadString("10221:onIfMatch = 310 -310"); - // force decay of f1(1285) to K0s, anti-K, pi+ - (AliPythia8::Instance())->ReadString("20223:onMode = off"); - (AliPythia8::Instance())->ReadString("20223:oneChannel = 1 1.00 0 310 -321 211"); - (AliPythia8::Instance())->ReadString("20223:onIfMatch = 310 -321 211"); - // force decay of f1(1420) to K0s, anti-K, pi+ - (AliPythia8::Instance())->ReadString("20333:onMode = off"); - (AliPythia8::Instance())->ReadString("20333:oneChannel = 1 1.00 0 310 -321 211"); - (AliPythia8::Instance())->ReadString("20333:onIfMatch = 310 -321 211"); - AliDecayerPythia *dec = new AliDecayerPythia; - AliPDG::AddParticlesToPdgDataBase(); + dec->SetForceDecay(kF0F1K0s); + AliDecayerPythia *dec2 = new AliDecayerPythia; // for f1 to K0s K anti-pi + dec2->SetForceDecay(kF1K0sK); + AliPDG::AddParticlesToPdgDataBase(); ctl->UseSingleInjectionPerEvent(); for (int idx = 0; idx < nParticles; ++idx) @@ -64,5 +37,11 @@ AliGenerator *GeneratorCustom() AliGenerator *inj = GeneratorParam(particleList[idx].n, particleList[idx].pdg, 1., particleList[idx].maxpt, -particleList[idx].maxy, particleList[idx].maxy, dec); ctl->AddGenerator(inj, (TString(particleList[idx].name) + " injector").Data(), 1.); } + // f1(1285) to K0s K anti-pi + AliGenerator *inj_20223 = GeneratorParam(1, 20223, 1., 20., -0.9, 0.9, dec2); + ctl->AddGenerator(inj_20223, "f1(1285) to K0,K,pi- injector", 1.); + // f1(1420) to K0s K anti-pi + AliGenerator *inj_20333 = GeneratorParam(1, 20333, 1., 20., -0.9, 0.9, dec2); + ctl->AddGenerator(inj_20333, "f1(1420) to K0,K,pi- injector", 1.); return ctl; }