#include <EgammaEcalExtractor.h>

Inheritance diagram for egammaisolation::EgammaEcalExtractor:

Public Member Functions
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &c) const override

reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const override

	EgammaEcalExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)

	EgammaEcalExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &iC)

void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks) override

	~EgammaEcalExtractor () override

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
	make single IsoDeposit based on a CandidateBaseRef as input More...

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const

virtual void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Attributes
edm::EDGetTokenT< reco::BasicClusterCollection >	basicClusterToken_

double	conesize_

double	etMin_

bool	scmatch_

edm::EDGetTokenT< reco::SuperClusterCollection >	superClusterToken_

Detailed Description

Definition at line 39 of file EgammaEcalExtractor.h.

Constructor & Destructor Documentation

◆ EgammaEcalExtractor() [1/2]

egammaisolation::EgammaEcalExtractor::EgammaEcalExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 41 of file EgammaEcalExtractor.h.

41 : EgammaEcalExtractor(par, iC) {}

◆ EgammaEcalExtractor() [2/2]

egammaisolation::EgammaEcalExtractor::EgammaEcalExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &	iC
	)

inline

Definition at line 42 of file EgammaEcalExtractor.h.

         : etMin_(par.getParameter<double>("etMin")),
           conesize_(par.getParameter<double>("extRadius")),
           scmatch_(par.getParameter<bool>("superClusterMatch")),
           basicClusterToken_(
               iC.consumes<reco::BasicClusterCollection>(par.getParameter<edm::InputTag>("basicClusters"))),
           superClusterToken_(
               iC.consumes<reco::SuperClusterCollection>(par.getParameter<edm::InputTag>("superClusters"))) {}

◆ ~EgammaEcalExtractor()

EgammaEcalExtractor::~EgammaEcalExtractor ( )

override

Definition at line 23 of file EgammaEcalExtractor.cc.

23 {}

Member Function Documentation

◆ deposit() [1/2]

reco::IsoDeposit EgammaEcalExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Candidate &	track
	)		const

overridevirtual

make single IsoDeposit based on a candidate as input purely virtual: have to implement in concrete implementations

Reimplemented from reco::isodeposit::IsoDepositExtractor.

Definition at line 25 of file EgammaEcalExtractor.cc.

                                                                                     {
   edm::Handle<reco::SuperClusterCollection> superClusterCollectionH;
   edm::Handle<reco::BasicClusterCollection> basicClusterCollectionH;
   ev.getByToken(superClusterToken_, superClusterCollectionH);
   ev.getByToken(basicClusterToken_, basicClusterCollectionH);
  
   reco::SuperClusterRef sc = candidate.get<reco::SuperClusterRef>();
   math::XYZPoint position = sc->position();
   // match the photon hybrid supercluster with those with Algo==0 (island)
   double delta1 = 1000.;
   double deltacur = 1000.;
   const reco::SuperCluster *matchedsupercluster = nullptr;
   bool MATCHEDSC = false;
  
   Direction candDir(position.eta(), position.phi());
   reco::IsoDeposit deposit(candDir);
   deposit.setVeto(reco::IsoDeposit::Veto(candDir, 0));  // no veto is needed for this deposit
   deposit.addCandEnergy(sc->energy() * sin(2 * atan(exp(-sc->eta()))));
  
   for (reco::SuperClusterCollection::const_iterator scItr = superClusterCollectionH->begin();
        scItr != superClusterCollectionH->end();
        ++scItr) {
     const reco::SuperCluster *supercluster = &(*scItr);
  
     if (supercluster->seed()->algo() == 0) {
       deltacur = ROOT::Math::VectorUtil::DeltaR(supercluster->position(), position);
       if (deltacur < delta1) {
         delta1 = deltacur;
         matchedsupercluster = supercluster;
         MATCHEDSC = true;
       }
     }
   }
  
   const reco::BasicCluster *cluster = nullptr;
  
   //loop over basic clusters
   for (reco::BasicClusterCollection::const_iterator cItr = basicClusterCollectionH->begin();
        cItr != basicClusterCollectionH->end();
        ++cItr) {
     cluster = &(*cItr);
     //    double ebc_bcchi2 = cluster->chi2();
     int ebc_bcalgo = cluster->algo();
     double ebc_bce = cluster->energy();
     double ebc_bceta = cluster->eta();
     double ebc_bcet = ebc_bce * sin(2 * atan(exp(ebc_bceta)));
     double newDelta = 0.;
  
     if (ebc_bcet > etMin_ && ebc_bcalgo == 0) {
       //      if (ebc_bcchi2 < 30.) {
  
       if (MATCHEDSC || !scmatch_) {  //skip selection if user wants to fill all superclusters
         bool inSuperCluster = false;
  
         if (scmatch_) {  // only try the matching if needed
           reco::CaloCluster_iterator theEclust = matchedsupercluster->clustersBegin();
           // loop over the basic clusters of the matched supercluster
           for (; theEclust != matchedsupercluster->clustersEnd(); ++theEclust) {
             if ((**theEclust) == (*cluster))
               inSuperCluster = true;
           }
         }
         if (!inSuperCluster || !scmatch_) {  //skip selection if user wants to fill all superclusters
           newDelta = ROOT::Math::VectorUtil::DeltaR(cluster->position(), position);
           if (newDelta < conesize_) {
             deposit.addDeposit(Direction(cluster->eta(), cluster->phi()), ebc_bcet);
           }
         }
       }
       //      } // matches ebc_bcchi2
     }  // matches ebc_bcet && ebc_bcalgo
   }
  
   //  std::cout << "Will return ecalIsol = " << ecalIsol << std::endl;
   return deposit;
 }

