Skip to content

Commit 5627e5f

Browse files
committed
Corrected error in accessing McCollision indices
1 parent dd8c138 commit 5627e5f

File tree

1 file changed

+6
-4
lines changed

1 file changed

+6
-4
lines changed

PWGLF/TableProducer/Strangeness/sigmaminustask.cxx

Lines changed: 6 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -114,7 +114,7 @@ struct sigmaminustask {
114114
}
115115
PROCESS_SWITCH(sigmaminustask, processData, "Data processing", true);
116116

117-
void processMC(CollisionsFullMC const& collisions, aod::KinkCands const& KinkCands, aod::McTrackLabels const& trackLabelsMC, aod::McParticles const& particlesMC, TracksFull const&)
117+
void processMC(CollisionsFullMC const& collisions, aod::KinkCands const& KinkCands, aod::McTrackLabels const& trackLabelsMC, aod::McParticles const& particlesMC, aod::McCollisions const&, TracksFull const&)
118118
{
119119
for (const auto& collision : collisions) {
120120
if (std::abs(collision.posZ()) > cutzvertex || !collision.sel8()) {
@@ -167,9 +167,11 @@ struct sigmaminustask {
167167
float decayRadiusMC = std::sqrt(deltaXMother * deltaXMother + deltaYMother * deltaYMother);
168168

169169
// Check coherence of MCcollision Id for daughter MCparticle and reconstructed collision
170-
auto mcCollision = mcTrackPiDau.template mcCollision_as<aod::McCollisions>();
171-
bool mcCollisionIdCheck = collision.mcCollisionId() == mcCollision.globalIndex();
172-
170+
bool mcCollisionIdCheck = false;
171+
if (collision.has_mcCollision()) {
172+
mcCollisionIdCheck = collision.mcCollision().globalIndex() == mcTrackPiDau.mcCollisionId();
173+
}
174+
173175
rSigmaMinus.fill(HIST("h2MassPtMCRec"), kinkCand.mothSign() * kinkCand.ptMoth(), kinkCand.mSigmaMinus());
174176
if (mcTrackSigma.pdgCode() > 0) {
175177
rSigmaMinus.fill(HIST("h2MassResolution_plus"), kinkCand.mSigmaMinus(), kinkCand.mSigmaMinus() - MotherMassMC);

0 commit comments

Comments
 (0)