---

EgammaHLTEcalRecIsolationProducer Class Reference

#include <RecoEgamma/EgammaHLTProducers/interface/EgammaHLTEcalRecIsolationProducer.h>

Inheritance diagram for EgammaHLTEcalRecIsolationProducer:

List of all members.

Public Member Functions
	EgammaHLTEcalRecIsolationProducer (const edm::ParameterSet &)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~EgammaHLTEcalRecIsolationProducer ()
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	recoEcalCandidateProducer_
bool	subtract_
bool	tryBoth_
bool	useIsolEt_

---

Detailed Description

Definition at line 33 of file EgammaHLTEcalRecIsolationProducer.h.

---

Constructor & Destructor Documentation

EgammaHLTEcalRecIsolationProducer::EgammaHLTEcalRecIsolationProducer

(

const edm::ParameterSet &

config

)

[explicit]

Definition at line 40 of file EgammaHLTEcalRecIsolationProducer.cc.

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

                                                                                                  : conf_(config)
{
  // use configuration file to setup input/output collection names
  //inputs
  recoEcalCandidateProducer_    = conf_.getParameter<edm::InputTag>("recoEcalCandidateProducer");
  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");

  //register your products
  produces < reco::RecoEcalCandidateIsolationMap >();
  
}

EgammaHLTEcalRecIsolationProducer::~EgammaHLTEcalRecIsolationProducer

(

)

Definition at line 71 of file EgammaHLTEcalRecIsolationProducer.cc.

{}

---

Member Function Documentation

void EgammaHLTEcalRecIsolationProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)			[virtual]

Implements edm::EDProducer.

Definition at line 75 of file EgammaHLTEcalRecIsolationProducer.cc.

References begin, DetId::Ecal, ecalBarrelRecHitCollection_, ecalBarrelRecHitProducer_, ecalEndcapRecHitCollection_, ecalEndcapRecHitProducer_, egIsoConeSizeInBarrel_, egIsoConeSizeInEndcap_, egIsoConeSizeOut_, egIsoJurassicWidth_, egIsoPtMinBarrel_, egIsoPtMinEndcap_, funct::exp(), edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByLabel(), EgammaRecHitIsolation::getEnergySum(), EgammaRecHitIsolation::getEtSum(), edm::AssociationMap< Tag >::insert(), edm::InputTag::label(), edm::ESHandle< T >::product(), edm::Event::put(), recoEcalCandidateProducer_, funct::sin(), subtract_, tryBoth_, and useIsolEt_.

                                                                                            {

  // Get the RecoEcalCandidate Collection
  edm::Handle<reco::RecoEcalCandidateCollection> recoecalcandHandle;
  iEvent.getByLabel(recoEcalCandidateProducer_,recoecalcandHandle);

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

  //create the meta hit collections inorder that we can pass them into the isolation objects

  EcalRecHitMetaCollection ecalBarrelHits(*ecalBarrelRecHitHandle);
  EcalRecHitMetaCollection ecalEndcapHits(*ecalEndcapRecHitHandle);

  //Get Calo Geometry
  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  const CaloGeometry* caloGeom = pG.product();

  //prepare product
  reco::RecoEcalCandidateIsolationMap isoMap;

  //create algorithm objects
  EgammaRecHitIsolation ecalBarrelIsol(egIsoConeSizeOut_,egIsoConeSizeInBarrel_,egIsoJurassicWidth_,egIsoPtMinBarrel_,caloGeom,&ecalBarrelHits,DetId::Ecal);
  EgammaRecHitIsolation ecalEndcapIsol(egIsoConeSizeOut_,egIsoConeSizeInEndcap_,egIsoJurassicWidth_,egIsoPtMinEndcap_,caloGeom,&ecalEndcapHits,DetId::Ecal);


  for (reco::RecoEcalCandidateCollection::const_iterator iRecoEcalCand= recoecalcandHandle->begin(); iRecoEcalCand!=recoecalcandHandle->end(); iRecoEcalCand++) {
    
    //create reference for storage in isolation map
    reco::RecoEcalCandidateRef recoecalcandref(reco::RecoEcalCandidateRef(recoecalcandHandle,iRecoEcalCand -recoecalcandHandle ->begin()));
    
    //ecal isolation is centered on supecluster
    reco::SuperClusterRef superClus = iRecoEcalCand->get<reco::SuperClusterRef>();

    //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
    float isol =0.;

    if(tryBoth_){ //barrel + endcap
      if(useIsolEt_) isol =  ecalBarrelIsol.getEtSum(&(*iRecoEcalCand)) + ecalEndcapIsol.getEtSum(&(*iRecoEcalCand));
      else           isol =  ecalBarrelIsol.getEnergySum(&(*iRecoEcalCand)) + ecalEndcapIsol.getEnergySum(&(*iRecoEcalCand));
    }
    else if( fabs(superClus->eta())<1.479) { //barrel
      if(useIsolEt_) isol =  ecalBarrelIsol.getEtSum(&(*iRecoEcalCand));
      else           isol =  ecalBarrelIsol.getEnergySum(&(*iRecoEcalCand));
    }
    else{ //endcap
      if(useIsolEt_) isol =  ecalEndcapIsol.getEtSum(&(*iRecoEcalCand));
      else           isol =  ecalEndcapIsol.getEnergySum(&(*iRecoEcalCand));
    }

    //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_) isol-= subtractVal;



    isoMap.insert(recoecalcandref, isol);

  }

  std::auto_ptr<reco::RecoEcalCandidateIsolationMap> isolMap(new reco::RecoEcalCandidateIsolationMap(isoMap));
  iEvent.put(isolMap);

}

