Public Member Functions
double	candIsolation (const reco::Candidate *c) const

template<typename T >
double	isolation (const T *t) const

bool	operator() (const reco::Candidate &z) const

	ZToMuMuIsoDepositSelector (const edm::ParameterSet &cfg, edm::ConsumesCollector &iC)

Private Attributes
double	alpha

double	beta

double	dREcal

double	dRHcal

double	dRTrk

double	dRVetoTrk

double	etEcalThreshold

double	etHcalThreshold

Isolator	isolator_

double	ptThreshold

bool	relativeIsolation

Detailed Description

template<typename Isolator>
class ZToMuMuIsoDepositSelector< Isolator >

Definition at line 58 of file ZToMuMuIsoDepositSelector.cc.

Constructor & Destructor Documentation

template<typename Isolator >

ZToMuMuIsoDepositSelector< Isolator >::ZToMuMuIsoDepositSelector	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

inline

Definition at line 60 of file ZToMuMuIsoDepositSelector.cc.

       : isolator_(cfg.template getParameter<double>("isoCut")),
         ptThreshold(cfg.getUntrackedParameter<double>("ptThreshold")),
         etEcalThreshold(cfg.getUntrackedParameter<double>("etEcalThreshold")),
         etHcalThreshold(cfg.getUntrackedParameter<double>("etHcalThreshold")),
         dRVetoTrk(cfg.getUntrackedParameter<double>("deltaRVetoTrk")),
         dRTrk(cfg.getUntrackedParameter<double>("deltaRTrk")),
         dREcal(cfg.getUntrackedParameter<double>("deltaREcal")),
         dRHcal(cfg.getUntrackedParameter<double>("deltaRHcal")),
         alpha(cfg.getUntrackedParameter<double>("alpha")),
         beta(cfg.getUntrackedParameter<double>("beta")),
         relativeIsolation(cfg.template getParameter<bool>("relativeIsolation")) {}

Member Function Documentation

template<typename Isolator >

double ZToMuMuIsoDepositSelector< Isolator >::candIsolation ( const reco::Candidate * c ) const

inline

Definition at line 104 of file ZToMuMuIsoDepositSelector.cc.

References HltBtagPostValidation_cff::c, Exception, edm::errors::InvalidReference, and amptDefaultParameters_cff::mu.

                                                      {
     const pat::Muon* mu = dynamic_cast<const pat::Muon*>(c);
     if (mu != nullptr)
       return isolation(mu);
     const pat::GenericParticle* trk = dynamic_cast<const pat::GenericParticle*>(c);
     if (trk != nullptr)
       return isolation(trk);
     throw edm::Exception(edm::errors::InvalidReference)
         << "Candidate daughter #0 is neither pat::Muons nor pat::GenericParticle\n";
     return -1;
   }

template<typename Isolator >

template<typename T >

double ZToMuMuIsoDepositSelector< Isolator >::isolation ( const T * t ) const

inline

Definition at line 74 of file ZToMuMuIsoDepositSelector.cc.

References zMuMuMuonUserData::alpha, zMuMuMuonUserData::beta, DeadROC_duringRun::dir, zMuMuMuonUserData::dREcal, zMuMuMuonUserData::dRHcal, zMuMuMuonUserData::dRTrk, zMuMuMuonUserData::dRVetoTrk, pat::EcalIso, zMuMuMuonUserData::etEcalThreshold, zMuMuMuonUserData::etHcalThreshold, pat::HcalIso, ewkMuLumiMonitorDQM_cfi::ptThreshold, ZMuMuCategoriesNtuples_cff::relativeIsolation, reco::IsoDeposit::sumWithin(), pat::TrackIso, and ZMuMuCategoriesSequences_cff::trkIso.

                                      {
     const pat::IsoDeposit* trkIso = t->isoDeposit(pat::TrackIso);
     const pat::IsoDeposit* ecalIso = t->isoDeposit(pat::EcalIso);
     const pat::IsoDeposit* hcalIso = t->isoDeposit(pat::HcalIso);
 
     Direction dir = Direction(t->eta(), t->phi());
 
     IsoDeposit::AbsVetos vetosTrk;
     vetosTrk.push_back(new ConeVeto(dir, dRVetoTrk));
     vetosTrk.push_back(new ThresholdVeto(ptThreshold));
 
     IsoDeposit::AbsVetos vetosEcal;
     vetosEcal.push_back(new ConeVeto(dir, 0.));
     vetosEcal.push_back(new ThresholdVeto(etEcalThreshold));
 
     IsoDeposit::AbsVetos vetosHcal;
     vetosHcal.push_back(new ConeVeto(dir, 0.));
     vetosHcal.push_back(new ThresholdVeto(etHcalThreshold));
 
     double isovalueTrk = (trkIso->sumWithin(dRTrk, vetosTrk));
     double isovalueEcal = (ecalIso->sumWithin(dREcal, vetosEcal));
     double isovalueHcal = (hcalIso->sumWithin(dRHcal, vetosHcal));
 
     double iso =
         alpha * (((1 + beta) / 2 * isovalueEcal) + ((1 - beta) / 2 * isovalueHcal)) + ((1 - alpha) * isovalueTrk);
     if (relativeIsolation)
       iso /= t->pt();
     return iso;
   }

template<typename Isolator >

bool ZToMuMuIsoDepositSelector< Isolator >::operator() ( const reco::Candidate & z ) const

inline

Definition at line 115 of file ZToMuMuIsoDepositSelector.cc.

References candIsolation(), reco::Candidate::daughter(), Exception, reco::Candidate::hasMasterClone(), edm::errors::InvalidReference, reco::Candidate::masterClone(), and reco::Candidate::numberOfDaughters().

                                                 {
     if (z.numberOfDaughters() != 2)
       throw edm::Exception(edm::errors::InvalidReference)
           << "Candidate has " << z.numberOfDaughters() << " daughters, 2 expected\n";
     const reco::Candidate* dau0 = z.daughter(0);
     const reco::Candidate* dau1 = z.daughter(1);
     if (!(dau0->hasMasterClone() && dau1->hasMasterClone()))
       throw edm::Exception(edm::errors::InvalidReference) << "Candidate daughters have no master clone\n";
     const reco::Candidate *m0 = &*dau0->masterClone(), *m1 = &*dau1->masterClone();
     return isolator_(candIsolation(m0), candIsolation(m1));
   }