Skip to content

Commit 91afd7d

Browse files
mcoquet642mcoquet642
committed
clang format
1 parent e5843f5 commit 91afd7d

1 file changed

Lines changed: 50 additions & 52 deletions

File tree

‎PWGDQ/Core/VarManager.h‎

Lines changed: 50 additions & 52 deletions
Original file line numberDiff line numberDiff line change
@@ -84,8 +84,6 @@
8484
#include <utility>
8585
#include <vector>
8686

87-
#include <cmath>
88-
8987
using SMatrix55 = ROOT::Math::SMatrix<double, 5, 5, ROOT::Math::MatRepSym<double, 5>>;
9088
using SMatrix5 = ROOT::Math::SVector<double, 5>;
9189
using Vec3D = ROOT::Math::SVector<double, 3>;
@@ -179,10 +177,10 @@ class VarManager : public TObject
179177

180178
enum BarrelTrackFilteringBits {
181179
kIsConversionLeg = 0, // electron from conversions
182-
kIsK0sLeg = 1, // pion from K0s
183-
kIsLambdaLeg = 2, // proton or pion from Lambda
184-
kIsALambdaLeg = 3, // proton or pion from anti-Lambda
185-
kIsOmegaLeg = 4, // kaon from Omega baryon decay
180+
kIsK0sLeg = 1, // pion from K0s
181+
kIsLambdaLeg = 2, // proton or pion from Lambda
182+
kIsALambdaLeg = 3, // proton or pion from anti-Lambda
183+
kIsOmegaLeg = 4, // kaon from Omega baryon decay
186184
kDalitzBits = 5, // first bit for Dalitz tagged tracks
187185
kBarrelUserCutsBits = 13, // first bit for user track cuts
188186
kIsTPCPostcalibrated = 63 // tracks were postcalibrated for the TPC PID
@@ -3530,7 +3528,7 @@ void VarManager::FillTrackCollisionMC(T1 const& track, const std::array<double,
35303528
float m = o2::constants::physics::MassJPsi;
35313529
if (massHyp != 0.f) {
35323530
m = massHyp;
3533-
}
3531+
}
35343532

35353533
// displaced vertex is computed with decay product (track) and vertex collPos
35363534
values[kMCVertexingLxy] = (collPos[0] - track.vx()) * (collPos[0] - track.vx()) +
@@ -3725,10 +3723,10 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
37253723
double arg = scalar / Ptot12;
37263724
if (arg > 1.) {
37273725
arg = 1.;
3728-
}
3726+
}
37293727
if (arg < -1) {
37303728
arg = -1;
3731-
}
3729+
}
37323730
values[kOpeningAngle] = TMath::ACos(arg);
37333731
}
37343732
}
@@ -3755,7 +3753,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
37553753
ROOT::Math::XYZVectorF xaxis_HE{(yaxis_HE.Cross(zaxis_HE)).Unit()};
37563754
if (fgUsedVars[kCosThetaHE]) {
37573755
values[kCosThetaHE] = zaxis_HE.Dot(v_CM);
3758-
}
3756+
}
37593757
if (fgUsedVars[kPhiHE]) {
37603758
values[kPhiHE] = TMath::ATan2(yaxis_HE.Dot(v_CM), xaxis_HE.Dot(v_CM));
37613759
if (values[kPhiHE] < 0) {
@@ -3787,7 +3785,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
37873785
ROOT::Math::XYZVectorF xaxis_CS{(yaxis_CS.Cross(zaxis_CS)).Unit()};
37883786
if (fgUsedVars[kCosThetaCS]) {
37893787
values[kCosThetaCS] = zaxis_CS.Dot(v_CM);
3790-
}
3788+
}
37913789
if (fgUsedVars[kPhiCS]) {
37923790
values[kPhiCS] = TMath::ATan2(yaxis_CS.Dot(v_CM), xaxis_CS.Dot(v_CM));
37933791
if (values[kPhiCS] < 0) {
@@ -3851,7 +3849,7 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
38513849
ROOT::Math::XYZVectorF zaxis_RM(randomCostheta, std::sqrt(1 - randomCostheta * randomCostheta) * std::cos(randomPhi), std::sqrt(1 - randomCostheta * randomCostheta) * std::sin(randomPhi));
38523850
if (fgUsedVars[kCosThetaRM]) {
38533851
values[kCosThetaRM] = zaxis_RM.Dot(v_CM);
3854-
}
3852+
}
38553853
}
38563854
}
38573855

