@@ -3502,18 +3502,18 @@ struct HfTrackIndexSkimCreatorLfCascades {
35023502 {
35033503 // pt
35043504 if (config.ptMinXicplusLfCasc > 0 .f ) {
3505- auto pt = RecoDecay::pt (pVecTrack0, pVecTrack1, pVecTrack2) + config.ptTolerance ; // add tolerance because of no reco decay vertex
3505+ const auto pt = RecoDecay::pt (pVecTrack0, pVecTrack1, pVecTrack2) + config.ptTolerance ; // add tolerance because of no reco decay vertex
35063506 if (pt < config.ptMinXicplusLfCasc ) {
35073507 return false ;
35083508 }
35093509 }
35103510
35113511 // invariant mass
3512- double invMassMin = config.massXiPiPiMin ;
3513- double invMassMax = config.massXiPiPiMax ;
3514- if (invMassMin >= 0 . && invMassMax > 0 .) {
3515- auto arrMom = std::array{pVecTrack0, pVecTrack1, pVecTrack2};
3516- auto invMass2 = RecoDecay::m2 (arrMom, arrMass3Prong[hf_cand_casc_lf::DecayType3Prong::XicplusToXiPiPi]);
3512+ if (config. massXiPiPiMin >= 0 . && config.massXiPiPiMax > 0 .) {
3513+ const double invMassMin = config.massXiPiPiMin ;
3514+ const double invMassMax = config. massXiPiPiMax ;
3515+ const auto arrMom = std::array{pVecTrack0, pVecTrack1, pVecTrack2};
3516+ const auto invMass2 = RecoDecay::m2 (arrMom, arrMass3Prong[hf_cand_casc_lf::DecayType3Prong::XicplusToXiPiPi]);
35173517 if (invMass2 < invMassMin * invMassMin || invMass2 >= invMassMax * invMassMax) {
35183518 return false ;
35193519 }
@@ -3532,23 +3532,23 @@ struct HfTrackIndexSkimCreatorLfCascades {
35323532 {
35333533 // pt
35343534 if (config.ptMinXicplusLfCasc > 0 .f ) {
3535- auto pt = RecoDecay::pt (pVecCand);
3535+ const auto pt = RecoDecay::pt (pVecCand);
35363536 if (pt < config.ptMinXicplusLfCasc ) {
35373537 return false ;
35383538 }
35393539 }
35403540
35413541 // CPA
35423542 if (config.xicCosPA > -1 .f ) {
3543- auto cpa = RecoDecay::cpa (primVtx, secVtx, pVecCand);
3543+ const auto cpa = RecoDecay::cpa (primVtx, secVtx, pVecCand);
35443544 if (cpa < config.xicCosPA ) {
35453545 return false ;
35463546 }
35473547 }
35483548
35493549 // decay length
35503550 if (config.decayLengthXicMin > 0 .f ) {
3551- auto decayLength = RecoDecay::distance (primVtx, secVtx);
3551+ const auto decayLength = RecoDecay::distance (primVtx, secVtx);
35523552 if (decayLength < config.decayLengthXicMin ) {
35533553 return false ;
35543554 }
@@ -3577,28 +3577,28 @@ struct HfTrackIndexSkimCreatorLfCascades {
35773577 for (const auto & collision : collisions) {
35783578
35793579 // set the magnetic field from CCDB
3580- auto bc = collision.bc_as <o2::aod::BCsWithTimestamps>();
3580+ const auto bc = collision.bc_as <o2::aod::BCsWithTimestamps>();
35813581 initCCDB (bc, runNumber, ccdb, config.isRun2 ? config.ccdbPathGrp : config.ccdbPathGrpMag , lut, config.isRun2 );
3582- auto magneticField = o2::base::Propagator::Instance ()->getNominalBz (); // z component
3582+ const auto magneticField = o2::base::Propagator::Instance ()->getNominalBz (); // z component
35833583
35843584 df2.setBz (magneticField);
35853585 df2.setRefitWithMatCorr (config.refitWithMatCorr );
35863586
35873587 // cascade loop
3588- auto thisCollId = collision.globalIndex ();
3589- auto groupedCascades = cascades.sliceBy (cascadesPerCollision, thisCollId);
3588+ const auto thisCollId = collision.globalIndex ();
3589+ const auto groupedCascades = cascades.sliceBy (cascadesPerCollision, thisCollId);
35903590
35913591 for (const auto & casc : groupedCascades) {
35923592
35933593 registry.fill (HIST (" hCandidateCounter" ), 0.5 ); // all cascade candidates
35943594
35953595 // ----------------accessing particles in the decay chain-------------
35963596 // cascade daughter - charged particle
3597- auto trackCascDauCharged = casc.bachelor_as <aod::TracksWCovDca>(); // meson <- xi track
3597+ const auto trackCascDauCharged = casc.bachelor_as <aod::TracksWCovDca>(); // meson <- xi track
35983598 // cascade daughter - V0
3599- auto trackV0PosDau = casc.posTrack_as <aod::TracksWCovDca>(); // p <- V0 track (positive track) 0
3599+ const auto trackV0PosDau = casc.posTrack_as <aod::TracksWCovDca>(); // p <- V0 track (positive track) 0
36003600 // V0 negative daughter
3601- auto trackV0NegDau = casc.negTrack_as <aod::TracksWCovDca>(); // pion <- V0 track (negative track) 1
3601+ const auto trackV0NegDau = casc.negTrack_as <aod::TracksWCovDca>(); // pion <- V0 track (negative track) 1
36023602
36033603 // check that particles come from the same collision
36043604 if (config.rejDiffCollTrack ) {
@@ -3618,8 +3618,8 @@ struct HfTrackIndexSkimCreatorLfCascades {
36183618 continue ;
36193619 }
36203620
3621- std::array<float , 3 > vertexCasc = {casc.x (), casc.y (), casc.z ()};
3622- std::array<float , 3 > pVecCasc = {casc.px (), casc.py (), casc.pz ()};
3621+ const std::array<float , 3 > vertexCasc = {casc.x (), casc.y (), casc.z ()};
3622+ const std::array<float , 3 > pVecCasc = {casc.px (), casc.py (), casc.pz ()};
36233623 std::array<float , 21 > covCasc = {0 .};
36243624 constexpr int MomInd[6 ] = {9 , 13 , 14 , 18 , 19 , 20 }; // cov matrix elements for momentum component
36253625 for (int i = 0 ; i < 6 ; i++) {
@@ -3643,14 +3643,14 @@ struct HfTrackIndexSkimCreatorLfCascades {
36433643 trackParCovCascOmega.setPID (o2::track::PID::OmegaMinus);
36443644
36453645 // --------------combining cascade and pion tracks--------------
3646- auto groupedBachTrackIndices = trackIndices.sliceBy (trackIndicesPerCollision, thisCollId);
3646+ const auto groupedBachTrackIndices = trackIndices.sliceBy (trackIndicesPerCollision, thisCollId);
36473647 for (auto trackIdCharmBachelor1 = groupedBachTrackIndices.begin (); trackIdCharmBachelor1 != groupedBachTrackIndices.end (); ++trackIdCharmBachelor1) {
36483648
36493649 hfFlag = 0 ;
36503650 isGoogForXi2Prong = true ;
36513651 isGoogForOmega2Prong = true ;
36523652
3653- auto trackCharmBachelor1 = trackIdCharmBachelor1.track_as <aod::TracksWCovDca>();
3653+ const auto trackCharmBachelor1 = trackIdCharmBachelor1.track_as <aod::TracksWCovDca>();
36543654
36553655 if ((config.rejDiffCollTrack ) && (trackCascDauCharged.collisionId () != trackCharmBachelor1.collisionId ())) {
36563656 continue ;
@@ -3667,7 +3667,7 @@ struct HfTrackIndexSkimCreatorLfCascades {
36673667 }
36683668
36693669 // primary pion track to be processed with DCAFitter
3670- auto trackParCovCharmBachelor1 = getTrackParCov (trackCharmBachelor1);
3670+ const auto trackParCovCharmBachelor1 = getTrackParCov (trackCharmBachelor1);
36713671
36723672 // find charm baryon decay using xi PID hypothesis (xi pi channel)
36733673 int nVtxFrom2ProngFitterXiHyp = 0 ;
@@ -3800,7 +3800,7 @@ struct HfTrackIndexSkimCreatorLfCascades {
38003800 continue ;
38013801 }
38023802
3803- auto trackCharmBachelor2 = trackIdCharmBachelor2.track_as <aod::TracksWCovDca>();
3803+ const auto trackCharmBachelor2 = trackIdCharmBachelor2.track_as <aod::TracksWCovDca>();
38043804
38053805 if ((config.rejDiffCollTrack ) && (trackCascDauCharged.collisionId () != trackCharmBachelor2.collisionId ())) {
38063806 continue ;
@@ -3841,7 +3841,7 @@ struct HfTrackIndexSkimCreatorLfCascades {
38413841 if (df2.isPropagateTracksToVertexDone ()) {
38423842 std::array<float , 3 > pVec1{0 .};
38433843 std::array<float , 3 > pVec2{0 .};
3844- std::array<float , 3 > pVec3{pVecCasc}; // Use the Xi track for the 3-prong calculations.
3844+ const std::array<float , 3 > pVec3{pVecCasc}; // Use the Xi track for the 3-prong calculations.
38453845 // get bachelor momenta at the Xic vertex
38463846 df2.getTrack (0 ).getPxPyPzGlo (pVec1);
38473847 df2.getTrack (1 ).getPxPyPzGlo (pVec2);
0 commit comments