From 009e79a13b86fe5f40d17921c1614efd6187dd07 Mon Sep 17 00:00:00 2001
From: Bong-Hwi Lim <bong-hwi.lim@cern.ch>
Date: Fri, 20 Jan 2023 15:12:37 +0100
Subject: [PATCH] Add f0f1 generator

---
 .../PWGLF/Pythia8_Monash2013_Rsn_f0f1.C       | 36 +++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100644 MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C

diff --git a/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C b/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C
new file mode 100644
index 00000000..e5909b1f
--- /dev/null
+++ 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;
+}