#include <IsolatorByNominalEfficiency.h>

Inheritance diagram for muonisolation::IsolatorByNominalEfficiency:

Classes
class	ConeSizes

Public Types
typedef MuIsoBaseIsolator::DepositContainer	DepositContainer

Public Types inherited from muonisolation::MuIsoBaseIsolator
typedef std::vector< DepositAndVetos >	DepositContainer

enum	ResultType { ISOL_INT_TYPE = 0, ISOL_FLOAT_TYPE, ISOL_BOOL_TYPE, ISOL_INVALID_TYPE }

typedef reco::IsoDeposit::Veto	Veto

typedef reco::IsoDeposit::Vetos	Vetos

Public Member Functions
Cuts	cuts (float nominalEfficiency) const

	IsolatorByNominalEfficiency (const std::string &thrFile, const std::vector< std::string > &ceff, const std::vector< double > &weights)
	Constructor. More...

	IsolatorByNominalEfficiency (const std::string &thrFile, const std::vector< std::string > &ceff, const std::vector< double > &weights, const std::vector< double > &thresh)

Result	result (const DepositContainer &deposits, const edm::Event *=nullptr) const override
	Compute the deposit within the cone and return the isolation result. More...

	~IsolatorByNominalEfficiency () override

Public Member Functions inherited from muonisolation::MuIsoBaseIsolator
virtual Result	result (const DepositContainer &deposits, const reco::Candidate &muon, const edm::Event *=nullptr) const
	Compute and return the isolation variable, with vetoes and the muon. More...

virtual Result	result (const DepositContainer &deposits, const reco::Track &muon, const edm::Event *=nullptr) const
	Compute and return the isolation variable, with vetoes and the muon. More...

virtual	~MuIsoBaseIsolator ()

Private Types
typedef std::multimap< float, int >	mapNomEff_Cone

Private Member Functions
int	bestConeForEfficiencyIndex (float eff_thr) const

mapNomEff_Cone	cones (const std::vector< std::string > &names)

std::string	findPath (const std::string &fileName)

ResultType	resultType () const override

virtual double	weightedSum (const DepositContainer &deposits, float dRcone) const

Private Attributes
mapNomEff_Cone	coneForEfficiency

ConeSizes	theConesInfo

std::vector< double >	theDepThresholds

std::vector< double >	theWeights

NominalEfficiencyThresholds *	thresholds

Detailed Description

Definition at line 22 of file IsolatorByNominalEfficiency.h.

Member Typedef Documentation

◆ DepositContainer

typedef MuIsoBaseIsolator::DepositContainer muonisolation::IsolatorByNominalEfficiency::DepositContainer

Definition at line 24 of file IsolatorByNominalEfficiency.h.

◆ mapNomEff_Cone

typedef std::multimap<float, int> muonisolation::IsolatorByNominalEfficiency::mapNomEff_Cone

private

Definition at line 61 of file IsolatorByNominalEfficiency.h.

Constructor & Destructor Documentation

◆ IsolatorByNominalEfficiency() [1/2]

muonisolation::IsolatorByNominalEfficiency::IsolatorByNominalEfficiency	(	const std::string &	thrFile,
		const std::vector< std::string > &	ceff,
		const std::vector< double > &	weights
	)

Constructor.

◆ IsolatorByNominalEfficiency() [2/2]

muonisolation::IsolatorByNominalEfficiency::IsolatorByNominalEfficiency	(	const std::string &	thrFile,
		const std::vector< std::string > &	ceff,
		const std::vector< double > &	weights,
		const std::vector< double > &	thresh
	)

◆ ~IsolatorByNominalEfficiency()

IsolatorByNominalEfficiency::~IsolatorByNominalEfficiency ( )

override

Definition at line 44 of file IsolatorByNominalEfficiency.cc.

44 { delete thresholds; }

References thresholds.

Member Function Documentation

◆ bestConeForEfficiencyIndex()

int IsolatorByNominalEfficiency::bestConeForEfficiencyIndex ( float eff_thr ) const

private

Definition at line 107 of file IsolatorByNominalEfficiency.cc.

                                                                                {
   //FIXME use upper_bound
   int best_cone;
   if (!coneForEfficiency.empty()) {
     best_cone = (--(coneForEfficiency.end()))->second;
   } else
     return 0;
  
   mapNomEff_Cone::const_reverse_iterator it;
   for (it = coneForEfficiency.rbegin(); it != coneForEfficiency.rend(); it++) {
     if (eff_thr <= (*it).first)
       best_cone = (*it).second;
   }
   return best_cone;
 }

References coneForEfficiency.

Referenced by cuts(), and result().

◆ cones()

IsolatorByNominalEfficiency::mapNomEff_Cone IsolatorByNominalEfficiency::cones ( const std::vector< std::string > & names )

private

Definition at line 46 of file IsolatorByNominalEfficiency.cc.

                                                                                                          {
   mapNomEff_Cone result;
   for (vector<string>::const_iterator is = usrVec.begin(); is != usrVec.end(); is++) {
     char* evp = nullptr;
     int cone = strtol((*is).c_str(), &evp, 10);
     float effic = strtod(evp + 1, &evp);
     result.insert(make_pair(effic, cone));
   }
   return result;
 }

References result().

◆ cuts()

Cuts IsolatorByNominalEfficiency::cuts ( float nominalEfficiency ) const

