EPOS LHC + injected (anti)nuclei (up to He3)

alisw · Dec 17, 2021 · baa8b2a · baa8b2a · alcaliva · Jan 10, 2022
1 parent 1fd6d79
commit baa8b2a
Showing 1 changed file with 54 additions and 0 deletions.
diff --git a/MC/CustomGenerators/PWGLF/EposLHC_Nuclei_Injected.C b/MC/CustomGenerators/PWGLF/EposLHC_Nuclei_Injected.C
@@ -0,0 +1,54 @@
+AliGenerator *GeneratorCustom()  {
+
+    AliGenCocktail *ctl = (AliGenCocktail*) GeneratorCocktail("EposLHC_Nuclei_Injected");
+
+    AliGenerator   *epos  = GeneratorEPOSLHC();
+    ctl->AddGenerator(epos,"EPOSLHC", 1.);
+
+    const Int_t pdgcodes[4] = {
+        2212,
+        1000010020,
+        1000020030,
+        1000010030
+    };
+
+    const Char_t *names[4] = {
+        "Proton",
+        "Deuteron",
+        "Helium-3",
+        "Triton"
+    };
+
+    int npart = 2;
+
+    for (Int_t ipart = 0; ipart < 4; ++ipart) {
+        for (Int_t negative = 0; negative < 2; negative++) {
+
+            AliGenBox *box = new AliGenBox(npart);
+            Int_t pdg = pdgcodes[ipart];
+            if (negative) pdg = -pdg;
+            box->SetPart(pdg);
+            box->SetPtRange(1.0,10.0);
+            box->SetPhiRange(0.0,360.0);
+            box->SetYRange(-0.5,0.5);
+            ctl->AddGenerator(box, Form("%s%i",names[ipart],negative), 1);
+        }
+    }
+
+    return ctl;    
+}
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+