Inheritance diagram for ZGlobalVsSAIsolationAnalyzer:

Public Types
typedef math::XYZVector	Vector

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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_

EDGetTokenT< CandidateView >	srcToken_

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

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Detailed Description

Definition at line 33 of file ZGlobalVsSAIsolationAnalyzer.cc.

Member Typedef Documentation

◆ Vector

typedef math::XYZVector ZGlobalVsSAIsolationAnalyzer::Vector

Definition at line 35 of file ZGlobalVsSAIsolationAnalyzer.cc.

Constructor & Destructor Documentation

◆ ZGlobalVsSAIsolationAnalyzer()

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

Definition at line 55 of file ZGlobalVsSAIsolationAnalyzer.cc.

     : srcToken_(consumes<CandidateView>(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

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 115 of file ZGlobalVsSAIsolationAnalyzer.cc.

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

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

◆ endJob()

void ZGlobalVsSAIsolationAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 125 of file ZGlobalVsSAIsolationAnalyzer.cc.

                                           {
   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;
 }

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

◆ evaluate()

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

private

Definition at line 91 of file ZGlobalVsSAIsolationAnalyzer.cc.

                                                                     {
   const pat::Muon* mu = dynamic_cast<const pat::Muon*>(&*dau->masterClone());
   if (mu == nullptr)
     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_++;
   }
 }

References DeadROC_duringRun::dir, pat::EcalIso, Exception, reco::Candidate::get(), pat::HcalIso, edm::errors::InvalidReference, isolated(), reco::Candidate::masterClone(), amptDefaultParameters_cff::mu, selGlobal_, selSA_, totGlobal_, totSA_, and pat::TrackIso.

Referenced by analyze().

◆ isolated()

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

private

Definition at line 72 of file ZGlobalVsSAIsolationAnalyzer.cc.

                                                                              {
   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_;
 }