From 07b1b45457772ea3d0d381215aa47803a4c2b4b1 Mon Sep 17 00:00:00 2001
From: mvolkl <38325723+mvolkl@users.noreply.github.com>
Date: Wed, 1 Mar 2023 16:58:50 +0100
Subject: [PATCH] HFe enhanced pythia8 generator, passed basic tests

---
 .../PWGHF/Pythia8_Monash2013_HFe.C            | 29 +++++++++++++++++++
 1 file changed, 29 insertions(+)
 create mode 100644 MC/CustomGenerators/PWGHF/Pythia8_Monash2013_HFe.C

diff --git a/MC/CustomGenerators/PWGHF/Pythia8_Monash2013_HFe.C b/MC/CustomGenerators/PWGHF/Pythia8_Monash2013_HFe.C
new file mode 100644
index 00000000..bd054d11
--- /dev/null
+++ b/MC/CustomGenerators/PWGHF/Pythia8_Monash2013_HFe.C
@@ -0,0 +1,29 @@
+AliGenerator *GeneratorCustom(TString opt = "")
+{
+
+  //Switches for prompt/nonprompt
+  Process_t process[2] = {kPyCharmppMNRwmi, kPyBeautyppMNRwmi}; //charm or beauty
+  Int_t iprocess = uidConfig % 2;
+
+  AliGenPythiaPlus* pyth = (AliGenPythiaPlus*)GeneratorPythia8(kPythia8Tune_Monash2013);
+  pyth->SetProcess(process[iprocess]);
+  
+  TVirtualMCDecayer* decayer = new AliDecayerPythia();
+  decayer->SetForceDecay(kAll);
+  decayer->Init();
+  gMC->SetExternalDecayer(decayer);
+
+  // Pt transfer of the hard scattering (set for all cases)
+  //pyth->SetPtHard(pthardminConfig, pthardmaxConfig);
+
+  // From overview
+  pyth->SetCutOnChild(1);
+  pyth->SetPdgCodeParticleforAcceptanceCut(11);
+  pyth->SetChildYRange(-1.2, 1.2);
+  pyth->SetChildPtRange(0, 10000.);
+   
+  // lambda_b from PDG 2019: tau0 = 1.471 ps = 441 m/c = 0.441 mm/c
+  (AliPythia8::Instance())->ReadString("5122:tau0 = 4.41000e-01");
+  
+  return pyth;
+}