#include <GEDGsfElectronFinalizer.h>

Inheritance diagram for GEDGsfElectronFinalizer:

Public Member Functions
	GEDGsfElectronFinalizer (const edm::ParameterSet &)

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

	~GEDGsfElectronFinalizer ()

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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

	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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
std::unique_ptr < ModifyObjectValueBase >	gedRegression_

unsigned	nDeps_

std::string	outputCollectionLabel_

edm::EDGetTokenT < reco::PFCandidateCollection >	pfCandidates_

edm::EDGetTokenT < reco::GsfElectronCollection >	previousGsfElectrons_

std::vector< edm::EDGetTokenT < edm::ValueMap< double > > >	tokenElectronIsoVals_

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

Definition at line 18 of file GEDGsfElectronFinalizer.h.

Constructor & Destructor Documentation

GEDGsfElectronFinalizer::GEDGsfElectronFinalizer ( const edm::ParameterSet & cfg )

explicit

Definition at line 17 of file GEDGsfElectronFinalizer.cc.

References edm::ParameterSet::existsAs(), reco::get(), edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), fwrapper::plugin, and AlCaHLTBitMon_QueryRunRegistry::string.

  {   
    previousGsfElectrons_ = consumes<reco::GsfElectronCollection>(cfg.getParameter<edm::InputTag>("previousGsfElectronsTag"));
    pfCandidates_ = consumes<reco::PFCandidateCollection>(cfg.getParameter<edm::InputTag>("pfCandidatesTag"));
    outputCollectionLabel_ = cfg.getParameter<std::string>("outputCollectionLabel");
    edm::ParameterSet pfIsoVals(cfg.getParameter<edm::ParameterSet> ("pfIsolationValues"));
    
    tokenElectronIsoVals_.push_back(consumes<edm::ValueMap<double> >(pfIsoVals.getParameter<edm::InputTag>("pfSumChargedHadronPt")));
    tokenElectronIsoVals_.push_back(consumes<edm::ValueMap<double> >(pfIsoVals.getParameter<edm::InputTag>("pfSumPhotonEt")));
    tokenElectronIsoVals_.push_back(consumes<edm::ValueMap<double> >(pfIsoVals.getParameter<edm::InputTag>("pfSumNeutralHadronEt")));
    tokenElectronIsoVals_.push_back(consumes<edm::ValueMap<double> >(pfIsoVals.getParameter<edm::InputTag>("pfSumPUPt")));
 //   std::vector<std::string> isoNames = pfIsoVals.getParameterNamesForType<edm::InputTag>();
 //   for(const std::string& name : isoNames) {
 //     edm::InputTag tag = 
 //       pfIsoVals.getParameter<edm::InputTag>(name);
 //     tokenElectronIsoVals_.push_back(consumes<edm::ValueMap<double> >(tag));   
 //   }
 
    nDeps_ =  tokenElectronIsoVals_.size();
 
    if( cfg.existsAs<edm::ParameterSet>("regressionConfig") ) {
      const edm::ParameterSet& iconf = cfg.getParameterSet("regressionConfig");
      const std::string& mname = iconf.getParameter<std::string>("modifierName");
      ModifyObjectValueBase* plugin = 
        ModifyObjectValueFactory::get()->create(mname,iconf);
      gedRegression_.reset(plugin);
      edm::ConsumesCollector sumes = consumesCollector();
      gedRegression_->setConsumes(sumes);
    } else {
      gedRegression_.reset(nullptr);
    }
 
    produces<reco::GsfElectronCollection> (outputCollectionLabel_);
 }

GEDGsfElectronFinalizer::~GEDGsfElectronFinalizer ( )

Definition at line 52 of file GEDGsfElectronFinalizer.cc.

53 {}

Member Function Documentation

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

virtual

Implements edm::stream::EDProducerBase.

Definition at line 56 of file GEDGsfElectronFinalizer.cc.

