122 {
123
124 uint64_t event = ctx.eventID().event_number();
125 ATH_MSG_DEBUG(
"********************************************************* EVENT NUMBER = " << event);
126
129 ATH_MSG_ERROR(
"Failed to retrieve the MuonDetectorManager conditions object");
130 return StatusCode::FAILURE;
131 }
132
137 if( !readMcEventCollection.isValid() ){
139 return StatusCode::FAILURE;
140 }
141 if(
m_doTruth) ptrMcEventCollection = readMcEventCollection.cptr();
143 if( !readMuonEntryLayer.isValid() ){
145 return StatusCode::FAILURE;
146 }
148 ptrMuonEntryLayer = readMuonEntryLayer.cptr();
149 MMLoadVariables
load = MMLoadVariables();
150 std::map<std::pair<uint64_t, unsigned int>,evInf_entry> Event_Info;
151 ATH_CHECK(
load.getTruthInfo(ctx, ptrMcEventCollection, ptrMuonEntryLayer, Event_Info));
152
153
154 for (const auto &it : Event_Info) {
155 m_trigger_trueEtaRange->push_back(
it.second.eta_ip);
156 m_trigger_truePtRange->push_back(
it.second.pt);
157 m_trigger_trueThe->push_back(
it.second.theta_ip);
158 m_trigger_truePhi->push_back(
it.second.phi_ip);
159 m_trigger_trueDth->push_back(
it.second.dtheta);
160 m_trigger_trueEtaPos->push_back(
it.second.eta_pos);
161 m_trigger_trueThePos->push_back(
it.second.theta_pos);
162 m_trigger_truePhiPos->push_back(
it.second.phi_pos);
163 m_trigger_trueEtaEnt->push_back(
it.second.eta_ent);
164 m_trigger_trueTheEnt->push_back(
it.second.theta_ent);
165 m_trigger_truePhiEnt->push_back(
it.second.phi_ent);
166 }
167 }
168 }
169
171 if( !readMmDigitContainer.isValid() ){
173 return StatusCode::FAILURE;
174 }
175
176 for (const MmDigitCollection* digitCollection : *readMmDigitContainer) {
177
178 std::vector<std::shared_ptr<MMT_Hit> > ev_hits;
179 for (const MmDigit* digit : *digitCollection) {
180 const Identifier
id =
digit->identify();
182
187 const int multiplet =
m_idHelperSvc->mmIdHelper().multilayer(
id);
189
191
192 const MuonGM::MMReadoutElement* readout = detManager->getMMReadoutElement(id);
193 if (channel < 1 or channel > (readout->
getDesign(
id))->totalStrips)
continue;
194
195 const float stripTime =
digit->stripResponseTime();
196
197 if (stripTime < 0.) continue;
198 const int BC = std::ceil(stripTime/25.);
199 ev_hits.emplace_back(std::make_shared<MMT_Hit>(
id, stationName, stationEta, stationPhi, sector, multiplet, gasGap, channel, stripTime, BC, detManager.
cptr()));
200
201 if (ev_hits.back()->infSlope()) {
203 ev_hits.pop_back();
204 continue;
205 }
207 m_trigger_VMM->push_back(ev_hits.back()->getVMM());
208 m_trigger_plane->push_back(ev_hits.back()->getPlane());
209 m_trigger_station->push_back(ev_hits.back()->getStationEta());
210 m_trigger_strip->push_back(ev_hits.back()->getChannel());
211 m_trigger_RZslopes->push_back(ev_hits.back()->getRZSlope());
212 }
213 }
214
215
216 if (do_MMDiamonds and ev_hits.size() >= (
m_diamond->getXthreshold()+
m_diamond->getUVthreshold())) {
217 const bool isLarge = (std::ranges::all_of(ev_hits, [](const auto &hit) { return hit->getSector() == 'L'; }));
218 const char sector = (isLarge) ? 'L' : 'S';
219 const char side = (std::ranges::all_of(ev_hits, [](
const auto &hit) {
return hit->getStationEta() < 0; })) ?
'C' :
'A';
220 const int sectorPhi = ev_hits[0]->getSectorPhi();
221 const bool allSectorPhi = (std::ranges::all_of(ev_hits, [&](const auto &hit) { return hit->getSectorPhi() == sectorPhi; }));
222 if (not allSectorPhi) {
223 ATH_MSG_ERROR(
"Available digits belongs to different sectors IDs, unable to assign an unique ID in output RDO");
224 return StatusCode::FAILURE;
225 }
226 const bool isEta1 = (std::ranges::all_of(ev_hits, [](const auto &hit) { return std::abs(hit->getStationEta()) == 1; }));
227
228
229 std::vector<MMT_Road> ev_roads;
230 m_diamond->createRoads(ev_roads, isLarge, isEta1);
231
232
233 std::vector<slope_t> diamondSlopes;
234 m_diamond->findDiamonds(ev_hits, ev_roads, diamondSlopes, sectorPhi);
235
236
237 if (not diamondSlopes.empty()) {
239 m_trigger_diamond_ntrig->push_back(diamondSlopes.size());
240 for (const auto &slope : diamondSlopes) {
241 m_trigger_diamond_sector->push_back(sector);
242 m_trigger_diamond_sectorPhi->push_back(sectorPhi);
243 m_trigger_diamond_bc->push_back(slope.BC);
244 m_trigger_diamond_totalCount->push_back(slope.totalCount);
245 m_trigger_diamond_xCount->push_back(slope.xCount);
246 m_trigger_diamond_uCount->push_back(slope.uCount);
247 m_trigger_diamond_iX->push_back(slope.iRoad);
248 m_trigger_diamond_iU->push_back(slope.iRoadu);
249 m_trigger_diamond_iV->push_back(slope.iRoadv);
250 m_trigger_diamond_age->push_back(slope.age);
251 m_trigger_diamond_mx->push_back(slope.mx);
252 m_trigger_diamond_my->push_back(slope.my);
253 m_trigger_diamond_Uavg->push_back(slope.uavg);
254 m_trigger_diamond_Vavg->push_back(slope.vavg);
255 m_trigger_diamond_mxl->push_back(slope.mxl);
256 m_trigger_diamond_theta->push_back(slope.theta);
257 m_trigger_diamond_eta->push_back(slope.eta);
258 m_trigger_diamond_dtheta->push_back(slope.dtheta);
259 m_trigger_diamond_phi->push_back(slope.phi);
260 m_trigger_diamond_phiShf->push_back(slope.phiShf);
261 }
262 }
263
264
265 std::vector<int> slopeBC;
266 for (const auto &slope : diamondSlopes) slopeBC.push_back(slope.BC);
267 std::ranges::sort(slopeBC);
268 slopeBC.erase(
std::unique(slopeBC.begin(), slopeBC.end()), slopeBC.end() );
269 for (const auto bc : slopeBC) {
270 Muon::NSW_TrigRawData* trigRawData = new Muon::NSW_TrigRawData(sectorPhi-1, side, bc);
271
272 for (const auto &slope : diamondSlopes) {
273 if (bc == slope.BC) {
274 Muon::NSW_TrigRawDataSegment* trigRawDataSegment = new Muon::NSW_TrigRawDataSegment();
275
276
279 else {
282 for (uint8_t i=0;
i<
nPhi;
i++) {
283 if ((slope.phi) < (
m_phiMin+i*phiSteps)) {
285 break;
286 }
287 }
289 }
290 if (
m_doNtuple) m_trigger_diamond_TP_phi_id->push_back(phi_id);
291
292
293 double extrapolatedR = 7824.46*std::abs(std::tan(slope.theta));
296 else {
299 for (uint8_t j=0;
j<nR;
j++) {
300 if (extrapolatedR < (
m_rMin+j*Rsteps)) {
302 break;
303 }
304 }
306 }
307 if (
m_doNtuple) m_trigger_diamond_TP_R_id->push_back(R_id);
308
309
312 else {
315 for (uint8_t k=0;
k<ndTheta;
k++) {
316 if ((slope.dtheta) < (
m_dThetaMin+k*dThetaSteps)) {
318 break;
319 }
320 }
322 }
323 if (
m_doNtuple) m_trigger_diamond_TP_dTheta_id->push_back(dTheta_id);
324
325
326 trigRawDataSegment->
setLowRes(slope.lowRes);
327
328 trigRawData->
push_back(trigRawDataSegment);
329 }
330 }
332 }
333 }
334 }
335 else {
336 ATH_MSG_DEBUG(
"Available hits are " << ev_hits.size() <<
", less than X+UV threshold, skipping digit collection");
337 }
338 }
339
340 return StatusCode::SUCCESS;
341 }
#define ATH_MSG_WARNING(x)
AtlasHitsVector< TrackRecord > TrackRecordCollection
value_type push_back(value_type pElem)
Add an element to the end of the collection.
const MuonChannelDesign * getDesign(const Identifier &id) const
returns the MuonChannelDesign class for the given identifier
void setLowRes(bool lowRes)
void setPhiIndex(uint8_t phiIndex)
void setDeltaTheta(uint8_t deltaTheta)
void setRIndex(uint8_t rIndex)
const_pointer_type cptr()
constexpr uint8_t stationPhi
station Phi 1 to 8
constexpr uint8_t stationEta
1 to 3
float j(const xAOD::IParticle &, const xAOD::TrackMeasurementValidation &hit, const Eigen::Matrix3d &jab_inv)
constexpr int nPhi
Default bin number of phi for vertex map.
load(f, use_proxy=1, key=None)
DataModel_detail::iterator< DVL > unique(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of unique for DataVector/List.