94 return StatusCode::SUCCESS;
146 for (std::pair<
const xAOD::IParticle*,
const std::vector<const CaloCell*>> particleCellsPair : muonCellsMap) {
149 const std::vector<const CaloCell*>&
cells = particleCellsPair.second;
151 std::vector< float > cells_energy;
152 std::vector< float > cells_et;
153 std::vector< float > cells_eta;
154 std::vector< float > cells_phi;
156 std::vector< int > cells_gain;
157 std::vector< bool > cells_bad;
158 std::vector< int > cells_sampling;
159 std::vector< float > cells_time;
160 std::vector< unsigned short > cells_quality;
162 std::vector< float > cells_sinTh;
163 std::vector< float > cells_cosTh;
164 std::vector< float > cells_cotTh;
165 std::vector< float > cells_x;
166 std::vector< float > cells_y;
167 std::vector< float > cells_z;
169 std::vector< float > cells_r;
170 std::vector< float > cells_dx;
171 std::vector< float > cells_dy;
172 std::vector< float > cells_dz;
173 std::vector< float > cells_dr;
174 std::vector< float > cells_volume;
175 std::vector< float > cells_dphi;
176 std::vector< float > cells_deta;
178 std::vector< int > cells_side;
179 std::vector< unsigned short > cells_section;
180 std::vector< unsigned short > cells_module;
181 std::vector< unsigned short > cells_tower;
182 std::vector< unsigned short > cells_sample;
184 std::vector< unsigned short > cells_pmt1_ros;
185 std::vector< unsigned short > cells_pmt2_ros;
186 std::vector< unsigned short > cells_pmt1_drawer;
187 std::vector< unsigned short > cells_pmt2_drawer;
188 std::vector< unsigned short > cells_pmt1_channel;
189 std::vector< unsigned short > cells_pmt2_channel;
191 std::vector< float > cells_pmt1_energy;
192 std::vector< float > cells_pmt2_energy;
193 std::vector< float > cells_pmt1_time;
194 std::vector< float > cells_pmt2_time;
196 std::vector< unsigned short > cells_pmt1_quality;
197 std::vector< unsigned short > cells_pmt2_quality;
198 std::vector< unsigned short > cells_pmt1_qbit;
199 std::vector< unsigned short > cells_pmt2_qbit;
201 std::vector< bool > cells_pmt1_bad;
202 std::vector< bool > cells_pmt2_bad;
203 std::vector< unsigned short > cells_pmt1_gain;
204 std::vector< unsigned short > cells_pmt2_gain;
212 cells_energy.push_back(
cell->energy() );
213 cells_et.push_back(
cell->et() );
214 cells_eta.push_back(
cell->eta() );
215 cells_phi.push_back(
cell->phi() );
218 cells_gain.push_back(
cell->gain() );
219 cells_bad.push_back(
cell->badcell() );
220 cells_sampling.push_back( cell_dde->
getSampling() );
221 cells_time.push_back(
cell->time() );
222 cells_quality.push_back(
cell->quality() );
225 cells_sinTh.push_back(
cell->sinTh() );
226 cells_cosTh.push_back(
cell->cosTh() );
227 cells_cotTh.push_back(
cell->cotTh() );
228 cells_x.push_back(
cell->x() );
229 cells_y.push_back(
cell->y() );
230 cells_z.push_back(
cell->z() );
233 cells_r.push_back( cell_dde->
r() );
234 cells_dx.push_back( cell_dde->
dx() );
235 cells_dy.push_back( cell_dde->
dy() );
236 cells_dz.push_back( cell_dde->
dz() );
237 cells_dr.push_back( cell_dde->
dr() );
238 cells_volume.push_back( cell_dde->
volume() );
239 cells_dphi.push_back( cell_dde->
dphi() );
240 cells_deta.push_back( cell_dde->
deta() );
273 cells_pmt1_ros.push_back( ros1 );
274 cells_pmt2_ros.push_back( ros2 );
275 cells_pmt1_drawer.push_back( drawer1 );
276 cells_pmt2_drawer.push_back( drawer2 );
277 cells_pmt1_channel.push_back( channel1 );
278 cells_pmt2_channel.push_back( channel2 );
280 cells_pmt1_energy.push_back( tile_cell->
ene1() );
281 cells_pmt2_energy.push_back( tile_cell->
ene2() );
282 cells_pmt1_time.push_back( tile_cell->
time1() );
283 cells_pmt2_time.push_back( tile_cell->
time2() );
284 cells_pmt1_quality.push_back( tile_cell->
qual1() );
285 cells_pmt2_quality.push_back( tile_cell->
qual2() );
286 cells_pmt1_qbit.push_back( tile_cell->
qbit1() );
287 cells_pmt2_qbit.push_back( tile_cell->
qbit2() );
288 cells_pmt1_bad.push_back( tile_cell->
badch1() );
289 cells_pmt2_bad.push_back( tile_cell->
badch2() );
290 cells_pmt1_gain.push_back( tile_cell->
gain1() );
291 cells_pmt2_gain.push_back( tile_cell->
gain2() );
294 cellsEnergy(*
particle) = std::move(cells_energy);
295 cellsEt(*
particle) = std::move(cells_et);
296 cellsEta(*
particle) = std::move(cells_eta);
297 cellsPhi(*
particle) = std::move(cells_phi);
298 cellsGain(*
particle) = std::move(cells_gain);
299 cellsBad(*
particle) = std::move(cells_bad);
300 cellsSampling(*
particle) = std::move(cells_sampling);
301 cellsTime(*
particle) = std::move(cells_time);
302 cellsQuality(*
particle) = std::move(cells_quality);
304 cellsSinTh(*
particle) = std::move(cells_sinTh);
305 cellsCosTh(*
particle) = std::move(cells_cosTh);
306 cellsCotTh(*
particle) = std::move(cells_cotTh);
307 cellsX(*
particle) = std::move(cells_x);
308 cellsY(*
particle) = std::move(cells_y);
309 cellsZ(*
particle) = std::move(cells_z);
310 cellsR(*
particle) = std::move(cells_r);
311 cellsDx(*
particle) = std::move(cells_dx);
312 cellsDy(*
particle) = std::move(cells_dy);
313 cellsDz(*
particle) = std::move(cells_dz);
314 cellsDr(*
particle) = std::move(cells_dr);
315 cellsVolume(*
particle) = std::move(cells_volume);
316 cellsDeta(*
particle) = std::move(cells_deta);
317 cellsDphi(*
particle) = std::move(cells_dphi);
318 cellsSide(*
particle) = std::move(cells_side);
319 cellsSection(*
particle) = std::move(cells_section);
320 cellsModule(*
particle) = std::move(cells_module);
321 cellsTower(*
particle) = std::move(cells_tower);
322 cellsSample(*
particle) = std::move(cells_sample);
324 cellsPmt1Ros(*
particle) = std::move(cells_pmt1_ros);
325 cellsPmt2Ros(*
particle) = std::move(cells_pmt2_ros);
326 cellsPmt1Drawer(*
particle) = std::move(cells_pmt1_drawer);
327 cellsPmt2Drawer(*
particle) = std::move(cells_pmt2_drawer);
328 cellsPmt1Channel(*
particle) = std::move(cells_pmt1_channel);
329 cellsPmt2Channel(*
particle) = std::move(cells_pmt2_channel);
330 cellsPmt1Energy(*
particle) = std::move(cells_pmt1_energy);
331 cellsPmt2Energy(*
particle) = std::move(cells_pmt2_energy);
332 cellsPmt1Time(*
particle) = std::move(cells_pmt1_time);
333 cellsPmt2Time(*
particle) = std::move(cells_pmt2_time);
334 cellsPmt1Quality(*
particle) = std::move(cells_pmt1_quality);
335 cellsPmt2Quality(*
particle) = std::move(cells_pmt2_quality);
336 cellsPmt1_Qbit(*
particle) = std::move(cells_pmt1_qbit);
337 cellsPmt2_Qbit(*
particle) = std::move(cells_pmt2_qbit);
338 cellsPmt1Bad(*
particle) = std::move(cells_pmt1_bad);
339 cellsPmt2Bad(*
particle) = std::move(cells_pmt2_bad);
340 cellsPmt1Gain(*
particle) = std::move(cells_pmt1_gain);
341 cellsPmt2Gain(*
particle) = std::move(cells_pmt2_gain);
345 return StatusCode::SUCCESS;