@@ -126,12 +126,12 @@ struct Derivedcascadeanalysis {
126126 Configurable<float > globalTracksCorrelpar0High{" globalTracksCorrelpar0High" , 226 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
127127 Configurable<float > globalTracksCorrelpar1High{" globalTracksCorrelpar1High" , -0.0181686 , " [0]*exp([1]*centrality)+[2], mean plus 3*sigma" };
128128 Configurable<float > globalTracksCorrelpar2High{" globalTracksCorrelpar2High" , -22 , " [0]*exp([1]*centrality)+[2], mean plus 3*sigma" };
129- Configurable<float > pvContribCorrelpar0Low{" pvcontribCorrelpar0Low " , 152 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
130- Configurable<float > pvContribCorrelpar1Low{" pvcontribCorrelpar1Low " , -0.0431016 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
131- Configurable<float > pvContribCorrelpar2Low{" pvcontribCorrelpar2Low " , -15.3776 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
132- Configurable<float > pvContribCorrelpar0High{" pvcontribCorrelpar0High " , 384.861 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
133- Configurable<float > pvContribCorrelpar1High{" pvcontribCorrelpar1High " , -0.0181686 , " [0]*exp([1]*centrality)+[2], mean plus 3*sigma" };
134- Configurable<float > pvContribCorrelpar2High{" pvcontribCorrelpar2High " , -39 , " [0]*exp([1]*centrality)+[2], mean plus 3*sigma" };
129+ Configurable<float > pvContribCorrelpar0Low{" pvContribCorrelpar0Low " , 152 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
130+ Configurable<float > pvContribCorrelpar1Low{" pvContribCorrelpar1Low " , -0.0431016 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
131+ Configurable<float > pvContribCorrelpar2Low{" pvContribCorrelpar2Low " , -15.3776 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
132+ Configurable<float > pvContribCorrelpar0High{" pvContribCorrelpar0High " , 384.861 , " [0]*exp([1]*centrality)+[2], mean minus 3*sigma" };
133+ Configurable<float > pvContribCorrelpar1High{" pvContribCorrelpar1High " , -0.0181686 , " [0]*exp([1]*centrality)+[2], mean plus 3*sigma" };
134+ Configurable<float > pvContribCorrelpar2High{" pvContribCorrelpar2High " , -39 , " [0]*exp([1]*centrality)+[2], mean plus 3*sigma" };
135135 } eventSelectionRun3Flags;
136136
137137 struct : ConfigurableGroup {
@@ -1145,11 +1145,11 @@ struct Derivedcascadeanalysis {
11451145 auto cascMC = casc.template cascMCCore_as <soa::Join<aod::CascMCCores, aod::CascMCCollRefs>>();
11461146 ptmc = RecoDecay::sqrtSumOfSquares (cascMC.pxMC (), cascMC.pyMC ());
11471147
1148- if (cascMC.isPhysicalPrimary () && ((isXi && std::abs (cascMC.pdgCode ()) == 3312 ) || (!isXi && std::abs (cascMC.pdgCode ()) == 3334 )))
1148+ if (cascMC.isPhysicalPrimary () && ((isXi && std::abs (cascMC.pdgCode ()) == PDG_t:: kXiMinus ) || (!isXi && std::abs (cascMC.pdgCode ()) == PDG_t:: kOmegaMinus )))
11491149 isTrueMCCascade = true ;
1150- if (isTrueMCCascade && ((isPositive && cascMC.pdgCodePositive () == 211 && cascMC.pdgCodeNegative () == - 2212 ) || (isNegative && cascMC.pdgCodePositive () == 2212 && cascMC.pdgCodeNegative () == - 211 )))
1150+ if (isTrueMCCascade && ((isPositive && cascMC.pdgCodePositive () == PDG_t:: kPiPlus && cascMC.pdgCodeNegative () == PDG_t:: kProtonBar ) || (isNegative && cascMC.pdgCodePositive () == PDG_t:: kProton && cascMC.pdgCodeNegative () == PDG_t:: kPiMinus )))
11511151 isCorrectLambdaDecay = true ;
1152- if (isTrueMCCascade && isCorrectLambdaDecay && ((isXi && std::abs (cascMC.pdgCodeBachelor ()) == 211 ) || (!isXi && std::abs (cascMC.pdgCodeBachelor ()) == 321 )))
1152+ if (isTrueMCCascade && isCorrectLambdaDecay && ((isXi && std::abs (cascMC.pdgCodeBachelor ()) == PDG_t:: kPiPlus ) || (!isXi && std::abs (cascMC.pdgCodeBachelor ()) == PDG_t:: kKPlus )))
11531153 isTrueMCCascadeDecay = true ;
11541154
11551155 if (qaFlags.doBefSelCheck && isTrueMCCascade) {
@@ -1458,9 +1458,9 @@ struct Derivedcascadeanalysis {
14581458
14591459 float ptmc = RecoDecay::sqrtSumOfSquares (cascMC.pxMC (), cascMC.pyMC ());
14601460 float ymc = 1e3 ;
1461- if (std::abs (cascMC.pdgCode ()) == 3312 )
1461+ if (std::abs (cascMC.pdgCode ()) == PDG_t:: kXiMinus )
14621462 ymc = RecoDecay::y (std::array{cascMC.pxMC (), cascMC.pyMC (), cascMC.pzMC ()}, o2::constants::physics::MassXiMinus);
1463- else if (std::abs (cascMC.pdgCode ()) == 3334 )
1463+ else if (std::abs (cascMC.pdgCode ()) == PDG_t:: kOmegaMinus )
14641464 ymc = RecoDecay::y (std::array{cascMC.pxMC (), cascMC.pyMC (), cascMC.pzMC ()}, o2::constants::physics::MassOmegaMinus);
14651465
14661466 if (ispO && (ymc > candidateSelectionValues.maxRapCut || ymc < candidateSelectionValues.minRapCut ))
@@ -1502,7 +1502,7 @@ struct Derivedcascadeanalysis {
15021502 nChEta1 = collision.multNTracksPVeta1 ();
15031503 }
15041504
1505- if (cascMC.pdgCode () == 3312 && isXi) {
1505+ if (cascMC.pdgCode () == PDG_t:: kXiMinus && isXi) {
15061506 histos.fill (HIST (" h2dGenXiMinus" ), centrality, ptmc);
15071507 histos.fill (HIST (" h2dGenXiMinusVsNch" ), nChEta1, ptmc);
15081508 histos.fill (HIST (" h2dGenXiMinusEta" ), RecoDecay::eta (std::array{cascMC.pxMC (), cascMC.pyMC (), cascMC.pzMC ()}));
@@ -1515,7 +1515,7 @@ struct Derivedcascadeanalysis {
15151515 histos.fill (HIST (" h2dGenXiMinusVsCentIR" ), ptmc, centrality, intRate);
15161516 histos.fill (HIST (" h2dGenXiMinusVsNchVsOccupancy" ), ptmc, nChEta1, occupancy);
15171517 }
1518- if (cascMC.pdgCode () == - 3312 && isXi) {
1518+ if (cascMC.pdgCode () == PDG_t:: kXiPlusBar && isXi) {
15191519 histos.fill (HIST (" h2dGenXiPlus" ), centrality, ptmc);
15201520 histos.fill (HIST (" h2dGenXiPlusVsNch" ), nChEta1, ptmc);
15211521 histos.fill (HIST (" h2dGenXiPlusVsMultMCVsCentrality" ), mcCollision.multMCNParticlesEta05 (), centrality, ptmc);
@@ -1524,7 +1524,7 @@ struct Derivedcascadeanalysis {
15241524 histos.fill (HIST (" h2dGenXiPlusVsNchVsOccupancy" ), ptmc, nChEta1, occupancy);
15251525 histos.fill (HIST (" h2dGenXiPlusVsCentIR" ), ptmc, centrality, intRate);
15261526 }
1527- if (cascMC.pdgCode () == 3334 && !isXi) {
1527+ if (cascMC.pdgCode () == PDG_t:: kOmegaMinus && !isXi) {
15281528 histos.fill (HIST (" h2dGenOmegaMinus" ), centrality, ptmc);
15291529 histos.fill (HIST (" h2dGenOmegaMinusVsNch" ), nChEta1, ptmc);
15301530 histos.fill (HIST (" h2dGenOmegaMinusEta" ), RecoDecay::eta (std::array{cascMC.pxMC (), cascMC.pyMC (), cascMC.pzMC ()}));
@@ -1537,7 +1537,7 @@ struct Derivedcascadeanalysis {
15371537 histos.fill (HIST (" h2dGenOmegaMinusVsNchVsOccupancy" ), ptmc, nChEta1, occupancy);
15381538 histos.fill (HIST (" h2dGenOmegaMinusVsCentIR" ), ptmc, centrality, intRate);
15391539 }
1540- if (cascMC.pdgCode () == - 3334 && !isXi) {
1540+ if (cascMC.pdgCode () == PDG_t:: kOmegaPlusBar && !isXi) {
15411541 histos.fill (HIST (" h2dGenOmegaPlus" ), centrality, ptmc);
15421542 histos.fill (HIST (" h2dGenOmegaPlusVsNch" ), nChEta1, ptmc);
15431543 histos.fill (HIST (" h2dGenOmegaPlusVsMultMCVsCentrality" ), mcCollision.multMCNParticlesEta05 (), centrality, ptmc);
0 commit comments