#include <RecoEgamma/EgammaElectronProducers/interface/SiStripElectronAssociator.h>

Inheritance diagram for SiStripElectronAssociator:

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

	SiStripElectronAssociator (const edm::ParameterSet &)

	~SiStripElectronAssociator () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
edm::InputTag	electronsLabel_

edm::EDGetTokenT< reco::SiStripElectronCollection >	siStripElectronCollection_

edm::EDGetTokenT< reco::TrackCollection >	trackCollection_

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>="">

Description: <one line="" class="" summary>="">

Usage: <usage>

Definition at line 42 of file SiStripElectronAssociator.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 50 of file SiStripElectronAssociator.cc.

References electronsLabel_, edm::ParameterSet::getParameter(), edm::InputTag::label(), siStripElectronCollection_, and trackCollection_.

 {
    //register your products
   electronsLabel_ = iConfig.getParameter<edm::InputTag>("electronsLabel");
   produces<reco::ElectronCollection>(electronsLabel_.label());
 
    //now do what ever other initialization is needed
   siStripElectronCollection_ = consumes<reco::SiStripElectronCollection>(iConfig.getParameter<edm::InputTag>("siStripElectronCollection"));
   trackCollection_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("trackCollection"));
 }

SiStripElectronAssociator::~SiStripElectronAssociator ( )

override

Definition at line 62 of file SiStripElectronAssociator.cc.

63 {}

Member Function Documentation

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

override

Definition at line 67 of file SiStripElectronAssociator.cc.

