Inheritance diagram for CaloParticleDebugger:

Public Member Functions
	CaloParticleDebugger (const edm::ParameterSet &)

	~CaloParticleDebugger () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	endJob () override

void	fillSimHits (std::map< int, float > &, const edm::Event &, const edm::EventSetup &)

Private Attributes
edm::InputTag	caloParticles_

edm::EDGetTokenT< std::vector< CaloParticle > >	caloParticlesToken_

std::vector< edm::InputTag >	collectionTags_

std::vector< edm::EDGetTokenT< std::vector< PCaloHit > > >	collectionTagsToken_

edm::InputTag	genParticles_

edm::EDGetTokenT< std::vector< reco::GenParticle > >	genParticlesToken_

int	geometryType_ = 0

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomToken_

edm::InputTag	simClusters_

edm::EDGetTokenT< std::vector< SimCluster > >	simClustersToken_

edm::InputTag	simTracks_

edm::EDGetTokenT< std::vector< SimTrack > >	simTracksToken_

edm::InputTag	simVertices_

edm::EDGetTokenT< std::vector< SimVertex > >	simVerticesToken_

edm::InputTag	trackingParticles_

edm::EDGetTokenT< std::vector< TrackingParticle > >	trackingParticlesToken_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 47 of file CaloParticleDebugger.cc.

Constructor & Destructor Documentation

◆ CaloParticleDebugger()

CaloParticleDebugger::CaloParticleDebugger ( const edm::ParameterSet & iConfig )

explicit

Definition at line 89 of file CaloParticleDebugger.cc.

References caloParticles_, caloParticlesToken_, collectionTags_, collectionTagsToken_, edm::ConsumesCollector::consumes(), edm::EDConsumerBase::consumesCollector(), edm::ConsumesCollector::esConsumes(), genParticles_, genParticlesToken_, geomToken_, simClusters_, simClustersToken_, simTracks_, simTracksToken_, simVertices_, simVerticesToken_, trackingParticles_, and trackingParticlesToken_.

     : simTracks_(iConfig.getParameter<edm::InputTag>("simTracks")),
       genParticles_(iConfig.getParameter<edm::InputTag>("genParticles")),
       simVertices_(iConfig.getParameter<edm::InputTag>("simVertices")),
       trackingParticles_(iConfig.getParameter<edm::InputTag>("trackingParticles")),
       caloParticles_(iConfig.getParameter<edm::InputTag>("caloParticles")),
       simClusters_(iConfig.getParameter<edm::InputTag>("simClusters")),
       collectionTags_(iConfig.getParameter<std::vector<edm::InputTag>>("collectionTags")) {
   edm::ConsumesCollector&& iC = consumesCollector();
   simTracksToken_ = iC.consumes<std::vector<SimTrack>>(simTracks_);
   genParticlesToken_ = iC.consumes<std::vector<reco::GenParticle>>(genParticles_);
   simVerticesToken_ = iC.consumes<std::vector<SimVertex>>(simVertices_);
   trackingParticlesToken_ = iC.consumes<std::vector<TrackingParticle>>(trackingParticles_);
   caloParticlesToken_ = iC.consumes<std::vector<CaloParticle>>(caloParticles_);
   simClustersToken_ = iC.consumes<std::vector<SimCluster>>(simClusters_);
   for (auto const& collectionTag : collectionTags_) {
     collectionTagsToken_.push_back(iC.consumes<std::vector<PCaloHit>>(collectionTag));
   }
   geomToken_ = iC.esConsumes<CaloGeometry, CaloGeometryRecord>();
 }

◆ ~CaloParticleDebugger()

CaloParticleDebugger::~CaloParticleDebugger ( )

override

Definition at line 110 of file CaloParticleDebugger.cc.

110 {}

Member Function Documentation

◆ analyze()

void CaloParticleDebugger::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 117 of file CaloParticleDebugger.cc.

