PrimaryVertexMonitor Class Reference

#include <PrimaryVertexMonitor.h>

Inheritance diagram for PrimaryVertexMonitor:

Classes
struct	IPMonitoring

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 *	dxy2
IPMonitoring	dxy_pt1
IPMonitoring	dxy_pt10
IPMonitoring	dz_pt1
IPMonitoring	dz_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 52 of file PrimaryVertexMonitor.h.

Constructor & Destructor Documentation

PrimaryVertexMonitor()

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

explicit

Definition at line 17 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),
      phi_pt1(nullptr),
      eta_pt1(nullptr),
      phi_pt10(nullptr),
      eta_pt10(nullptr),
      dxy2(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 327 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(), vertexToken_, and run3scouting_cff::vx.

                                                                                      {
  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 53 of file PrimaryVertexMonitor.cc.

References AlignmentLabel_, beamSpotInputTag_, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), PrimaryVertexMonitor::IPMonitoring::bookIPMonitor(), dqm::implementation::IBooker::bookProfile(), bsBeamWidthX, bsBeamWidthY, bsDxdz, bsDydz, bsSigmaZ, bsType, bsX, bsY, bsZ, chi2ndf, chi2prob, conf_, dqmLabel, dxy2, dxy_pt1, dxy_pt10, beam_dqm_sourceclient-live_cfg::DxyBin, beam_dqm_sourceclient-live_cfg::DxyMax, beam_dqm_sourceclient-live_cfg::DxyMin, dz_pt1, dz_pt10, eta_pt1, eta_pt10, MonitorTrackInnerTrackMuons_cff::EtaBin, 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, MonitorTrackSTAMuons_cfi::PhiMax, MonitorTrackSTAMuons_cfi::PhiMin, PrimaryVertexMonitor::IPMonitoring::pTcut_, score, dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, sumpt, TopFolderName_, trksWeight, PrimaryVertexMonitor::IPMonitoring::varname_, 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 PhiBin = conf_.getParameter<int>("PhiBin");
  double PhiMin = conf_.getParameter<double>("PhiMin");
  double PhiMax = conf_.getParameter<double>("PhiMax");

  int EtaBin = conf_.getParameter<int>("EtaBin");
  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.);

  dxy2 = iBooker.book1D("dxyzoom", "PV tracks (p_{T} > 1 GeV) d_{xy} (#mum)", DxyBin, DxyMin / 5., DxyMax / 5.);

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

  // initialize and book the monitors;
  dxy_pt1.varname_ = "xy";
  dxy_pt1.pTcut_ = 1.f;
  dxy_pt1.bookIPMonitor(iBooker, conf_);

  dxy_pt10.varname_ = "xy";
  dxy_pt10.pTcut_ = 10.f;
  dxy_pt10.bookIPMonitor(iBooker, conf_);

  dz_pt1.varname_ = "z";
  dz_pt1.pTcut_ = 1.f;
  dz_pt1.bookIPMonitor(iBooker, conf_);

  dz_pt10.varname_ = "z";
  dz_pt10.pTcut_ = 10.f;
  dz_pt10.bookIPMonitor(iBooker, conf_);
}

pvTracksPlots()

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

private

Definition at line 406 of file PrimaryVertexMonitor.cc.