References funct::abs(), reco::GsfElectron::gsfTrack(), i, reco::GsfElectron::mvaOutput(), reco::GsfElectron::setMvaOutput(), reco::GsfElectron::setPassPflowPreselection(), reco::GsfElectron::setPfIsolationVariables(), reco::GsfElectron::MvaOutput::status, reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::PflowIsolationVariables::sumPUPt.

  {
    
    // Output collection
    std::auto_ptr<reco::GsfElectronCollection> outputElectrons_p(new reco::GsfElectronCollection);
    
    if( gedRegression_ ) {
      gedRegression_->setEvent(event);
      gedRegression_->setEventContent(setup);
    }
 
    // read input collections
    // electrons
    edm::Handle<reco::GsfElectronCollection> gedElectronHandle;
    event.getByToken(previousGsfElectrons_,gedElectronHandle);
 
    // PFCandidates
    edm::Handle<reco::PFCandidateCollection> pfCandidateHandle;
    event.getByToken(pfCandidates_,pfCandidateHandle);
    // value maps
    std::vector< edm::Handle< edm::ValueMap<double> > > isolationValueMaps(nDeps_);
 
    for(unsigned i=0; i < nDeps_ ; ++i) {
      event.getByToken(tokenElectronIsoVals_[i],isolationValueMaps[i]);
    }
 
    // prepare a map of PFCandidates having a valid GsfTrackRef to save time
    std::map<reco::GsfTrackRef, const reco::PFCandidate* > gsfPFMap;
    reco::PFCandidateCollection::const_iterator it = pfCandidateHandle->begin();
    reco::PFCandidateCollection::const_iterator itend = pfCandidateHandle->end() ;
    for(;it!=itend;++it) {
      // First check that the GsfTrack is non null
      if( it->gsfTrackRef().isNonnull()) {
        if(abs(it->pdgId())==11) // consider only the electrons 
      gsfPFMap[it->gsfTrackRef()]=&(*it);
      }
    }
 
    
    // Now loop on the electrons
    unsigned nele=gedElectronHandle->size();
    for(unsigned iele=0; iele<nele;++iele) {
      reco::GsfElectronRef myElectronRef(gedElectronHandle,iele);
      
      reco::GsfElectron newElectron(*myElectronRef);
      reco::GsfElectron::PflowIsolationVariables isoVariables;
      isoVariables.sumChargedHadronPt = (*(isolationValueMaps)[0])[myElectronRef];
      isoVariables.sumPhotonEt = (*(isolationValueMaps)[1])[myElectronRef];
      isoVariables.sumNeutralHadronEt = (*(isolationValueMaps)[2])[myElectronRef];
      isoVariables.sumPUPt = (*(isolationValueMaps)[3])[myElectronRef];
      newElectron.setPfIsolationVariables(isoVariables);
 
      // now set a status if not already done (in GEDGsfElectronProducer.cc)
      //     std::cout << " previous status " << newElectron.mvaOutput().status << std::endl;
      if(newElectron.mvaOutput().status<=0) { 
        std::map<reco::GsfTrackRef, const  reco::PFCandidate * >::const_iterator itcheck=gsfPFMap.find(newElectron.gsfTrack());
        reco::GsfElectron::MvaOutput myMvaOutput(newElectron.mvaOutput());
        if(itcheck!=gsfPFMap.end()) {
      // it means that there is a PFCandidate with the same GsfTrack
      myMvaOutput.status = 3; //as defined in PFCandidateEGammaExtra.h
      newElectron.setPassPflowPreselection(true); //this is currently fully redundant with mvaOutput.stats so candidate for removal
        }
        else{
      myMvaOutput.status = 4 ; //
      newElectron.setPassPflowPreselection(false);//this is currently fully redundant with mvaOutput.stats so candidate for removal   
        }
        newElectron.setMvaOutput(myMvaOutput);
      }
      
      if( gedRegression_ ) {
        gedRegression_->modifyObject(newElectron);
      }
      outputElectrons_p->push_back(newElectron);
    }
    
    event.put(outputElectrons_p,outputCollectionLabel_);
  }

Member Data Documentation

std::unique_ptr<ModifyObjectValueBase> GEDGsfElectronFinalizer::gedRegression_

private

Definition at line 31 of file GEDGsfElectronFinalizer.h.

unsigned GEDGsfElectronFinalizer::nDeps_

private

Definition at line 32 of file GEDGsfElectronFinalizer.h.

std::string GEDGsfElectronFinalizer::outputCollectionLabel_

private

Definition at line 29 of file GEDGsfElectronFinalizer.h.

edm::EDGetTokenT<reco::PFCandidateCollection> GEDGsfElectronFinalizer::pfCandidates_

private

Definition at line 28 of file GEDGsfElectronFinalizer.h.

edm::EDGetTokenT<reco::GsfElectronCollection> GEDGsfElectronFinalizer::previousGsfElectrons_

private

Definition at line 27 of file GEDGsfElectronFinalizer.h.

std::vector<edm::EDGetTokenT<edm::ValueMap<double> > > GEDGsfElectronFinalizer::tokenElectronIsoVals_

private

Definition at line 30 of file GEDGsfElectronFinalizer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation