Skip to content

Update decay table for b-hadron to open charm MC #2472

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

Merged
stefanopolitano merged 3 commits into from
Sep 24, 2026
+27 −69
Merged

Update decay table for b-hadron to open charm MC #2472

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
a2e65c9
Update decay table for b-hadron to open charm MC
fgrosa Sep 24, 2026
6093c1f
Update test macro
fgrosa Sep 24, 2026
f1cb404
Fix typo
fgrosa Sep 24, 2026
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
13 changes: 5 additions & 8 deletions MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Bforced_gap5_Mode2.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,17 +6,14 @@ int External() {

std::vector<int> checkPdgHadron{411, 421, 431, 4122, 4132, 4232, 4332, 511, 521, 531, 5122};
std::map<int, std::vector<std::vector<int>>> checkHadronDecays{ // sorted pdg of daughters
{411, {{-321, 211, 211}, {-313, 211}, {211, 311}, {211, 333}}}, // D+
{411, {{-321, 211, 211}, {-313, 211}, {211, 333}}}, // D+
{421, {{-321, 211}, {-321, 111, 211}}}, // D0
{431, {{211, 333}, {-313, 321}}}, // Ds+
{4122, {{-313, 2212}, {-321, 2224}, {211, 102134}, {-321, 211, 2212}, {311, 2212}}}, // Lc+
{4132, {{211, 3312}}}, // Xic0
{4232, {{-313, 2212}, {-321, 3324}, {211, 211, 3312}, {-321, 211, 2212}}}, // Xic+
{4332, {{211, 3334}}}, // Omegac+
{511, {{-411, 211}, {-413, 211}, {-411, 213}, {431, -211}}}, // B0
{521, {{-421, 211}, {-423, 211}, {-421, 213}}}, // B+
{4122, {{-313, 2212}, {-321, 2224}, {211, 102134}, {-321, 211, 2212}}}, // Lc+
{511, {{-411, 211}, {-413, 211}, {-411, 213}, {-411, 321}, {-413, 213}, {-413, 321}}}, // B0
{521, {{-421, 211}, {-423, 211}, {-421, 213}, {-421, 321}}}, // B+
{531, {{-431, 211}, {-433, 211}, {-431, 213}}}, // Bs0
{5122, {{4122, -211}, {4122, -213}, {4122, 211, -211, -211}, {4212, -211}}} // Lb0
{5122, {{4122, -211}, {4122, -213}, {14122, -211}, {4124, -211}, {4212, -211}}} // Lb0
};

TFile file(path.c_str(), "READ");
Expand Down
83 changes: 22 additions & 61 deletions MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand All @@ -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
Expand Down Expand Up @@ -131,60 +106,46 @@ 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:oneChannel = 1 0.35 0 -411 211
511:addChannel = 1 0.35 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 -413 213
511:addChannel = 1 0.05 0 -411 321
511:addChannel = 1 0.05 0 -413 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
Expand All @@ -202,9 +163,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
Loading