#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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

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

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

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

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, particleFlow_cfi::dEta, particleFlow_cfi::dPhi, 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_) ;
  }