References funct::abs(), caloParticlesToken_, haddnano::cl, mps_fire::end, fillSimHits(), AJJGenJetFilter_cfi::genParticles, runTauDisplay::genParticlesH, genParticlesToken_, runTauDisplay::gp, mps_fire::i, heavyIonCSV_trainingSettings::idx, iEvent, dqmiolumiharvest::j, edm::Handle< T >::product(), SimTracksters_cff::simclusters, simClustersToken_, simTracksToken_, simVerticesToken_, jetUpdater_cfi::sort, submitPVValidationJobs::t, cmsswSequenceInfo::tp, trackingParticlesToken_, pwdgSkimBPark_cfi::tracks, findQualityFiles::v, and AlignmentTracksFromVertexSelector_cfi::vertices.

                                                                                       {
   using namespace edm;
   using std::begin;
   using std::end;
   using std::iota;
   using std::sort;
 
   edm::Handle<std::vector<SimTrack>> simTracksH;
   edm::Handle<std::vector<reco::GenParticle>> genParticlesH;
   edm::Handle<std::vector<SimVertex>> simVerticesH;
   edm::Handle<std::vector<TrackingParticle>> trackingParticlesH;
   edm::Handle<std::vector<CaloParticle>> caloParticlesH;
   edm::Handle<std::vector<SimCluster>> simClustersH;
 
   iEvent.getByToken(simTracksToken_, simTracksH);
   auto const& tracks = *simTracksH.product();
   std::vector<int> sorted_tracks_idx(tracks.size());
   iota(begin(sorted_tracks_idx), end(sorted_tracks_idx), 0);
   sort(begin(sorted_tracks_idx), end(sorted_tracks_idx), [&tracks](int i, int j) {
     return tracks[i].momentum().eta() < tracks[j].momentum().eta();
   });
 
   iEvent.getByToken(genParticlesToken_, genParticlesH);
   auto const& genParticles = *genParticlesH.product();
   std::vector<int> sorted_genParticles_idx(genParticles.size());
   iota(begin(sorted_genParticles_idx), end(sorted_genParticles_idx), 0);
   sort(begin(sorted_genParticles_idx), end(sorted_genParticles_idx), [&genParticles](int i, int j) {
     return genParticles[i].momentum().eta() < genParticles[j].momentum().eta();
   });
 
   iEvent.getByToken(simVerticesToken_, simVerticesH);
   auto const& vertices = *simVerticesH.product();
   std::vector<int> sorted_vertices_idx(vertices.size());
   iota(begin(sorted_vertices_idx), end(sorted_vertices_idx), 0);
   sort(begin(sorted_vertices_idx), end(sorted_vertices_idx), [&vertices](int i, int j) {
     return vertices[i].vertexId() < vertices[j].vertexId();
   });
 
   iEvent.getByToken(trackingParticlesToken_, trackingParticlesH);
   auto const& trackingpart = *trackingParticlesH.product();
   std::vector<int> sorted_tp_idx(trackingpart.size());
   iota(begin(sorted_tp_idx), end(sorted_tp_idx), 0);
   sort(begin(sorted_tp_idx), end(sorted_tp_idx), [&trackingpart](int i, int j) {
     return trackingpart[i].eta() < trackingpart[j].eta();
   });
 
   iEvent.getByToken(caloParticlesToken_, caloParticlesH);
   auto const& calopart = *caloParticlesH.product();
   std::vector<int> sorted_cp_idx(calopart.size());
   iota(begin(sorted_cp_idx), end(sorted_cp_idx), 0);
   sort(begin(sorted_cp_idx), end(sorted_cp_idx), [&calopart](int i, int j) {
     return calopart[i].eta() < calopart[j].eta();
   });
 
   iEvent.getByToken(simClustersToken_, simClustersH);
   auto const& simclusters = *simClustersH.product();
   std::vector<int> sorted_simcl_idx(simclusters.size());
   iota(begin(sorted_simcl_idx), end(sorted_simcl_idx), 0);
   sort(begin(sorted_simcl_idx), end(sorted_simcl_idx), [&simclusters](int i, int j) {
     return simclusters[i].eta() < simclusters[j].eta();
   });
 
   // Let's first fill in hits information
   std::map<int, float> detIdToTotalSimEnergy;
   fillSimHits(detIdToTotalSimEnergy, iEvent, iSetup);
 
   int idx = 0;
 
   std::map<int, int> trackid_to_track_index;
   LogVerbatim("CaloParticleDebuggerSimTracks") << "\n\n**Printing SimTracks information **";
   LogVerbatim("CaloParticleDebuggerSimTracks") << "IDX\tTrackId\tPDGID\tMOMENTUM(x,y,z,E)\tVertexIdx\tGenPartIdx";
   for (auto i : sorted_tracks_idx) {
     auto const& t = tracks[i];
     LogVerbatim("CaloParticleDebuggerSimTracks") << idx << "\t" << t.trackId() << "\t" << t;
     LogVerbatim("CaloParticleDebuggerSimTracks")
         << "Crossed  Boundary: " << t.crossedBoundary() << "  Position Boundary: " << t.getPositionAtBoundary()
         << "  Momentum Boundary: " << t.getMomentumAtBoundary() << "  Vtx:               " << t.vertIndex()
         << "  Momemtum Origin:   " << t.momentum();
     trackid_to_track_index[t.trackId()] = idx;
     idx++;
   }
 
   LogVerbatim("CaloParticleDebuggerGenParticles") << "\n\n**Printing GenParticles information **";
   LogVerbatim("CaloParticleDebuggerGenParticles") << "IDX\tPDGID\tMOMENTUM(x,y,z)\tVertex(x,y,z)";
   for (auto i : sorted_genParticles_idx) {
     auto const& gp = genParticles[i];
     LogVerbatim("CaloParticleDebuggerGenParticles")
         << i << "\t" << gp.pdgId() << "\t" << gp.momentum() << "\t" << gp.vertex();
   }
 
   LogVerbatim("CaloParticleDebuggerSimVertices") << "\n\n**Printing SimVertex information **";
   LogVerbatim("CaloParticleDebuggerSimVertices") << "IDX\tPOSITION(x,y,z)\tPARENT_INDEX\tVERTEX_ID";
   for (auto i : sorted_vertices_idx) {
     auto const& v = vertices[i];
     LogVerbatim("CaloParticleDebuggerSimVertices") << i << "\t" << v;
   }
 
   LogVerbatim("CaloParticleDebuggerTrackingParticles") << "\n\n**Printing TrackingParticles information **";
   for (auto i : sorted_tp_idx) {
     auto const& tp = trackingpart[i];
     LogVerbatim("CaloParticleDebuggerTrackingParticles") << i << "\t" << tp;
   }
 
   LogVerbatim("CaloParticleDebuggerCaloParticles") << "\n\n**Printing CaloParticles information **";
   idx = 0;
   for (auto i : sorted_cp_idx) {
     auto const& cp = calopart[i];
     LogVerbatim("CaloParticleDebuggerCaloParticles")
         << "\n\n"
         << idx++ << " |Eta|: " << std::abs(cp.momentum().eta()) << "\tType: " << cp.pdgId()
         << "\tEnergy: " << cp.energy() << "\tIdx: " << cp.g4Tracks()[0].trackId();  // << cp ;
     double total_sim_energy = 0.;
     double total_cp_energy = 0.;
     LogVerbatim("CaloParticleDebuggerCaloParticles") << "--> Overall simclusters in CP: " << cp.simClusters().size();
     // All the next mess just to print the simClusters ordered
     auto const& simcs = cp.simClusters();
     std::vector<int> sorted_sc_idx(simcs.size());
     iota(begin(sorted_sc_idx), end(sorted_sc_idx), 0);
     sort(begin(sorted_sc_idx), end(sorted_sc_idx), [&simcs](int i, int j) {
       return simcs[i]->momentum().eta() < simcs[j]->momentum().eta();
     });
     for (auto i : sorted_sc_idx) {
       LogVerbatim("CaloParticleDebuggerCaloParticles") << *(simcs[i]);
     }
 
     for (auto const& sc : cp.simClusters()) {
       for (auto const& cl : sc->hits_and_fractions()) {
         total_sim_energy += detIdToTotalSimEnergy[cl.first] * cl.second;
         total_cp_energy += cp.energy() * cl.second;
       }
     }
     LogVerbatim("CaloParticleDebuggerCaloParticles")
         << "--> Overall SC energy (sum using sim energies): " << total_sim_energy;
     LogVerbatim("CaloParticleDebuggerCaloParticles")
         << "--> Overall SC energy (sum using CaloP energies): " << total_cp_energy;
   }
 
   idx = 0;
   LogVerbatim("CaloParticleDebuggerSimClusters") << "\n\n**Printing SimClusters information **";
   for (auto i : sorted_simcl_idx) {
     auto const& simcl = simclusters[i];
     LogVerbatim("CaloParticleDebuggerSimClusters")
         << "\n\n"
         << idx++ << " |Eta|: " << std::abs(simcl.momentum().eta()) << "\tType: " << simcl.pdgId()
         << "\tEnergy: " << simcl.energy() << "\tKey: " << i;  // << simcl ;
     auto const& simtrack = simcl.g4Tracks()[0];
     LogVerbatim("CaloParticleDebuggerSimClusters") << "  Crossed  Boundary: " << simtrack.crossedBoundary()
                                                    << "  Position Boundary: " << simtrack.getPositionAtBoundary()
                                                    << "  Momentum Boundary: " << simtrack.getMomentumAtBoundary();
     double total_sim_energy = 0.;
     LogVerbatim("CaloParticleDebuggerSimClusters") << "--> Overall simclusters's size: " << simcl.numberOfRecHits();
     for (auto const& cl : simcl.hits_and_fractions()) {
       total_sim_energy += detIdToTotalSimEnergy[cl.first] * cl.second;
     }
     LogVerbatim("CaloParticleDebuggerSimClusters") << simcl;
     LogVerbatim("CaloParticleDebuggerSimClusters")
         << "--> Overall SimCluster energy (sum using sim energies): " << total_sim_energy;
   }
 }

