#include <EgammaEcalRecHitIsolationProducer.h>

Inheritance diagram for EgammaEcalRecHitIsolationProducer:

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

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

	~EgammaEcalRecHitIsolationProducer () override

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometrytoken_

edm::ParameterSet	conf_

edm::InputTag	ecalBarrelRecHitCollection_

edm::InputTag	ecalBarrelRecHitProducer_

edm::InputTag	ecalEndcapRecHitCollection_

edm::InputTag	ecalEndcapRecHitProducer_

double	egIsoConeSizeInBarrel_

double	egIsoConeSizeInEndcap_

double	egIsoConeSizeOut_

double	egIsoEMinBarrel_

double	egIsoEMinEndcap_

double	egIsoJurassicWidth_

double	egIsoPtMinBarrel_

double	egIsoPtMinEndcap_

edm::InputTag	emObjectProducer_

edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	sevLvToken_

bool	subtract_

bool	tryBoth_

bool	useIsolEt_

bool	useNumCrystals_

bool	vetoClustered_

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 40 of file EgammaEcalRecHitIsolationProducer.h.

Constructor & Destructor Documentation

◆ EgammaEcalRecHitIsolationProducer()

EgammaEcalRecHitIsolationProducer::EgammaEcalRecHitIsolationProducer ( const edm::ParameterSet & config )

explicit

Definition at line 24 of file EgammaEcalRecHitIsolationProducer.cc.

                                                                                                   : conf_(config) {
   // use configuration file to setup input/output collection names
   //inputs
   emObjectProducer_ = conf_.getParameter<edm::InputTag>("emObjectProducer");
   ecalBarrelRecHitProducer_ = conf_.getParameter<edm::InputTag>("ecalBarrelRecHitProducer");
   ecalBarrelRecHitCollection_ = conf_.getParameter<edm::InputTag>("ecalBarrelRecHitCollection");
   ecalEndcapRecHitProducer_ = conf_.getParameter<edm::InputTag>("ecalEndcapRecHitProducer");
   ecalEndcapRecHitCollection_ = conf_.getParameter<edm::InputTag>("ecalEndcapRecHitCollection");
  
   //vetos
   egIsoPtMinBarrel_ = conf_.getParameter<double>("etMinBarrel");
   egIsoEMinBarrel_ = conf_.getParameter<double>("eMinBarrel");
   egIsoPtMinEndcap_ = conf_.getParameter<double>("etMinEndcap");
   egIsoEMinEndcap_ = conf_.getParameter<double>("eMinEndcap");
   egIsoConeSizeInBarrel_ = conf_.getParameter<double>("intRadiusBarrel");
   egIsoConeSizeInEndcap_ = conf_.getParameter<double>("intRadiusEndcap");
   egIsoConeSizeOut_ = conf_.getParameter<double>("extRadius");
   egIsoJurassicWidth_ = conf_.getParameter<double>("jurassicWidth");
  
   // options
   useIsolEt_ = conf_.getParameter<bool>("useIsolEt");
   tryBoth_ = conf_.getParameter<bool>("tryBoth");
   subtract_ = conf_.getParameter<bool>("subtract");
   useNumCrystals_ = conf_.getParameter<bool>("useNumCrystals");
   vetoClustered_ = conf_.getParameter<bool>("vetoClustered");
  
   //EventSetup Tokens
   sevLvToken_ = esConsumes();
   caloGeometrytoken_ = esConsumes();
  
   //register your products
   produces<edm::ValueMap<double>>();
 }

References caloGeometrytoken_, conf_, ecalBarrelRecHitCollection_, ecalBarrelRecHitProducer_, ecalEndcapRecHitCollection_, ecalEndcapRecHitProducer_, egIsoConeSizeInBarrel_, egIsoConeSizeInEndcap_, egIsoConeSizeOut_, egIsoEMinBarrel_, egIsoEMinEndcap_, egIsoJurassicWidth_, egIsoPtMinBarrel_, egIsoPtMinEndcap_, emObjectProducer_, edm::ParameterSet::getParameter(), sevLvToken_, subtract_, tryBoth_, useIsolEt_, useNumCrystals_, and vetoClustered_.

◆ ~EgammaEcalRecHitIsolationProducer()

EgammaEcalRecHitIsolationProducer::~EgammaEcalRecHitIsolationProducer ( )

override

Definition at line 58 of file EgammaEcalRecHitIsolationProducer.cc.

58 {}

Member Function Documentation

◆ produce()

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

override

