279 if (trigger_cells_h.
isValid()) {
280 for (
auto tc_itr = trigger_cells.
begin(0); tc_itr != trigger_cells.
end(0); tc_itr++) {
284 histograms.h_tc_mipPt_->Fill(tc_itr->mipPt());
285 histograms.h_tc_energy_->Fill(tc_itr->energy());
288 histograms.h_tc_x_->Fill(tc_itr->position().x());
289 histograms.h_tc_y_->Fill(tc_itr->position().y());
290 histograms.h_tc_z_->Fill(tc_itr->position().z());
302 for (
auto cl_itr =
clusters.begin(0); cl_itr !=
clusters.end(0); cl_itr++) {
304 histograms.h_cl_mipPt_->Fill(cl_itr->mipPt());
306 histograms.h_cl_energy_->Fill(cl_itr->energy());
310 histograms.h_cl_cells_n_->Fill(cl_itr->constituents().size());
320 if (multiclusters_h.
isValid()) {
321 for (
auto cl3d_itr = multiclusters.
begin(0); cl3d_itr != multiclusters.
end(0); cl3d_itr++) {
324 histograms.h_cl3d_energy_->Fill(cl3d_itr->energy());
325 histograms.h_cl3d_eta_->Fill(cl3d_itr->eta());
326 histograms.h_cl3d_phi_->Fill(cl3d_itr->phi());
327 histograms.h_cl3d_clusters_n_->Fill(cl3d_itr->constituents().size());
329 histograms.h_cl3d_showerlength_->Fill(cl3d_itr->showerLength());
330 histograms.h_cl3d_coreshowerlength_->Fill(cl3d_itr->coreShowerLength());
331 histograms.h_cl3d_firstlayer_->Fill(cl3d_itr->firstLayer());
332 histograms.h_cl3d_maxlayer_->Fill(cl3d_itr->maxLayer());
333 histograms.h_cl3d_seetot_->Fill(cl3d_itr->sigmaEtaEtaTot());
334 histograms.h_cl3d_seemax_->Fill(cl3d_itr->sigmaEtaEtaMax());
335 histograms.h_cl3d_spptot_->Fill(cl3d_itr->sigmaPhiPhiTot());
336 histograms.h_cl3d_sppmax_->Fill(cl3d_itr->sigmaPhiPhiMax());
337 histograms.h_cl3d_szz_->Fill(cl3d_itr->sigmaZZ());
338 histograms.h_cl3d_srrtot_->Fill(cl3d_itr->sigmaRRTot());
339 histograms.h_cl3d_srrmax_->Fill(cl3d_itr->sigmaRRMax());
340 histograms.h_cl3d_srrmean_->Fill(cl3d_itr->sigmaRRMean());
341 histograms.h_cl3d_emaxe_->Fill(cl3d_itr->eMax() / cl3d_itr->energy());
342 histograms.h_cl3d_hoe_->Fill(cl3d_itr->hOverE());
343 histograms.h_cl3d_meanz_->Fill(cl3d_itr->zBarycenter());
344 histograms.h_cl3d_layer10_->Fill(cl3d_itr->layer10percent());
345 histograms.h_cl3d_layer50_->Fill(cl3d_itr->layer50percent());
346 histograms.h_cl3d_layer90_->Fill(cl3d_itr->layer90percent());
347 histograms.h_cl3d_first1layers_->Fill(cl3d_itr->first1layers());
348 histograms.h_cl3d_first3layers_->Fill(cl3d_itr->first3layers());
349 histograms.h_cl3d_first5layers_->Fill(cl3d_itr->first5layers());
350 histograms.h_cl3d_firstHcal1layers_->Fill(cl3d_itr->firstHcal1layers());
351 histograms.h_cl3d_firstHcal3layers_->Fill(cl3d_itr->firstHcal3layers());
352 histograms.h_cl3d_firstHcal5layers_->Fill(cl3d_itr->firstHcal5layers());
353 histograms.h_cl3d_last1layers_->Fill(cl3d_itr->last1layers());
354 histograms.h_cl3d_last3layers_->Fill(cl3d_itr->last3layers());
355 histograms.h_cl3d_last5layers_->Fill(cl3d_itr->last5layers());
356 histograms.h_cl3d_emax1layers_->Fill(cl3d_itr->emax1layers());
357 histograms.h_cl3d_emax3layers_->Fill(cl3d_itr->emax3layers());
358 histograms.h_cl3d_emax5layers_->Fill(cl3d_itr->emax5layers());
359 histograms.h_cl3d_eot_->Fill(cl3d_itr->eot());
360 histograms.h_cl3d_ebm0_->Fill(cl3d_itr->ebm0());
361 histograms.h_cl3d_ebm1_->Fill(cl3d_itr->ebm1());
362 histograms.h_cl3d_hbm_->Fill(cl3d_itr->hbm());
363 histograms.h_cl3d_ntc67_->Fill(cl3d_itr->triggerCells67percent());
364 histograms.h_cl3d_ntc90_->Fill(cl3d_itr->triggerCells90percent());
366 histograms.h_cl3d_quality_->Fill(cl3d_itr->hwQual());
377 for (
auto tower_itr =
towers.begin(0); tower_itr !=
towers.end(0); tower_itr++) {
379 histograms.h_tower_pt_->Fill(tower_itr->pt());
380 histograms.h_tower_energy_->Fill(tower_itr->energy());
381 histograms.h_tower_eta_->Fill(tower_itr->eta());
382 histograms.h_tower_phi_->Fill(tower_itr->phi());
383 histograms.h_tower_etEm_->Fill(tower_itr->etEm());
384 histograms.h_tower_etHad_->Fill(tower_itr->etHad());
385 histograms.h_tower_iEta_->Fill(tower_itr->id().iEta());
386 histograms.h_tower_iPhi_->Fill(tower_itr->id().iPhi());
const edm::EDGetToken clusters_token_
const_iterator begin(int bx) const
std::unique_ptr< HGCalTriggerClusterIdentificationBase > id_
const edm::EDGetToken towers_token_
const edm::EDGetToken multiclusters_token_
const_iterator end(int bx) const
std::shared_ptr< HGCalTriggerTools > triggerTools_
const edm::EDGetToken trigger_cells_token_