@@ -3981,7 +3979,7 @@ void VarManager::FillPairRotation(T1 const& t1, T2 const& t2, float* values)
39813979

39823980
if (rotationphi2 > 2. * TMath::Pi()) {
39833981
rotationphi2 -= 2. * TMath::Pi();
3984-
}
3982+
}
39853983

39863984
values[kCharge] = t1.sign() + t2.sign();
39873985
values[kCharge1] = t1.sign();
@@ -4270,7 +4268,7 @@ void VarManager::FillPairME(T1 const& t1, T2 const& t2, float* values)
42704268
ROOT::Math::XYZVectorF xaxis_HE{(yaxis_HE.Cross(zaxis_HE)).Unit()};
42714269
if (fgUsedVars[kCosThetaHE]) {
42724270
values[kCosThetaHE] = zaxis_HE.Dot(v_CM);
4273-
}
4271+
}
42744272
if (fgUsedVars[kPhiHE]) {
42754273
values[kPhiHE] = TMath::ATan2(yaxis_HE.Dot(v_CM), xaxis_HE.Dot(v_CM));
42764274
if (values[kPhiHE] < 0) {
@@ -4302,7 +4300,7 @@ void VarManager::FillPairME(T1 const& t1, T2 const& t2, float* values)
43024300
ROOT::Math::XYZVectorF xaxis_CS{(yaxis_CS.Cross(zaxis_CS)).Unit()};
43034301
if (fgUsedVars[kCosThetaCS]) {
43044302
values[kCosThetaCS] = zaxis_CS.Dot(v_CM);
4305-
}
4303+
}
43064304
if (fgUsedVars[kPhiCS]) {
43074305
values[kPhiCS] = TMath::ATan2(yaxis_CS.Dot(v_CM), xaxis_CS.Dot(v_CM));
43084306
if (values[kPhiCS] < 0) {
@@ -4366,7 +4364,7 @@ void VarManager::FillPairME(T1 const& t1, T2 const& t2, float* values)
43664364
ROOT::Math::XYZVectorF zaxis_RM(randomCostheta, std::sqrt(1 - randomCostheta * randomCostheta) * std::cos(randomPhi), std::sqrt(1 - randomCostheta * randomCostheta) * std::sin(randomPhi));
43674365
if (fgUsedVars[kCosThetaRM]) {
43684366
values[kCosThetaRM] = zaxis_RM.Dot(v_CM);
4369-
}
4367+
}
43704368
}
43714369
}
43724370

@@ -4544,7 +4542,7 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
45444542
ROOT::Math::XYZVectorF xaxis_HE{(yaxis_HE.Cross(zaxis_HE)).Unit()};
45454543
if (fgUsedVars[kMCCosThetaHE]) {
45464544
values[kMCCosThetaHE] = zaxis_HE.Dot(v_CM);
4547-
}
4545+
}
45484546
if (fgUsedVars[kMCPhiHE]) {
45494547
values[kMCPhiHE] = TMath::ATan2(yaxis_HE.Dot(v_CM), xaxis_HE.Dot(v_CM));
45504548
if (values[kMCPhiHE] < 0) {
@@ -4576,7 +4574,7 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
45764574
ROOT::Math::XYZVectorF xaxis_CS{(yaxis_CS.Cross(zaxis_CS)).Unit()};
45774575
if (fgUsedVars[kMCCosThetaCS]) {
45784576
values[kMCCosThetaCS] = zaxis_CS.Dot(v_CM);
4579-
}
4577+
}
45804578
if (fgUsedVars[kMCPhiCS]) {
45814579
values[kMCPhiCS] = TMath::ATan2(yaxis_CS.Dot(v_CM), xaxis_CS.Dot(v_CM));
45824580
if (values[kMCPhiCS] < 0) {
@@ -4640,7 +4638,7 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
46404638
ROOT::Math::XYZVectorF zaxis_RM(randomCostheta, std::sqrt(1 - randomCostheta * randomCostheta) * std::cos(randomPhi), std::sqrt(1 - randomCostheta * randomCostheta) * std::sin(randomPhi));
46414639
if (fgUsedVars[kMCCosThetaRM]) {
46424640
values[kMCCosThetaRM] = zaxis_RM.Dot(v_CM);
4643-
}
4641+
}
46444642
}
46454643
}
46464644

