128 {
129
130 uint64_t event = ctx.eventID().event_number();
131 ATH_MSG_DEBUG(
"********************************************************* EVENT NUMBER = " << event);
132
135 ATH_MSG_ERROR(
"Failed to retrieve the MuonDetectorManager conditions object");
136 return StatusCode::FAILURE;
137 }
138
143 if( !readMcEventCollection.isValid() ){
145 return StatusCode::FAILURE;
146 }
147 if(
m_doTruth) ptrMcEventCollection = readMcEventCollection.cptr();
149 if( !readMuonEntryLayer.isValid() ){
151 return StatusCode::FAILURE;
152 }
154 ptrMuonEntryLayer = readMuonEntryLayer.cptr();
155 MMLoadVariables
load = MMLoadVariables();
156 std::map<std::pair<uint64_t, unsigned int>,evInf_entry> Event_Info;
157 ATH_CHECK(
load.getTruthInfo(ctx, ptrMcEventCollection, ptrMuonEntryLayer, Event_Info));
158
159
160 for (const auto &it : Event_Info) {
161 m_trigger_trueEtaRange->push_back(
it.second.eta_ip);
162 m_trigger_truePtRange->push_back(
it.second.pt);
163 m_trigger_trueThe->push_back(
it.second.theta_ip);
164 m_trigger_truePhi->push_back(
it.second.phi_ip);
165 m_trigger_trueDth->push_back(
it.second.dtheta);
166 m_trigger_trueEtaPos->push_back(
it.second.eta_pos);
167 m_trigger_trueThePos->push_back(
it.second.theta_pos);
168 m_trigger_truePhiPos->push_back(
it.second.phi_pos);
169 m_trigger_trueEtaEnt->push_back(
it.second.eta_ent);
170 m_trigger_trueTheEnt->push_back(
it.second.theta_ent);
171 m_trigger_truePhiEnt->push_back(
it.second.phi_ent);
172 }
173 }
174 }
175
177 if( !readMmDigitContainer.isValid() ){
179 return StatusCode::FAILURE;
180 }
181
182 for (const MmDigitCollection* digitCollection : *readMmDigitContainer) {
183
184 std::vector<std::shared_ptr<MMT_Hit> > ev_hits;
185 for (const MmDigit* digit : *digitCollection) {
186 const Identifier
id =
digit->identify();
188
193 const int multiplet =
m_idHelperSvc->mmIdHelper().multilayer(
id);
195
197
198 const MuonGM::MMReadoutElement* readout = detManager->getMMReadoutElement(id);
199 if (channel < 1 or channel > (readout->
getDesign(
id))->totalStrips)
continue;
200
201 const float stripTime =
digit->stripResponseTime();
202
203 if (stripTime < 0.) continue;
204 const int BC = std::ceil(stripTime/25.);
205 ev_hits.emplace_back(std::make_shared<MMT_Hit>(
id, stationName, stationEta, stationPhi, sector, multiplet, gasGap, channel, stripTime, BC, detManager.
cptr()));
206
207 if (ev_hits.back()->infSlope()) {
209 ev_hits.pop_back();
210 continue;
211 }
213 m_trigger_VMM->push_back(ev_hits.back()->getVMM());
214 m_trigger_plane->push_back(ev_hits.back()->getPlane());
215 m_trigger_station->push_back(ev_hits.back()->getStationEta());
216 m_trigger_strip->push_back(ev_hits.back()->getChannel());
217 m_trigger_RZslopes->push_back(ev_hits.back()->getRZSlope());
218 }
219 }
220
221
222 if (do_MMDiamonds and ev_hits.size() >= (
m_diamond->getXthreshold()+
m_diamond->getUVthreshold())) {
223 const bool isLarge = (std::all_of(ev_hits.begin(), ev_hits.end(), [] (const auto &hit) { return hit->getSector() == 'L'; }));
224 const char sector = (isLarge) ? 'L' : 'S';
225 const char side = (std::all_of(ev_hits.begin(), ev_hits.end(), [] (
const auto &hit) { return hit->getStationEta() < 0; })) ?
'C' :
'A';
226 const int sectorPhi = ev_hits[0]->getSectorPhi();
227 const bool allSectorPhi = (std::all_of(ev_hits.begin(), ev_hits.end(), [&] (const auto &hit) { return hit->getSectorPhi() == sectorPhi; }));
228 if (not allSectorPhi) {
229 ATH_MSG_ERROR(
"Available digits belongs to different sectors IDs, unable to assign an unique ID in output RDO");
230 return StatusCode::FAILURE;
231 }
232
233
234 std::vector<MMT_Road> ev_roads;
235 m_diamond->createRoads(ev_roads, isLarge);
236
237
238 std::vector<slope_t> diamondSlopes;
239 m_diamond->findDiamonds(ev_hits, ev_roads, diamondSlopes, sectorPhi);
240
241
242 if (not diamondSlopes.empty()) {
244 m_trigger_diamond_ntrig->push_back(diamondSlopes.size());
245 for (const auto &slope : diamondSlopes) {
246 m_trigger_diamond_sector->push_back(sector);
247 m_trigger_diamond_sectorPhi->push_back(sectorPhi);
248 m_trigger_diamond_bc->push_back(slope.BC);
249 m_trigger_diamond_totalCount->push_back(slope.totalCount);
250 m_trigger_diamond_realCount->push_back(slope.realCount);
251 m_trigger_diamond_XbkgCount->push_back(slope.xbkg);
252 m_trigger_diamond_UVbkgCount->push_back(slope.uvbkg);
253 m_trigger_diamond_XmuonCount->push_back(slope.xmuon);
254 m_trigger_diamond_UVmuonCount->push_back(slope.uvmuon);
255 m_trigger_diamond_iX->push_back(slope.iRoad);
256 m_trigger_diamond_iU->push_back(slope.iRoadu);
257 m_trigger_diamond_iV->push_back(slope.iRoadv);
258 m_trigger_diamond_age->push_back(slope.age);
259 m_trigger_diamond_mx->push_back(slope.mx);
260 m_trigger_diamond_my->push_back(slope.my);
261 m_trigger_diamond_Uavg->push_back(slope.uavg);
262 m_trigger_diamond_Vavg->push_back(slope.vavg);
263 m_trigger_diamond_mxl->push_back(slope.mxl);
264 m_trigger_diamond_theta->push_back(slope.theta);
265 m_trigger_diamond_eta->push_back(slope.eta);
266 m_trigger_diamond_dtheta->push_back(slope.dtheta);
267 m_trigger_diamond_phi->push_back(slope.phi);
268 m_trigger_diamond_phiShf->push_back(slope.phiShf);
269 }
270 }
271
272
273 std::vector<int> slopeBC;
274 for (const auto &slope : diamondSlopes) slopeBC.push_back(slope.BC);
275 std::sort(slopeBC.begin(), slopeBC.end());
276 slopeBC.erase(
std::unique(slopeBC.begin(), slopeBC.end()), slopeBC.end() );
277 for (const auto &bc : slopeBC) {
278 Muon::NSW_TrigRawData* trigRawData = new Muon::NSW_TrigRawData(sectorPhi-1, side, bc);
279
280 for (const auto &slope : diamondSlopes) {
281 if (bc == slope.BC) {
282 Muon::NSW_TrigRawDataSegment* trigRawDataSegment = new Muon::NSW_TrigRawDataSegment();
283
284
287 else {
290 for (uint8_t i=0;
i<
nPhi;
i++) {
291 if ((slope.phi) < (
m_phiMin+i*phiSteps)) {
293 break;
294 }
295 }
297 }
298 if (
m_doNtuple) m_trigger_diamond_TP_phi_id->push_back(phi_id);
299
300
301 double extrapolatedR = 7824.46*std::abs(std::tan(slope.theta));
304 else {
307 for (uint8_t j=0; j<nR; j++) {
308 if (extrapolatedR < (
m_rMin+j*Rsteps)) {
309 R_id = j;
310 break;
311 }
312 }
314 }
315 if (
m_doNtuple) m_trigger_diamond_TP_R_id->push_back(R_id);
316
317
320 else {
323 for (uint8_t k=0;
k<ndTheta;
k++) {
324 if ((slope.dtheta) < (
m_dThetaMin+k*dThetaSteps)) {
326 break;
327 }
328 }
330 }
331 if (
m_doNtuple) m_trigger_diamond_TP_dTheta_id->push_back(dTheta_id);
332
333
334 trigRawDataSegment->
setLowRes(slope.lowRes);
335
336 trigRawData->
push_back(trigRawDataSegment);
337 }
338 }
340 }
341 }
342 }
343 else {
344 ATH_MSG_DEBUG(
"Available hits are " << ev_hits.size() <<
", less than X+UV threshold, skipping digit collection");
345 }
346 }
347
348 return StatusCode::SUCCESS;
349 }
#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
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.
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.