#include <RecoMuonValidator.h>

Inheritance diagram for RecoMuonValidator:

Classes
struct	CommonME

struct	HistoDimensions

struct	MuonME

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

virtual int	countMuonHits (const reco::Track &track) const

virtual int	countTrackerHits (const reco::Track &track) const

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &eventSetup) override

void	dqmEndRun (edm::Run const &, edm::EventSetup const &) override

	RecoMuonValidator (const edm::ParameterSet &pset)

	~RecoMuonValidator () override

Public Member Functions inherited from DQMOneEDAnalyzer<>
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) 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

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)

Protected Attributes
edm::InputTag	beamspotLabel_

edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_

CommonME *	commonME_

DQMStore *	dbe_

bool	doAbsEta_

bool	doAssoc_

HistoDimensions	hDim

edm::InputTag	muAssocLabel_

edm::EDGetTokenT< reco::MuonToTrackingParticleAssociator >	muAssocToken_

edm::InputTag	muonLabel_

MuonME *	muonME_

std::string	muonSelection_

edm::EDGetTokenT< edm::View< reco::Muon > >	muonToken_

std::string	outputFileName_

edm::InputTag	primvertexLabel_

edm::EDGetTokenT< reco::VertexCollection >	primvertexToken_

edm::ParameterSet	pset

edm::InputTag	simLabel_

edm::EDGetTokenT< TrackingParticleCollection >	simToken_

std::string	subDir_

std::string	subsystemname_

TrackingParticleSelector	tpSelector_

reco::MuonTrackType	trackType_

bool	usePFMuon_

unsigned int	verbose_

Protected Attributes inherited from DQMOneEDAnalyzer<>
edm::EDPutTokenT< DQMToken >	runToken_

Private Attributes
StringCutObjectSelector< reco::Muon >	selector_

bool	wantTightMuon_

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer<>
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMOneEDAnalyzer<>
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)=0

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

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

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 29 of file RecoMuonValidator.h.

Constructor & Destructor Documentation

◆ RecoMuonValidator()

RecoMuonValidator::RecoMuonValidator ( const edm::ParameterSet & pset )

