|
13 | 13 | // Task performing forward track DCA computation |
14 | 14 | // |
15 | 15 |
|
16 | | -#include "Common/DataModel/TrackSelectionTables.h" |
17 | 16 | #include "Common/Core/fwdtrackUtilities.h" |
| 17 | +#include "Common/DataModel/TrackSelectionTables.h" |
| 18 | + |
18 | 19 | #include "CCDB/BasicCCDBManager.h" |
19 | 20 | #include "DataFormatsParameters/GRPMagField.h" |
20 | 21 | #include "DetectorsBase/GeometryManager.h" |
21 | 22 | #include "DetectorsBase/Propagator.h" |
22 | 23 | #include "GlobalTracking/MatchGlobalFwd.h" |
23 | | - |
24 | 24 | #include <Framework/AnalysisDataModel.h> |
25 | 25 | #include <Framework/AnalysisHelpers.h> |
26 | 26 | #include <Framework/AnalysisTask.h> |
@@ -81,21 +81,21 @@ struct FwdTrackExtension { |
81 | 81 | } |
82 | 82 | const float zField = grpmag->getNominalL3Field(); |
83 | 83 | for (auto& track : tracks) { |
84 | | - const auto trackType = track.trackType(); |
85 | | - o2::dataformats::GlobalFwdTrack fwdtrack = o2::aod::fwdtrackutils::getTrackParCovFwd(track,track); |
86 | | - if (fRefitGlobalMuon && (trackType == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack || trackType == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalForwardTrack)) { |
87 | | - auto muontrack = track.template matchMCHTrack_as<MuonsWithCov>(); |
88 | | - auto mfttrack = track.template matchMFTTrack_as<aod::MFTTracks>(); |
89 | | - o2::dataformats::GlobalFwdTrack propmuon = o2::aod::fwdtrackutils::propagateMuon(muontrack, muontrack, collision, o2::aod::fwdtrackutils::propagationPoint::kToVertex, 0.f, zField); |
90 | | - SMatrix5 tpars(mfttrack.x(), mfttrack.y(), mfttrack.phi(), mfttrack.tgl(), mfttrack.signed1Pt()); |
91 | | - SMatrix55 tcovs{}; |
92 | | - o2::track::TrackParCovFwd mft{mfttrack.z(), tpars, tcovs, mfttrack.chi2()}; |
93 | | - fwdtrack = o2::aod::fwdtrackutils::refitGlobalMuonCov(propmuon,mft); |
94 | | - } |
95 | | - const auto proptrack = o2::aod::fwdtrackutils::propagateTrackParCovFwd(fwdtrack, trackType, collision, o2::aod::fwdtrackutils::propagationPoint::kToDCA, 0.f, zField); |
96 | | - const float dcaX = (proptrack.getX() - collision.posX()); |
97 | | - const float dcaY = (proptrack.getY() - collision.posY()); |
98 | | - fwdDCA(dcaX, dcaY); |
| 84 | + const auto trackType = track.trackType(); |
| 85 | + o2::dataformats::GlobalFwdTrack fwdtrack = o2::aod::fwdtrackutils::getTrackParCovFwd(track, track); |
| 86 | + if (fRefitGlobalMuon && (trackType == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalMuonTrack || trackType == o2::aod::fwdtrack::ForwardTrackTypeEnum::GlobalForwardTrack)) { |
| 87 | + auto muontrack = track.template matchMCHTrack_as<MuonsWithCov>(); |
| 88 | + auto mfttrack = track.template matchMFTTrack_as<aod::MFTTracks>(); |
| 89 | + o2::dataformats::GlobalFwdTrack propmuon = o2::aod::fwdtrackutils::propagateMuon(muontrack, muontrack, collision, o2::aod::fwdtrackutils::propagationPoint::kToVertex, 0.f, zField); |
| 90 | + SMatrix5 tpars(mfttrack.x(), mfttrack.y(), mfttrack.phi(), mfttrack.tgl(), mfttrack.signed1Pt()); |
| 91 | + SMatrix55 tcovs{}; |
| 92 | + o2::track::TrackParCovFwd mft{mfttrack.z(), tpars, tcovs, mfttrack.chi2()}; |
| 93 | + fwdtrack = o2::aod::fwdtrackutils::refitGlobalMuonCov(propmuon, mft); |
| 94 | + } |
| 95 | + const auto proptrack = o2::aod::fwdtrackutils::propagateTrackParCovFwd(fwdtrack, trackType, collision, o2::aod::fwdtrackutils::propagationPoint::kToDCA, 0.f, zField); |
| 96 | + const float dcaX = (proptrack.getX() - collision.posX()); |
| 97 | + const float dcaY = (proptrack.getY() - collision.posY()); |
| 98 | + fwdDCA(dcaX, dcaY); |
99 | 99 | } |
100 | 100 | } |
101 | 101 | }; |
|
0 commit comments