ElectronIDValueMapProducer Class Reference

Inheritance diagram for ElectronIDValueMapProducer:

Public Member Functions
	ElectronIDValueMapProducer (const edm::ParameterSet &)
	~ElectronIDValueMapProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override
void	writeValueMap (edm::Event &iEvent, const edm::Handle< edm::View< reco::GsfElectron > > &handle, const std::vector< float > &values, const std::string &label) const

Private Attributes
edm::EDGetTokenT< EcalRecHitCollection >	ebReducedRecHitCollection_
edm::EDGetTokenT< EcalRecHitCollection >	ebReducedRecHitCollectionMiniAOD_
edm::EDGetTokenT< EcalRecHitCollection >	eeReducedRecHitCollection_
edm::EDGetTokenT< EcalRecHitCollection >	eeReducedRecHitCollectionMiniAOD_
edm::EDGetTokenT< EcalRecHitCollection >	esReducedRecHitCollection_
edm::EDGetTokenT< EcalRecHitCollection >	esReducedRecHitCollectionMiniAOD_
noZS::EcalClusterLazyTools *	lazyToolnoZS
edm::EDGetToken	src_
edm::EDGetToken	srcMiniAOD_

Static Private Attributes
static char	eleFull5x5Circularity_ [] = "eleFull5x5Circularity"
static char	eleFull5x5E1x5_ [] = "eleFull5x5E1x5"
static char	eleFull5x5E2x5_ [] = "eleFull5x5E2x5"
static char	eleFull5x5E5x5_ [] = "eleFull5x5E5x5"
static char	eleFull5x5R9_ [] = "eleFull5x5R9"
static char	eleFull5x5SigmaIEtaIEta_ [] = "eleFull5x5SigmaIEtaIEta"
static char	eleFull5x5SigmaIEtaIPhi_ [] = "eleFull5x5SigmaIEtaIPhi"

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

Definition at line 21 of file ElectronIDValueMapProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 70 of file ElectronIDValueMapProducer.cc.

References ebReducedRecHitCollection_, ebReducedRecHitCollectionMiniAOD_, eeReducedRecHitCollection_, eeReducedRecHitCollectionMiniAOD_, eleFull5x5Circularity_, eleFull5x5E1x5_, eleFull5x5E2x5_, eleFull5x5E5x5_, eleFull5x5R9_, eleFull5x5SigmaIEtaIEta_, eleFull5x5SigmaIEtaIPhi_, esReducedRecHitCollection_, esReducedRecHitCollectionMiniAOD_, edm::ParameterSet::getParameter(), src_, and srcMiniAOD_.

                                                                                     {

  //
  // Declare consummables, handle both AOD and miniAOD case
  //
  ebReducedRecHitCollection_        = mayConsume<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>
                                       ("ebReducedRecHitCollection"));
  ebReducedRecHitCollectionMiniAOD_ = mayConsume<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>
                                       ("ebReducedRecHitCollectionMiniAOD"));

  eeReducedRecHitCollection_        = mayConsume<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>
                                       ("eeReducedRecHitCollection"));
  eeReducedRecHitCollectionMiniAOD_ = mayConsume<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>
                                       ("eeReducedRecHitCollectionMiniAOD"));

  esReducedRecHitCollection_        = mayConsume<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>
                                       ("esReducedRecHitCollection"));
  esReducedRecHitCollectionMiniAOD_ = mayConsume<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>
                                       ("esReducedRecHitCollectionMiniAOD"));

  src_        = mayConsume<edm::View<reco::GsfElectron> >(iConfig.getParameter<edm::InputTag>("src"));
  srcMiniAOD_ = mayConsume<edm::View<reco::GsfElectron> >(iConfig.getParameter<edm::InputTag>("srcMiniAOD"));

  produces<edm::ValueMap<float> >(eleFull5x5SigmaIEtaIEta_);  
  produces<edm::ValueMap<float> >(eleFull5x5SigmaIEtaIPhi_); 
  produces<edm::ValueMap<float> >(eleFull5x5E1x5_);
  produces<edm::ValueMap<float> >(eleFull5x5E2x5_);
  produces<edm::ValueMap<float> >(eleFull5x5E5x5_);
  produces<edm::ValueMap<float> >(eleFull5x5R9_);  
  produces<edm::ValueMap<float> >(eleFull5x5Circularity_);  

}

