#include <Alignment/GeneralPurposeVertexAnalyzer/plugins/GeneralPurposeVertexAnalyzer.cc>

Inheritance diagram for GeneralPurposeVertexAnalyzer:

Public Member Functions
	GeneralPurposeVertexAnalyzer (const edm::ParameterSet &)

	~GeneralPurposeVertexAnalyzer () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

void	beginJob () override

template<typename T , typename... Args>
T *	book (const Args &...args) const

void	pvTracksPlots (const reco::Vertex &v)

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

Private Attributes
const edm::InputTag	beamSpotInputTag_

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

TH1D *	bsBeamWidthX

TH1D *	bsBeamWidthY

TH1D *	bsDxdz

TH1D *	bsDydz

TH1D *	bsSigmaZ

TH1D *	bsType

TH1D *	bsX

TH1D *	bsY

TH1D *	bsZ

TH1D *	chi2ndf

TH1D *	chi2prob

TH1D *	dxy

TH1D *	dxy2

const int	dxyBin_

TH1D *	dxyErr

const double	dxyMax_

const double	dxyMin_

TProfile *	dxyVsEta_pt1

TProfile *	dxyVsEta_pt10

TProfile2D *	dxyVsEtaVsPhi_pt1

TProfile2D *	dxyVsEtaVsPhi_pt10

TProfile *	dxyVsPhi_pt1

TProfile *	dxyVsPhi_pt10

TH1D *	dz

const int	dzBin_

TH1D *	dzErr

const double	dzMax_

const double	dzMin_

TProfile *	dzVsEta_pt1

TProfile *	dzVsEta_pt10

TProfile2D *	dzVsEtaVsPhi_pt1

TProfile2D *	dzVsEtaVsPhi_pt10

TProfile *	dzVsPhi_pt1

TProfile *	dzVsPhi_pt10

bool	errorPrinted_

TH1D *	eta_pt1

TH1D *	eta_pt10

const int	etaBin2D_

const int	etaBin_

const double	etaMax_

const double	etaMin_

edm::Service< TFileService >	fs_

TH1D *	nans [2]

TH1I *	nbgvtx

TH1D *	nbtksinvtx [2]

TH1I *	nbvtx

const int	ndof_

TH1D *	ntracks

TH1D *	ntracksVsZ [2]

TH1D *	phi_pt1

TH1D *	phi_pt10

const int	phiBin2D_

const int	phiBin_

const double	phiMax_

const double	phiMin_

TH1D *	score [2]

const edm::EDGetTokenT< VertexScore >	scoreToken_

TH1D *	sumpt

const int	tkNoBin_

const double	tkNoMax_

const double	tkNoMin_

TH1D *	trksWeight [2]

TH1D *	tt [2]

TH1D *	type [2]

const edm::InputTag	vertexInputTag_

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

const double	vposx_

const double	vposy_

TH1D *	vtxchi2 [2]

TH1D *	vtxndf [2]

TH1D *	vtxprob [2]

TH1D *	weight

TH1D *	xDiff [2]

TH1D *	xerr [2]

TH2D *	xerrVsTrks [2]

TH1D *	xrec [2]

TH1D *	yDiff [2]

TH1D *	yerr [2]

TH2D *	yerrVsTrks [2]

TH1D *	yrec [2]

TH1D *	zerr [2]

TH2D *	zerrVsTrks [2]

TH1D *	zrec [2]

Static Private Attributes
static constexpr int	cmToUm = 10000

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

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

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

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

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

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

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

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

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: monitor vertex properties for alignment purposes, largely copied from DQMOffline/RecoB/plugins/PrimaryVertexMonitor.cc

Definition at line 55 of file GeneralPurposeVertexAnalyzer.cc.

Member Typedef Documentation

◆ VertexScore

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

private

Definition at line 78 of file GeneralPurposeVertexAnalyzer.cc.

Constructor & Destructor Documentation

◆ GeneralPurposeVertexAnalyzer()

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

explicit

Definition at line 134 of file GeneralPurposeVertexAnalyzer.cc.

References TFileService::kSharedResource.

     : ndof_(iConfig.getParameter<int>("ndof")),
       errorPrinted_(false),
       vertexInputTag_(iConfig.getParameter<edm::InputTag>("vertexLabel")),
       beamSpotInputTag_(iConfig.getParameter<edm::InputTag>("beamSpotLabel")),
       vertexToken_(consumes<reco::VertexCollection>(vertexInputTag_)),
       scoreToken_(consumes<VertexScore>(vertexInputTag_)),
       beamspotToken_(consumes<reco::BeamSpot>(beamSpotInputTag_)),
       // to be configured for each year...
       vposx_(iConfig.getParameter<double>("Xpos")),
       vposy_(iConfig.getParameter<double>("Ypos")),
       tkNoBin_(iConfig.getParameter<int>("TkSizeBin")),
       tkNoMin_(iConfig.getParameter<double>("TkSizeMin")),
       tkNoMax_(iConfig.getParameter<double>("TkSizeMax")),
       dxyBin_(iConfig.getParameter<int>("DxyBin")),
       dxyMin_(iConfig.getParameter<double>("DxyMin")),
       dxyMax_(iConfig.getParameter<double>("DxyMax")),
       dzBin_(iConfig.getParameter<int>("DzBin")),
       dzMin_(iConfig.getParameter<double>("DzMin")),
       dzMax_(iConfig.getParameter<double>("DzMax")),
       phiBin_(iConfig.getParameter<int>("PhiBin")),
       phiBin2D_(iConfig.getParameter<int>("PhiBin2D")),
       phiMin_(iConfig.getParameter<double>("PhiMin")),
       phiMax_(iConfig.getParameter<double>("PhiMax")),
       etaBin_(iConfig.getParameter<int>("EtaBin")),
       etaBin2D_(iConfig.getParameter<int>("EtaBin2D")),
       etaMin_(iConfig.getParameter<double>("EtaMin")),
       etaMax_(iConfig.getParameter<double>("EtaMax")),
       // histograms
       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) {
   usesResource(TFileService::kSharedResource);
 }