Definition at line 144 of file IsolatorByNominalEfficiency.cc.

                                                                     {
   vector<double> etaBounds = thresholds->bins();
   vector<double> coneSizes;
   vector<double> cutvalues;
   for (vector<double>::const_iterator it = etaBounds.begin(), itEnd = etaBounds.end(); it < itEnd; ++it) {
     float eta = (*it);
     int icone = bestConeForEfficiencyIndex(nominalEfficiency);
     coneSizes.push_back(theConesInfo.size(icone));
     NominalEfficiencyThresholds::ThresholdLocation location = {eta - 1.e-3f, icone};
     cutvalues.push_back(thresholds->thresholdValueForEfficiency(location, nominalEfficiency));
   }
   return Cuts(etaBounds, coneSizes, cutvalues);
 }

References bestConeForEfficiencyIndex(), muonisolation::NominalEfficiencyThresholds::bins(), PVValHelper::eta, f, EcalCondDBWriter_cfi::location, muonisolation::IsolatorByNominalEfficiency::ConeSizes::size(), theConesInfo, thresholds, and muonisolation::NominalEfficiencyThresholds::thresholdValueForEfficiency().

Referenced by L3NominalEfficiencyConfigurator::cuts().

◆ findPath()

string IsolatorByNominalEfficiency::findPath ( const std::string & fileName )

private

Definition at line 57 of file IsolatorByNominalEfficiency.cc.

                                                                    {
   // FIXME
   edm::FileInPath f(fileName);
   return f.fullPath();
 }

References f, and MillePedeFileConverter_cfg::fileName.

◆ result()

MuIsoBaseIsolator::Result IsolatorByNominalEfficiency::result	(	const DepositContainer &	deposits,
		const edm::Event *	= `nullptr`
	)		const

overridevirtual

Compute the deposit within the cone and return the isolation result.

Implements muonisolation::MuIsoBaseIsolator.

Definition at line 63 of file IsolatorByNominalEfficiency.cc.

                                                                                      {
   if (deposits.empty()) {
     cout << "IsolatorByNominalEfficiency: no deposit" << endl;
     return Result(resultType());  //FIXME
   }
  
   // To determine the threshold, the direction of the cone of the first
   // set of deposits is used.
   // For algorithms where different cone axis definitions are used
   // for different types deposits (eg. HCAL and ECAL deposits for
   // calorimeter isolation), the first one is used to determine the threshold
   // value!
   float theEta = deposits.back().dep->eta();
  
   // Try descending efficiency values to find the point where the candidate
   // becomes non isolated
  
   float nominalEfficiency = 1.;
   const float deltaeff = 0.005;
   const float mineff = deltaeff;
   for (float eff = .995; eff > mineff; eff -= deltaeff) {
     int cone = bestConeForEfficiencyIndex(eff);
     float coneSize = theConesInfo.size(cone);
     NominalEfficiencyThresholds::ThresholdLocation location = {theEta, cone};
     float thres = thresholds->thresholdValueForEfficiency(location, eff);
     float sumDep = weightedSum(deposits, coneSize);
     //      cout << "   Eff=" << eff
     //         << " eta=" << theEta
     //         << " cone=" << cone
     //         << " dR=" << coneSize
     //         << " thres=" << thres
     //         << " deposit=" << sumDep
     //         << " isolated=" << (sumDep < thres)
     //         << endl;
     if (sumDep > thres)
       break;
     nominalEfficiency = eff;
   }
   Result res(resultType());
   res.valFloat = nominalEfficiency;
   return res;
 }

References bestConeForEfficiencyIndex(), HLT_2018_cff::coneSize, gather_cfg::cout, CandIsolatorFromDeposits_cfi::deposits, EcalCondDBWriter_cfi::location, resultType(), muonisolation::IsolatorByNominalEfficiency::ConeSizes::size(), theConesInfo, thresholds, muonisolation::NominalEfficiencyThresholds::thresholdValueForEfficiency(), and weightedSum().

Referenced by cones().

◆ resultType()

ResultType muonisolation::IsolatorByNominalEfficiency::resultType ( ) const

inlineoverrideprivatevirtual

Implements muonisolation::MuIsoBaseIsolator.

Definition at line 66 of file IsolatorByNominalEfficiency.h.

66 { return ISOL_FLOAT_TYPE; }

References muonisolation::MuIsoBaseIsolator::ISOL_FLOAT_TYPE.

Referenced by result().

◆ weightedSum()

double IsolatorByNominalEfficiency::weightedSum	(	const DepositContainer &	deposits,
		float	dRcone
	)		const

privatevirtual

Definition at line 123 of file IsolatorByNominalEfficiency.cc.

                                                                                                     {
   double sumDep = 0;
  
   assert(deposits.size() == theWeights.size());
  
   vector<double>::const_iterator w = theWeights.begin();
   vector<double>::const_iterator dThresh = theDepThresholds.begin();
   for (DepositContainer::const_iterator dep = deposits.begin(); dep != deposits.end(); dep++) {
     if (dep->vetos != nullptr) {
       sumDep += dep->dep->depositAndCountWithin(dRcone, *dep->vetos, *dThresh).first * (*w);
     } else {
       sumDep += dep->dep->depositAndCountWithin(dRcone, Vetos(), *dThresh).first * (*w);
     }
     if (sumDep < 0.)
       sumDep = 0.;
     w++;
     dThresh++;
   }
   return sumDep;
 }