@@ -4866,7 +4864,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
48664864
v12 = v1 + v2;
48674865
if (fgPVrecalKF) {
48684866
primaryVertexNew = RecalculatePrimaryVertex(t1, t2, collision);
4869-
}
4867+
}
48704868

48714869
} else if constexpr (pairType == kDecayToMuMu && muonHasCov) {
48724870
// Get pca candidate from forward DCA fitter
@@ -4933,7 +4931,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
49334931
values[kVertexingTauxyProjectedNs] = values[kVertexingTauxyProjected] / o2::constants::physics::LightSpeedCm2NS;
49344932
if (fgPVrecalKF) {
49354933
values[kVertexingTauxyProjectedPoleJPsiMassRecalculatePV] = values[kVertexingLxyProjectedRecalculatePV] * o2::constants::physics::MassJPsi / (v12.Pt());
4936-
}
4934+
}
49374935
values[kVertexingTauzProjected] = values[kVertexingLzProjected] * v12.M() / TMath::Abs(v12.Pz());
49384936
values[kVertexingTauxyzProjected] = values[kVertexingLxyzProjected] * v12.M() / (v12.P());
49394937
}
@@ -4977,7 +4975,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
49774975
values[kKFMass] = mass;
49784976
} else {
49794977
values[kKFMass] = -999.;
4980-
}
4978+
}
49814979
}
49824980

49834981
if constexpr (eventHasVtxCov) {
@@ -4996,15 +4994,15 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
49964994
values[kVertexingLxyzErr] = (KFPV.GetCovariance(0) + KFGeoTwoProng.GetCovariance(0)) * dxPair2PV * dxPair2PV + (KFPV.GetCovariance(2) + KFGeoTwoProng.GetCovariance(2)) * dyPair2PV * dyPair2PV + (KFPV.GetCovariance(5) + KFGeoTwoProng.GetCovariance(5)) * dzPair2PV * dzPair2PV + 2 * ((KFPV.GetCovariance(1) + KFGeoTwoProng.GetCovariance(1)) * dxPair2PV * dyPair2PV + (KFPV.GetCovariance(3) + KFGeoTwoProng.GetCovariance(3)) * dxPair2PV * dzPair2PV + (KFPV.GetCovariance(4) + KFGeoTwoProng.GetCovariance(4)) * dyPair2PV * dzPair2PV);
49974995
if (fabs(values[kVertexingLxy]) < 1.e-8f) {
49984996
values[kVertexingLxy] = 1.e-8f;
4999-
}
4997+
}
50004998
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
50014999
if (fabs(values[kVertexingLz]) < 1.e-8f) {
50025000
values[kVertexingLz] = 1.e-8f;
5003-
}
5001+
}
50045002
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
50055003
if (fabs(values[kVertexingLxyz]) < 1.e-8f) {
50065004
values[kVertexingLxyz] = 1.e-8f;
5007-
}
5005+
}
50085006
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
50095007
values[kVertexingTauxy] = KFGeoTwoProng.GetPseudoProperDecayTime(KFPV, KFGeoTwoProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
50105008
values[kVertexingTauz] = -1 * dzPair2PV * KFGeoTwoProng.GetMass() / (TMath::Abs(KFGeoTwoProng.GetPz()) * o2::constants::physics::LightSpeedCm2NS);
@@ -5038,10 +5036,10 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
50385036

50395037
if (fgUsedVars[kKFChi2OverNDFGeo]) {
50405038
values[kKFChi2OverNDFGeo] = KFGeoTwoProng.GetChi2() / KFGeoTwoProng.GetNDF();
5041-
}
5039+
}
50425040
if (fgUsedVars[kKFCosPA]) {
50435041
values[kKFCosPA] = calculateCosPA(KFGeoTwoProng, KFPV);
5044-
}
5042+
}
50455043