◆ ~GeneralPurposeVertexAnalyzer()

GeneralPurposeVertexAnalyzer::~GeneralPurposeVertexAnalyzer ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 207 of file GeneralPurposeVertexAnalyzer.cc.

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

                                                                                               {
   using namespace edm;
 
   const auto &recVtxs = iEvent.getHandle(vertexToken_);
   const auto &scores = iEvent.getHandle(scoreToken_);
   const auto &beamSpotHandle = iEvent.getHandle(beamspotToken_);
   reco::BeamSpot beamSpot = *beamSpotHandle;
 
   // check for absent products and simply "return" in that case
   if (recVtxs.isValid() == false || beamSpotHandle.isValid() == false) {
     edm::LogWarning("GeneralPurposeVertexAnalyzer")
         << " 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("GeneralPurposeVertexAnalyzer")
               << "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;
   }
 
   nbvtx->Fill(recVtxs->size());
   int ng = 0;
   for (auto const &vx : (*recVtxs)) {
     if (vx.isValid() && !vx.isFake() && vx.ndof() >= ndof_) {
       ++ng;
     }
   }
   nbgvtx->Fill(ng);
 
   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());
 }

◆ beginJob()

void GeneralPurposeVertexAnalyzer::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 419 of file GeneralPurposeVertexAnalyzer.cc.

