Skip to content

Commit f6e527d

Browse files
axti98axti98
committed
Fix signs
1 parent 990e3ba commit f6e527d

2 files changed

Lines changed: 36 additions & 35 deletions

File tree

‎PWGDQ/Core/DQConstants.h‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -161,6 +161,7 @@ inline constexpr int NumTrkCovParsNoX = 5;
161161
inline constexpr int NumTrkComponents = 3;
162162

163163
// Vertexing defaults
164+
inline constexpr float InvalidVertexingError = 1.e8f;
164165
inline constexpr float VertexingLxyzTolerance = 1.0e-8f;
165166
inline constexpr float InitialShift = 0.0f;
166167

‎PWGDQ/Core/VarManager.h‎

Lines changed: 35 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -1462,7 +1462,7 @@ class VarManager : public TObject
14621462

14631463
static float getDeltaPsiInRange(float psi1, float psi2, unsigned int harmonic)
14641464
{
1465-
return RecoDecay::constrainAngle(psi1 - psi2, o2::constants::math::PI / harmonic, harmonic);
1465+
return RecoDecay::constrainAngle(psi1 - psi2, -o2::constants::math::PI / harmonic, harmonic);
14661466
}
14671467
template <typename T, typename T1>
14681468
static o2::dataformats::VertexBase RecalculatePrimaryVertex(T const& track0, T const& track1, const T1& collision);
@@ -2079,7 +2079,7 @@ void VarManager::FillEvent(T const& event, float* values)
20792079
}
20802080

20812081
if (fgUsedVars[kRandomPsi2]) {
2082-
values[kRandomPsi2] = gRandom->Uniform(o2::constants::math::PIHalf, o2::constants::math::PIHalf);
2082+
values[kRandomPsi2] = gRandom->Uniform(-o2::constants::math::PIHalf, o2::constants::math::PIHalf);
20832083
}
20842084

