PWGEM: LMee cocktail tutorial (#2463)

* python script

* plotting macro

* config files

* parametrizations

* decaytable

* res/eff/mult/DCA files

* update plotting macro

Author: DanielSamitz

Tutorials/PWGEM/Cocktail/configs/ana_config.json

@@ -0,0 +1,18 @@
+{
+    "em-lmee-lf-cocktail": {
+        "cfgResFileName": "PATH/Cocktail/files/resolution/Resolution_PbPb5TeV.root",
+        "cfgEffFileName": "PATH/Cocktail/files/resolution/EffMult.root",
+        "cfgMultFileName": "PATH/Cocktail/files/resolution/EffMult.root",
+        "cfgDCAFileName": "PATH/Cocktail/files/resolution/DCA.root",
+        "cfgWriteTTree": "false",
+        "cfgDoPairing": "false",
+        "cfgMaxEta": "0.8",
+        "cfgMinPt": "0.2",
+        "cfgMaxPt": "10.0",
+        "cfgMinOpAng": "0",
+        "cfgALTweight": "1",
+        "cfgDoVirtPh": "false",
+        "cfgPhotonPtFileName": "",
+        "cfgPhotonPtDirName": ""
+    }
+}

Tutorials/PWGEM/Cocktail/configs/ini_example.ini

@@ -0,0 +1,3 @@
+[GeneratorExternal]
+fileName = ${O2DPG_ROOT}/MC/config/PWGEM/external/generator/GeneratorEMCocktailV2.C
+funcName=GenerateEMCocktail(400,0,3,63,"${O2DPG_ROOT}/MC/config/PWGEM/parametrizations/PbPb5TeV_central.json","5TeV_0005_wRatio_etatest",350,0.0,30.0,10000,1,1,0,0,"",0,1.1,"${O2DPG_ROOT}/MC/config/PWGEM/decaytables/decaytable_LMee.dat",0)

Tutorials/PWGEM/Cocktail/configs/sim_config.json

@@ -0,0 +1,21 @@
+{
+    "collisionSystem": "kPbPb",
+    "centrality": "k0005",
+    "selectedMothers": 63,
+    "paramFile": "PATH/Cocktail/files/parametrizations/PbPb_5TeV.json",
+    "paramFileDir": "5TeV_0005_wRatio",
+    "numberOfParticles": 350,
+    "minPt": 0.0,
+    "maxPt": 30.0,
+    "yGenRange": 1.1,
+    "externalDecayer": 1,
+    "decayMode": 3,
+    "dynamicalPtRange": 0,
+    "useLMeeDecaytable": "PATH/Cocktail/files/decaytables/decaytable_LMee.dat",
+    "weightingMode": 0,
+    "pythiaErrorTolerance": 10000,
+    "paramV2FileDir": "",
+    "useYWeights": 0,
+    "decayLongLived": 1,
+    "toFixEP": 0
+}

Tutorials/PWGEM/Cocktail/files/decaytables/decaytable_LMee.dat

@@ -0,0 +1,46 @@
+{
+    "13TeV_default": {
+        "histoMtScaleFactor": {
+            "331": 0.4,
+            "333": 0.13,
+            "443": 0.054,
+            "3212": 0.49,
+            "310": 1.0,
+            "2224": 1.0,
+            "2214": 1.0,
+            "1114": 1.0,
+            "2114": 1.0,
+            "213": 1.0,
+            "-213": 1.0,
+            "313": 1.0,
+            "130": 1.0,
+            "3122": 1.0
+        },
+        "111_pt": "((x*(66.3458/TMath::Power(TMath::Exp((-(0.218056*x))-(0.0842573*(x*x)))+(x/0.741916),6.22709)))/(1634.39/TMath::Power(TMath::Exp((-(0.346984*x))-(0.0798057*(x*x)))+(x/0.648328),6.04003)))*(2098.79/TMath::Power(TMath::Exp((-(0.391494*x))-(0.0443212*(x*x)))+(x/0.575932),6.08223))",
+        "221_pt": "(((x*(66.3458/TMath::Power(TMath::Exp((-(0.218056*x))-(0.0842573*(x*x)))+(x/0.741916),6.22709)))/(1634.39/TMath::Power(TMath::Exp((-(0.346984*x))-(0.0798057*(x*x)))+(x/0.648328),6.04003)))*(2098.79/TMath::Power(TMath::Exp((-(0.391494*x))-(0.0443212*(x*x)))+(x/0.575932),6.08223)))*(((0.504625*exp((((-(157384*x))-sqrt((x*x)+(0.5478530000*0.5478530000)))/sqrt(1-(-(157384*-(157384)))))/0.1))+(42.0992*(-(0.0112628*TMath::Power(1+((x/1.39355)*(x/1.39355)),-(0.919045))))))/(exp((((-(157384*x))-sqrt((x*x)+(0.1349770000*0.1349770000)))/sqrt(1-(-(157384*-(157384)))))/0.1)+(-(0.0112628*TMath::Power(1+((x/1.39355)*(x/1.39355)),-(0.919045))))))",
+        "113_pt": "(((x*(66.3458/TMath::Power(TMath::Exp((-(0.218056*x))-(0.0842573*(x*x)))+(x/0.741916),6.22709)))/(1634.39/TMath::Power(TMath::Exp((-(0.346984*x))-(0.0798057*(x*x)))+(x/0.648328),6.04003)))*(2098.79/TMath::Power(TMath::Exp((-(0.391494*x))-(0.0443212*(x*x)))+(x/0.575932),6.08223)))*((((0.37788/(1+exp(-((x-0.265203)/0.274547))))+(-(0.393597*TMath::Gaus(x,0.324349,0.819888))))+(-(0.230459*TMath::Gaus(x,1.23251,1.24114))))+0.412584)",
+        "223_pt": "(((x*(66.3458/TMath::Power(TMath::Exp((-(0.218056*x))-(0.0842573*(x*x)))+(x/0.741916),6.22709)))/(1634.39/TMath::Power(TMath::Exp((-(0.346984*x))-(0.0798057*(x*x)))+(x/0.648328),6.04003)))*(2098.79/TMath::Power(TMath::Exp((-(0.391494*x))-(0.0443212*(x*x)))+(x/0.575932),6.08223)))*((((0.435069/(1+exp(-((x-0.302485)/0.27614))))+(-(0.258853*TMath::Gaus(x,0.42493,0.693856))))+(-(0.202117*TMath::Gaus(x,1.07795,1.09857))))+0.242026)"
+    },
+    "13TeV_PreliminaryChargedPions": {
+        "histoMtScaleFactor": {
+            "331": 0.4,
+            "333": 0.13,
+            "443": 0.054,
+            "3212": 0.49,
+            "310": 1.0,
+            "2224": 1.0,
+            "2214": 1.0,
+            "1114": 1.0,
+            "2114": 1.0,
+            "213": 1.0,
+            "-213": 1.0,
+            "313": 1.0,
+            "130": 1.0,
+            "3122": 1.0
+        },
+        "111_pt": "x*(123.313*(pow((exp((0.171351*x)-(abs(0.00418378)*(x*x)))+(x/0.919801)),-6.29677)))",
+        "221_pt": "(x*(123.313*(pow((exp((0.171351*x)-(abs(0.00418378)*(x*x)))+(x/0.919801)),-(6.29677)))))*(((0.539894*exp((((0.994515*x)-sqrt((x*x)+(0.5478530000*0.5478530000)))/sqrt(1-(0.994515*0.994515)))/1.69543))+(-(0.0125445*(0.608582*TMath::Power(1+((x/22.1715)*(x/22.1715)),-(229.876))))))/(exp((((0.994515*x)-sqrt((x*x)+(0.1349770000*0.1349770000)))/sqrt(1-(0.994515*0.994515)))/1.69543)+(0.608582*TMath::Power(1+((x/22.1715)*(x/22.1715)),-(229.876)))))",
+        "113_pt": "(x*(123.313*(pow((exp((0.171351*x)-(abs(0.00418378)*(x*x)))+(x/0.919801)),-(6.29677)))))*((((0.37788/(1+exp(-((x-0.265203)/0.274547))))+(-(0.393597*TMath::Gaus(x,0.324349,0.819888))))+(-(0.230459*TMath::Gaus(x,1.23251,1.24114))))+0.412584)",
+        "223_pt": "(x*(123.313*(pow((exp((0.171351*x)-(abs(0.00418378)*(x*x)))+(x/0.919801)),-(6.29677)))))*((((0.435069/(1+exp(-((x-0.302485)/0.27614))))+(-(0.258853*TMath::Gaus(x,0.42493,0.693856))))+(-(0.202117*TMath::Gaus(x,1.07795,1.09857))))+0.242026)"
+    }
+}