Skip to content

Commit 98d2c7c

Browse files
fgrosafgrosa
committed
Add Xic0 excited states
1 parent 1a62026 commit 98d2c7c

File tree

4 files changed

+49
-17
lines changed

4 files changed

+49
-17
lines changed

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

Lines changed: 10 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -32,20 +32,24 @@ public:
3232
mHadronPdgList = hadronPdgList;
3333
mPartPdgToReplaceList = partPdgToReplaceList;
3434
mFreqReplaceList = freqReplaceList;
35-
// Ds1*(2700), Ds1*(2860), Ds3*(2860), Xic(3055), Xic(3080)
36-
mCustomPartPdgs = {30433, 40433, 437, 4325, 4326};
35+
// Ds1*(2700), Ds1*(2860), Ds3*(2860), Xic(3055)+, Xic(3080)+, Xic(3055)0, Xic(3080)0
36+
mCustomPartPdgs = {30433, 40433, 437, 4315, 4316, 4325, 4326};
3737
mCustomPartMasses[30433] = 2.714f;
3838
mCustomPartMasses[40433] = 2.859f;
3939
mCustomPartMasses[437] = 2.860f;
40+
mCustomPartMasses[4315] = 3.0590f;
41+
mCustomPartMasses[4316] = 3.0799f;
4042
mCustomPartMasses[4325] = 3.0559f;
4143
mCustomPartMasses[4326] = 3.0772f;
4244
mCustomPartWidths[30433] = 0.122f;
4345
mCustomPartWidths[40433] = 0.160f;
4446
mCustomPartWidths[437] = 0.053f;
47+
mCustomPartWidths[4315] = 0.0064f;
48+
mCustomPartWidths[4316] = 0.0056f;
4549
mCustomPartWidths[4325] = 0.0078f;
4650
mCustomPartWidths[4326] = 0.0036f;
4751
Print();
48-
}
52+
}
4953

5054
/// Destructor
5155
~GeneratorPythia8GapTriggeredHF() = default;
@@ -285,10 +289,12 @@ protected:
285289
}
286290
float energy = std::sqrt(px*px + py*py + pz*pz + mass*mass);
287291

292+
// we remove particle to replace and its daughters
293+
mPythia.event[iPartToReplace].undoDecay();
288294
mPythia.event.remove(iPartToReplace, iPartToReplace, true); // we remove the original particle
295+
289296
mPythia.event.append(charge * pdgCodeNew, status, mothers[0], mothers[1], 0, 0, 0, 0, px, py, pz, energy, mass);
290297
mPythia.moreDecays(); // we need to decay the new particle
291-
292298
return true;
293299
}
294300

‎MC/config/PWGHF/ini/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.ini‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
### The external generator derives from GeneratorPythia8.
22
[GeneratorExternal]
33
fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_pythia8_gaptriggered_hf.C
4-
funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5, -1.5, 1.5, {10433, 435}, {{10433, 30433}, {10433, 437}, {435, 4325}, {435, 4326}}, {0.1, 0.1, 0.1, 0.1})
4+
funcName=GeneratorPythia8GapTriggeredCharmAndBeauty(5, -1.5, 1.5, -1.5, 1.5, {10433, 435}, {{10433, 30433}, {10433, 437}, {435, 4325}, {435, 4326}, {425, 4315}, {425, 4316}}, {0.1, 0.1, 0.1, 0.1, 0.5, 0.5})
55

66
[GeneratorPythia8]
77
config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_DResoTrigger.cfg

‎MC/config/PWGHF/ini/tests/GeneratorHF_D2H_ccbar_and_bbbar_gap5_DResoTrigger.C‎

Lines changed: 9 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -4,19 +4,22 @@ int External() {
44
int checkPdgQuarkOne{4};
55
int checkPdgQuarkTwo{5};
66
float ratioTrigger = 1./5.; // one event triggered out of 5
7-
std::array<std::array<int, 2>, 4> pdgReplParticles = {{10433, 30433}, {10433, 437}, {435, 4325}, {435, 4326}};
8-
std::array<std::array<int, 2>, 4> pdgReplPartCounters = {{0, 0}, {0, 0}, {0, 0}, {0, 0}};
9-
std::array<float, 4> freqRepl = {0.1, 0.1, 0.1, 0.1}; // one event triggered out of 5
7+
std::array<std::array<int, 2>, 6> pdgReplParticles = {{10433, 30433}, {10433, 437}, {435, 4325}, {435, 4326}, {425, 4315}, {425, 4316}};
8+
std::array<std::array<int, 2>, 6> pdgReplPartCounters = {{0, 0}, {0, 0}, {0, 0}, {0, 0}, {0, 0}, {0, 0}};
9+
std::array<float, 4> freqRepl = {0.1, 0.1, 0.1, 0.1, 0.5, 0.5}; // one event triggered out of 5
1010

11-
std::vector<int> checkPdgHadron{10433, 30433, 435, 437, 4325, 4326};
11+
std::vector<int> checkPdgHadron{10433, 30433, 435, 437, 4325, 4326, 4315, 4316};
1212
std::map<int, std::vector<std::vector<int>>> checkHadronDecays{ // sorted pdg of daughters
1313
{411, {{-321, 211, 211}, {-313, 211}, {211, 311}, {211, 333}}}, // D+
14+
{421, {{-321, 211}, {-321, 211, 111}}}, // D0
1415
{435, {{311, 413}, {311, 411}}}, // Ds2*(2573)
1516
{10433, {{311, 413}}}, // Ds1(2536)
1617
{30433, {{311, 413}}}, // Ds1*(2700)
1718
{437, {{311, 413}}}, // Ds3*(2860)
18-
{4325, {{411, 3122}}}, // Xic(3055)
19-
{4326, {{411, 3122}}}, // Xic(3080)
19+
{4325, {{411, 3122}}}, // Xic(3055)+
20+
{4326, {{411, 3122}}}, // Xic(3080)+
21+
{4315, {{421, 3122}}}, // Xic(3055)+
22+
{4316, {{421, 3122}}}, // Xic(3080)+
2023
{531, {{-435, -11, 12}, {-10433, -11, 12}, {-435, -13, 14}, {-10433, -13, 14}, {-435, -15, 16}, {-10433, -15, 16}, {-435, 211}}}// Bs0
2124
};
2225

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

Lines changed: 29 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -31,22 +31,32 @@ StringFlav:mesonCL1S1J2= 5
3131
437:all = D*_s3(2860)+ D*_s3(2860)- 7 3 0 2.860 0.053 2.51 10. 0.
3232
437:mayDecay = on
3333

34-
### Xic(3055)
34+
### Xic(3055)+
3535
4325:all = Xi_c(3055)+ Xi_cbar(3055)- 4 3 0 3.0559 0.0078 2.986 10. 0.
3636
4325:mayDecay = on
3737

38-
### Xic(3080)
38+
### Xic(3080)+
3939
4326:all = Xi_c(3080)+ Xi_cbar(3080)- 2 3 0 3.0772 0.0036 2.986 10. 0.
4040
4326:mayDecay = on
4141

42+
### Xic(3055)0
43+
4315:all = Xi_c(3055)0 Xi_cbar(3055)0 4 0 0 3.0590 0.0064 2.981 10. 0.
44+
4315:mayDecay = on
45+
46+
### Xic(3080)0
47+
4316:all = Xi_c(3080)0 Xi_cbar(3080)0 2 0 0 3.0799 0.0056 2.981 10. 0.
48+
4316:mayDecay = on
49+
4250
### turn off all charm resonances decays
4351
10433:onMode = off # Ds1(2536)
4452
435:onMode = off # Ds2*(2573)
4553
30433:onMode = off # Ds1*(2700)
4654
40433:onMode = off # Ds1*(2860)
4755
437:onMode = off # Ds3*(2860)
48-
4325:onMode = off # Xic(3055)
49-
4326:onMode = off # Xic(3080)
56+
4325:onMode = off # Xic(3055)+
57+
4326:onMode = off # Xic(3080)+
58+
4315:onMode = off # Xic(3055)0
59+
4316:onMode = off # Xic(3080)0
5060

5161
### turn off all beauty hadron decays
5262
531:onMode = off
@@ -74,14 +84,22 @@ StringFlav:mesonCL1S1J2= 5
7484
437:oneChannel = 1 1 0 413 311
7585
437:onIfMatch = 413 311
7686

77-
### add Xic(3055)
87+
### add Xic(3055)+
7888
4325:oneChannel = 1 1 0 411 3122
7989
4325:onIfMatch = 411 3122
8090

81-
### add Xic(3080)
91+
### add Xic(3080)+
8292
4326:oneChannel = 1 1 0 411 3122
8393
4326:onIfMatch = 411 3122
8494

95+
### add Xic(3055)0
96+
4315:oneChannel = 1 1 0 421 3122
97+
4315:onIfMatch = 421 3122
98+
99+
### add Xic(3080)0
100+
4316:oneChannel = 1 1 0 421 3122
101+
4316:onIfMatch = 421 3122
102+
85103
### add Bs0
86104
531:oneChannel = 1 0.0070000 0 12 -11 -435
87105
531:addChannel = 1 0.0070000 0 12 -11 -10433
@@ -113,6 +131,9 @@ StringFlav:mesonCL1S1J2= 5
113131
411:addChannel = 1 0.0104 0 -313 211
114132
411:addChannel = 1 0.0156 0 311 211
115133
411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi
134+
### set D0 BRs
135+
421:oneChannel = 1 0.9 0 -321 211
136+
421:addChannel = 1 0.1 0 -321 211 111
116137

117138
### K* -> K pi
118139
313:onMode = off
@@ -127,6 +148,8 @@ StringFlav:mesonCL1S1J2= 5
127148

128149
### D0 -> K pi
129150
421:onIfMatch = 321 211
151+
### D0 -> K pi pi0
152+
421:onIfMatch = 321 211 111
130153

131154
### D+/- -> K pi pi
132155
411:onIfMatch = 321 211 211

0 commit comments

Comments
 (0)