GsfElectronProducer Class Reference

#include <GsfElectronProducer.h>

Inheritance diagram for GsfElectronProducer:

Public Member Functions
	GsfElectronProducer (const edm::ParameterSet &, const gsfAlgoHelpers::HeavyObjectCache *)
void	produce (edm::Event &, const edm::EventSetup &) override
Public Member Functions inherited from GsfElectronBaseProducer
	GsfElectronBaseProducer (const edm::ParameterSet &, const gsfAlgoHelpers::HeavyObjectCache *)
	~GsfElectronBaseProducer () override
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< gsfAlgoHelpers::HeavyObjectCache > >
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Protected Member Functions
void	beginEvent (edm::Event &, const edm::EventSetup &)
Protected Member Functions inherited from GsfElectronBaseProducer
void	beginEvent (edm::Event &, const edm::EventSetup &)
void	fillEvent (reco::GsfElectronCollection &electrons, edm::Event &event)
const edm::OrphanHandle< reco::GsfElectronCollection > &	orphanHandle () const

Private Member Functions
void	addPflowInfo (reco::GsfElectronCollection &electrons, edm::Event const &event) const
void	checkPfTranslatorParameters (edm::ParameterSet const &)
reco::GsfElectronCollection	clonePreviousElectrons (edm::Event const &event) const
void	setPflowPreselectionFlag (reco::GsfElectron &ele) const

Private Attributes
bool	pfTranslatorParametersChecked_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< gsfAlgoHelpers::HeavyObjectCache > >
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
Static Public Member Functions inherited from GsfElectronBaseProducer
static void	fillDescriptions (edm::ConfigurationDescriptions &)
static void	globalEndJob (gsfAlgoHelpers::HeavyObjectCache const *)
static std::unique_ptr< gsfAlgoHelpers::HeavyObjectCache >	initializeGlobalCache (const edm::ParameterSet &conf)
Protected Attributes inherited from GsfElectronBaseProducer
GsfElectronAlgo *	algo_
const GsfElectronAlgo::CutsConfiguration	cutsCfg_
const GsfElectronAlgo::CutsConfiguration	cutsCfgPflow_
ElectronHcalHelper::Configuration	hcalCfg_
ElectronHcalHelper::Configuration	hcalCfgPflow_
GsfElectronAlgo::InputTagsConfiguration	inputCfg_
edm::EDGetTokenT< edm::ValueMap< float > >	pfMVA_
GsfElectronAlgo::StrategyConfiguration	strategyCfg_

Detailed Description

Definition at line 7 of file GsfElectronProducer.h.

Constructor & Destructor Documentation

GsfElectronProducer::GsfElectronProducer ( const edm::ParameterSet &  cfg, const gsfAlgoHelpers::HeavyObjectCache *  hoc  )

explicit

Definition at line 26 of file GsfElectronProducer.cc.

  : GsfElectronBaseProducer(cfg,hoc), pfTranslatorParametersChecked_(false)
 {}

Member Function Documentation

void GsfElectronProducer::addPflowInfo ( reco::GsfElectronCollection &  electrons, edm::Event const &  event  ) const

private

Definition at line 116 of file GsfElectronProducer.cc.

