|
22 | 22 | #include "CommonConstants/LHCConstants.h" |
23 | 23 | #include "Framework/AnalysisDataModel.h" |
24 | 24 |
|
25 | | -#include "TLorentzVector.h" |
26 | | - |
27 | 25 | using BCsWithBcSels = o2::soa::Join<o2::aod::BCs, o2::aod::BcSels>; |
28 | 26 |
|
29 | 27 | using ForwardTracks = o2::soa::Join<o2::aod::FwdTracks, o2::aod::FwdTracksCov>; |
@@ -96,13 +94,15 @@ struct FITInfo { |
96 | 94 | int32_t distClosestBcT0A = 999; |
97 | 95 | }; |
98 | 96 |
|
| 97 | +constexpr double AbsorberMid = 26.5; |
| 98 | + |
99 | 99 | template <typename T, typename TSelectorsArray> |
100 | 100 | void applyFwdCuts(UPCCutparHolder& upcCuts, const T& track, TSelectorsArray& fwdSelectors) |
101 | 101 | { |
102 | 102 | fwdSelectors[kFwdSelPt] = track.pt() > upcCuts.getFwdPtLow() && track.pt() < upcCuts.getFwdPtHigh(); // check pt |
103 | 103 | fwdSelectors[kFwdSelEta] = track.eta() > upcCuts.getFwdEtaLow() && track.eta() < upcCuts.getFwdEtaHigh(); // check pseudorapidity |
104 | 104 | fwdSelectors[kFwdSelRabs] = track.rAtAbsorberEnd() > upcCuts.getMuonRAtAbsorberEndLow() && track.rAtAbsorberEnd() < upcCuts.getMuonRAtAbsorberEndHigh(); // check muon R |
105 | | - fwdSelectors[kFwdSelpDCA] = track.rAtAbsorberEnd() < 26.5 ? track.pDca() < upcCuts.getMuonPDcaHighFirst() : track.pDca() < upcCuts.getMuonPDcaHighSecond(); // check pDCA |
| 105 | + fwdSelectors[kFwdSelpDCA] = track.rAtAbsorberEnd() < AbsorberMid ? track.pDca() < upcCuts.getMuonPDcaHighFirst() : track.pDca() < upcCuts.getMuonPDcaHighSecond(); // check pDCA |
106 | 106 | fwdSelectors[kFwdSelChi2] = track.chi2() > upcCuts.getFwdChi2Low() && track.chi2() < upcCuts.getFwdChi2High(); // check chi2 |
107 | 107 | } |
108 | 108 |
|
@@ -131,7 +131,7 @@ void applyBarrelCuts(UPCCutparHolder& upcCuts, const T& track, TSelectorsArray& |
131 | 131 |
|
132 | 132 | if (upcCuts.getCheckMaxDcaXY()) { |
133 | 133 | float dca = track.dcaXY(); |
134 | | - float maxDCA = 0.0105f + 0.0350f / pow(track.pt(), 1.1f); |
| 134 | + float maxDCA = 0.0105f + 0.0350f / std::pow(track.pt(), 1.1f); |
135 | 135 | barrelSelectors[kBarrelSelDCAXY] = dca < maxDCA; |
136 | 136 | } else { |
137 | 137 | barrelSelectors[kBarrelSelDCAXY] = true; |
|
0 commit comments