#include <PrimaryVertexMonitor.h>

Inheritance diagram for PrimaryVertexMonitor:

Public Member Functions
virtual void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)

virtual void	endJob ()

	PrimaryVertexMonitor (const edm::ParameterSet &pSet)

	~PrimaryVertexMonitor ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Private Attributes
edm::InputTag	beamSpotLabel

MonitorElement *	bsBeamWidthX

MonitorElement *	bsBeamWidthY

MonitorElement *	bsDxdz

MonitorElement *	bsDydz

MonitorElement *	bsSigmaZ

MonitorElement *	bsType

MonitorElement *	bsX

MonitorElement *	bsY

MonitorElement *	bsZ

std::string	dqmLabel

DQMStore *	dqmStore_

edm::InputTag	moduleLabel

MonitorElement *	nans [2]

MonitorElement *	nbtksinvtx [2]

MonitorElement *	nbvtx

MonitorElement *	trksWeight [2]

MonitorElement *	tt [2]

MonitorElement *	type [2]

MonitorElement *	vtxchi2 [2]

MonitorElement *	vtxndf [2]

MonitorElement *	vtxprob [2]

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]

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

CurrentProcessingContext const *	currentContext () const

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

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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)

Detailed Description

Definition at line 23 of file PrimaryVertexMonitor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file PrimaryVertexMonitor.cc.

