Set the D+ -> φ π+ to 25% and average the rest, move the D+ -> K*0(892) pi+ to where the right final state is

Author: wuctlby

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

@@ -42,7 +42,7 @@ BeamRemnants:saturation 5
 4332:tau0 = 0.0803
 
 ### Force golden charm hadrons decay modes for D2H studies
-### HF decays
+# HF decays
 ### BR are set to yield 50% of signal in golden channel and uniform abundance of corr. bkg channels (weighted by BR from PDG) 
 ###
 ### D0 decays
@@ -63,20 +63,20 @@ BeamRemnants:saturation 5
 
 ### D+ decays
 ### D+ -> K- π+ π+ (50%)
-411:oneChannel = 1 0.44120 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
-411:addChannel = 1 0.05880 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
+411:oneChannel = 1 0.40189 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
+411:addChannel = 1 0.05356 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
+411:addChannel = 1 0.04455 0 -313 211              ### D+ -> K*0(892) pi+ 1.04%
 ### D+ -> K- π+ π+ π0 (small, 3%)
 411:addChannel = 1 0.03000 0 -321 211 211 111      ### D+ -> K- π+ π+ π0 6.25%
-### D+ -> K- K+ π+ (25.00%)
-411:addChannel = 1 0.06243 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
-411:addChannel = 1 0.00764 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
-411:addChannel = 1 0.04553 0 -313 211              ### D+ -> K*0(892) pi+
-411:addChannel = 1 0.06757 0  333 211              ### D+ -> φ π+ 2.69e-3
-411:addChannel = 1 0.06683 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
-### D+ -> π- π+ π+ (22.00%)
-411:addChannel = 1 0.05824 0  113 211               ### D+ -> rho0 π+ 8.4e-4
-411:addChannel = 1 0.03235 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
-411:addChannel = 1 0.12941 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
+### D+ -> K- K+ π+ (36.00%, set 25% for D+ -> φ π+, 11% for the rest)
+411:addChannel = 1 0.25000 0  333 211              ### D+ -> φ π+ 2.69e-3 !needed for signal
+411:addChannel = 1 0.03929 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
+411:addChannel = 1 0.02865 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
+411:addChannel = 1 0.04206 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
+### D+ -> π- π+ π+ (11.00%)
+411:addChannel = 1 0.02911 0  113 211               ### D+ -> rho0 π+ 8.4e-4
+411:addChannel = 1 0.01618 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
+411:addChannel = 1 0.06471 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
 
 ### Ds+ decays
 ### Ds+ -> K- K+ π+ (50%)
@@ -175,7 +175,7 @@ BeamRemnants:saturation 5
 3334:onMode = off
 3334:onIfAll = 3122 -321
 
-### switch off all decay channels
+# switch off all decay channels
 411:onMode = off
 421:onMode = off
 431:onMode = off

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

@@ -46,7 +46,7 @@ BeamRemnants:saturation 5
 4332:tau0 = 0.0803
 
 ### Force golden charm hadrons decay modes for D2H studies
-### HF decays
+# HF decays
 ### BR are set to yield 50% of signal in golden channel and uniform abundance of corr. bkg channels (weighted by BR from PDG) 
 ###
 ### D0 decays
@@ -67,20 +67,20 @@ BeamRemnants:saturation 5
 
 ### D+ decays
 ### D+ -> K- π+ π+ (50%)
-411:oneChannel = 1 0.44120 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
-411:addChannel = 1 0.05880 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
+411:oneChannel = 1 0.40189 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
+411:addChannel = 1 0.05356 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
+411:addChannel = 1 0.04455 0 -313 211              ### D+ -> K*0(892) pi+ 1.04%
 ### D+ -> K- π+ π+ π0 (small, 3%)
 411:addChannel = 1 0.03000 0 -321 211 211 111      ### D+ -> K- π+ π+ π0 6.25%