References bsBeamWidthX, bsBeamWidthY, bsDxdz, bsDydz, bsSigmaZ, bsType, bsX, bsY, bsZ, chi2ndf, chi2prob, dxy, dxy2, dxyBin_, dxyErr, dxyMax_, dxyMin_, dxyVsEta_pt1, dxyVsEta_pt10, dxyVsEtaVsPhi_pt1, dxyVsEtaVsPhi_pt10, dxyVsPhi_pt1, dxyVsPhi_pt10, dz, dzBin_, dzErr, dzMax_, dzMin_, dzVsEta_pt1, dzVsEta_pt10, dzVsEtaVsPhi_pt1, dzVsEtaVsPhi_pt10, dzVsPhi_pt1, dzVsPhi_pt10, eta_pt1, eta_pt10, etaBin2D_, etaBin_, etaMax_, etaMin_, mps_fire::i, nans, nbgvtx, nbtksinvtx, nbvtx, ntracks, ntracksVsZ, phi_pt1, phi_pt10, phiBin2D_, phiBin_, phiMax_, phiMin_, score, AlCaHLTBitMon_QueryRunRegistry::string, sumpt, tkNoBin_, tkNoMax_, tkNoMin_, trksWeight, vposx_, vposy_, vtxchi2, vtxndf, vtxprob, xDiff, xerr, xerrVsTrks, xrec, yDiff, yerr, yerrVsTrks, yrec, zerr, zerrVsTrks, and zrec.

                                             {
   nbvtx = book<TH1I>("vtxNbr", "Reconstructed Vertices in Event", 80, -0.5, 79.5);
   nbgvtx = book<TH1I>("goodvtxNbr", "Reconstructed Good Vertices in Event", 80, -0.5, 79.5);
 
   nbtksinvtx[0] = book<TH1D>("otherVtxTrksNbr", "Reconstructed Tracks in Vertex (other Vtx)", 40, -0.5, 99.5);
   ntracksVsZ[0] =
       book<TProfile>("otherVtxTrksVsZ", "Reconstructed Tracks in Vertex (other Vtx) vs Z", 80, -20., 20., 0., 100., "");
   ntracksVsZ[0]->SetXTitle("z-bs");
   ntracksVsZ[0]->SetYTitle("#tracks");
 
   score[0] = book<TH1D>("otherVtxScore", "sqrt(score) (other Vtx)", 100, 0., 400.);
   trksWeight[0] = book<TH1D>("otherVtxTrksWeight", "Total weight of Tracks in Vertex (other Vtx)", 40, 0, 100.);
   vtxchi2[0] = book<TH1D>("otherVtxChi2", "#chi^{2} (other Vtx)", 100, 0., 200.);
   vtxndf[0] = book<TH1D>("otherVtxNdf", "ndof (other Vtx)", 100, 0., 200.);
   vtxprob[0] = book<TH1D>("otherVtxProb", "#chi^{2} probability (other Vtx)", 100, 0., 1.);
   nans[0] = book<TH1D>("otherVtxNans", "Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (other Vtx)", 9, 0.5, 9.5);
 
   nbtksinvtx[1] = book<TH1D>("tagVtxTrksNbr", "Reconstructed Tracks in Vertex (tagged Vtx)", 100, -0.5, 99.5);
   ntracksVsZ[1] =
       book<TProfile>("tagVtxTrksVsZ", "Reconstructed Tracks in Vertex (tagged Vtx) vs Z", 80, -20., 20., 0., 100., "");
   ntracksVsZ[1]->SetXTitle("z-bs");
   ntracksVsZ[1]->SetYTitle("#tracks");
 
   score[1] = book<TH1D>("tagVtxScore", "sqrt(score) (tagged Vtx)", 100, 0., 400.);
   trksWeight[1] = book<TH1D>("tagVtxTrksWeight", "Total weight of Tracks in Vertex (tagged Vtx)", 100, 0, 100.);
   vtxchi2[1] = book<TH1D>("tagVtxChi2", "#chi^{2} (tagged Vtx)", 100, 0., 200.);
   vtxndf[1] = book<TH1D>("tagVtxNdf", "ndof (tagged Vtx)", 100, 0., 200.);
   vtxprob[1] = book<TH1D>("tagVtxProb", "#chi^{2} probability (tagged Vtx)", 100, 0., 1.);
   nans[1] = book<TH1D>("tagVtxNans", "Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (tagged Vtx)", 9, 0.5, 9.5);
 
   xrec[0] = book<TH1D>("otherPosX", "Position x Coordinate (other Vtx)", 100, vposx_ - 0.1, vposx_ + 0.1);
   yrec[0] = book<TH1D>("otherPosY", "Position y Coordinate (other Vtx)", 100, vposy_ - 0.1, vposy_ + 0.1);
   zrec[0] = book<TH1D>("otherPosZ", "Position z Coordinate (other Vtx)", 100, -20., 20.);
   xDiff[0] = book<TH1D>("otherDiffX", "X distance from BeamSpot (other Vtx)", 100, -500, 500);
   yDiff[0] = book<TH1D>("otherDiffY", "Y distance from BeamSpot (other Vtx)", 100, -500, 500);
   xerr[0] = book<TH1D>("otherErrX", "Uncertainty x Coordinate (other Vtx)", 100, 0., 100);
   yerr[0] = book<TH1D>("otherErrY", "Uncertainty y Coordinate (other Vtx)", 100, 0., 100);
   zerr[0] = book<TH1D>("otherErrZ", "Uncertainty z Coordinate (other Vtx)", 100, 0., 100);
   xerrVsTrks[0] = book<TH2D>(
       "otherErrVsWeightX", "Uncertainty x Coordinate vs. track weight (other Vtx)", 100, 0, 100., 100, 0., 100);
   yerrVsTrks[0] = book<TH2D>(
       "otherErrVsWeightY", "Uncertainty y Coordinate vs. track weight (other Vtx)", 100, 0, 100., 100, 0., 100);
   zerrVsTrks[0] = book<TH2D>(
       "otherErrVsWeightZ", "Uncertainty z Coordinate vs. track weight (other Vtx)", 100, 0, 100., 100, 0., 100);
 
   xrec[1] = book<TH1D>("tagPosX", "Position x Coordinate (tagged Vtx)", 100, vposx_ - 0.1, vposx_ + 0.1);
   yrec[1] = book<TH1D>("tagPosY", "Position y Coordinate (tagged Vtx)", 100, vposy_ - 0.1, vposy_ + 0.1);
   zrec[1] = book<TH1D>("tagPosZ", "Position z Coordinate (tagged Vtx)", 100, -20., 20.);
   xDiff[1] = book<TH1D>("tagDiffX", "X distance from BeamSpot (tagged Vtx)", 100, -500, 500);
   yDiff[1] = book<TH1D>("tagDiffY", "Y distance from BeamSpot (tagged Vtx)", 100, -500, 500);
   xerr[1] = book<TH1D>("tagErrX", "Uncertainty x Coordinate (tagged Vtx)", 100, 0., 100);
   yerr[1] = book<TH1D>("tagErrY", "Uncertainty y Coordinate (tagged Vtx)", 100, 0., 100);
   zerr[1] = book<TH1D>("tagErrZ", "Uncertainty z Coordinate (tagged Vtx)", 100, 0., 100);
   xerrVsTrks[1] = book<TH2D>(
       "tagErrVsWeightX", "Uncertainty x Coordinate vs. track weight (tagged Vtx)", 100, 0, 100., 100, 0., 100);
   yerrVsTrks[1] = book<TH2D>(
       "tagErrVsWeightY", "Uncertainty y Coordinate vs. track weight (tagged Vtx)", 100, 0, 100., 100, 0., 100);
   zerrVsTrks[1] = book<TH2D>(
       "tagErrVsWeightZ", "Uncertainty z Coordinate vs. track weight (tagged Vtx)", 100, 0, 100., 100, 0., 100);
 
   type[0] = book<TH1D>("otherType", "Vertex type (other Vtx)", 3, -0.5, 2.5);
   type[1] = book<TH1D>("tagType", "Vertex type (tagged Vtx)", 3, -0.5, 2.5);
 
   for (int i = 0; i < 2; ++i) {
     type[i]->GetXaxis()->SetBinLabel(1, "Valid, real");
     type[i]->GetXaxis()->SetBinLabel(2, "Valid, fake");
     type[i]->GetXaxis()->SetBinLabel(3, "Invalid");
   }
 
   bsX = book<TH1D>("bsX", "BeamSpot x0", 100, vposx_ - 0.1, vposx_ + 0.1);
   bsY = book<TH1D>("bsY", "BeamSpot y0", 100, vposy_ - 0.1, vposy_ + 0.1);
   bsZ = book<TH1D>("bsZ", "BeamSpot z0", 100, -2., 2.);
   bsSigmaZ = book<TH1D>("bsSigmaZ", "BeamSpot sigmaZ", 100, 0., 10.);
   bsDxdz = book<TH1D>("bsDxdz", "BeamSpot dxdz", 100, -0.0003, 0.0003);
   bsDydz = book<TH1D>("bsDydz", "BeamSpot dydz", 100, -0.0003, 0.0003);
   bsBeamWidthX = book<TH1D>("bsBeamWidthX", "BeamSpot BeamWidthX", 100, 0., 100.);
   bsBeamWidthY = book<TH1D>("bsBeamWidthY", "BeamSpot BeamWidthY", 100, 0., 100.);
   bsType = book<TH1D>("bsType", "BeamSpot type", 4, -1.5, 2.5);
   bsType->GetXaxis()->SetBinLabel(1, "Unknown");
   bsType->GetXaxis()->SetBinLabel(2, "Fake");
   bsType->GetXaxis()->SetBinLabel(3, "LHC");
   bsType->GetXaxis()->SetBinLabel(4, "Tracker");
 
   // repeated strings in titles
   std::string s_1 = "PV Tracks (p_{T} > 1 GeV)";
   std::string s_10 = "PV Tracks (p_{T} > 10 GeV)";
 
   ntracks = book<TH1D>("ntracks", fmt::sprintf("number of %s", s_1).c_str(), tkNoBin_, tkNoMin_, tkNoMax_);
   ntracks->SetXTitle(fmt::sprintf("Number of %s per Event", s_1).c_str());
   ntracks->SetYTitle("Number of Events");
 
   weight = book<TH1D>("weight", fmt::sprintf("weight of %s", s_1).c_str(), 100, 0., 1.);
   weight->SetXTitle(fmt::sprintf("weight of %s per Event", s_1).c_str());
   weight->SetYTitle("Number of Events");
 
   sumpt = book<TH1D>("sumpt", fmt::sprintf("#Sum p_{T} of %s", s_1).c_str(), 100, -0.5, 249.5);
   chi2ndf = book<TH1D>("chi2ndf", fmt::sprintf("%s #chi^{2}/ndof", s_1).c_str(), 100, 0., 20.);
   chi2prob = book<TH1D>("chi2prob", fmt::sprintf("%s #chi^{2} probability", s_1).c_str(), 100, 0., 1.);
   dxy = book<TH1D>("dxy", fmt::sprintf("%s d_{xy} (#mum)", s_1).c_str(), dxyBin_, dxyMin_, dxyMax_);
   dxy2 = book<TH1D>("dxyzoom", fmt::sprintf("%s d_{xy} (#mum)", s_1).c_str(), dxyBin_, dxyMin_ / 5., dxyMax_ / 5.);
   dxyErr = book<TH1D>("dxyErr", fmt::sprintf("%s d_{xy} error (#mum)", s_1).c_str(), 100, 0., 2000.);
   dz = book<TH1D>("dz", fmt::sprintf("%s d_{z} (#mum)", s_1).c_str(), dzBin_, dzMin_, dzMax_);
   dzErr = book<TH1D>("dzErr", fmt::sprintf("%s d_{z} error(#mum)", s_1).c_str(), 100, 0., 10000.);
 
   phi_pt1 =
       book<TH1D>("phi_pt1", fmt::sprintf("%s #phi; PV tracks #phi;#tracks", s_1).c_str(), phiBin_, phiMin_, phiMax_);
   eta_pt1 =
       book<TH1D>("eta_pt1", fmt::sprintf("%s #eta; PV tracks #eta;#tracks", s_1).c_str(), etaBin_, etaMin_, etaMax_);
   phi_pt10 =
       book<TH1D>("phi_pt10", fmt::sprintf("%s #phi; PV tracks #phi;#tracks", s_10).c_str(), phiBin_, phiMin_, phiMax_);
   eta_pt10 =
       book<TH1D>("eta_pt10", fmt::sprintf("%s #phi; PV tracks #eta;#tracks", s_10).c_str(), etaBin_, etaMin_, etaMax_);
 
   dxyVsPhi_pt1 = book<TProfile>("dxyVsPhi_pt1",
                                 fmt::sprintf("%s d_{xy} (#mum) VS track #phi", s_1).c_str(),
                                 phiBin_,
                                 phiMin_,
                                 phiMax_,
                                 dxyMin_,
                                 dxyMax_);
   dxyVsPhi_pt1->SetXTitle("PV track (p_{T} > 1 GeV) #phi");
   dxyVsPhi_pt1->SetYTitle("PV track (p_{T} > 1 GeV) d_{xy} (#mum)");
 
   dzVsPhi_pt1 = book<TProfile>("dzVsPhi_pt1",
                                fmt::sprintf("%s d_{z} (#mum) VS track #phi", s_1).c_str(),
                                phiBin_,
                                phiMin_,
                                phiMax_,
                                dzMin_,
                                dzMax_);
   dzVsPhi_pt1->SetXTitle("PV track (p_{T} > 1 GeV) #phi");
   dzVsPhi_pt1->SetYTitle("PV track (p_{T} > 1 GeV) d_{z} (#mum)");
 
   dxyVsEta_pt1 = book<TProfile>("dxyVsEta_pt1",
                                 fmt::sprintf("%s d_{xy} (#mum) VS track #eta", s_1).c_str(),
                                 etaBin_,
                                 etaMin_,
                                 etaMax_,
                                 dxyMin_,
                                 dxyMax_);
   dxyVsEta_pt1->SetXTitle("PV track (p_{T} > 1 GeV) #eta");
   dxyVsEta_pt1->SetYTitle("PV track (p_{T} > 1 GeV) d_{xy} (#mum)");
 
   dzVsEta_pt1 = book<TProfile>("dzVsEta_pt1",
                                fmt::sprintf("%s d_{z} (#mum) VS track #eta", s_1).c_str(),
                                etaBin_,
                                etaMin_,
                                etaMax_,
                                dzMin_,
                                dzMax_);
   dzVsEta_pt1->SetXTitle("PV track (p_{T} > 1 GeV) #eta");
   dzVsEta_pt1->SetYTitle("PV track (p_{T} > 1 GeV) d_{z} (#mum)");
 
   dxyVsEtaVsPhi_pt1 = book<TProfile2D>("dxyVsEtaVsPhi_pt1",
                                        fmt::sprintf("%s d_{xy} (#mum) VS track #eta VS track #phi", s_1).c_str(),
                                        etaBin2D_,
                                        etaMin_,
                                        etaMax_,
                                        phiBin2D_,
                                        phiMin_,
                                        phiMax_,
                                        dxyMin_,
                                        dxyMax_);
   dxyVsEtaVsPhi_pt1->SetXTitle("PV track (p_{T} > 1 GeV) #eta");
   dxyVsEtaVsPhi_pt1->SetYTitle("PV track (p_{T} > 1 GeV) #phi");
   dxyVsEtaVsPhi_pt1->SetZTitle("PV track (p_{T} > 1 GeV) d_{xy} (#mum)");
 
   dzVsEtaVsPhi_pt1 = book<TProfile2D>("dzVsEtaVsPhi_pt1",
                                       fmt::sprintf("%s d_{z} (#mum) VS track #eta VS track #phi", s_1).c_str(),
                                       etaBin2D_,
                                       etaMin_,
                                       etaMax_,
                                       phiBin2D_,
                                       phiMin_,
                                       phiMax_,
                                       dzMin_,
                                       dzMax_);
   dzVsEtaVsPhi_pt1->SetXTitle("PV track (p_{T} > 1 GeV) #eta");
   dzVsEtaVsPhi_pt1->SetYTitle("PV track (p_{T} > 1 GeV) #phi");
   dzVsEtaVsPhi_pt1->SetZTitle("PV track (p_{T} > 1 GeV) d_{z} (#mum)");
 
   dxyVsPhi_pt10 = book<TProfile>("dxyVsPhi_pt10",
                                  fmt::sprintf("%s d_{xy} (#mum) VS track #phi", s_10).c_str(),
                                  phiBin_,
                                  phiMin_,
                                  phiMax_,
                                  dxyMin_,
                                  dxyMax_);
   dxyVsPhi_pt10->SetXTitle("PV track (p_{T} > 10 GeV) #phi");
   dxyVsPhi_pt10->SetYTitle("PV track (p_{T} > 10 GeV) d_{xy} (#mum)");
 
   dzVsPhi_pt10 = book<TProfile>("dzVsPhi_pt10",
                                 fmt::sprintf("%s d_{z} (#mum) VS track #phi", s_10).c_str(),
                                 phiBin_,
                                 phiMin_,
                                 phiMax_,
                                 dzMin_,
                                 dzMax_);
   dzVsPhi_pt10->SetXTitle("PV track (p_{T} > 10 GeV) #phi");
   dzVsPhi_pt10->SetYTitle("PV track (p_{T} > 10 GeV) d_{z} (#mum)");
 
   dxyVsEta_pt10 = book<TProfile>("dxyVsEta_pt10",
                                  fmt::sprintf("%s d_{xy} (#mum) VS track #eta", s_10).c_str(),
                                  etaBin_,
                                  etaMin_,
                                  etaMax_,
                                  dxyMin_,
                                  dxyMax_);
   dxyVsEta_pt10->SetXTitle("PV track (p_{T} > 10 GeV) #eta");
   dxyVsEta_pt10->SetYTitle("PV track (p_{T} > 10 GeV) d_{xy} (#mum)");
 
   dzVsEta_pt10 = book<TProfile>("dzVsEta_pt10",
                                 fmt::sprintf("%s d_{z} (#mum) VS track #eta", s_10).c_str(),
                                 etaBin_,
                                 etaMin_,
                                 etaMax_,
                                 dzMin_,
                                 dzMax_);
   dzVsEta_pt10->SetXTitle("PV track (p_{T} > 10 GeV) #eta");
   dzVsEta_pt10->SetYTitle("PV track (p_{T} > 10 GeV) d_{z} (#mum)");
 
   dxyVsEtaVsPhi_pt10 = book<TProfile2D>("dxyVsEtaVsPhi_pt10",
                                         fmt::sprintf("%s d_{xy} (#mum) VS track #eta VS track #phi", s_10).c_str(),
                                         etaBin2D_,
                                         etaMin_,
                                         etaMax_,
                                         phiBin2D_,
                                         phiMin_,
                                         phiMax_,
                                         dxyMin_,
                                         dxyMax_);
   dxyVsEtaVsPhi_pt10->SetXTitle("PV track (p_{T} > 10 GeV) #eta");
   dxyVsEtaVsPhi_pt10->SetYTitle("PV track (p_{T} > 10 GeV) #phi");
   dxyVsEtaVsPhi_pt10->SetZTitle("PV track (p_{T} > 10 GeV) d_{xy} (#mum)");
 
   dzVsEtaVsPhi_pt10 = book<TProfile2D>("dzVsEtaVsPhi_pt10",
                                        fmt::sprintf("%s d_{z} (#mum) VS track #eta VS track #phi", s_10).c_str(),
                                        etaBin2D_,
                                        etaMin_,
                                        etaMax_,
                                        phiBin2D_,
                                        phiMin_,
                                        phiMax_,
                                        dzMin_,
                                        dzMax_);
   dzVsEtaVsPhi_pt10->SetXTitle("PV track (p_{T} > 10 GeV) #eta");
   dzVsEtaVsPhi_pt10->SetYTitle("PV track (p_{T} > 10 GeV) #phi");
   dzVsEtaVsPhi_pt10->SetZTitle("PV track (p_{T} > 10 GeV) d_{z} (#mum)");
 }