Definition at line 555 of file RecoMuonValidator.cc.

     : selector_(pset.getParameter<std::string>("selection")) {
   //  pset=ps;
   verbose_ = pset.getUntrackedParameter<unsigned int>("verbose", 0);
  
   outputFileName_ = pset.getUntrackedParameter<string>("outputFileName", "");
  
   wantTightMuon_ = pset.getParameter<bool>("wantTightMuon");
   beamspotLabel_ = pset.getParameter<edm::InputTag>("beamSpot");
   primvertexLabel_ = pset.getParameter<edm::InputTag>("primaryVertex");
   beamspotToken_ = consumes<reco::BeamSpot>(beamspotLabel_);
   primvertexToken_ = consumes<reco::VertexCollection>(primvertexLabel_);
  
   // Set histogram dimensions from config
  
   hDim.nBinP = pset.getUntrackedParameter<unsigned int>("nBinP");
   hDim.minP = pset.getUntrackedParameter<double>("minP");
   hDim.maxP = pset.getUntrackedParameter<double>("maxP");
  
   hDim.nBinPt = pset.getUntrackedParameter<unsigned int>("nBinPt");
   hDim.minPt = pset.getUntrackedParameter<double>("minPt");
   hDim.maxPt = pset.getUntrackedParameter<double>("maxPt");
  
   doAbsEta_ = pset.getUntrackedParameter<bool>("doAbsEta");
   hDim.doAbsEta = doAbsEta_;
   hDim.nBinEta = pset.getUntrackedParameter<unsigned int>("nBinEta");
   hDim.minEta = pset.getUntrackedParameter<double>("minEta");
   hDim.maxEta = pset.getUntrackedParameter<double>("maxEta");
  
   hDim.nBinDxy = pset.getUntrackedParameter<unsigned int>("nBinDxy");
   hDim.minDxy = pset.getUntrackedParameter<double>("minDxy");
   hDim.maxDxy = pset.getUntrackedParameter<double>("maxDxy");
  
   hDim.nBinDz = pset.getUntrackedParameter<unsigned int>("nBinDz");
   hDim.minDz = pset.getUntrackedParameter<double>("minDz");
   hDim.maxDz = pset.getUntrackedParameter<double>("maxDz");
  
   hDim.nBinPhi = pset.getUntrackedParameter<unsigned int>("nBinPhi");
   hDim.minPhi = pset.getUntrackedParameter<double>("minPhi", -TMath::Pi());
   hDim.maxPhi = pset.getUntrackedParameter<double>("maxPhi", TMath::Pi());
  
   hDim.nBinErr = pset.getUntrackedParameter<unsigned int>("nBinErr");
   hDim.nBinPull = pset.getUntrackedParameter<unsigned int>("nBinPull");
  
   hDim.wPull = pset.getUntrackedParameter<double>("wPull");
  
   hDim.minErrP = pset.getUntrackedParameter<double>("minErrP");
   hDim.maxErrP = pset.getUntrackedParameter<double>("maxErrP");
  
   hDim.minErrPt = pset.getUntrackedParameter<double>("minErrPt");
   hDim.maxErrPt = pset.getUntrackedParameter<double>("maxErrPt");
  
   hDim.minErrQPt = pset.getUntrackedParameter<double>("minErrQPt");
   hDim.maxErrQPt = pset.getUntrackedParameter<double>("maxErrQPt");
  
   hDim.minErrEta = pset.getUntrackedParameter<double>("minErrEta");
   hDim.maxErrEta = pset.getUntrackedParameter<double>("maxErrEta");
  
   hDim.minErrPhi = pset.getUntrackedParameter<double>("minErrPhi");
   hDim.maxErrPhi = pset.getUntrackedParameter<double>("maxErrPhi");
  
   hDim.minErrDxy = pset.getUntrackedParameter<double>("minErrDxy");
   hDim.maxErrDxy = pset.getUntrackedParameter<double>("maxErrDxy");
  
   hDim.minErrDz = pset.getUntrackedParameter<double>("minErrDz");
   hDim.maxErrDz = pset.getUntrackedParameter<double>("maxErrDz");
  
   hDim.nTrks = pset.getUntrackedParameter<unsigned int>("nTrks");
   hDim.nAssoc = pset.getUntrackedParameter<unsigned int>("nAssoc");
   hDim.nDof = pset.getUntrackedParameter<unsigned int>("nDof", 55);
  
   // Labels for simulation and reconstruction tracks
   simLabel_ = pset.getParameter<InputTag>("simLabel");
   muonLabel_ = pset.getParameter<InputTag>("muonLabel");
   simToken_ = consumes<TrackingParticleCollection>(simLabel_);
   muonToken_ = consumes<edm::View<reco::Muon> >(muonLabel_);
  
   // Labels for sim-reco association
   doAssoc_ = pset.getUntrackedParameter<bool>("doAssoc", true);
   muAssocLabel_ = pset.getParameter<InputTag>("muAssocLabel");
   if (doAssoc_) {
     muAssocToken_ = consumes<reco::MuonToTrackingParticleAssociator>(muAssocLabel_);
   }
  
   // Different momentum assignment and additional histos in case of PF muons
   usePFMuon_ = pset.getUntrackedParameter<bool>("usePFMuon");
   hDim.usePFMuon = usePFMuon_;
  
   //type of track
   std::string trackType = pset.getParameter<std::string>("trackType");
   if (trackType == "inner")
     trackType_ = reco::InnerTk;
   else if (trackType == "outer")
     trackType_ = reco::OuterTk;
   else if (trackType == "global")
     trackType_ = reco::GlobalTk;
   else if (trackType == "segments")
     trackType_ = reco::Segments;
   else
     throw cms::Exception("Configuration") << "Track type '" << trackType << "' not supported.\n";
  
   //  seedPropagatorName_ = pset.getParameter<string>("SeedPropagator");
  
   ParameterSet tpset = pset.getParameter<ParameterSet>("tpSelector");
   tpSelector_ = TrackingParticleSelector(tpset.getParameter<double>("ptMin"),
                                          tpset.getParameter<double>("ptMax"),
                                          tpset.getParameter<double>("minRapidity"),
                                          tpset.getParameter<double>("maxRapidity"),
                                          tpset.getParameter<double>("tip"),
                                          tpset.getParameter<double>("lip"),
                                          tpset.getParameter<int>("minHit"),
                                          tpset.getParameter<bool>("signalOnly"),
                                          tpset.getParameter<bool>("intimeOnly"),
                                          tpset.getParameter<bool>("chargedOnly"),
                                          tpset.getParameter<bool>("stableOnly"),
                                          tpset.getParameter<std::vector<int> >("pdgId"));
  
   // the service parameters
   ParameterSet serviceParameters = pset.getParameter<ParameterSet>("ServiceParameters");
  
   // retrieve the instance of DQMService
   dbe_ = Service<DQMStore>().operator->();
   subsystemname_ = pset.getUntrackedParameter<std::string>("subSystemFolder", "YourSubsystem");
  
   subDir_ = pset.getUntrackedParameter<string>("subDir");
   if (subDir_.empty())
     subDir_ = "RecoMuonV";
   if (subDir_[subDir_.size() - 1] == '/')
     subDir_.erase(subDir_.size() - 1);
 }

