egammaisolation::EgammaEcalExtractor Class Reference

Inheritance diagram for egammaisolation::EgammaEcalExtractor:

Public Member Functions
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const override
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &c) 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::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...
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 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 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 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 37 of file EgammaEcalExtractor.cc.

Constructor & Destructor Documentation

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

inline

Definition at line 39 of file EgammaEcalExtractor.cc.

: EgammaEcalExtractor(par, iC) {}

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

inline

Definition at line 40 of file EgammaEcalExtractor.cc.

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

override

Definition at line 85 of file EgammaEcalExtractor.cc.

{}

Member Function Documentation

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 53 of file EgammaEcalExtractor.cc.

References Exception, and mergeVDriftHistosByStation::name.

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

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 87 of file EgammaEcalExtractor.cc.

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

                                                                                    {
  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 (auto const& scItr : ev.get(superClusterToken_)) {
    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 (auto const& cItr : ev.get(basicClusterToken_)) {
    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;
}

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 51 of file EgammaEcalExtractor.cc.

                                                                                                                 {
    }

Member Data Documentation

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

private

Definition at line 72 of file EgammaEcalExtractor.cc.

double egammaisolation::EgammaEcalExtractor::conesize_

private

Definition at line 69 of file EgammaEcalExtractor.cc.

double egammaisolation::EgammaEcalExtractor::etMin_

private

Definition at line 68 of file EgammaEcalExtractor.cc.

bool egammaisolation::EgammaEcalExtractor::scmatch_

private

Definition at line 70 of file EgammaEcalExtractor.cc.

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

private

Definition at line 73 of file EgammaEcalExtractor.cc.