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

bool	hasAbilityToProduceInRuns () 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 34 of file EgammaEcalRecHitIsolationProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 31 of file EgammaEcalRecHitIsolationProducer.cc.

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

                                                                                                   : 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> >();
 }

EgammaEcalRecHitIsolationProducer::~EgammaEcalRecHitIsolationProducer ( )

override

Definition at line 65 of file EgammaEcalRecHitIsolationProducer.cc.

65 {}

Member Function Documentation

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

override

Definition at line 74 of file EgammaEcalRecHitIsolationProducer.cc.

References DetId::Ecal, ecalBarrelRecHitCollection_, ecalBarrelRecHitProducer_, ecalEndcapRecHitCollection_, ecalEndcapRecHitProducer_, egIsoConeSizeInBarrel_, egIsoConeSizeInEndcap_, egIsoConeSizeOut_, egIsoEMinBarrel_, egIsoEMinEndcap_, egIsoJurassicWidth_, egIsoPtMinBarrel_, egIsoPtMinEndcap_, emObjectProducer_, JetChargeProducer_cfi::exp, objects.autophobj::filler, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), EgammaRecHitIsolation::getEnergySum(), EgammaRecHitIsolation::getEtSum(), mps_fire::i, edm::InputTag::label(), eostools::move(), edm::ESHandle< T >::product(), edm::Event::put(), EgammaRecHitIsolation::setUseNumCrystals(), EgammaRecHitIsolation::setVetoClustered(), funct::sin(), subtract_, tryBoth_, useIsolEt_, useNumCrystals_, and vetoClustered_.

 {
 
 
   // 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));
 
 }