◆ book()

template<typename T , typename... Args>

T * GeneralPurposeVertexAnalyzer::book ( const Args &... args ) const

private

Definition at line 413 of file GeneralPurposeVertexAnalyzer.cc.

References writedatasetfile::args, fs_, TFileService::make(), and submitPVValidationJobs::t.

                                                                {
   T *t = fs_->make<T>(args...);
   return t;
 }

◆ fillDescriptions()

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

static

Definition at line 670 of file GeneralPurposeVertexAnalyzer.cc.

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

                                                                                               {
   edm::ParameterSetDescription desc;
   desc.add<int>("ndof", 4);
   desc.add<edm::InputTag>("vertexLabel", edm::InputTag("offlinePrimaryVertices"));
   desc.add<edm::InputTag>("beamSpotLabel", edm::InputTag("offlineBeamSpot"));
   desc.add<double>("Xpos", 0.1);
   desc.add<double>("Ypos", 0.0);
   desc.add<int>("TkSizeBin", 100);
   desc.add<double>("TkSizeMin", 499.5);
   desc.add<double>("TkSizeMax", -0.5);
   desc.add<int>("DxyBin", 100);
   desc.add<double>("DxyMin", 5000.);
   desc.add<double>("DxyMax", -5000.);
   desc.add<int>("DzBin", 100);
   desc.add<double>("DzMin", -2000.0);
   desc.add<double>("DzMax", 2000.0);
   desc.add<int>("PhiBin", 32);
   desc.add<int>("PhiBin2D", 12);
   desc.add<double>("PhiMin", -M_PI);
   desc.add<double>("PhiMax", M_PI);
   desc.add<int>("EtaBin", 26);
   desc.add<int>("EtaBin2D", 8);
   desc.add<double>("EtaMin", -2.7);
   desc.add<double>("EtaMax", 2.7);
   descriptions.addWithDefaultLabel(desc);
 }

