185 return StatusCode::SUCCESS;
237 for (std::pair<
const xAOD::IParticle*,
const std::vector<const CaloCell*>> particleCellsPair : muonCellsMap) {
240 const std::vector<const CaloCell*>&
cells = particleCellsPair.second;
242 std::vector< float > cells_energy;
243 std::vector< float > cells_et;
244 std::vector< float > cells_eta;
245 std::vector< float > cells_phi;
247 std::vector< int > cells_gain;
248 std::vector< bool > cells_bad;
249 std::vector< int > cells_sampling;
250 std::vector< float > cells_time;
251 std::vector< unsigned short > cells_quality;
253 std::vector< float > cells_sinTh;
254 std::vector< float > cells_cosTh;
255 std::vector< float > cells_cotTh;
256 std::vector< float > cells_x;
257 std::vector< float > cells_y;
258 std::vector< float > cells_z;
260 std::vector< float > cells_r;
261 std::vector< float > cells_dx;
262 std::vector< float > cells_dy;
263 std::vector< float > cells_dz;
264 std::vector< float > cells_dr;
265 std::vector< float > cells_volume;
266 std::vector< float > cells_dphi;
267 std::vector< float > cells_deta;
269 std::vector< int > cells_side;
270 std::vector< unsigned short > cells_section;
271 std::vector< unsigned short > cells_module;
272 std::vector< unsigned short > cells_tower;
273 std::vector< unsigned short > cells_sample;
275 std::vector< unsigned short > cells_pmt1_ros;
276 std::vector< unsigned short > cells_pmt2_ros;
277 std::vector< unsigned short > cells_pmt1_drawer;
278 std::vector< unsigned short > cells_pmt2_drawer;
279 std::vector< unsigned short > cells_pmt1_channel;
280 std::vector< unsigned short > cells_pmt2_channel;
282 std::vector< float > cells_pmt1_energy;
283 std::vector< float > cells_pmt2_energy;
284 std::vector< float > cells_pmt1_time;
285 std::vector< float > cells_pmt2_time;
287 std::vector< unsigned short > cells_pmt1_quality;
288 std::vector< unsigned short > cells_pmt2_quality;
289 std::vector< unsigned short > cells_pmt1_qbit;
290 std::vector< unsigned short > cells_pmt2_qbit;
292 std::vector< bool > cells_pmt1_bad;
293 std::vector< bool > cells_pmt2_bad;
294 std::vector< unsigned short > cells_pmt1_gain;
295 std::vector< unsigned short > cells_pmt2_gain;
303 cells_energy.push_back(
cell->energy() );
304 cells_et.push_back(
cell->et() );
305 cells_eta.push_back(
cell->eta() );
306 cells_phi.push_back(
cell->phi() );
309 cells_gain.push_back(
cell->gain() );
310 cells_bad.push_back(
cell->badcell() );
311 cells_sampling.push_back( cell_dde->
getSampling() );
312 cells_time.push_back(
cell->time() );
313 cells_quality.push_back(
cell->quality() );
316 cells_sinTh.push_back(
cell->sinTh() );
317 cells_cosTh.push_back(
cell->cosTh() );
318 cells_cotTh.push_back(
cell->cotTh() );
319 cells_x.push_back(
cell->x() );
320 cells_y.push_back(
cell->y() );
321 cells_z.push_back(
cell->z() );
324 cells_r.push_back( cell_dde->
r() );
325 cells_dx.push_back( cell_dde->
dx() );
326 cells_dy.push_back( cell_dde->
dy() );
327 cells_dz.push_back( cell_dde->
dz() );
328 cells_dr.push_back( cell_dde->
dr() );
329 cells_volume.push_back( cell_dde->
volume() );
330 cells_dphi.push_back( cell_dde->
dphi() );
331 cells_deta.push_back( cell_dde->
deta() );
364 cells_pmt1_ros.push_back( ros1 );
365 cells_pmt2_ros.push_back( ros2 );
366 cells_pmt1_drawer.push_back( drawer1 );
367 cells_pmt2_drawer.push_back( drawer2 );
368 cells_pmt1_channel.push_back( channel1 );
369 cells_pmt2_channel.push_back( channel2 );
371 cells_pmt1_energy.push_back( tile_cell->
ene1() );
372 cells_pmt2_energy.push_back( tile_cell->
ene2() );
373 cells_pmt1_time.push_back( tile_cell->
time1() );
374 cells_pmt2_time.push_back( tile_cell->
time2() );
375 cells_pmt1_quality.push_back( tile_cell->
qual1() );
376 cells_pmt2_quality.push_back( tile_cell->
qual2() );
377 cells_pmt1_qbit.push_back( tile_cell->
qbit1() );
378 cells_pmt2_qbit.push_back( tile_cell->
qbit2() );
379 cells_pmt1_bad.push_back( tile_cell->
badch1() );
380 cells_pmt2_bad.push_back( tile_cell->
badch2() );
381 cells_pmt1_gain.push_back( tile_cell->
gain1() );
382 cells_pmt2_gain.push_back( tile_cell->
gain2() );
385 cellsEnergy(*
particle) = std::move(cells_energy);
386 cellsEt(*
particle) = std::move(cells_et);
387 cellsEta(*
particle) = std::move(cells_eta);
388 cellsPhi(*
particle) = std::move(cells_phi);
389 cellsGain(*
particle) = std::move(cells_gain);
390 cellsBad(*
particle) = std::move(cells_bad);
391 cellsSampling(*
particle) = std::move(cells_sampling);
392 cellsTime(*
particle) = std::move(cells_time);
393 cellsQuality(*
particle) = std::move(cells_quality);
395 cellsSinTh(*
particle) = std::move(cells_sinTh);
396 cellsCosTh(*
particle) = std::move(cells_cosTh);
397 cellsCotTh(*
particle) = std::move(cells_cotTh);
398 cellsX(*
particle) = std::move(cells_x);
399 cellsY(*
particle) = std::move(cells_y);
400 cellsZ(*
particle) = std::move(cells_z);
401 cellsR(*
particle) = std::move(cells_r);
402 cellsDx(*
particle) = std::move(cells_dx);
403 cellsDy(*
particle) = std::move(cells_dy);
404 cellsDz(*
particle) = std::move(cells_dz);
405 cellsDr(*
particle) = std::move(cells_dr);
406 cellsVolume(*
particle) = std::move(cells_volume);
407 cellsDeta(*
particle) = std::move(cells_deta);
408 cellsDphi(*
particle) = std::move(cells_dphi);
409 cellsSide(*
particle) = std::move(cells_side);
410 cellsSection(*
particle) = std::move(cells_section);
411 cellsModule(*
particle) = std::move(cells_module);
412 cellsTower(*
particle) = std::move(cells_tower);
413 cellsSample(*
particle) = std::move(cells_sample);
415 cellsPmt1Ros(*
particle) = std::move(cells_pmt1_ros);
416 cellsPmt2Ros(*
particle) = std::move(cells_pmt2_ros);
417 cellsPmt1Drawer(*
particle) = std::move(cells_pmt1_drawer);
418 cellsPmt2Drawer(*
particle) = std::move(cells_pmt2_drawer);
419 cellsPmt1Channel(*
particle) = std::move(cells_pmt1_channel);
420 cellsPmt2Channel(*
particle) = std::move(cells_pmt2_channel);
421 cellsPmt1Energy(*
particle) = std::move(cells_pmt1_energy);
422 cellsPmt2Energy(*
particle) = std::move(cells_pmt2_energy);
423 cellsPmt1Time(*
particle) = std::move(cells_pmt1_time);
424 cellsPmt2Time(*
particle) = std::move(cells_pmt2_time);
425 cellsPmt1Quality(*
particle) = std::move(cells_pmt1_quality);
426 cellsPmt2Quality(*
particle) = std::move(cells_pmt2_quality);
427 cellsPmt1_Qbit(*
particle) = std::move(cells_pmt1_qbit);
428 cellsPmt2_Qbit(*
particle) = std::move(cells_pmt2_qbit);
429 cellsPmt1Bad(*
particle) = std::move(cells_pmt1_bad);
430 cellsPmt2Bad(*
particle) = std::move(cells_pmt2_bad);
431 cellsPmt1Gain(*
particle) = std::move(cells_pmt1_gain);
432 cellsPmt2Gain(*
particle) = std::move(cells_pmt2_gain);
436 return StatusCode::SUCCESS;