◆ ~RecoMuonValidator()

RecoMuonValidator::~RecoMuonValidator ( )

override

Definition at line 794 of file RecoMuonValidator.cc.

794 {}

Member Function Documentation

◆ analyze()

void RecoMuonValidator::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 813 of file RecoMuonValidator.cc.

                                                                                 {
   // Look for the Primary Vertex (and use the BeamSpot instead, if you can't find it):
   reco::Vertex::Point posVtx;
   reco::Vertex::Error errVtx;
   edm::Handle<reco::VertexCollection> recVtxs;
   event.getByToken(primvertexToken_, recVtxs);
   unsigned int theIndexOfThePrimaryVertex = 999.;
   for (unsigned int ind = 0; ind < recVtxs->size(); ++ind) {
     if ((*recVtxs)[ind].isValid() && !((*recVtxs)[ind].isFake())) {
       theIndexOfThePrimaryVertex = ind;
       break;
     }
   }
   if (theIndexOfThePrimaryVertex < 100) {
     posVtx = ((*recVtxs)[theIndexOfThePrimaryVertex]).position();
     errVtx = ((*recVtxs)[theIndexOfThePrimaryVertex]).error();
   } else {
     LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
     edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
     event.getByToken(beamspotToken_, recoBeamSpotHandle);
     reco::BeamSpot bs = *recoBeamSpotHandle;
     posVtx = bs.position();
     errVtx(0, 0) = bs.BeamWidthX();
     errVtx(1, 1) = bs.BeamWidthY();
     errVtx(2, 2) = bs.sigmaZ();
   }
   const reco::Vertex thePrimaryVertex(posVtx, errVtx);
  
   // Get TrackingParticles
   Handle<TrackingParticleCollection> simHandle;
   event.getByToken(simToken_, simHandle);
   const TrackingParticleCollection simColl = *(simHandle.product());
  
   // Get Muons
   Handle<edm::View<Muon> > muonHandle;
   event.getByToken(muonToken_, muonHandle);
   View<Muon> muonColl = *(muonHandle.product());
  
   reco::MuonToTrackingParticleAssociator const* assoByHits = nullptr;
   if (doAssoc_) {
     edm::Handle<reco::MuonToTrackingParticleAssociator> associatorBase;
     event.getByToken(muAssocToken_, associatorBase);
     assoByHits = associatorBase.product();
   }
  
   const TrackingParticleCollection::size_type nSim = simColl.size();
  
   edm::RefToBaseVector<reco::Muon> Muons;
   for (size_t i = 0; i < muonHandle->size(); ++i) {
     Muons.push_back(muonHandle->refAt(i));
   }
  
   edm::RefVector<TrackingParticleCollection> allTPs;
   for (size_t i = 0; i < nSim; ++i) {
     allTPs.push_back(TrackingParticleRef(simHandle, i));
   }
  
   muonME_->hNSim_->Fill(nSim);
   muonME_->hNMuon_->Fill(muonColl.size());
  
   reco::MuonToSimCollection muonToSimColl;
   reco::SimToMuonCollection simToMuonColl;
  
   if (doAssoc_) {
     assoByHits->associateMuons(muonToSimColl, simToMuonColl, Muons, trackType_, allTPs);
   } else {
     /*
     // SimToMuon associations
     Handle<reco::RecoToSimCollection> simToTrkMuHandle;
     event.getByLabel(trkMuAssocLabel_, simToTrkMuHandle);
     trkSimRecColl = *(simToTrkMuHandle.product());
  
     Handle<reco::RecoToSimCollection> simToStaMuHandle;
     event.getByLabel(staMuAssocLabel_, simToStaMuHandle);
     staSimRecColl = *(simToStaMuHandle.product());
  
     Handle<reco::RecoToSimCollection> simToGlbMuHandle;
     event.getByLabel(glbMuAssocLabel_, simToGlbMuHandle);
     glbSimRecColl = *(simToGlbMuHandle.product());
  
     // MuonToSim associations
     Handle<reco::SimToRecoCollection> trkMuToSimHandle;
     event.getByLabel(trkMuAssocLabel_, trkMuToSimHandle);
     trkRecSimColl = *(trkMuToSimHandle.product());
  
     Handle<reco::SimToRecoCollection> staMuToSimHandle;
     event.getByLabel(staMuAssocLabel_, staMuToSimHandle);
     staRecSimColl = *(staMuToSimHandle.product());
  
     Handle<reco::SimToRecoCollection> glbMuToSimHandle;
     event.getByLabel(glbMuAssocLabel_, glbMuToSimHandle);
     glbRecSimColl = *(glbMuToSimHandle.product());
 */
   }
  
   int glbNTrackerHits = 0;
   int trkNTrackerHits = 0;
   int glbNMuonHits = 0;
   int staNMuonHits = 0;
   int NTrackerHits = 0;
   int NMuonHits = 0;
  
   // Analyzer reco::Muon
   for (View<Muon>::const_iterator iMuon = muonColl.begin(); iMuon != muonColl.end(); ++iMuon) {
     double muonP, muonPt, muonEta, muonPhi;
     if (usePFMuon_) {
       muonP = iMuon->pfP4().P();
       muonPt = iMuon->pfP4().Pt();
       muonEta = iMuon->pfP4().Eta();
       muonPhi = iMuon->pfP4().Phi();
     } else {
       muonP = iMuon->p();
       muonPt = iMuon->pt();
       muonEta = iMuon->eta();
       muonPhi = iMuon->phi();
     }
  
     //histograms for fractions
     commonME_->hMuonAllP_->Fill(muonP);
     commonME_->hMuonAllPt_->Fill(muonPt);
     commonME_->hMuonAllEta_->Fill(muonEta);
     commonME_->hMuonAllPhi_->Fill(muonPhi);
  
     if (!selector_(*iMuon))
       continue;
     if (wantTightMuon_) {
       if (!muon::isTightMuon(*iMuon, thePrimaryVertex))
         continue;
     }
  
     TrackRef Track = iMuon->track();
  
     if (Track.isNonnull()) {
       commonME_->hMuonTrackP_->Fill(Track->p());
       commonME_->hMuonTrackPt_->Fill(Track->pt());
       commonME_->hMuonTrackEta_->Fill(Track->eta());
       commonME_->hMuonTrackPhi_->Fill(Track->phi());
  
       //ip histograms
       commonME_->hMuonTrackDxy_->Fill(Track->dxy());
       commonME_->hMuonTrackDz_->Fill(Track->dz());
     }
  
     if (iMuon->isGlobalMuon()) {
       Track = iMuon->combinedMuon();
       glbNTrackerHits = countTrackerHits(*Track);
       glbNMuonHits = countMuonHits(*Track);
     } else if (iMuon->isTrackerMuon()) {
       Track = iMuon->track();
       trkNTrackerHits = countTrackerHits(*Track);
     } else {
       Track = iMuon->standAloneMuon();
     }
  
     NTrackerHits = countTrackerHits(*Track);
     muonME_->hNTrackerHits_->Fill(NTrackerHits);
     muonME_->hNTrackerHits_vs_Pt_->Fill(Track->pt(), NTrackerHits);
     muonME_->hNTrackerHits_vs_Eta_->Fill(Track->eta(), NTrackerHits);
  
     NMuonHits = countMuonHits(*Track);
     muonME_->hNMuonHits_->Fill(NMuonHits);
     muonME_->hNMuonHits_vs_Pt_->Fill(Track->pt(), NMuonHits);
     muonME_->hNMuonHits_vs_Eta_->Fill(Track->eta(), NMuonHits);
  
     //list of histos for each type
  
     //      muonME_->hNTrks_->Fill();
     muonME_->hNTrksEta_->Fill(Track->eta());
     muonME_->hNTrksPt_->Fill(Track->pt());
  
     commonME_->hMuonP_->Fill(muonP);
     commonME_->hMuonPt_->Fill(muonPt);
     commonME_->hMuonEta_->Fill(muonEta);
     commonME_->hMuonPhi_->Fill(muonPhi);
  
     if (iMuon->isGlobalMuon()) {
       double gtHitPat = iMuon->globalTrack()->hitPattern().numberOfAllHits(HitPattern::TRACK_HITS) -
                         iMuon->globalTrack()->hitPattern().numberOfValidHits();
  
       double itHitPat = iMuon->innerTrack()->hitPattern().numberOfAllHits(HitPattern::TRACK_HITS) -
                         iMuon->innerTrack()->hitPattern().numberOfValidHits();
  
       double otHitPat = iMuon->outerTrack()->hitPattern().numberOfAllHits(HitPattern::TRACK_HITS) -
                         iMuon->outerTrack()->hitPattern().numberOfValidHits();
  
       commonME_->hNInvalidHitsGTHitPattern_->Fill(gtHitPat);
       commonME_->hNInvalidHitsITHitPattern_->Fill(itHitPat);
       commonME_->hNInvalidHitsOTHitPattern_->Fill(otHitPat);
       commonME_->hNDeltaInvalidHitsHitPattern_->Fill(gtHitPat - itHitPat - otHitPat);
  
       //must be global and standalone
       if (iMuon->isStandAloneMuon()) {
         commonME_->hStaToGlbDiffNMuonHitsEta_->Fill(Track->eta(), staNMuonHits - glbNMuonHits);
         commonME_->hStaToGlbDiffNMuonHitsPt_->Fill(Track->pt(), staNMuonHits - glbNMuonHits);
         commonME_->hStaToGlbDiffNMuonHits_->Fill(staNMuonHits - glbNMuonHits);
       }
  
       //must be global and tracker
       if (iMuon->isTrackerMuon()) {
         commonME_->hTrkToGlbDiffNTrackerHitsEta_->Fill(Track->eta(), trkNTrackerHits - glbNTrackerHits);
         commonME_->hTrkToGlbDiffNTrackerHitsPt_->Fill(Track->pt(), trkNTrackerHits - glbNTrackerHits);
         commonME_->hTrkToGlbDiffNTrackerHits_->Fill(trkNTrackerHits - glbNTrackerHits);
       }
     }
  
   }  //end of reco muon loop
  
   // Associate by hits
   for (TrackingParticleCollection::size_type i = 0; i < nSim; i++) {
     TrackingParticleRef simRef(simHandle, i);
     const TrackingParticle* simTP = simRef.get();
     if (!tpSelector_(*simTP))
       continue;
  
     //denominators for efficiency plots
     const double simP = simRef->p();
     const double simPt = simRef->pt();
     const double simEta = doAbsEta_ ? fabs(simRef->eta()) : simRef->eta();
     const double simPhi = simRef->phi();
  
     GlobalPoint simVtx(simRef->vertex().x(), simRef->vertex().y(), simRef->vertex().z());
     GlobalVector simMom(simRef->momentum().x(), simRef->momentum().y(), simRef->momentum().z());
     const double simDxy = -simVtx.x() * sin(simPhi) + simVtx.y() * cos(simPhi);
     const double simDz = simVtx.z() - (simVtx.x() * simMom.x() + simVtx.y() * simMom.y()) * simMom.z() / simMom.perp2();
  
     const unsigned int nSimHits = simRef->numberOfHits();
  
     muonME_->hSimP_->Fill(simP);
     muonME_->hSimPt_->Fill(simPt);
     muonME_->hSimEta_->Fill(simEta);
     muonME_->hSimPhi_->Fill(simPhi);
     muonME_->hSimDxy_->Fill(simDxy);
     muonME_->hSimDz_->Fill(simDz);
     muonME_->hNSimHits_->Fill(nSimHits);
  
     // Get sim-reco association for a simRef
     vector<pair<RefToBase<Muon>, double> > MuRefV;
     if (simToMuonColl.find(simRef) != simToMuonColl.end()) {
       MuRefV = simToMuonColl[simRef];
  
       if (!MuRefV.empty()) {
         muonME_->hNSimToReco_->Fill(MuRefV.size());
         const Muon* Mu = MuRefV.begin()->first.get();
         if (!selector_(*Mu))
           continue;
         if (wantTightMuon_) {
           if (!muon::isTightMuon(*Mu, thePrimaryVertex))
             continue;
         }
  
         muonME_->fill(&*simTP, Mu);
       }
     }
   }
 }

