#include <PrimaryVertexMonitor.h>

Inheritance diagram for PrimaryVertexMonitor:

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

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

	PrimaryVertexMonitor (const edm::ParameterSet &pSet)

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

Private Types
using	VertexScore = edm::ValueMap< float >

Private Member Functions
void	pvTracksPlots (const reco::Vertex &v)

void	vertexPlots (const reco::Vertex &v, const reco::BeamSpot &beamSpot, int i)

Private Attributes
std::string	AlignmentLabel_

edm::InputTag	beamSpotInputTag_

edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_

MonitorElement *	bsBeamWidthX

MonitorElement *	bsBeamWidthY

MonitorElement *	bsDxdz

MonitorElement *	bsDydz

MonitorElement *	bsSigmaZ

MonitorElement *	bsType

MonitorElement *	bsX

MonitorElement *	bsY

MonitorElement *	bsZ

MonitorElement *	chi2ndf

MonitorElement *	chi2prob

edm::ParameterSet	conf_

std::string	dqmLabel

MonitorElement *	dxy

MonitorElement *	dxy2

MonitorElement *	dxyErr

MonitorElement *	dxyVsEta_pt1

MonitorElement *	dxyVsEta_pt10

MonitorElement *	dxyVsEtaVsPhi_pt1

MonitorElement *	dxyVsEtaVsPhi_pt10

MonitorElement *	dxyVsPhi_pt1

MonitorElement *	dxyVsPhi_pt10

MonitorElement *	dz

MonitorElement *	dzErr

MonitorElement *	dzVsEta_pt1

MonitorElement *	dzVsEta_pt10

MonitorElement *	dzVsEtaVsPhi_pt1

MonitorElement *	dzVsEtaVsPhi_pt10

MonitorElement *	dzVsPhi_pt1

MonitorElement *	dzVsPhi_pt10

bool	errorPrinted_

MonitorElement *	eta_pt1

MonitorElement *	eta_pt10

MonitorElement *	nans [2]

MonitorElement *	nbgvtx

MonitorElement *	nbtksinvtx [2]

MonitorElement *	nbvtx

int	ndof_

MonitorElement *	ntracks

MonitorElement *	ntracksVsZ [2]

MonitorElement *	phi_pt1

MonitorElement *	phi_pt10

MonitorElement *	score [2]

edm::EDGetTokenT< VertexScore >	scoreToken_

MonitorElement *	sumpt

std::string	TopFolderName_

MonitorElement *	trksWeight [2]

MonitorElement *	tt [2]

MonitorElement *	type [2]

bool	useHPfoAlignmentPlots_

edm::InputTag	vertexInputTag_

edm::EDGetTokenT< reco::VertexCollection >	vertexToken_

MonitorElement *	vtxchi2 [2]

MonitorElement *	vtxndf [2]

MonitorElement *	vtxprob [2]

MonitorElement *	weight

MonitorElement *	xDiff [2]

MonitorElement *	xerr [2]

MonitorElement *	xerrVsTrks [2]

MonitorElement *	xrec [2]

MonitorElement *	yDiff [2]

MonitorElement *	yerr [2]

MonitorElement *	yerrVsTrks [2]

MonitorElement *	yrec [2]

MonitorElement *	zerr [2]

MonitorElement *	zerrVsTrks [2]

MonitorElement *	zrec [2]

Static Private Attributes
static constexpr int	cmToUm = 10000

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

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

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

Definition at line 26 of file PrimaryVertexMonitor.h.

Member Typedef Documentation

◆ VertexScore

using PrimaryVertexMonitor::VertexScore = edm::ValueMap<float>

private

Definition at line 41 of file PrimaryVertexMonitor.h.

Constructor & Destructor Documentation

◆ PrimaryVertexMonitor()

PrimaryVertexMonitor::PrimaryVertexMonitor ( const edm::ParameterSet & pSet )

explicit

Definition at line 15 of file PrimaryVertexMonitor.cc.

