MuonIDTableProducer Class Reference

Inheritance diagram for MuonIDTableProducer:

Public Member Functions
	MuonIDTableProducer (const edm::ParameterSet &iConfig)
	~MuonIDTableProducer () override
Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
	~EDProducerBase () 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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::StreamID, edm::Event &, edm::EventSetup const &) const override

Static Private Member Functions
static bool	isMediumMuonHIP (const pat::Muon &muon)
static bool	isSoftMuonHIP (const pat::Muon &muon, const reco::Vertex &vtx)
static bool	isTrackerHighPt (const pat::Muon &mu, const reco::Vertex &primaryVertex)

Private Attributes
std::string	name_
edm::EDGetTokenT< std::vector< pat::Muon > >	src_
edm::EDGetTokenT< std::vector< reco::Vertex > >	srcVtx_

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	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
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 20 of file MuonIDTableProducer.cc.

Constructor & Destructor Documentation

MuonIDTableProducer::MuonIDTableProducer ( const edm::ParameterSet &  iConfig )

inlineexplicit

Definition at line 22 of file MuonIDTableProducer.cc.

                                                                     :
            name_(iConfig.getParameter<std::string>("name")),
            src_(consumes<std::vector<pat::Muon>>(iConfig.getParameter<edm::InputTag>("muons"))),
            srcVtx_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("vertices")))
        {
            produces<nanoaod::FlatTable>();
        }

MuonIDTableProducer::~MuonIDTableProducer ( )

inlineoverride

Definition at line 30 of file MuonIDTableProducer.cc.

{};

Member Function Documentation

static void MuonIDTableProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

inlinestatic

