#include <RecoDisplacedMuonValidator.h>

Inheritance diagram for RecoDisplacedMuonValidator:

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

	RecoDisplacedMuonValidator (const edm::ParameterSet &pset)

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

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

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

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > 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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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_

bool	tpRefVector

edm::EDGetTokenT< TrackingParticleRefVector >	tpRefVectorToken_

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
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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 edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
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)

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

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

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

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

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

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

ProducesCollector	producesCollector ()

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

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

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

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

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

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

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

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 29 of file RecoDisplacedMuonValidator.h.

Constructor & Destructor Documentation

◆ RecoDisplacedMuonValidator()

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

Definition at line 555 of file RecoDisplacedMuonValidator.cc.

     : selector_(pset.getParameter<std::string>("selection")) {
   // dump cfg parameters
   edm::LogVerbatim("RecoDisplacedMuonValidator") << "constructing RecoDisplacedMuonValidator: " << pset.dump();
   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");
   tpRefVector = pset.getParameter<bool>("tpRefVector");
   if (tpRefVector)
     tpRefVectorToken_ = consumes<TrackingParticleRefVector>(simLabel_);
   else
     simToken_ = consumes<TrackingParticleCollection>(simLabel_);
 
   muonLabel_ = pset.getParameter<InputTag>("muonLabel");
   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_ = "RecoDisplacedMuonV";
   if (subDir_[subDir_.size() - 1] == '/')
     subDir_.erase(subDir_.size() - 1);
 }

◆ ~RecoDisplacedMuonValidator()

RecoDisplacedMuonValidator::~RecoDisplacedMuonValidator ( )

override

Definition at line 800 of file RecoDisplacedMuonValidator.cc.

800 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 819 of file RecoDisplacedMuonValidator.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("RecoDisplacedMuonValidator") << "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
   TrackingParticleRefVector TPrefV;
   const TrackingParticleRefVector* ptr_TPrefV = nullptr;
   Handle<TrackingParticleCollection> simHandle;
   Handle<TrackingParticleRefVector> TPCollectionRefVector_H;
 
   if (tpRefVector) {
     event.getByToken(tpRefVectorToken_, TPCollectionRefVector_H);
     ptr_TPrefV = TPCollectionRefVector_H.product();
     TPrefV = *ptr_TPrefV;
   } else {
     event.getByToken(simToken_, simHandle);
     size_t nTP = simHandle->size();
     for (size_t i = 0; i < nTP; ++i) {
       TPrefV.push_back(TrackingParticleRef(simHandle, i));
     }
     ptr_TPrefV = &TPrefV;
   }
 
   // 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 size_t nSim = ptr_TPrefV->size();
 
   edm::RefToBaseVector<reco::Muon> Muons;
   for (size_t i = 0; i < muonHandle->size(); ++i) {
     Muons.push_back(muonHandle->refAt(i));
   }
 
   muonME_->hNSim_->Fill(nSim);
   muonME_->hNMuon_->Fill(muonColl.size());
 
   reco::MuonToSimCollection muonToSimColl;
   reco::SimToMuonCollection simToMuonColl;
 
   if (doAssoc_) {
     edm::LogVerbatim("RecoDisplacedMuonValidator")
         << "\n >>> MuonToSim association : " << muAssocLabel_ << " <<< \n"
         << "     muon collection : " << muonLabel_ << " (size = " << muonHandle->size() << ") \n"
         << "     TrackingParticle collection : " << simLabel_ << " (size = " << nSim << ")";
 
     assoByHits->associateMuons(muonToSimColl, simToMuonColl, Muons, trackType_, TPrefV);
   } 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 (trackType_ == reco::InnerTk) {
       Track = iMuon->track();
       trkNTrackerHits = countTrackerHits(*Track);
     } else if (trackType_ == reco::OuterTk) {
       Track = iMuon->standAloneMuon();
     } else if (trackType_ == reco::GlobalTk) {
       Track = iMuon->combinedMuon();
       glbNTrackerHits = countTrackerHits(*Track);
       glbNMuonHits = countMuonHits(*Track);
     }
 
     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());
 
     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 (size_t i = 0; i < nSim; i++) {
     TrackingParticleRef simRef = TPrefV[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 RecoDisplacedMuonValidator::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMOneEDAnalyzer<>.

Definition at line 692 of file RecoDisplacedMuonValidator.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 RecoDisplacedMuonValidator::countMuonHits ( const reco::Track & track ) const

virtual

Definition at line 1086 of file RecoDisplacedMuonValidator.cc.

References submitPVResolutionJobs::count, DetId::det(), DetId::Muon, mps_fire::result, and HLT_2024v14_cff::track.

Referenced by analyze().

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

◆ countTrackerHits()

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

virtual

Definition at line 1101 of file RecoDisplacedMuonValidator.cc.

References submitPVResolutionJobs::count, DetId::det(), mps_fire::result, HLT_2024v14_cff::track, and DetId::Tracker.

Referenced by analyze().

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

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 806 of file RecoDisplacedMuonValidator.cc.

806 {}

◆ dqmEndRun()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer<>.

Definition at line 811 of file RecoDisplacedMuonValidator.cc.

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

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