muonisolation::IsolatorByNominalEfficiency Class Reference

#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.

References thresholds.

{ delete thresholds; }

Member Function Documentation

bestConeForEfficiencyIndex()

int IsolatorByNominalEfficiency::bestConeForEfficiencyIndex ( float  eff_thr ) const

private

Definition at line 107 of file IsolatorByNominalEfficiency.cc.

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

Referenced by cuts(), and result().

                                                                               {
  //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;
}

cones()

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

private

Definition at line 46 of file IsolatorByNominalEfficiency.cc.

References hltPFPuppi_cfi::cone, and result().

                                                                                                         {
  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;
}

cuts()

Cuts IsolatorByNominalEfficiency::cuts

(

float

nominalEfficiency

)

const

Definition at line 144 of file IsolatorByNominalEfficiency.cc.

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

Referenced by L3NominalEfficiencyConfigurator::cuts().

                                                                    {
  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);
}

findPath()

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

private

Definition at line 57 of file IsolatorByNominalEfficiency.cc.

References f, and MillePedeFileConverter_cfg::fileName.

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

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.

References bestConeForEfficiencyIndex(), hltPFPuppi_cfi::cone, HLT_2023v12_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().

                                                                                     {
  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;
}

resultType()

ResultType muonisolation::IsolatorByNominalEfficiency::resultType ( ) const

inlineoverrideprivatevirtual

Implements muonisolation::MuIsoBaseIsolator.

Definition at line 66 of file IsolatorByNominalEfficiency.h.

References muonisolation::MuIsoBaseIsolator::ISOL_FLOAT_TYPE.

Referenced by result().

{ return ISOL_FLOAT_TYPE; }

weightedSum()

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

privatevirtual

Definition at line 123 of file IsolatorByNominalEfficiency.cc.

References cms::cuda::assert(), CandIsolatorFromDeposits_cfi::deposits, JetPlusTrackCorrections_cff::dRcone, theDepThresholds, theWeights, and w().

Referenced by result().

                                                                                                    {
  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;
}

Member Data Documentation

coneForEfficiency

mapNomEff_Cone muonisolation::IsolatorByNominalEfficiency::coneForEfficiency

private

Definition at line 69 of file IsolatorByNominalEfficiency.h.

Referenced by bestConeForEfficiencyIndex().

theConesInfo

ConeSizes muonisolation::IsolatorByNominalEfficiency::theConesInfo

private

Definition at line 73 of file IsolatorByNominalEfficiency.h.

Referenced by cuts(), and result().

theDepThresholds

std::vector<double> muonisolation::IsolatorByNominalEfficiency::theDepThresholds

private

Definition at line 72 of file IsolatorByNominalEfficiency.h.

Referenced by weightedSum().

theWeights

std::vector<double> muonisolation::IsolatorByNominalEfficiency::theWeights

private

Definition at line 71 of file IsolatorByNominalEfficiency.h.

Referenced by weightedSum().

thresholds

NominalEfficiencyThresholds* muonisolation::IsolatorByNominalEfficiency::thresholds

private

Definition at line 70 of file IsolatorByNominalEfficiency.h.

Referenced by cuts(), result(), and ~IsolatorByNominalEfficiency().