References beamSpotInputTag_, beamspotToken_, edm::ParameterSet::getParameter(), scoreToken_, vertexInputTag_, and vertexToken_.

     : conf_(pSet),
       TopFolderName_(pSet.getParameter<std::string>("TopFolderName")),
       AlignmentLabel_(pSet.getParameter<std::string>("AlignmentLabel")),
       ndof_(pSet.getParameter<int>("ndof")),
       useHPfoAlignmentPlots_(pSet.getParameter<bool>("useHPforAlignmentPlots")),
       errorPrinted_(false),
       nbvtx(nullptr),
       bsX(nullptr),
       bsY(nullptr),
       bsZ(nullptr),
       bsSigmaZ(nullptr),
       bsDxdz(nullptr),
       bsDydz(nullptr),
       bsBeamWidthX(nullptr),
       bsBeamWidthY(nullptr),
       bsType(nullptr),
       sumpt(nullptr),
       ntracks(nullptr),
       weight(nullptr),
       chi2ndf(nullptr),
       chi2prob(nullptr),
       dxy(nullptr),
       dxy2(nullptr),
       dz(nullptr),
       dxyErr(nullptr),
       dzErr(nullptr),
       phi_pt1(nullptr),
       eta_pt1(nullptr),
       phi_pt10(nullptr),
       eta_pt10(nullptr),
       dxyVsPhi_pt1(nullptr),
       dzVsPhi_pt1(nullptr),
       dxyVsEta_pt1(nullptr),
       dzVsEta_pt1(nullptr),
       dxyVsEtaVsPhi_pt1(nullptr),
       dzVsEtaVsPhi_pt1(nullptr),
       dxyVsPhi_pt10(nullptr),
       dzVsPhi_pt10(nullptr),
       dxyVsEta_pt10(nullptr),
       dzVsEta_pt10(nullptr),
       dxyVsEtaVsPhi_pt10(nullptr),
       dzVsEtaVsPhi_pt10(nullptr) {
   vertexInputTag_ = pSet.getParameter<InputTag>("vertexLabel");
   beamSpotInputTag_ = pSet.getParameter<InputTag>("beamSpotLabel");
   vertexToken_ = consumes<reco::VertexCollection>(vertexInputTag_);
   scoreToken_ = consumes<VertexScore>(vertexInputTag_);
   beamspotToken_ = consumes<reco::BeamSpot>(beamSpotInputTag_);
 }

◆ ~PrimaryVertexMonitor()

PrimaryVertexMonitor::~PrimaryVertexMonitor ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 381 of file PrimaryVertexMonitor.cc.

References pwdgSkimBPark_cfi::beamSpot, beamSpotInputTag_, beamspotToken_, bsBeamWidthX, bsBeamWidthY, bsDxdz, bsDydz, bsSigmaZ, bsType, bsX, bsY, bsZ, cmToUm, errorPrinted_, dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), mps_fire::i, iEvent, edm::HandleBase::isValid(), nbgvtx, nbvtx, ndof_, convertSQLiteXML::ok, pvTracksPlots(), score, scoreToken_, mathSSE::sqrt(), findQualityFiles::v, vertexInputTag_, vertexPlots(), and vertexToken_.

                                                                                       {
   Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(vertexToken_, recVtxs);
 
   Handle<VertexScore> scores;
   iEvent.getByToken(scoreToken_, scores);
 
   edm::Handle<reco::BeamSpot> beamSpotHandle;
   iEvent.getByToken(beamspotToken_, beamSpotHandle);
 
   //
   // check for absent products and simply "return" in that case
   //
   if (recVtxs.isValid() == false || beamSpotHandle.isValid() == false) {
     edm::LogWarning("PrimaryVertexMonitor")
         << " Some products not available in the event: VertexCollection " << vertexInputTag_ << " " << recVtxs.isValid()
         << " BeamSpot " << beamSpotInputTag_ << " " << beamSpotHandle.isValid() << ". Skipping plots for this event";
     return;
   }
 
   // check upfront that refs to track are (likely) to be valid
   {
     bool ok = true;
     for (const auto& v : *recVtxs) {
       if (v.tracksSize() > 0) {
         const auto& ref = v.trackRefAt(0);
         if (ref.isNull() || !ref.isAvailable()) {
           if (!errorPrinted_)
             edm::LogWarning("PrimaryVertexMonitor")
                 << "Skipping vertex collection: " << vertexInputTag_
                 << " since likely the track collection the vertex has refs pointing to is missing (at least the first "
                    "TrackBaseRef is null or not available)";
           else
             errorPrinted_ = true;
           ok = false;
         }
       }
     }
     if (!ok)
       return;
   }
 
   BeamSpot beamSpot = *beamSpotHandle;
 
   nbvtx->Fill(recVtxs->size() * 1.);
   int ng = 0;
   for (auto const& vx : (*recVtxs))
     if (vx.isValid() && !vx.isFake() && vx.ndof() >= ndof_)
       ++ng;
   nbgvtx->Fill(ng * 1.);
 
   if (scores.isValid() && !(*scores).empty()) {
     auto pvScore = (*scores).get(0);
     score[1]->Fill(std::sqrt(pvScore));
     for (unsigned int i = 1; i < (*scores).size(); ++i)
       score[0]->Fill(std::sqrt((*scores).get(i)));
   }
 
   // fill PV tracks MEs (as now, for alignment)
   if (!recVtxs->empty()) {
     vertexPlots(recVtxs->front(), beamSpot, 1);
     pvTracksPlots(recVtxs->front());
 
     for (reco::VertexCollection::const_iterator v = recVtxs->begin() + 1; v != recVtxs->end(); ++v)
       vertexPlots(*v, beamSpot, 0);
   }
 
   // Beamline plots:
   bsX->Fill(beamSpot.x0());
   bsY->Fill(beamSpot.y0());
   bsZ->Fill(beamSpot.z0());
   bsSigmaZ->Fill(beamSpot.sigmaZ());
   bsDxdz->Fill(beamSpot.dxdz());
   bsDydz->Fill(beamSpot.dydz());
   bsBeamWidthX->Fill(beamSpot.BeamWidthX() * cmToUm);
   bsBeamWidthY->Fill(beamSpot.BeamWidthY() * cmToUm);
   bsType->Fill(beamSpot.type());
 }

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 67 of file PrimaryVertexMonitor.cc.