References reco::CaloCluster::algo(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), electronAnalyzer_cfi::DeltaR, reco::CaloCluster::energy(), reco::CaloCluster::eta(), ev, JetChargeProducer_cfi::exp, reco::Candidate::get(), reco::CaloCluster::phi(), reco::CaloCluster::position(), position, reco::SuperCluster::seed(), and funct::sin().

◆ deposit() [2/2]

reco::IsoDeposit egammaisolation::EgammaEcalExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

inlineoverridevirtual

make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 55 of file EgammaEcalExtractor.h.

                                                                     {
       throw cms::Exception("Configuration Error")
           << "This extractor " << (typeid(this).name()) << " is not made for tracks";
     }

References Exception, and Skims_PA_cff::name.

◆ fillVetos()

void egammaisolation::EgammaEcalExtractor::fillVetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::TrackCollection &	tracks
	)

inlineoverridevirtual

fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 53 of file EgammaEcalExtractor.h.

53 {

54 }

Member Data Documentation

◆ basicClusterToken_

edm::EDGetTokenT<reco::BasicClusterCollection> egammaisolation::EgammaEcalExtractor::basicClusterToken_

private

Definition at line 74 of file EgammaEcalExtractor.h.

◆ conesize_

double egammaisolation::EgammaEcalExtractor::conesize_

private

Definition at line 71 of file EgammaEcalExtractor.h.

◆ etMin_

double egammaisolation::EgammaEcalExtractor::etMin_

private

Definition at line 70 of file EgammaEcalExtractor.h.

◆ scmatch_

bool egammaisolation::EgammaEcalExtractor::scmatch_

private

Definition at line 72 of file EgammaEcalExtractor.h.

◆ superClusterToken_

edm::EDGetTokenT<reco::SuperClusterCollection> egammaisolation::EgammaEcalExtractor::superClusterToken_

private

Definition at line 75 of file EgammaEcalExtractor.h.

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EgammaEcalExtractor() [1/2]

◆ EgammaEcalExtractor() [2/2]

◆ ~EgammaEcalExtractor()

Member Function Documentation

◆ deposit() [1/2]

◆ deposit() [2/2]

◆ fillVetos()

Member Data Documentation

◆ basicClusterToken_

◆ conesize_

◆ etMin_

◆ scmatch_

◆ superClusterToken_