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

Inheritance diagram for SiStripElectronProducer:

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

	SiStripElectronProducer (const edm::ParameterSet &)

	~SiStripElectronProducer () 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)

Private Attributes
SiStripElectronAlgo *	algo_p

edm::EDGetTokenT< SiStripMatchedRecHit2DCollection >	matched_sistrips2dtag_

edm::EDGetTokenT< SiStripRecHit2DCollection >	rphi_sistrips2dtag_

std::string	siHitProducer_

std::string	siMatchedHitCollection_

std::string	siRphiHitCollection_

std::string	siStereoHitCollection_

std::string	siStripElectronsLabel_

edm::EDGetTokenT< SiStripRecHit2DCollection >	stereo_sistrips2dtag_

std::string	superClusterCollection_

std::string	superClusterProducer_

edm::EDGetTokenT< reco::SuperClusterCollection >	superClustertag_

std::string	trackCandidatesLabel_

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

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

Detailed Description

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

Usage: <usage>

Definition at line 35 of file SiStripElectronProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 46 of file SiStripElectronProducer.cc.

References algo_p, edm::ParameterSet::getParameter(), LogDebug, matched_sistrips2dtag_, rphi_sistrips2dtag_, siHitProducer_, siMatchedHitCollection_, siRphiHitCollection_, siStereoHitCollection_, siStripElectronsLabel_, stereo_sistrips2dtag_, AlCaHLTBitMon_QueryRunRegistry::string, superClusterCollection_, superClusterProducer_, superClustertag_, and trackCandidatesLabel_.

 {
    // register your products
    siStripElectronsLabel_ = iConfig.getParameter<std::string>("siStripElectronsLabel");
    trackCandidatesLabel_ = iConfig.getParameter<std::string>("trackCandidatesLabel");
    produces<reco::SiStripElectronCollection>(siStripElectronsLabel_);
    produces<TrackCandidateCollection>(trackCandidatesLabel_);
 
    // get parameters
    siHitProducer_ = iConfig.getParameter<std::string>("siHitProducer");
    siRphiHitCollection_ = iConfig.getParameter<std::string>("siRphiHitCollection");
    siStereoHitCollection_ = iConfig.getParameter<std::string>("siStereoHitCollection");
    siMatchedHitCollection_ = iConfig.getParameter<std::string>("siMatchedHitCollection");
 
    superClusterProducer_ = iConfig.getParameter<std::string>("superClusterProducer");
    superClusterCollection_ = iConfig.getParameter<std::string>("superClusterCollection");
    
    rphi_sistrips2dtag_ = 
      consumes<SiStripRecHit2DCollection>(edm::InputTag(siHitProducer_,siRphiHitCollection_));
    stereo_sistrips2dtag_ = 
      consumes<SiStripRecHit2DCollection>(edm::InputTag(siHitProducer_,siStereoHitCollection_));
    matched_sistrips2dtag_ = 
      consumes<SiStripMatchedRecHit2DCollection>(edm::InputTag(siHitProducer_,siMatchedHitCollection_));
    superClustertag_ = 
      consumes<reco::SuperClusterCollection>(edm::InputTag(superClusterProducer_,superClusterCollection_));
 
    algo_p = new SiStripElectronAlgo(
       iConfig.getParameter<int32_t>("maxHitsOnDetId"),
       iConfig.getParameter<double>("originUncertainty"),
       iConfig.getParameter<double>("phiBandWidth"),      // this is in radians
       iConfig.getParameter<double>("maxNormResid"),
       iConfig.getParameter<int32_t>("minHits"),
       iConfig.getParameter<double>("maxReducedChi2"));
 
    LogDebug("") << " Welcome to SiStripElectronProducer " ;
 
 }

SiStripElectronProducer::~SiStripElectronProducer ( )

override

Definition at line 90 of file SiStripElectronProducer.cc.

References algo_p.

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

Member Function Documentation

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

override

Definition at line 115 of file SiStripElectronProducer.cc.

References algo_p, reco::CaloCluster::energy(), SiStripElectronAlgo::findElectron(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, LogDebug, matched_sistrips2dtag_, eostools::move(), SiStripElectronAlgo::prepareEvent(), edm::Handle< T >::product(), edm::Event::put(), rphi_sistrips2dtag_, SimDataFormats::CaloAnalysis::sc, siStripElectronsLabel_, stereo_sistrips2dtag_, harvestTrackValidationPlots::str, superClustertag_, and trackCandidatesLabel_.

 {
    // Extract data from the event
    edm::ESHandle<TrackerGeometry> trackerHandle;
    iSetup.get<TrackerDigiGeometryRecord>().get(trackerHandle);
 
    edm::Handle<SiStripRecHit2DCollection> rphiHitsHandle;
    iEvent.getByToken(rphi_sistrips2dtag_, rphiHitsHandle);
 
    edm::Handle<SiStripRecHit2DCollection> stereoHitsHandle;
    iEvent.getByToken(stereo_sistrips2dtag_, stereoHitsHandle);
 
    edm::Handle<SiStripMatchedRecHit2DCollection> matchedHitsHandle;
    iEvent.getByToken(matched_sistrips2dtag_, matchedHitsHandle);
 
    edm::ESHandle<MagneticField> magneticFieldHandle;
    iSetup.get<IdealMagneticFieldRecord>().get(magneticFieldHandle);
 
    edm::Handle<reco::SuperClusterCollection> superClusterHandle;
    iEvent.getByToken(superClustertag_, superClusterHandle);
 
    // Set up SiStripElectronAlgo for this event
    algo_p->prepareEvent(trackerHandle, rphiHitsHandle, stereoHitsHandle, matchedHitsHandle, magneticFieldHandle);
 
    // Prepare the output electron candidates and clouds to be filled
    auto electronOut = std::make_unique<reco::SiStripElectronCollection>();
    auto trackCandidateOut = std::make_unique<TrackCandidateCollection>();
 
    //Retrieve tracker topology from geometry
    edm::ESHandle<TrackerTopology> tTopoHand;
    iSetup.get<TrackerTopologyRcd>().get(tTopoHand);
    const TrackerTopology *tTopo=tTopoHand.product();
 
    // counter for electron candidates
    int siStripElectCands = 0 ;
 
    std::ostringstream str;
 
 
    // Loop over clusters
    str << "Starting loop over superclusters."<< "\n" << std::endl;
    for (unsigned int i = 0;  i < superClusterHandle.product()->size();  i++) {
       const reco::SuperCluster* sc = &(*reco::SuperClusterRef(superClusterHandle, i));
       double energy = sc->energy();
 
       if (algo_p->findElectron(*electronOut, *trackCandidateOut, reco::SuperClusterRef(superClusterHandle, i),tTopo)) {
     str << "Supercluster energy: " << energy << ", FOUND an electron." << "\n" << std::endl;
      ++siStripElectCands ;
       }
       else {
      str << "Supercluster energy: " << energy << ", DID NOT FIND an electron."<< "\n" << std::endl;
       }
    }
    str << "Ending loop over superclusters." << "\n" << std::endl;
    
    str << " Found " << siStripElectCands 
             << " SiStripElectron Candidates before track fit " 
             << "\n" << std::endl ;
 
    LogDebug("SiStripElectronProducer") << str.str();
 
    // Put the electron candidates and the tracking trajectories into the event
    iEvent.put(std::move(electronOut), siStripElectronsLabel_);
    iEvent.put(std::move(trackCandidateOut), trackCandidatesLabel_);
 }