References AlignmentLabel_, beamSpotInputTag_, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), dqm::implementation::IBooker::bookProfile2D(), bsBeamWidthX, bsBeamWidthY, bsDxdz, bsDydz, bsSigmaZ, bsType, bsX, bsY, bsZ, chi2ndf, chi2prob, conf_, dqmLabel, dxy, dxy2, beam_dqm_sourceclient-live_cfg::DxyBin, dxyErr, beam_dqm_sourceclient-live_cfg::DxyMax, beam_dqm_sourceclient-live_cfg::DxyMin, dxyVsEta_pt1, dxyVsEta_pt10, dxyVsEtaVsPhi_pt1, dxyVsEtaVsPhi_pt10, dxyVsPhi_pt1, dxyVsPhi_pt10, dz, PrimaryVertexMonitor_cff::DzBin, dzErr, PrimaryVertexMonitor_cff::DzMax, PrimaryVertexMonitor_cff::DzMin, dzVsEta_pt1, dzVsEta_pt10, dzVsEtaVsPhi_pt1, dzVsEtaVsPhi_pt10, dzVsPhi_pt1, dzVsPhi_pt10, eta_pt1, eta_pt10, MonitorTrackInnerTrackMuons_cff::EtaBin, PrimaryVertexMonitor_cff::EtaBin2D, MonitorTrackInnerTrackMuons_cff::EtaMax, MonitorTrackInnerTrackMuons_cff::EtaMin, edm::ParameterSet::getParameter(), mps_fire::i, edm::InputTag::label(), nans, nbgvtx, nbtksinvtx, nbvtx, ntracks, ntracksVsZ, phi_pt1, phi_pt10, MonitorTrackSTAMuons_cfi::PhiBin, PrimaryVertexMonitor_cff::PhiBin2D, MonitorTrackSTAMuons_cfi::PhiMax, MonitorTrackSTAMuons_cfi::PhiMin, score, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, sumpt, TopFolderName_, trksWeight, vertexInputTag_, vtxchi2, vtxndf, vtxprob, xDiff, xerr, xerrVsTrks, xrec, yDiff, yerr, yerrVsTrks, yrec, zerr, zerrVsTrks, and zrec.

                                                                                                        {
   std::string dqmLabel = "";
 
   //
   // Book all histograms.
   //
 
   //  get the store
   dqmLabel = TopFolderName_ + "/" + vertexInputTag_.label();
   iBooker.setCurrentFolder(dqmLabel);
 
   //   xPos = iBooker.book1D ("xPos","x Coordinate" ,100, -0.1, 0.1);
 
   nbvtx = iBooker.book1D("vtxNbr", "Reconstructed Vertices in Event", 80, -0.5, 79.5);
   nbgvtx = iBooker.book1D("goodvtxNbr", "Reconstructed Good Vertices in Event", 80, -0.5, 79.5);
 
   // to be configured each year...
   auto vposx = conf_.getParameter<double>("Xpos");
   auto vposy = conf_.getParameter<double>("Ypos");
 
   nbtksinvtx[0] = iBooker.book1D("otherVtxTrksNbr", "Reconstructed Tracks in Vertex (other Vtx)", 40, -0.5, 99.5);
   ntracksVsZ[0] = iBooker.bookProfile(
       "otherVtxTrksVsZ", "Reconstructed Tracks in Vertex (other Vtx) vs Z", 80, -20., 20., 50, 0, 100, "");
   ntracksVsZ[0]->setAxisTitle("z-bs", 1);
   ntracksVsZ[0]->setAxisTitle("#tracks", 2);
 
   score[0] = iBooker.book1D("otherVtxScore", "sqrt(score) (other Vtx)", 100, 0., 400.);
   trksWeight[0] = iBooker.book1D("otherVtxTrksWeight", "Total weight of Tracks in Vertex (other Vtx)", 40, 0, 100.);
   vtxchi2[0] = iBooker.book1D("otherVtxChi2", "#chi^{2} (other Vtx)", 100, 0., 200.);
   vtxndf[0] = iBooker.book1D("otherVtxNdf", "ndof (other Vtx)", 100, 0., 200.);
   vtxprob[0] = iBooker.book1D("otherVtxProb", "#chi^{2} probability (other Vtx)", 100, 0., 1.);
   nans[0] = iBooker.book1D("otherVtxNans", "Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (other Vtx)", 9, 0.5, 9.5);
 
   nbtksinvtx[1] = iBooker.book1D("tagVtxTrksNbr", "Reconstructed Tracks in Vertex (tagged Vtx)", 100, -0.5, 99.5);
   ntracksVsZ[1] = iBooker.bookProfile(
       "tagVtxTrksVsZ", "Reconstructed Tracks in Vertex (tagged Vtx) vs Z", 80, -20., 20., 50, 0, 100, "");
   ntracksVsZ[1]->setAxisTitle("z-bs", 1);
   ntracksVsZ[1]->setAxisTitle("#tracks", 2);
 
   score[1] = iBooker.book1D("tagVtxScore", "sqrt(score) (tagged Vtx)", 100, 0., 400.);
   trksWeight[1] = iBooker.book1D("tagVtxTrksWeight", "Total weight of Tracks in Vertex (tagged Vtx)", 100, 0, 100.);
   vtxchi2[1] = iBooker.book1D("tagVtxChi2", "#chi^{2} (tagged Vtx)", 100, 0., 200.);
   vtxndf[1] = iBooker.book1D("tagVtxNdf", "ndof (tagged Vtx)", 100, 0., 200.);
   vtxprob[1] = iBooker.book1D("tagVtxProb", "#chi^{2} probability (tagged Vtx)", 100, 0., 1.);
   nans[1] = iBooker.book1D("tagVtxNans", "Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (tagged Vtx)", 9, 0.5, 9.5);
 
   xrec[0] = iBooker.book1D("otherPosX", "Position x Coordinate (other Vtx)", 100, vposx - 0.1, vposx + 0.1);
   yrec[0] = iBooker.book1D("otherPosY", "Position y Coordinate (other Vtx)", 100, vposy - 0.1, vposy + 0.1);
   zrec[0] = iBooker.book1D("otherPosZ", "Position z Coordinate (other Vtx)", 100, -20., 20.);
   xDiff[0] = iBooker.book1D("otherDiffX", "X distance from BeamSpot (other Vtx)", 100, -500, 500);
   yDiff[0] = iBooker.book1D("otherDiffY", "Y distance from BeamSpot (other Vtx)", 100, -500, 500);
   xerr[0] = iBooker.book1D("otherErrX", "Uncertainty x Coordinate (other Vtx)", 100, 0., 100);
   yerr[0] = iBooker.book1D("otherErrY", "Uncertainty y Coordinate (other Vtx)", 100, 0., 100);
   zerr[0] = iBooker.book1D("otherErrZ", "Uncertainty z Coordinate (other Vtx)", 100, 0., 100);
   xerrVsTrks[0] = iBooker.book2D(
       "otherErrVsWeightX", "Uncertainty x Coordinate vs. track weight (other Vtx)", 100, 0, 100., 100, 0., 100);
   yerrVsTrks[0] = iBooker.book2D(
       "otherErrVsWeightY", "Uncertainty y Coordinate vs. track weight (other Vtx)", 100, 0, 100., 100, 0., 100);
   zerrVsTrks[0] = iBooker.book2D(
       "otherErrVsWeightZ", "Uncertainty z Coordinate vs. track weight (other Vtx)", 100, 0, 100., 100, 0., 100);
 
   xrec[1] = iBooker.book1D("tagPosX", "Position x Coordinate (tagged Vtx)", 100, vposx - 0.1, vposx + 0.1);
   yrec[1] = iBooker.book1D("tagPosY", "Position y Coordinate (tagged Vtx)", 100, vposy - 0.1, vposy + 0.1);
   zrec[1] = iBooker.book1D("tagPosZ", "Position z Coordinate (tagged Vtx)", 100, -20., 20.);
   xDiff[1] = iBooker.book1D("tagDiffX", "X distance from BeamSpot (tagged Vtx)", 100, -500, 500);
   yDiff[1] = iBooker.book1D("tagDiffY", "Y distance from BeamSpot (tagged Vtx)", 100, -500, 500);
   xerr[1] = iBooker.book1D("tagErrX", "Uncertainty x Coordinate (tagged Vtx)", 100, 0., 100);
   yerr[1] = iBooker.book1D("tagErrY", "Uncertainty y Coordinate (tagged Vtx)", 100, 0., 100);
   zerr[1] = iBooker.book1D("tagErrZ", "Uncertainty z Coordinate (tagged Vtx)", 100, 0., 100);
   xerrVsTrks[1] = iBooker.book2D(
       "tagErrVsWeightX", "Uncertainty x Coordinate vs. track weight (tagged Vtx)", 100, 0, 100., 100, 0., 100);
   yerrVsTrks[1] = iBooker.book2D(
       "tagErrVsWeightY", "Uncertainty y Coordinate vs. track weight (tagged Vtx)", 100, 0, 100., 100, 0., 100);
   zerrVsTrks[1] = iBooker.book2D(
       "tagErrVsWeightZ", "Uncertainty z Coordinate vs. track weight (tagged Vtx)", 100, 0, 100., 100, 0., 100);
 
   type[0] = iBooker.book1D("otherType", "Vertex type (other Vtx)", 3, -0.5, 2.5);
   type[1] = iBooker.book1D("tagType", "Vertex type (tagged Vtx)", 3, -0.5, 2.5);
   for (int i = 0; i < 2; ++i) {
     type[i]->setBinLabel(1, "Valid, real");
     type[i]->setBinLabel(2, "Valid, fake");
     type[i]->setBinLabel(3, "Invalid");
   }
 
   //  get the store
   dqmLabel = TopFolderName_ + "/" + beamSpotInputTag_.label();
   iBooker.setCurrentFolder(dqmLabel);
 
   bsX = iBooker.book1D("bsX", "BeamSpot x0", 100, vposx - 0.1, vposx + 0.1);
   bsY = iBooker.book1D("bsY", "BeamSpot y0", 100, vposy - 0.1, vposy + 0.1);
   bsZ = iBooker.book1D("bsZ", "BeamSpot z0", 100, -2., 2.);
   bsSigmaZ = iBooker.book1D("bsSigmaZ", "BeamSpot sigmaZ", 100, 0., 10.);
   bsDxdz = iBooker.book1D("bsDxdz", "BeamSpot dxdz", 100, -0.0003, 0.0003);
   bsDydz = iBooker.book1D("bsDydz", "BeamSpot dydz", 100, -0.0003, 0.0003);
   bsBeamWidthX = iBooker.book1D("bsBeamWidthX", "BeamSpot BeamWidthX", 100, 0., 100.);
   bsBeamWidthY = iBooker.book1D("bsBeamWidthY", "BeamSpot BeamWidthY", 100, 0., 100.);
   bsType = iBooker.book1D("bsType", "BeamSpot type", 4, -1.5, 2.5);
   bsType->setBinLabel(1, "Unknown");
   bsType->setBinLabel(2, "Fake");
   bsType->setBinLabel(3, "LHC");
   bsType->setBinLabel(4, "Tracker");
 
   //  get the store
   dqmLabel = TopFolderName_ + "/" + AlignmentLabel_;
   iBooker.setCurrentFolder(dqmLabel);
 
   int TKNoBin = conf_.getParameter<int>("TkSizeBin");
   double TKNoMin = conf_.getParameter<double>("TkSizeMin");
   double TKNoMax = conf_.getParameter<double>("TkSizeMax");
 
   int DxyBin = conf_.getParameter<int>("DxyBin");
   double DxyMin = conf_.getParameter<double>("DxyMin");
   double DxyMax = conf_.getParameter<double>("DxyMax");
 
   int DzBin = conf_.getParameter<int>("DzBin");
   double DzMin = conf_.getParameter<double>("DzMin");
   double DzMax = conf_.getParameter<double>("DzMax");
 
   int PhiBin = conf_.getParameter<int>("PhiBin");
   int PhiBin2D = conf_.getParameter<int>("PhiBin2D");
   double PhiMin = conf_.getParameter<double>("PhiMin");
   double PhiMax = conf_.getParameter<double>("PhiMax");
 
   int EtaBin = conf_.getParameter<int>("EtaBin");
   int EtaBin2D = conf_.getParameter<int>("EtaBin2D");
   double EtaMin = conf_.getParameter<double>("EtaMin");
   double EtaMax = conf_.getParameter<double>("EtaMax");
 
   ntracks = iBooker.book1D("ntracks", "number of PV tracks (p_{T} > 1 GeV)", TKNoBin, TKNoMin, TKNoMax);
   ntracks->setAxisTitle("Number of PV Tracks (p_{T} > 1 GeV) per Event", 1);
   ntracks->setAxisTitle("Number of Event", 2);
 
   weight = iBooker.book1D("weight", "weight of PV tracks (p_{T} > 1 GeV)", 100, 0., 1.);
   weight->setAxisTitle("weight of PV Tracks (p_{T} > 1 GeV) per Event", 1);
   weight->setAxisTitle("Number of Event", 2);
 
   sumpt = iBooker.book1D("sumpt", "#Sum p_{T} of PV tracks (p_{T} > 1 GeV)", 100, -0.5, 249.5);
   chi2ndf = iBooker.book1D("chi2ndf", "PV tracks (p_{T} > 1 GeV) #chi^{2}/ndof", 100, 0., 20.);
   chi2prob = iBooker.book1D("chi2prob", "PV tracks (p_{T} > 1 GeV) #chi^{2} probability", 100, 0., 1.);
 
   dxy = iBooker.book1D("dxy", "PV tracks (p_{T} > 1 GeV) d_{xy} (#mum)", DxyBin, DxyMin, DxyMax);
   dxy2 = iBooker.book1D("dxyzoom", "PV tracks (p_{T} > 1 GeV) d_{xy} (#mum)", DxyBin, DxyMin / 5., DxyMax / 5.);
   dxyErr = iBooker.book1D("dxyErr", "PV tracks (p_{T} > 1 GeV) d_{xy} error (#mum)", 100, 0., 2000.);
   dz = iBooker.book1D("dz", "PV tracks (p_{T} > 1 GeV) d_{z} (#mum)", DzBin, DzMin, DzMax);
   dzErr = iBooker.book1D("dzErr", "PV tracks (p_{T} > 1 GeV) d_{z} error(#mum)", 100, 0., 10000.);
 
   phi_pt1 = iBooker.book1D("phi_pt1", "PV tracks (p_{T} > 1 GeV) #phi; PV tracks #phi;#tracks", PhiBin, PhiMin, PhiMax);
   eta_pt1 = iBooker.book1D("eta_pt1", "PV tracks (p_{T} > 1 GeV) #eta; PV tracks #eta;#tracks", EtaBin, EtaMin, EtaMax);
   phi_pt10 =
       iBooker.book1D("phi_pt10", "PV tracks (p_{T} > 10 GeV) #phi; PV tracks #phi;#tracks", PhiBin, PhiMin, PhiMax);
   eta_pt10 =
       iBooker.book1D("eta_pt10", "PV tracks (p_{T} > 10 GeV) #phi; PV tracks #eta;#tracks", EtaBin, EtaMin, EtaMax);
 
   dxyVsPhi_pt1 = iBooker.bookProfile("dxyVsPhi_pt1",
                                      "PV tracks (p_{T} > 1 GeV) d_{xy} (#mum) VS track #phi",
                                      PhiBin,
                                      PhiMin,
                                      PhiMax,
                                      DxyBin,
                                      DxyMin,
                                      DxyMax,
                                      "");
   dxyVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 1);
   dxyVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy} (#mum)", 2);
 
   dzVsPhi_pt1 = iBooker.bookProfile("dzVsPhi_pt1",
                                     "PV tracks (p_{T} > 1 GeV) d_{z} (#mum) VS track #phi",
                                     PhiBin,
                                     PhiMin,
                                     PhiMax,
                                     DzBin,
                                     DzMin,
                                     DzMax,
                                     "");
   dzVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 1);
   dzVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z} (#mum)", 2);
 
   dxyVsEta_pt1 = iBooker.bookProfile("dxyVsEta_pt1",
                                      "PV tracks (p_{T} > 1 GeV) d_{xy} (#mum) VS track #eta",
                                      EtaBin,
                                      EtaMin,
                                      EtaMax,
                                      DxyBin,
                                      DxyMin,
                                      DxyMax,
                                      "");
   dxyVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
   dxyVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy} (#mum)", 2);
 
   dzVsEta_pt1 = iBooker.bookProfile("dzVsEta_pt1",
                                     "PV tracks (p_{T} > 1 GeV) d_{z} (#mum) VS track #eta",
                                     EtaBin,
                                     EtaMin,
                                     EtaMax,
                                     DzBin,
                                     DzMin,
                                     DzMax,
                                     "");
   dzVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
   dzVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z} (#mum)", 2);
 
   dxyVsEtaVsPhi_pt1 = iBooker.bookProfile2D("dxyVsEtaVsPhi_pt1",
                                             "PV tracks (p_{T} > 1 GeV) d_{xy} (#mum) VS track #eta VS track #phi",
                                             EtaBin2D,
                                             EtaMin,
                                             EtaMax,
                                             PhiBin2D,
                                             PhiMin,
                                             PhiMax,
                                             DxyBin,
                                             DxyMin,
                                             DxyMax,
                                             "");
   dxyVsEtaVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
   dxyVsEtaVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 2);
   dxyVsEtaVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy} (#mum)", 3);
 
   dzVsEtaVsPhi_pt1 = iBooker.bookProfile2D("dzVsEtaVsPhi_pt1",
                                            "PV tracks (p_{T} > 1 GeV) d_{z} (#mum) VS track #eta VS track #phi",
                                            EtaBin2D,
                                            EtaMin,
                                            EtaMax,
                                            PhiBin2D,
                                            PhiMin,
                                            PhiMax,
                                            DzBin,
                                            DzMin,
                                            DzMax,
                                            "");
   dzVsEtaVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
   dzVsEtaVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 2);
   dzVsEtaVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z} (#mum)", 3);
 
   dxyVsPhi_pt10 = iBooker.bookProfile("dxyVsPhi_pt10",
                                       "PV tracks (p_{T} > 10 GeV) d_{xy} (#mum) VS track #phi",
                                       PhiBin,
                                       PhiMin,
                                       PhiMax,
                                       DxyBin,
                                       DxyMin,
                                       DxyMax,
                                       "");
   dxyVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 1);
   dxyVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy} (#mum)", 2);
 
   dzVsPhi_pt10 = iBooker.bookProfile("dzVsPhi_pt10",
                                      "PV tracks (p_{T} > 10 GeV) d_{z} (#mum) VS track #phi",
                                      PhiBin,
                                      PhiMin,
                                      PhiMax,
                                      DzBin,
                                      DzMin,
                                      DzMax,
                                      "");
   dzVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 1);
   dzVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z} (#mum)", 2);
 
   dxyVsEta_pt10 = iBooker.bookProfile("dxyVsEta_pt10",
                                       "PV tracks (p_{T} > 10 GeV) d_{xy} (#mum) VS track #eta",
                                       EtaBin,
                                       EtaMin,
                                       EtaMax,
                                       DxyBin,
                                       DxyMin,
                                       DxyMax,
                                       "");
   dxyVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
   dxyVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy} (#mum)", 2);
 
   dzVsEta_pt10 = iBooker.bookProfile("dzVsEta_pt10",
                                      "PV tracks (p_{T} > 10 GeV) d_{z} (#mum) VS track #eta",
                                      EtaBin,
                                      EtaMin,
                                      EtaMax,
                                      DzBin,
                                      DzMin,
                                      DzMax,
                                      "");
   dzVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
   dzVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z} (#mum)", 2);
 
   dxyVsEtaVsPhi_pt10 = iBooker.bookProfile2D("dxyVsEtaVsPhi_pt10",
                                              "PV tracks (p_{T} > 10 GeV) d_{xy} (#mum) VS track #eta VS track #phi",
                                              EtaBin2D,
                                              EtaMin,
                                              EtaMax,
                                              PhiBin2D,
                                              PhiMin,
                                              PhiMax,
                                              DxyBin,
                                              DxyMin,
                                              DxyMax,
                                              "");
   dxyVsEtaVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
   dxyVsEtaVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 2);
   dxyVsEtaVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy} (#mum)", 3);
 
   dzVsEtaVsPhi_pt10 = iBooker.bookProfile2D("dzVsEtaVsPhi_pt10",
                                             "PV tracks (p_{T} > 10 GeV) d_{z} (#mum) VS track #eta VS track #phi",
                                             EtaBin2D,
                                             EtaMin,
                                             EtaMax,
                                             PhiBin2D,
                                             PhiMin,
                                             PhiMax,
                                             DzBin,
                                             DzMin,
                                             DzMax,
                                             "");
   dzVsEtaVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
   dzVsEtaVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 2);
   dzVsEtaVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z} (#mum)", 3);
 }

