@@ -160,6 +160,8 @@ struct QuarkoniaToHyperons {
160160 Configurable<bool > requireNoCollInTimeRangeStd{" requireNoCollInTimeRangeStd" , true , " reject collisions corrupted by the cannibalism, with other collisions within +/- 10 microseconds" };
161161 Configurable<bool > requireNoCollInTimeRangeNarrow{" requireNoCollInTimeRangeNarrow" , false , " reject collisions corrupted by the cannibalism, with other collisions within +/- 10 microseconds" };
162162
163+ Configurable<float > maxZVtxPosition{" maxZVtxPosition" , 10 ., " max Z vtx position" };
164+
163165 Configurable<bool > buildK0sK0sPairs{" buildK0sK0sPairs" , false , " Build K0s K0s from charmonia decay" };
164166 Configurable<bool > buildLaLaBarPairs{" buildLaLaBarPairs" , false , " Build Lambda antiLambda from charmonia decay" };
165167 Configurable<bool > buildXiXiBarPairs{" buildXiXiBarPairs" , false , " Build Xi antiXi from charmonia decay" };
@@ -1238,7 +1240,7 @@ struct QuarkoniaToHyperons {
12381240 if (fillHists)
12391241 histos.fill (HIST (" hEventSelection" ), 4 /* Not at TF border */ );
12401242
1241- if (std::abs (collision.posZ ()) > 10 . f ) {
1243+ if (std::abs (collision.posZ ()) > maxZVtxPosition ) {
12421244 return false ;
12431245 }
12441246 if (fillHists)
@@ -1712,18 +1714,18 @@ struct QuarkoniaToHyperons {
17121714
17131715 if (isXi) {
17141716 if (casc.sign () < 0 ) {
1715- if (cascMC.pdgCode () != 3312 || cascMC.pdgCodePositive () != 2212 || cascMC.pdgCodeNegative () != - 211 || cascMC.pdgCodeBachelor () != -211 )
1717+ if (cascMC.pdgCode () != PDG_t:: kXiMinus || cascMC.pdgCodePositive () != PDG_t:: kProton || cascMC.pdgCodeNegative () != PDG_t:: kPiMinus || cascMC.pdgCodeBachelor () != -PDG_t:: kPiMinus )
17161718 return false ;
17171719 } else {
1718- if (cascMC.pdgCode () != - 3312 || cascMC.pdgCodePositive () != 211 || cascMC.pdgCodeNegative () != - 2212 || cascMC.pdgCodeBachelor () != 211 )
1720+ if (cascMC.pdgCode () != PDG_t:: kXiPlus || cascMC.pdgCodePositive () != PDG_t:: kPiPlus || cascMC.pdgCodeNegative () != PDG_t:: kProtonBar || cascMC.pdgCodeBachelor () != PDG_t:: kPiPlus )
17191721 return false ;
17201722 }
17211723 } else {
17221724 if (casc.sign () < 0 ) {
1723- if (cascMC.pdgCode () != 3334 || cascMC.pdgCodePositive () != 2212 || cascMC.pdgCodeNegative () != - 211 || cascMC.pdgCodeBachelor () != - 321 )
1725+ if (cascMC.pdgCode () != PDG_t:: kOmegaMinus || cascMC.pdgCodePositive () != PDG_t:: kProton || cascMC.pdgCodeNegative () != PDG_t:: kPiMinus || cascMC.pdgCodeBachelor () != PDG_t:: kKMinus )
17241726 return false ;
17251727 } else {
1726- if (cascMC.pdgCode () != - 3334 || cascMC.pdgCodePositive () != 211 || cascMC.pdgCodeNegative () != - 2212 || cascMC.pdgCodeBachelor () != 321 )
1728+ if (cascMC.pdgCode () != PDG_t:: kOmegaPlus || cascMC.pdgCodePositive () != PDG_t:: kPiPlus || cascMC.pdgCodeNegative () != PDG_t:: kProtonBar || cascMC.pdgCodeBachelor () != PDG_t:: kKPlus )
17271729 return false ;
17281730 }
17291731 }
@@ -1740,17 +1742,17 @@ struct QuarkoniaToHyperons {
17401742 uint64_t bitMap = 0 ;
17411743 // check for specific particle species
17421744
1743- if (v0.pdgCode () == 310 && v0.pdgCodePositive () == 211 && v0.pdgCodeNegative () == - 211 ) {
1745+ if (v0.pdgCode () == PDG_t:: kK0Short && v0.pdgCodePositive () == PDG_t:: kPiPlus && v0.pdgCodeNegative () == PDG_t:: kPiMinus ) {
17441746 BITSET (bitMap, selConsiderK0Short);
17451747 if (v0.isPhysicalPrimary ())
17461748 BITSET (bitMap, selPhysPrimK0Short);
17471749 }
1748- if (v0.pdgCode () == 3122 && v0.pdgCodePositive () == 2212 && v0.pdgCodeNegative () == - 211 ) {
1750+ if (v0.pdgCode () == PDG_t:: kLambda0 && v0.pdgCodePositive () == PDG_t:: kProton && v0.pdgCodeNegative () == PDG_t:: kPiMinus ) {
17491751 BITSET (bitMap, selConsiderLambda);
17501752 if (v0.isPhysicalPrimary ())
17511753 BITSET (bitMap, selPhysPrimLambda);
17521754 }
1753- if (v0.pdgCode () == - 3122 && v0.pdgCodePositive () == 211 && v0.pdgCodeNegative () == - 2212 ) {
1755+ if (v0.pdgCode () == PDG_t:: kLambda0Bar && v0.pdgCodePositive () == PDG_t:: kPiPlus && v0.pdgCodeNegative () == PDG_t:: kProtonBar ) {
17541756 BITSET (bitMap, selConsiderAntiLambda);
17551757 if (v0.isPhysicalPrimary ())
17561758 BITSET (bitMap, selPhysPrimAntiLambda);
@@ -2730,9 +2732,9 @@ struct QuarkoniaToHyperons {
27302732
27312733 float ptmc = RecoDecay::sqrtSumOfSquares (v0MC.pxPosMC () + v0MC.pxNegMC (), v0MC.pyPosMC () + v0MC.pyNegMC ());
27322734 float ymc = 1e-3 ;
2733- if (v0MC.pdgCode () == 310 )
2735+ if (v0MC.pdgCode () == PDG_t:: kK0Short )
27342736 ymc = RecoDecay::y (std::array{v0MC.pxPosMC () + v0MC.pxNegMC (), v0MC.pyPosMC () + v0MC.pyNegMC (), v0MC.pzPosMC () + v0MC.pzNegMC ()}, o2::constants::physics::MassKaonNeutral);
2735- else if (std::fabs (v0MC.pdgCode ()) == 3122 )
2737+ else if (std::fabs (v0MC.pdgCode ()) == PDG_t:: kLambda0 )
27362738 ymc = RecoDecay::y (std::array{v0MC.pxPosMC () + v0MC.pxNegMC (), v0MC.pyPosMC () + v0MC.pyNegMC (), v0MC.pzPosMC () + v0MC.pzNegMC ()}, o2::constants::physics::MassLambda);
27372739
27382740 selMap[v0sGrouped[collision.globalIndex ()][i]] = computeReconstructionBitmap (v0, collision, ymc, ymc, ptmc);
@@ -2804,9 +2806,9 @@ struct QuarkoniaToHyperons {
28042806 auto cascadeMC = cascade.cascMCCore_as <soa::Join<aod::CascMCCores, aod::CascMCCollRefs>>();
28052807
28062808 float ymc = 1e-3 ;
2807- if (std::fabs (cascadeMC.pdgCode ()) == 3312 )
2809+ if (std::fabs (cascadeMC.pdgCode ()) == PDG_t:: kXiMinus )
28082810 ymc = RecoDecay::y (std::array{cascadeMC.pxMC (), cascadeMC.pyMC (), cascadeMC.pzMC ()}, o2::constants::physics::MassXiMinus);
2809- else if (std::fabs (cascadeMC.pdgCode ()) == 3334 )
2811+ else if (std::fabs (cascadeMC.pdgCode ()) == PDG_t:: kOmegaMinus )
28102812 ymc = RecoDecay::y (std::array{cascadeMC.pxMC (), cascadeMC.pyMC (), cascadeMC.pzMC ()}, o2::constants::physics::MassOmegaMinus);
28112813
28122814 if (buildXiXiBarPairs) {
0 commit comments