GsfElectronFull5x5Filler Class Reference

#include <GsfElectronFull5x5Filler.h>

Inheritance diagram for GsfElectronFull5x5Filler:

Public Member Functions
void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
void	calculateShowerShape_full5x5 (const reco::SuperClusterRef &theClus, bool pflow, reco::GsfElectron::ShowerShape &showerShape)
	GsfElectronFull5x5Filler (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~GsfElectronFull5x5Filler () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Private Attributes
edm::Handle< EcalRecHitCollection >	_ebRecHits
edm::EDGetTokenT< EcalRecHitCollection >	_ebRecHitsToken
edm::Handle< EcalRecHitCollection >	_eeRecHits
edm::EDGetTokenT< EcalRecHitCollection >	_eeRecHitsToken
const CaloGeometry *	_geometry
std::unique_ptr< ElectronHcalHelper >	_hcalHelper
std::unique_ptr< ElectronHcalHelper >	_hcalHelperPflow
edm::EDGetTokenT< reco::GsfElectronCollection >	_source
const CaloTopology *	_topology

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 13 of file GsfElectronFull5x5Filler.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 38 of file GsfElectronFull5x5Filler.cc.

References edm::ParameterSet::getParameter(), ElectronHcalHelper::Configuration::hcalTowers, ElectronHcalHelper::Configuration::hOverEConeSize, ElectronHcalHelper::Configuration::hOverEPtMin, funct::true, and ElectronHcalHelper::Configuration::useTowers.

 {
   _source = consumes<reco::GsfElectronCollection>(cfg.getParameter<edm::InputTag>("source"));
   
   ElectronHcalHelper::Configuration hcalCfg, hcalCfgPflow;
   
   hcalCfg.hOverEConeSize = cfg.getParameter<double>("hOverEConeSize") ;
   if (hcalCfg.hOverEConeSize>0)
     {
       hcalCfg.useTowers = true ;
       hcalCfg.hcalTowers = 
     consumes<CaloTowerCollection>(cfg.getParameter<edm::InputTag>("hcalTowers")) ;
       hcalCfg.hOverEPtMin = cfg.getParameter<double>("hOverEPtMin") ;
     }
   hcalCfgPflow.hOverEConeSize = cfg.getParameter<double>("hOverEConeSizePflow") ;
   if (hcalCfgPflow.hOverEConeSize>0)
     {
       hcalCfgPflow.useTowers = true ;
       hcalCfgPflow.hcalTowers = 
     consumes<CaloTowerCollection>(cfg.getParameter<edm::InputTag>("hcalTowers")) ;
       hcalCfgPflow.hOverEPtMin = cfg.getParameter<double>("hOverEPtMinPflow") ;
     }
   _hcalHelper.reset(new ElectronHcalHelper(hcalCfg));
   _hcalHelperPflow.reset(new ElectronHcalHelper(hcalCfg));
   
   _ebRecHitsToken = consumes<EcalRecHitCollection>(cfg.getParameter<edm::InputTag>("barrelRecHitCollectionTag"));
   _eeRecHitsToken = consumes<EcalRecHitCollection>(cfg.getParameter<edm::InputTag>("endcapRecHitCollectionTag"));
   
   produces<reco::GsfElectronCollection >();
 }

GsfElectronFull5x5Filler::~GsfElectronFull5x5Filler ( )

override

Definition at line 69 of file GsfElectronFull5x5Filler.cc.

 {}

Member Function Documentation

void GsfElectronFull5x5Filler::beginLuminosityBlock ( edm::LuminosityBlock const &  lb, edm::EventSetup const &  es  )

override

Definition at line 94 of file GsfElectronFull5x5Filler.cc.

References edm::EventSetup::get(), and edm::ESHandle< T >::product().

                                                  {
  edm::ESHandle<CaloGeometry> caloGeom ;
  edm::ESHandle<CaloTopology> caloTopo ;
  es.get<CaloGeometryRecord>().get(caloGeom);
  es.get<CaloTopologyRecord>().get(caloTopo);
  _geometry = caloGeom.product();
  _topology = caloTopo.product();
}

void GsfElectronFull5x5Filler::calculateShowerShape_full5x5	(	const reco::SuperClusterRef &	theClus,
		bool	pflow,
		reco::GsfElectron::ShowerShape &	showerShape
	)

Definition at line 104 of file GsfElectronFull5x5Filler.cc.

References EcalClusterToolsT< noZS >::covariances(), gamEcalExtractorBlocks_cff::detector, EcalClusterToolsT< noZS >::e1x5(), reco::GsfElectron::ShowerShape::e1x5, EcalClusterToolsT< noZS >::e2x5Max(), reco::GsfElectron::ShowerShape::e2x5Max, EcalClusterToolsT< noZS >::e3x3(), EcalClusterToolsT< noZS >::e5x5(), reco::GsfElectron::ShowerShape::e5x5, EcalBarrel, reco::GsfElectron::ShowerShape::hcalDepth1OverEcal, reco::GsfElectron::ShowerShape::hcalDepth1OverEcalBc, reco::GsfElectron::ShowerShape::hcalDepth2OverEcal, reco::GsfElectron::ShowerShape::hcalDepth2OverEcalBc, reco::GsfElectron::ShowerShape::hcalTowersBehindClusters, reco::CaloCluster::hitsAndFractions(), edm::isNotFinite(), EcalClusterToolsT< noZS >::localCovariances(), nullptr, reco::GsfElectron::ShowerShape::r9, reco::GsfElectron::ShowerShape::sigmaEtaEta, reco::GsfElectron::ShowerShape::sigmaIetaIeta, reco::GsfElectron::ShowerShape::sigmaIphiIphi, mathSSE::sqrt(), and DetId::subdetId().

 {
  const reco::CaloCluster & seedCluster = *(theClus->seed()) ;
  // temporary, till CaloCluster->seed() is made available
  DetId seedXtalId = seedCluster.hitsAndFractions()[0].first ;
  int detector = seedXtalId.subdetId() ;
  
  const EcalRecHitCollection * recHits = nullptr ;
  if (detector==EcalBarrel)
   {
    recHits = _ebRecHits.product() ;
   }
  else
   {
    recHits = _eeRecHits.product() ;
   }

  std::vector<float> covariances = noZS::EcalClusterTools::covariances(seedCluster,recHits,_topology,_geometry) ;
  std::vector<float> localCovariances = noZS::EcalClusterTools::localCovariances(seedCluster,recHits,_topology) ;
  showerShape.sigmaEtaEta = sqrt(covariances[0]) ;
  showerShape.sigmaIetaIeta = sqrt(localCovariances[0]) ;
  if (!edm::isNotFinite(localCovariances[2])) showerShape.sigmaIphiIphi = sqrt(localCovariances[2]) ;
  showerShape.e1x5 = noZS::EcalClusterTools::e1x5(seedCluster,recHits,_topology)  ;
  showerShape.e2x5Max = noZS::EcalClusterTools::e2x5Max(seedCluster,recHits,_topology)  ;
  showerShape.e5x5 = noZS::EcalClusterTools::e5x5(seedCluster,recHits,_topology) ;
  showerShape.r9 = noZS::EcalClusterTools::e3x3(seedCluster,recHits,_topology)/theClus->rawEnergy() ;

  if (pflow)
   {
    showerShape.hcalDepth1OverEcal = _hcalHelperPflow->hcalESumDepth1(*theClus)/theClus->energy() ;
    showerShape.hcalDepth2OverEcal = _hcalHelperPflow->hcalESumDepth2(*theClus)/theClus->energy() ;
    showerShape.hcalTowersBehindClusters = _hcalHelperPflow->hcalTowersBehindClusters(*theClus) ;
    showerShape.hcalDepth1OverEcalBc = _hcalHelperPflow->hcalESumDepth1BehindClusters(showerShape.hcalTowersBehindClusters)/showerShape.e5x5 ;
    showerShape.hcalDepth2OverEcalBc = _hcalHelperPflow->hcalESumDepth2BehindClusters(showerShape.hcalTowersBehindClusters)/showerShape.e5x5 ;
   }
  else
   {
    showerShape.hcalDepth1OverEcal = _hcalHelper->hcalESumDepth1(*theClus)/theClus->energy() ;
    showerShape.hcalDepth2OverEcal = _hcalHelper->hcalESumDepth2(*theClus)/theClus->energy() ;
    showerShape.hcalTowersBehindClusters = _hcalHelper->hcalTowersBehindClusters(*theClus) ;
    showerShape.hcalDepth1OverEcalBc = _hcalHelper->hcalESumDepth1BehindClusters(showerShape.hcalTowersBehindClusters)/showerShape.e5x5 ;
    showerShape.hcalDepth2OverEcalBc = _hcalHelper->hcalESumDepth2BehindClusters(showerShape.hcalTowersBehindClusters)/showerShape.e5x5 ;
   }
 }

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

override

Definition at line 73 of file GsfElectronFull5x5Filler.cc.

References reco::GsfElectron::full5x5_setShowerShape(), eostools::move(), MillePedeFileConverter_cfg::out, reco::GsfElectron::superCluster(), and groupFilesInBlocks::temp.

 {
   auto out = std::make_unique<reco::GsfElectronCollection>();
   
   edm::Handle<reco::GsfElectronCollection> eles;
   event.getByToken(_source,eles);
   
   event.getByToken(_ebRecHitsToken,_ebRecHits);
   event.getByToken(_eeRecHitsToken,_eeRecHits);

   for( const auto& ele : *eles ) {
     reco::GsfElectron temp(ele);
     reco::GsfElectron::ShowerShape full5x5_ss;
     calculateShowerShape_full5x5(temp.superCluster(),true,full5x5_ss);
     temp.full5x5_setShowerShape(full5x5_ss);
     out->push_back(temp);
   }
   
   event.put(std::move(out));
 }

Member Data Documentation

edm::Handle<EcalRecHitCollection> GsfElectronFull5x5Filler::_ebRecHits

private

Definition at line 31 of file GsfElectronFull5x5Filler.h.

edm::EDGetTokenT<EcalRecHitCollection> GsfElectronFull5x5Filler::_ebRecHitsToken

private

Definition at line 29 of file GsfElectronFull5x5Filler.h.

edm::Handle<EcalRecHitCollection> GsfElectronFull5x5Filler::_eeRecHits

private

Definition at line 32 of file GsfElectronFull5x5Filler.h.

edm::EDGetTokenT<EcalRecHitCollection> GsfElectronFull5x5Filler::_eeRecHitsToken

private

Definition at line 29 of file GsfElectronFull5x5Filler.h.

const CaloGeometry* GsfElectronFull5x5Filler::_geometry

private

Definition at line 34 of file GsfElectronFull5x5Filler.h.

std::unique_ptr<ElectronHcalHelper> GsfElectronFull5x5Filler::_hcalHelper

private

Definition at line 30 of file GsfElectronFull5x5Filler.h.

std::unique_ptr<ElectronHcalHelper> GsfElectronFull5x5Filler::_hcalHelperPflow

private

Definition at line 30 of file GsfElectronFull5x5Filler.h.

edm::EDGetTokenT<reco::GsfElectronCollection> GsfElectronFull5x5Filler::_source

private

Definition at line 28 of file GsfElectronFull5x5Filler.h.

const CaloTopology* GsfElectronFull5x5Filler::_topology

private

Definition at line 33 of file GsfElectronFull5x5Filler.h.