Skip to content

Commit fd508bf

Browse files
fgrosafgrosa
authored
Fix HF tests (#2208)
* Fix HF tests * Fix another test file * Revert wrong change * Few more fixes * Even more fixes * Fix typo and slightly lower threshold
1 parent 22da986 commit fd508bf

6 files changed

+121
-46
lines changed

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

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -68,7 +68,7 @@ int External() {
6868
for (int j{track.getFirstDaughterTrackId()}; j <= track.getLastDaughterTrackId(); ++j) {
6969
auto pdgDau = tracks->at(j).GetPdgCode();
7070
pdgsDecay.push_back(pdgDau);
71-
if (pdgDau != 333) { // phi is antiparticle of itself
71+
if (pdgDau != 333 && pdgDau != 111) { // phi is antiparticle of itself
7272
pdgsDecayAntiPart.push_back(-pdgDau);
7373
} else {
7474
pdgsDecayAntiPart.push_back(pdgDau);
@@ -129,7 +129,7 @@ int External() {
129129
return 1;
130130
}
131131

132-
if (averagePt < 8.) { // by testing locally it should be around 8.5 GeV/c with pthard bin 20-200 (contrary to 2-2.5 GeV/c of SoftQCD)
132+
if (averagePt < 7.) { // by testing locally it should be around 8.5 GeV/c with pthard bin 20-200 (contrary to 2-2.5 GeV/c of SoftQCD)
133133
std::cerr << "Average pT of charmed hadrons " << averagePt << " lower than expected\n";
134134
return 1;
135135
}

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

Lines changed: 61 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -7,13 +7,65 @@ int External() {
77

88
std::vector<int> checkPdgHadron{411, 421, 431, 4122, 4132, 4232, 4332};
99
std::map<int, std::vector<std::vector<int>>> checkHadronDecays{ // sorted pdg of daughters
10-
{411, {{-321, 211, 211}, {-313, 211}, {211, 311}, {211, 333}}}, // D+
11-
{421, {{-321, 211}, {-321, 111, 211}}}, // D0
12-
{431, {{211, 333}, {-313, 321}}}, // Ds+
13-
{4122, {{-313, 2212}, {-321, 2224}, {211, 102134}, {-321, 211, 2212}, {311, 2212}}}, // Lc+
14-
{4132, {{211, 3312}}}, // Xic0
15-
{4232, {{-313, 2212}, {-321, 3324}, {211, 211, 3312}, {-321, 211, 2212}}}, // Xic+
16-
{4332, {{211, 3334}}} // Omegac+
10+
{411, {
11+
{-321, 211, 211}, // K- π+ π+ (non-resonant)
12+
{-321, 111, 211, 211}, // K- π+ π+ π0 (non-resonant)
13+
{-313, 321}, // K*0(892) K+
14+
{-10311, 321}, // K*0(1430) K+
15+
{211, 333}, // φ π+
16+
{-321, 211, 321}, // K- K+ π+ (non-resonant)
17+
{113, 211}, // ρ0 π+
18+
{211, 225}, // f2(1270) π+
19+
{-211, 211, 211} // π- π+ π+ (non-resonant)
20+
}},
21+
{421, {
22+
{-321, 211}, // K- π+ (non-resonant)
23+
{-321, 111, 211}, // K- π+ π0
24+
{-321, 213}, // ρ+ K-
25+
{-313, 111}, // antiK*0(892) π0
26+
{-323, 211}, // K*-(892) π+
27+
{-211, 211}, // π- π+
28+
{-211, 213}, // ρ+ π-
29+
{-211, 111, 211}, // π- π+ π0
30+
{-321, 321} // K- K+
31+
}},
32+
{431, {
33+
{211, 333}, // φ π+
34+
{-313, 321}, // antiK*(892) K+
35+
{213, 333}, // φ ρ
36+
{113, 211}, // ρ π+
37+
{211, 225}, // f2(1270) π+
38+
{-211, 211, 211}, // π- π+ π+ (s-wave)
39+
{211, 313}, // K*(892)0 π+
40+
{321, 10221}, // f0(1370) K+
41+
{113, 321}, // ρ0 K+
42+
{-211, 211, 321}, // π- K+ π+ (non-resonant)
43+
{211, 221} // η π+
44+
}},
45+
{4122, {
46+
{-321, 211, 2212}, // p K- π+ (non-resonant)
47+
{-313, 2212}, // p K*0(892)
48+
{-321, 2224}, // Δ++ K-
49+
{211, 102134}, // Λ(1520) π+
50+
{-321, 111, 211, 2212}, // p K- π+ π0
51+
{-211, 211, 2212}, // p π- π+
52+
{333, 2212} // p φ
53+
}},
54+
{4232, {
55+
{-321, 211, 2212}, // Xic+ -> p, K-, pi+
56+
{-313, 2212}, // Xic+ -> p, Kbar^*(892)0
57+
{211, 211, 3312}, // Xic+ -> Xi-, pi+, pi+
58+
{333, 2212}, // Xic+ -> p, phi(1020)0
59+
{-211, 211, 3222}, // Xic+ -> Sigma+, pi-, pi+
60+
{211, 3324} // Xic+ -> Xi(1530)0, pi+
61+
}},
62+
{4132, {
63+
{211, 3312}, // Xic0 -> Xi-, pi+
64+
}},
65+
{4332, {
66+
{211, 3334}, // Omegac0 -> Omega-, pi+
67+
{211, 3312} // Omegac0 -> Xi-, pi+
68+
}}
1769
};
1870

1971
TFile file(path.c_str(), "READ");
@@ -66,7 +118,7 @@ int External() {
66118
for (int j{track.getFirstDaughterTrackId()}; j <= track.getLastDaughterTrackId(); ++j) {
67119
auto pdgDau = tracks->at(j).GetPdgCode();
68120
pdgsDecay.push_back(pdgDau);
69-
if (pdgDau != 333) { // phi is antiparticle of itself
121+
if (pdgDau != 333 && pdgDau != 111 && pdgDau != 221 && pdgDau != 113 && pdgDau != 225) { // phi is antiparticle of itself
70122
pdgsDecayAntiPart.push_back(-pdgDau);
71123
} else {
72124
pdgsDecayAntiPart.push_back(pdgDau);
@@ -119,7 +171,7 @@ int External() {
119171
}
120172

121173
float fracForcedDecays = float(nSignalGoodDecay) / nSignals;
122-
if (fracForcedDecays < 0.9) { // we put some tolerance (e.g. due to oscillations which might change the final state)
174+
if (fracForcedDecays < 0.85) { // we put some tolerance (e.g. due to oscillations which might change the final state)
123175
std::cerr << "Fraction of signals decaying into the correct channel " << fracForcedDecays << " lower than expected\n";
124176
return 1;
125177
}

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

Lines changed: 18 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -5,11 +5,25 @@ int External() {
55
int checkPdgQuarkTwo{5};
66
float ratioTrigger = 1./5; // one event triggered out of 5
77

8-
std::vector<int> checkPdgHadron{4132, 4332, 4232};
8+
std::vector<int> checkPdgHadron{4132, 4232, 4332};
99
std::map<int, std::vector<std::vector<int>>> checkHadronDecays{ // sorted pdg of daughters
10-
{4132, {{211, 3312}}}, // Xic0
11-
{4332, {{211, 3334}, {211, 3312}}}, // Omegac0
12-
{4232, {{313, 2212}, {211, 321, 2212}, {211, 3324}, {211, 211, 3312}}} //Xic+
10+
{4232, {
11+
{-321, 211, 2212}, // Xic+ -> p, K-, pi+
12+
{-313, 2212}, // Xic+ -> p, Kbar^*(892)0
13+
{211, 211, 3312}, // Xic+ -> Xi-, pi+, pi+
14+
{333, 2212}, // Xic+ -> p, phi(1020)0
15+
{-211, 211, 3222}, // Xic+ -> Sigma+, pi-, pi+
16+
{211, 3324}, // Xic+ -> Xi(1530)0, pi+
17+
}},
18+
19+
{4132, {
20+
{211, 3312}, // Xic0 -> Xi-, pi+
21+
}},
22+
23+
{4332, {
24+
{211, 3334}, // Omegac0 -> Omega-, pi+
25+
{211, 3312}, // Omegac0 -> Xi-, pi+
26+
}},
1327
};
1428

1529
TFile file(path.c_str(), "READ");

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

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -57,7 +57,6 @@ int External() {
5757
{2212, -211, 211}, // Lambdac+ -> p, pi-, pi+
5858
{2212, 333}, // Lambdac+ -> p, phi(1020)0
5959
}},
60-
6160
{4232, {
6261
{2212, -321, 211}, // Xic+ -> p, K-, pi+
6362
{2212, -313}, // Xic+ -> p, Kbar^*(892)0

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

Lines changed: 38 additions & 28 deletions
Original file line numberDiff line numberDiff line change
@@ -5,57 +5,67 @@ int External() {
55
int checkPdgQuarkTwo{5};
66
float ratioTrigger = 1./5; // one event triggered out of 5
77

8-
std::vector<int> checkPdgHadron{411, 421, 431, 4122, 4232};
8+
std::vector<int> checkPdgHadron{411, 421, 431, 4122, 4132, 4232, 4332};
99
std::map<int, std::vector<std::vector<int>>> checkHadronDecays{ // sorted pdg of daughters
1010
{411, {
1111
{-321, 211, 211}, // K- π+ π+ (non-resonant)
12+
{-321, 111, 211, 211}, // K- π+ π+ π0 (non-resonant)
1213
{-313, 321}, // K*0(892) K+
1314
{-10311, 321}, // K*0(1430) K+
1415
{211, 333}, // φ π+
15-
{-321, 321, 211}, // K- K+ π+ (non-resonant)
16+
{-321, 211, 321}, // K- K+ π+ (non-resonant)
1617
{113, 211}, // ρ0 π+
17-
{225, 211}, // f2(1270) π+
18-
{-211, 211, 211} // π- π+ π+ (non-resonant)
18+
{211, 225}, // f2(1270) π+
19+
{-211, 211, 211} // π- π+ π+ (non-resonant)
1920
}},
2021
{421, {
2122
{-321, 211}, // K- π+ (non-resonant)
2223
{-321, 111, 211}, // K- π+ π0
23-
{213, -321}, // ρ+ K-
24+
{-321, 213}, // ρ+ K-
2425
{-313, 111}, // antiK*0(892) π0
2526
{-323, 211}, // K*-(892) π+
2627
{-211, 211}, // π- π+
27-
{213, -211}, // ρ+ π-
28-
{-211, 211, 111}, // π- π+ π0
29-
{-321, 321} // K- K+
28+
{-211, 213}, // ρ+ π-
29+
{-211, 111, 211}, // π- π+ π0
30+
{-321, 321} // K- K+
3031
}},
3132
{431, {
3233
{211, 333}, // φ π+
3334
{-313, 321}, // antiK*(892) K+
34-
{333, 213}, // φ ρ
35+
{213, 333}, // φ ρ
3536
{113, 211}, // ρ π+
36-
{225, 211}, // f2(1270) π+
37+
{211, 225}, // f2(1270) π+
3738
{-211, 211, 211}, // π- π+ π+ (s-wave)
38-
{313, 211}, // K*(892)0 π+
39-
{10221, 321}, // f0(1370) K+
39+
{211, 313}, // K*(892)0 π+
40+
{321, 10221}, // f0(1370) K+
4041
{113, 321}, // ρ0 K+
41-
{-211, 321, 211}, // π- K+ π+ (non-resonant)
42-
{221, 211} // η π+
42+
{-211, 211, 321}, // π- K+ π+ (non-resonant)
43+
{211, 221} // η π+
4344
}},
4445
{4122, {
45-
{2212, -321, 211}, // p K- π+ (non-resonant)
46-
{2212, -313}, // p K*0(892)
47-
{2224, -321}, // Δ++ K-
48-
{102134, 211}, // Λ(1520) K-
49-
{2212, -321, 211, 111}, // p K- π+ π0
50-
{2212, -211, 211}, // p π- π+
51-
{2212, 333} // p φ
52-
}},
53-
{4232, {
54-
{-313, 2212}, // antiK*0(892) p
55-
{2212, -321, 211}, // p K- π+
56-
{2212, 333}, // p φ
57-
{3222, -211, 211} // Σ+ π- π+
46+
{-321, 211, 2212}, // p K- π+ (non-resonant)
47+
{-313, 2212}, // p K*0(892)
48+
{-321, 2224}, // Δ++ K-
49+
{211, 102134}, // Λ(1520) π+
50+
{-321, 111, 211, 2212}, // p K- π+ π0
51+
{-211, 211, 2212}, // p π- π+
52+
{333, 2212} // p φ
5853
}},
54+
{4232, {
55+
{-321, 211, 2212}, // Xic+ -> p, K-, pi+
56+
{-313, 2212}, // Xic+ -> p, Kbar^*(892)0
57+
{211, 211, 3312}, // Xic+ -> Xi-, pi+, pi+
58+
{333, 2212}, // Xic+ -> p, phi(1020)0
59+
{-211, 211, 3222}, // Xic+ -> Sigma+, pi-, pi+
60+
{211, 3324} // Xic+ -> Xi(1530)0, pi+
61+
}},
62+
{4132, {
63+
{211, 3312}, // Xic0 -> Xi-, pi+
64+
}},
65+
{4332, {
66+
{211, 3334}, // Omegac0 -> Omega-, pi+
67+
{211, 3312} // Omegac0 -> Xi-, pi+
68+
}}
5969
};
6070

6171
TFile file(path.c_str(), "READ");
@@ -108,7 +118,7 @@ int External() {
108118
for (int j{track.getFirstDaughterTrackId()}; j <= track.getLastDaughterTrackId(); ++j) {
109119
auto pdgDau = tracks->at(j).GetPdgCode();
110120
pdgsDecay.push_back(pdgDau);
111-
if (pdgDau != 333) { // phi is antiparticle of itself
121+
if (pdgDau != 333 && pdgDau != 111 && pdgDau != 221 && pdgDau != 113 && pdgDau != 225) { // phi is antiparticle of itself
112122
pdgsDecayAntiPart.push_back(-pdgDau);
113123
} else {
114124
pdgsDecayAntiPart.push_back(pdgDau);

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

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -5,7 +5,7 @@
55
#include <iostream>
66

77
int External() {
8-
std::string path{"/home/mattia/Documenti/cernbox/Documents/PostDoc/D2H/MC/corrBkgSigmaC/tf1/genevents_Kine.root"};
8+
std::string path{"o2sim_Kine.root"};
99

1010
int checkPdgQuarkOne{4};
1111
int checkPdgQuarkTwo{5};
@@ -177,7 +177,7 @@ int External() {
177177
}
178178

179179
float fracForcedDecays = float(nSignalGoodDecay) / nSignals;
180-
if (fracForcedDecays < 0.9) { // we put some tolerance (e.g. due to oscillations which might change the final state)
180+
if (fracForcedDecays < 0.5) { // we put some tolerance (e.g. due to oscillations which might change the final state)
181181
std::cerr << "Fraction of signals decaying into the correct channel " << fracForcedDecays << " lower than expected\n";
182182
return 1;
183183
}

0 commit comments

Comments
 (0)