@@ -591,7 +591,8 @@ struct cascadeFlow {
591591 double ComputeEPResolutionwShifts (const TCollision& coll, double psiT0C, double psiV0A, double psiT0A, double psiTPCA, double psiTPCC, TProfile3D* shiftprofileA, TProfile3D* shiftprofileB, TProfile3D* shiftprofileC, TProfile3D* shiftprofileD, TProfile3D* shiftprofileE)
592592 {
593593 float collcentrality = coll.centFT0C ();
594- if (isCollisionCentrality == 1 ) collcentrality = coll.centFT0M ();
594+ if (isCollisionCentrality == 1 )
595+ collcentrality = coll.centFT0M ();
595596 int nmode = 2 ;
596597 auto deltapsiFT0C = 0.0 ;
597598 auto deltapsiFV0A = 0.0 ;
@@ -806,15 +807,15 @@ struct cascadeFlow {
806807 invMassLambda = v0.mLambda ();
807808 double ctauLambda = v0.distovertotmom (coll.posX (), coll.posY (), coll.posZ ()) * o2::constants::physics::MassLambda0;
808809 analysisLambdaSample (coll.centFT0C (),
809- coll.centFT0M (),
810+ coll.centFT0M (),
810811 hasEventPlane,
811812 hasSpectatorPlane,
812813 chargeIndex,
813814 v0.pt (),
814815 v0.phi (),
815816 v0.eta (),
816817 invMassLambda,
817- v0.mK0Short (),
818+ v0.mK0Short (),
818819 ctauLambda,
819820 v0.v0radius (),
820821 v0.dcapostopv (),
@@ -2753,15 +2754,17 @@ struct cascadeFlow {
27532754 theta1 = o2::constants::math::PI + theta; // pi/2 < theta1 < pi --> pi/4 < theta1/2 < pi/2 --> 1 < tan (theta1/2) --> negative eta
27542755
27552756 float cascMCeta = -std::log (std::tan (theta1 / 2 ));
2756-
2757+
27572758 float pxLambda = cascmc.pxPosMC () + cascmc.pxNegMC ();
27582759 float pyLambda = cascmc.pyPosMC () + cascmc.pyNegMC ();
27592760 float pzLambda = cascmc.pzPosMC () + cascmc.pzNegMC ();
27602761 float ptmcLambda = RecoDecay::sqrtSumOfSquares (pxLambda, pyLambda);
27612762 float thetaLambda = std::atan (ptmcLambda / pzLambda);
27622763 float theta1Lambda = 0 ;
2763- if (thetaLambda > 0 ) theta1Lambda = thetaLambda;
2764- else theta1Lambda = o2::constants::math::PI + thetaLambda;
2764+ if (thetaLambda > 0 )
2765+ theta1Lambda = thetaLambda;
2766+ else
2767+ theta1Lambda = o2::constants::math::PI + thetaLambda;
27652768 float lambdaMCeta = -std::log (std::tan (theta1Lambda / 2 ));
27662769
27672770 float cascMCy = 0 ;
@@ -2774,10 +2777,10 @@ struct cascadeFlow {
27742777 if (std::abs (cascMCy) < yCascMCGen)
27752778 histosMCGen.fill (HIST (" h2DGenXiY05" ), centrality, ptmc);
27762779 histosMCGen.fill (HIST (" hGenXiY" ), cascMCy);
2777- if (std::abs (cascMCy) < yXiToLambdaMCGen && std::abs (lambdaMCeta) < etaLambdaFromXiMCGen){
2778- histosMCGen.fill (HIST (" h2DGenXiVsPtLambda" ), centrality, ptmcLambda); // to compute secondary lambda efficiency
2779- histosMCGen.fill (HIST (" h2DGenLambdaFromXiPtMatrix" ), ptmc, ptmcLambda);
2780- }
2780+ if (std::abs (cascMCy) < yXiToLambdaMCGen && std::abs (lambdaMCeta) < etaLambdaFromXiMCGen) {
2781+ histosMCGen.fill (HIST (" h2DGenXiVsPtLambda" ), centrality, ptmcLambda); // to compute secondary lambda efficiency
2782+ histosMCGen.fill (HIST (" h2DGenLambdaFromXiPtMatrix" ), ptmc, ptmcLambda);
2783+ }
27812784 } else if (std::abs (cascmc.pdgCode ()) == PDG_t::kOmegaMinus ) {
27822785 cascMCy = RecoDecay::y (std::array{cascmc.pxMC (), cascmc.pyMC (), cascmc.pzMC ()}, constants::physics::MassOmegaMinus);
27832786 if (std::abs (cascMCeta) < etaCascMCGen) {
@@ -2899,9 +2902,9 @@ struct cascadeFlow {
28992902 } else {
29002903 histos.fill (HIST (" hCentvsPtvsPrimaryFracLambda" ), collisionCentrality, v0.pt (), 1 );
29012904 histos.fill (HIST (" hCentvsPrimaryFracLambda" ), collisionCentrality, 1 );
2902- histosMCReco.fill (HIST (" h2DRecoTrueLambdaSec" ), collisionCentrality, ptmc);
2903- if (v0MC.pdgCodeMother () == PDG_t::kXiMinus )
2904- histosMCReco.fill (HIST (" h2DRecoTrueLambdaFromXi" ), collisionCentrality, ptmc);
2905+ histosMCReco.fill (HIST (" h2DRecoTrueLambdaSec" ), collisionCentrality, ptmc);
2906+ if (v0MC.pdgCodeMother () == PDG_t::kXiMinus )
2907+ histosMCReco.fill (HIST (" h2DRecoTrueLambdaFromXi" ), collisionCentrality, ptmc);
29052908 }
29062909 } else if (isTrueALambda) {
29072910 if (isPrimary) {
@@ -2911,9 +2914,9 @@ struct cascadeFlow {
29112914 } else {
29122915 histos.fill (HIST (" hCentvsPtvsPrimaryFracLambda" ), collisionCentrality, v0.pt (), 3 );
29132916 histos.fill (HIST (" hCentvsPrimaryFracLambda" ), collisionCentrality, 3 );
2914- histosMCReco.fill (HIST (" h2DRecoTrueAntiLambdaSec" ), collisionCentrality, ptmc);
2915- if (v0MC.pdgCodeMother () == -PDG_t::kXiMinus )
2916- histosMCReco.fill (HIST (" h2DRecoTrueAntiLambdaFromXi" ), collisionCentrality, ptmc);
2917+ histosMCReco.fill (HIST (" h2DRecoTrueAntiLambdaSec" ), collisionCentrality, ptmc);
2918+ if (v0MC.pdgCodeMother () == -PDG_t::kXiMinus )
2919+ histosMCReco.fill (HIST (" h2DRecoTrueAntiLambdaFromXi" ), collisionCentrality, ptmc);
29172920 }
29182921 }
29192922 }
0 commit comments