References reco::TrackBase::charge(), MillePedeFileConverter_cfg::e, metsig::electron, electronsLabel_, Exception, edm::Event::getByToken(), mps_fire::i, edm::InputTag::label(), LogDebug, mag(), eostools::move(), convertSQLitetoXML_cfg::output, reco::TrackBase::p(), edm::Handle< T >::product(), edm::Event::put(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), reco::Track::recHits(), reco::Electron::setSuperCluster(), reco::Electron::setTrack(), siStripElectronCollection_, siStripElectrons_cfi::siStripElectrons, trackCollection_, l1t::tracks, reco::TrackBase::vx(), reco::TrackBase::vy(), and reco::TrackBase::vz().

 {
 
   // position tolerance for equality of 2 hits set to 10 microns
   static const double positionTol = 1e-3 ; 
 
    edm::Handle<reco::SiStripElectronCollection> siStripElectrons;
    iEvent.getByToken(siStripElectronCollection_, siStripElectrons);
 
    edm::Handle<reco::TrackCollection> tracks;
    iEvent.getByToken(trackCollection_, tracks);
 
    std::map<const reco::SiStripElectron*, bool> alreadySeen;
    for (reco::SiStripElectronCollection::const_iterator strippyIter = siStripElectrons->begin();  strippyIter != siStripElectrons->end();  ++strippyIter) {
       alreadySeen[&(*strippyIter)] = false;
    }
 
    // Output the high-level Electrons
    auto output = std::make_unique<reco::ElectronCollection>();
 
    LogDebug("SiStripElectronAssociator") << " About to loop over tracks " << std::endl ;
    LogDebug("SiStripElectronAssociator") << " Number of tracks in Associator " << tracks.product()->size() ;
    LogDebug("SiStripElectronAssociator") << " Number of SiStripElectron Candidates " << siStripElectrons->size();
 
    // The reco::Track's hits are a (improper?) subset of the reco::SiStripElectron's
    // countSiElFit counts electron candidates that return w/ a good track fit
    int countSiElFit = 0 ;
    for (unsigned int i = 0;  i < tracks.product()->size();  i++) {
       const reco::Track* trackPtr = &(*reco::TrackRef(tracks, i));
       
       // If the reco::Track and the reco::SiStripElectron share even
       // one hit in common, they belong to each other.  (Disjoint sets
       // of hits are assigned to electrons.)  So let's look at one hit.
 
       // But first, make sure the track's hit list is not empty.
       if (trackPtr->recHits().empty()) { continue; }
 
       // Detector id is not enough to completely specify a hit
       uint32_t id = (*trackPtr->recHits().begin())->geographicalId().rawId();
       LocalPoint pos = (*trackPtr->recHits().begin())->localPosition();
       
       LogDebug("SiStripElectronAssociator") << " New Track Candidate " << i
                                             << " DetId " << id
                                             << " pos " << pos << "\n";
 
       // Find the electron with that hit!
       bool foundElectron = false;
       for (reco::SiStripElectronCollection::const_iterator strippyIter = siStripElectrons->begin();  strippyIter != siStripElectrons->end();  ++strippyIter) {
         if (!alreadySeen[&(*strippyIter)]) {
           
           bool hitInCommon = false;
           LogDebug("SiStripElectronAssociator") << " Looping over Mono hits " << "\n" ;
           
           for (std::vector<SiStripRecHit2D>::const_iterator hitIter = strippyIter->rphiRecHits().begin();  hitIter != strippyIter->rphiRecHits().end();  ++hitIter) {
             
             LogDebug("SiStripElectronAssociator") << " SiStripCand " 
                                                   << " DetId " << hitIter->geographicalId().rawId()
                                                   << " localPos " << hitIter->localPosition()
                                                   << " deltasPos " << (hitIter->localPosition() - pos).mag() ;
             
             if (hitIter->geographicalId().rawId() == id   &&
                 (hitIter->localPosition() - pos).mag() < positionTol ) {
               hitInCommon = true;
               LogDebug("SiStripElectronAssociator") << " hitInCommon True " << "\n" ;
               break;
             } else {
               LogDebug("SiStripElectronAssociator") << " hitInCommon False " << "\n" ;
             }
           } // end loop over rphi hits
           
           if(!hitInCommon) {
             LogDebug("SiStripElectronAssociator") << " Looping over Stereo hits " << "\n" ;
 
             for (std::vector<SiStripRecHit2D>::const_iterator hitIter = strippyIter->stereoRecHits().begin();  hitIter != strippyIter->stereoRecHits().end();  ++hitIter) {
               
                LogDebug("SiStripElectronAssociator") << " SiStripCand " 
                                                      << " DetId " << hitIter->geographicalId().rawId()
                                                      << " localPos " << hitIter->localPosition()
                                                      << " deltasPos " << (hitIter->localPosition() - pos).mag() ;
                
                if (hitIter->geographicalId().rawId() == id   &&
                    (hitIter->localPosition() - pos).mag() < positionTol) {
                  hitInCommon = true;
                   LogDebug("SiStripElectronAssociator") << " hitInCommon True " << "\n"  ;
                  break;
                } else {
                  LogDebug("SiStripElectronAssociator") << " hitInCommon False " << "\n" ;
                }
                
             } // end loop over stereo hits
           } // end of hitInCommon check for loop over stereo hits
           
           if (hitInCommon) {
              LogDebug("SiStripElectronAssociator") << " Hit in Common Found \n" ;
              ++countSiElFit ;
              foundElectron = true;
              alreadySeen[&(*strippyIter)] = true;
              
              reco::Electron electron((trackPtr->charge() > 0 ? 1 : -1),
                                      math::XYZTLorentzVector(trackPtr->px(),
                                                              trackPtr->py(),
                                                              trackPtr->pz(),
                                                              trackPtr->p()),
                                      math::XYZPoint(trackPtr->vx(),
                                                     trackPtr->vy(),
                                                     trackPtr->vz()));
              electron.setSuperCluster(strippyIter->superCluster());
              electron.setTrack(reco::TrackRef(tracks, i));
              
              output->push_back(electron);
              break ; // no need to check other SiStripElectrons if this one is good 
              
           } else {
              LogDebug("SiStripElectronAssociator") << "Hit in Common NOT found \n"  ;
           }
           // endif this electron belongs to this track
           
           
         } // endif we haven't seen this electron before
         LogDebug("SiStripElectronAssociator") << "Done with this electron " << "\n\n\n";
       } // end loop over electrons
       
       LogDebug("SiStripElectronAssociator") << "Testing if foundElectron " << foundElectron << std::endl;
       
       if (!foundElectron) {
         throw cms::Exception("Configuration")
           << " Inconsistent track collection!"; 
       }
       
       LogDebug("SiStripElectronAssociator") << "At end of track loop \n" << std::endl; 
       
    } // end loop over tracks
    
    
    LogDebug("SiStripElectronAssociator") << " Number of SiStripElectrons returned with a good fit " 
                                          << countSiElFit << "\n"<<  std::endl ;
    iEvent.put(std::move(output), electronsLabel_.label());
 }