◆ pvTracksPlots()

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

private

Definition at line 285 of file GeneralPurposeVertexAnalyzer.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, f, reco::TrackBase::highPurity, ntracks, phi, phi_pt1, phi_pt10, DiDispStaMuonMonitor_cfi::pt, HLT_2023v12_cff::pt2, sumpt, submitPVValidationJobs::t, 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());
 
   float sumPT = 0.;
   for (const auto &t : v.tracks()) {
     const bool isHighPurity = t->quality(reco::TrackBase::highPurity);
     if (!isHighPurity) {
       continue;
     }
 
     const float pt = t->pt();
     if (pt < 1.f) {
       continue;
     }
 
     const float pt2 = pt * pt;
     const float eta = t->eta();
     const float phi = t->phi();
     const float w = v.trackWeight(t);
     const float chi2NDF = t->normalizedChi2();
     const float chi2Prob = TMath::Prob(t->chi2(), static_cast<int>(t->ndof()));
     const float Dxy = t->dxy(myVertex) * cmToUm;
     const float Dz = t->dz(myVertex) * cmToUm;
     const float DxyErr = t->dxyError() * cmToUm;
     const float DzErr = t->dzError() * cmToUm;
 
     sumPT += pt2;
 
     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.f) {
       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(static_cast<float>(v.tracks().size()));
   sumpt->Fill(sumPT);
 }