References GsfElectronAlgo::InputTagsConfiguration::edIsoVals, relativeConstraints::empty, edm::ParameterSet::empty(), Exception, runEdmFileComparison::found, edm::ParameterSet::getParameter(), GsfElectronBaseProducer::inputCfg_, gedGsfElectronFinalizer_cfi::pfIsolationValues, GsfElectronAlgo::InputTagsConfiguration::pfIsoVals, GsfElectronAlgo::InputTagsConfiguration::pflowGsfElectronsTag, GsfElectronAlgo::InputTagsConfiguration::previousGsfElectrons, setPflowPreselectionFlag(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, and reco::GsfElectron::PflowIsolationVariables::sumPhotonEt.

Referenced by produce().

{
  //Isolation Value Maps for PF and EcalDriven electrons
  typedef std::vector<edm::Handle<edm::ValueMap<double> > > IsolationValueMaps;
  IsolationValueMaps pfIsolationValues;
  IsolationValueMaps edIsolationValues;

  //Fill in the Isolation Value Maps for PF and EcalDriven electrons
  std::vector<edm::InputTag> inputTagIsoVals;
  if (!inputCfg_.pfIsoVals.empty()) {
    inputTagIsoVals.push_back(inputCfg_.pfIsoVals.getParameter<edm::InputTag>("pfSumChargedHadronPt"));
    inputTagIsoVals.push_back(inputCfg_.pfIsoVals.getParameter<edm::InputTag>("pfSumPhotonEt"));
    inputTagIsoVals.push_back(inputCfg_.pfIsoVals.getParameter<edm::InputTag>("pfSumNeutralHadronEt"));

    pfIsolationValues.resize(inputTagIsoVals.size());

    for (size_t j = 0; j < inputTagIsoVals.size(); ++j) {
      event.getByLabel(inputTagIsoVals[j], pfIsolationValues[j]);
    }
  }

  if (!inputCfg_.edIsoVals.empty()) {
    inputTagIsoVals.clear();
    inputTagIsoVals.push_back(inputCfg_.edIsoVals.getParameter<edm::InputTag>("edSumChargedHadronPt"));
    inputTagIsoVals.push_back(inputCfg_.edIsoVals.getParameter<edm::InputTag>("edSumPhotonEt"));
    inputTagIsoVals.push_back(inputCfg_.edIsoVals.getParameter<edm::InputTag>("edSumNeutralHadronEt"));

    edIsolationValues.resize(inputTagIsoVals.size());

    for (size_t j = 0; j < inputTagIsoVals.size(); ++j) {
      event.getByLabel(inputTagIsoVals[j], edIsolationValues[j]);
    }
  }

  bool found;
  auto edElectrons = event.getHandle(inputCfg_.previousGsfElectrons);
  auto pfElectrons = event.getHandle(inputCfg_.pflowGsfElectronsTag);
  reco::GsfElectronCollection::const_iterator pfElectron, edElectron;
  unsigned int edIndex, pfIndex;

  for (auto& el : electrons) {
    // Retreive info from pflow electrons
    found = false;
    for (pfIndex = 0, pfElectron = pfElectrons->begin(); pfElectron != pfElectrons->end(); pfIndex++, pfElectron++) {
      if (pfElectron->gsfTrack() == el.gsfTrack()) {
        if (found) {
          edm::LogWarning("GsfElectronProducer") << "associated pfGsfElectron already found";
        } else {
          found = true;

          // Isolation Values
          if (!(pfIsolationValues).empty()) {
            reco::GsfElectronRef pfElectronRef(pfElectrons, pfIndex);
            reco::GsfElectron::PflowIsolationVariables isoVariables;
            isoVariables.sumChargedHadronPt = (*(pfIsolationValues)[0])[pfElectronRef];
            isoVariables.sumPhotonEt = (*(pfIsolationValues)[1])[pfElectronRef];
            isoVariables.sumNeutralHadronEt = (*(pfIsolationValues)[2])[pfElectronRef];
            el.setPfIsolationVariables(isoVariables);
          }

          // el.setPfIsolationVariables(pfElectron->pfIsolationVariables()) ;
          el.setMvaInput(pfElectron->mvaInput());
          el.setMvaOutput(pfElectron->mvaOutput());
          if (el.ecalDrivenSeed()) {
            el.setP4(GsfElectron::P4_PFLOW_COMBINATION, pfElectron->p4(GsfElectron::P4_PFLOW_COMBINATION),
                     pfElectron->p4Error(GsfElectron::P4_PFLOW_COMBINATION), false);
          } else {
            el.setP4(GsfElectron::P4_PFLOW_COMBINATION, pfElectron->p4(GsfElectron::P4_PFLOW_COMBINATION),
                     pfElectron->p4Error(GsfElectron::P4_PFLOW_COMBINATION), true);
          }
          double noCutMin = -999999999.;
          if (el.mva_e_pi() < noCutMin) {
            throw cms::Exception("GsfElectronAlgo|UnexpectedMvaValue") << "unexpected MVA value: " << el.mva_e_pi();
          }
        }
      }
    }

    // Isolation Values
    // Retreive not found info from ed electrons
    if (!(edIsolationValues).empty()) {
      edIndex = 0, edElectron = edElectrons->begin();
      while ((found == false) && (edElectron != edElectrons->end())) {
        if (edElectron->gsfTrack() == el.gsfTrack()) {
          found = true;

          // CONSTRUCTION D UNE REF dans le handle previousElectrons avec l'indice edIndex,
          // puis recuperation dans la ValueMap ED

          reco::GsfElectronRef edElectronRef(edElectrons, edIndex);
          reco::GsfElectron::PflowIsolationVariables isoVariables;
          isoVariables.sumChargedHadronPt = (*(edIsolationValues)[0])[edElectronRef];
          isoVariables.sumPhotonEt = (*(edIsolationValues)[1])[edElectronRef];
          isoVariables.sumNeutralHadronEt = (*(edIsolationValues)[2])[edElectronRef];
          el.setPfIsolationVariables(isoVariables);
        }

        edIndex++;
        edElectron++;
      }
    }

    // Preselection
    setPflowPreselectionFlag(el);
  }
}

void GsfElectronProducer::beginEvent ( edm::Event &  event, const edm::EventSetup &  setup  )

protected

Definition at line 64 of file GsfElectronProducer.cc.

References GsfElectronBaseProducer::beginEvent(), checkPfTranslatorParameters(), edm::parameterSet(), GsfElectronBaseProducer::pfMVA_, pfTranslatorParametersChecked_, edm::HandleBase::provenance(), and funct::true.

 {
  // extra configuration checks
  if (!pfTranslatorParametersChecked_)
   {
    pfTranslatorParametersChecked_ = true ;
    edm::Handle<edm::ValueMap<float> > pfMva ;
    event.getByToken(pfMVA_,pfMva) ;
    checkPfTranslatorParameters(edm::parameterSet(*pfMva.provenance())) ;
   }

  // call to base class
  return GsfElectronBaseProducer::beginEvent(event,setup) ;
 }

void GsfElectronProducer::checkPfTranslatorParameters ( edm::ParameterSet const &  pset )

private

Definition at line 79 of file GsfElectronProducer.cc.

References GsfElectronAlgo::StrategyConfiguration::applyPreselection, GsfElectronBaseProducer::cutsCfg_, GsfElectronBaseProducer::cutsCfgPflow_, Exception, edm::ParameterSet::getParameter(), GsfElectronAlgo::CutsConfiguration::minMVA, and GsfElectronBaseProducer::strategyCfg_.

Referenced by beginEvent().

 {
  edm::ParameterSet mvaBlock = pset.getParameter<edm::ParameterSet>("MVACutBlock") ;
  double pfTranslatorMinMva = mvaBlock.getParameter<double>("MVACut") ;
  double pfTranslatorUndefined = -99. ;
  if (strategyCfg_.applyPreselection&&(cutsCfgPflow_.minMVA<pfTranslatorMinMva))
   {
    // For pure tracker seeded electrons, if MVA is under translatorMinMva, there is no supercluster
    // of any kind available, so GsfElectronCoreProducer has already discarded the electron.
    edm::LogWarning("GsfElectronAlgo|MvaCutTooLow")
      <<"Parameter minMVAPflow ("<<cutsCfgPflow_.minMVA<<") will have no effect on purely tracker seeded electrons."
      <<" It is inferior to the cut already applied by PFlow translator ("<<pfTranslatorMinMva<<")." ;
   }
  if (strategyCfg_.applyPreselection&&(cutsCfg_.minMVA<pfTranslatorMinMva))
   {
    // For ecal seeded electrons, there is a cluster and GsfElectronCoreProducer has kept all electrons,
    // but when MVA is under translatorMinMva, the translator has not stored the supercluster and
    // forced the MVA value to translatorUndefined
    if (cutsCfg_.minMVA>pfTranslatorUndefined)
     {
      edm::LogWarning("GsfElectronAlgo|IncompletePflowInformation")
        <<"Parameter minMVA  ("<<cutsCfg_.minMVA<<")is inferior to the cut applied by PFlow translator ("<<pfTranslatorMinMva<<")."
        <<" Some ecal (and eventually tracker) seeded electrons may lack their MVA value and PFlow supercluster." ;
     }
    else
     {
      // the MVA value has been forced to translatorUndefined, inferior minMVAPflow
      // so the cut actually applied is the PFlow one
      throw cms::Exception("GsfElectronAlgo|BadMvaCut")
        <<"Parameter minMVA is inferior to the lowest possible value."
        <<" Every electron will be blessed whatever other criteria." ;
     }
   }
 }

reco::GsfElectronCollection GsfElectronProducer::clonePreviousElectrons ( edm::Event const &  event ) const

private

Definition at line 31 of file GsfElectronProducer.cc.

References lowPtElectrons_cff::electrons, edm::Event::get(), GsfElectronAlgo::InputTagsConfiguration::gsfElectronCores, GsfElectronBaseProducer::inputCfg_, and GsfElectronAlgo::InputTagsConfiguration::previousGsfElectrons.

Referenced by produce().

{
  reco::GsfElectronCollection electrons;

  auto coreElectrons = event.getHandle(inputCfg_.gsfElectronCores);
  const GsfElectronCoreCollection* newCores = coreElectrons.product();

  for (auto const& oldElectron : event.get(inputCfg_.previousGsfElectrons)) {
    const GsfElectronCoreRef oldCoreRef = oldElectron.core();
    const GsfTrackRef oldElectronGsfTrackRef = oldCoreRef->gsfTrack();
    unsigned int icore;
    for (icore = 0; icore < newCores->size(); ++icore) {
      if (oldElectronGsfTrackRef == (*newCores)[icore].gsfTrack()) {
        const GsfElectronCoreRef coreRef = edm::Ref<GsfElectronCoreCollection>(coreElectrons, icore);
        electrons.emplace_back(oldElectron, coreRef);
        break;
      }
    }
  }
  return electrons;
}

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

override

Definition at line 54 of file GsfElectronProducer.cc.

References GsfElectronAlgo::StrategyConfiguration::addPflowElectrons, addPflowInfo(), GsfElectronBaseProducer::algo_, clonePreviousElectrons(), GsfElectronAlgo::completeElectrons(), lowPtElectrons_cff::electrons, GsfElectronBaseProducer::fillEvent(), and GsfElectronBaseProducer::strategyCfg_.

 {
  auto electrons = clonePreviousElectrons(event) ;
  // don't add pflow only electrons if one so wish
  if (strategyCfg_.addPflowElectrons)
    { algo_->completeElectrons(electrons, event, setup, globalCache()) ; }
  addPflowInfo(electrons, event) ;
  fillEvent(electrons, event) ;
 }

void GsfElectronProducer::setPflowPreselectionFlag ( reco::GsfElectron &  ele ) const

private

Definition at line 224 of file GsfElectronProducer.cc.

References reco::GsfElectron::core(), GsfElectronBaseProducer::cutsCfg_, GsfElectronBaseProducer::cutsCfgPflow_, LogTrace, GsfElectronAlgo::CutsConfiguration::minMVA, reco::GsfElectron::MvaOutput::mva_e_pi, reco::GsfElectron::mvaOutput(), reco::GsfElectron::passingMvaPreselection(), reco::GsfElectron::setPassMvaPreselection(), and reco::GsfElectron::setPassPflowPreselection().

Referenced by addPflowInfo().

{
  ele.setPassMvaPreselection(false);

  if (ele.core()->ecalDrivenSeed()) {
    if (ele.mvaOutput().mva_e_pi >= cutsCfg_.minMVA)
      ele.setPassMvaPreselection(true);
  } else {
    if (ele.mvaOutput().mva_e_pi >= cutsCfgPflow_.minMVA)
      ele.setPassMvaPreselection(true);
  }

  if (ele.passingMvaPreselection()) {
    LogTrace("GsfElectronAlgo") << "Main mva criterion is satisfied";
  }

  ele.setPassPflowPreselection(ele.passingMvaPreselection());
}

Member Data Documentation

bool GsfElectronProducer::pfTranslatorParametersChecked_

private

Definition at line 25 of file GsfElectronProducer.h.

Referenced by beginEvent().