◆ beginJob()

void CaloParticleDebugger::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 278 of file CaloParticleDebugger.cc.

278 {}

◆ endJob()

void CaloParticleDebugger::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 281 of file CaloParticleDebugger.cc.

281 {}

◆ fillDescriptions()

void CaloParticleDebugger::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 360 of file CaloParticleDebugger.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_2022v15_cff::InputTag.

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("simTracks", edm::InputTag("g4SimHits"));
   desc.add<edm::InputTag>("genParticles", edm::InputTag("genParticles"));
   desc.add<edm::InputTag>("simVertices", edm::InputTag("g4SimHits"));
   desc.add<edm::InputTag>("trackingParticles", edm::InputTag("mix", "MergedTrackTruth"));
   desc.add<edm::InputTag>("caloParticles", edm::InputTag("mix", "MergedCaloTruth"));
   desc.add<edm::InputTag>("simClusters", edm::InputTag("mix", "MergedCaloTruth"));
   desc.add<std::vector<edm::InputTag>>("collectionTags",
                                        {edm::InputTag("g4SimHits", "HGCHitsEE"),
                                         edm::InputTag("g4SimHits", "HGCHitsHEfront"),
                                         edm::InputTag("g4SimHits", "HcalHits")});
   descriptions.add("caloParticleDebugger", desc);
 }