ElectronIDValueMapProducer::~ElectronIDValueMapProducer ( )

override

Definition at line 103 of file ElectronIDValueMapProducer.cc.

                                                        {
}

Member Function Documentation

void ElectronIDValueMapProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 187 of file ElectronIDValueMapProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

                                                                                            {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

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

overrideprivate

Definition at line 106 of file ElectronIDValueMapProducer.cc.

References heepElectronID_HEEPV50_CSA14_25ns_cff::e1x5, EcalClusterLazyToolsT< EcalClusterToolsImpl >::e1x5(), heepElectronID_HEEPV50_CSA14_25ns_cff::e2x5, EcalClusterLazyToolsT< EcalClusterToolsImpl >::e2x5Max(), EcalClusterLazyToolsT< EcalClusterToolsImpl >::e3x3(), heepElectronID_HEEPV50_CSA14_25ns_cff::e5x5, EcalClusterLazyToolsT< EcalClusterToolsImpl >::e5x5(), ebReducedRecHitCollection_, ebReducedRecHitCollectionMiniAOD_, eeReducedRecHitCollection_, eeReducedRecHitCollectionMiniAOD_, eleFull5x5Circularity_, eleFull5x5E1x5_, eleFull5x5E2x5_, eleFull5x5E5x5_, eleFull5x5R9_, eleFull5x5SigmaIEtaIEta_, eleFull5x5SigmaIEtaIPhi_, esReducedRecHitCollection_, esReducedRecHitCollectionMiniAOD_, edm::Event::getByToken(), mps_fire::i, edm::detail::isnan(), edm::HandleBase::isValid(), lazyToolnoZS, EcalClusterLazyToolsT< EcalClusterToolsImpl >::localCovariances(), mathSSE::sqrt(), TrackRefitter_38T_cff::src, src_, srcMiniAOD_, and writeValueMap().

                                                                                      {

  using namespace edm;
  
  edm::Handle<edm::View<reco::GsfElectron> > src;

  // Retrieve the collection of electrons from the event.
  // If we fail to retrieve the collection with the standard AOD
  // name, we next look for the one with the stndard miniAOD name.
  bool isAOD = true;
  iEvent.getByToken(src_, src);

  if( !src.isValid() ){
    isAOD = false;
    iEvent.getByToken(srcMiniAOD_,src);
  }

  if( isAOD )
    lazyToolnoZS = new noZS::EcalClusterLazyTools(iEvent, iSetup, 
                          ebReducedRecHitCollection_, 
                          eeReducedRecHitCollection_,
                          esReducedRecHitCollection_ );
  else
    lazyToolnoZS = new noZS::EcalClusterLazyTools(iEvent, iSetup, 
                          ebReducedRecHitCollectionMiniAOD_, 
                          eeReducedRecHitCollectionMiniAOD_,
                          esReducedRecHitCollectionMiniAOD_ );
 
  // size_t n = src->size();
  std::vector<float> eleFull5x5SigmaIEtaIEta, eleFull5x5SigmaIEtaIPhi;
  std::vector<float> eleFull5x5R9, eleFull5x5Circularity;
  std::vector<float> eleFull5x5E1x5,eleFull5x5E2x5,eleFull5x5E5x5;
  
  // reco::GsfElectron::superCluster() is virtual so we can exploit polymorphism
  for (size_t i = 0; i < src->size(); ++i){
    auto iEle = src->ptrAt(i);
    const auto& theseed = *(iEle->superCluster()->seed());

    std::vector<float> vCov = lazyToolnoZS->localCovariances( theseed );
    const float see = (isnan(vCov[0]) ? 0. : sqrt(vCov[0]));
    const float sep = vCov[1];
    eleFull5x5SigmaIEtaIEta.push_back(see);
    eleFull5x5SigmaIEtaIPhi.push_back(sep);
    eleFull5x5R9.push_back(lazyToolnoZS->e3x3( theseed ) / iEle->superCluster()->rawEnergy() );    
    
    const float e1x5 = lazyToolnoZS->e1x5( theseed );
    const float e2x5 = lazyToolnoZS->e2x5Max( theseed );
    const float e5x5 = lazyToolnoZS->e5x5( theseed );
    const float circularity = (e5x5 != 0.) ? 1.-e1x5/e5x5 : -1;
    
    eleFull5x5E1x5.push_back(e1x5); 
    eleFull5x5E2x5.push_back(e2x5);
    eleFull5x5E5x5.push_back(e5x5);
    eleFull5x5Circularity.push_back(circularity);
  }
  
  writeValueMap(iEvent, src, eleFull5x5SigmaIEtaIEta, eleFull5x5SigmaIEtaIEta_);  
  writeValueMap(iEvent, src, eleFull5x5SigmaIEtaIPhi, eleFull5x5SigmaIEtaIPhi_);  
  writeValueMap(iEvent, src, eleFull5x5R9, eleFull5x5R9_);  
  writeValueMap(iEvent, src, eleFull5x5E1x5, eleFull5x5E1x5_);  
  writeValueMap(iEvent, src, eleFull5x5E2x5, eleFull5x5E2x5_);   
  writeValueMap(iEvent, src, eleFull5x5E5x5, eleFull5x5E5x5_);  
  writeValueMap(iEvent, src, eleFull5x5Circularity, eleFull5x5Circularity_);  
  
  delete lazyToolnoZS;
}

void ElectronIDValueMapProducer::writeValueMap ( edm::Event &  iEvent, const edm::Handle< edm::View< reco::GsfElectron > > &  handle, const std::vector< float > &  values, const std::string &  label  ) const

private

Definition at line 173 of file ElectronIDValueMapProducer.cc.

References objects.autophobj::filler, cmsBatch::handle, eostools::move(), and edm::Event::put().

Referenced by produce().

{
  using namespace edm; 
  using namespace std;
  auto valMap = std::make_unique<ValueMap<float>>();
  edm::ValueMap<float>::Filler filler(*valMap);
  filler.insert(handle, values.begin(), values.end());
  filler.fill();
  iEvent.put(std::move(valMap), label);
}

Member Data Documentation

edm::EDGetTokenT<EcalRecHitCollection> ElectronIDValueMapProducer::ebReducedRecHitCollection_

private

Definition at line 42 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> ElectronIDValueMapProducer::ebReducedRecHitCollectionMiniAOD_

private

Definition at line 48 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> ElectronIDValueMapProducer::eeReducedRecHitCollection_

private

Definition at line 43 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> ElectronIDValueMapProducer::eeReducedRecHitCollectionMiniAOD_

private

Definition at line 49 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5Circularity_ = "eleFull5x5Circularity"

staticprivate

Definition at line 59 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5E1x5_ = "eleFull5x5E1x5"

staticprivate

Definition at line 55 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5E2x5_ = "eleFull5x5E2x5"

staticprivate

Definition at line 56 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5E5x5_ = "eleFull5x5E5x5"

staticprivate

Definition at line 57 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5R9_ = "eleFull5x5R9"

staticprivate

Definition at line 58 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5SigmaIEtaIEta_ = "eleFull5x5SigmaIEtaIEta"

staticprivate

Definition at line 53 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

char ElectronIDValueMapProducer::eleFull5x5SigmaIEtaIPhi_ = "eleFull5x5SigmaIEtaIPhi"

staticprivate

Definition at line 54 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> ElectronIDValueMapProducer::esReducedRecHitCollection_

private

Definition at line 44 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> ElectronIDValueMapProducer::esReducedRecHitCollectionMiniAOD_

private

Definition at line 50 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

noZS::EcalClusterLazyTools* ElectronIDValueMapProducer::lazyToolnoZS

private

Definition at line 39 of file ElectronIDValueMapProducer.cc.

Referenced by produce().

edm::EDGetToken ElectronIDValueMapProducer::src_

private

Definition at line 45 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().

edm::EDGetToken ElectronIDValueMapProducer::srcMiniAOD_

private

Definition at line 51 of file ElectronIDValueMapProducer.cc.

Referenced by ElectronIDValueMapProducer(), and produce().