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 Member Functions
bool	isNumber (const std::string &str) const

double	toNumber (const std::string &str) const

Private Attributes
reco::isodeposit::AbsVetos	barrelVetos_

double	deltaR_

reco::isodeposit::AbsVetos	endcapVetos_

reco::isodeposit::EventDependentAbsVetos	evdepVetos_

edm::Handle< reco::IsoDepositMap >	hDeps_

Mode	mode_

bool	skipDefaultVeto_

edm::EDGetTokenT< reco::IsoDepositMap >	srcToken_

bool	usePivotForBarrelEndcaps_

bool	usesFunction_

double	weight_

StringObjectFunction< reco::Candidate >	weightExpr_

Detailed Description

Definition at line 43 of file PFCandIsolatorFromDeposit.cc.

Constructor & Destructor Documentation

◆ SingleDeposit()

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

Definition at line 83 of file PFCandIsolatorFromDeposit.cc.

     : 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")),
       usePivotForBarrelEndcaps_(iConfig.getParameter<bool>("PivotCoordinatesForEBEE"))
 //,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
     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 = nullptr;
   static const std::regex ecalSwitch("^Ecal(Barrel|Endcaps):(.*)");
 
   for (vstring::const_iterator it = vetos.begin(), ed = vetos.end(); it != ed; ++it) {
     std::cmatch match;
     // in that case, make two series of vetoes
     if (usePivotForBarrelEndcaps_) {
       if (regex_match(it->c_str(), match, ecalSwitch)) {
         if (match[1] == "Barrel") {
           //        std::cout << " Adding Barrel veto " << std::string(match[2]) << std::endl;
           barrelVetos_.push_back(
               IsoDepositVetoFactory::make(std::string(match[2]).c_str(), evdep, iC));  // I don't know a better syntax
         }
         if (match[1] == "Endcaps") {
           //        std::cout << " Adding Endcap veto " << std::string(match[2]) << std::endl;
           endcapVetos_.push_back(IsoDepositVetoFactory::make(std::string(match[2]).c_str(), evdep, iC));
         }
       } else {
         barrelVetos_.push_back(IsoDepositVetoFactory::make(it->c_str(), evdep, iC));
         endcapVetos_.push_back(IsoDepositVetoFactory::make(it->c_str(), evdep, iC));
       }
     } else {
       //only one serie of vetoes, just barrel
       barrelVetos_.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 << "PFCandIsolatorFromDeposits::SingleDeposit::SingleDeposit: Total of " << vetos_.size() << " vetos" << std::endl;
 }

Member Function Documentation

◆ cleanup()

void PFCandIsolatorFromDeposits::SingleDeposit::cleanup ( )

Definition at line 155 of file PFCandIsolatorFromDeposit.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

Referenced by esMonitoring.FDJsonServer::handle_close().

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

◆ compute()

double PFCandIsolatorFromDeposits::SingleDeposit::compute ( const reco::CandidateBaseRef & cand )

Definition at line 174 of file PFCandIsolatorFromDeposit.cc.

                                                                                         {
   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
   bool barrel = true;
   if (usePivotForBarrelEndcaps_) {
     const reco::PFCandidate *myPFCand = dynamic_cast<const reco::PFCandidate *>(&(*cand));
     if (myPFCand) {
       // exact barrel boundary
       barrel = fabs(myPFCand->positionAtECALEntrance().eta()) < 1.479;
     } else {
       const reco::RecoCandidate *myRecoCand = dynamic_cast<const reco::RecoCandidate *>(&(*cand));
       if (myRecoCand) {
         // not optimal. isEB should be used.
         barrel = (fabs(myRecoCand->superCluster()->eta()) < 1.479);
       }
     }
   }
   // if ! usePivotForBarrelEndcaps_ only the barrel series is used, which does not prevent the vetoes do be different in barrel & endcaps
   reco::isodeposit::AbsVetos *vetos = (barrel) ? &barrelVetos_ : &endcapVetos_;
 
   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();
   }
   throw cms::Exception("Logic error") << "Should not happen at " << __FILE__ << ", line "
                                       << __LINE__;  // avoid gcc warning
 }

◆ isNumber()

bool PFCandIsolatorFromDeposits::SingleDeposit::isNumber ( const std::string & str ) const

private

Definition at line 77 of file PFCandIsolatorFromDeposit.cc.

References str.

Referenced by SingleDeposit().

                                                                                  {
   static const std::regex re("^[+-]?(\\d+\\.?|\\d*\\.\\d*)$");
   return regex_match(str.c_str(), re);
 }

◆ map()

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

inline

Definition at line 49 of file PFCandIsolatorFromDeposit.cc.

References hDeps_.

49 { return *hDeps_; }

PFCandIsolatorFromDeposits::SingleDeposit::hDeps_

edm::Handle< reco::IsoDepositMap > hDeps_

Definition: PFCandIsolatorFromDeposit.cc:64

◆ open()

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

Definition at line 167 of file PFCandIsolatorFromDeposit.cc.

References iEvent, and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

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