◆ fillSimHits()

void CaloParticleDebugger::fillSimHits	(	std::map< int, float > &	detIdToTotalSimEnergy,
		const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 283 of file CaloParticleDebugger.cc.

References cms::cuda::assert(), collectionTags_, collectionTagsToken_, HGCalTopology::dddConstants(), HcalTopology::dddConstants(), EcalMappingRecord_cfi::eegeom, DetId::Forward, ForwardEmpty, relativeConstraints::geom, geometryType_, geomToken_, edm::EventSetup::getHandle(), DetId::Hcal, HcalEndcap, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCEE, HGCHEF, mps_fire::i, l1ctLayer2EG_cff::id, iEvent, pixelTopology::layer, HcalHitRelabeller::relabel(), fileinputsource_cfi::sec, rpcPointValidation_cfi::simHit, HGCalDDDConstants::simToReco(), HcalDetId::subdet(), unpackBuffers-CaloStage2::token, HcalGeometry::topology(), and HGCalTestNumbering::unpackHexagonIndex().

Referenced by analyze().

                                                                     {
   // Taken needed quantities from the EventSetup
   const auto& geom = iSetup.getHandle(geomToken_);
   const HGCalGeometry *eegeom = nullptr, *fhgeom = nullptr, *bhgeomnew = nullptr;
   const HcalGeometry* bhgeom = nullptr;
   const HGCalDDDConstants* hgddd[3];
   const HGCalTopology* hgtopo[3];
   const HcalDDDRecConstants* hcddd = nullptr;
   eegeom =
       static_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::HGCalEE, ForwardSubdetector::ForwardEmpty));
   //check if it's the new geometry
   if (eegeom) {
     geometryType_ = 1;
     fhgeom = static_cast<const HGCalGeometry*>(
         geom->getSubdetectorGeometry(DetId::HGCalHSi, ForwardSubdetector::ForwardEmpty));
     bhgeomnew = static_cast<const HGCalGeometry*>(
         geom->getSubdetectorGeometry(DetId::HGCalHSc, ForwardSubdetector::ForwardEmpty));
   } else {
     geometryType_ = 0;
     eegeom = static_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::Forward, HGCEE));
     fhgeom = static_cast<const HGCalGeometry*>(geom->getSubdetectorGeometry(DetId::Forward, HGCHEF));
     bhgeom = static_cast<const HcalGeometry*>(geom->getSubdetectorGeometry(DetId::Hcal, HcalEndcap));
   }
   hgtopo[0] = &(eegeom->topology());
   hgtopo[1] = &(fhgeom->topology());
   if (bhgeomnew)
     hgtopo[2] = &(bhgeomnew->topology());
 
   for (unsigned i = 0; i < 3; ++i) {
     if (hgtopo[i])
       hgddd[i] = &(hgtopo[i]->dddConstants());
   }
 
   if (bhgeom)
     hcddd = bhgeom->topology().dddConstants();
 
   // loop over the collections
   int token = 0;
   for (auto const& collectionTag : collectionTags_) {
     edm::Handle<std::vector<PCaloHit>> hSimHits;
     const bool isHcal = (collectionTag.instance().find("HcalHits") != std::string::npos);
     iEvent.getByToken(collectionTagsToken_[token++], hSimHits);
     for (auto const& simHit : *hSimHits) {
       DetId id(0);
       const uint32_t simId = simHit.id();
       if (geometryType_ == 1) {
         //no test numbering in new geometry
         id = simId;
       } else if (isHcal) {
         assert(hcddd);
         HcalDetId hid = HcalHitRelabeller::relabel(simId, hcddd);
         if (hid.subdet() == HcalEndcap)
           id = hid;
       } else {
         int subdet, layer, cell, sec, subsec, zp;
         HGCalTestNumbering::unpackHexagonIndex(simId, subdet, zp, layer, sec, subsec, cell);
         const HGCalDDDConstants* ddd = hgddd[subdet - 3];
         std::pair<int, int> recoLayerCell = ddd->simToReco(cell, layer, sec, hgtopo[subdet - 3]->detectorType());
         cell = recoLayerCell.first;
         layer = recoLayerCell.second;
         // skip simhits with bad barcodes or non-existant layers
         if (layer == -1 || simHit.geantTrackId() == 0)
           continue;
         id = HGCalDetId((ForwardSubdetector)subdet, zp, layer, subsec, sec, cell);
       }
 
       if (DetId(0) == id)
         continue;
 
       detIdToTotalSimEnergy[id.rawId()] += simHit.energy();
     }
   }  // end of loop over InputTags
 }

