73{
75
76
77 std::vector<const xAOD::Muon*> candidateSAmuons;
79 bool isSA = (
muon->muonType() == xAOD::Muon::MuonStandAlone);
80 bool isCalo = (
muon->muonType() == xAOD::Muon::CaloTagged);
81 bool isSegment = (
muon->muonType() == xAOD::Muon::SegmentTagged);
82 bool isSiForward = (
muon->muonType() == xAOD::Muon::SiliconAssociatedForwardMuon);
83
85 continue;
86 }
87
89 if (!MuSAMSTP) {
90 if (isCalo || isSegment || isSiForward) {
91 ATH_MSG_VERBOSE(
"Skipping non-SA, non-Combined muon type in validation mode!");
92 } else {
94 }
95 continue;
96 }
97
98
101 continue;
102 }
103
104 if (MuSAMSTP->
pt() > 13000000) {
105 ATH_MSG_DEBUG(
"Skipping SA muon with pT " << (MuSAMSTP->
pt() / 1000.) <<
" GeV!");
106 continue;
107 }
108
109
111 muon->parameter(spectrometerFieldIntegral, xAOD::Muon::spectrometerFieldIntegral);
112 if (spectrometerFieldIntegral < 0.1) {
113 ATH_MSG_DEBUG(
"Skipping SA muon with spectrometerFieldIntegral " <<
muon->spectrometerFieldIntegral <<
" T*m!");
114 continue;
115 }
116
117 candidateSAmuons.push_back(muon);
118 }
119
120
121 if (candidateSAmuons.size() < 2) {
123 return StatusCode::SUCCESS;
124 }
125
126
127 std::vector<std::shared_ptr<xAOD::TrackParticle>> extrapolatedMuSATracks;
128 std::map<const xAOD::TrackParticle*, const xAOD::Muon*> muonToExtrapolatedTrackMap;
129 for (const auto muon : candidateSAmuons) {
131 auto extrapolatedMuSATrack =
extrapolateMuSA(*MuSAMSTP, eventInfo, ctx);
132 if (extrapolatedMuSATrack->definingParameters()[
Trk::d0] > 8000) {
133 ATH_MSG_DEBUG(
"Failed to extrapolate MuSA track, skipping!");
134 continue;
135 }
136
137 extrapolatedMuSATracks.emplace_back(extrapolatedMuSATrack.release());
138 muonToExtrapolatedTrackMap[extrapolatedMuSATracks.back().get()] =
muon;
139 ATH_MSG_VERBOSE(
"Extrapolated MuSA track! Total extrapolated so far: " << extrapolatedMuSATracks.size());
140 }
141
142 if (extrapolatedMuSATracks.size() < 2) {
144 return StatusCode::SUCCESS;
145 }
146
147 std::unique_ptr<Trk::IVKalState> state =
m_vertexFitter->makeState(ctx);
148 std::vector<const xAOD::TrackParticle*> tracksToFit(2);
149 std::vector<const xAOD::NeutralParticle*> dummyNeutrals;
150
151 for (
unsigned int i = 0;
i < extrapolatedMuSATracks.size();
i++) {
152 for (
unsigned int j = i+1;
j < extrapolatedMuSATracks.size();
j++) {
153 MuSAVtxFitterTool::WrkVrt MuSACandidate;
154 tracksToFit[0] = extrapolatedMuSATracks[
i].get();
155 tracksToFit[1] = extrapolatedMuSATracks[
j].get();
157 MuSACandidate.pos, MuSACandidate.mom,
158 MuSACandidate.charge, MuSACandidate.cov,
159 MuSACandidate.chi2PerTrk, MuSACandidate.trkAtVrt,
160 MuSACandidate.chi2, *state);
161 if (
res.isSuccess()) {
163 MuSACandidate.newExtrapolatedTracks.push_back(extrapolatedMuSATracks[i]);
164 MuSACandidate.newExtrapolatedTracks.push_back(extrapolatedMuSATracks[j]);
165 MuSACandidate.muonCandidates = {
166 muonToExtrapolatedTrackMap[tracksToFit[0]],
167 muonToExtrapolatedTrackMap[tracksToFit[1]]
168 };
169 MuSACandidate.minOpAng = muonToExtrapolatedTrackMap[tracksToFit[0]]->p4().DeltaR(
170 muonToExtrapolatedTrackMap[tracksToFit[1]]->p4());
171 workVerticesContainer.emplace_back(std::move(MuSACandidate));
172 } else {
174 }
175 }
176 }
177
178
179 return StatusCode::SUCCESS;
180}
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
std::pair< std::vector< unsigned int >, bool > res
xAOD::MuonContainer * muonContainer
virtual double pt() const override final
The transverse momentum ( ) of the particle.
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
::StatusCode StatusCode
StatusCode definition for legacy code.
@ spectrometerFieldIntegral
Discriminators and further variables.
float j(const xAOD::IParticle &, const xAOD::TrackMeasurementValidation &hit, const Eigen::Matrix3d &jab_inv)
TrackParticle_v1 TrackParticle
Reference the current persistent version: