Merge branch 'alisw:master' into master

alisw · Apr 11, 2023 · 82d2d2a · 82d2d2a
2 parents c8e2a25 + d8e8f25
commit 82d2d2a
Show file tree

Hide file tree

Showing 7 changed files with 430 additions and 17 deletions.
diff --git a/MC/CustomGenerators/PWGDQ/Pythia8_Monash_HighMassDimuon.C b/MC/CustomGenerators/PWGDQ/Pythia8_Monash_HighMassDimuon.C
@@ -0,0 +1,18 @@
+AliGenerator *GeneratorCustom(TString opt = "")
+
+{
+
+    AliGenPythiaPlus* pyth = (AliGenPythiaPlus*)GeneratorPythia8(kPythia8Tune_Monash2013);
+
+    pyth->SetProcess(kPyHeavyFlavppMNRwmi);
+
+    pyth->SetCutOnChild(1);
+    pyth->SetPdgCodeParticleforAcceptanceCut(13);
+    pyth->SetChildYRange(-4.1, -2.4);
+    pyth->SetChildPtRange(0, 10000.);
+
+    pyth->SetStackFillOpt(AliGenPythiaPlus::kHeavyFlavor);
+
+    return pyth;
+
+}
diff --git a/MC/CustomGenerators/PWGLF/Dpmjet_Rsn_pPb003.C b/MC/CustomGenerators/PWGLF/Dpmjet_Rsn_pPb003.C
@@ -0,0 +1,47 @@
+AliGenerator *GeneratorCustom(TString opt = "")
+{
+  // Custom generator to inject f0, f1 and f2
+  // Options:
+  // pPb:  flat pT distibution in a range 0-20 GeV/c; 4 particles per event;
+  // rapidity range: -0.6 < y < 0.1 : Anchored to period LHC16q
+
+  Int_t ninj = 4; // 4 injected particles per species
+
+  // Set pT range (pT_L, pT_H)
+  Double_t  pT_L = 0.0;
+  Double_t  pT_H = 20.0;
+
+  // Set Rapidity
+  Double_t y_L = -0.6;
+  Double_t y_H = 0.1;
+
+  AliGenCocktail *ctl  = (AliGenCocktail*) GeneratorCocktail("DPMJET_RsnInject");
+  AliGenDPMjet *dpmjet  = (AliGenDPMjet*) GeneratorPhojet();
+  ctl->AddGenerator(dpmjet,"Dpmjet",           1.);
+
+  Double_t randm = gRandom->Rndm();
+
+  if ( randm < 0.2) {
+    AliGenerator *injf0_1710_pPb     = GeneratorInjector(ninj, 10331, pT_L, pT_H,  y_L, y_H); // f0(1710)
+    ctl->AddGenerator(injf0_1710_pPb,"Injector(Rsnf0_1710)", 1.);
+  }
+  else if ( randm >= 0.2 && randm < 0.4) {
+    AliGenerator *injf2_1525_pPb     = GeneratorInjector(ninj,    335, pT_L, pT_H,  y_L, y_H); // f2(1525)
+    ctl->AddGenerator(injf2_1525_pPb,"Injector(Rsnf2_1525)", 1.);
+  }
+  else if ( randm >= 0.4 && randm < 0.6) {
+    AliGenerator *injf1_1285_pPb     = GeneratorInjector(ninj,    20223, pT_L, pT_H,  y_L, y_H); // f1(1285)
+    ctl->AddGenerator(injf1_1285_pPb,"Injector(Rsnf1_1285)", 1.);
+  }
+  else if ( randm >= 0.6 && randm < 0.8) {
+    AliGenerator *injf1_1420_pPb     = GeneratorInjector(ninj,    20333, pT_L, pT_H,  y_L, y_H); // f1(1420)
+    ctl->AddGenerator(injf1_1420_pPb,"Injector(Rsnf1_1420)", 1.);
+  }
+  else {
+    AliGenerator *injf0_1500_pPb     = GeneratorInjector(ninj,    90302, pT_L, pT_H,  y_L, y_H); // f0(1500)
+    ctl->AddGenerator(injf0_1500_pPb,"Injector(Rsnf2)", 1.);
+  }
+
+  return ctl;
+
+}
diff --git a/MC/CustomGenerators/PWGLF/Hijing_Nuclex006.C b/MC/CustomGenerators/PWGLF/Hijing_Nuclex006.C
@@ -5,10 +5,11 @@ GeneratorCustom()
   AliGenCocktail *ctl   = (AliGenCocktail*) GeneratorCocktail("Hijing_Nuclex006");
   AliGenerator   *hij   = GeneratorHijing();
   ctl->AddGenerator(hij,  "Hijing", 1.);
