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

Inheritance diagram for SiStripElectronProducer:

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

	SiStripElectronProducer (const edm::ParameterSet &)

	~SiStripElectronProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

Definition at line 89 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
	)

virtual

Implements edm::EDProducer.

Definition at line 114 of file SiStripElectronProducer.cc.

References algo_p, relval_parameters_module::energy, reco::CaloCluster::energy(), SiStripElectronAlgo::findElectron(), edm::EventSetup::get(), edm::Event::getByToken(), i, LogDebug, matched_sistrips2dtag_, SiStripElectronAlgo::prepareEvent(), edm::Handle< T >::product(), edm::ESHandle< class >::product(), edm::Event::put(), rphi_sistrips2dtag_, siStripElectronsLabel_, stereo_sistrips2dtag_, 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
    std::auto_ptr<reco::SiStripElectronCollection> electronOut(new reco::SiStripElectronCollection);
    std::auto_ptr<TrackCandidateCollection> trackCandidateOut(new TrackCandidateCollection);
 
    //Retrieve tracker topology from geometry
    edm::ESHandle<TrackerTopology> tTopoHand;
    iSetup.get<IdealGeometryRecord>().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(electronOut, siStripElectronsLabel_);
    iEvent.put(trackCandidateOut, trackCandidatesLabel_);
 }