◆ bookHistograms()

void RecoMuonValidator::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

override

Definition at line 686 of file RecoMuonValidator.cc.

                                                                {
   // book histograms
   ibooker.cd();
  
   ibooker.setCurrentFolder(subDir_);
  
   commonME_ = new CommonME;
   muonME_ = new MuonME;
  
   //commonME
   const int nHits = 100;
  
   // - diffs
   commonME_->hTrkToGlbDiffNTrackerHits_ = ibooker.book1D("TrkGlbDiffNTrackerHits",
                                                          "Difference of number of tracker hits (tkMuon - globalMuon)",
                                                          2 * nHits + 1,
                                                          -nHits - 0.5,
                                                          nHits + 0.5);
   commonME_->hStaToGlbDiffNMuonHits_ = ibooker.book1D("StaGlbDiffNMuonHits",
                                                       "Difference of number of muon hits (staMuon - globalMuon)",
                                                       2 * nHits + 1,
                                                       -nHits - 0.5,
                                                       nHits + 0.5);
  
   commonME_->hTrkToGlbDiffNTrackerHitsEta_ =
       ibooker.book2D("TrkGlbDiffNTrackerHitsEta",
                      "Difference of number of tracker hits (tkMuon - globalMuon)",
                      hDim.nBinEta,
                      hDim.minEta,
                      hDim.maxEta,
                      2 * nHits + 1,
                      -nHits - 0.5,
                      nHits + 0.5);
   commonME_->hStaToGlbDiffNMuonHitsEta_ = ibooker.book2D("StaGlbDiffNMuonHitsEta",
                                                          "Difference of number of muon hits (staMuon - globalMuon)",
                                                          hDim.nBinEta,
                                                          hDim.minEta,
                                                          hDim.maxEta,
                                                          2 * nHits + 1,
                                                          -nHits - 0.5,
                                                          nHits + 0.5);
  
   commonME_->hTrkToGlbDiffNTrackerHitsPt_ = ibooker.book2D("TrkGlbDiffNTrackerHitsPt",
                                                            "Difference of number of tracker hits (tkMuon - globalMuon)",
                                                            hDim.nBinPt,
                                                            hDim.minPt,
                                                            hDim.maxPt,
                                                            2 * nHits + 1,
                                                            -nHits - 0.5,
                                                            nHits + 0.5);
   commonME_->hStaToGlbDiffNMuonHitsPt_ = ibooker.book2D("StaGlbDiffNMuonHitsPt",
                                                         "Difference of number of muon hits (staMuon - globalMuon)",
                                                         hDim.nBinPt,
                                                         hDim.minPt,
                                                         hDim.maxPt,
                                                         2 * nHits + 1,
                                                         -nHits - 0.5,
                                                         nHits + 0.5);
  
   // -global muon hit pattern
   commonME_->hNInvalidHitsGTHitPattern_ = ibooker.book1D(
       "NInvalidHitsGTHitPattern", "Number of invalid hits on a global track", nHits + 1, -0.5, nHits + 0.5);
   commonME_->hNInvalidHitsITHitPattern_ = ibooker.book1D(
       "NInvalidHitsITHitPattern", "Number of invalid hits on an inner track", nHits + 1, -0.5, nHits + 0.5);
   commonME_->hNInvalidHitsOTHitPattern_ = ibooker.book1D(
       "NInvalidHitsOTHitPattern", "Number of invalid hits on an outer track", nHits + 1, -0.5, nHits + 0.5);
   commonME_->hNDeltaInvalidHitsHitPattern_ =
       ibooker.book1D("hNDeltaInvalidHitsHitPattern",
                      "The discrepancy for Number of invalid hits on an global track and inner and outer tracks",
                      2 * nHits + 1,
                      -nHits - 0.5,
                      nHits + 0.5);
  
   //muon based kinematics
   commonME_->hMuonP_ = ibooker.book1D("PMuon", "p of muon", hDim.nBinP, hDim.minP, hDim.maxP);
   commonME_->hMuonPt_ = ibooker.book1D("PtMuon", "p_{T} of muon", hDim.nBinPt, hDim.minPt, hDim.maxPt);
   commonME_->hMuonEta_ = ibooker.book1D("EtaMuon", "#eta of muon", hDim.nBinEta, hDim.minEta, hDim.maxEta);
   commonME_->hMuonPhi_ = ibooker.book1D("PhiMuon", "#phi of muon", hDim.nBinPhi, hDim.minPhi, hDim.maxPhi);
   //track based kinematics
   commonME_->hMuonTrackP_ = ibooker.book1D("PMuonTrack", "p of reco muon track", hDim.nBinP, hDim.minP, hDim.maxP);
   commonME_->hMuonTrackPt_ =
       ibooker.book1D("PtMuonTrack", "p_{T} of reco muon track", hDim.nBinPt, hDim.minPt, hDim.maxPt);
   commonME_->hMuonTrackEta_ =
       ibooker.book1D("EtaMuonTrack", "#eta of reco muon track", hDim.nBinEta, hDim.minEta, hDim.maxEta);
   commonME_->hMuonTrackPhi_ =
       ibooker.book1D("PhiMuonTrack", "#phi of reco muon track", hDim.nBinPhi, hDim.minPhi, hDim.maxPhi);
   commonME_->hMuonTrackDxy_ =
       ibooker.book1D("DxyMuonTrack", "Dxy of reco muon track", hDim.nBinDxy, hDim.minDxy, hDim.maxDxy);
   commonME_->hMuonTrackDz_ =
       ibooker.book1D("DzMuonTrack", "Dz of reco muon track", hDim.nBinDz, hDim.minDz, hDim.maxDz);
  
   //histograms for fractions
   commonME_->hMuonAllP_ = ibooker.book1D("PMuonAll", "p of muons of all types", hDim.nBinP, hDim.minP, hDim.maxP);
   commonME_->hMuonAllPt_ =
       ibooker.book1D("PtMuonAll", "p_{T} of muon of all types", hDim.nBinPt, hDim.minPt, hDim.maxPt);
   commonME_->hMuonAllEta_ =
       ibooker.book1D("EtaMuonAll", "#eta of muon of all types", hDim.nBinEta, hDim.minEta, hDim.maxEta);
   commonME_->hMuonAllPhi_ =
       ibooker.book1D("PhiMuonAll", "#phi of muon of all types", hDim.nBinPhi, hDim.minPhi, hDim.maxPhi);
  
   muonME_->bookHistos(ibooker, subDir_, hDim);
 }

