Public Member Functions
bool	operator() (const reco::Candidate &z) const

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

Private Types
typedef float(pat::Lepton< reco::Muon >::*	LeptonIsolationType) () const

typedef float(pat::GenericParticle::*	TrackIsolationType) () const

Private Attributes
Isolator	isolator_

LeptonIsolationType	leptonIsolation_

TrackIsolationType	trackIsolation_

Detailed Description

template<typename Isolator>
class ZToMuMuIsolationSelector< Isolator >

Definition at line 50 of file ZToMuMuIsolationSelector.cc.

Member Typedef Documentation

template<typename Isolator >

typedef float(pat::Lepton<reco::Muon>::* ZToMuMuIsolationSelector< Isolator >::LeptonIsolationType) () const

private

Definition at line 110 of file ZToMuMuIsolationSelector.cc.

template<typename Isolator >

typedef float(pat::GenericParticle::* ZToMuMuIsolationSelector< Isolator >::TrackIsolationType) () const

private

Definition at line 111 of file ZToMuMuIsolationSelector.cc.

Constructor & Destructor Documentation

template<typename Isolator >

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

inline

Definition at line 52 of file ZToMuMuIsolationSelector.cc.

References pat::GenericParticle::caloIso(), edm::errors::Configuration, pat::GenericParticle::ecalIso(), Exception, pat::GenericParticle::hcalIso(), AlCaHLTBitMon_QueryRunRegistry::string, and pat::GenericParticle::trackIso().

       : isolator_(cfg.template getParameter<double>("isoCut")) {
     std::string iso = cfg.template getParameter<std::string>("isolationType");
     if (iso == "track") {
       leptonIsolation_ = &pat::Lepton<reco::Muon>::trackIso;
       trackIsolation_ = &pat::GenericParticle::trackIso;
     } else if (iso == "ecal") {
       leptonIsolation_ = &pat::Lepton<reco::Muon>::ecalIso;
       trackIsolation_ = &pat::GenericParticle::ecalIso;
     } else if (iso == "hcal") {
       leptonIsolation_ = &pat::Lepton<reco::Muon>::hcalIso;
       trackIsolation_ = &pat::GenericParticle::hcalIso;
     } else if (iso == "calo") {
       leptonIsolation_ = &pat::Lepton<reco::Muon>::caloIso;
       trackIsolation_ = &pat::GenericParticle::caloIso;
     } else
       throw edm::Exception(edm::errors::Configuration) << "Invalid isolation type: " << iso << ". Valid types are:"
                                                        << "'track', 'ecal', 'hcal', 'calo'\n";
   }

Member Function Documentation

template<typename Isolator >

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

inline

Definition at line 71 of file ZToMuMuIsolationSelector.cc.

References 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();
     double iso0 = -1, iso1 = -1;
     const pat::Muon *mu0 = dynamic_cast<const pat::Muon *>(m0);
     if (mu0 != nullptr) {
       iso0 = ((*mu0).*(leptonIsolation_))();
     } else {
       const pat::GenericParticle *trk0 = dynamic_cast<const pat::GenericParticle *>(m0);
       if (trk0 != nullptr) {
         iso0 = ((*trk0).*(trackIsolation_))();
       } else {
         throw edm::Exception(edm::errors::InvalidReference)
             << "Candidate daughter #0 is neither pat::Muons nor pat::GenericParticle\n";
       }
     }
     const pat::Muon *mu1 = dynamic_cast<const pat::Muon *>(m1);
     if (mu1 != nullptr) {
       iso1 = ((*mu1).*(leptonIsolation_))();
     } else {
       const pat::GenericParticle *trk1 = dynamic_cast<const pat::GenericParticle *>(m1);
       if (trk1 != nullptr) {
         iso1 = ((*trk1).*(trackIsolation_))();
       } else {
         throw edm::Exception(edm::errors::InvalidReference)
             << "Candidate daughter #1 is neither pat::Muons nor pat::GenericParticle\n";
       }
     }
     bool pass = isolator_(iso0, iso1);
     return pass;
   }