Add f0f1 generator

alisw · Jan 20, 2023 · 009e79a · 009e79a
1 parent 71d8d77
commit 009e79a
Showing 1 changed file with 36 additions and 0 deletions.
diff --git a/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C b/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C
@@ -0,0 +1,36 @@
+struct particle_inj {
+  int n;
+  char name[32];
+  int pdg;
+  double maxpt;
+  double maxy;
+};
+
+AliGenerator* GeneratorCustom()
+{
+  AliGenCocktail *ctl = (AliGenCocktail*) GeneratorCocktail("Monash2013_Rsn_f2f0");
+  // pythia8
+  AliGenerator   *py8 = GeneratorPythia8(kPythia8Tune_Monash2013);
+  ctl->AddGenerator(py8, "Pythia8 (Monash2013)", 1.);
+
+  // randomly injected particles
+  const int nParticles = 2;
+  particle_inj particleList[2] = { // {name,pdgcode,maxpt,maxy},
+    {1,"f0(1710)",10331,10.,0.9},
+    {1,"f2(1525)",335,10.,0.9},
+    {1,"f1(1285)",20223,10.,0.9},
+    {1,"f1(1420)",20333,10.,0.9},
+    {1,"f1(1420)",20333,10.,0.9},
+    {1,"f0(1500)",90302,10.,0.9}
+  };
+
+  AliDecayerPythia *dec = new AliDecayerPythia;
+  AliPDG::AddParticlesToPdgDataBase();
+
+  ctl->UseSingleInjectionPerEvent();
+  for (int idx = 0; idx < nParticles; ++idx) {
+    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.);
+  }
+  return ctl;
+}