#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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

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, randomXiThetaGunProducer_cfi::energy, 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_, 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_);
 }