#include <DiMuonVertexMonitor.h>

Inheritance diagram for DiMuonVertexMonitor:

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

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

	DiMuonVertexMonitor (const edm::ParameterSet &)

	~DiMuonVertexMonitor () override=default

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

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

void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Private Member Functions
const reco::Vertex *	findClosestVertex (const TransientVertex aTransVtx, const reco::VertexCollection *vertices) const

Private Attributes
MonitorElement *	hCosPhi3D_

MonitorElement *	hCosPhi_

MonitorElement *	hCosPhiInv3D_

MonitorElement *	hCosPhiInv_

MonitorElement *	hCutFlow_

MonitorElement *	hdxy_

MonitorElement *	hdxyErr_

MonitorElement *	hdz_

MonitorElement *	hdzErr_

MonitorElement *	hInvMass_

MonitorElement *	hIP2d_

MonitorElement *	hIP2dsig_

MonitorElement *	hIP3d_

MonitorElement *	hIP3dsig_

MonitorElement *	hSVChi2_

MonitorElement *	hSVCompatibility3D_

MonitorElement *	hSVCompatibility_

MonitorElement *	hSVDist3D_

MonitorElement *	hSVDist3DErr_

MonitorElement *	hSVDist3DSig_

MonitorElement *	hSVDist_

MonitorElement *	hSVDistErr_

MonitorElement *	hSVDistSig_

MonitorElement *	hSVNormChi2_

MonitorElement *	hSVProb_

std::pair< float, float >	massLimits_

const float	maxSVdist_

const std::string	MEFolderName_

const std::string	motherName_

const edm::EDGetTokenT< reco::TrackCollection >	tracksToken_

const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord >	ttbESToken_

const bool	useClosestVertex_

const edm::EDGetTokenT< reco::VertexCollection >	vertexToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQM/TrackerMonitorTrack/src/DiMuonVertexMonitor.cc Monitoring quantities related to the DiMuon vertex during Tracker Alignment AlCaReco Production

Definition at line 32 of file DiMuonVertexMonitor.h.

Constructor & Destructor Documentation

◆ DiMuonVertexMonitor()

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

explicit

Definition at line 29 of file DiMuonVertexMonitor.cc.