-  //	hypertriton, hyperhydrogen-4, hyperhelium-4, double-hyperhydrogen-4
-  AliGenerator   *nu1a  = Generator_Nuclex(0x10C010, kFALSE, 40);
-  AliGenerator   *nu1b  = Generator_Nuclex(0x10C010, kTRUE, 40);
+  //	hypertriton, hyperhydrogen-4, hyperhelium-4, double-hyperhydrogen-4, 4(Xi-)He
+  AliGenerator   *nu1a  = Generator_Nuclex((0x10 | 0x4000 | 0x8000 | 0x100000 | 0x200000), kFALSE, 40, 10., 1.);
+  AliGenerator   *nu1b  = Generator_Nuclex((0x10 | 0x4000 | 0x8000 | 0x100000 | 0x200000), kTRUE,  40, 10., 1.);
   ctl->AddGenerator(nu1a,  "Nuclex1a", 1.);
   ctl->AddGenerator(nu1b,  "Nuclex1b", 1.);
+
   return ctl;
 }
diff --git a/MC/CustomGenerators/PWGLF/Hijing_Nuclex007.C b/MC/CustomGenerators/PWGLF/Hijing_Nuclex007.C
@@ -1,20 +1,34 @@
-// Requires AliRoot from v5-09-59b --> fix for A=4 hypernuclei masses
+// Requires latest AliRoot --> fix for A=4 hypernuclei masses
 AliGenerator *
 GeneratorCustom()
 {
   AliGenCocktail *ctl   = (AliGenCocktail*) GeneratorCocktail("Hijing_Nuclex007");
   AliGenerator   *hij   = GeneratorHijing();
   ctl->AddGenerator(hij,  "Hijing", 1.);
-  //	hypertriton, hyperhydrogen-4, hyperhelium-4, double-hyperhydrogen-4
-  AliGenerator   *nu1a  = Generator_Nuclex((0x10 | 0x4000 | 0x8000 | 0x100000), kFALSE, 40);
-  AliGenerator   *nu1b  = Generator_Nuclex((0x10 | 0x4000 | 0x8000 | 0x100000), kTRUE,  40);
+
+  //	Exotica + hypertriton, hyperhydrogen-4, hyperhelium-4, double-hyperhydrogen-4, 4(Xi-)He
+  AliGenerator   *nu1a  = Generator_Nuclex((0x10 | 0x20 | 0x40 | 0x80 | 0x100 | 0x200 | 0x400 | 0x800 | 0x1000 | 0x2000 | 0x4000 | 0x8000 | 0x100000 | 0x200000), kFALSE, 40, 10., 1.);
+  AliGenerator   *nu1b  = Generator_Nuclex((0x10 | 0x20 | 0x40 | 0x80 | 0x100 | 0x200 | 0x400 | 0x800 | 0x1000 | 0x2000 | 0x4000 | 0x8000 | 0x100000 | 0x200000), kTRUE,  40, 10., 1.);
   ctl->AddGenerator(nu1a,  "Nuclex1a", 1.);
   ctl->AddGenerator(nu1b,  "Nuclex1b", 1.);
   //	helium3, triton, alpha
-  AliGenerator   *nu2a  = Generator_Nuclex((0x2 | 0x4 | 0x8), kFALSE, 10);
-  AliGenerator   *nu2b  = Generator_Nuclex((0x2 | 0x4 | 0x8), kTRUE,  10);
+  AliGenerator   *nu2a  = Generator_Nuclex((0x2 | 0x4 | 0x8), kFALSE, 10, 10., 1.);
+  AliGenerator   *nu2b  = Generator_Nuclex((0x2 | 0x4 | 0x8), kTRUE,  10, 10., 1.);
   ctl->AddGenerator(nu2a,  "Nuclex2a", 1.);
   ctl->AddGenerator(nu2b,  "Nuclex2b", 1.);
 
+  AliGenerator   *ilam  = GeneratorInjector(20, 3122,  0., 10., -1., 1.);
+  AliGenerator   *ixi   = GeneratorInjector(20, 3312,  0., 10., -1., 1.);
+  AliGenerator   *iom   = GeneratorInjector(20, 3334,  0., 10., -1., 1.);
+  AliGenerator   *ialam = GeneratorInjector(20, -3122, 0., 10., -1., 1.);
+  AliGenerator   *iaxi  = GeneratorInjector(20, -3312, 0., 10., -1., 1.);
+  AliGenerator   *iaom  = GeneratorInjector(20, -3334, 0., 10., -1., 1.);
+  ctl->AddGenerator(ilam,  "Injector (Lambda)", 1.);
+  ctl->AddGenerator(ixi,   "Injector (Xi)",     1.);
+  ctl->AddGenerator(iom,   "Injector (Omega)",  1.);
+  ctl->AddGenerator(ialam, "Injector (Lambda)", 1.);
+  ctl->AddGenerator(iaxi,  "Injector (Xi)",     1.);
+  ctl->AddGenerator(iaom,  "Injector (Omega)",  1.);
+
   return ctl;
 }
diff --git a/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C b/MC/CustomGenerators/PWGLF/Pythia8_Monash2013_Rsn_f0f1.C
@@ -0,0 +1,35 @@
+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 = 5;
+  particle_inj particleList[5] = { // {name,pdgcode,maxpt,maxy},
+    {1,"f0(1710)",10331,20.,0.9},
+    {1,"f2(1525)",335,20.,0.9},
+    {1,"f1(1285)",20223,20.,0.9},
+    {1,"f1(1420)",20333,20.,0.9},
+    {1,"f0(1500)",90302,20.,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;
+}