Skip to content

Add configurations for MC production for SigmaC background studies. #2063

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
mfaggin wants to merge 1 commit into from
Closed

Add configurations for MC production for SigmaC background studies. #2063

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
8 changes: 6 additions & 2 deletions MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -33,22 +33,26 @@ public:
mHadronPdgList = hadronPdgList;
mPartPdgToReplaceList = partPdgToReplaceList;
mFreqReplaceList = freqReplaceList;
// Ds1*(2700), Ds1*(2860), Ds3*(2860), Xic(3055)+, Xic(3080)+, Xic(3055)0, Xic(3080)0
mCustomPartPdgs = {30433, 40433, 437, 4315, 4316, 4325, 4326};
// Ds1*(2700), Ds1*(2860), Ds3*(2860), Xic(3055)+, Xic(3080)+, Xic(3055)0, Xic(3080)0, LambdaC(2595), LambdaC(2625)
mCustomPartPdgs = {30433, 40433, 437, 4315, 4316, 4325, 4326, 4124, 14122};
mCustomPartMasses[30433] = 2.714f;
mCustomPartMasses[40433] = 2.859f;
mCustomPartMasses[437] = 2.860f;
mCustomPartMasses[4315] = 3.0590f;
mCustomPartMasses[4316] = 3.0799f;
mCustomPartMasses[4325] = 3.0559f;
mCustomPartMasses[4326] = 3.0772f;
mCustomPartMasses[4124] = 2.62810f;
mCustomPartMasses[14122] = 2.59225f;
mCustomPartWidths[30433] = 0.122f;
mCustomPartWidths[40433] = 0.160f;
mCustomPartWidths[437] = 0.053f;
mCustomPartWidths[4315] = 0.0064f;
mCustomPartWidths[4316] = 0.0056f;
mCustomPartWidths[4325] = 0.0078f;
mCustomPartWidths[4326] = 0.0036f;
mCustomPartWidths[4124] = 0.00052f;
mCustomPartWidths[14122] = 0.0026f;
Print();
}

Expand Down
10 changes: 10 additions & 0 deletions MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkgSigmaC.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,10 @@

### The external generator derives from GeneratorPythia8.
[GeneratorExternal]
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C
funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5, -1.5, 1.5, {413}, {{413, 14122}, {413, 4124}}, {0.5, 0.5})

[GeneratorPythia8]
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_CorrBkgSigmaC.cfg
includePartonEvent=false
### not needed for jet studies, hence no need to keep parton event
84 changes: 84 additions & 0 deletions MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_CorrBkgSigmaC.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,84 @@
### authors: Mattia Faggin (mattia.faggin@cern.ch)
### last update: July 2025

### beams
Beams:idA 2212 # proton
Beams:idB 2212 # proton
Beams:eCM 13600. # GeV

### processes
SoftQCD:inelastic on # all inelastic processes

### decays
ParticleDecays:limitTau0 on
ParticleDecays:tau0Max 10.

### switching on Pythia Mode2
ColourReconnection:mode 1
ColourReconnection:allowDoubleJunRem off
ColourReconnection:m0 0.3
ColourReconnection:allowJunctions on
ColourReconnection:junctionCorrection 1.20
ColourReconnection:timeDilationMode 2
ColourReconnection:timeDilationPar 0.18
StringPT:sigma 0.335
StringZ:aLund 0.36
StringZ:bLund 0.56
StringFlav:probQQtoQ 0.078
StringFlav:ProbStoUD 0.2
StringFlav:probQQ1toQQ0join 0.0275,0.0275,0.0275,0.0275
MultiPartonInteractions:pT0Ref 2.15
BeamRemnants:remnantMode 1
BeamRemnants:saturation 5

# switch off charm-hadron decays
4122:onMode = off
14122:onMode = off ### Lambda_c(2595)+ already present in PYTHIA
4124:onMode = off ### Lambda_c(2625)+ already present in PYTHIA

## Λc decays
### Λc+ -> p K- π+
4122:oneChannel = 1 0.0350 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5%
4122:addChannel = 1 0.0196 100 2212 -313 ### Λc+ -> p K*0(892) 1.96%
4122:addChannel = 1 0.0108 100 2224 -321 ### Λc+ -> Delta++ K- 1.08%
4122:addChannel = 1 0.0022 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3
4122:onIfMatch = 2212 321 211
4122:onIfMatch = 2212 313
4122:onIfMatch = 2224 321
4122:onIfMatch = 102134 211
### Λc(2595)+
14122:oneChannel = 1 0.24 0 4222 -211 ### Λc(2595)+ -> Σc(2455)++ π- (24%)
14122:addChannel = 1 0.24 0 4112 211 ### Λc(2595)+ -> Σc(2455)0 π+ (24%)
14122:addChannel = 1 0.18 0 4122 211 -211 ### Λc(2595)+ -> Λc+ π+ π- (18%)
14122:onIfMatch = 4222 211 ### Λc(2595)+ -> Σc(2455)++ π-
14122:onIfMatch = 4112 211 ### Λc(2595)+ -> Σc(2455)0 π+
14122:onIfMatch = 4122 211 211 ### Λc(2595)+ -> Λc+ π+ π-
### Λc(2625)+
4124:oneChannel = 1 0.24 0 4222 -211 ### Λc(2592)+ -> Σc(2455)++ π- (24%)
4124:addChannel = 1 0.24 0 4112 211 ### Λc(2592)+ -> Σc(2455)0 π+ (24%)
4124:addChannel = 1 0.18 0 4122 211 -211 ### Λc(2592)+ -> Λc+ π+ π- (18%)
4124:onIfMatch = 4222 211 ### Λc(2625)+ -> Σc(2455)++ π-
4124:onIfMatch = 4112 211 ### Λc(2625)+ -> Σc(2455)0 π+
4124:onIfMatch = 4122 211 211 ### Λc(2625)+ -> Λc+ π+ π-

# Correct decay lengths (wrong in PYTHIA8 decay table)
# Lb
5122:tau0 = 0.4390
# Xic0
4132:tau0 = 0.0455
# OmegaC
4332:tau0 = 0.0803

# Correct Breit-Wigner width
4124:mWidth = 0.00052 # <0.52 MeV

# Force the decay of resonances in the decay chain
### K*0(892) -> K- π+
313:onMode = off
313:onIfAll = 321 211
### for Λc -> Delta++ K-
2224:onMode = off
2224:onIfAll = 2212 211
### for Λc -> Lambda(1520) K-
102134:onMode = off
102134:onIfAll = 2212 321