50465044
// in principle, they should be in FillTrack
50475045
if (fgUsedVars[kKFTrack0DCAxyz] || fgUsedVars[kKFTrack1DCAxyz]) {
@@ -5054,10 +5052,10 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
50545052
}
50555053
if (fgUsedVars[kKFDCAxyzBetweenProngs]) {
50565054
values[kKFDCAxyzBetweenProngs] = trk0KF.GetDistanceFromParticle(trk1KF);
5057-
}
5055+
}
50585056
if (fgUsedVars[kKFDCAxyBetweenProngs]) {
50595057
values[kKFDCAxyBetweenProngs] = trk0KF.GetDistanceFromParticleXY(trk1KF);
5060-
}
5058+
}
50615059

50625060
if (fgUsedVars[kKFTracksDCAxyzMax]) {
50635061
values[kKFTracksDCAxyzMax] = values[kKFTrack0DCAxyz] > values[kKFTrack1DCAxyz] ? values[kKFTrack0DCAxyz] : values[kKFTrack1DCAxyz];
@@ -5092,7 +5090,7 @@ void VarManager::FillPairVertexing(C const& collision, T const& t1, T const& t2,
50925090
values[kKFMassGeoTop] = mass;
50935091
} else {
50945092
values[kKFMassGeoTop] = -999.;
5095-
}
5093+
}
50965094
}
50975095
if (propToSV) {
50985096
if constexpr ((pairType == kDecayToMuMu) && muonHasCov) {
@@ -5190,7 +5188,7 @@ void VarManager::FillPairVertexingRecomputePV(C const& /*collision*/, T const& t
51905188

51915189
if (fgFitterTwoProngBarrel.getNCandidates() == 0) {
51925190
return;
5193-
}
5191+
}
51945192
Vec3D secondaryVertex;
51955193

51965194
if (!fgUsedKF) { // to be updated when seconday vertex is computed with KF
@@ -5393,7 +5391,7 @@ void VarManager::FillTripletVertexing(C const& collision, T const& t1, T const&
53935391
values[kKFMass] = mass;
53945392
} else {
53955393
values[kKFMass] = -999.;
5396-
}
5394+
}
53975395
}
53985396

53995397
if constexpr (eventHasVtxCov) {
@@ -5413,15 +5411,15 @@ void VarManager::FillTripletVertexing(C const& collision, T const& t1, T const&
54135411
values[kVertexingLxyzErr] = (KFPV.GetCovariance(0) + KFGeoThreeProng.GetCovariance(0)) * dxTriplet2PV * dxTriplet2PV + (KFPV.GetCovariance(2) + KFGeoThreeProng.GetCovariance(2)) * dyTriplet2PV * dyTriplet2PV + (KFPV.GetCovariance(5) + KFGeoThreeProng.GetCovariance(5)) * dzTriplet2PV * dzTriplet2PV + 2 * ((KFPV.GetCovariance(1) + KFGeoThreeProng.GetCovariance(1)) * dxTriplet2PV * dyTriplet2PV + (KFPV.GetCovariance(3) + KFGeoThreeProng.GetCovariance(3)) * dxTriplet2PV * dzTriplet2PV + (KFPV.GetCovariance(4) + KFGeoThreeProng.GetCovariance(4)) * dyTriplet2PV * dzTriplet2PV);
54145412
if (fabs(values[kVertexingLxy]) < 1.e-8f) {
54155413
values[kVertexingLxy] = 1.e-8f;
5416-
}
5414+
}
54175415
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
54185416
if (fabs(values[kVertexingLz]) < 1.e-8f) {
54195417
values[kVertexingLz] = 1.e-8f;
5420-
}
5418+
}
54215419
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
54225420
if (fabs(values[kVertexingLxyz]) < 1.e-8f) {
54235421
values[kVertexingLxyz] = 1.e-8f;
5424-
}
5422+
}
54255423
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
54265424

