Update decay table for b-hadron to open charm MC

Author: fgrosa

MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg

@@ -48,24 +48,12 @@ BeamRemnants:saturation 5
 ### add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other
 411:oneChannel = 1 0.0752 0 -321 211 211
 411:addChannel = 1 0.0104 0 -313 211
-411:addChannel = 1 0.0156 0 311 211
-411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi
+411:addChannel = 1 0.0100 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi
 ## add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other
 4122:oneChannel = 1 0.0196 100 2212 -313
 4122:addChannel = 1 0.0108 100 2224 -321
 4122:addChannel = 1 0.022 100 102134 211
 4122:addChannel = 1 0.035 0 2212 -321 211
-4122:addChannel = 1 0.0159 0 2212 311
-### add Xic+ decays absent in PYTHIA8 decay table
-4232:addChannel = 1 0.2 0 2212 -313
-4232:addChannel = 1 0.2 0 2212 -321 211
-4232:addChannel = 1 0.2 0 3324 211
-4232:addChannel = 1 0.2 0 3312 211 211
-### add Xic0 decays absent in PYTHIA8 decay table
-4132:addChannel = 1 0.0143 0 3312 211
-### add OmegaC decays absent in PYTHIA8 decay table
-4332:addChannel = 1 0.5 0 3334 211
-4332:addChannel = 1 0.5 0 3312 211
 
 ### K* -> K pi
 313:onMode = off
@@ -82,25 +70,12 @@ BeamRemnants:saturation 5
 ### for Lambda_c -> Lambda(1520) K-
 102134:onMode = off
 102134:onIfAll = 2212 321
-### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p 
-### and Omega_c -> pi Xi -> pi pi Lambda -> pi pi pi p
-3312:onMode = off
-3312:onIfAll = 3122 -211
-3122:onMode = off
-3122:onIfAll = 2212 -211
-### for Omega_c -> pi Omega -> pi K Lambda -> pi K pi p
-3334:onMode = off
-3334:onIfAll = 3122 -321
 
 ### switch off all decay channels
 411:onMode = off
 421:onMode = off
 431:onMode = off
 4122:onMode = off
-4232:onMode = off
-4132:onMode = off
-443:onMode = off
-4332:onMode = off
 511:onMode = off
 521:onMode = off
 531:onMode = off
@@ -131,60 +106,44 @@ BeamRemnants:saturation 5
 4122:onIfMatch = 102134 211
 ### Lambda_c -> p K pi
 4122:onIfMatch = 2212 321 211
-### Lambda_c -> pK0s
-4122:onIfMatch = 2212 311 
-
-### Xic+ -> pK*0
-4232:onIfMatch = 2212 313
-### Xic+ -> p K- pi+
-4232:onIfMatch = 2212 321 211
-### Xic+ -> Xi*0 pi+, Xi*->Xi- pi+
-4232:onIfMatch = 3324 211
-### Xic+ -> Xi- pi+ pi+
-4232:onIfMatch = 3312 211 211
-    
-### Xic0 -> Xi- pi+
-4132:onIfMatch = 3312 211
-
-### Omega_c -> Omega pi
-4332:onIfMatch = 3334 211
-### Omega_c -> Xi pi
-4332:onIfMatch = 3312 211
 
 ### Force also golden beauty hadrons decay modes for D2H studies
 ### add B0 decays
 511:oneChannel = 1 0.4 0 -411 211
-511:addChannel = 1 0.3 0 -413 211
-511:addChannel = 1 0.2 0 431 -211
+511:addChannel = 1 0.4 0 -413 211
 ### add B+ decays
-521:oneChannel = 1 0.8 0 -421 211
+521:oneChannel = 1 0.7 0 -421 211
 ### add Bs0 decays
-531:oneChannel = 1 0.8 0 -431 211
+531:oneChannel = 1 0.7 0 -431 211
 ### add Lb decays
-5122:oneChannel = 1 0.8 0 4122 -211
+5122:oneChannel = 1 0.6 0 4122 -211
 
 ### we also add channels useful for studies of partly reconstructed decays / correlated backgrounds
 ### add B0 decays
 511:addChannel = 1 0.1 0 -411 213
+511:addChannel = 1 0.1 0 -411 321
 ### add B+ decays
 521:addChannel = 1 0.1 0 -421 213
 521:addChannel = 1 0.1 0 -423 211
+521:addChannel = 1 0.1 0 -421 321
 ### add Bs0 decays
 531:addChannel = 1 0.1 0 -431 213
 531:addChannel = 1 0.1 0 -433 211
+531:addChannel = 1 0.1 0 -433 321
 ### add Lb decays
-5122:addChannel = 1 0.1 0 4122 -213
-5122:addChannel = 1 0.05 0 4122 211 -211 -211
-5122:addChannel = 1 0.05 0 4212 -211
+5122:addChannel = 1 0.10 0 14122 -211
+5122:addChannel = 1 0.10 0 4124 -211
+5122:addChannel = 1 0.10 0 4212 -211
+5122:addChannel = 1 0.10 0 4122 -321
 
 ### B0 -> D pi
 511:onIfMatch = 411 211
 ### B0 -> D* pi
 511:onIfMatch = 413 211
-### B0 -> Ds pi
-511:onIfMatch = 431 211
 ### B0 -> D rho
 511:onIfMatch = 411 213
+### B0 -> D rho
+511:onIfMatch = 411 321
 
 ### B+ -> D0 pi
 521:onIfMatch = 421 211
@@ -202,9 +161,9 @@ BeamRemnants:saturation 5
 
 ### Lb -> Lc pi
 5122:onIfMatch = 4122 211
-### Lb -> Lc rho
-5122:onIfMatch = 4122 213
-### Lb -> Lc pi pi pi
-5122:onIfMatch = 4122 211 211 211
-### Lb -> Sc pi
+### Lb -> Lc(2595)+ pi
+5122:onIfMatch = 14122 211
+### Lb -> Lc(2625)+ pi
+5122:onIfMatch = 4124 211
+### Lb -> Sc(2455)+ pi
 5122:onIfMatch = 4212 211