Skip to content

[PWGHF] Updating config for pp 13.6, 5.36 TeV and PbPb 5.36 TeV with correlated bkg #1986

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
stefanopolitano wants to merge 1 commit into from
Closed

[PWGHF] Updating config for pp 13.6, 5.36 TeV and PbPb 5.36 TeV with correlated bkg #1986

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
132 changes: 107 additions & 25 deletions MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -41,45 +41,64 @@ BeamRemnants:saturation 5
# OmegaC
4332:tau0 = 0.0803

### Force golden charm hadrons decay modes for D2H studies
### add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other
421:oneChannel = 1 0.0389 0 -321 211
421:addChannel = 1 0.00389 0 -321 211 111
### 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
## 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
### HF decays
### BR are set to yield 50% the golden channel and 50% the correlated bkg ones
###
### D0 decays
421:oneChannel = 1 0.500 0 -321 211 ### D0 -> K- π+
421:addChannel = 1 0.100 0 -321 211 111 ### D0 -> K- π+ π0
421:addChannel = 1 0.100 0 -211 211 ### D0 -> π- π+
421:addChannel = 1 0.100 0 -211 211 111 ### D0 -> π- π+ π0
421:addChannel = 1 0.100 0 -321 321 ### D0 -> K- K+
421:addChannel = 1 0.100 0 -321 321 111 ### D0 -> K- K+ π0

### D+ decays
411:oneChannel = 1 0.50 0 -321 211 211 ### D+ -> K- π+ π+
411:addChannel = 1 0.17 0 -211 211 211 ### D+ -> π- π+ π+
411:addChannel = 1 0.16 0 313 211 ### D+ -> K* π+ -> K- π+ π+
411:addChannel = 1 0.17 0 333 211 ### D+ -> φ π+ -> K- K+ π+

### Ds+ decays
431:oneChannel = 1 0.500 0 333 211 ### Ds+ -> φ π+ -> K- K+ π+
431:addChannel = 1 0.125 0 313 321 ### Ds+ -> K* K+ -> K- K+ π+
431:addChannel = 1 0.125 0 -211 211 211 ### Ds+ -> π- π+ π+
431:addChannel = 1 0.125 0 -211 321 211 ### Ds+ -> π- K+ π+
431:addChannel = 1 0.125 0 -321 321 321 ### Ds+ -> K- K+ K+

## Lc decays
4122:oneChannel = 1 0.2500 2212 -321 211 ### Λc+ -> p K- π+
4122:addChannel = 1 0.0560 100 2212 -313 ### Λc+ -> p K* -> p K- π+
4122:addChannel = 1 0.0550 100 2224 -321 ### Λc+ -> Delta++ K- -> p K- π+
4122:addChannel = 1 0.0550 100 102134 211 ### Λc+ -> Lambda(1520) K- -> p K- π+
4122:addChannel = 1 0.1670 0 2212 211 211 ### Λc+ -> p π+ π+
4122:addChannel = 1 0.2500 0 2212 311 ### Λc+ -> p K0s -> p π+ π-
4122:addChannel = 1 0.1670 0 2212 -321 321 ### Λc+ -> p K- K+

### Ξc+ decays
4232:oneChannel = 1 0.50 0 2212 -321 211 ### Ξc+ -> p K- π+
4232:addChannel = 1 0.25 0 -3312 211 211 ### Ξc+ -> Ξ- π+ π+
4232:addChannel = 1 0.25 0 2212 333 ### Ξc+ -> p φ -> p K- K+

### 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
# Allow the decay of resonances in the decay chain
### K* -> K π
313:onMode = off
313:onIfAll = 321 211
### for Ds -> Phi pi+
### for Ds -> φ π+
333:onMode = off
333:onIfAll = 321 321
### for D0 -> rho0 pi+ k-
### for D0 -> rho0 π+ k-
113:onMode = off
113:onIfAll = 211 211
### for Lambda_c -> Delta++ K-
### for Λc -> Delta++ K-
2224:onMode = off
2224:onIfAll = 2212 211
### for Lambda_c -> Lambda(1520) K-
### for Λc -> Lambda(1520) K-
102134:onMode = off
102134:onIfAll = 2212 321
### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
Expand All @@ -92,6 +111,69 @@ BeamRemnants:saturation 5
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