◆ countMuonHits()

int RecoMuonValidator::countMuonHits ( const reco::Track & track ) const

virtual

Definition at line 1069 of file RecoMuonValidator.cc.

                                                                  {
   TransientTrackingRecHit::ConstRecHitContainer result;
  
   int count = 0;
  
   for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
     if ((*hit)->isValid()) {
       DetId recoid = (*hit)->geographicalId();
       if (recoid.det() == DetId::Muon)
         count++;
     }
   }
   return count;
 }

References KineDebug3::count(), DetId::det(), DetId::Muon, mps_fire::result, and HLT_2018_cff::track.

Referenced by analyze().

◆ countTrackerHits()

int RecoMuonValidator::countTrackerHits ( const reco::Track & track ) const

virtual

Definition at line 1084 of file RecoMuonValidator.cc.

                                                                     {
   TransientTrackingRecHit::ConstRecHitContainer result;
  
   int count = 0;
  
   for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
     if ((*hit)->isValid()) {
       DetId recoid = (*hit)->geographicalId();
       if (recoid.det() == DetId::Tracker)
         count++;
     }
   }
   return count;
 }

References KineDebug3::count(), DetId::det(), mps_fire::result, HLT_2018_cff::track, and DetId::Tracker.

Referenced by analyze().

◆ dqmBeginRun()