References massLimits_, and motherName_.

     : ttbESToken_(esConsumes(edm::ESInputTag("", "TransientTrackBuilder"))),
       tracksToken_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("muonTracks"))),
       vertexToken_(consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertices"))),
       motherName_(iConfig.getParameter<std::string>("decayMotherName")),
       MEFolderName_(iConfig.getParameter<std::string>("FolderName")),
       useClosestVertex_(iConfig.getParameter<bool>("useClosestVertex")),
       maxSVdist_(iConfig.getParameter<double>("maxSVdist")) {
   if (motherName_.find('Z') != std::string::npos) {
     massLimits_ = std::make_pair(50., 120);
   } else if (motherName_.find("J/#psi") != std::string::npos) {
     massLimits_ = std::make_pair(2.7, 3.4);
   } else if (motherName_.find("#Upsilon") != std::string::npos) {
     massLimits_ = std::make_pair(8.9, 9.9);
   } else {
     edm::LogError("DiMuonVertexMonitor") << " unrecognized decay mother particle: " << motherName_
                                          << " setting the default for the Z->mm (50.,120.)" << std::endl;
     massLimits_ = std::make_pair(50., 120);
   }
 }

◆ ~DiMuonVertexMonitor()

DiMuonVertexMonitor::~DiMuonVertexMonitor ( )

overridedefault

Member Function Documentation

◆ analyze()

void DiMuonVertexMonitor::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 106 of file DiMuonVertexMonitor.cc.

References TransientTrackBuilder::build(), cmToum, VertexDistanceXY::compatibility(), VertexDistance3D::compatibility(), TransientVertex::degreesOfFreedom(), DeadROC_duringRun::dir, VertexDistanceXY::distance(), VertexDistance3D::distance(), HLT_2023v12_cff::distance, Measurement1D::error(), dqm::impl::MonitorElement::Fill(), findClosestVertex(), edm::EventSetup::getData(), hCosPhi3D_, hCosPhi_, hCosPhiInv3D_, hCosPhiInv_, hdxy_, hdxyErr_, hdz_, hdzErr_, hInvMass_, hIP2d_, hIP2dsig_, hIP3d_, hIP3dsig_, hSVChi2_, hSVCompatibility3D_, hSVCompatibility_, hSVDist3D_, hSVDist3DErr_, hSVDist3DSig_, hSVDist_, hSVDistErr_, hSVDistSig_, hSVNormChi2_, hSVProb_, iEvent, createfilelist::int, electrons_cff::ip3d, edm::HandleBase::isValid(), reco::Vertex::isValid(), TransientVertex::isValid(), LogTrace, maxSVdist_, mumass2, reco::Vertex::position(), TransientVertex::position(), edm::Handle< T >::product(), IPTools::signedImpactParameter3D(), IPTools::signedTransverseImpactParameter(), mathSSE::sqrt(), FrontierCondition_GT_autoExpress_cfi::t0, RandomServiceHelper::t1, TransientVertex::totalChiSquared(), HLT_2023v12_cff::track, tracksToken_, ttbESToken_, useClosestVertex_, Measurement1D::value(), KalmanVertexFitter::vertex(), vertexToken_, AlignmentTracksFromVertexSelector_cfi::vertices, cms::Exception::what(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                      {
   std::vector<const reco::Track*> myTracks;
   const auto trackHandle = iEvent.getHandle(tracksToken_);
   if (!trackHandle.isValid()) {
     edm::LogError("DiMuonVertexMonitor") << "invalid track collection encountered!";
     return;
   }
 
   for (const auto& muonTrk : *trackHandle) {
     myTracks.emplace_back(&muonTrk);
   }
 
   const TransientTrackBuilder* theB = &iSetup.getData(ttbESToken_);
   TransientVertex mumuTransientVtx;
   std::vector<reco::TransientTrack> tks;
 
   if (myTracks.size() != 2) {
     edm::LogWarning("DiMuonVertexMonitor") << "There are not enough tracks to monitor!";
     return;
   }
 
   const auto& t1 = myTracks[1]->momentum();
   const auto& t0 = myTracks[0]->momentum();
   const auto& ditrack = t1 + t0;
 
   const auto& tplus = myTracks[0]->charge() > 0 ? myTracks[0] : myTracks[1];
   const auto& tminus = myTracks[0]->charge() < 0 ? myTracks[0] : myTracks[1];
 
   TLorentzVector p4_tplus(tplus->px(), tplus->py(), tplus->pz(), sqrt((tplus->p() * tplus->p()) + mumass2));
   TLorentzVector p4_tminus(tminus->px(), tminus->py(), tminus->pz(), sqrt((tminus->p() * tminus->p()) + mumass2));
 
   const auto& Zp4 = p4_tplus + p4_tminus;
   float track_invMass = Zp4.M();
   hInvMass_->Fill(track_invMass);
 
   // creat the pair of TLorentVectors used to make the plos
   std::pair<TLorentzVector, TLorentzVector> tktk_p4 = std::make_pair(p4_tplus, p4_tminus);
 
   math::XYZPoint ZpT(ditrack.x(), ditrack.y(), 0);
   math::XYZPoint Zp(ditrack.x(), ditrack.y(), ditrack.z());
 
   for (const auto& track : myTracks) {
     reco::TransientTrack trajectory = theB->build(track);
     tks.push_back(trajectory);
   }
 
   KalmanVertexFitter kalman(true);
   mumuTransientVtx = kalman.vertex(tks);
 
   double SVProb = TMath::Prob(mumuTransientVtx.totalChiSquared(), (int)mumuTransientVtx.degreesOfFreedom());
   hSVProb_->Fill(SVProb);
   hSVChi2_->Fill(mumuTransientVtx.totalChiSquared());
   hSVNormChi2_->Fill(mumuTransientVtx.totalChiSquared() / (int)mumuTransientVtx.degreesOfFreedom());
 
   if (!mumuTransientVtx.isValid())
     return;
 
   const reco::Vertex* theClosestVertex;
   // get collection of reconstructed vertices from event
   edm::Handle<reco::VertexCollection> vertexHandle = iEvent.getHandle(vertexToken_);
   if (vertexHandle.isValid()) {
     const reco::VertexCollection* vertices = vertexHandle.product();
     theClosestVertex = this->findClosestVertex(mumuTransientVtx, vertices);
   } else {
     edm::LogWarning("DiMuonVertexMonitor") << "invalid vertex collection encountered Skipping event!";
     return;
   }
 
   reco::Vertex theMainVtx;
   if (!useClosestVertex_ || theClosestVertex == nullptr) {
     theMainVtx = *theClosestVertex;
   } else {
     theMainVtx = vertexHandle.product()->front();
   }
 
   const math::XYZPoint theMainVtxPos(theMainVtx.position().x(), theMainVtx.position().y(), theMainVtx.position().z());
   const math::XYZPoint myVertex(
       mumuTransientVtx.position().x(), mumuTransientVtx.position().y(), mumuTransientVtx.position().z());
   const math::XYZPoint deltaVtx(
       theMainVtxPos.x() - myVertex.x(), theMainVtxPos.y() - myVertex.y(), theMainVtxPos.z() - myVertex.z());
 
   if (theMainVtx.isValid()) {
     // fill the impact parameter plots
     for (const auto& track : myTracks) {
       hdxy_->Fill(track->dxy(theMainVtxPos) * cmToum);
       hdz_->Fill(track->dz(theMainVtxPos) * cmToum);
       hdxyErr_->Fill(track->dxyError() * cmToum);
       hdzErr_->Fill(track->dzError() * cmToum);
 
       const auto& ttrk = theB->build(track);
       Global3DVector dir(track->px(), track->py(), track->pz());
       const auto& ip2d = IPTools::signedTransverseImpactParameter(ttrk, dir, theMainVtx);
       const auto& ip3d = IPTools::signedImpactParameter3D(ttrk, dir, theMainVtx);
 
       hIP2d_->Fill(ip2d.second.value() * cmToum);
       hIP3d_->Fill(ip3d.second.value() * cmToum);
       hIP2dsig_->Fill(ip2d.second.significance());
       hIP3dsig_->Fill(ip3d.second.significance());
     }
 
     // Z Vertex distance in the xy plane
     VertexDistanceXY vertTool;
     double distance = vertTool.distance(mumuTransientVtx, theMainVtx).value();
     double dist_err = vertTool.distance(mumuTransientVtx, theMainVtx).error();
     float compatibility = 0.;
 
     try {
       compatibility = vertTool.compatibility(mumuTransientVtx, theMainVtx);
     } catch (cms::Exception& er) {
       LogTrace("DiMuonVertexMonitor") << "caught std::exception " << er.what() << std::endl;
     }
 
     hSVCompatibility_->Fill(compatibility);
     hSVDist_->Fill(distance * cmToum);
     hSVDistErr_->Fill(dist_err * cmToum);
     hSVDistSig_->Fill(distance / dist_err);
 
     // Z Vertex distance in 3D
     VertexDistance3D vertTool3D;
     double distance3D = vertTool3D.distance(mumuTransientVtx, theMainVtx).value();
     double dist3D_err = vertTool3D.distance(mumuTransientVtx, theMainVtx).error();
     float compatibility3D = 0.;
 
     try {
       compatibility3D = vertTool3D.compatibility(mumuTransientVtx, theMainVtx);
     } catch (cms::Exception& er) {
       LogTrace("DiMuonVertexMonitor") << "caught std::exception " << er.what() << std::endl;
     }
 
     hSVCompatibility3D_->Fill(compatibility3D);
     hSVDist3D_->Fill(distance3D * cmToum);
     hSVDist3DErr_->Fill(dist3D_err * cmToum);
     hSVDist3DSig_->Fill(distance3D / dist3D_err);
 
     // cut on the PV - SV distance
     if (distance * cmToum < maxSVdist_) {
       double cosphi = (ZpT.x() * deltaVtx.x() + ZpT.y() * deltaVtx.y()) /
                       (sqrt(ZpT.x() * ZpT.x() + ZpT.y() * ZpT.y()) *
                        sqrt(deltaVtx.x() * deltaVtx.x() + deltaVtx.y() * deltaVtx.y()));
 
       double cosphi3D = (Zp.x() * deltaVtx.x() + Zp.y() * deltaVtx.y() + Zp.z() * deltaVtx.z()) /
                         (sqrt(Zp.x() * Zp.x() + Zp.y() * Zp.y() + Zp.z() * Zp.z()) *
                          sqrt(deltaVtx.x() * deltaVtx.x() + deltaVtx.y() * deltaVtx.y() + deltaVtx.z() * deltaVtx.z()));
 
       hCosPhi_->Fill(cosphi);
       hCosPhi3D_->Fill(cosphi3D);
       // inverted
       hCosPhiInv_->Fill(-cosphi);
       hCosPhiInv3D_->Fill(-cosphi3D);
     }
   } else {
     edm::LogWarning("DiMuonVertexMonitor") << "hardest primary vertex in the event is not valid!";
   }
 }

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 50 of file DiMuonVertexMonitor.cc.

References dqm::implementation::IBooker::book1D(), hCosPhi3D_, hCosPhi_, hCosPhiInv3D_, hCosPhiInv_, hdxy_, hdxyErr_, hdz_, hdzErr_, hInvMass_, hIP2d_, hIP2dsig_, hIP3d_, hIP3dsig_, hSVChi2_, hSVCompatibility3D_, hSVCompatibility_, hSVDist3D_, hSVDist3DErr_, hSVDist3DSig_, hSVDist_, hSVDistErr_, hSVDistSig_, hSVNormChi2_, hSVProb_, massLimits_, MEFolderName_, motherName_, dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                       {
   iBooker.setCurrentFolder(MEFolderName_ + "/DiMuonVertexMonitor");
 
   // clang-format off
   std::string ts = fmt::sprintf(";%s vertex probability;N(#mu#mu pairs)", motherName_);
   std::string ps = "N(#mu#mu pairs)";
   hSVProb_ = iBooker.book1D("VtxProb", ts, 100, 0., 1.);
 
   ts = fmt::sprintf("#chi^{2} of the %s vertex; #chi^{2} of the %s vertex; %s", motherName_, motherName_, ps);
   hSVChi2_ = iBooker.book1D("VtxChi2", ts, 200, 0., 200.);
 
   ts = fmt::sprintf("#chi^{2}/ndf of the %s vertex; #chi^{2}/ndf of %s vertex; %s", motherName_, motherName_, ps);
   hSVNormChi2_ = iBooker.book1D("VtxNormChi2", ts, 100, 0., 20.);
 
   std::string histTit = motherName_ + " #rightarrow #mu^{+}#mu^{-}";
   ts = fmt::sprintf("%s;PV- %sV xy distance [#mum];%s", histTit, motherName_, ps);
   hSVDist_ = iBooker.book1D("VtxDist", ts, 100, 0., 300.);
 
   ts = fmt::sprintf("%s;PV-%sV xy distance error [#mum];%s", histTit, motherName_, ps);
   hSVDistErr_ = iBooker.book1D("VtxDistErr", ts, 100, 0., 1000.);
 
   ts = fmt::sprintf("%s;PV-%sV xy distance signficance;%s", histTit, motherName_, ps);
   hSVDistSig_ = iBooker.book1D("VtxDistSig", ts, 100, 0., 5.);
 
   ts = fmt::sprintf("compatibility of %s vertex; compatibility of %s vertex; %s", motherName_, motherName_, ps);
   hSVCompatibility_ = iBooker.book1D("VtxCompatibility", ts, 100, 0., 100.);
 
   ts = fmt::sprintf("%s;PV-%sV 3D distance [#mum];%s", histTit, motherName_, ps);
   hSVDist3D_ = iBooker.book1D("VtxDist3D", ts, 100, 0., 300.);
 
   ts = fmt::sprintf("%s;PV-%sV 3D distance error [#mum];%s", histTit, motherName_, ps);
   hSVDist3DErr_ = iBooker.book1D("VtxDist3DErr", ts, 100, 0., 1000.);
 
   ts = fmt::sprintf("%s;PV-%sV 3D distance signficance;%s", histTit, motherName_, ps);
   hSVDist3DSig_ = iBooker.book1D("VtxDist3DSig", ts, 100, 0., 5.);
 
   ts = fmt::sprintf("3D compatibility of %s vertex;3D compatibility of %s vertex; %s", motherName_, motherName_, ps);
   hSVCompatibility3D_ = iBooker.book1D("VtxCompatibility3D", ts, 100, 0., 100.);
 
   hInvMass_ = iBooker.book1D("InvMass", fmt::sprintf("%s;M(#mu,#mu) [GeV];%s", histTit, ps), 70., massLimits_.first, massLimits_.second);
   hCosPhi_ = iBooker.book1D("CosPhi", fmt::sprintf("%s;cos(#phi_{xy});%s", histTit, ps), 50, -1., 1.);
   hCosPhi3D_ = iBooker.book1D("CosPhi3D", fmt::sprintf("%s;cos(#phi_{3D});%s", histTit, ps), 50, -1., 1.);
   hCosPhiInv_ = iBooker.book1D("CosPhiInv", fmt::sprintf("%s;inverted cos(#phi_{xy});%s", histTit, ps), 50, -1., 1.);
   hCosPhiInv3D_ = iBooker.book1D("CosPhiInv3D", fmt::sprintf("%s;inverted cos(#phi_{3D});%s", histTit, ps), 50, -1., 1.);
 
   hdxy_ = iBooker.book1D("dxy", fmt::sprintf("%s;muon track d_{xy}(PV) [#mum];muon tracks", histTit), 150, -300, 300);
   hdz_ = iBooker.book1D("dz", fmt::sprintf("%s;muon track d_{z}(PV) [#mum];muon tracks", histTit), 150, -300, 300);
   hdxyErr_ = iBooker.book1D("dxyErr", fmt::sprintf("%s;muon track err_{dxy} [#mum];muon tracks", histTit), 250, 0., 500.);
   hdzErr_ = iBooker.book1D("dzErr", fmt::sprintf("%s;muon track err_{dz} [#mum];muon tracks", histTit), 250, 0., 500.);
   hIP2d_ = iBooker.book1D("IP2d", fmt::sprintf("%s;muon track IP_{2D} [#mum];muon tracks", histTit), 150, -300, 300);
   hIP3d_ = iBooker.book1D("IP3d", fmt::sprintf("%s;muon track IP_{3D} [#mum];muon tracks", histTit), 150, -300, 300);
   hIP2dsig_ = iBooker.book1D("IP2Dsig", fmt::sprintf("%s;muon track IP_{2D} significance;muon tracks", histTit), 100, 0., 5.);
   hIP3dsig_ = iBooker.book1D("IP3Dsig", fmt::sprintf("%s;muon track IP_{3D} significance;muon tracks", histTit), 100, 0., 5.);
   // clang-format on
 }

◆ fillDescriptions()

void DiMuonVertexMonitor::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 290 of file DiMuonVertexMonitor.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                      {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("muonTracks", edm::InputTag("ALCARECOTkAlDiMuon"));
   desc.add<edm::InputTag>("vertices", edm::InputTag("offlinePrimaryVertices"));
   desc.add<std::string>("FolderName", "DiMuonVertexMonitor");
   desc.add<std::string>("decayMotherName", "Z");
   desc.add<bool>("useClosestVertex", true);
   desc.add<double>("maxSVdist", 50.);
   descriptions.addWithDefaultLabel(desc);
 }

◆ findClosestVertex()

const reco::Vertex * DiMuonVertexMonitor::findClosestVertex	(	const TransientVertex	aTransVtx,
		const reco::VertexCollection *	vertices
	)		const

private

Definition at line 262 of file DiMuonVertexMonitor.cc.

References counter, VertexDistance3D::distance(), TransientVertex::isValid(), Measurement1D::value(), AlignmentTracksFromVertexSelector_cfi::vertices, and extraflags_cff::vtx.

Referenced by analyze().

                                                                                                        {
   reco::Vertex* defaultVtx = nullptr;
 
   if (!aTransVtx.isValid())
     return defaultVtx;
 
   // find the closest vertex to the secondary vertex in 3D
   VertexDistance3D vertTool3D;
   float minD = 9999.;
   int closestVtxIndex = 0;
   int counter = 0;
   for (const auto& vtx : *vertices) {
     double dist3D = vertTool3D.distance(aTransVtx, vtx).value();
     if (dist3D < minD) {
       minD = dist3D;
       closestVtxIndex = counter;
     }
     counter++;
   }
 
   if ((*vertices).at(closestVtxIndex).isValid()) {
     return &(vertices->at(closestVtxIndex));
   } else {
     return defaultVtx;
   }
 }

Member Data Documentation

◆ hCosPhi3D_

MonitorElement* DiMuonVertexMonitor::hCosPhi3D_

private

Definition at line 76 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hCosPhi_

MonitorElement* DiMuonVertexMonitor::hCosPhi_

private

Definition at line 75 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hCosPhiInv3D_

MonitorElement* DiMuonVertexMonitor::hCosPhiInv3D_

private

Definition at line 78 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hCosPhiInv_

MonitorElement* DiMuonVertexMonitor::hCosPhiInv_

private

Definition at line 77 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hCutFlow_

MonitorElement* DiMuonVertexMonitor::hCutFlow_

private

Definition at line 80 of file DiMuonVertexMonitor.h.

◆ hdxy_

MonitorElement* DiMuonVertexMonitor::hdxy_

private

Definition at line 83 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hdxyErr_

MonitorElement* DiMuonVertexMonitor::hdxyErr_

private

Definition at line 85 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hdz_

MonitorElement* DiMuonVertexMonitor::hdz_

private

Definition at line 84 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hdzErr_

MonitorElement* DiMuonVertexMonitor::hdzErr_

private

Definition at line 86 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hInvMass_

MonitorElement* DiMuonVertexMonitor::hInvMass_

private

Definition at line 79 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hIP2d_

MonitorElement* DiMuonVertexMonitor::hIP2d_

private

Definition at line 87 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hIP2dsig_

MonitorElement* DiMuonVertexMonitor::hIP2dsig_

private

Definition at line 89 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hIP3d_

MonitorElement* DiMuonVertexMonitor::hIP3d_

private

Definition at line 88 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hIP3dsig_

MonitorElement* DiMuonVertexMonitor::hIP3dsig_

private

Definition at line 90 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVChi2_

MonitorElement* DiMuonVertexMonitor::hSVChi2_

private

Definition at line 62 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVCompatibility3D_

MonitorElement* DiMuonVertexMonitor::hSVCompatibility3D_

private

Definition at line 73 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVCompatibility_

MonitorElement* DiMuonVertexMonitor::hSVCompatibility_

private

Definition at line 68 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVDist3D_

MonitorElement* DiMuonVertexMonitor::hSVDist3D_

private

Definition at line 70 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVDist3DErr_

MonitorElement* DiMuonVertexMonitor::hSVDist3DErr_

private

Definition at line 71 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVDist3DSig_

MonitorElement* DiMuonVertexMonitor::hSVDist3DSig_

private

Definition at line 72 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVDist_

MonitorElement* DiMuonVertexMonitor::hSVDist_

private

Definition at line 65 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVDistErr_

MonitorElement* DiMuonVertexMonitor::hSVDistErr_

private

Definition at line 66 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVDistSig_

MonitorElement* DiMuonVertexMonitor::hSVDistSig_

private

Definition at line 67 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVNormChi2_

MonitorElement* DiMuonVertexMonitor::hSVNormChi2_

private

Definition at line 63 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ hSVProb_

MonitorElement* DiMuonVertexMonitor::hSVProb_

private

Definition at line 61 of file DiMuonVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

◆ massLimits_

std::pair<float, float> DiMuonVertexMonitor::massLimits_

private

Definition at line 57 of file DiMuonVertexMonitor.h.

Referenced by bookHistograms(), and DiMuonVertexMonitor().

◆ maxSVdist_

const float DiMuonVertexMonitor::maxSVdist_

private

Definition at line 58 of file DiMuonVertexMonitor.h.

Referenced by analyze().

◆ MEFolderName_

const std::string DiMuonVertexMonitor::MEFolderName_

private

Definition at line 54 of file DiMuonVertexMonitor.h.

Referenced by bookHistograms().

◆ motherName_

const std::string DiMuonVertexMonitor::motherName_

private

Definition at line 53 of file DiMuonVertexMonitor.h.

Referenced by bookHistograms(), and DiMuonVertexMonitor().

◆ tracksToken_

const edm::EDGetTokenT<reco::TrackCollection> DiMuonVertexMonitor::tracksToken_

private

Definition at line 49 of file DiMuonVertexMonitor.h.

Referenced by analyze().

◆ ttbESToken_

const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> DiMuonVertexMonitor::ttbESToken_

private

Definition at line 46 of file DiMuonVertexMonitor.h.

Referenced by analyze().

◆ useClosestVertex_

const bool DiMuonVertexMonitor::useClosestVertex_

private

Definition at line 55 of file DiMuonVertexMonitor.h.

Referenced by analyze().

◆ vertexToken_

const edm::EDGetTokenT<reco::VertexCollection> DiMuonVertexMonitor::vertexToken_

private

Definition at line 51 of file DiMuonVertexMonitor.h.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DiMuonVertexMonitor()

◆ ~DiMuonVertexMonitor()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ fillDescriptions()

◆ findClosestVertex()

Member Data Documentation

◆ hCosPhi3D_

◆ hCosPhi_

◆ hCosPhiInv3D_

◆ hCosPhiInv_

◆ hCutFlow_

◆ hdxy_

◆ hdxyErr_

◆ hdz_

◆ hdzErr_

◆ hInvMass_

◆ hIP2d_

◆ hIP2dsig_

◆ hIP3d_

◆ hIP3dsig_

◆ hSVChi2_

◆ hSVCompatibility3D_

◆ hSVCompatibility_

◆ hSVDist3D_

◆ hSVDist3DErr_

◆ hSVDist3DSig_

◆ hSVDist_

◆ hSVDistErr_

◆ hSVDistSig_

◆ hSVNormChi2_

◆ hSVProb_

◆ massLimits_

◆ maxSVdist_

◆ MEFolderName_

◆ motherName_

◆ tracksToken_

◆ ttbESToken_

◆ useClosestVertex_

◆ vertexToken_