MC/PWGEM: no forced b-decays for HF Cocktail (#1711)

* no forced b-decays
* check for forced charm decays

Author: DanielSamitz

MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini

@@ -1,6 +1,6 @@
 [GeneratorExternal]
 fileName = ${O2DPG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C
-funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332")
+funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;")
 
 [GeneratorPythia8]
 config = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail.cfg
@@ -10,4 +10,4 @@ includePartonEvent=true
 
 [DecayerPythia8]
 config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
-config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
+config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg

MC/config/PWGEM/ini/Pythia8_Charm_Cocktail.ini

@@ -10,4 +10,4 @@ includePartonEvent=true
 
 [DecayerPythia8]
 config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg
-config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg
+config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg

MC/config/PWGEM/ini/tests/Pythia8_Beauty_Cocktail.C

@@ -1,6 +1,6 @@
 int External()
 {
-    	
+    
     std::string path{"o2sim_Kine.root"};
     TFile file(path.c_str(), "READ");
     if (file.IsZombie()) {
@@ -24,11 +24,11 @@ int External()
             if (pdg == -5){
                 found_bbar = true;
             }
-            if ( abs(pdg) == 511 || abs(pdg) == 521 || abs(pdg) == 531 || abs(pdg) == 5122 || abs(pdg) == 5132 || abs(pdg) == 5232 || abs(pdg) == 5332){
+            if ( abs(pdg) == 411 || abs(pdg) == 421 || abs(pdg) == 431 || abs(pdg) == 4122 || abs(pdg) == 4232 || abs(pdg) == 4132 || abs(pdg) == 4332){
                 int ifirstdaughter = track.getFirstDaughterTrackId();
                 int ilastdaughter = track.getLastDaughterTrackId();
                 if (ifirstdaughter == -1 || ilastdaughter == -1){
-                    std::cerr << "Found beauty hadron that does not have daughters" << "\n";
+                    std::cerr << "Found charm hadron that does not have daughters" << "\n";
                     return 1;
                 }
                 bool found_electron = false;
@@ -47,7 +47,7 @@ int External()
                     }
                 }
                 if (!found_electron){
-                    std::cerr << "Found beauty hadron that does not decay to electron" << "\n";
+                    std::cerr << "Found charm hadron that does not decay to electron" << "\n";
                     return 1;
                 }
             }
@@ -56,6 +56,6 @@ int External()
             std::cerr << "Found event without b-bbar pair" << "\n";
             return 1;
         }
-	}
+    }
     return 0;
 }

…M/pythia8/decayer/force_semileptonic.cfg …a8/decayer/force_semileptonic_beauty.cfgMC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg renamed to MC/config/PWGEM/pythia8/decayer/force_semileptonic_beauty.cfg MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg renamed to MC/config/PWGEM/pythia8/decayer/force_semileptonic_beauty.cfg

@@ -1,25 +1,3 @@
-### only semileptonic decays for charm
-### D+
-411:oneChannel = 1 0.087 0 -311 -11 12
-411:addChannel = 1 0.040 0 -321 211 -11 12
-411:addChannel = 1 0.037 0 -313 -11 12
-### D0
-421:oneChannel = 1 0.035 0 -321 -11 12
-421:addChannel = 1 0.022 0 -323 -11 12
-421:addChannel = 1 0.016 0 -321 111 -11 12
-421:addChannel = 1 0.014 0 -311 -211 -11 12
-### Ds
-431:oneChannel = 1 0.025 0 333 -11 12
-431:addChannel = 1 0.027 0 221 -11 12
-### Lambdac
-4122:oneChannel = 1 0.036 0 3122 -11 12
-### chi_{c}^{+}
-4232:oneChannel = 1 0.07 0 3322 -11 12
-### chi_{c}^{0}
-4132:oneChannel = 1 0.014 0 3312 -11 12
-### Omega_{c}
-4332:oneChannel = 1 0.01224 0 3334 -11 12
-
 ### only semileptonic decays for beauty
 ### B0
 511:oneChannel = 1 0.0207000 0 12 -11 -411

MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg

@@ -0,0 +1,26 @@
+### only semileptonic decays for charm
+### D+
+411:oneChannel = 1 0.087 0 -311 -11 12
+411:addChannel = 1 0.040 0 -321 211 -11 12
+411:addChannel = 1 0.037 0 -313 -11 12
+### D0
+421:oneChannel = 1 0.035 0 -321 -11 12
+421:addChannel = 1 0.022 0 -323 -11 12
+421:addChannel = 1 0.016 0 -321 111 -11 12
+421:addChannel = 1 0.014 0 -311 -211 -11 12
+### Ds
+431:oneChannel = 1 0.025 0 333 -11 12
+431:addChannel = 1 0.027 0 221 -11 12
+### Lambdac
+4122:oneChannel = 1 0.036 0 3122 -11 12
+### chi_{c}^{+}
+4232:oneChannel = 1 0.07 0 3322 -11 12
+### chi_{c}^{0}
+4132:oneChannel = 1 0.014 0 3312 -11 12
+### Omega_{c}
+4332:oneChannel = 1 0.01224 0 3334 -11 12
+
+# Correct OmegaC decay length (wrong in PYTHIA8 decay table) (mm/c)
+4332:tau0 = 0.08000000000
+# Correct Lb decay length (wrong in PYTHIA8 decay table)
+5122:tau0 = 4.41000e-01