#include <Producers/InterestingTrackEcalDetIdProducer/src/InterestingTrackEcalDetIdProducer.cc>

Inheritance diagram for InterestingTrackEcalDetIdProducer:

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

	~InterestingTrackEcalDetIdProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
void	beginRun (edm::Run const &, const edm::EventSetup &) override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const CaloTopology *	caloTopology_

double	minTrackPt_

TrackDetectorAssociator	trackAssociator_

TrackAssociatorParameters	trackAssociatorParameters_

edm::ParameterSet	trackAssociatorPS_

edm::InputTag	trackCollection_

edm::EDGetTokenT< reco::TrackCollection >	trackCollectionToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 47 of file InterestingTrackEcalDetIdProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 79 of file InterestingTrackEcalDetIdProducer.cc.

References edm::ParameterSet::getParameter(), TrackAssociatorParameters::loadParameters(), trackAssociator_, trackAssociatorParameters_, trackAssociatorPS_, trackCollectionToken_, and TrackDetectorAssociator::useDefaultPropagator().

     :
 
       trackAssociatorPS_(iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters")),
       minTrackPt_(iConfig.getParameter<double>("MinTrackPt"))
 
 {
   trackCollectionToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("TrackCollection"));
   trackAssociator_.useDefaultPropagator();
   edm::ConsumesCollector iC = consumesCollector();
   trackAssociatorParameters_.loadParameters(trackAssociatorPS_, iC);
 
   produces<DetIdCollection>();
 }

InterestingTrackEcalDetIdProducer::~InterestingTrackEcalDetIdProducer ( )

override

Definition at line 94 of file InterestingTrackEcalDetIdProducer.cc.

                                                                       {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void InterestingTrackEcalDetIdProducer::beginRun	(	edm::Run const &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 139 of file InterestingTrackEcalDetIdProducer.cc.

References caloTopology_, DEFINE_FWK_MODULE, and edm::EventSetup::get().

                                                                                                {
   edm::ESHandle<CaloTopology> theCaloTopology;
   iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
   caloTopology_ = &(*theCaloTopology);
 }

void InterestingTrackEcalDetIdProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 104 of file InterestingTrackEcalDetIdProducer.cc.

References TrackDetectorAssociator::associate(), caloTopology_, DetId::Ecal, TrackDetMatchInfo::EcalRecHits, spr::find(), TrackDetMatchInfo::findMaxDeposition(), edm::Event::getByToken(), TrackDetectorAssociator::getFreeTrajectoryState(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), photons_cff::ids, info(), interestingDetIdCollectionProducer_cfi::interestingDetIdCollection, minTrackPt_, eostools::move(), edm::Event::put(), DetId::rawId(), DetId::subdetId(), ecaldqm::topology(), trackAssociator_, trackAssociatorParameters_, trackCollectionToken_, and PDWG_EXOHSCP_cff::tracks.

                                                                                              {
   using namespace edm;
 
   auto interestingDetIdCollection = std::make_unique<DetIdCollection>();
 
   // Get tracks from event
   edm::Handle<reco::TrackCollection> tracks;
   iEvent.getByToken(trackCollectionToken_, tracks);
 
   // Loop over tracks
   for (reco::TrackCollection::const_iterator tkItr = tracks->begin(); tkItr != tracks->end(); ++tkItr) {
     if (tkItr->pt() < minTrackPt_)
       continue;
 
     TrackDetMatchInfo info = trackAssociator_.associate(
         iEvent, iSetup, trackAssociator_.getFreeTrajectoryState(iSetup, *tkItr), trackAssociatorParameters_);
 
     DetId centerId = info.findMaxDeposition(TrackDetMatchInfo::EcalRecHits);
 
     if (centerId.rawId() == 0)
       continue;
 
     // Find 5x5 around max
     const CaloSubdetectorTopology* topology = caloTopology_->getSubdetectorTopology(DetId::Ecal, centerId.subdetId());
     const std::vector<DetId>& ids = topology->getWindow(centerId, 5, 5);
     for (std::vector<DetId>::const_iterator idItr = ids.begin(); idItr != ids.end(); ++idItr) {
       if (std::find(interestingDetIdCollection->begin(), interestingDetIdCollection->end(), *idItr) ==
           interestingDetIdCollection->end())
         interestingDetIdCollection->push_back(*idItr);
     }
   }
 
   iEvent.put(std::move(interestingDetIdCollection));
 }