#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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
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_

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

Constructor & Destructor Documentation

◆ EgammaEcalRecHitIsolationProducer()

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

explicit

Definition at line 30 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");
  
   //register your products
   produces<edm::ValueMap<double>>();
 }

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

◆ ~EgammaEcalRecHitIsolationProducer()

EgammaEcalRecHitIsolationProducer::~EgammaEcalRecHitIsolationProducer ( )

override

Definition at line 60 of file EgammaEcalRecHitIsolationProducer.cc.

60 {}

Member Function Documentation

◆ produce()

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

override

Definition at line 67 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.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
   const EcalSeverityLevelAlgo* sevLevel = sevlv.product();
  
   //Get Calo Geometry
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   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));
 }