◆ vertexPlots()

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

private

Definition at line 356 of file GeneralPurposeVertexAnalyzer.cc.

References pwdgSkimBPark_cfi::beamSpot, ChiSquaredProbability(), cmToUm, 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

◆ beamSpotInputTag_

const edm::InputTag GeneralPurposeVertexAnalyzer::beamSpotInputTag_

private

Definition at line 76 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ beamspotToken_

const edm::EDGetTokenT<reco::BeamSpot> GeneralPurposeVertexAnalyzer::beamspotToken_

private

Definition at line 80 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ bsBeamWidthX

TH1D * GeneralPurposeVertexAnalyzer::bsBeamWidthX

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsBeamWidthY

TH1D * GeneralPurposeVertexAnalyzer::bsBeamWidthY

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsDxdz

TH1D * GeneralPurposeVertexAnalyzer::bsDxdz

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsDydz

TH1D * GeneralPurposeVertexAnalyzer::bsDydz

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsSigmaZ

TH1D * GeneralPurposeVertexAnalyzer::bsSigmaZ

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsType

TH1D * GeneralPurposeVertexAnalyzer::bsType

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsX

TH1D* GeneralPurposeVertexAnalyzer::bsX

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsY

TH1D * GeneralPurposeVertexAnalyzer::bsY

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ bsZ

TH1D * GeneralPurposeVertexAnalyzer::bsZ

private

Definition at line 117 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ chi2ndf

TH1D * GeneralPurposeVertexAnalyzer::chi2ndf

private

Definition at line 119 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ chi2prob

TH1D * GeneralPurposeVertexAnalyzer::chi2prob

private

Definition at line 119 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ cmToUm

constexpr int GeneralPurposeVertexAnalyzer::cmToUm = 10000

staticprivate

Definition at line 82 of file GeneralPurposeVertexAnalyzer.cc.

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

◆ dxy

TH1D* GeneralPurposeVertexAnalyzer::dxy

private

Definition at line 120 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), Electron.Electron::cutBasedId(), ntupleDataFormat.Track::dxyPull(), and pvTracksPlots().

◆ dxy2

TH1D * GeneralPurposeVertexAnalyzer::dxy2

private

