Inheritance diagram for ZGlobalVsSAIsolationAnalyzer:

Public Types
typedef math::XYZVector	Vector

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Member Functions
	ZGlobalVsSAIsolationAnalyzer (const edm::ParameterSet &cfg)

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

virtual void	endJob () override

void	evaluate (const reco::Candidate *dau)

bool	isolated (const Direction &dir, const pat::IsoDeposit trkIsoDep, const pat::IsoDeposit ecalIsoDep, const pat::IsoDeposit *hcalIsoDep)

Private Attributes
double	alpha

double	beta

double	dREcal

double	dRHcal

double	dRTrk

double	dRVeto

double	etEcalThreshold

double	etHcalThreshold

double	isoCut_

double	ptThreshold

unsigned long	selGlobal_

unsigned long	selSA_

InputTag	src_

unsigned long	totGlobal_

unsigned long	totSA_

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 33 of file ZGlobalVsSAIsolationAnalyzer.cc.

Member Typedef Documentation

typedef math::XYZVector ZGlobalVsSAIsolationAnalyzer::Vector

Definition at line 35 of file ZGlobalVsSAIsolationAnalyzer.cc.

Constructor & Destructor Documentation

ZGlobalVsSAIsolationAnalyzer::ZGlobalVsSAIsolationAnalyzer ( const edm::ParameterSet & cfg )

Definition at line 52 of file ZGlobalVsSAIsolationAnalyzer.cc.

                                                                                  :
   src_(cfg.getParameter<InputTag>("src")),
   dRVeto(cfg.getParameter<double>("veto")),
   dRTrk(cfg.getParameter<double>("deltaRTrk")),
   dREcal(cfg.getParameter<double>("deltaREcal")),
   dRHcal(cfg.getParameter<double>("deltaRHcal")),
   ptThreshold(cfg.getParameter<double>("ptThreshold")),
   etEcalThreshold(cfg.getParameter<double>("etEcalThreshold")),
   etHcalThreshold(cfg.getParameter<double>("etHcalThreshold")),
   alpha(cfg.getParameter<double>("alpha")),
   beta(cfg.getParameter<double>("beta")),
   isoCut_(cfg.getParameter<double>("isoCut")),
   selGlobal_(0), selSA_(0), totGlobal_(0), totSA_(0) {
 }

Member Function Documentation

void ZGlobalVsSAIsolationAnalyzer::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 105 of file ZGlobalVsSAIsolationAnalyzer.cc.

References reco::Candidate::daughter(), WZMuSkim_cff::dimuons, evaluate(), i, and src_.

                                                                                             {
   Handle<CandidateView> dimuons;
   event.getByLabel(src_, dimuons);
   for(unsigned int i=0; i< dimuons->size(); ++ i) {
     const Candidate & zmm = (* dimuons)[i];
     evaluate(zmm.daughter(0));
     evaluate(zmm.daughter(1));
   }    
 }

void ZGlobalVsSAIsolationAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 115 of file ZGlobalVsSAIsolationAnalyzer.cc.

References gather_cfg::cout, selGlobal_, selSA_, mathSSE::sqrt(), totGlobal_, and totSA_.

                                           {
   cout << "Isolation efficiency report:" << endl;
   double eff, err;
   eff =  double(selGlobal_)/double(totGlobal_);
   err = sqrt(eff*(1.-eff)/double(totGlobal_));
   cout <<"Global: " << selGlobal_ << "/" << totGlobal_ << " = " <<eff <<"+/-" << err<< endl;
   eff = double(selSA_)/double(totSA_);
   err = sqrt(eff*(1.-eff)/double(totSA_));  
   cout <<"St.Al.: " << selSA_ << "/" << totSA_ << " = " << eff <<"+/-" << err << endl;
 }

void ZGlobalVsSAIsolationAnalyzer::evaluate ( const reco::Candidate * dau )

private

Definition at line 84 of file ZGlobalVsSAIsolationAnalyzer.cc.

References dir, pat::EcalIso, reco::LeafCandidate::eta(), edm::hlt::Exception, reco::Candidate::get(), pat::HcalIso, edm::errors::InvalidReference, pat::Lepton< LeptonType >::isoDeposit(), isolated(), reco::Candidate::masterClone(), RPCpg::mu, reco::LeafCandidate::phi(), selGlobal_, selSA_, totGlobal_, totSA_, and pat::TrackIso.

Referenced by analyze().

                                                                     {
   const pat::Muon * mu = dynamic_cast<const pat::Muon *>(&*dau->masterClone());
   if(mu == 0) throw Exception(errors::InvalidReference) << "Daughter is not a muon!\n";
   const pat::IsoDeposit * trkIsoDep = mu->isoDeposit(pat::TrackIso);
   const pat::IsoDeposit * ecalIsoDep = mu->isoDeposit(pat::EcalIso);
   const pat::IsoDeposit * hcalIsoDep = mu->isoDeposit(pat::HcalIso);
   // global muon
   {
     Direction dir = Direction(mu->eta(), mu->phi());
     if(isolated(dir, trkIsoDep, ecalIsoDep, hcalIsoDep)) selGlobal_++;
     totGlobal_++;
   }
   // stand-alone
   {
     TrackRef sa = dau->get<TrackRef,reco::StandAloneMuonTag>();
     Direction dir = Direction(sa->eta(), sa->phi());
     if(isolated(dir, trkIsoDep, ecalIsoDep, hcalIsoDep)) selSA_++;
     totSA_++;   
   }
 }

bool ZGlobalVsSAIsolationAnalyzer::isolated	(	const Direction &	dir,
		const pat::IsoDeposit *	trkIsoDep,
		const pat::IsoDeposit *	ecalIsoDep,
		const pat::IsoDeposit *	hcalIsoDep
	)

private

Definition at line 67 of file ZGlobalVsSAIsolationAnalyzer.cc.

References alpha, beta, dREcal, dRHcal, dRTrk, dRVeto, etEcalThreshold, etHcalThreshold, isoCut_, ptThreshold, and reco::IsoDeposit::sumWithin().

Referenced by evaluate().

                                                                                             {
   IsoDeposit::AbsVetos vetoTrk, vetoEcal, vetoHcal;
   vetoTrk.push_back(new ConeVeto(dir, dRVeto));
   vetoTrk.push_back(new ThresholdVeto(ptThreshold));
   vetoEcal.push_back(new ConeVeto(dir, 0.));
   vetoEcal.push_back(new ThresholdVeto(etEcalThreshold));
   vetoHcal.push_back(new ConeVeto(dir, 0.));
   vetoHcal.push_back(new ThresholdVeto(etHcalThreshold));
   
   double trkIso = trkIsoDep->sumWithin(dir, dRTrk, vetoTrk);
   double ecalIso = ecalIsoDep->sumWithin(dir, dREcal, vetoEcal);
   double hcalIso = hcalIsoDep->sumWithin(dir, dRHcal, vetoHcal);
   double iso = alpha*((0.5*(1+beta)*ecalIso) + (0.5*(1-beta)*hcalIso)) + (1-alpha)*trkIso;
   return iso < isoCut_;      
 }