#include <GsfElectronCoreBaseProducer.h>

Inheritance diagram for GsfElectronCoreBaseProducer:

Public Member Functions
	GsfElectronCoreBaseProducer (const edm::ParameterSet &conf)

	~GsfElectronCoreBaseProducer () override

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

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::vector< ModuleDescription const * > &modules, 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescription (edm::ParameterSetDescription &)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	fillElectronCore (reco::GsfElectronCore *)

void	initEvent (edm::Event &event, const edm::EventSetup &setup)

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

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 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)

Protected Attributes
edm::Handle< reco::TrackCollection >	ctfTracksH_

edm::Handle< reco::GsfPFRecTrackCollection >	gsfPfRecTracksH_

edm::Handle< reco::GsfTrackCollection >	gsfTracksH_

bool	useGsfPfRecTracks_

Private Member Functions
std::pair< reco::TrackRef, float >	getCtfTrackRef (const reco::GsfTrackRef &)

Private Attributes
edm::EDGetTokenT< reco::TrackCollection >	ctfTracksTag_

edm::EDGetTokenT< reco::GsfPFRecTrackCollection >	gsfPfRecTracksTag_

edm::EDGetTokenT< reco::GsfTrackCollection >	gsfTracksTag_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 31 of file GsfElectronCoreBaseProducer.h.

Constructor & Destructor Documentation

GsfElectronCoreBaseProducer::GsfElectronCoreBaseProducer ( const edm::ParameterSet & conf )

explicit

Definition at line 28 of file GsfElectronCoreBaseProducer.cc.

References edm::ParameterSet::getParameter().

  {
   produces<GsfElectronCoreCollection>() ;
   gsfPfRecTracksTag_ = mayConsume<reco::GsfPFRecTrackCollection>(config.getParameter<edm::InputTag>("gsfPfRecTracks")) ;
   gsfTracksTag_ = consumes<reco::GsfTrackCollection>(config.getParameter<edm::InputTag>("gsfTracks"));
   ctfTracksTag_ = consumes<reco::TrackCollection>(config.getParameter<edm::InputTag>("ctfTracks"));
   useGsfPfRecTracks_ = config.getParameter<bool>("useGsfPfRecTracks") ;
  }

GsfElectronCoreBaseProducer::~GsfElectronCoreBaseProducer ( )

override

Definition at line 37 of file GsfElectronCoreBaseProducer.cc.

38 {}

Member Function Documentation

void GsfElectronCoreBaseProducer::fillDescription ( edm::ParameterSetDescription & desc )

static

Definition at line 20 of file GsfElectronCoreBaseProducer.cc.

References edm::ParameterSetDescription::add().

Referenced by GEDGsfElectronCoreProducer::fillDescriptions().

  {
   desc.add<edm::InputTag>("gsfPfRecTracks",edm::InputTag("pfTrackElec")) ;
   desc.add<edm::InputTag>("gsfTracks",edm::InputTag("electronGsfTracks")) ;
   desc.add<edm::InputTag>("ctfTracks",edm::InputTag("generalTracks")) ;
   desc.add<bool>("useGsfPfRecTracks",true) ;
  }

void GsfElectronCoreBaseProducer::fillElectronCore ( reco::GsfElectronCore * eleCore )

protected

Definition at line 54 of file GsfElectronCoreBaseProducer.cc.

References getCtfTrackRef(), reco::GsfElectronCore::gsfTrack(), and reco::GsfElectronCore::setCtfTrack().

Referenced by GsfElectronCoreEcalDrivenProducer::produceEcalDrivenCore(), GEDGsfElectronCoreProducer::produceElectronCore(), and GsfElectronCoreProducer::produceTrackerDrivenCore().

  {
   const GsfTrackRef & gsfTrackRef = eleCore->gsfTrack() ;
 
   std::pair<TrackRef,float> ctfpair = getCtfTrackRef(gsfTrackRef) ;
   eleCore->setCtfTrack(ctfpair.first,ctfpair.second) ;
  }

std::pair< TrackRef, float > GsfElectronCoreBaseProducer::getCtfTrackRef ( const reco::GsfTrackRef & gsfTrackRef )

private

Definition at line 68 of file GsfElectronCoreBaseProducer.cc.

References funct::abs(), TrackingRecHit::all, counter, reco::HitPattern::getHitPattern(), reco::HitPattern::hitPattern, min(), pi, reco::HitPattern::pixelHitFilter(), mathSSE::sqrt(), reco::HitPattern::stripTIBHitFilter(), reco::HitPattern::stripTIDHitFilter(), and reco::HitPattern::TRACK_HITS.

Referenced by fillElectronCore().

  {
   float maxFracShared = 0;
   TrackRef ctfTrackRef = TrackRef() ;
   const TrackCollection * ctfTrackCollection = ctfTracksH_.product() ;
 
   // get the Hit Pattern for the gsfTrack
   const HitPattern& gsfHitPattern = gsfTrackRef->hitPattern();
 
   unsigned int counter ;
   TrackCollection::const_iterator ctfTkIter ;
   for ( ctfTkIter = ctfTrackCollection->begin() , counter = 0 ;
         ctfTkIter != ctfTrackCollection->end() ; ctfTkIter++, counter++ )
    {
 
     double dEta = gsfTrackRef->eta() - ctfTkIter->eta();
     double dPhi = gsfTrackRef->phi() - ctfTkIter->phi();
     double pi = acos(-1.);
     if(std::abs(dPhi) > pi) dPhi = 2*pi - std::abs(dPhi);
 
     // dont want to look at every single track in the event!
     if(sqrt(dEta*dEta + dPhi*dPhi) > 0.3) continue;
 
     unsigned int shared = 0 ;
     int gsfHitCounter = 0 ;
     int numGsfInnerHits = 0 ;
     int numCtfInnerHits = 0 ;
     // get the CTF Track Hit Pattern
     const HitPattern& ctfHitPattern = ctfTkIter->hitPattern();
 
     trackingRecHit_iterator elHitsIt ;
     for ( elHitsIt = gsfTrackRef->recHitsBegin() ;
           elHitsIt != gsfTrackRef->recHitsEnd() ;
           elHitsIt++, gsfHitCounter++ )
      {
       if(!((**elHitsIt).isValid()))  //count only valid Hits
        { continue ; }
 
       // look only in the pixels/TIB/TID
       uint32_t gsfHit = gsfHitPattern.getHitPattern(HitPattern::TRACK_HITS, gsfHitCounter);
       if (!(HitPattern::pixelHitFilter(gsfHit) ||
           HitPattern::stripTIBHitFilter(gsfHit) ||
           HitPattern::stripTIDHitFilter(gsfHit)))
        { continue ; }
 
       numGsfInnerHits++ ;
 
       int ctfHitsCounter = 0 ;
       numCtfInnerHits = 0 ;
       trackingRecHit_iterator ctfHitsIt ;
       for ( ctfHitsIt = ctfTkIter->recHitsBegin() ;
             ctfHitsIt != ctfTkIter->recHitsEnd() ;
             ctfHitsIt++, ctfHitsCounter++ )
        {
         if(!((**ctfHitsIt).isValid())) //count only valid Hits!
          { continue ; }
 
       uint32_t ctfHit = ctfHitPattern.getHitPattern(HitPattern::TRACK_HITS, ctfHitsCounter);
       if(!(HitPattern::pixelHitFilter(ctfHit) ||
             HitPattern::stripTIBHitFilter(ctfHit) ||
             HitPattern::stripTIDHitFilter(ctfHit)))
        { continue ; }
 
       numCtfInnerHits++ ;
 
         if( (**elHitsIt).sharesInput(&(**ctfHitsIt),TrackingRecHit::all) )
          {
           shared++ ;
           break ;
          }
 
        } //ctfHits iterator
 
      } //gsfHits iterator
 
     if ((numGsfInnerHits==0)||(numCtfInnerHits==0))
      { continue ; }
 
     if ( static_cast<float>(shared)/std::min(numGsfInnerHits,numCtfInnerHits) > maxFracShared )
      {
       maxFracShared = static_cast<float>(shared)/std::min(numGsfInnerHits, numCtfInnerHits);
       ctfTrackRef = TrackRef(ctfTracksH_,counter);
      }
 
    } //ctfTrack iterator
 
   return make_pair(ctfTrackRef,maxFracShared) ;
  }

void GsfElectronCoreBaseProducer::initEvent	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

protected

Definition at line 46 of file GsfElectronCoreBaseProducer.cc.

Referenced by GsfElectronCoreEcalDrivenProducer::produce(), GEDGsfElectronCoreProducer::produce(), and GsfElectronCoreProducer::produce().

  {
   if (useGsfPfRecTracks_)
    { event.getByToken(gsfPfRecTracksTag_,gsfPfRecTracksH_) ; }
   event.getByToken(gsfTracksTag_,gsfTracksH_) ;
   event.getByToken(ctfTracksTag_,ctfTracksH_) ;
  }