References chi2ndf, chi2prob, cmToUm, dxy2, dxy_pt1, dxy_pt10, dz_pt1, dz_pt10, PVValHelper::eta, eta_pt1, eta_pt10, dqm::impl::MonitorElement::Fill(), reco::TrackBase::highPurity, createfilelist::int, PrimaryVertexMonitor::IPMonitoring::IP_, PrimaryVertexMonitor::IPMonitoring::IPErr_, PrimaryVertexMonitor::IPMonitoring::IPErrVsEta_, PrimaryVertexMonitor::IPMonitoring::IPErrVsEtaVsPhi_, PrimaryVertexMonitor::IPMonitoring::IPErrVsPhi_, PrimaryVertexMonitor::IPMonitoring::IPVsEta_, PrimaryVertexMonitor::IPMonitoring::IPVsEtaVsPhi_, PrimaryVertexMonitor::IPMonitoring::IPVsPhi_, 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
    phi_pt1->Fill(phi);
    eta_pt1->Fill(eta);

    weight->Fill(w);
    chi2ndf->Fill(chi2NDF);
    chi2prob->Fill(chi2Prob);
    dxy2->Fill(Dxy);

    dxy_pt1.IP_->Fill(Dxy);
    dxy_pt1.IPVsPhi_->Fill(phi, Dxy);
    dxy_pt1.IPVsEta_->Fill(eta, Dxy);
    dxy_pt1.IPVsEtaVsPhi_->Fill(eta, phi, Dxy);

    dxy_pt1.IPErr_->Fill(DxyErr);
    dxy_pt1.IPErrVsPhi_->Fill(phi, DxyErr);
    dxy_pt1.IPErrVsEta_->Fill(eta, DxyErr);
    dxy_pt1.IPErrVsEtaVsPhi_->Fill(eta, phi, DxyErr);

    dz_pt1.IP_->Fill(Dz);
    dz_pt1.IPVsPhi_->Fill(phi, Dz);
    dz_pt1.IPVsEta_->Fill(eta, Dz);
    dz_pt1.IPVsEtaVsPhi_->Fill(eta, phi, Dz);

    dz_pt1.IPErr_->Fill(DzErr);
    dz_pt1.IPErrVsPhi_->Fill(phi, DzErr);
    dz_pt1.IPErrVsEta_->Fill(eta, DzErr);
    dz_pt1.IPErrVsEtaVsPhi_->Fill(eta, phi, DzErr);

    if (pt < 10.)
      continue;

    phi_pt10->Fill(phi);
    eta_pt10->Fill(eta);

    dxy_pt10.IP_->Fill(Dxy);
    dxy_pt10.IPVsPhi_->Fill(phi, Dxy);
    dxy_pt10.IPVsEta_->Fill(eta, Dxy);
    dxy_pt10.IPVsEtaVsPhi_->Fill(eta, phi, Dxy);

    dxy_pt10.IPErr_->Fill(DxyErr);
    dxy_pt10.IPErrVsPhi_->Fill(phi, DxyErr);
    dxy_pt10.IPErrVsEta_->Fill(eta, DxyErr);
    dxy_pt10.IPErrVsEtaVsPhi_->Fill(eta, phi, DxyErr);

    dz_pt10.IP_->Fill(Dz);
    dz_pt10.IPVsPhi_->Fill(phi, Dz);
    dz_pt10.IPVsEta_->Fill(eta, Dz);
    dz_pt10.IPVsEtaVsPhi_->Fill(eta, phi, Dz);

    dz_pt10.IPErr_->Fill(DzErr);
    dz_pt10.IPErrVsPhi_->Fill(phi, DzErr);
    dz_pt10.IPErrVsEta_->Fill(eta, DzErr);
    dz_pt10.IPErrVsEtaVsPhi_->Fill(eta, phi, DzErr);
  }
  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 505 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.);
        }
      }
    }
  }
}

Member Data Documentation

AlignmentLabel_

std::string PrimaryVertexMonitor::AlignmentLabel_

private

Definition at line 62 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

beamSpotInputTag_

edm::InputTag PrimaryVertexMonitor::beamSpotInputTag_

private

Definition at line 55 of file PrimaryVertexMonitor.h.

Referenced by analyze(), bookHistograms(), and PrimaryVertexMonitor().

beamspotToken_

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

private

Definition at line 51 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and PrimaryVertexMonitor().

bsBeamWidthX

MonitorElement * PrimaryVertexMonitor::bsBeamWidthX

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsBeamWidthY

MonitorElement * PrimaryVertexMonitor::bsBeamWidthY

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsDxdz

