diff --git a/PWGEM/Dilepton/Utils/SemiCharmTag.h b/PWGEM/Dilepton/Utils/SemiCharmTag.h index 82c965b0ccb..8f3610e1e23 100644 --- a/PWGEM/Dilepton/Utils/SemiCharmTag.h +++ b/PWGEM/Dilepton/Utils/SemiCharmTag.h @@ -175,8 +175,8 @@ LHPair makePairLeptonTrack(TFitter& fitter, TCollision const& collision, TLepton pair.mass = v12.M(); pair.pt = v12.Pt(); pair.p = v12.P(); - pair.deta = v1.Eta() - v2.Eta(); // lepton - hadron - pair.dphi = RecoDecay::constrainAngle(RecoDecay::constrainAngle(v1.Phi(), 0, 1U) - RecoDecay::constrainAngle(v2.Phi(), 0, 1U), -M_PI/2, 1U); // lepton - hadron + pair.deta = v1.Eta() - v2.Eta(); // lepton - hadron + pair.dphi = RecoDecay::constrainAngle(RecoDecay::constrainAngle(v1.Phi(), 0, 1U) - RecoDecay::constrainAngle(v2.Phi(), 0, 1U), -M_PI / 2, 1U); // lepton - hadron pair.isOK = true; return pair; } @@ -302,8 +302,8 @@ LHPair makePairLeptonV0(TFitter& fitter, TCollision const& collision, TLepton co pair.mass = v12.M(); pair.pt = v12.Pt(); pair.p = v12.P(); - pair.deta = v1.Eta() - v2.Eta(); // lepton - hadron - pair.dphi = RecoDecay::constrainAngle(RecoDecay::constrainAngle(v1.Phi(), 0, 1U) - RecoDecay::constrainAngle(v2.Phi(), 0, 1U), -M_PI/2, 1U); // lepton - hadron + pair.deta = v1.Eta() - v2.Eta(); // lepton - hadron + pair.dphi = RecoDecay::constrainAngle(RecoDecay::constrainAngle(v1.Phi(), 0, 1U) - RecoDecay::constrainAngle(v2.Phi(), 0, 1U), -M_PI / 2, 1U); // lepton - hadron pair.isOK = true; return pair; @@ -434,8 +434,8 @@ LHPair makePairLeptonCascade(TFitter& fitter, TCollision const& collision, TLept pair.mass = v12.M(); pair.pt = v12.Pt(); pair.p = v12.P(); - pair.deta = v1.Eta() - v2.Eta(); // lepton - hadron - pair.dphi = RecoDecay::constrainAngle(RecoDecay::constrainAngle(v1.Phi(), 0, 1U) - RecoDecay::constrainAngle(v2.Phi(), 0, 1U), -M_PI/2, 1U); // lepton - hadron + pair.deta = v1.Eta() - v2.Eta(); // lepton - hadron + pair.dphi = RecoDecay::constrainAngle(RecoDecay::constrainAngle(v1.Phi(), 0, 1U) - RecoDecay::constrainAngle(v2.Phi(), 0, 1U), -M_PI / 2, 1U); // lepton - hadron pair.isOK = true; return pair;