CandIsolatorFromDeposits::SingleDeposit Class Reference

Public Member Functions
void	cleanup ()
double	compute (const reco::CandidateBaseRef &cand)
const reco::IsoDepositMap &	map ()
void	open (const edm::Event &iEvent, const edm::EventSetup &iSetup)
	SingleDeposit (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

Private Attributes
double	deltaR_
reco::isodeposit::EventDependentAbsVetos	evdepVetos_
edm::Handle< reco::IsoDepositMap >	hDeps_
Mode	mode_
bool	skipDefaultVeto_
edm::EDGetTokenT< reco::IsoDepositMap >	srcToken_
bool	usesFunction_
reco::isodeposit::AbsVetos	vetos_
double	weight_
StringObjectFunction< reco::Candidate >	weightExpr_

Detailed Description

Definition at line 39 of file CandIsolatorFromDeposits.h.

Constructor & Destructor Documentation

CandIsolatorFromDeposits::SingleDeposit::SingleDeposit	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 36 of file CandIsolatorFromDeposits.cc.

References CandIsolatorFromDeposits::Count, evdepVetos_, Exception, edm::ParameterSet::getParameter(), isNumber(), IsoDepositVetoFactory::make(), CandIsolatorFromDeposits::Max, CandIsolatorFromDeposits::MaxRelative, CandIsolatorFromDeposits::MeanDR, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, mode_, CandIsolatorFromDeposits::NearestDR, AlCaHLTBitMon_QueryRunRegistry::string, CandIsolatorFromDeposits::Sum, CandIsolatorFromDeposits::Sum2, CandIsolatorFromDeposits::Sum2Relative, CandIsolatorFromDeposits::SumDR, CandIsolatorFromDeposits::SumRelative, toNumber(), usesFunction_, boostedElectronIsolation_cff::vetos, vetos_, and weight_.

Referenced by CandIsolatorFromDeposits::CandIsolatorFromDeposits().

    : srcToken_(iC.consumes<reco::IsoDepositMap>(iConfig.getParameter<edm::InputTag>("src"))),
      deltaR_(iConfig.getParameter<double>("deltaR")),
      weightExpr_(iConfig.getParameter<std::string>("weight")),
      skipDefaultVeto_(iConfig.getParameter<bool>("skipDefaultVeto"))
//,vetos_(new AbsVetos())
{
  std::string mode = iConfig.getParameter<std::string>("mode");
  if (mode == "sum")
    mode_ = Sum;
  else if (mode == "sumRelative")
    mode_ = SumRelative;
  else if (mode == "sum2")
    mode_ = Sum2;
  else if (mode == "sum2Relative")
    mode_ = Sum2Relative;
  else if (mode == "max")
    mode_ = Max;
  else if (mode == "maxRelative")
    mode_ = MaxRelative;
  else if (mode == "nearestDR")
    mode_ = NearestDR;
  else if (mode == "count")
    mode_ = Count;
  else if (mode == "meanDR")
    mode_ = MeanDR;
  else if (mode == "sumDR")
    mode_ = SumDR;
  else
    throw cms::Exception("Not Implemented") << "Mode '" << mode << "' not implemented. "
                                            << "Supported modes are 'sum', 'sumRelative', 'count'." <<
        //"Supported modes are 'sum', 'sumRelative', 'max', 'maxRelative', 'count'." << // TODO: on request only
        "New methods can be easily implemented if requested.";
  typedef std::vector<std::string> vstring;
  vstring vetos = iConfig.getParameter<vstring>("vetos");
  reco::isodeposit::EventDependentAbsVeto *evdep;
  for (vstring::const_iterator it = vetos.begin(), ed = vetos.end(); it != ed; ++it) {
    vetos_.push_back(IsoDepositVetoFactory::make(it->c_str(), evdep, iC));
    if (evdep)
      evdepVetos_.push_back(evdep);
  }
  std::string weight = iConfig.getParameter<std::string>("weight");
  if (isNumber(weight)) {
    //std::cout << "Weight is a simple number, " << toNumber(weight) << std::endl;
    weight_ = toNumber(weight);
    usesFunction_ = false;
  } else {
    usesFunction_ = true;
    //std::cout << "Weight is a function, this might slow you down... " << std::endl;
  }
  //std::cout << "CandIsolatorFromDeposits::SingleDeposit::SingleDeposit: Total of " << vetos_.size() << " vetos" << std::endl;
}

Member Function Documentation

void CandIsolatorFromDeposits::SingleDeposit::cleanup

(

)

Definition at line 88 of file CandIsolatorFromDeposits.cc.

References evdepVetos_, and vetos_.

                                                    {
  for (AbsVetos::iterator it = vetos_.begin(), ed = vetos_.end(); it != ed; ++it) {
    delete *it;
  }
  vetos_.clear();
  // NOTE: we DON'T have to delete the evdepVetos_, they have already been deleted above. We just clear the vectors
  evdepVetos_.clear();
}

double CandIsolatorFromDeposits::SingleDeposit::compute

(

const reco::CandidateBaseRef &

cand

)

Definition at line 103 of file CandIsolatorFromDeposits.cc.

References reco::IsoDeposit::algoWithin(), reco::IsoDeposit::candEnergy(), CandIsolatorFromDeposits::Count, reco::IsoDeposit::countWithin(), deltaR_, PVValHelper::eta, reco::IsoDeposit::eta(), Exception, CandIsolatorFromDeposits::Max, CandIsolatorFromDeposits::MaxRelative, reco::IsoDeposit::maxWithin(), CandIsolatorFromDeposits::MeanDR, mode_, CandIsolatorFromDeposits::NearestDR, reco::IsoDeposit::nearestDR(), phi, reco::IsoDeposit::phi(), skipDefaultVeto_, CandIsolatorFromDeposits::Sum, CandIsolatorFromDeposits::Sum2, CandIsolatorFromDeposits::Sum2Relative, reco::IsoDeposit::sum2Within(), CandIsolatorFromDeposits::SumDR, CandIsolatorFromDeposits::SumRelative, reco::IsoDeposit::sumWithin(), usesFunction_, vetos_, weight_, and weightExpr_.

                                                                                      {
  const IsoDeposit &dep = (*hDeps_)[cand];
  double eta = dep.eta(), phi = dep.phi();  // better to center on the deposit direction
                                            // that could be, e.g., the impact point at calo
  for (AbsVetos::iterator it = vetos_.begin(), ed = vetos_.end(); it != ed; ++it) {
    (*it)->centerOn(eta, phi);
  }
  double weight = (usesFunction_ ? weightExpr_(*cand) : weight_);
  switch (mode_) {
    case Count:
      return weight * dep.countWithin(deltaR_, vetos_, skipDefaultVeto_);
    case Sum:
      return weight * dep.sumWithin(deltaR_, vetos_, skipDefaultVeto_);
    case SumRelative:
      return weight * dep.sumWithin(deltaR_, vetos_, skipDefaultVeto_) / dep.candEnergy();
    case Sum2:
      return weight * dep.sum2Within(deltaR_, vetos_, skipDefaultVeto_);
    case Sum2Relative:
      return weight * dep.sum2Within(deltaR_, vetos_, skipDefaultVeto_) / (dep.candEnergy() * dep.candEnergy());
    case Max:
      return weight * dep.maxWithin(deltaR_, vetos_, skipDefaultVeto_);
    case NearestDR:
      return weight * dep.nearestDR(deltaR_, vetos_, skipDefaultVeto_);
    case MaxRelative:
      return weight * dep.maxWithin(deltaR_, vetos_, skipDefaultVeto_) / dep.candEnergy();
    case MeanDR:
      return weight * dep.algoWithin<reco::IsoDeposit::MeanDRAlgo>(deltaR_, vetos_, skipDefaultVeto_);
    case SumDR:
      return weight * dep.algoWithin<reco::IsoDeposit::SumDRAlgo>(deltaR_, vetos_, skipDefaultVeto_);
  }
  throw cms::Exception("Logic error") << "Should not happen at " << __FILE__ << ", line "
                                      << __LINE__;  // avoid gcc warning
}

const reco::IsoDepositMap& CandIsolatorFromDeposits::SingleDeposit::map ( )

inline

Definition at line 45 of file CandIsolatorFromDeposits.h.

Referenced by CandIsolatorFromDeposits::produce().

{ return *hDeps_; }

void CandIsolatorFromDeposits::SingleDeposit::open	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 96 of file CandIsolatorFromDeposits.cc.

References evdepVetos_, edm::Event::getByToken(), hDeps_, and srcToken_.

                                                                                                    {
  iEvent.getByToken(srcToken_, hDeps_);
  for (EventDependentAbsVetos::iterator it = evdepVetos_.begin(), ed = evdepVetos_.end(); it != ed; ++it) {
    (*it)->setEvent(iEvent, iSetup);
  }
}

Member Data Documentation

double CandIsolatorFromDeposits::SingleDeposit::deltaR_

private

Definition at line 50 of file CandIsolatorFromDeposits.h.

Referenced by compute().

reco::isodeposit::EventDependentAbsVetos CandIsolatorFromDeposits::SingleDeposit::evdepVetos_

private

Definition at line 55 of file CandIsolatorFromDeposits.h.

Referenced by cleanup(), open(), and SingleDeposit().

edm::Handle<reco::IsoDepositMap> CandIsolatorFromDeposits::SingleDeposit::hDeps_

private

Definition at line 57 of file CandIsolatorFromDeposits.h.

Referenced by open().

Mode CandIsolatorFromDeposits::SingleDeposit::mode_

private

Definition at line 48 of file CandIsolatorFromDeposits.h.

Referenced by compute(), and SingleDeposit().

bool CandIsolatorFromDeposits::SingleDeposit::skipDefaultVeto_

private

Definition at line 56 of file CandIsolatorFromDeposits.h.

Referenced by compute().

edm::EDGetTokenT<reco::IsoDepositMap> CandIsolatorFromDeposits::SingleDeposit::srcToken_

private

Definition at line 49 of file CandIsolatorFromDeposits.h.

Referenced by open().

bool CandIsolatorFromDeposits::SingleDeposit::usesFunction_

private

Definition at line 51 of file CandIsolatorFromDeposits.h.

Referenced by compute(), and SingleDeposit().

reco::isodeposit::AbsVetos CandIsolatorFromDeposits::SingleDeposit::vetos_

private

Definition at line 54 of file CandIsolatorFromDeposits.h.

Referenced by cleanup(), compute(), and SingleDeposit().

double CandIsolatorFromDeposits::SingleDeposit::weight_

private

Definition at line 52 of file CandIsolatorFromDeposits.h.

Referenced by compute(), and SingleDeposit().

StringObjectFunction<reco::Candidate> CandIsolatorFromDeposits::SingleDeposit::weightExpr_

private

Definition at line 53 of file CandIsolatorFromDeposits.h.

Referenced by compute().