Definition at line 120 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dxyBin_

const int GeneralPurposeVertexAnalyzer::dxyBin_

private

Definition at line 90 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ dxyErr

TH1D * GeneralPurposeVertexAnalyzer::dxyErr

private

Definition at line 120 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), ntupleDataFormat.Track::dxyPull(), and pvTracksPlots().

◆ dxyMax_

const double GeneralPurposeVertexAnalyzer::dxyMax_

private

Definition at line 92 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ dxyMin_

const double GeneralPurposeVertexAnalyzer::dxyMin_

private

Definition at line 91 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ dxyVsEta_pt1

TProfile* GeneralPurposeVertexAnalyzer::dxyVsEta_pt1

private

Definition at line 124 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dxyVsEta_pt10

TProfile* GeneralPurposeVertexAnalyzer::dxyVsEta_pt10

private

Definition at line 127 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dxyVsEtaVsPhi_pt1

TProfile2D* GeneralPurposeVertexAnalyzer::dxyVsEtaVsPhi_pt1

private

Definition at line 125 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dxyVsEtaVsPhi_pt10

TProfile2D* GeneralPurposeVertexAnalyzer::dxyVsEtaVsPhi_pt10

private

Definition at line 128 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dxyVsPhi_pt1

TProfile* GeneralPurposeVertexAnalyzer::dxyVsPhi_pt1

private

Definition at line 123 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dxyVsPhi_pt10

TProfile* GeneralPurposeVertexAnalyzer::dxyVsPhi_pt10

private

Definition at line 126 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dz

TH1D * GeneralPurposeVertexAnalyzer::dz

private

Definition at line 120 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), Electron.Electron::cutBasedId(), ntupleDataFormat.Track::dzPull(), and pvTracksPlots().

◆ dzBin_

const int GeneralPurposeVertexAnalyzer::dzBin_

private

Definition at line 94 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ dzErr

TH1D * GeneralPurposeVertexAnalyzer::dzErr

private

Definition at line 120 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), ntupleDataFormat.Track::dzPull(), and pvTracksPlots().

◆ dzMax_

const double GeneralPurposeVertexAnalyzer::dzMax_

private

Definition at line 96 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ dzMin_

const double GeneralPurposeVertexAnalyzer::dzMin_

private

Definition at line 95 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ dzVsEta_pt1

TProfile * GeneralPurposeVertexAnalyzer::dzVsEta_pt1

private

Definition at line 124 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dzVsEta_pt10

TProfile * GeneralPurposeVertexAnalyzer::dzVsEta_pt10

private

Definition at line 127 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dzVsEtaVsPhi_pt1

TProfile2D * GeneralPurposeVertexAnalyzer::dzVsEtaVsPhi_pt1

private

Definition at line 125 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dzVsEtaVsPhi_pt10

TProfile2D * GeneralPurposeVertexAnalyzer::dzVsEtaVsPhi_pt10

private

Definition at line 128 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dzVsPhi_pt1

TProfile * GeneralPurposeVertexAnalyzer::dzVsPhi_pt1

private

Definition at line 123 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ dzVsPhi_pt10

TProfile * GeneralPurposeVertexAnalyzer::dzVsPhi_pt10

private

Definition at line 126 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ errorPrinted_

bool GeneralPurposeVertexAnalyzer::errorPrinted_

private

Definition at line 74 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ eta_pt1

TH1D * GeneralPurposeVertexAnalyzer::eta_pt1

private

Definition at line 121 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ eta_pt10

TH1D * GeneralPurposeVertexAnalyzer::eta_pt10

private

Definition at line 122 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ etaBin2D_

const int GeneralPurposeVertexAnalyzer::etaBin2D_

private

Definition at line 104 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ etaBin_

const int GeneralPurposeVertexAnalyzer::etaBin_

private

Definition at line 103 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ etaMax_

const double GeneralPurposeVertexAnalyzer::etaMax_

private

Definition at line 106 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ etaMin_

const double GeneralPurposeVertexAnalyzer::etaMin_

private

Definition at line 105 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ fs_

edm::Service<TFileService> GeneralPurposeVertexAnalyzer::fs_

private

Definition at line 71 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by book().

◆ nans

TH1D * GeneralPurposeVertexAnalyzer::nans[2]

private

Definition at line 115 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ nbgvtx

TH1I * GeneralPurposeVertexAnalyzer::nbgvtx

private

Definition at line 109 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ nbtksinvtx

TH1D* GeneralPurposeVertexAnalyzer::nbtksinvtx[2]

private

Definition at line 110 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ nbvtx

TH1I* GeneralPurposeVertexAnalyzer::nbvtx

private

Definition at line 109 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ ndof_

const int GeneralPurposeVertexAnalyzer::ndof_

private

Definition at line 73 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ ntracks

TH1D * GeneralPurposeVertexAnalyzer::ntracks

private

Definition at line 119 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ ntracksVsZ

TH1D* GeneralPurposeVertexAnalyzer::ntracksVsZ[2]

private

Definition at line 114 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ phi_pt1

TH1D* GeneralPurposeVertexAnalyzer::phi_pt1

private

Definition at line 121 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ phi_pt10

TH1D* GeneralPurposeVertexAnalyzer::phi_pt10

private

Definition at line 122 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and pvTracksPlots().

◆ phiBin2D_

const int GeneralPurposeVertexAnalyzer::phiBin2D_

private

Definition at line 99 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ phiBin_

const int GeneralPurposeVertexAnalyzer::phiBin_

private

Definition at line 98 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ phiMax_

const double GeneralPurposeVertexAnalyzer::phiMax_

private

Definition at line 101 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ phiMin_

const double GeneralPurposeVertexAnalyzer::phiMin_

private

Definition at line 100 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ score

TH1D * GeneralPurposeVertexAnalyzer::score[2]

private

Definition at line 110 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ scoreToken_

const edm::EDGetTokenT<VertexScore> GeneralPurposeVertexAnalyzer::scoreToken_

private

Definition at line 79 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ sumpt

TH1D* GeneralPurposeVertexAnalyzer::sumpt

private

Definition at line 119 of file GeneralPurposeVertexAnalyzer.cc.

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

◆ tkNoBin_

const int GeneralPurposeVertexAnalyzer::tkNoBin_

private

Definition at line 86 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ tkNoMax_

const double GeneralPurposeVertexAnalyzer::tkNoMax_

private

Definition at line 88 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ tkNoMin_

const double GeneralPurposeVertexAnalyzer::tkNoMin_

private

Definition at line 87 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ trksWeight

TH1D * GeneralPurposeVertexAnalyzer::trksWeight[2]

private

Definition at line 110 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ tt

TH1D* GeneralPurposeVertexAnalyzer::tt[2]

private

Definition at line 111 of file GeneralPurposeVertexAnalyzer.cc.

◆ type

TH1D* GeneralPurposeVertexAnalyzer::type[2]

private

Definition at line 116 of file GeneralPurposeVertexAnalyzer.cc.

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

◆ vertexInputTag_

const edm::InputTag GeneralPurposeVertexAnalyzer::vertexInputTag_

private

Definition at line 76 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ vertexToken_

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

private

Definition at line 77 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by analyze().

◆ vposx_

const double GeneralPurposeVertexAnalyzer::vposx_

private

Definition at line 84 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ vposy_

const double GeneralPurposeVertexAnalyzer::vposy_

private

Definition at line 85 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob().

◆ vtxchi2

TH1D* GeneralPurposeVertexAnalyzer::vtxchi2[2]

private

Definition at line 115 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ vtxndf

TH1D * GeneralPurposeVertexAnalyzer::vtxndf[2]

private

Definition at line 115 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ vtxprob

TH1D * GeneralPurposeVertexAnalyzer::vtxprob[2]

private

Definition at line 115 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ weight

TH1D * GeneralPurposeVertexAnalyzer::weight

private

Definition at line 119 of file GeneralPurposeVertexAnalyzer.cc.

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

TH1D * GeneralPurposeVertexAnalyzer::xDiff[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ xerr

TH1D * GeneralPurposeVertexAnalyzer::xerr[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

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

◆ xerrVsTrks

TH2D* GeneralPurposeVertexAnalyzer::xerrVsTrks[2]

private

Definition at line 113 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ xrec

TH1D* GeneralPurposeVertexAnalyzer::xrec[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ yDiff

TH1D * GeneralPurposeVertexAnalyzer::yDiff[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ yerr

TH1D * GeneralPurposeVertexAnalyzer::yerr[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), python.rootplot.utilities.Hist::av_yerr(), beginJob(), 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

TH2D * GeneralPurposeVertexAnalyzer::yerrVsTrks[2]

private

Definition at line 113 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ yrec

TH1D * GeneralPurposeVertexAnalyzer::yrec[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ zerr

TH1D * GeneralPurposeVertexAnalyzer::zerr[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ zerrVsTrks

TH2D * GeneralPurposeVertexAnalyzer::zerrVsTrks[2]

private

Definition at line 113 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

◆ zrec

TH1D * GeneralPurposeVertexAnalyzer::zrec[2]

private

Definition at line 112 of file GeneralPurposeVertexAnalyzer.cc.

Referenced by beginJob(), and vertexPlots().

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ VertexScore

Constructor & Destructor Documentation

◆ GeneralPurposeVertexAnalyzer()

◆ ~GeneralPurposeVertexAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ book()

◆ fillDescriptions()

◆ pvTracksPlots()

◆ vertexPlots()

Member Data Documentation

◆ beamSpotInputTag_

◆ beamspotToken_

◆ bsBeamWidthX

◆ bsBeamWidthY

◆ bsDxdz

◆ bsDydz

◆ bsSigmaZ

◆ bsType

◆ bsX

◆ bsY

◆ bsZ

◆ chi2ndf

◆ chi2prob

◆ cmToUm

◆ dxy

◆ dxy2

◆ dxyBin_

◆ dxyErr

◆ dxyMax_

◆ dxyMin_

◆ dxyVsEta_pt1

◆ dxyVsEta_pt10

◆ dxyVsEtaVsPhi_pt1

◆ dxyVsEtaVsPhi_pt10

◆ dxyVsPhi_pt1

◆ dxyVsPhi_pt10

◆ dz

◆ dzBin_

◆ dzErr

◆ dzMax_

◆ dzMin_

◆ dzVsEta_pt1

◆ dzVsEta_pt10

◆ dzVsEtaVsPhi_pt1

◆ dzVsEtaVsPhi_pt10

◆ dzVsPhi_pt1

◆ dzVsPhi_pt10

◆ errorPrinted_

◆ eta_pt1

◆ eta_pt10

◆ etaBin2D_

◆ etaBin_

◆ etaMax_

◆ etaMin_

◆ fs_

◆ nans

◆ nbgvtx

◆ nbtksinvtx

◆ nbvtx

◆ ndof_

◆ ntracks

◆ ntracksVsZ

◆ phi_pt1

◆ phi_pt10

◆ phiBin2D_

◆ phiBin_

◆ phiMax_

◆ phiMin_

◆ score