Implementing Marcello comments

Author: stefanopolitano

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

@@ -1,4 +1,4 @@
-### authors: Stefano Politano (stefano.politano@cern.ch), Mattia Faggin (mattia.faggin@cern.ch), Fabrizio Grosa (fabrizio.grosa@cern.ch)
+### authors: Stefano Politano (stefano.politano@cern.ch), Mattia Faggin (mattia.faggin@cern.ch), Fabrizio Grosa (fabrizio.grosa@cern.ch), Marcello Di Costanzo (marcello.di.costanzo@cern.ch)
 ### last update:  April 2025
 
 ### beams
@@ -46,14 +46,15 @@ BeamRemnants:saturation 5
 ### D0 -> K- π+ (20%)
 421:oneChannel = 1 0.20 0 -321 211              ### D0 -> K- π+ 3.94% 
 ### D0 -> K- π+ π0 (20%)
-421:addChannel = 1 0.04 0 -321 211 111          ### D0 -> K- π+ π0 (non-resonant) 1.15% (e.g. 115/(115+231+195)*0.2 = 0.04)
-421:addChannel = 1 0.09 0 -313 211              ### D0 -> K*-(892) π+ 2.31%
-421:addChannel = 1 0.07 0 -323 111              ### D0 -> antiK*0(892) π0 1.95%
-### D0 -> π- π+ π0 (20%)
-421:addChannel = 1 0.18 0 -211 211              ### D0 -> π- π+ π0 (non-resonant) 1.0e-4
-421:addChannel = 1 0.02 0 -211 211              ### D0 -> rho+ π- 1.01%
+421:addChannel = 1 0.01 0 -321 211 111          ### D0 -> K- π+ π0 (non-resonant) 1.15% (e.g. 115/(115+231+195)*0.2 = 0.04)
+421:addChannel = 1 0.14 0  213 -211             ### D0 -> rho+ K- 11.2%
+421:addChannel = 1 0.03 0 -313 111              ### D0 -> antiK*0(892) π0 1.95%
+421:addChannel = 1 0.03 0 -323 211              ### D0 -> K*-(892) π+ 2.31%
 ### D0 -> π- π+ (20%)
-421:addChannel = 1 0.20 0 -211 211 111          ### D0 -> π- π+ (non-resonant) 1.45e-3
+421:addChannel = 1 0.20 0 -211 211              ### D0 -> π- π+ (non-resonant) 1.0e-4
+### D0 -> π- π+ π0 (20%)
+421:addChannel = 1 0.14 0  213 -211             ### D0 -> rho+ π- 1.01%
+421:addChannel = 1 0.06 0 -211 211 111          ### D0 -> π- π+ π0 (non-resonant) 1.3e-4
 ### D0 -> K- K+ (20%)
 421:addChannel = 1 0.20 0 -321 321              ### D0 -> K- K+ (non-resonant) 4.08e-3
 
@@ -109,11 +110,11 @@ BeamRemnants:saturation 5
 4232:oneChannel = 1 0.34 0 2212 -321 211        ### Ξc+ -> p K- π+
 ### Ξc+ -> p K- K+ (33%)
 4232:addChannel = 1 0.33 0 2212 333             ### Ξc+ -> p φ
-### Ξc+ -> sigma π+ π- (33%)
-4232:addChannel = 1 0.33 0 3222 -211 211        ### Ξc+ -> sigma π π
+### Ξc+ -> sigma+ π+ π- (33%)
+4232:addChannel = 1 0.33 0 3222 -211 211        ### Ξc+ -> sigma+ π- π+
 
 # Allow the decay of resonances in the decay chain
-### K*(892)0 -> K- π+
+### K*0(892) -> K- π+
 313:onMode = off
 313:onIfAll = 321 211
 ### K*(892)+ -> K- π0
@@ -122,7 +123,7 @@ BeamRemnants:saturation 5
 ### K*(1430)0 -> K- π0
 10311:onMode = off
 10311:onIfAll = 321 211
-### rho+ -> π- π0
+### rho+ -> π+ π0
 213:onMode = off
 213:onIfAll = 211 111
 ### φ -> K+ K-
@@ -159,8 +160,16 @@ BeamRemnants:saturation 5
 421:onIfMatch = 321 211
 ### D0 -> K π π0
 421:onIfMatch = 321 211 111
+### D0 -> rho K
+421:onIfMatch = 213 211
+### D0 -> antiK*0 π0
+421:onIfMatch = 313 111
+### D0 -> K*- π+
+421:onIfMatch = 323 211 
 ### D0 -> π π
 421:onIfMatch = 211 211
+### D0 -> rho+ π-
+421:onIfMatch = 213 -211
 ### D0 -> π π π0
 421:onIfMatch = 211 211 111
 ### D0 -> K K
@@ -170,8 +179,6 @@ BeamRemnants:saturation 5
 411:onIfMatch = 321 211 211
 ### D+/- -> K π π 
 411:onIfMatch = 321 211 211 111
-### D+/- -> K* π
-411:onIfMatch = 313 211
 ### D+/- -> K* K
 411:onIfMatch = 313 321
 ### D+/- -> antiK* π
@@ -226,6 +233,6 @@ BeamRemnants:saturation 5
 ### Ξc+ -> p K- π+
 4232:onIfMatch = 2212 321 211
 ### Ξc+ -> Ξ- π+ π+
-4232:onIfMatch = 3312 211 211
+4232:onIfMatch = 3222 211 211
 ### Ξc+ -> p φ
 4232:onIfMatch = 2212 333

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

@@ -1,6 +1,7 @@
 ### authors: Fabrizio Grosa (fabrizio.grosa@cern.ch)
 ###          Mattia Faggin (mattia.faggin@cern.ch)
-###          Stefano Politano (stefano.politano@cern.ch)       
+###          Stefano Politano (stefano.politano@cern.ch)
+###50        Marcello Di Costanzo (marcello.di.costanzo@cern.ch)
 ### last update:  June 2024
 
 ### beams
@@ -53,14 +54,15 @@ BeamRemnants:saturation 5
 ### D0 -> K- π+ (20%)
 421:oneChannel = 1 0.20 0 -321 211              ### D0 -> K- π+ 3.94% 
 ### D0 -> K- π+ π0 (20%)
-421:addChannel = 1 0.04 0 -321 211 111          ### D0 -> K- π+ π0 (non-resonant) 1.15% (e.g. 115/(115+231+195)*0.2 = 0.04)
-421:addChannel = 1 0.09 0 -313 211              ### D0 -> K*-(892) π+ 2.31%
-421:addChannel = 1 0.07 0 -323 111              ### D0 -> antiK*0(892) π0 1.95%
-### D0 -> π- π+ π0 (20%)
-421:addChannel = 1 0.18 0 -211 211              ### D0 -> π- π+ π0 (non-resonant) 1.0e-4
-421:addChannel = 1 0.02 0 -211 211              ### D0 -> rho+ π- 1.01%
+421:addChannel = 1 0.01 0 -321 211 111          ### D0 -> K- π+ π0 (non-resonant) 1.15% (e.g. 115/(115+231+195)*0.2 = 0.04)
+421:addChannel = 1 0.14 0  213 -211             ### D0 -> rho+ K- 11.2%
+421:addChannel = 1 0.03 0 -313 111              ### D0 -> antiK*0(892) π0 1.95%
+421:addChannel = 1 0.03 0 -323 211              ### D0 -> K*-(892) π+ 2.31%
 ### D0 -> π- π+ (20%)
-421:addChannel = 1 0.20 0 -211 211 111          ### D0 -> π- π+ (non-resonant) 1.45e-3
+421:addChannel = 1 0.20 0 -211 211              ### D0 -> π- π+ (non-resonant) 1.0e-4
+### D0 -> π- π+ π0 (20%)
+421:addChannel = 1 0.14 0  213 -211             ### D0 -> rho+ π- 1.01%
+421:addChannel = 1 0.06 0 -211 211 111          ### D0 -> π- π+ π0 (non-resonant) 1.3e-4
 ### D0 -> K- K+ (20%)
 421:addChannel = 1 0.20 0 -321 321              ### D0 -> K- K+ (non-resonant) 4.08e-3
 
@@ -116,11 +118,11 @@ BeamRemnants:saturation 5
 4232:oneChannel = 1 0.34 0 2212 -321 211        ### Ξc+ -> p K- π+
 ### Ξc+ -> p K- K+ (33%)
 4232:addChannel = 1 0.33 0 2212 333             ### Ξc+ -> p φ
-### Ξc+ -> sigma π+ π- (33%)
-4232:addChannel = 1 0.33 0 3222 -211 211        ### Ξc+ -> sigma π π
+### Ξc+ -> sigma+ π+ π- (33%)
+4232:addChannel = 1 0.33 0 3222 -211 211        ### Ξc+ -> sigma+ π- π+
 
 # Allow the decay of resonances in the decay chain
-### K*(892)0 -> K- π+
+### K*0(892) -> K- π+
 313:onMode = off
 313:onIfAll = 321 211
 ### K*(892)+ -> K- π0
@@ -129,7 +131,7 @@ BeamRemnants:saturation 5
 ### K*(1430)0 -> K- π0
 10311:onMode = off
 10311:onIfAll = 321 211
-### rho+ -> π- π0
+### rho+ -> π+ π0
 213:onMode = off
 213:onIfAll = 211 111
 ### φ -> K+ K-
@@ -166,8 +168,16 @@ BeamRemnants:saturation 5
 421:onIfMatch = 321 211
 ### D0 -> K π π0
 421:onIfMatch = 321 211 111
+### D0 -> rho K
+421:onIfMatch = 213 211
+### D0 -> antiK*0 π0
+421:onIfMatch = 313 111
+### D0 -> K*- π+
+421:onIfMatch = 323 211 
 ### D0 -> π π
 421:onIfMatch = 211 211
+### D0 -> rho+ π-
+421:onIfMatch = 213 -211
 ### D0 -> π π π0
 421:onIfMatch = 211 211 111
 ### D0 -> K K
@@ -177,8 +187,6 @@ BeamRemnants:saturation 5
 411:onIfMatch = 321 211 211
 ### D+/- -> K π π 
 411:onIfMatch = 321 211 211 111
-### D+/- -> K* π
-411:onIfMatch = 313 211
 ### D+/- -> K* K
 411:onIfMatch = 313 321
 ### D+/- -> antiK* π
@@ -233,6 +241,6 @@ BeamRemnants:saturation 5
 ### Ξc+ -> p K- π+
 4232:onIfMatch = 2212 321 211
 ### Ξc+ -> Ξ- π+ π+
-4232:onIfMatch = 3312 211 211
+4232:onIfMatch = 3222 211 211
 ### Ξc+ -> p φ
-4232:onIfMatch = 2212 333
+4232:onIfMatch = 2212 333