# Allow the decay of HF
### D0 -> K π
421:onIfMatch = 321 211
### D0 -> K π π0
421:onIfMatch = 321 211 111
### D0 -> π π
421:onIfMatch = 211 211
### D0 -> π π π0
421:onIfMatch = 211 211 111
### D0 -> K K
421:onIfMatch = 321 321
### D0 -> K K π0
421:onIfMatch = 321 321 111

### D+/- -> K π π
411:onIfMatch = 321 211 211
### D+/- -> K* π
411:onIfMatch = 313 211
### D+/- -> φ π
411:onIfMatch = 333 211
### D+/- -> π π π
411:onIfMatch = 211 211 211

### Ds -> φ π
431:onIfMatch = 333 211
### Ds -> K* K
431:onIfMatch = 321 313
### Ds -> π π π
431:onIfMatch = 211 211 211
### Ds -> K π π
431:onIfMatch = 321 211 211
### Ds -> K K K
431:onIfMatch = 321 321 321

### Λc -> p K π
4122:onIfMatch = 2212 321 211
### Λc -> p K*
4122:onIfMatch = 2212 313
### Λc -> Delta++ K
4122:onIfMatch = 2224 321
### Λc -> Lambda(1520) π
4122:onIfMatch = 102134 211
### Λc -> p π π
4122:onIfMatch = 2212 211 211
### Λc -> pK0s
4122:onIfMatch = 2212 311
### Λc -> p K K
4122:onIfMatch = 2212 321 321

### Ξc+ -> p K- π+
4232:onIfMatch = 2212 321 211
### Ξc+ -> Ξ- π+ π+
4232:onIfMatch = 3312 211 211
### Ξc+ -> p φ
4232:onIfMatch = 2212 333

### switch off all decay channels
411:onMode = off
421:onMode = off
Expand Down
132 changes: 107 additions & 25 deletions MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_hardQCD_5TeV.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,45 +45,64 @@ BeamRemnants:saturation 5
# OmegaC
4332:tau0 = 0.0803

### Force golden charm hadrons decay modes for D2H studies
### add D0 decays absent in PYTHIA8 decay table and set BRs from PDG for other
421:oneChannel = 1 0.0389 0 -321 211
421:addChannel = 1 0.00389 0 -321 211 111
### 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
## 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
### HF decays
### BR are set to yield 50% the golden channel and 50% the correlated bkg ones
###
### D0 decays
421:oneChannel = 1 0.500 0 -321 211 ### D0 -> K- π+
421:addChannel = 1 0.100 0 -321 211 111 ### D0 -> K- π+ π0
421:addChannel = 1 0.100 0 -211 211 ### D0 -> π- π+
421:addChannel = 1 0.100 0 -211 211 111 ### D0 -> π- π+ π0
421:addChannel = 1 0.100 0 -321 321 ### D0 -> K- K+
421:addChannel = 1 0.100 0 -321 321 111 ### D0 -> K- K+ π0

### D+ decays
411:oneChannel = 1 0.50 0 -321 211 211 ### D+ -> K- π+ π+
411:addChannel = 1 0.17 0 -211 211 211 ### D+ -> π- π+ π+
411:addChannel = 1 0.16 0 313 211 ### D+ -> K* π+ -> K- π+ π+
411:addChannel = 1 0.17 0 333 211 ### D+ -> φ π+ -> K- K+ π+

### Ds+ decays
431:oneChannel = 1 0.500 0 333 211 ### Ds+ -> φ π+ -> K- K+ π+
431:addChannel = 1 0.125 0 313 321 ### Ds+ -> K* K+ -> K- K+ π+
431:addChannel = 1 0.125 0 -211 211 211 ### Ds+ -> π- π+ π+
431:addChannel = 1 0.125 0 -211 321 211 ### Ds+ -> π- K+ π+
431:addChannel = 1 0.125 0 -321 321 321 ### Ds+ -> K- K+ K+