References edm::ParameterSet::getParameter(), i, cppFunctionSkipper::operator, and spu::zerr().

 {
   moduleLabel = pSet.getParameter<InputTag>("vertexLabel");
   beamSpotLabel = pSet.getParameter<InputTag>("beamSpotLabel");
 
   //
   // Book all histograms.
   //
 
   //  get the store
   dqmStore_ = edm::Service<DQMStore>().operator->();
   dqmLabel = "OfflinePV/"+moduleLabel.label();
   dqmStore_->setCurrentFolder(dqmLabel);
 
 //   xPos = dqmStore_->book1D ("xPos","x Coordinate" ,100, -0.1, 0.1);
 
   nbvtx      = dqmStore_->book1D("vtxNbr","Reconstructed Vertices in Event",50,-0.5,49.5);
 
   nbtksinvtx[0] = dqmStore_->book1D("otherVtxTrksNbr","Reconstructed Tracks in Vertex (other Vtx)",40,-0.5,99.5); 
   trksWeight[0] = dqmStore_->book1D("otherVtxTrksWeight","Total weight of Tracks in Vertex (other Vtx)",40,0,100.); 
   vtxchi2[0]    = dqmStore_->book1D("otherVtxChi2","#chi^{2} (other Vtx)",100,0.,200.);
   vtxndf[0]     = dqmStore_->book1D("otherVtxNdf","ndof (other Vtx)",100,0.,200.);
   vtxprob[0]    = dqmStore_->book1D("otherVtxProb","#chi^{2} probability (other Vtx)",100,0.,1.);
   nans[0]       = dqmStore_->book1D("otherVtxNans","Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (other Vtx)",9,0.5,9.5);
 
   nbtksinvtx[1] = dqmStore_->book1D("tagVtxTrksNbr","Reconstructed Tracks in Vertex (tagged Vtx)",100,-0.5,99.5); 
   trksWeight[1] = dqmStore_->book1D("tagVtxTrksWeight","Total weight of Tracks in Vertex (tagged Vtx)",100,0,100.); 
   vtxchi2[1]    = dqmStore_->book1D("tagVtxChi2","#chi^{2} (tagged Vtx)",100,0.,200.);
   vtxndf[1]     = dqmStore_->book1D("tagVtxNdf","ndof (tagged Vtx)",100,0.,200.);
   vtxprob[1]    = dqmStore_->book1D("tagVtxProb","#chi^{2} probability (tagged Vtx)",100,0.,1.);
   nans[1]       = dqmStore_->book1D("tagVtxNans","Illegal values for x,y,z,xx,xy,xz,yy,yz,zz (tagged Vtx)",9,0.5,9.5);
 
   xrec[0]    = dqmStore_->book1D("otherPosX","Position x Coordinate (other Vtx)",100,-0.1,0.1);
   yrec[0]    = dqmStore_->book1D("otherPosY","Position y Coordinate (other Vtx)",100,-0.1,0.1);
   zrec[0]        = dqmStore_->book1D("otherPosZ","Position z Coordinate (other Vtx)",100,-20.,20.);
   xDiff[0]   = dqmStore_->book1D("otherDiffX","X distance from BeamSpot (other Vtx)",100,-500,500);
   yDiff[0]   = dqmStore_->book1D("otherDiffY","Y distance from BeamSpot (other Vtx)",100,-500,500);
   xerr[0]    = dqmStore_->book1D("otherErrX","Uncertainty x Coordinate (other Vtx)",100,-0.1,0.1);
   yerr[0]    = dqmStore_->book1D("otherErrY","Uncertainty y Coordinate (other Vtx)",100,-0.1,0.1);
   zerr[0]        = dqmStore_->book1D("otherErrZ","Uncertainty z Coordinate (other Vtx)",100,-20.,20.);
   xerrVsTrks[0]  = dqmStore_->book2D("otherErrVsWeightX","Uncertainty x Coordinate vs. track weight (other Vtx)",100,0,100.,100,-0.1,0.1);
   yerrVsTrks[0]  = dqmStore_->book2D("otherErrVsWeightY","Uncertainty y Coordinate vs. track weight (other Vtx)",100,0,100.,100,-0.1,0.1);
   zerrVsTrks[0]  = dqmStore_->book2D("otherErrVsWeightZ","Uncertainty z Coordinate vs. track weight (other Vtx)",100,0,100.,100,-0.1,0.1);
 
 
   xrec[1]     = dqmStore_->book1D("tagPosX","Position x Coordinate (tagged Vtx)",100,-0.1,0.1);
   yrec[1]     = dqmStore_->book1D("tagPosY","Position y Coordinate (tagged Vtx)",100,-0.1,0.1);
   zrec[1]     = dqmStore_->book1D("tagPosZ","Position z Coordinate (tagged Vtx)",100,-20.,20.);
   xDiff[1]    = dqmStore_->book1D("tagDiffX","X distance from BeamSpot (tagged Vtx)",100,-500, 500);
   yDiff[1]    = dqmStore_->book1D("tagDiffY","Y distance from BeamSpot (tagged Vtx)",100,-500, 500);
   xerr[1]     = dqmStore_->book1D("tagErrX","Uncertainty x Coordinate (tagged Vtx)",100,0.,100);
   yerr[1]     = dqmStore_->book1D("tagErrY","Uncertainty y Coordinate (tagged Vtx)",100,0.,100);
   zerr[1]     = dqmStore_->book1D("tagErrZ","Uncertainty z Coordinate (tagged Vtx)",100,0.,100);
   xerrVsTrks[1]  = dqmStore_->book2D("tagErrVsWeightX","Uncertainty x Coordinate vs. track weight (tagged Vtx)",100,0,100.,100,0.,100);
   yerrVsTrks[1]  = dqmStore_->book2D("tagErrVsWeightY","Uncertainty y Coordinate vs. track weight (tagged Vtx)",100,0,100.,100,0.,100);
   zerrVsTrks[1]  = dqmStore_->book2D("tagErrVsWeightZ","Uncertainty z Coordinate vs. track weight (tagged Vtx)",100,0,100.,100,0.,100);
 
   type[0] = dqmStore_->book1D("otherType","Vertex type (other Vtx)",3,-0.5,2.5);
   type[1] = dqmStore_->book1D("tagType","Vertex type (tagged Vtx)",3,-0.5,2.5);
   for (int i=0;i<2;++i){
     type[i]->getTH1F()->GetXaxis()->SetBinLabel(1,"Valid, real");
     type[i]->getTH1F()->GetXaxis()->SetBinLabel(2,"Valid, fake");
     type[i]->getTH1F()->GetXaxis()->SetBinLabel(3,"Invalid");
   }
 
   bsX       = dqmStore_->book1D("bsX", "BeamSpot x0", 100,-0.1,0.1);
   bsY       = dqmStore_->book1D("bsY", "BeamSpot y0", 100,-0.1,0.1);
   bsZ       = dqmStore_->book1D("bsZ", "BeamSpot z0", 100,-2.,2.);
   bsSigmaZ  = dqmStore_->book1D("bsSigmaZ", "BeamSpot sigmaZ", 100, 0., 10. );
   bsDxdz    = dqmStore_->book1D("bsDxdz", "BeamSpot dxdz", 100, -0.0003, 0.0003);
   bsDydz    = dqmStore_->book1D("bsDydz", "BeamSpot dydz", 100, -0.0003, 0.0003);
   bsBeamWidthX  = dqmStore_->book1D("bsBeamWidthX", "BeamSpot BeamWidthX", 100, 0., 100.);
   bsBeamWidthY  = dqmStore_->book1D("bsBeamWidthY", "BeamSpot BeamWidthY", 100, 0., 100.);
   bsType    = dqmStore_->book1D("bsType", "BeamSpot type", 4, -1.5, 2.5);
   bsType->getTH1F()->GetXaxis()->SetBinLabel(1,"Unknown");
   bsType->getTH1F()->GetXaxis()->SetBinLabel(2,"Fake");
   bsType->getTH1F()->GetXaxis()->SetBinLabel(3,"LHC");
   bsType->getTH1F()->GetXaxis()->SetBinLabel(4,"Tracker");
 }

PrimaryVertexMonitor::~PrimaryVertexMonitor ( )