Definition at line 32 of file MuonIDTableProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), produce(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                  {
            edm::ParameterSetDescription desc;
            desc.add<edm::InputTag>("muons")->setComment("input muon collection");
            desc.add<edm::InputTag>("vertices", edm::InputTag("offlineSlimmedPrimaryVertices"))->setComment("input vertex collection, for dxy/dz");
            desc.add<std::string>("name")->setComment("name of the muon nanoaod::FlatTable we are extending with IDs");
            descriptions.add("muonIDTable", desc);
        }

bool MuonIDTableProducer::isMediumMuonHIP ( const pat::Muon &  muon )

staticprivate

Definition at line 85 of file MuonIDTableProducer.cc.

References reco::MuonQuality::chi2LocalPosition, reco::Muon::combinedQuality(), pat::Muon::globalTrack(), pat::Muon::innerTrack(), reco::Muon::isGlobalMuon(), muon::isLooseMuon(), muon::segmentCompatibility(), and reco::MuonQuality::trkKink.

Referenced by produce().

                                                           {
      bool goodGlob = mu.isGlobalMuon() && 
                      mu.globalTrack()->normalizedChi2() < 3 && 
                      mu.combinedQuality().chi2LocalPosition < 12 && 
                      mu.combinedQuality().trkKink < 20; 
      bool isMedium = muon::isLooseMuon(mu) && 
                      mu.innerTrack()->validFraction() > 0.49 && 
                      muon::segmentCompatibility(mu) > (goodGlob ? 0.303 : 0.451); 
      return isMedium; 
}

bool MuonIDTableProducer::isSoftMuonHIP ( const pat::Muon &  muon, const reco::Vertex &  vtx  )

staticprivate

Definition at line 96 of file MuonIDTableProducer.cc.

References funct::abs(), pat::Muon::innerTrack(), muon::isGoodMuon(), reco::Vertex::position(), and muon::TMOneStationTight.

Referenced by produce().

                                                                                {
    return muon::isGoodMuon(mu, muon::TMOneStationTight) && 
           mu.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 &&
           mu.innerTrack()->hitPattern().pixelLayersWithMeasurement() > 0 && 
           std::abs(mu.innerTrack()->dxy(vtx.position())) < 0.3 && 
           std::abs(mu.innerTrack()->dz(vtx.position())) < 20.;
}

bool MuonIDTableProducer::isTrackerHighPt ( const pat::Muon &  mu, const reco::Vertex &  primaryVertex  )

staticprivate

Definition at line 104 of file MuonIDTableProducer.cc.

References funct::abs(), DEFINE_FWK_MODULE, pat::Muon::innerTrack(), pat::Muon::muonBestTrack(), reco::Muon::numberOfMatchedStations(), and reco::Vertex::position().

Referenced by produce().

                                                                                              {
        return ( mu.numberOfMatchedStations() > 1 
                         && (mu.muonBestTrack()->ptError()/mu.muonBestTrack()->pt()) < 0.3 
                         && std::abs(mu.muonBestTrack()->dxy(primaryVertex.position())) < 0.2 
                         && std::abs(mu.muonBestTrack()->dz(primaryVertex.position())) < 0.5 
                         && mu.innerTrack()->hitPattern().numberOfValidPixelHits() > 0 
                         && mu.innerTrack()->hitPattern().trackerLayersWithMeasurement() > 5 );
}

void MuonIDTableProducer::produce ( edm::StreamID  , edm::Event &  , edm::EventSetup const &    ) const

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 54 of file MuonIDTableProducer.cc.

References nanoaod::FlatTable::BoolColumn, edm::Event::getByToken(), mps_fire::i, muon::isHighPtMuon(), muon::isMediumMuon(), isMediumMuonHIP(), muon::isSoftMuon(), isSoftMuonHIP(), muon::isTightMuon(), isTrackerHighPt(), eostools::move(), RPCpg::mu, extraflags_cff::muons, name_, edm::Event::put(), MetAnalyzer::pv(), edm::Event::run(), src_, srcVtx_, nanoaod::FlatTable::UInt8Column, and electrons_cff::vertices.

Referenced by fillDescriptions().

{

    edm::Handle<std::vector<pat::Muon>> muons;
    iEvent.getByToken(src_, muons);
    edm::Handle<std::vector<reco::Vertex>> vertices;
    iEvent.getByToken(srcVtx_, vertices);

    unsigned int ncand = muons->size();

    const reco::Vertex & pv = vertices->front(); // consistent with IP information in slimmedLeptons.

    bool isRun2016BCDEF = (272007 <= iEvent.run() && iEvent.run() <= 278808);
    std::vector<uint8_t> tight(ncand, 0), highPt(ncand, 0), soft(ncand, 0), medium(ncand, 0);
    for (unsigned int i = 0; i < ncand; ++i) {
        const pat::Muon & mu = (*muons)[i];
        tight[i] = muon::isTightMuon(mu, pv);
        highPt[i] = muon::isHighPtMuon(mu, pv) ? 2 : isTrackerHighPt(mu,pv);
        soft[i] = isRun2016BCDEF ? isSoftMuonHIP(mu,pv) : muon::isSoftMuon(mu, pv);
        medium[i] = isRun2016BCDEF ? isMediumMuonHIP(mu) : muon::isMediumMuon(mu);
    }

    auto tab  = std::make_unique<nanoaod::FlatTable>(ncand, name_, false, true);
    tab->addColumn<uint8_t>("tightId", tight, "POG Tight muon ID", nanoaod::FlatTable::BoolColumn);
    tab->addColumn<uint8_t>("highPtId", highPt, "POG highPt muon ID (1 = tracker high pT, 2 = global high pT, which includes tracker high pT)", nanoaod::FlatTable::UInt8Column);
    tab->addColumn<uint8_t>("softId", soft, "POG Soft muon ID (using the relaxed cuts in the data Run 2016 B-F periods, and standard cuts elsewhere)", nanoaod::FlatTable::BoolColumn);
    tab->addColumn<uint8_t>("mediumId", medium, "POG Medium muon ID (using the relaxed cuts in the data Run 2016 B-F periods, and standard cuts elsewhere)", nanoaod::FlatTable::BoolColumn);

    iEvent.put(std::move(tab));
}

Member Data Documentation

std::string MuonIDTableProducer::name_

private

Definition at line 43 of file MuonIDTableProducer.cc.

Referenced by produce().

edm::EDGetTokenT<std::vector<pat::Muon> > MuonIDTableProducer::src_

private

Definition at line 44 of file MuonIDTableProducer.cc.

Referenced by produce().

edm::EDGetTokenT<std::vector<reco::Vertex> > MuonIDTableProducer::srcVtx_

private

Definition at line 45 of file MuonIDTableProducer.cc.

Referenced by produce().