Skip to content

Commit 68247df

Browse files
author
Mattia Faggin
committed
Add configurations for MC production for SigmaC background studies.
1 parent 3448606 commit 68247df

File tree

3 files changed

+100
-2
lines changed

3 files changed

+100
-2
lines changed

‎MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C‎

Lines changed: 6 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -33,22 +33,26 @@ public:
3333
mHadronPdgList = hadronPdgList;
3434
mPartPdgToReplaceList = partPdgToReplaceList;
3535
mFreqReplaceList = freqReplaceList;
36-
// Ds1*(2700), Ds1*(2860), Ds3*(2860), Xic(3055)+, Xic(3080)+, Xic(3055)0, Xic(3080)0
37-
mCustomPartPdgs = {30433, 40433, 437, 4315, 4316, 4325, 4326};
36+
// Ds1*(2700), Ds1*(2860), Ds3*(2860), Xic(3055)+, Xic(3080)+, Xic(3055)0, Xic(3080)0, LambdaC(2595), LambdaC(2625)
37+
mCustomPartPdgs = {30433, 40433, 437, 4315, 4316, 4325, 4326, 4124, 14122};
3838
mCustomPartMasses[30433] = 2.714f;
3939
mCustomPartMasses[40433] = 2.859f;
4040
mCustomPartMasses[437] = 2.860f;
4141
mCustomPartMasses[4315] = 3.0590f;
4242
mCustomPartMasses[4316] = 3.0799f;
4343
mCustomPartMasses[4325] = 3.0559f;
4444
mCustomPartMasses[4326] = 3.0772f;
45+
mCustomPartMasses[4124] = 2.62810f;
46+
mCustomPartMasses[14122] = 2.59225f;
4547
mCustomPartWidths[30433] = 0.122f;
4648
mCustomPartWidths[40433] = 0.160f;
4749
mCustomPartWidths[437] = 0.053f;
4850
mCustomPartWidths[4315] = 0.0064f;
4951
mCustomPartWidths[4316] = 0.0056f;
5052
mCustomPartWidths[4325] = 0.0078f;
5153
mCustomPartWidths[4326] = 0.0036f;
54+
mCustomPartWidths[4124] = 0.00052f;
55+
mCustomPartWidths[14122] = 0.0026f;
5256
Print();
5357
}
5458

‎MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_Mode2_corrBkgSigmaC.ini‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,10 @@
1+
2+
### The external generator derives from GeneratorPythia8.
3+
[GeneratorExternal]
4+
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C
5+
funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5, -1.5, 1.5, {413}, {{413, 14122}, {413, 4124}}, {0.5, 0.5})
6+
7+
[GeneratorPythia8]
8+
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_CorrBkgSigmaC.cfg
9+
includePartonEvent=false
10+
### not needed for jet studies, hence no need to keep parton event

‎MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_CorrBkgSigmaC.cfg‎