void RecoMuonValidator::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &	eventSetup
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 800 of file RecoMuonValidator.cc.

800 {}

◆ dqmEndRun()

void RecoMuonValidator::dqmEndRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 805 of file RecoMuonValidator.cc.

                                                                      {
   if (dbe_ && !outputFileName_.empty())
     dbe_->save(outputFileName_);
 }

References dbe_, outputFileName_, and dqm::implementation::DQMStore::save().

Member Data Documentation

◆ beamspotLabel_

edm::InputTag RecoMuonValidator::beamspotLabel_

protected

Definition at line 53 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ beamspotToken_

edm::EDGetTokenT<reco::BeamSpot> RecoMuonValidator::beamspotToken_

protected

Definition at line 55 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ commonME_

CommonME* RecoMuonValidator::commonME_

protected

Definition at line 77 of file RecoMuonValidator.h.

Referenced by analyze(), and bookHistograms().

◆ dbe_

DQMStore* RecoMuonValidator::dbe_

protected

Definition at line 63 of file RecoMuonValidator.h.

Referenced by dqmEndRun(), and RecoMuonValidator().

◆ doAbsEta_

bool RecoMuonValidator::doAbsEta_

protected

Definition at line 65 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ doAssoc_

bool RecoMuonValidator::doAssoc_