◆ pvTracksPlots()

void PrimaryVertexMonitor::pvTracksPlots ( const reco::Vertex & v )

private

Definition at line 460 of file PrimaryVertexMonitor.cc.

References chi2ndf, chi2prob, cmToUm, dxy, dxy2, dxyErr, dxyVsEta_pt1, dxyVsEta_pt10, dxyVsEtaVsPhi_pt1, dxyVsEtaVsPhi_pt10, dxyVsPhi_pt1, dxyVsPhi_pt10, dz, dzErr, dzVsEta_pt1, dzVsEta_pt10, dzVsEtaVsPhi_pt1, dzVsEtaVsPhi_pt10, dzVsPhi_pt1, dzVsPhi_pt10, PVValHelper::eta, eta_pt1, eta_pt10, dqm::impl::MonitorElement::Fill(), reco::TrackBase::highPurity, createfilelist::int, BeamSpotPI::nTracks, ntracks, phi, phi_pt1, phi_pt10, DiDispStaMuonMonitor_cfi::pt, sumpt, submitPVValidationJobs::t, useHPfoAlignmentPlots_, findQualityFiles::v, and w().

Referenced by analyze().

                                                         {
   if (!v.isValid())
     return;
   if (v.isFake())
     return;
 
   if (v.tracksSize() == 0) {
     ntracks->Fill(0);
     return;
   }
 
   const math::XYZPoint myVertex(v.position().x(), v.position().y(), v.position().z());
 
   size_t nTracks = 0;
   float sumPT = 0.;
 
   for (reco::Vertex::trackRef_iterator t = v.tracks_begin(); t != v.tracks_end(); t++) {
     bool isHighPurity = (**t).quality(reco::TrackBase::highPurity);
     if (!isHighPurity && useHPfoAlignmentPlots_)
       continue;
 
     float pt = (**t).pt();
     if (pt < 1.)
       continue;
 
     nTracks++;
 
     float eta = (**t).eta();
     float phi = (**t).phi();
 
     float w = v.trackWeight(*t);
     float chi2NDF = (**t).normalizedChi2();
     float chi2Prob = TMath::Prob((**t).chi2(), (int)(**t).ndof());
     float Dxy = (**t).dxy(myVertex) * cmToUm;
     float Dz = (**t).dz(myVertex) * cmToUm;
     float DxyErr = (**t).dxyError() * cmToUm;
     float DzErr = (**t).dzError() * cmToUm;
 
     sumPT += pt * pt;
 
     // fill MEs
     weight->Fill(w);
     chi2ndf->Fill(chi2NDF);
     chi2prob->Fill(chi2Prob);
     dxy->Fill(Dxy);
     dxy2->Fill(Dxy);
     dz->Fill(Dz);
     dxyErr->Fill(DxyErr);
     dzErr->Fill(DzErr);
     phi_pt1->Fill(phi);
     eta_pt1->Fill(eta);
 
     dxyVsPhi_pt1->Fill(phi, Dxy);
     dzVsPhi_pt1->Fill(phi, Dz);
     dxyVsEta_pt1->Fill(eta, Dxy);
     dzVsEta_pt1->Fill(eta, Dz);
     dxyVsEtaVsPhi_pt1->Fill(eta, phi, Dxy);
     dzVsEtaVsPhi_pt1->Fill(eta, phi, Dz);
 
     if (pt < 10.)
       continue;
 
     phi_pt10->Fill(phi);
     eta_pt10->Fill(eta);
     dxyVsPhi_pt10->Fill(phi, Dxy);
     dzVsPhi_pt10->Fill(phi, Dz);
     dxyVsEta_pt10->Fill(eta, Dxy);
     dzVsEta_pt10->Fill(eta, Dz);
     dxyVsEtaVsPhi_pt10->Fill(eta, phi, Dxy);
     dzVsEtaVsPhi_pt10->Fill(eta, phi, Dz);
   }
   ntracks->Fill(float(nTracks));
   sumpt->Fill(sumPT);
 }

