|
46 | 46 |
|
47 | 47 | #include <TF1.h> |
48 | 48 | #include <TH2.h> |
49 | | -#include <TH3.h> |
50 | 49 | #include <THnSparse.h> |
| 50 | +#include <TPDGCode.h> |
| 51 | +#include <TString.h> |
51 | 52 |
|
52 | 53 | #include <array> |
53 | 54 | #include <cmath> |
| 55 | +#include <cstdint> |
54 | 56 | #include <numeric> |
55 | 57 | #include <string> |
56 | 58 | #include <vector> |
@@ -81,23 +83,23 @@ struct lambda1405candidate { |
81 | 83 | float phi = -1; // Phi of the Lambda(1405) candidate |
82 | 84 | float eta = -1; // Eta of the Lambda(1405) candidate |
83 | 85 |
|
84 | | - bool isSigmaPlus = false; // True if compatible with Sigma+ |
85 | | - bool isSigmaMinus = false; // True if compatible with Sigma- |
86 | | - float sigmaMinusMass = -1; // Invariant mass of the Sigma- candidate |
87 | | - float sigmaPlusMass = -1; // Invariant mass of the Sigma+ candidate |
88 | | - float xiMinusMass = -1; // Invariant mass of the Xi- candidate |
89 | | - int sigmaSign = 0; // Sign of the Sigma candidate: 1 for matter, -1 for antimatter |
90 | | - float sigmaPt = -1; // pT of the Sigma daughter |
91 | | - float sigmaAlphaAP = -1; // Alpha of the Sigma |
92 | | - float sigmaQtAP = -1; // qT of the Sigma |
93 | | - float sigmaRadius = -1; // Radius of the Sigma decay vertex |
94 | | - float dcaSigmaToPv = -1; // DCA of the Sigma candidate to the primary vertex |
95 | | - float kinkPt = -1; // pT of the kink daughter |
96 | | - float kinkPiNSigTpc = -1; // Number of sigmas for the pion candidate from Sigma kink in Tpc |
97 | | - float kinkPiNSigTof = -1; // Number of sigmas for the pion candidate from Sigma kink in Tof |
98 | | - float kinkPrNSigTpc = -1; // Number of sigmas for the proton candidate from Sigma kink in Tpc |
99 | | - float kinkPrNSigTof = -1; // Number of sigmas for the proton candidate from Sigma kink in Tof |
100 | | - float kinkDcaDauToPv = -1; // DCA of the kink daughter to the primary vertex |
| 86 | + bool isSigmaPlus = false; // True if compatible with Sigma+ |
| 87 | + bool isSigmaMinus = false; // True if compatible with Sigma- |
| 88 | + float sigmaMinusMass = -1; // Invariant mass of the Sigma- candidate |
| 89 | + float sigmaPlusMass = -1; // Invariant mass of the Sigma+ candidate |
| 90 | + float xiMinusMass = -1; // Invariant mass of the Xi- candidate |
| 91 | + int sigmaSign = 0; // Sign of the Sigma candidate: 1 for matter, -1 for antimatter |
| 92 | + float sigmaPt = -1; // pT of the Sigma daughter |
| 93 | + float sigmaAlphaAP = -1; // Alpha of the Sigma |
| 94 | + float sigmaQtAP = -1; // qT of the Sigma |
| 95 | + float sigmaRadius = -1; // Radius of the Sigma decay vertex |
| 96 | + float dcaSigmaToPv = -1; // DCA of the Sigma candidate to the primary vertex |
| 97 | + float kinkPt = -1; // pT of the kink daughter |
| 98 | + float kinkPiNSigTpc = -1; // Number of sigmas for the pion candidate from Sigma kink in Tpc |
| 99 | + float kinkPiNSigTof = -1; // Number of sigmas for the pion candidate from Sigma kink in Tof |
| 100 | + float kinkPrNSigTpc = -1; // Number of sigmas for the proton candidate from Sigma kink in Tpc |
| 101 | + float kinkPrNSigTof = -1; // Number of sigmas for the proton candidate from Sigma kink in Tof |
| 102 | + float kinkDcaDauToPv = -1; // DCA of the kink daughter to the primary vertex |
101 | 103 |
|
102 | 104 | float bachPiPt = -1; // pT of the pion daughter |
103 | 105 | float bachPiNSigTpc = -1; // Number of sigmas for the pion candidate |
@@ -1011,7 +1013,7 @@ struct lambda1405analysis { |
1011 | 1013 | if (!selectPiBach(piTrack, primaryVertex, lambda1405Cand.centMult)) { |
1012 | 1014 | continue; |
1013 | 1015 | } |
1014 | | - rSelections.fill(HIST("hSelectionsL1405"), 2); // Bach Pi selection |
| 1016 | + rSelections.fill(HIST("hSelectionsL1405"), 2); // Bach Pi selection |
1015 | 1017 |
|
1016 | 1018 | auto piMom = std::array{piTrack.px(), piTrack.py(), piTrack.pz()}; |
1017 | 1019 | float invMass{-1.f}; |
@@ -1337,7 +1339,7 @@ struct lambda1405analysis { |
1337 | 1339 | auto labelBachPi = trackLabelsMC.rawIteratorAt(lambda1405Cand.bachPiId); |
1338 | 1340 | if (!labelBachPi.has_mcParticle()) { |
1339 | 1341 | rSelections.fill(HIST("hRecoNotMatchedCounter"), 3); // Bach pion not matched |
1340 | | - continue; // Skip if no valid MC association |
| 1342 | + continue; // Skip if no valid MC association |
1341 | 1343 | } |
1342 | 1344 |
|
1343 | 1345 | auto genBachPi = labelBachPi.template mcParticle_as<aod::McParticles>(); |
|
0 commit comments