@@ -536,7 +536,7 @@ struct lambdaspincorrderived {
536536 }
537537 proton2 = ROOT::Math::PtEtaPhiMVector (v02.protonPt (), v02.protonEta (), v02.protonPhi (), o2::constants::physics::MassProton);
538538 lambda2 = ROOT::Math::PtEtaPhiMVector (v02.lambdaPt (), v02.lambdaEta (), v02.lambdaPhi (), v02.lambdaMass ());
539- histos.fill (HIST (" deltaPhiSame" ), std::abs ( RecoDecay::constrainAngle (v0.lambdaPhi (), 0 . 0F , harmonic) - RecoDecay::constrainAngle ( v02.lambdaPhi (), 0 .0F , harmonic) ));
539+ histos.fill (HIST (" deltaPhiSame" ), RecoDecay::constrainAngle (v0.lambdaPhi () - v02.lambdaPhi (), 0 .0F , harmonicDphi ));
540540 if (v0.v0Status () == 0 && v02.v0Status () == 0 ) {
541541 fillHistograms (0 , 0 , lambda, lambda2, proton, proton2, 0 , 1.0 );
542542 }
@@ -608,7 +608,7 @@ struct lambdaspincorrderived {
608608 lambda = ROOT::Math::PtEtaPhiMVector (t3.lambdaPt (), t3.lambdaEta (), t3.lambdaPhi (), t3.lambdaMass ());
609609 proton2 = ROOT::Math::PtEtaPhiMVector (t2.protonPt (), t2.protonEta (), t2.protonPhi (), o2::constants::physics::MassProton);
610610 lambda2 = ROOT::Math::PtEtaPhiMVector (t2.lambdaPt (), t2.lambdaEta (), t2.lambdaPhi (), t2.lambdaMass ());
611- histos.fill (HIST (" deltaPhiMix" ), std::abs ( RecoDecay::constrainAngle (t3.lambdaPhi (), 0 . 0F , harmonic) - RecoDecay::constrainAngle ( t2.lambdaPhi (), 0 .0F , harmonic) ));
611+ histos.fill (HIST (" deltaPhiMix" ), RecoDecay::constrainAngle (t3.lambdaPhi () - t2.lambdaPhi (), 0 .0F , harmonicDphi ));
612612 if (t3.v0Status () == 0 && t2.v0Status () == 0 ) {
613613 fillHistograms (0 , 0 , lambda, lambda2, proton, proton2, 1 , 1.0 );
614614 }
@@ -1030,9 +1030,7 @@ struct lambdaspincorrderived {
10301030 auto proton2 = ROOT::Math::PtEtaPhiMVector (t2.protonPt (), t2.protonEta (), t2.protonPhi (), o2::constants::physics::MassProton);
10311031 auto lambda2 = ROOT::Math::PtEtaPhiMVector (t2.lambdaPt (), t2.lambdaEta (), t2.lambdaPhi (), t2.lambdaMass ());
10321032
1033- const float dPhi = std::fabs (
1034- RecoDecay::constrainAngle (lambda.Phi (), 0 .0F , harmonic) -
1035- RecoDecay::constrainAngle (lambda2.Phi (), 0 .0F , harmonic));
1033+ const float dPhi = std::fabs (RecoDecay::constrainAngle (lambda.Phi () - lambda2.Phi (), 0 .0F , harmonicDphi));
10361034 histos.fill (HIST (" deltaPhiMix" ), dPhi, wBase);
10371035 fillHistograms (tX.v0Status (), t2.v0Status (), lambda, lambda2, proton, proton2, 1 , wBase);
10381036 }
0 commit comments