protected

Definition at line 66 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ hDim

HistoDimensions RecoMuonValidator::hDim

protected

Definition at line 123 of file RecoMuonValidator.h.

Referenced by bookHistograms(), and RecoMuonValidator().

◆ muAssocLabel_

edm::InputTag RecoMuonValidator::muAssocLabel_

protected

Definition at line 50 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ muAssocToken_

edm::EDGetTokenT<reco::MuonToTrackingParticleAssociator> RecoMuonValidator::muAssocToken_

protected

Definition at line 51 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ muonLabel_

edm::InputTag RecoMuonValidator::muonLabel_

protected

Definition at line 45 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ muonME_

MuonME* RecoMuonValidator::muonME_

protected

Definition at line 74 of file RecoMuonValidator.h.

Referenced by analyze(), and bookHistograms().

◆ muonSelection_

std::string RecoMuonValidator::muonSelection_

protected

Definition at line 46 of file RecoMuonValidator.h.

◆ muonToken_

edm::EDGetTokenT<edm::View<reco::Muon> > RecoMuonValidator::muonToken_

protected

Definition at line 48 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ outputFileName_

std::string RecoMuonValidator::outputFileName_

protected

Definition at line 58 of file RecoMuonValidator.h.

Referenced by dqmEndRun(), and RecoMuonValidator().