Definition at line 95 of file PrimaryVertexMonitor.cc.

96 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 98 of file PrimaryVertexMonitor.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), edm::Event::getByLabel(), edm::HandleBase::isValid(), reco::BeamSpot::sigmaZ(), findQualityFiles::v, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

 {
   Handle<reco::VertexCollection> recVtxs;
   iEvent.getByLabel(moduleLabel, recVtxs);
 
   edm::Handle<reco::BeamSpot> beamSpotHandle;
   iEvent.getByLabel(beamSpotLabel,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 "
       <<moduleLabel<<" " 
       <<recVtxs.isValid() <<" BeamSpot "
       <<beamSpotLabel<<" "
       <<beamSpotHandle.isValid()<<". Skipping plots for this event";
     return;
   }
 
   BeamSpot beamSpot = *beamSpotHandle;
 
   nbvtx->Fill(recVtxs->size()*1.);
 
   vertexPlots(recVtxs->front(), beamSpot, 1);
 
   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()*10000);
   bsBeamWidthY->Fill(beamSpot.BeamWidthY()*10000);
   // bsType->Fill(beamSpot.type());
 
 }

void PrimaryVertexMonitor::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 196 of file PrimaryVertexMonitor.cc.

197 {

198 }

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

private

Definition at line 142 of file PrimaryVertexMonitor.cc.

References reco::Vertex::chi2(), ChiSquaredProbability(), reco::Vertex::covariance(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), HcalObjRepresent::Fill(), i, getHLTprescales::index, reco::Vertex::isFake(), edm::isNotFinite(), reco::Vertex::isValid(), j, gen::k, reco::Vertex::ndof(), reco::Vertex::position(), lumiQTWidget::t, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), reco::Vertex::trackWeight(), histoStyle::weight, reco::BeamSpot::x0(), reco::Vertex::xError(), reco::BeamSpot::y0(), reco::Vertex::yError(), reco::BeamSpot::z0(), spu::zerr(), and reco::Vertex::zError().

 {
 
     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());
 
       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)*10000);
       yDiff[i]->Fill((v.position().y() - yb)*10000);
 
       xerr[i]->Fill(v.xError()*10000);
       yerr[i]->Fill(v.yError()*10000);
       zerr[i]->Fill(v.zError()*10000);
       xerrVsTrks[i]->Fill(weight, v.xError()*10000);
       yerrVsTrks[i]->Fill(weight, v.yError()*10000);
       zerrVsTrks[i]->Fill(weight, v.zError()*10000);
 
       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

edm::InputTag PrimaryVertexMonitor::beamSpotLabel

private

Definition at line 37 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsBeamWidthX

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsBeamWidthY

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsDxdz

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsDydz

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsSigmaZ

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsType

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement* PrimaryVertexMonitor::bsX

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsY

private

Definition at line 49 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::bsZ

private

Definition at line 49 of file PrimaryVertexMonitor.h.

std::string PrimaryVertexMonitor::dqmLabel

private

Definition at line 40 of file PrimaryVertexMonitor.h.

DQMStore* PrimaryVertexMonitor::dqmStore_

private

Definition at line 39 of file PrimaryVertexMonitor.h.

edm::InputTag PrimaryVertexMonitor::moduleLabel

private

Definition at line 37 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::nans[2]

private

Definition at line 47 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::nbtksinvtx[2]

private

Definition at line 43 of file PrimaryVertexMonitor.h.

MonitorElement* PrimaryVertexMonitor::nbvtx

private

Definition at line 43 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::trksWeight[2]

private

Definition at line 43 of file PrimaryVertexMonitor.h.

MonitorElement* PrimaryVertexMonitor::tt[2]

private

Definition at line 44 of file PrimaryVertexMonitor.h.

MonitorElement* PrimaryVertexMonitor::type[2]

private

Definition at line 48 of file PrimaryVertexMonitor.h.

Referenced by cuy.ValElement::__init__(), argparse.Action::__init__(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties().

MonitorElement* PrimaryVertexMonitor::vtxchi2[2]

private

Definition at line 47 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::vtxndf[2]

private

Definition at line 47 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::vtxprob[2]

private

Definition at line 47 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::xDiff[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::xerr[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

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

MonitorElement* PrimaryVertexMonitor::xerrVsTrks[2]

private

Definition at line 46 of file PrimaryVertexMonitor.h.

MonitorElement* PrimaryVertexMonitor::xrec[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::yDiff[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::yerr[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

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

MonitorElement * PrimaryVertexMonitor::yerrVsTrks[2]

private

Definition at line 46 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::yrec[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::zerr[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::zerrVsTrks[2]

private

Definition at line 46 of file PrimaryVertexMonitor.h.

MonitorElement * PrimaryVertexMonitor::zrec[2]

private

Definition at line 45 of file PrimaryVertexMonitor.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation