#include <ZtoMMEventSelector.h>

Inheritance diagram for ZtoMMEventSelector:

Public Member Functions
bool	filter (edm::Event &, edm::EventSetup const &) override

	ZtoMMEventSelector (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
const edm::InputTag	bsTag_

const edm::EDGetTokenT< reco::BeamSpot >	bsToken_

const double	maxD0_

const double	maxDz_

const double	maxEta_

const double	maxInvMass_

const double	maxIso_

const double	maxNormChi2_

const int	minChambers_

const double	minInvMass_

const int	minMatchedStations_

const int	minMatches_

const int	minPixelHits_

const double	minPt_

const double	minPtHighest_

const int	minStripHits_

const edm::InputTag	muonTag_

const edm::EDGetTokenT< reco::MuonCollection >	muonToken_

bool	verbose_

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDFilterBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 10 of file ZtoMMEventSelector.h.

Constructor & Destructor Documentation

ZtoMMEventSelector::ZtoMMEventSelector ( const edm::ParameterSet & ps )

explicit

Definition at line 19 of file ZtoMMEventSelector.cc.

                                                                : 
   verbose_(ps.getUntrackedParameter<bool>("verbose", false)),
   muonTag_(ps.getUntrackedParameter<edm::InputTag>("muonInputTag", edm::InputTag("muons"))),
   bsTag_(ps.getUntrackedParameter<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot"))),
   muonToken_(consumes<reco::MuonCollection>(muonTag_)),
   bsToken_(consumes<reco::BeamSpot>(bsTag_)),
   maxEta_(ps.getUntrackedParameter<double>("maxEta", 2.1)),
   minPt_(ps.getUntrackedParameter<double>("minPt", 5)),
   maxNormChi2_(ps.getUntrackedParameter<double>("maxNormChi2", 10)),
   maxD0_(ps.getUntrackedParameter<double>("maxD0", 0.02)),
   maxDz_(ps.getUntrackedParameter<double>("maxDz", 20.)),
   minPixelHits_(ps.getUntrackedParameter<uint32_t>("minPixelHits", 1)),
   minStripHits_(ps.getUntrackedParameter<uint32_t>("minStripHits", 8)),
   minChambers_(ps.getUntrackedParameter<uint32_t>("minChambers", 2)),
   minMatches_(ps.getUntrackedParameter<uint32_t>("minMatches", 2)),
   minMatchedStations_(ps.getUntrackedParameter<double>("minMatchedStations", 2)),
   maxIso_(ps.getUntrackedParameter<double>("maxIso", 0.3)),
   minPtHighest_(ps.getUntrackedParameter<double>("minPtHighest", 24)),
   minInvMass_(ps.getUntrackedParameter<double>("minInvMass", 60)),
   maxInvMass_(ps.getUntrackedParameter<double>("maxInvMass", 120))
 {
 }

Member Function Documentation

bool ZtoMMEventSelector::filter	(	edm::Event &	iEvent,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Implements edm::stream::EDFilterBase.

Definition at line 41 of file ZtoMMEventSelector.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, bsTag_, bsToken_, vertices_cff::chi2, DEFINE_FWK_MODULE, edm::Event::getByToken(), muon::GlobalMuonPromptTight, muon::isGoodMuon(), edm::HandleBase::isValid(), list(), ResonanceBuilder::mass, SiStripPI::max, maxD0_, maxDz_, maxEta_, maxInvMass_, maxIso_, maxNormChi2_, minChambers_, minMatchedStations_, minMatches_, minPixelHits_, minPt_, minPtHighest_, minStripHits_, RPCpg::mu, muonTag_, muonToken_, ndof, reco::HitPattern::numberOfValidPixelHits(), reco::HitPattern::numberOfValidStripHits(), reco::BeamSpot::position(), reco::tau::disc::Pt(), reco::MuonPFIsolation::sumChargedParticlePt, reco::MuonPFIsolation::sumNeutralHadronEt, reco::MuonPFIsolation::sumPhotonEt, and reco::MuonPFIsolation::sumPUPt.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

                                                                              {
   // Read Muon Collection
   edm::Handle<reco::MuonCollection> muonColl;
   iEvent.getByToken(muonToken_, muonColl);
 
   // and the beamspot
   edm::Handle<reco::BeamSpot> beamSpot;
   iEvent.getByToken(bsToken_, beamSpot);
 
   std::vector<TLorentzVector> list; 
   if (muonColl.isValid()) {
     for (auto const& mu : *muonColl) {
       if (!mu.isGlobalMuon()) continue;
       if (!mu.isPFMuon()) continue;
       if (std::fabs(mu.eta()) >= maxEta_) continue;
       if (mu.pt() < minPt_) continue;
 
       reco::TrackRef gtk = mu.globalTrack();
       double chi2 = gtk->chi2();
       double ndof = gtk->ndof();
       double chbyndof = (ndof > 0) ? chi2/ndof : 0;
       if (chbyndof >= maxNormChi2_) continue;
 
       reco::TrackRef tk = mu.innerTrack();
       double trkd0 = tk->d0();
       double trkdz = tk->dz();
       if (beamSpot.isValid()) {
     trkd0 = -(tk->dxy(beamSpot->position()));
     if (std::fabs(trkd0) >= maxD0_) continue;
     trkdz = tk->dz(beamSpot->position()); 
     if (std::fabs(trkdz) >= maxDz_) continue;
       }
       else {
     edm::LogError("ZtoMMEventSelector") << "Error >> Failed to get BeamSpot for label: "
                         << bsTag_;
       }
 
       const reco::HitPattern& hitp = gtk->hitPattern(); 
       if (hitp.numberOfValidPixelHits() < minPixelHits_) continue;
       if (hitp.numberOfValidStripHits() < minStripHits_) continue;
 
       // Hits/section in the muon chamber
       if (mu.numberOfChambers() < minChambers_) continue;
       if (mu.numberOfMatches() < minMatches_) continue;
       if (mu.numberOfMatchedStations() < minMatchedStations_) continue;
       if (!muon::isGoodMuon(mu, muon::GlobalMuonPromptTight)) continue;
 
       // PF Isolation      
       const reco::MuonPFIsolation& pfIso04 = mu.pfIsolationR04();
       double absiso = pfIso04.sumChargedParticlePt + std::max(0.0, pfIso04.sumNeutralHadronEt + pfIso04.sumPhotonEt - 0.5 * pfIso04.sumPUPt);
       if (absiso/mu.pt() > maxIso_) continue;
 
       TLorentzVector lv;
       lv.SetPtEtaPhiE(mu.pt(), mu.eta(), mu.phi(), mu.energy());
       list.push_back(lv); 
     }
   }
   else {
     edm::LogError("ZtoMMEventSelector") << "Error >> Failed to get MuonCollection for label: " << muonTag_;
     return false;
   }
 
   if (list.size() < 2) return false;
   if (list[0].Pt() < minPtHighest_) return false;  
   TLorentzVector zv = list[0] + list[1];
   double mass = zv.M();
   if (mass < minInvMass_ || mass > maxInvMass_) return false;
     
   return true;
 }

Member Data Documentation

const edm::InputTag ZtoMMEventSelector::bsTag_

private

Definition at line 19 of file ZtoMMEventSelector.h.

Referenced by filter().

const edm::EDGetTokenT<reco::BeamSpot> ZtoMMEventSelector::bsToken_

private

Definition at line 21 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::maxD0_

private

Definition at line 26 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::maxDz_

private

Definition at line 27 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::maxEta_

private

Definition at line 23 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::maxInvMass_

private

Definition at line 36 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::maxIso_

private

Definition at line 33 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::maxNormChi2_

private

Definition at line 25 of file ZtoMMEventSelector.h.

Referenced by filter().

const int ZtoMMEventSelector::minChambers_

private

Definition at line 30 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::minInvMass_

private

Definition at line 35 of file ZtoMMEventSelector.h.

const int ZtoMMEventSelector::minMatchedStations_

private

Definition at line 32 of file ZtoMMEventSelector.h.

Referenced by filter().

const int ZtoMMEventSelector::minMatches_

private

Definition at line 31 of file ZtoMMEventSelector.h.

Referenced by filter().

const int ZtoMMEventSelector::minPixelHits_

private

Definition at line 28 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::minPt_

private

Definition at line 24 of file ZtoMMEventSelector.h.

Referenced by filter().

const double ZtoMMEventSelector::minPtHighest_

private

Definition at line 34 of file ZtoMMEventSelector.h.

Referenced by filter().

const int ZtoMMEventSelector::minStripHits_

private

Definition at line 29 of file ZtoMMEventSelector.h.

Referenced by filter().

const edm::InputTag ZtoMMEventSelector::muonTag_

private

Definition at line 18 of file ZtoMMEventSelector.h.

Referenced by filter().

const edm::EDGetTokenT<reco::MuonCollection> ZtoMMEventSelector::muonToken_

private

Definition at line 20 of file ZtoMMEventSelector.h.

Referenced by filter().

bool ZtoMMEventSelector::verbose_

private

Definition at line 17 of file ZtoMMEventSelector.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation