@@ -1544,9 +1544,9 @@ struct lambda1405analysis {
15441544
15451545 // Recompute the sigma momentum
15461546 bool success{false };
1547- float sigmaPRecalc = recalcSigmaMomReco (success, isSigmaMinusKink,
1548- sigmaCand.pxMoth (), sigmaCand.pyMoth (), sigmaCand.pzMoth (),
1549- sigmaCand.pxDaug (), sigmaCand.pyDaug (), sigmaCand.pzDaug ());
1547+ float sigmaPRecalc = recalcSigmaMom (success, isSigmaMinusKink,
1548+ sigmaCand.pxMoth (), sigmaCand.pyMoth (), sigmaCand.pzMoth (),
1549+ sigmaCand.pxDaug (), sigmaCand.pyDaug (), sigmaCand.pzDaug ());
15501550 if (!success && skipSigmasFailedRecompMom) {
15511551 return ;
15521552 }
@@ -1561,7 +1561,7 @@ struct lambda1405analysis {
15611561 sigmaMomReco[2 ] *= scale;
15621562
15631563 sigmaPt = std::sqrt (sigmaMomReco[0 ] * sigmaMomReco[0 ] + sigmaMomReco[1 ] * sigmaMomReco[1 ]);
1564- if (lambda1405Cand. isSigmaMinus ) {
1564+ if (isSigmaMinusKink ) {
15651565 rSigmaMinus.fill (HIST (" hDeltaPxRecalcSigmaMinus" ), sigmaMomReco[0 ] - sigmaCand.pxMoth (), sigmaCand.pxMoth ());
15661566 rSigmaMinus.fill (HIST (" hDeltaPyRecalcSigmaMinus" ), sigmaMomReco[1 ] - sigmaCand.pyMoth (), sigmaCand.pyMoth ());
15671567 rSigmaMinus.fill (HIST (" hDeltaPzRecalcSigmaMinus" ), sigmaMomReco[2 ] - sigmaCand.pzMoth (), sigmaCand.pzMoth ());
@@ -1589,35 +1589,35 @@ struct lambda1405analysis {
15891589 float genSigmaQtAP = qtAP (sigmaMomGen, kinkMomGen);
15901590 if (isSigmaMinusKink) {
15911591 rSigmaMinus.fill (HIST (" hSparseGenSigmaMinus" ), genMassSigma, genSigma.pt (), genSigmaAlphaAP, genSigmaQtAP, genCentMult, genNumContrib, genOcc);
1592- } else if (isSigmaPlusKink ) {
1592+ } else if (isSigmaPlusToPiKink || isSigmaPlusToPrKink ) {
15931593 rSigmaPlus.fill (HIST (" hSparseGenSigmaPlus" ), genMassSigma, genSigma.pt (), genSigmaAlphaAP, genSigmaQtAP, genCentMult, genNumContrib, genOcc);
15941594 }
15951595 }
15961596
15971597 // Fill table with sigma properties for efficiency studies
15981598 outputSigmaEffMC (
1599- sigmaCand.pxMoth , sigmaCand.pxMoth - genSigma.px (),
1600- sigmaCand.pyMoth , sigmaCand.pyMoth - genSigma.py (),
1601- sigmaCand.pzMoth , sigmaCand.pzMoth - genSigma.pz (),
1602- sigmaPt, sigmaPt - genSigma.pz (),
1599+ sigmaCand.pxMoth () , sigmaCand.pxMoth () - genSigma.px (),
1600+ sigmaCand.pyMoth () , sigmaCand.pyMoth () - genSigma.py (),
1601+ sigmaCand.pzMoth () , sigmaCand.pzMoth () - genSigma.pz (),
1602+ sigmaCand. pt (), sigmaCand. pt () - genSigma.pt (),
16031603 massSigma, massSigma - genMassSigma,
1604- sigmaCand.pxMoth - sigmaMomReco[0 ],
1605- sigmaCand.pyMoth - sigmaMomReco[1 ],
1606- sigmaCand.pzMoth - sigmaMomReco[2 ],
1604+ sigmaCand.pxMoth () - sigmaMomReco[0 ],
1605+ sigmaCand.pyMoth () - sigmaMomReco[1 ],
1606+ sigmaCand.pzMoth () - sigmaMomReco[2 ],
16071607 genSigma.phi (),
16081608 genSigma.eta (),
1609- sigmaCand.pxDaug , sigmaCand.pxDaug - genKinkDaug.px (),
1610- sigmaCand.pyDaug , sigmaCand.pyDaug - genKinkDaug.py (),
1611- sigmaCand.pzDaug , sigmaCand.pzDaug - genKinkDaug.pz (),
1612- kinkPt, kinkPt - genKinkDaug.pz (),
1609+ sigmaCand.pxDaug () , sigmaCand.pxDaug () - genKinkDaug.px (),
1610+ sigmaCand.pyDaug () , sigmaCand.pyDaug () - genKinkDaug.py (),
1611+ sigmaCand.pzDaug () , sigmaCand.pzDaug () - genKinkDaug.pz (),
1612+ kinkPt, kinkPt - genKinkDaug.pt (),
16131613 genKinkDaug.phi (),
16141614 genKinkDaug.eta (),
1615- sigmaCand.xDecVtx , sigmaCand.xDecVtx - genKinkDaug.vx ,
1616- sigmaCand.yDecVtx , sigmaCand.yDecVtx - genKinkDaug.vy ,
1617- sigmaCand.zDecVtx , sigmaCand.zDecVtx - genKinkDaug.vz ,
1615+ sigmaCand.xDecVtx () , sigmaCand.xDecVtx () - genKinkDaug.vx () ,
1616+ sigmaCand.yDecVtx () , sigmaCand.yDecVtx () - genKinkDaug.vy () ,
1617+ sigmaCand.zDecVtx () , sigmaCand.zDecVtx () - genKinkDaug.vz () ,
16181618 recoSigmaAlphaAP, recoSigmaQtAP,
1619- recoRecalcPtSigmaQtAP , recoRecalcPtSigmaQtAP,
1620- sigmaCand.dcaDaugPv , sigmaCand.dcaMothPv ,
1619+ recoRecalcPtSigmaAlphaAP , recoRecalcPtSigmaQtAP,
1620+ sigmaCand.dcaDaugPv () , sigmaCand.dcaMothPv () ,
16211621 genSigma.pdgCode (), genKinkDaug.pdgCode (),
16221622 centMult, collision.ft0cOccupancyInTimeRange (), collision.numContrib ());
16231623 }
0 commit comments