Member Data Documentation

◆ caloParticles_

edm::InputTag CaloParticleDebugger::caloParticles_

private

Definition at line 63 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger().

◆ caloParticlesToken_

edm::EDGetTokenT<std::vector<CaloParticle> > CaloParticleDebugger::caloParticlesToken_

private

Definition at line 70 of file CaloParticleDebugger.cc.

Referenced by analyze(), and CaloParticleDebugger().

◆ collectionTags_

std::vector<edm::InputTag> CaloParticleDebugger::collectionTags_

private

Definition at line 65 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger(), and fillSimHits().

◆ collectionTagsToken_

std::vector<edm::EDGetTokenT<std::vector<PCaloHit> > > CaloParticleDebugger::collectionTagsToken_

private

Definition at line 72 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger(), and fillSimHits().

◆ genParticles_

edm::InputTag CaloParticleDebugger::genParticles_

private

Definition at line 60 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger().

◆ genParticlesToken_

edm::EDGetTokenT<std::vector<reco::GenParticle> > CaloParticleDebugger::genParticlesToken_

private

Definition at line 67 of file CaloParticleDebugger.cc.

Referenced by analyze(), and CaloParticleDebugger().

◆ geometryType_

int CaloParticleDebugger::geometryType_ = 0

private