-### D+ -> K- K+ π+ (25.00%)
-411:addChannel = 1 0.06243 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
-411:addChannel = 1 0.00764 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
-411:addChannel = 1 0.04553 0 -313 211              ### D+ -> K*0(892) pi+
-411:addChannel = 1 0.06757 0  333 211              ### D+ -> φ π+ 2.69e-3
-411:addChannel = 1 0.06683 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
-### D+ -> π- π+ π+ (22.00%)
-411:addChannel = 1 0.05824 0  113 211               ### D+ -> rho0 π+ 8.4e-4
-411:addChannel = 1 0.03235 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
-411:addChannel = 1 0.12941 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
+### D+ -> K- K+ π+ (36.00%, set 25% for D+ -> φ π+, 11% for the rest)
+411:addChannel = 1 0.25000 0  333 211              ### D+ -> φ π+ 2.69e-3 !needed for signal
+411:addChannel = 1 0.03929 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
+411:addChannel = 1 0.02865 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
+411:addChannel = 1 0.04206 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
+### D+ -> π- π+ π+ (11.00%)
+411:addChannel = 1 0.02911 0  113 211               ### D+ -> rho0 π+ 8.4e-4
+411:addChannel = 1 0.01618 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
+411:addChannel = 1 0.06471 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
 
 ### Ds+ decays
 ### Ds+ -> K- K+ π+ (50%)
@@ -179,7 +179,7 @@ BeamRemnants:saturation 5
 3334:onMode = off
 3334:onIfAll = 3122 -321
 
-### switch off all decay channels
+# Switch off all decay channels
 411:onMode = off
 421:onMode = off
 431:onMode = off

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

@@ -42,7 +42,7 @@ BeamRemnants:saturation 5
 4332:tau0 = 0.0803
 
 ### Force golden charm hadrons decay modes for D2H studies
-### HF decays
+# HF decays
 ### BR are set to yield 50% of signal in golden channel and uniform abundance of corr. bkg channels (weighted by BR from PDG) 
 ###
 ### D0 decays
@@ -63,20 +63,20 @@ BeamRemnants:saturation 5
 
 ### D+ decays
 ### D+ -> K- π+ π+ (50%)
-411:oneChannel = 1 0.44120 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
-411:addChannel = 1 0.05880 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
+411:oneChannel = 1 0.40189 0 -321 211 211          ### D+ -> K- π+ π+ 9.38%
+411:addChannel = 1 0.05356 0 -10311 211            ### D+ -> antiK*0(1430) π+ 1.25%
+411:addChannel = 1 0.04455 0 -313 211              ### D+ -> K*0(892) pi+ 1.04%
 ### D+ -> K- π+ π+ π0 (small, 3%)
 411:addChannel = 1 0.03000 0 -321 211 211 111      ### D+ -> K- π+ π+ π0 6.25%
-### D+ -> K- K+ π+ (25.00%)
-411:addChannel = 1 0.06243 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
-411:addChannel = 1 0.00764 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
-411:addChannel = 1 0.04553 0 -313 211              ### D+ -> K*0(892) pi+
-411:addChannel = 1 0.06757 0  333 211              ### D+ -> φ π+ 2.69e-3
-411:addChannel = 1 0.06683 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
-### D+ -> π- π+ π+ (22.00%)
-411:addChannel = 1 0.05824 0  113 211               ### D+ -> rho0 π+ 8.4e-4
-411:addChannel = 1 0.03235 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
-411:addChannel = 1 0.12941 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
+### D+ -> K- K+ π+ (36.00%, set 25% for D+ -> φ π+, 11% for the rest)
+411:addChannel = 1 0.25000 0  333 211              ### D+ -> φ π+ 2.69e-3 !needed for signal
+411:addChannel = 1 0.03929 0 -313 321              ### D+ -> K*0(892) K+ 2.49e-3
+411:addChannel = 1 0.02865 0 -10311 321            ### D+ -> antiK*0(1430) K+ 1.82e-3
+411:addChannel = 1 0.04206 0 -321 321 211          ### D+ -> K- K+ π+ (non-resonant) 2.68e-3 
+### D+ -> π- π+ π+ (11.00%)
+411:addChannel = 1 0.02911 0  113 211               ### D+ -> rho0 π+ 8.4e-4
+411:addChannel = 1 0.01618 0  225 211               ### D+ -> f2(1270) π+ 4.6e-4
+411:addChannel = 1 0.06471 0 -211 211 211           ### D+ -> π- π+ π+ (non-resonant) 1.0e-4
 
 ### Ds+ decays
 ### Ds+ -> K- K+ π+ (50%)
@@ -175,7 +175,7 @@ BeamRemnants:saturation 5
 3334:onMode = off
 3334:onIfAll = 3122 -321
 
-### switch off all decay channels
+# switch off all decay channels
 411:onMode = off
 421:onMode = off
 431:onMode = off