#include <CutsIsolatorWithCorrection.h>

Inheritance diagram for CutsIsolatorWithCorrection:

Public Member Functions
	CutsIsolatorWithCorrection (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)

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

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

ResultType	resultType () const 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	~MuIsoBaseIsolator ()

Private Member Functions
double	depSum (const DepositContainer &deposits, double dr, double corr) const

Private Attributes
bool	theAndOrCuts

bool	theCutAbsIso

bool	theCutRelativeIso

muonisolation::Cuts	theCuts

muonisolation::Cuts	theCutsRel

double	theEffAreaSFBarrel

double	theEffAreaSFEndcap

bool	theReturnAbsoluteSum

bool	theReturnRelativeSum

double	theRhoMax

double	theRhoScaleBarrel

double	theRhoScaleEndcap

edm::EDGetTokenT< double >	theRhoToken

bool	theUseRhoCorrection

Additional Inherited Members
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

Detailed Description

Definition at line 9 of file CutsIsolatorWithCorrection.h.

Constructor & Destructor Documentation

◆ CutsIsolatorWithCorrection()

CutsIsolatorWithCorrection::CutsIsolatorWithCorrection	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &&	iC
	)

Definition at line 8 of file CutsIsolatorWithCorrection.cc.

References Exception, theCutAbsIso, and theCutRelativeIso.

     : theCuts(par.getParameter<std::vector<double> >("EtaBounds"),
               par.getParameter<std::vector<double> >("ConeSizes"),
               par.getParameter<std::vector<double> >("Thresholds")),
       theCutsRel(par.getParameter<std::vector<double> >("EtaBoundsRel"),
                  par.getParameter<std::vector<double> >("ConeSizesRel"),
                  par.getParameter<std::vector<double> >("ThresholdsRel")),
       theCutAbsIso(par.getParameter<bool>("CutAbsoluteIso")),
       theCutRelativeIso(par.getParameter<bool>("CutRelativeIso")),
       theUseRhoCorrection(par.getParameter<bool>("UseRhoCorrection")),
       theRhoToken(iC.consumes<double>(par.getParameter<edm::InputTag>("RhoSrc"))),
       theRhoMax(par.getParameter<double>("RhoMax")),
       theRhoScaleBarrel(par.getParameter<double>("RhoScaleBarrel")),
       theRhoScaleEndcap(par.getParameter<double>("RhoScaleEndcap")),
       theEffAreaSFBarrel(par.getParameter<double>("EffAreaSFBarrel")),
       theEffAreaSFEndcap(par.getParameter<double>("EffAreaSFEndcap")),
       theReturnAbsoluteSum(par.getParameter<bool>("ReturnAbsoluteSum")),
       theReturnRelativeSum(par.getParameter<bool>("ReturnRelativeSum")),
       theAndOrCuts(par.getParameter<bool>("AndOrCuts")) {
   if (!(theCutAbsIso || theCutRelativeIso))
     throw cms::Exception("BadConfiguration")
         << "Something has to be cut: set either CutAbsoluteIso or CutRelativeIso to true";
 }

Member Function Documentation

◆ depSum()

double CutsIsolatorWithCorrection::depSum	(	const DepositContainer &	deposits,
		double	dr,
		double	corr
	)		const

private

Definition at line 32 of file CutsIsolatorWithCorrection.cc.

References alignCSCRings::corr, CandIsolatorFromDeposits_cfi::deposits, and l1ctLayer1_cff::dr.

Referenced by result().

                                                                                                         {
   double dephlt = -corr;
   unsigned int nDeps = deposits.size();
   //  edm::LogWarning("CutsIsolatorWithCorrection::depSumIn")
   //    << "add nDeposit "<< nDeps<< " \t dr "<<dr<<" \t corr "<<corr;
   for (unsigned int iDep = 0; iDep < nDeps; ++iDep) {
     double lDep = deposits[iDep].dep->depositWithin(dr);
     dephlt += lDep;
     //    edm::LogWarning("CutsIsolatorWithCorrection::depSumIDep")
     //      <<"dep "<<iDep<<" \t added "<<lDep<<" \t sumnow "<<dephlt;
   }
 
   return dephlt;
 }

◆ result() [1/2]

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

inlineoverridevirtual

Compute and return the isolation variable.

Implements muonisolation::MuIsoBaseIsolator.

Definition at line 15 of file CutsIsolatorWithCorrection.h.

References muonisolation::MuIsoBaseIsolator::ISOL_BOOL_TYPE, and muonisolation::MuIsoBaseIsolator::Result::valBool.

                                                                                             {
     Result answer(ISOL_BOOL_TYPE);
     answer.valBool = false;
     // fail miserably...
     return answer;
   }

◆ result() [2/2]

MuIsoBaseIsolator::Result CutsIsolatorWithCorrection::result	(	const DepositContainer &	deposits,
		const reco::Track &	muon,
		const edm::Event *	= `nullptr`
	)		const

overridevirtual

Compute and return the isolation variable, with vetoes and the muon.

Reimplemented from muonisolation::MuIsoBaseIsolator.

Definition at line 47 of file CutsIsolatorWithCorrection.cc.

References L1NNTauProducer_cff::conesize, muonisolation::Cuts::CutSpec::conesize, CandIsolatorFromDeposits_cfi::deposits, depSum(), reco::TrackBase::eta(), makeMEIFBenchmarkPlots::ev, nano_mu_digi_cff::float, muonisolation::MuIsoBaseIsolator::ISOL_BOOL_TYPE, pi, edm::Handle< T >::product(), reco::TrackBase::pt(), rho, HLT_2023v12_cff::rhoScale, theAndOrCuts, theCutAbsIso, theCutRelativeIso, theCuts, theCutsRel, theEffAreaSFBarrel, theEffAreaSFEndcap, theReturnAbsoluteSum, theReturnRelativeSum, theRhoMax, theRhoScaleBarrel, theRhoScaleEndcap, theRhoToken, theUseRhoCorrection, muonisolation::Cuts::CutSpec::threshold, muonisolation::MuIsoBaseIsolator::Result::valBool, and muonisolation::MuIsoBaseIsolator::Result::valFloat.

                                                                                        {
   Result answer(ISOL_BOOL_TYPE);
 
   bool absDecision = false;
   bool relDecision = false;
 
   double rho = 0.0;
   double effAreaSF = 1.0;
 
   static const double pi = 3.14159265358979323846;
 
   //  edm::LogWarning("CutsIsolatorWithCorrection::resultIn")
   //    <<"Start tk.pt "<<tk.pt()<<" \t tk.eta "<<tk.eta()<<" \t tk.phi "<<tk.phi();
 
   if (theUseRhoCorrection) {
     edm::Handle<double> rhoHandle;
     ev->getByToken(theRhoToken, rhoHandle);
     rho = *(rhoHandle.product());
     if (rho < 0.0)
       rho = 0.0;
     double rhoScale = fabs(tk.eta()) > 1.442 ? theRhoScaleEndcap : theRhoScaleBarrel;
     effAreaSF = fabs(tk.eta()) > 1.442 ? theEffAreaSFEndcap : theEffAreaSFBarrel;
     //    edm::LogWarning("CutsIsolatorWithCorrection::resultInRho")
     //      << "got rho "<<rho<<" vs max "<<theRhoMax<<" will scale by "<<rhoScale;
     if (rho > theRhoMax) {
       rho = theRhoMax;
     }
     rho = rho * rhoScale;
     //    edm::LogWarning("CutsIsolatorWithCorrection::resultOutRho")<<" final rho "<<rho;
   }
 
   if (theCutAbsIso) {
     muonisolation::Cuts::CutSpec cuts_here = theCuts(tk.eta());
     double conesize = cuts_here.conesize;
     double dephlt = depSum(deposits, conesize, rho * conesize * conesize * pi * effAreaSF);
     if (theReturnAbsoluteSum)
       answer.valFloat = (float)dephlt;
     if (dephlt < cuts_here.threshold) {
       absDecision = true;
     } else {
       absDecision = false;
     }
     //    edm::LogWarning("CutsIsolatorWithCorrection::resultOutAbsIso")
     //      <<"compared dephlt "<<dephlt<<" \t with "<<cuts_here.threshold;
   } else
     absDecision = true;
 
   if (theCutRelativeIso) {
     muonisolation::Cuts::CutSpec cuts_here = theCutsRel(tk.eta());
     double conesize = cuts_here.conesize;
     double dephlt = depSum(deposits, conesize, rho * conesize * conesize * pi * effAreaSF) / tk.pt();
     if (theReturnRelativeSum)
       answer.valFloat = (float)dephlt;
     if (dephlt < cuts_here.threshold) {
       relDecision = true;
     } else {
       relDecision = false;
     }
     //    edm::LogWarning("CutsIsolatorWithCorrection::resultOutRelIso")
     //      <<"compared dephlt "<<dephlt<<" \t with "<<cuts_here.threshold;
   } else
     relDecision = true;
 
   if (theAndOrCuts) {
     answer.valBool = absDecision && relDecision;
   } else {
     answer.valBool = absDecision || relDecision;
   }
 
   //  edm::LogWarning("CutsIsolatorWithCorrection::result")
   //    <<"isAbsIsolated "<<absDecision<<" \t isRelIsolated  "<<relDecision
   //    <<" \t combined with AND "<<theAndOrCuts
   //    <<" = "  << answer.valBool;
 
   return answer;
 }