#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 49 of file InterestingTrackEcalDetIdProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 85 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 101 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 157 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 116 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 l1t::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));
 
 }