---

Member Data Documentation

edm::ParameterSet EgammaHLTEcalRecIsolationProducer::conf_ [private]

Definition at line 62 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer().

edm::InputTag EgammaHLTEcalRecIsolationProducer::ecalBarrelRecHitCollection_ [private]

Definition at line 45 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

edm::InputTag EgammaHLTEcalRecIsolationProducer::ecalBarrelRecHitProducer_ [private]

Definition at line 44 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

edm::InputTag EgammaHLTEcalRecIsolationProducer::ecalEndcapRecHitCollection_ [private]

Definition at line 47 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

edm::InputTag EgammaHLTEcalRecIsolationProducer::ecalEndcapRecHitProducer_ [private]

Definition at line 46 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

double EgammaHLTEcalRecIsolationProducer::egIsoConeSizeInBarrel_ [private]

Definition at line 54 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

double EgammaHLTEcalRecIsolationProducer::egIsoConeSizeInEndcap_ [private]

Definition at line 55 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

double EgammaHLTEcalRecIsolationProducer::egIsoConeSizeOut_ [private]

Definition at line 53 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

double EgammaHLTEcalRecIsolationProducer::egIsoEMinBarrel_ [private]

Definition at line 50 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer().

double EgammaHLTEcalRecIsolationProducer::egIsoEMinEndcap_ [private]

Definition at line 52 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer().

double EgammaHLTEcalRecIsolationProducer::egIsoJurassicWidth_ [private]

Definition at line 56 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

double EgammaHLTEcalRecIsolationProducer::egIsoPtMinBarrel_ [private]

Definition at line 49 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

double EgammaHLTEcalRecIsolationProducer::egIsoPtMinEndcap_ [private]

Definition at line 51 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

edm::InputTag EgammaHLTEcalRecIsolationProducer::recoEcalCandidateProducer_ [private]

Definition at line 43 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

bool EgammaHLTEcalRecIsolationProducer::subtract_ [private]

Definition at line 60 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

bool EgammaHLTEcalRecIsolationProducer::tryBoth_ [private]

Definition at line 59 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

bool EgammaHLTEcalRecIsolationProducer::useIsolEt_ [private]

Definition at line 58 of file EgammaHLTEcalRecIsolationProducer.h.

Referenced by EgammaHLTEcalRecIsolationProducer(), and produce().

---

The documentation for this class was generated from the following files:

• RecoEgamma/EgammaHLTProducers/interface/EgammaHLTEcalRecIsolationProducer.h
• RecoEgamma/EgammaHLTProducers/src/EgammaHLTEcalRecIsolationProducer.cc

---