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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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

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

std::vector< ESResolverIndex > 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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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	genBarcodes_

edm::EDGetTokenT< std::vector< int > >	genBarcodesToken_

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 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 91 of file CaloParticleDebugger.cc.

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

     : simTracks_(iConfig.getParameter<edm::InputTag>("simTracks")),
       genParticles_(iConfig.getParameter<edm::InputTag>("genParticles")),
       genBarcodes_(iConfig.getParameter<edm::InputTag>("genBarcodes")),
       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_);
   genBarcodesToken_ = iC.consumes<std::vector<int>>(genBarcodes_);
   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 114 of file CaloParticleDebugger.cc.

114 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 121 of file CaloParticleDebugger.cc.

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

                                                                                       {
   using namespace edm;
   using std::begin;
   using std::end;
 
   edm::Handle<std::vector<SimTrack>> simTracksH;
   edm::Handle<std::vector<reco::GenParticle>> genParticlesH;
   edm::Handle<std::vector<int>> genBarcodesH;
   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();
 
   iEvent.getByToken(genParticlesToken_, genParticlesH);
   auto const& genParticles = *genParticlesH.product();
 
   iEvent.getByToken(genBarcodesToken_, genBarcodesH);
   auto const& genBarcodes = *genBarcodesH.product();
 
   iEvent.getByToken(simVerticesToken_, simVerticesH);
   auto const& vertices = *simVerticesH.product();
 
   iEvent.getByToken(trackingParticlesToken_, trackingParticlesH);
   auto const& trackingpart = *trackingParticlesH.product();
 
   iEvent.getByToken(caloParticlesToken_, caloParticlesH);
   auto const& calopart = *caloParticlesH.product();
 
   iEvent.getByToken(simClustersToken_, simClustersH);
   auto const& simclusters = *simClustersH.product();
 
   // Let's first fill in hits information
   std::map<int, float> detIdToTotalSimEnergy;
   fillSimHits(detIdToTotalSimEnergy, iEvent, iSetup);
 
   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 (size_t i = 0; i < tracks.size(); ++i) {
     auto const& t = tracks[i];
     LogVerbatim("CaloParticleDebuggerSimTracks")
         << i << "\t" << t.trackId() << "\t" << t << "  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()] = i;
   }
 
   LogVerbatim("CaloParticleDebuggerGenParticles") << "\n\n**Printing GenParticles information **";
   LogVerbatim("CaloParticleDebuggerGenParticles") << "BARCODE\tPDGID\tMOMENTUM(x,y,z)\tVertex(x,y,z)";
   for (size_t i = 0; i < genParticles.size(); ++i) {
     auto const& gp = genParticles[i];
     LogVerbatim("CaloParticleDebuggerGenParticles")
         << genBarcodes[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 (size_t i = 0; i < vertices.size(); ++i) {
     auto const& v = vertices[i];
     LogVerbatim("CaloParticleDebuggerSimVertices") << i << "\t" << v;
   }
 
   LogVerbatim("CaloParticleDebuggerTrackingParticles") << "\n\n**Printing TrackingParticles information **";
   for (size_t i = 0; i < trackingpart.size(); ++i) {
     auto const& tp = trackingpart[i];
     LogVerbatim("CaloParticleDebuggerTrackingParticles") << i << "\t" << tp;
   }
 
   LogVerbatim("CaloParticleDebuggerCaloParticles") << "\n\n**Printing CaloParticles information **";
   for (size_t i = 0; i < calopart.size(); ++i) {
     auto const& cp = calopart[i];
     LogVerbatim("CaloParticleDebuggerCaloParticles") << "\n\n"
                                                      << i << "\tType: " << cp.pdgId() << "\tEnergy: " << cp.energy()
                                                      << "\tBarcode: " << cp.g4Tracks()[0].genpartIndex()
                                                      << " G4_trackID: " << 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();
     for (size_t j = 0; j < simcs.size(); ++j) {
       LogVerbatim("CaloParticleDebuggerCaloParticles") << *(simcs[j]);
     }
 
     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;
   }
 
   LogVerbatim("CaloParticleDebuggerSimClusters") << "\n\n**Printing SimClusters information **";
   for (size_t i = 0; i < simclusters.size(); ++i) {
     auto const& simcl = simclusters[i];
     LogVerbatim("CaloParticleDebuggerSimClusters")
         << "\n\n"
         << i << "\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 242 of file CaloParticleDebugger.cc.

242 {}

◆ endJob()

void CaloParticleDebugger::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 245 of file CaloParticleDebugger.cc.

245 {}

◆ fillDescriptions()

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

static

Definition at line 326 of file CaloParticleDebugger.cc.

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

                                                                                       {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("simTracks", edm::InputTag("g4SimHits"));
   desc.add<edm::InputTag>("genParticles", edm::InputTag("genParticles"));
   desc.add<edm::InputTag>("genBarcodes", 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 247 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, nano_mu_digi_cff::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());
   else
     hgtopo[2] = nullptr;
 
   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
 }