◆ primvertexLabel_

edm::InputTag RecoMuonValidator::primvertexLabel_

protected

Definition at line 54 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ primvertexToken_

edm::EDGetTokenT<reco::VertexCollection> RecoMuonValidator::primvertexToken_

protected

Definition at line 56 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ pset

edm::ParameterSet RecoMuonValidator::pset

protected

Definition at line 61 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ selector_

StringCutObjectSelector<reco::Muon> RecoMuonValidator::selector_

private

Definition at line 126 of file RecoMuonValidator.h.

Referenced by analyze().

◆ simLabel_

edm::InputTag RecoMuonValidator::simLabel_

protected

Definition at line 44 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ simToken_

edm::EDGetTokenT<TrackingParticleCollection> RecoMuonValidator::simToken_

protected

Definition at line 47 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ subDir_

std::string RecoMuonValidator::subDir_

protected

Definition at line 59 of file RecoMuonValidator.h.

Referenced by bookHistograms(), and RecoMuonValidator().

◆ subsystemname_

std::string RecoMuonValidator::subsystemname_

protected

Definition at line 60 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ tpSelector_

TrackingParticleSelector RecoMuonValidator::tpSelector_

protected

Definition at line 69 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ trackType_

reco::MuonTrackType RecoMuonValidator::trackType_

protected

Definition at line 72 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ usePFMuon_

bool RecoMuonValidator::usePFMuon_

protected

Definition at line 67 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

◆ verbose_

unsigned int RecoMuonValidator::verbose_

protected

Definition at line 42 of file RecoMuonValidator.h.

Referenced by RecoMuonValidator().

◆ wantTightMuon_

bool RecoMuonValidator::wantTightMuon_

private

Definition at line 127 of file RecoMuonValidator.h.

Referenced by analyze(), and RecoMuonValidator().

Classes

Public Member Functions

Protected Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RecoMuonValidator()

◆ ~RecoMuonValidator()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ countMuonHits()

◆ countTrackerHits()

◆ dqmBeginRun()

◆ dqmEndRun()

Member Data Documentation

◆ beamspotLabel_

◆ beamspotToken_

◆ commonME_

◆ dbe_

◆ doAbsEta_

◆ doAssoc_

◆ hDim

◆ muAssocLabel_

◆ muAssocToken_

◆ muonLabel_

◆ muonME_

◆ muonSelection_

◆ muonToken_

◆ outputFileName_

◆ primvertexLabel_

◆ primvertexToken_

◆ pset

◆ selector_

◆ simLabel_

◆ simToken_

◆ subDir_

◆ subsystemname_

◆ tpSelector_

◆ trackType_

◆ usePFMuon_

◆ verbose_

◆ wantTightMuon_