◆ vertexPlots()

void PrimaryVertexMonitor::vertexPlots	(	const reco::Vertex &	v,
		const reco::BeamSpot &	beamSpot,
		int	i
	)

private

Definition at line 535 of file PrimaryVertexMonitor.cc.

References pwdgSkimBPark_cfi::beamSpot, ChiSquaredProbability(), cmToUm, dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), mps_fire::i, edm::isNotFinite(), dqmiolumiharvest::j, dqmdumpme::k, nans, nbtksinvtx, ntracksVsZ, submitPVValidationJobs::t, trksWeight, findQualityFiles::v, vtxchi2, vtxndf, vtxprob, xDiff, xerr, xerrVsTrks, xrec, yDiff, yerr, yerrVsTrks, yrec, zerr, zerrVsTrks, and zrec.

Referenced by analyze().

                                                                                        {
   if (i < 0 || i > 1)
     return;
   if (!v.isValid())
     type[i]->Fill(2.);
   else if (v.isFake())
     type[i]->Fill(1.);
   else
     type[i]->Fill(0.);
 
   if (v.isValid() && !v.isFake()) {
     float weight = 0;
     for (reco::Vertex::trackRef_iterator t = v.tracks_begin(); t != v.tracks_end(); t++)
       weight += v.trackWeight(*t);
     trksWeight[i]->Fill(weight);
     nbtksinvtx[i]->Fill(v.tracksSize());
     ntracksVsZ[i]->Fill(v.position().z() - beamSpot.z0(), v.tracksSize());
 
     vtxchi2[i]->Fill(v.chi2());
     vtxndf[i]->Fill(v.ndof());
     vtxprob[i]->Fill(ChiSquaredProbability(v.chi2(), v.ndof()));
 
     xrec[i]->Fill(v.position().x());
     yrec[i]->Fill(v.position().y());
     zrec[i]->Fill(v.position().z());
 
     float xb = beamSpot.x0() + beamSpot.dxdz() * (v.position().z() - beamSpot.z0());
     float yb = beamSpot.y0() + beamSpot.dydz() * (v.position().z() - beamSpot.z0());
     xDiff[i]->Fill((v.position().x() - xb) * cmToUm);
     yDiff[i]->Fill((v.position().y() - yb) * cmToUm);
 
     xerr[i]->Fill(v.xError() * cmToUm);
     yerr[i]->Fill(v.yError() * cmToUm);
     zerr[i]->Fill(v.zError() * cmToUm);
     xerrVsTrks[i]->Fill(weight, v.xError() * cmToUm);
     yerrVsTrks[i]->Fill(weight, v.yError() * cmToUm);
     zerrVsTrks[i]->Fill(weight, v.zError() * cmToUm);
 
     nans[i]->Fill(1., edm::isNotFinite(v.position().x()) * 1.);
     nans[i]->Fill(2., edm::isNotFinite(v.position().y()) * 1.);
     nans[i]->Fill(3., edm::isNotFinite(v.position().z()) * 1.);
 
     int index = 3;
     for (int k = 0; k != 3; k++) {
       for (int j = k; j != 3; j++) {
         index++;
         nans[i]->Fill(index * 1., edm::isNotFinite(v.covariance(k, j)) * 1.);
         // in addition, diagonal element must be positive
         if (j == k && v.covariance(k, j) < 0) {
           nans[i]->Fill(index * 1., 1.);
         }
       }
     }
   }
 }