Lines changed: 84 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,84 @@
1+
### authors: Mattia Faggin (mattia.faggin@cern.ch)
2+
### last update: July 2025
3+
4+
### beams
5+
Beams:idA 2212 # proton
6+
Beams:idB 2212 # proton
7+
Beams:eCM 13600. # GeV
8+
9+
### processes
10+
SoftQCD:inelastic on # all inelastic processes
11+
12+
### decays
13+
ParticleDecays:limitTau0 on
14+
ParticleDecays:tau0Max 10.
15+
16+
### switching on Pythia Mode2
17+
ColourReconnection:mode 1
18+
ColourReconnection:allowDoubleJunRem off
19+
ColourReconnection:m0 0.3
20+
ColourReconnection:allowJunctions on
21+
ColourReconnection:junctionCorrection 1.20
22+
ColourReconnection:timeDilationMode 2
23+
ColourReconnection:timeDilationPar 0.18
24+
StringPT:sigma 0.335
25+
StringZ:aLund 0.36
26+
StringZ:bLund 0.56
27+
StringFlav:probQQtoQ 0.078
28+
StringFlav:ProbStoUD 0.2
29+
StringFlav:probQQ1toQQ0join 0.0275,0.0275,0.0275,0.0275
30+
MultiPartonInteractions:pT0Ref 2.15
31+
BeamRemnants:remnantMode 1
32+
BeamRemnants:saturation 5
33+
34+
# switch off charm-hadron decays
35+
4122:onMode = off
36+
14122:onMode = off ### Lambda_c(2595)+ already present in PYTHIA
37+
4124:onMode = off ### Lambda_c(2625)+ already present in PYTHIA
38+
39+
## Λc decays
40+
### Λc+ -> p K- π+
41+
4122:oneChannel = 1 0.0350 0 2212 -321 211 ### Λc+ -> p K- π+ (non-resonant) 3.5%
42+
4122:addChannel = 1 0.0196 100 2212 -313 ### Λc+ -> p K*0(892) 1.96%
43+
4122:addChannel = 1 0.0108 100 2224 -321 ### Λc+ -> Delta++ K- 1.08%
44+
4122:addChannel = 1 0.0022 100 102134 211 ### Λc+ -> Lambda(1520) K- 2.20e-3
45+
4122:onIfMatch = 2212 321 211
46+
4122:onIfMatch = 2212 313
47+
4122:onIfMatch = 2224 321
48+
4122:onIfMatch = 102134 211
49+
### Λc(2595)+
50+
14122:oneChannel = 1 0.24 0 4222 -211 ### Λc(2595)+ -> Σc(2455)++ π- (24%)
51+
14122:addChannel = 1 0.24 0 4112 211 ### Λc(2595)+ -> Σc(2455)0 π+ (24%)
52+
14122:addChannel = 1 0.18 0 4122 211 -211 ### Λc(2595)+ -> Λc+ π+ π- (18%)
53+
14122:onIfMatch = 4222 211 ### Λc(2595)+ -> Σc(2455)++ π-
54+
14122:onIfMatch = 4112 211 ### Λc(2595)+ -> Σc(2455)0 π+
55+
14122:onIfMatch = 4122 211 211 ### Λc(2595)+ -> Λc+ π+ π-
56+
### Λc(2625)+
57+
4124:oneChannel = 1 0.24 0 4222 -211 ### Λc(2592)+ -> Σc(2455)++ π- (24%)
58+
4124:addChannel = 1 0.24 0 4112 211 ### Λc(2592)+ -> Σc(2455)0 π+ (24%)
59+
4124:addChannel = 1 0.18 0 4122 211 -211 ### Λc(2592)+ -> Λc+ π+ π- (18%)
60+
4124:onIfMatch = 4222 211 ### Λc(2625)+ -> Σc(2455)++ π-
61+
4124:onIfMatch = 4112 211 ### Λc(2625)+ -> Σc(2455)0 π+
62+
4124:onIfMatch = 4122 211 211 ### Λc(2625)+ -> Λc+ π+ π-
63+
64+
# Correct decay lengths (wrong in PYTHIA8 decay table)
65+
# Lb
66+
5122:tau0 = 0.4390
67+
# Xic0
68+
4132:tau0 = 0.0455
69+
# OmegaC
70+
4332:tau0 = 0.0803
71+
72+
# Correct Breit-Wigner width
73+
4124:mWidth = 0.00052 # <0.52 MeV
74+
75+
# Force the decay of resonances in the decay chain
76+
### K*0(892) -> K- π+
77+
313:onMode = off
78+
313:onIfAll = 321 211
79+
### for Λc -> Delta++ K-
80+
2224:onMode = off
81+
2224:onIfAll = 2212 211
82+
### for Λc -> Lambda(1520) K-
83+
102134:onMode = off
84+
102134:onIfAll = 2212 321

0 commit comments

Comments
 (0)