54275425
values[kVertexingTauxy] = KFGeoThreeProng.GetPseudoProperDecayTime(KFPV, KFGeoThreeProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
@@ -5647,15 +5645,15 @@ void VarManager::FillDileptonTrackVertexing(C const& collision, T1 const& lepton
56475645
// Quantities between 3rd prong and candidate
56485646
if (fgUsedVars[kKFDCAxyzBetweenProngs]) {
56495647
values[kKFDCAxyzBetweenProngs] = KFGeoTwoLeptons.GetDistanceFromParticle(hadronKF);
5650-
}
5648+
}
56515649

56525650
KFGeoThreeProng.SetConstructMethod(2);
56535651
KFGeoThreeProng.AddDaughter(KFGeoTwoLeptons);
56545652
KFGeoThreeProng.AddDaughter(hadronKF);
56555653

56565654
if (fgUsedVars[kKFMass]) {
56575655
values[kKFMass] = KFGeoThreeProng.GetMass();
5658-
}
5656+
}
56595657

56605658
if constexpr (eventHasVtxCov) {
56615659
KFPVertex kfpVertex = createKFPVertexFromCollision(collision);
@@ -5673,30 +5671,30 @@ void VarManager::FillDileptonTrackVertexing(C const& collision, T1 const& lepton
56735671
values[kVertexingLxyzErr] = (KFPV.GetCovariance(0) + KFGeoThreeProng.GetCovariance(0)) * dxTriplet3PV * dxTriplet3PV + (KFPV.GetCovariance(2) + KFGeoThreeProng.GetCovariance(2)) * dyTriplet3PV * dyTriplet3PV + (KFPV.GetCovariance(5) + KFGeoThreeProng.GetCovariance(5)) * dzTriplet3PV * dzTriplet3PV + 2 * ((KFPV.GetCovariance(1) + KFGeoThreeProng.GetCovariance(1)) * dxTriplet3PV * dyTriplet3PV + (KFPV.GetCovariance(3) + KFGeoThreeProng.GetCovariance(3)) * dxTriplet3PV * dzTriplet3PV + (KFPV.GetCovariance(4) + KFGeoThreeProng.GetCovariance(4)) * dyTriplet3PV * dzTriplet3PV);
56745672
if (fabs(values[kVertexingLxy]) < 1.e-8f) {
56755673
values[kVertexingLxy] = 1.e-8f;
5676-
}
5674+
}
56775675
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
56785676
if (fabs(values[kVertexingLz]) < 1.e-8f) {
56795677
values[kVertexingLz] = 1.e-8f;
5680-
}
5678+
}
56815679
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
56825680
if (fabs(values[kVertexingLxyz]) < 1.e-8f) {
56835681
values[kVertexingLxyz] = 1.e-8f;
5684-
}
5682+
}
56855683
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
56865684

56875685
if (fgUsedVars[kVertexingTauxy]) {
56885686
values[kVertexingTauxy] = KFGeoThreeProng.GetPseudoProperDecayTime(KFPV, KFGeoThreeProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
5689-
}
5687+
}
56905688
if (fgUsedVars[kVertexingTauxyErr]) {
56915689
values[kVertexingTauxyErr] = values[kVertexingLxyErr] * KFGeoThreeProng.GetMass() / (KFGeoThreeProng.GetPt() * o2::constants::physics::LightSpeedCm2NS);
5692-
}
5690+
}
56935691

56945692
if (fgUsedVars[kCosPointingAngle]) {
56955693
values[VarManager::kCosPointingAngle] = (dxTriplet3PV * KFGeoThreeProng.GetPx() +
56965694
dyTriplet3PV * KFGeoThreeProng.GetPy() +
56975695
dzTriplet3PV * KFGeoThreeProng.GetPz()) /
56985696
(KFGeoThreeProng.GetP() * values[VarManager::kVertexingLxyz]);
5699-
}
5697+
}
57005698
} // end calculate vertex variables
57015699

57025700
// As defined in Run 2 (projected onto momentum)
@@ -6834,7 +6832,7 @@ void VarManager::FillDileptonTrackTrackVertexing(C const& collision, T1 const& l
68346832
values[kKFMass] = mass;
68356833
} else {
68366834
values[kKFMass] = -999.;
6837-
}
6835+
}
68386836
}
68396837