Definition at line 65 of file EgammaEcalRecHitIsolationProducer.cc.

                                                                                              {
   // Get the  filtered objects
   edm::Handle<edm::View<reco::Candidate>> emObjectHandle;
   iEvent.getByLabel(emObjectProducer_, emObjectHandle);
  
   // Next get Ecal hits barrel
   edm::Handle<EcalRecHitCollection> ecalBarrelRecHitHandle;  //EcalRecHitCollection is a typedef to
   iEvent.getByLabel(ecalBarrelRecHitProducer_.label(), ecalBarrelRecHitCollection_.label(), ecalBarrelRecHitHandle);
  
   // Next get Ecal hits endcap
   edm::Handle<EcalRecHitCollection> ecalEndcapRecHitHandle;
   iEvent.getByLabel(ecalEndcapRecHitProducer_.label(), ecalEndcapRecHitCollection_.label(), ecalEndcapRecHitHandle);
  
   edm::ESHandle<EcalSeverityLevelAlgo> sevlv = iSetup.getHandle(sevLvToken_);
   const EcalSeverityLevelAlgo* sevLevel = sevlv.product();
  
   //Get Calo Geometry
   edm::ESHandle<CaloGeometry> pG = iSetup.getHandle(caloGeometrytoken_);
   const CaloGeometry* caloGeom = pG.product();
  
   //reco::CandViewDoubleAssociations* isoMap = new reco::CandViewDoubleAssociations( reco::CandidateBaseRefProd( emObjectHandle ) );
   auto isoMap = std::make_unique<edm::ValueMap<double>>();
   edm::ValueMap<double>::Filler filler(*isoMap);
   std::vector<double> retV(emObjectHandle->size(), 0);
  
   EgammaRecHitIsolation ecalBarrelIsol(egIsoConeSizeOut_,
                                        egIsoConeSizeInBarrel_,
                                        egIsoJurassicWidth_,
                                        egIsoPtMinBarrel_,
                                        egIsoEMinBarrel_,
                                        caloGeom,
                                        *ecalBarrelRecHitHandle,
                                        sevLevel,
                                        DetId::Ecal);
   ecalBarrelIsol.setUseNumCrystals(useNumCrystals_);
   ecalBarrelIsol.setVetoClustered(vetoClustered_);
  
   EgammaRecHitIsolation ecalEndcapIsol(egIsoConeSizeOut_,
                                        egIsoConeSizeInEndcap_,
                                        egIsoJurassicWidth_,
                                        egIsoPtMinEndcap_,
                                        egIsoEMinEndcap_,
                                        caloGeom,
                                        *ecalEndcapRecHitHandle,
                                        sevLevel,
                                        DetId::Ecal);
   ecalEndcapIsol.setUseNumCrystals(useNumCrystals_);
   ecalEndcapIsol.setVetoClustered(vetoClustered_);
  
   for (size_t i = 0; i < emObjectHandle->size(); ++i) {
     //i need to know if its in the barrel/endcap so I get the supercluster handle to find out the detector eta
     //this might not be the best way, are we guaranteed that eta<1.5 is barrel
     //this can be safely replaced by another method which determines where the emobject is
     //then we either get the isolation Et or isolation Energy depending on user selection
     double isoValue = 0.;
  
     reco::SuperClusterRef superClus = emObjectHandle->at(i).get<reco::SuperClusterRef>();
  
     if (tryBoth_) {  //barrel + endcap
       if (useIsolEt_)
         isoValue =
             ecalBarrelIsol.getEtSum(&(emObjectHandle->at(i))) + ecalEndcapIsol.getEtSum(&(emObjectHandle->at(i)));
       else
         isoValue = ecalBarrelIsol.getEnergySum(&(emObjectHandle->at(i))) +
                    ecalEndcapIsol.getEnergySum(&(emObjectHandle->at(i)));
     } else if (fabs(superClus->eta()) < 1.479) {  //barrel
       if (useIsolEt_)
         isoValue = ecalBarrelIsol.getEtSum(&(emObjectHandle->at(i)));
       else
         isoValue = ecalBarrelIsol.getEnergySum(&(emObjectHandle->at(i)));
     } else {  //endcap
       if (useIsolEt_)
         isoValue = ecalEndcapIsol.getEtSum(&(emObjectHandle->at(i)));
       else
         isoValue = ecalEndcapIsol.getEnergySum(&(emObjectHandle->at(i)));
     }
  
     //we subtract off the electron energy here as well
     double subtractVal = 0;
  
     if (useIsolEt_)
       subtractVal = superClus.get()->rawEnergy() * sin(2 * atan(exp(-superClus.get()->eta())));
     else
       subtractVal = superClus.get()->rawEnergy();
  
     if (subtract_)
       isoValue -= subtractVal;
  
     retV[i] = isoValue;
     //all done, isolation is now in the map
  
   }  //end of loop over em objects
  
   filler.insert(emObjectHandle, retV.begin(), retV.end());
   filler.fill();
  
   iEvent.put(std::move(isoMap));
 }