Extended Box_Nuclex_001, added hypertriton, pt_max=10 GeV, npart=10

alisw · Sep 10, 2021 · bd8f1b8 · bd8f1b8
1 parent 7d3c5c7
commit bd8f1b8
Showing 1 changed file with 41 additions and 0 deletions.
diff --git a/MC/CustomGenerators/PWGLF/Box_Nuclex_001_Extended.C b/MC/CustomGenerators/PWGLF/Box_Nuclex_001_Extended.C
@@ -0,0 +1,41 @@
+AliGenerator *GeneratorCustom()  {
+
+    AliGenCocktail *ctl = (AliGenCocktail*) GeneratorCocktail("Box_Nuclex_001_Extended");
+
+    const Int_t pdgcodes[6] = {
+        2212,
+        1000010020,
+        1000020030,
+        1000010030,
+        1000020040,
+        1010010030
+    };
+
+    const Char_t *names[6] = {
+        "Proton",
+        "Deuteron",
+        "Helium-3",
+        "Triton",
+        "Helium-4",
+        "Hypertriton"
+    };
+
+    int npart = 10;
+
+    double maxpt[6] = {10.0,10.0,10.0,10.0,10.0,10.0};
+
+    for (Int_t ipart = 0; ipart < 6; ++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(0., maxpt[ipart]);
+            box->SetPhiRange(0., 360.);
+            box->SetYRange(-0.8,0.8);
+            ctl->AddGenerator(box, Form("%s%i",names[ipart],negative), 1);
+        }
+    }
+  return ctl;
+}