## Lc decays
4122:oneChannel = 1 0.2500 2212 -321 211 ### Λc+ -> p K- π+
4122:addChannel = 1 0.0560 100 2212 -313 ### Λc+ -> p K* -> p K- π+
4122:addChannel = 1 0.0550 100 2224 -321 ### Λc+ -> Delta++ K- -> p K- π+
4122:addChannel = 1 0.0550 100 102134 211 ### Λc+ -> Lambda(1520) K- -> p K- π+
4122:addChannel = 1 0.1670 0 2212 211 211 ### Λc+ -> p π+ π+
4122:addChannel = 1 0.2500 0 2212 311 ### Λc+ -> p K0s -> p π+ π-
4122:addChannel = 1 0.1670 0 2212 -321 321 ### Λc+ -> p K- K+

### Ξc+ decays
4232:oneChannel = 1 0.50 0 2212 -321 211 ### Ξc+ -> p K- π+
4232:addChannel = 1 0.25 0 -3312 211 211 ### Ξc+ -> Ξ- π+ π+
4232:addChannel = 1 0.25 0 2212 333 ### Ξc+ -> p φ -> p K- K+

### 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
# Allow the decay of resonances in the decay chain
### K* -> K π
313:onMode = off
313:onIfAll = 321 211
### for Ds -> Phi pi+
### for Ds -> φ π+
333:onMode = off
333:onIfAll = 321 321
### for D0 -> rho0 pi+ k-
### for D0 -> rho0 π+ k-
113:onMode = off
113:onIfAll = 211 211
### for Lambda_c -> Delta++ K-
### for Λc -> Delta++ K-
2224:onMode = off
2224:onIfAll = 2212 211
### for Lambda_c -> Lambda(1520) K-
### for Λc -> Lambda(1520) K-
102134:onMode = off
102134:onIfAll = 2212 321
### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p
Expand All @@ -96,6 +115,69 @@ BeamRemnants:saturation 5
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

# Allow the decay of HF
### D0 -> K π
421:onIfMatch = 321 211
### D0 -> K π π0
421:onIfMatch = 321 211 111
### D0 -> π π
421:onIfMatch = 211 211
### D0 -> π π π0
421:onIfMatch = 211 211 111
### D0 -> K K
421:onIfMatch = 321 321
### D0 -> K K π0
421:onIfMatch = 321 321 111

### D+/- -> K π π
411:onIfMatch = 321 211 211
### D+/- -> K* π
411:onIfMatch = 313 211
### D+/- -> φ π
411:onIfMatch = 333 211
### D+/- -> π π π
411:onIfMatch = 211 211 211

### Ds -> φ π
431:onIfMatch = 333 211
### Ds -> K* K
431:onIfMatch = 321 313
### Ds -> π π π
431:onIfMatch = 211 211 211
### Ds -> K π π
431:onIfMatch = 321 211 211
### Ds -> K K K
431:onIfMatch = 321 321 321

### Λc -> p K π
4122:onIfMatch = 2212 321 211
### Λc -> p K*
4122:onIfMatch = 2212 313
### Λc -> Delta++ K
4122:onIfMatch = 2224 321
### Λc -> Lambda(1520) π
4122:onIfMatch = 102134 211
### Λc -> p π π
4122:onIfMatch = 2212 211 211
### Λc -> pK0s
4122:onIfMatch = 2212 311
### Λc -> p K K
4122:onIfMatch = 2212 321 321

### Ξc+ -> p K- π+
4232:onIfMatch = 2212 321 211
### Ξc+ -> Ξ- π+ π+
4232:onIfMatch = 3312 211 211
### Ξc+ -> p φ
4232:onIfMatch = 2212 333

### switch off all decay channels
411:onMode = off
421:onMode = off
Expand Down
Loading