MonitorElement * PrimaryVertexMonitor::bsDxdz

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsDydz

MonitorElement * PrimaryVertexMonitor::bsDydz

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsSigmaZ

MonitorElement * PrimaryVertexMonitor::bsSigmaZ

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsType

MonitorElement * PrimaryVertexMonitor::bsType

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsX

MonitorElement* PrimaryVertexMonitor::bsX

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsY

MonitorElement * PrimaryVertexMonitor::bsY

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

bsZ

MonitorElement * PrimaryVertexMonitor::bsZ

private

Definition at line 77 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

chi2ndf

MonitorElement * PrimaryVertexMonitor::chi2ndf

private

Definition at line 79 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

chi2prob

MonitorElement * PrimaryVertexMonitor::chi2prob

private

Definition at line 79 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

cmToUm

constexpr int PrimaryVertexMonitor::cmToUm = 10000

staticprivate

Definition at line 67 of file PrimaryVertexMonitor.h.

Referenced by analyze(), pvTracksPlots(), and vertexPlots().

conf_

edm::ParameterSet PrimaryVertexMonitor::conf_

private

Definition at line 57 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

dqmLabel

std::string PrimaryVertexMonitor::dqmLabel

private

Definition at line 59 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

dxy2

MonitorElement* PrimaryVertexMonitor::dxy2

private

Definition at line 83 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

dxy_pt1

IPMonitoring PrimaryVertexMonitor::dxy_pt1

private

Definition at line 86 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

dxy_pt10

IPMonitoring PrimaryVertexMonitor::dxy_pt10

private

Definition at line 87 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

dz_pt1

IPMonitoring PrimaryVertexMonitor::dz_pt1

private

Definition at line 89 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

dz_pt10

IPMonitoring PrimaryVertexMonitor::dz_pt10

private

Definition at line 90 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

errorPrinted_

bool PrimaryVertexMonitor::errorPrinted_

private

Definition at line 65 of file PrimaryVertexMonitor.h.

Referenced by analyze().

eta_pt1

MonitorElement * PrimaryVertexMonitor::eta_pt1

private

Definition at line 80 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

eta_pt10

MonitorElement * PrimaryVertexMonitor::eta_pt10

private

Definition at line 81 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

nans

MonitorElement * PrimaryVertexMonitor::nans[2]

private

Definition at line 75 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

nbgvtx

MonitorElement * PrimaryVertexMonitor::nbgvtx

private

Definition at line 70 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

nbtksinvtx

MonitorElement * PrimaryVertexMonitor::nbtksinvtx[2]

private

Definition at line 70 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

nbvtx

MonitorElement* PrimaryVertexMonitor::nbvtx

private

Definition at line 70 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

ndof_

int PrimaryVertexMonitor::ndof_

private

Definition at line 63 of file PrimaryVertexMonitor.h.

Referenced by analyze().

ntracks

MonitorElement * PrimaryVertexMonitor::ntracks

private

Definition at line 79 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

ntracksVsZ

MonitorElement* PrimaryVertexMonitor::ntracksVsZ[2]

private

Definition at line 74 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

phi_pt1

MonitorElement* PrimaryVertexMonitor::phi_pt1

private

Definition at line 80 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

phi_pt10

MonitorElement* PrimaryVertexMonitor::phi_pt10

private

Definition at line 81 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

score

MonitorElement * PrimaryVertexMonitor::score[2]

private

Definition at line 70 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

scoreToken_

edm::EDGetTokenT<VertexScore> PrimaryVertexMonitor::scoreToken_

private

Definition at line 53 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and PrimaryVertexMonitor().

sumpt

MonitorElement* PrimaryVertexMonitor::sumpt

private

Definition at line 79 of file PrimaryVertexMonitor.h.

Referenced by ZMuMuRochCorAnalyzer.DiMuon::__str__(), bookHistograms(), and pvTracksPlots().