20852085
if constexpr ((fillMap & Collision) > 0) {
@@ -3752,7 +3752,7 @@ void VarManager::FillEnergyCorrelatorsMC(T const& track, T1 const& t1, float* va
37523752
MassHadron = std::sqrt(t1.e() * t1.e() - t1.p() * t1.p());
37533753
}
37543754
ROOT::Math::PtEtaPhiMVector v1(track.pt(), track.eta(), track.phi(), o2::constants::physics::MassJPsi);
3755-
float deltaphi = RecoDecay::constrainAngle(track.phi() - t1.phi(), o2::constants::math::PIHalf);
3755+
float deltaphi = RecoDecay::constrainAngle(track.phi() - t1.phi(), -o2::constants::math::PIHalf);
37563756
float deltaeta = t1.eta() - track.eta();
37573757
ROOT::Math::PtEtaPhiMVector v2(t1.pt(), t1.eta(), t1.phi(), MassHadron);
37583758
float E_boost = LorentzTransformJpsihadroncosChi("weight_boost", v1, v2);
@@ -3767,7 +3767,7 @@ void VarManager::FillEnergyCorrelatorsMC(T const& track, T1 const& t1, float* va
37673767
values[kMCdeltaeta] = deltaeta;
37683768
values[kMCHadronPt] = t1.pt();
37693769
values[kMCHadronEta] = t1.eta();
3770-
values[kMCHadronPhi] = RecoDecay::constrainAngle(t1.phi(), o2::constants::math::PIHalf);
3770+
values[kMCHadronPhi] = RecoDecay::constrainAngle(t1.phi(), -o2::constants::math::PIHalf);
37713771
values[kMCHadronPdgCode] = t1.pdgCode();
37723772
values[kMCWeight] = E_boost / o2::constants::physics::MassJPsi * Accweight;
37733773

@@ -3783,7 +3783,7 @@ void VarManager::FillEnergyCorrelatorsMC(T const& track, T1 const& t1, float* va
37833783
float randomPhi_toward = o2::constants::math::PIHalf;
37843784
float randomPhi_away = o2::constants::math::PIHalf;
37853785

3786-
float deltaphitrans = RecoDecay::constrainAngle(track.phi() - t1.phi(), o2::constants::math::PI);
3786+
float deltaphitrans = RecoDecay::constrainAngle(track.phi() - t1.phi(), -o2::constants::math::PI);
37873787
if ((deltaphitrans > -Transhigh * o2::constants::math::PI && deltaphitrans < -Translow * o2::constants::math::PI) || (deltaphitrans > Translow * o2::constants::math::PI && deltaphitrans < Transhigh * o2::constants::math::PI)) {
37883788
randomPhi_trans = gRandom->Uniform(-o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
37893789
randomPhi_toward = gRandom->Uniform(-o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
@@ -3898,7 +3898,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
38983898
values[kPhi2] = t2.phi();
38993899

39003900
if (fgUsedVars[kDeltaPhiPair2]) {
3901-
values[kDeltaPhiPair2] = RecoDecay::constrainAngle(v1.Phi() - v2.Phi(), o2::constants::math::PIHalf);
3901+
values[kDeltaPhiPair2] = RecoDecay::constrainAngle(v1.Phi() - v2.Phi(), -o2::constants::math::PIHalf);
39023902
}
39033903

39043904
if (fgUsedVars[kDeltaEtaPair2]) {
@@ -4482,7 +4482,7 @@ void VarManager::FillPairME(T1 const& t1, T2 const& t2, float* values)
44824482
values[kPhi2] = t2.phi();
44834483

44844484
if (fgUsedVars[kDeltaPhiPair2]) {
4485-
values[kDeltaPhiPair2] = RecoDecay::constrainAngle(v1.Phi() - v2.Phi(), o2::constants::math::PIHalf);
4485+
values[kDeltaPhiPair2] = RecoDecay::constrainAngle(v1.Phi() - v2.Phi(), -o2::constants::math::PIHalf);
44864486
}
44874487

44884488
if (fgUsedVars[kDeltaEtaPair2]) {
@@ -4700,7 +4700,7 @@ void VarManager::FillPairME(T1 const& t1, T2 const& t2, float* values)
47004700
float R2EP = t1.sign() > 0 ? values[kTwoR2EP1] : values[kTwoR2EP2];
47014701

47024702
// reaction plane
4703-
float phi = RecoDecay::constrainAngle(v_daughter.Phi(), o2::constants::math::PI);
4703+
float phi = RecoDecay::constrainAngle(v_daughter.Phi(), -o2::constants::math::PI);
47044704
values[kDeltaPhiME_TPC] = phi > Psi2A ? phi - Psi2A : Psi2A - phi;
47054705
values[kDeltaPhiME_TPC] = values[kDeltaPhiME_TPC] > o2::constants::math::PI ? o2::constants::math::TwoPI - values[kDeltaPhiME_TPC] : values[kDeltaPhiME_TPC];
47064706
values[kDeltaPhiME_FT0A] = phi > Psi2B ? phi - Psi2B : Psi2B - phi;
@@ -4845,7 +4845,7 @@ void VarManager::FillPairMEAcrossTFs(T const& t1, T const& t2, float* values)
48454845
float Psi2C = values[kPsi2C];
48464846

48474847
// reaction plane
4848-
float phi = RecoDecay::constrainAngle(v_daughter.Phi(), o2::constants::math::PI);
4848+
float phi = RecoDecay::constrainAngle(v_daughter.Phi(), -o2::constants::math::PI);
48494849
values[kDeltaPhi_TPC] = phi > Psi2A ? phi - Psi2A : Psi2A - phi;
48504850
values[kDeltaPhi_TPC] = values[kDeltaPhi_TPC] > o2::constants::math::PI ? o2::constants::math::TwoPI - values[kDeltaPhi_TPC] : values[kDeltaPhi_TPC];
48514851
values[kDeltaPhi_FT0A] = phi > Psi2B ? phi - Psi2B : Psi2B - phi;
@@ -5097,7 +5097,7 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
50975097
ROOT::Math::PtEtaPhiMVector v_daughter = boostv12(t1.pdgCode() > 0 ? v1 : v2);
50985098

50995099
// reaction plane
5100-
float phi = RecoDecay::constrainAngle(v_daughter.Phi(), o2::constants::math::PI);
5100+
float phi = RecoDecay::constrainAngle(v_daughter.Phi(), -o2::constants::math::PI);
51015101
values[kDeltaPhi_Random] = phi - values[kRandomPsi2];
51025102
values[kDeltaPhi_Random] = values[kDeltaPhi_Random] > o2::constants::math::PI ? o2::constants::math::TwoPI - values[kDeltaPhi_Random] : values[kDeltaPhi_Random];
51035103
values[kDeltaPhi_MC] = phi - values[kMCEventPlaneAngle];
@@ -5395,15 +5395,15 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
53955395
if (std::fabs(values[kVertexingLxy]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
53965396
values[kVertexingLxy] = o2::dqconstants::defaults::VertexingLxyzTolerance;
53975397
}
5398-
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
5398+
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
53995399
if (std::fabs(values[kVertexingLz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
54005400
values[kVertexingLz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
54015401
}
5402-
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
5402+
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
54035403
if (std::fabs(values[kVertexingLxyz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
54045404
values[kVertexingLxyz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
54055405
}
5406-
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
5406+
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
54075407
values[kVertexingTauxy] = KFGeoTwoProng.GetPseudoProperDecayTime(KFPV, KFGeoTwoProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
54085408
values[kVertexingTauz] = -1 * dzPair2PV * KFGeoTwoProng.GetMass() / (std::abs(KFGeoTwoProng.GetPz()) * o2::constants::physics::LightSpeedCm2NS);
54095409
values[kVertexingPz] = std::abs(KFGeoTwoProng.GetPz());
@@ -5773,15 +5773,15 @@ void VarManager::FillTripletVertexing(C const& collision, T const& t1, T const&
57735773
if (std::fabs(values[kVertexingLxy]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
57745774
values[kVertexingLxy] = o2::dqconstants::defaults::VertexingLxyzTolerance;
57755775
}
5776-
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
5776+
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
57775777
if (std::fabs(values[kVertexingLz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
57785778
values[kVertexingLz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
57795779
}
5780-
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
5780+
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
57815781
if (std::fabs(values[kVertexingLxyz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
57825782
values[kVertexingLxyz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
57835783
}
5784-
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
5784+
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
57855785

57865786
values[kVertexingTauxy] = KFGeoThreeProng.GetPseudoProperDecayTime(KFPV, KFGeoThreeProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
57875787
values[kVertexingTauz] = -1 * dzTriplet2PV * KFGeoThreeProng.GetMass() / (std::abs(KFGeoThreeProng.GetPz()) * o2::constants::physics::LightSpeedCm2NS);
@@ -6033,15 +6033,15 @@ void VarManager::FillDileptonTrackVertexing(C const& collision, T1 const& lepton
60336033
if (std::fabs(values[kVertexingLxy]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
60346034
values[kVertexingLxy] = o2::dqconstants::defaults::VertexingLxyzTolerance;
60356035
}
6036-
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
6036+
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
60376037
if (std::fabs(values[kVertexingLz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
60386038
values[kVertexingLz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
60396039
}
6040-
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
6040+
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
60416041
if (std::fabs(values[kVertexingLxyz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
60426042
values[kVertexingLxyz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
60436043
}
6044-
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
6044+
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
60456045

60466046
if (fgUsedVars[kVertexingTauxy]) {
60476047
values[kVertexingTauxy] = KFGeoThreeProng.GetPseudoProperDecayTime(KFPV, KFGeoThreeProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
@@ -6675,7 +6675,7 @@ void VarManager::FillDileptonHadron(T1 const& dilepton, T2 const& hadron, float*
66756675
}
66766676

66776677
if (fgUsedVars[kDeltaPhi]) {
6678-
values[kDeltaPhi] = RecoDecay::constrainAngle(dilepton.phi() - hadron.phi(), o2::constants::math::PIHalf);
6678+
values[kDeltaPhi] = RecoDecay::constrainAngle(dilepton.phi() - hadron.phi(), -o2::constants::math::PIHalf);
66796679
}
66806680
if (fgUsedVars[kDeltaPhiSym]) {
66816681
double delta = std::abs(dilepton.phi() - hadron.phi());
@@ -6727,10 +6727,10 @@ void VarManager::FillEnergyCorrelatorTriple(C const& collision, T1 const& lepton
67276727
values[kWeight_before] = v2.Pt() / v1.M() * weight;
67286728
values[kPtDau] = v2.pt();
67296729
values[kEtaDau] = v2.eta();
6730-
values[kPhiDau] = RecoDecay::constrainAngle(v2.phi(), o2::constants::math::PIHalf);
6730+
values[kPhiDau] = RecoDecay::constrainAngle(v2.phi(), -o2::constants::math::PIHalf);
67316731

6732-
float deltaphi = RecoDecay::constrainAngle(v1.phi() - v2.phi(), o2::constants::math::PI);
6733-
values[kDeltaPhi] = RecoDecay::constrainAngle(v1.phi() - v2.phi(), o2::constants::math::PIHalf);
6732+
float deltaphi = RecoDecay::constrainAngle(v1.phi() - v2.phi(), -o2::constants::math::PI);
6733+
values[kDeltaPhi] = RecoDecay::constrainAngle(v1.phi() - v2.phi(), -o2::constants::math::PIHalf);
67346734
values[kDeltaEta] = v1.eta() - v2.eta();
67356735
values[kCosChi_randomPhi_trans] = o2::dqconstants::InvalidValue;
67366736
values[kCosChi_randomPhi_toward] = o2::dqconstants::InvalidValue;
@@ -6745,9 +6745,9 @@ void VarManager::FillEnergyCorrelatorTriple(C const& collision, T1 const& lepton
67456745
float randomPhi_away = o2::constants::math::PIHalf;
67466746

67476747
if ((deltaphi > -Transhigh * o2::constants::math::PI && deltaphi < -Translow * o2::constants::math::PI) || (deltaphi > Translow * o2::constants::math::PI && deltaphi < Transhigh * o2::constants::math::PI)) {
6748-
randomPhi_trans = gRandom->Uniform(o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
6749-
randomPhi_toward = gRandom->Uniform(o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
6750-
randomPhi_away = gRandom->Uniform(o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
6748+
randomPhi_trans = gRandom->Uniform(-o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
6749+
randomPhi_toward = gRandom->Uniform(-o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
6750+
randomPhi_away = gRandom->Uniform(-o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
67516751
values[kPtDau_randomPhi_trans] = v2.pt();
67526752
ROOT::Math::PtEtaPhiMVector v2_randomPhi_trans(v2.pt(), v2.eta(), randomPhi_trans, o2::constants::physics::MassPionCharged);
67536753
values[kCosChi_randomPhi_trans] = LorentzTransformJpsihadroncosChi("coschi", v1, v2_randomPhi_trans);
@@ -6760,9 +6760,9 @@ void VarManager::FillEnergyCorrelatorTriple(C const& collision, T1 const& lepton
67606760
values[kCosChi_randomPhi_away] = LorentzTransformJpsihadroncosChi("coschi", v1, v2_randomPhi_away);
67616761
values[kWeight_randomPhi_away] = LorentzTransformJpsihadroncosChi("weight_boost", v1, v2_randomPhi_away) / v1.M() * weight;
67626762

6763-
values[kdeltaphi_randomPhi_trans] = RecoDecay::constrainAngle(v1.phi() - randomPhi_trans, o2::constants::math::PIHalf);
6764-
values[kdeltaphi_randomPhi_toward] = RecoDecay::constrainAngle(v1.phi() - randomPhi_toward, o2::constants::math::PIHalf);
6765-
values[kdeltaphi_randomPhi_away] = RecoDecay::constrainAngle(v1.phi() - randomPhi_away, o2::constants::math::PIHalf);
6763+
values[kdeltaphi_randomPhi_trans] = RecoDecay::constrainAngle(v1.phi() - randomPhi_trans, -o2::constants::math::PIHalf);
6764+
values[kdeltaphi_randomPhi_toward] = RecoDecay::constrainAngle(v1.phi() - randomPhi_toward, -o2::constants::math::PIHalf);
6765+
values[kdeltaphi_randomPhi_away] = RecoDecay::constrainAngle(v1.phi() - randomPhi_away, -o2::constants::math::PIHalf);
67666766
}
67676767
}
67686768
}
@@ -6804,10 +6804,10 @@ void VarManager::FillEnergyCorrelatorsUnfoldingTriple(T1 const& lepton1, T2 cons
68046804

68056805
values[kMCCosChi_randomPhi_trans_rec] = o2::dqconstants::InvalidValue;
68066806
values[kMCCosChi_randomPhi_trans_gen] = o2::dqconstants::InvalidValue;
6807-
float deltaphi_rec = RecoDecay::constrainAngle(dilepton.phi() - hadron.phi(), o2::constants::math::PIHalf);
6807+
float deltaphi_rec = RecoDecay::constrainAngle(dilepton.phi() - hadron.phi(), -o2::constants::math::PIHalf);
68086808

68096809
if ((deltaphi_rec > -Transhigh * o2::constants::math::PI && deltaphi_rec < -Translow * o2::constants::math::PI) || (deltaphi_rec > Translow * o2::constants::math::PI && deltaphi_rec < Transhigh * o2::constants::math::PI)) {
6810-
float randomPhi_trans_rec = gRandom->Uniform(o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
6810+
float randomPhi_trans_rec = gRandom->Uniform(-o2::constants::math::PIHalf, 3. * o2::constants::math::PIHalf);
68116811
ROOT::Math::PtEtaPhiMVector v2_randomPhi_trans_rec(v2_rec.pt(), v2_rec.eta(), randomPhi_trans_rec, o2::constants::physics::MassPionCharged);
68126812
values[kMCCosChi_randomPhi_trans_rec] = LorentzTransformJpsihadroncosChi("coschi", v1_rec, v2_randomPhi_trans_rec);
68136813
values[kMCWeight_randomPhi_trans_rec] = LorentzTransformJpsihadroncosChi("weight_boost", v1_rec, v2_randomPhi_trans_rec) / v1_rec.M() * Effweight_rec;
@@ -7196,15 +7196,15 @@ void VarManager::FillDileptonTrackTrackVertexing(C const& collision, T1 const& l
71967196
if (std::fabs(values[kVertexingLxy]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
71977197
values[kVertexingLxy] = o2::dqconstants::defaults::VertexingLxyzTolerance;
71987198
}
7199-
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
7199+
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
72007200
if (std::fabs(values[kVertexingLz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
72017201
values[kVertexingLz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
72027202
}
7203-
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
7203+
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
72047204
if (std::fabs(values[kVertexingLxyz]) < o2::dqconstants::defaults::VertexingLxyzTolerance) {
72057205
values[kVertexingLxyz] = o2::dqconstants::defaults::VertexingLxyzTolerance;
72067206
}
7207-
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::VertexingLxyzTolerance : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
7207+
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? o2::dqconstants::defaults::InvalidVertexingError : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
72087208

72097209
values[kVertexingTauxy] = KFGeoFourProng.GetPseudoProperDecayTime(KFPV, KFGeoFourProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
72107210
values[kVertexingTauz] = -1 * dzQuadlet2PV * KFGeoFourProng.GetMass() / (std::abs(KFGeoFourProng.GetPz()) * o2::constants::physics::LightSpeedCm2NS);
@@ -7772,7 +7772,7 @@ void VarManager::FillPairAlice3(T1 const& t1, T2 const& t2, float* values)
77727772
values[kPhi2] = t2.phi();
77737773

77747774
if (fgUsedVars[kDeltaPhiPair2]) {
7775-
values[kDeltaPhiPair2] = RecoDecay::constrainAngle(v1.Phi() - v2.Phi(), o2::constants::math::PIHalf);
7775+
values[kDeltaPhiPair2] = RecoDecay::constrainAngle(v1.Phi() - v2.Phi(), -o2::constants::math::PIHalf);
77767776
}
77777777

77787778
if (fgUsedVars[kDeltaEtaPair2]) {

0 commit comments

Comments
 (0)