@@ -30,17 +30,17 @@ using namespace constants::math;
3030
3131class UPCPairCuts
3232{
33- public:
34- enum Particle { Photon = 0 ,
33+ public:
34+ enum Particle { Photon = 0 ,
3535 K0 ,
3636 Lambda,
3737 Phi,
3838 Rho,
3939 ParticlesLastEntry };
4040
41- void SetHistogramRegistry (HistogramRegistry* registry) { histogramRegistry = registry; }
41+ void SetHistogramRegistry (HistogramRegistry* registry) { histogramRegistry = registry; }
4242
43- void SetPairCut (Particle particle, float cut)
43+ void SetPairCut (Particle particle, float cut)
4444 {
4545 LOGF (info, " Enabled pair cut for %d with value %f" , static_cast <int >(particle), cut);
4646 mCuts [particle] = cut;
@@ -66,11 +66,12 @@ class UPCPairCuts
6666
6767 template <typename T>
6868 bool twoTrackCut (T const & track1, T const & track2);
69- protected:
69+
70+ protected:
7071 float mCuts [ParticlesLastEntry] = {-1 };
7172 float mTwoTrackDistance = -1 ; // distance below which the pair is flagged as to be removed
7273 float mTwoTrackRadius = 0 .8f ; // radius at which the two track cuts are applied
73- int magField = 5 ; // magField: B field in kG
74+ int magField = 5 ; // magField: B field in kG
7475
7576 HistogramRegistry* histogramRegistry = nullptr ; // if set, control histograms are stored here
7677
@@ -149,7 +150,7 @@ bool UPCPairCuts::twoTrackCut(T const& track1, T const& track2)
149150 }
150151
151152 if (dphistarminabs < mTwoTrackDistance && std::fabs (deta) < mTwoTrackDistance ) {
152- // LOGF(debug, "Removed track pair %ld %ld with %f %f %f %f %d %f %f %d %d", track1.index(), track2.index(), deta, dphistarminabs, track1.phi2(), track1.pt(), track1.sign(), track2.phi2(), track2.pt(), track2.sign(), magField);
153+ // LOGF(debug, "Removed track pair %ld %ld with %f %f %f %f %d %f %f %d %d", track1.index(), track2.index(), deta, dphistarminabs, track1.phi2(), track1.pt(), track1.sign(), track2.phi2(), track2.pt(), track2.sign(), magField);
153154 return true ;
154155 }
155156
@@ -165,7 +166,7 @@ bool UPCPairCuts::twoTrackCut(T const& track1, T const& track2)
165166template <typename T>
166167bool UPCPairCuts::conversionCut (T const & track1, T const & track2, Particle conv, double cut)
167168{
168- // LOGF(info, "pt is %f %f", track1.pt(), track2.pt());
169+ // LOGF(info, "pt is %f %f", track1.pt(), track2.pt());
169170
170171 if (cut < 0 ) {
171172 return false ;
@@ -306,7 +307,7 @@ double UPCPairCuts::getInvMassSquaredFast(T const& track1, double m0_1, T const&
306307
307308 double mass2 = m0_1 * m0_1 + m0_2 * m0_2 + 2 .0f * (std::sqrt (e1squ * e2squ) - (pt1 * pt2 * (cosDeltaPhi + 1 .0f / tantheta1 / tantheta2)));
308309
309- // LOGF(debug, "%f %f %f %f %f %f %f %f %f", pt1, eta1, phi1, pt2, eta2, phi2, m0_1, m0_2, mass2);
310+ // LOGF(debug, "%f %f %f %f %f %f %f %f %f", pt1, eta1, phi1, pt2, eta2, phi2, m0_1, m0_2, mass2);
310311
311312 return mass2;
312313}
@@ -341,4 +342,4 @@ float UPCPairCuts::getDPhiStar(T const& track1, T const& track2, float radius, i
341342 return dphistar;
342343}
343344
344- #endif
345+ #endif
0 commit comments