68406838
KFPVertex kfpVertex = createKFPVertexFromCollision(collision);
@@ -6854,15 +6852,15 @@ void VarManager::FillDileptonTrackTrackVertexing(C const& collision, T1 const& l
68546852

68556853
if (fabs(values[kVertexingLxy]) < 1.e-8f) {
68566854
values[kVertexingLxy] = 1.e-8f;
6857-
}
6855+
}
68586856
values[kVertexingLxyErr] = values[kVertexingLxyErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyErr]) / values[kVertexingLxy];
68596857
if (fabs(values[kVertexingLz]) < 1.e-8f) {
68606858
values[kVertexingLz] = 1.e-8f;
6861-
}
6859+
}
68626860
values[kVertexingLzErr] = values[kVertexingLzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLzErr]) / values[kVertexingLz];
68636861
if (fabs(values[kVertexingLxyz]) < 1.e-8f) {
68646862
values[kVertexingLxyz] = 1.e-8f;
6865-
}
6863+
}
68666864
values[kVertexingLxyzErr] = values[kVertexingLxyzErr] < 0. ? 1.e8f : std::sqrt(values[kVertexingLxyzErr]) / values[kVertexingLxyz];
68676865

68686866
values[kVertexingTauxy] = KFGeoFourProng.GetPseudoProperDecayTime(KFPV, KFGeoFourProng.GetMass()) / (o2::constants::physics::LightSpeedCm2NS);
@@ -7447,10 +7445,10 @@ void VarManager::FillPairAlice3(T1 const& t1, T2 const& t2, float* values)
74477445
double arg = scalar / Ptot12;
74487446
if (arg > 1.) {
74497447
arg = 1.;
7450-
}
7448+
}
74517449
if (arg < -1) {
74527450
arg = -1;
7453-
}
7451+
}
74547452
values[kOpeningAngle] = TMath::ACos(arg);
74557453
}
74567454
}
@@ -7477,7 +7475,7 @@ void VarManager::FillPairAlice3(T1 const& t1, T2 const& t2, float* values)
74777475
ROOT::Math::XYZVectorF xaxis_HE{(yaxis_HE.Cross(zaxis_HE)).Unit()};
74787476
if (fgUsedVars[kCosThetaHE]) {
74797477
values[kCosThetaHE] = zaxis_HE.Dot(v_CM);
7480-
}
7478+
}
74817479
if (fgUsedVars[kPhiHE]) {
74827480
values[kPhiHE] = TMath::ATan2(yaxis_HE.Dot(v_CM), xaxis_HE.Dot(v_CM));
74837481
if (values[kPhiHE] < 0) {
@@ -7509,7 +7507,7 @@ void VarManager::FillPairAlice3(T1 const& t1, T2 const& t2, float* values)
75097507
ROOT::Math::XYZVectorF xaxis_CS{(yaxis_CS.Cross(zaxis_CS)).Unit()};
75107508
if (fgUsedVars[kCosThetaCS]) {
75117509
values[kCosThetaCS] = zaxis_CS.Dot(v_CM);
7512-
}
7510+
}
75137511
if (fgUsedVars[kPhiCS]) {
75147512
values[kPhiCS] = TMath::ATan2(yaxis_CS.Dot(v_CM), xaxis_CS.Dot(v_CM));
75157513
if (values[kPhiCS] < 0) {
@@ -7573,7 +7571,7 @@ void VarManager::FillPairAlice3(T1 const& t1, T2 const& t2, float* values)
75737571
ROOT::Math::XYZVectorF zaxis_RM(randomCostheta, std::sqrt(1 - randomCostheta * randomCostheta) * std::cos(randomPhi), std::sqrt(1 - randomCostheta * randomCostheta) * std::sin(randomPhi));
75747572
if (fgUsedVars[kCosThetaRM]) {
75757573
values[kCosThetaRM] = zaxis_RM.Dot(v_CM);
7576-
}
7574+
}
75777575
}
75787576
}
75797577

0 commit comments

Comments
 (0)