@@ -5,7 +5,7 @@ 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)
@@ -21,13 +21,13 @@ int External() {
2121 {421, {
2222 {-321, 211}, // K- π+ (non-resonant)
2323 {-321, 111, 211}, // K- π+ π0
24- {213, -321}, // ρ+ K-
24+ {-321, 213 }, // ρ+ K-
2525 {-313, 111}, // antiK*0(892) π0
2626 {-323, 211}, // K*-(892) π+
2727 {-211, 211}, // π- π+
28- {213, -211}, // ρ+ π-
29- {-211, 211, 111 }, // π- π+ π0
30- {-321, 321} // K- K+
28+ {-211, 213 }, // ρ+ π-
29+ {-211, 111, 211 }, // π- π+ π0
30+ {-321, 321} // K- K+
3131 }},
3232 {431, {
3333 {211, 333}, // φ π+
@@ -46,17 +46,26 @@ int External() {
4646 {-321, 211, 2212}, // p K- π+ (non-resonant)
4747 {-313, 2212}, // p K*0(892)
4848 {-321, 2224}, // Δ++ K-
49- {211, 102134}, // Λ(1520) K-
49+ {211, 102134}, // Λ(1520) π+
5050 {-321, 111, 211, 2212}, // p K- π+ π0
5151 {-211, 211, 2212}, // p π- π+
5252 {333, 2212} // p φ
5353 }},
54- {4232, {
55- {-313, 2212}, // antiK*0(892) p
56- {-321, 211, 2212}, // p K- π+
57- {333, 2212}, // p φ
58- {-211, 211, 3222} // Σ+ π- π+
59- }}
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+ }}
6069 };
6170
6271 TFile file(path.c_str(), "READ");
@@ -109,7 +118,7 @@ int External() {
109118 for (int j{track.getFirstDaughterTrackId()}; j <= track.getLastDaughterTrackId(); ++j) {
110119 auto pdgDau = tracks->at(j).GetPdgCode();
111120 pdgsDecay.push_back(pdgDau);
112- 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
113122 pdgsDecayAntiPart.push_back(-pdgDau);
114123 } else {
115124 pdgsDecayAntiPart.push_back(pdgDau);
0 commit comments