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

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

virtual	~IsolatorByNominalEfficiency ()

Public Member Functions inherited from muonisolation::MuIsoBaseIsolator
virtual Result	result (const DepositContainer &deposits, const reco::Candidate &muon, const edm::Event *=0) 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 *=0) 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)

virtual ResultType	resultType () const

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 21 of file IsolatorByNominalEfficiency.h.

Member Typedef Documentation

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

Definition at line 23 of file IsolatorByNominalEfficiency.h.

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

private

Definition at line 60 of file IsolatorByNominalEfficiency.h.

Constructor & Destructor Documentation

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

Constructor.

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

virtual

Definition at line 36 of file IsolatorByNominalEfficiency.cc.

References thresholds.

 {
   delete thresholds;
 }

Member Function Documentation

int IsolatorByNominalEfficiency::bestConeForEfficiencyIndex ( float eff_thr ) const

private

Definition at line 105 of file IsolatorByNominalEfficiency.cc.

References coneForEfficiency.

Referenced by cuts(), and result().

 {
 
   //FIXME use upper_bound
   int best_cone;
   if (coneForEfficiency.size() != 0) {
     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;
 }

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

private

Definition at line 42 of file IsolatorByNominalEfficiency.cc.

References result().

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

Cuts IsolatorByNominalEfficiency::cuts ( float nominalEfficiency ) const

Definition at line 145 of file IsolatorByNominalEfficiency.cc.

References bestConeForEfficiencyIndex(), muonisolation::NominalEfficiencyThresholds::bins(), eta(), f, 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);
 }

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

private

Definition at line 54 of file IsolatorByNominalEfficiency.cc.

References f, and edm::FileInPath::fullPath().

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

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

virtual

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

Implements muonisolation::MuIsoBaseIsolator.

Definition at line 62 of file IsolatorByNominalEfficiency.cc.

References bestConeForEfficiencyIndex(), gather_cfg::cout, resultType(), muonisolation::IsolatorByNominalEfficiency::ConeSizes::size(), theConesInfo, thresholds, muonisolation::NominalEfficiencyThresholds::thresholdValueForEfficiency(), muonisolation::MuIsoBaseIsolator::Result::valFloat, and weightedSum().

Referenced by cones().

                                                                                           {
 
   if (deposits.size()==0) {
     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;
 }

virtual ResultType muonisolation::IsolatorByNominalEfficiency::resultType ( ) const

inlineprivatevirtual

Implements muonisolation::MuIsoBaseIsolator.

Definition at line 65 of file IsolatorByNominalEfficiency.h.

References muonisolation::MuIsoBaseIsolator::ISOL_FLOAT_TYPE.

Referenced by result().

65 { return ISOL_FLOAT_TYPE;}

muonisolation::MuIsoBaseIsolator::ISOL_FLOAT_TYPE

Definition: MuIsoBaseIsolator.h:29

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

privatevirtual

Definition at line 123 of file IsolatorByNominalEfficiency.cc.

References 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 != 0){
       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;
 }