TopFolderName_

std::string PrimaryVertexMonitor::TopFolderName_

private

Definition at line 61 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

trksWeight

MonitorElement * PrimaryVertexMonitor::trksWeight[2]

private

Definition at line 70 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

tt

MonitorElement* PrimaryVertexMonitor::tt[2]

private

Definition at line 71 of file PrimaryVertexMonitor.h.

type

MonitorElement* PrimaryVertexMonitor::type[2]

private

Definition at line 76 of file PrimaryVertexMonitor.h.

Referenced by core.autovars.NTupleVariable::makeBranch(), DQM.DQMReader::read_objects(), and core.AutoHandle.AutoHandle::ReallyLoad().

useHPfoAlignmentPlots_

bool PrimaryVertexMonitor::useHPfoAlignmentPlots_

private

Definition at line 64 of file PrimaryVertexMonitor.h.

Referenced by pvTracksPlots().

vertexInputTag_

edm::InputTag PrimaryVertexMonitor::vertexInputTag_

private

Definition at line 55 of file PrimaryVertexMonitor.h.

Referenced by analyze(), bookHistograms(), and PrimaryVertexMonitor().

vertexToken_

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

private

Definition at line 50 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and PrimaryVertexMonitor().

vtxchi2

MonitorElement* PrimaryVertexMonitor::vtxchi2[2]

private

Definition at line 75 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

vtxndf

MonitorElement * PrimaryVertexMonitor::vtxndf[2]

private

Definition at line 75 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

vtxprob

MonitorElement * PrimaryVertexMonitor::vtxprob[2]

private

Definition at line 75 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

weight

MonitorElement * PrimaryVertexMonitor::weight

private

Definition at line 79 of file PrimaryVertexMonitor.h.

Referenced by dirstructure.Directory::calcStats(), dirstructure.Directory::get_fail_rate(), dirstructure.Directory::get_null_rate(), dirstructure.Directory::get_skiped_rate(), dirstructure.Directory::get_success_rate(), and dirstructure.Directory::print_report().

xDiff

MonitorElement * PrimaryVertexMonitor::xDiff[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

xerr

MonitorElement * PrimaryVertexMonitor::xerr[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), python.rootplot.utilities.Hist::av_xerr(), bookHistograms(), python.rootplot.utilities.Hist::delete_bin(), python.rootplot.root2matplotlib.Hist::errorbar(), python.rootplot.root2matplotlib.Hist::errorbarh(), python.rootplot.utilities.Hist::TGraph(), and vertexPlots().

xerrVsTrks

MonitorElement* PrimaryVertexMonitor::xerrVsTrks[2]

private

Definition at line 73 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

xrec

MonitorElement* PrimaryVertexMonitor::xrec[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

yDiff

MonitorElement * PrimaryVertexMonitor::yDiff[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

yerr

MonitorElement * PrimaryVertexMonitor::yerr[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), python.rootplot.utilities.Hist::av_yerr(), bookHistograms(), python.rootplot.utilities.Hist::delete_bin(), python.rootplot.root2matplotlib.Hist::errorbar(), python.rootplot.root2matplotlib.Hist::errorbarh(), python.rootplot.utilities.Hist::min(), python.rootplot.utilities.Hist::scale(), python.rootplot.utilities.Hist::TGraph(), python.rootplot.utilities.Hist::TH1F(), and vertexPlots().

yerrVsTrks

MonitorElement * PrimaryVertexMonitor::yerrVsTrks[2]

private

Definition at line 73 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

yrec

MonitorElement * PrimaryVertexMonitor::yrec[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

zerr

MonitorElement * PrimaryVertexMonitor::zerr[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

zerrVsTrks

MonitorElement * PrimaryVertexMonitor::zerrVsTrks[2]

private

Definition at line 73 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

zrec

MonitorElement * PrimaryVertexMonitor::zrec[2]

private

Definition at line 72 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().