[PWGLF]: add cfg and ini for deuteron production (#1881)

* add cfg and ini for deuteron production

* Add macro for ini

* Update pythia8_inel_136tev_wDeCoal.cfg

---------

Co-authored-by: Nicolò Jacazio <njacazio@users.noreply.github.com>

Author: ercolessi

MC/config/PWGLF/ini/GeneratorLF_pp1360_wDeCoal.ini

@@ -0,0 +1,6 @@
+[GeneratorExternal]
+fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_ALICE3.C
+funcName=generator_pythia8_ALICE3()
+
+[GeneratorPythia8]
+config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGLF/pythia8/generator/pythia8_inel_136tev_wDeCoal.cfg

MC/config/PWGLF/ini/GeneratorLF_pp1360_wDeInel.ini

@@ -0,0 +1,6 @@
+[GeneratorExternal]
+fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/ALICE3/pythia8/generator_pythia8_ALICE3.C
+funcName=generator_pythia8_ALICE3()
+
+[GeneratorPythia8]
+config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGLF/pythia8/generator/pythia8_inel_136tev_wDeInel.cfg

MC/config/PWGLF/ini/tests/GeneratorLF_pp1360_wDeCoal.C

@@ -0,0 +1,28 @@
+int External()
+{
+    std::string path{"o2sim_Kine.root"};
+
+    TFile file(path.c_str(), "READ");
+    if (file.IsZombie())
+    {
+        std::cerr << "Cannot open ROOT file " << path << "\n";
+        return 1;
+    }
+
+    auto tree = (TTree *)file.Get("o2sim");
+    if (!tree)
+    {
+        std::cerr << "Cannot find tree o2sim in file " << path << "\n";
+        return 1;
+    }
+    std::vector<o2::MCTrack> *tracks{};
+    tree->SetBranchAddress("MCTrack", &tracks);
+
+    auto nEvents = tree->GetEntries();
+    if (nEvents < 1)
+    {
+        std::cerr << "No events actually generated: not OK!";
+        return 1;
+    }
+    return 0;
+}

MC/config/PWGLF/ini/tests/GeneratorLF_pp1360_wDeInel.C

@@ -0,0 +1,28 @@
+int External()
+{
+    std::string path{"o2sim_Kine.root"};
+
+    TFile file(path.c_str(), "READ");
+    if (file.IsZombie())
+    {
+        std::cerr << "Cannot open ROOT file " << path << "\n";
+        return 1;
+    }
+
+    auto tree = (TTree *)file.Get("o2sim");
+    if (!tree)
+    {
+        std::cerr << "Cannot find tree o2sim in file " << path << "\n";
+        return 1;
+    }
+    std::vector<o2::MCTrack> *tracks{};
+    tree->SetBranchAddress("MCTrack", &tracks);
+
+    auto nEvents = tree->GetEntries();
+    if (nEvents < 1)
+    {
+        std::cerr << "No events actually generated: not OK!";
+        return 1;
+    }
+    return 0;
+}

MC/config/PWGLF/pythia8/generator/pythia8_inel_136tev_wDeCoal.cfg

@@ -0,0 +1,23 @@
+### Specify beams
+Beams:idA = 2212
+Beams:idB = 2212
+Beams:eCM = 13600.              ### energy
+
+Beams:frameType = 1
+ParticleDecays:limitTau0 = on
+ParticleDecays:tau0Max = 10.      ### match alice: 1cm/c = 10.0mm/c
+
+### processes
+SoftQCD:inelastic = on # all inelastic processes
+
+# default: do nothing, Monash 2013 will do its thing
+Tune:pp = 14
+
+Random:setSeed = on
+
+# enable deuteron production by coalescence collisions
+HadronLevel:DeuteronProduction = on
+DeuteronProduction:channels = {2212 2112 > 22}
+DeuteronProduction:models = {0}
+DeuteronProduction:norm = 1
+DeuteronProduction:parms = {0.195 1}

MC/config/PWGLF/pythia8/generator/pythia8_inel_136tev_wDeInel.cfg

@@ -0,0 +1,19 @@
+### Specify beams
+Beams:idA = 2212
+Beams:idB = 2212
+Beams:eCM = 13600.              ### energy
+
+Beams:frameType = 1
+ParticleDecays:limitTau0 = on
+ParticleDecays:tau0Max = 10.      ### match alice: 1cm/c = 10.0mm/c
+
+### processes
+SoftQCD:inelastic = on # all inelastic processes
+
+# default: do nothing, Monash 2013 will do its thing
+Tune:pp = 14
+
+Random:setSeed = on
+
+# enable deuteron production by inelastic collisions
+HadronLevel:DeuteronProduction = on