Definition at line 74 of file CaloParticleDebugger.cc.

Referenced by fillSimHits().

◆ geomToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> CaloParticleDebugger::geomToken_

private

Definition at line 73 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger(), and fillSimHits().

◆ simClusters_

edm::InputTag CaloParticleDebugger::simClusters_

private

Definition at line 64 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger().

◆ simClustersToken_

edm::EDGetTokenT<std::vector<SimCluster> > CaloParticleDebugger::simClustersToken_

private

Definition at line 71 of file CaloParticleDebugger.cc.

Referenced by analyze(), and CaloParticleDebugger().

◆ simTracks_

edm::InputTag CaloParticleDebugger::simTracks_

private

Definition at line 59 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger().

◆ simTracksToken_

edm::EDGetTokenT<std::vector<SimTrack> > CaloParticleDebugger::simTracksToken_

private

Definition at line 66 of file CaloParticleDebugger.cc.

Referenced by analyze(), and CaloParticleDebugger().

◆ simVertices_

edm::InputTag CaloParticleDebugger::simVertices_

private

Definition at line 61 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger().

◆ simVerticesToken_

edm::EDGetTokenT<std::vector<SimVertex> > CaloParticleDebugger::simVerticesToken_

private

Definition at line 68 of file CaloParticleDebugger.cc.

Referenced by analyze(), and CaloParticleDebugger().

◆ trackingParticles_

edm::InputTag CaloParticleDebugger::trackingParticles_

private

Definition at line 62 of file CaloParticleDebugger.cc.

Referenced by CaloParticleDebugger().

◆ trackingParticlesToken_

edm::EDGetTokenT<std::vector<TrackingParticle> > CaloParticleDebugger::trackingParticlesToken_

private

Definition at line 69 of file CaloParticleDebugger.cc.

Referenced by analyze(), and CaloParticleDebugger().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CaloParticleDebugger()

◆ ~CaloParticleDebugger()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

◆ fillSimHits()

Member Data Documentation

◆ caloParticles_

◆ caloParticlesToken_

◆ collectionTags_

◆ collectionTagsToken_

◆ genParticles_

◆ genParticlesToken_

◆ geometryType_

◆ geomToken_

◆ simClusters_

◆ simClustersToken_

◆ simTracks_

◆ simTracksToken_

◆ simVertices_

◆ simVerticesToken_

◆ trackingParticles_

◆ trackingParticlesToken_