PrimaryVertexMonitor Class Reference

#include <PrimaryVertexMonitor.h>

Inheritance diagram for PrimaryVertexMonitor:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	PrimaryVertexMonitor (const edm::ParameterSet &pSet)
	~PrimaryVertexMonitor ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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

Private Attributes
std::string	AlignmentLabel_
edm::InputTag	beamSpotInputTag_
edm::EDGetTokenT< reco::BeamSpot >	beamspotToken_
MonitorElement *	bsBeamWidthX
MonitorElement *	bsBeamWidthY
MonitorElement *	bsDxdz
MonitorElement *	bsDydz
MonitorElement *	bsSigmaZ
MonitorElement *	bsType
MonitorElement *	bsX
MonitorElement *	bsY
MonitorElement *	bsZ
MonitorElement *	chi2ndf
MonitorElement *	chi2prob
edm::ParameterSet	conf_
std::string	dqmLabel
MonitorElement *	dxy
MonitorElement *	dxy2
MonitorElement *	dxyErr
MonitorElement *	dxyVsEta_pt1
MonitorElement *	dxyVsEta_pt10
MonitorElement *	dxyVsPhi_pt1
MonitorElement *	dxyVsPhi_pt10
MonitorElement *	dz
MonitorElement *	dzErr
MonitorElement *	dzVsEta_pt1
MonitorElement *	dzVsEta_pt10
MonitorElement *	dzVsPhi_pt1
MonitorElement *	dzVsPhi_pt10
MonitorElement *	nans [2]
MonitorElement *	nbgvtx
MonitorElement *	nbtksinvtx [2]
MonitorElement *	nbvtx
MonitorElement *	ntracks
MonitorElement *	sumpt
std::string	TopFolderName_
MonitorElement *	trksWeight [2]
MonitorElement *	tt [2]
MonitorElement *	type [2]
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]

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 25 of file PrimaryVertexMonitor.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file PrimaryVertexMonitor.cc.

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

  : conf_          ( pSet )
  , TopFolderName_ ( pSet.getParameter<std::string>("TopFolderName") )
  , AlignmentLabel_( pSet.getParameter<std::string>("AlignmentLabel"))
  , nbvtx(NULL)
  , bsX(NULL)
  , bsY(NULL)
  , bsZ(NULL)
  , bsSigmaZ(NULL)
  , bsDxdz(NULL)
  , bsDydz(NULL)
  , bsBeamWidthX(NULL)
  , bsBeamWidthY(NULL)
  , bsType(NULL)
  , sumpt(NULL)
  , ntracks(NULL)
  , weight(NULL)
  , chi2ndf(NULL)
  , chi2prob(NULL)
  , dxy(NULL)
  , dxy2(NULL)
  , dz(NULL)
  , dxyErr(NULL)
  , dzErr(NULL)
  , dxyVsPhi_pt1(NULL)
  , dzVsPhi_pt1(NULL)
  , dxyVsEta_pt1(NULL)
  , dzVsEta_pt1(NULL)
  , dxyVsPhi_pt10(NULL)
  , dzVsPhi_pt10(NULL)
  , dxyVsEta_pt10(NULL)
  , dzVsEta_pt10(NULL)
{
  //  dqmStore_ = edm::Service<DQMStore>().operator->();


  vertexInputTag_   = pSet.getParameter<InputTag>("vertexLabel");
  beamSpotInputTag_ = pSet.getParameter<InputTag>("beamSpotLabel");
  vertexToken_   = consumes<reco::VertexCollection>(vertexInputTag_);
  beamspotToken_ = consumes<reco::BeamSpot>        (beamSpotInputTag_);

}

PrimaryVertexMonitor::~PrimaryVertexMonitor

(

)

Definition at line 224 of file PrimaryVertexMonitor.cc.

{}

Member Function Documentation

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

overridevirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 227 of file PrimaryVertexMonitor.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, beamSpotInputTag_, beamspotToken_, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), bsBeamWidthX, bsBeamWidthY, bsDxdz, bsDydz, bsSigmaZ, bsX, bsY, bsZ, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), MonitorElement::Fill(), edm::Event::getByToken(), edm::HandleBase::isValid(), nbgvtx, nbvtx, pvTracksPlots(), reco::BeamSpot::sigmaZ(), findQualityFiles::v, vertexInputTag_, vertexPlots(), vertexToken_, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

{
  Handle<reco::VertexCollection> recVtxs;
  iEvent.getByToken(vertexToken_, recVtxs);

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

  BeamSpot beamSpot = *beamSpotHandle;

  nbvtx->Fill(recVtxs->size()*1.);
  int ng=0;
  for (auto const & vx : (*recVtxs) )
    if (vx.isValid() && !vx.isFake()  && vx.ndof()>=4.) ++ng;
  nbgvtx->Fill(ng*1.);

  vertexPlots(recVtxs->front(), beamSpot, 1);

  // fill PV tracks MEs (as now, for alignment)
  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()*10000);
  bsBeamWidthY->Fill(beamSpot.BeamWidthY()*10000);
  // bsType->Fill(beamSpot.type());

}

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 62 of file PrimaryVertexMonitor.cc.

References AlignmentLabel_, beamSpotInputTag_, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), bsBeamWidthX, bsBeamWidthY, bsDxdz, bsDydz, bsSigmaZ, bsType, bsX, bsY, bsZ, chi2ndf, chi2prob, conf_, dqmLabel, dxy, dxy2, dxyErr, dxyVsEta_pt1, dxyVsEta_pt10, dxyVsPhi_pt1, dxyVsPhi_pt10, dz, dzErr, dzVsEta_pt1, dzVsEta_pt10, dzVsPhi_pt1, dzVsPhi_pt10, HLT_25ns14e33_v3_cff::EtaMax, Gflash::EtaMin, edm::ParameterSet::getParameter(), MonitorElement::getTH1F(), i, edm::InputTag::label(), nans, nbgvtx, nbtksinvtx, nbvtx, ntracks, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, sumpt, TopFolderName_, trksWeight, vertexInputTag_, vtxchi2, vtxndf, vtxprob, weight, 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",50,-0.5,49.5);
  nbgvtx      = iBooker.book1D("goodvtxNbr","Reconstructed Good Vertices in Event",50,-0.5,49.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); 
  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); 
  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,vposx-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]->getTH1F()->GetXaxis()->SetBinLabel(1,"Valid, real");
    type[i]->getTH1F()->GetXaxis()->SetBinLabel(2,"Valid, fake");
    type[i]->getTH1F()->GetXaxis()->SetBinLabel(3,"Invalid");
  }


  //  get the store
  dqmLabel = TopFolderName_+"/"+beamSpotInputTag_.label();
  iBooker.setCurrentFolder(dqmLabel);
  
  bsX       = iBooker.book1D("bsX", "BeamSpot x0", 100,-0.1,0.1);
  bsY       = iBooker.book1D("bsY", "BeamSpot y0", 100,-0.1,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->getTH1F()->GetXaxis()->SetBinLabel(1,"Unknown");
  bsType->getTH1F()->GetXaxis()->SetBinLabel(2,"Fake");
  bsType->getTH1F()->GetXaxis()->SetBinLabel(3,"LHC");
  bsType->getTH1F()->GetXaxis()->SetBinLabel(4,"Tracker");

  
  //  get the store
  dqmLabel = TopFolderName_+"/"+AlignmentLabel_;
  iBooker.setCurrentFolder(dqmLabel);

  int    TKNoBin    = conf_.getParameter<int>(   "TkSizeBin");
  double TKNoMin    = conf_.getParameter<double>("TkSizeMin");
  double TKNoMax    = conf_.getParameter<double>("TkSizeMax");

  int    DxyBin     = conf_.getParameter<int>(   "DxyBin");
  double DxyMin     = conf_.getParameter<double>("DxyMin");
  double DxyMax     = conf_.getParameter<double>("DxyMax");
  
  int    DzBin      = conf_.getParameter<int>(   "DzBin");
  double DzMin      = conf_.getParameter<double>("DzMin");
  double DzMax      = conf_.getParameter<double>("DzMax");
  
  int    PhiBin     = conf_.getParameter<int>(   "PhiBin");
  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. );
  dxy      = iBooker.book1D("dxy",     "PV tracks (p_{T} > 1 GeV) d_{xy} (cm)",         DxyBin, DxyMin, DxyMax);
  dxy2     = iBooker.book1D("dxyzoom", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm)",         DxyBin, DxyMin/5., DxyMax/5.);
  dz       = iBooker.book1D("dz",      "PV tracks (p_{T} > 1 GeV) d_{z} (cm)",          DzBin,  DzMin,  DzMax );
  dxyErr   = iBooker.book1D("dxyErr",  "PV tracks (p_{T} > 1 GeV) d_{xy} error (cm)",   100, 0.,   0.2 );
  dzErr    = iBooker.book1D("dzErr",   "PV tracks (p_{T} > 1 GeV) d_{z} error(cm)",     100, 0.,   1. );

  dxyVsPhi_pt1 = iBooker.bookProfile("dxyVsPhi_pt1", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm) VS track #phi",PhiBin, PhiMin, PhiMax, DxyBin, DxyMin, DxyMax,"");
  dxyVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi",  1);
  dxyVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy}",2);
  dzVsPhi_pt1  = iBooker.bookProfile("dzVsPhi_pt1",  "PV tracks (p_{T} > 1 GeV) d_{z} (cm) VS track #phi", PhiBin, PhiMin, PhiMax, DzBin,  DzMin,  DzMax, "");  
  dzVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #phi", 1);
  dzVsPhi_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z}",2);

  dxyVsEta_pt1 = iBooker.bookProfile("dxyVsEta_pt1", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm) VS track #eta",EtaBin, EtaMin, EtaMax, DxyBin, DxyMin, DxyMax,"");
  dxyVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta",  1);
  dxyVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{xy}",2);
  dzVsEta_pt1  = iBooker.bookProfile("dzVsEta_pt1",  "PV tracks (p_{T} > 1 GeV) d_{z} (cm) VS track #eta", EtaBin, EtaMin, EtaMax, DzBin,  DzMin,  DzMax, "");
  dzVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) #eta", 1);
  dzVsEta_pt1->setAxisTitle("PV track (p_{T} > 1 GeV) d_{z}",2);

  dxyVsPhi_pt10 = iBooker.bookProfile("dxyVsPhi_pt10", "PV tracks (p_{T} > 1 GeV) d_{xy} (cm) VS track #phi",PhiBin, PhiMin, PhiMax, DxyBin, DxyMin, DxyMax,"");
  dxyVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi",  1);
  dxyVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy}",2);
  dzVsPhi_pt10  = iBooker.bookProfile("dzVsPhi_pt10",  "PV tracks (p_{T} > 10 GeV) d_{z} (cm) VS track #phi", PhiBin, PhiMin, PhiMax, DzBin,  DzMin,  DzMax, "");  
  dzVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #phi", 1);
  dzVsPhi_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z}",2);

  dxyVsEta_pt10 = iBooker.bookProfile("dxyVsEta_pt10", "PV tracks (p_{T} > 10 GeV) d_{xy} (cm) VS track #eta",EtaBin, EtaMin, EtaMax, DxyBin, DxyMin, DxyMax,"");
  dxyVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta",  1);
  dxyVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{xy}",2);
  dzVsEta_pt10  = iBooker.bookProfile("dzVsEta_pt10",  "PV tracks (p_{T} > 10 GeV) d_{z} (cm) VS track #eta", EtaBin, EtaMin, EtaMax, DzBin,  DzMin,  DzMax, "");
  dzVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) #eta", 1);
  dzVsEta_pt10->setAxisTitle("PV track (p_{T} > 10 GeV) d_{z}",2);

}

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

private

Definition at line 279 of file PrimaryVertexMonitor.cc.

References chi2ndf, chi2prob, dxy, dxy2, dxyErr, dxyVsEta_pt1, dxyVsEta_pt10, dxyVsPhi_pt1, dxyVsPhi_pt10, dz, dzErr, dzVsEta_pt1, dzVsEta_pt10, dzVsPhi_pt1, dzVsPhi_pt10, eta, HcalObjRepresent::Fill(), reco::TrackBase::highPurity, reco::Vertex::isFake(), reco::Vertex::isValid(), ntracks, phi, reco::Vertex::position(), EnergyCorrector::pt, sumpt, tree::t, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), reco::Vertex::trackWeight(), w, and weight.

Referenced by analyze().

{

  const math::XYZPoint myVertex(v.position().x(),v.position().y(),v.position().z());

  if ( !v.isValid() ) return;
  if (  v.isFake()  ) return;

  if ( v.tracksSize() == 0 ) {
    ntracks -> Fill ( 0 );
    return;
  }

  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 ) 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);
    float Dz       = (**t).dz(myVertex);
    float DxyErr   = (**t).dxyError();
    float DzErr    = (**t).dzError();

    sumPT += pt*pt;

    
    // fill MEs
    weight   -> Fill (w);
    chi2ndf  -> Fill (chi2NDF);
    chi2prob -> Fill (chi2Prob);
    dxy      -> Fill (Dxy);
    dxy2     -> Fill (Dxy); 
    dz       -> Fill (Dz);
    dxyErr   -> Fill (DxyErr);
    dzErr    -> Fill (DzErr);
    
    dxyVsPhi_pt1 -> Fill (phi,Dxy);
    dzVsPhi_pt1  -> Fill (phi,Dz);
    dxyVsEta_pt1 -> Fill (eta,Dxy);
    dzVsEta_pt1  -> Fill (eta,Dz);

    if ( pt < 10. ) continue;
    dxyVsPhi_pt10 -> Fill (phi,Dxy);
    dzVsPhi_pt10  -> Fill (phi,Dz);
    dxyVsEta_pt10 -> Fill (eta,Dxy);
    dzVsEta_pt10  -> Fill (eta,Dz);
  }
  ntracks -> Fill (float(nTracks));
  sumpt   -> Fill (sumPT);

}

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

private

Definition at line 344 of file PrimaryVertexMonitor.cc.

References reco::Vertex::chi2(), ChiSquaredProbability(), reco::Vertex::covariance(), reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), MonitorElement::Fill(), HcalObjRepresent::Fill(), i, cmsHarvester::index, reco::Vertex::isFake(), edm::isNotFinite(), reco::Vertex::isValid(), j, relval_steps::k, nans, nbtksinvtx, reco::Vertex::ndof(), reco::Vertex::position(), tree::t, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), reco::Vertex::tracksSize(), reco::Vertex::trackWeight(), trksWeight, vtxchi2, vtxndf, vtxprob, weight, reco::BeamSpot::x0(), xDiff, xerr, reco::Vertex::xError(), xerrVsTrks, xrec, reco::BeamSpot::y0(), yDiff, yerr, reco::Vertex::yError(), yerrVsTrks, yrec, reco::BeamSpot::z0(), zerr, reco::Vertex::zError(), 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());

      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

std::string PrimaryVertexMonitor::AlignmentLabel_

private

Definition at line 49 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

edm::InputTag PrimaryVertexMonitor::beamSpotInputTag_

private

Definition at line 42 of file PrimaryVertexMonitor.h.

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

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

private

Definition at line 40 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and PrimaryVertexMonitor().

MonitorElement * PrimaryVertexMonitor::bsBeamWidthX

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsBeamWidthY

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsDxdz

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsDydz

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsSigmaZ

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsType

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

MonitorElement* PrimaryVertexMonitor::bsX

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsY

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::bsZ

private

Definition at line 58 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::chi2ndf

private

Definition at line 60 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::chi2prob

private

Definition at line 60 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

edm::ParameterSet PrimaryVertexMonitor::conf_

private

Definition at line 44 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

std::string PrimaryVertexMonitor::dqmLabel

private

Definition at line 46 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

MonitorElement* PrimaryVertexMonitor::dxy

private

Definition at line 61 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dxy2

private

Definition at line 61 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dxyErr

private

Definition at line 61 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement* PrimaryVertexMonitor::dxyVsEta_pt1

private

Definition at line 63 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement* PrimaryVertexMonitor::dxyVsEta_pt10

private

Definition at line 65 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement* PrimaryVertexMonitor::dxyVsPhi_pt1

private

Definition at line 62 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement* PrimaryVertexMonitor::dxyVsPhi_pt10

private

Definition at line 64 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dz

private

Definition at line 61 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dzErr

private

Definition at line 61 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dzVsEta_pt1

private

Definition at line 63 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dzVsEta_pt10

private

Definition at line 65 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dzVsPhi_pt1

private

Definition at line 62 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::dzVsPhi_pt10

private

Definition at line 64 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement * PrimaryVertexMonitor::nans[2]

private

Definition at line 56 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::nbgvtx

private

Definition at line 52 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::nbtksinvtx[2]

private

Definition at line 52 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement* PrimaryVertexMonitor::nbvtx

private

Definition at line 52 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement * PrimaryVertexMonitor::ntracks

private

Definition at line 60 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and pvTracksPlots().

MonitorElement* PrimaryVertexMonitor::sumpt

private

Definition at line 60 of file PrimaryVertexMonitor.h.

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

std::string PrimaryVertexMonitor::TopFolderName_

private

Definition at line 48 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms().

MonitorElement * PrimaryVertexMonitor::trksWeight[2]

private

Definition at line 52 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement* PrimaryVertexMonitor::tt[2]

private

Definition at line 53 of file PrimaryVertexMonitor.h.

MonitorElement* PrimaryVertexMonitor::type[2]

private

Definition at line 57 of file PrimaryVertexMonitor.h.

Referenced by cuy.ValElement::__init__(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), core.autovars.NTupleVariable::makeBranch(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and core.AutoHandle.AutoHandle::ReallyLoad().

edm::InputTag PrimaryVertexMonitor::vertexInputTag_

private

Definition at line 42 of file PrimaryVertexMonitor.h.

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

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

private

Definition at line 39 of file PrimaryVertexMonitor.h.

Referenced by analyze(), and PrimaryVertexMonitor().

MonitorElement* PrimaryVertexMonitor::vtxchi2[2]

private

Definition at line 56 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::vtxndf[2]

private

Definition at line 56 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::vtxprob[2]

private

Definition at line 56 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::weight

private

Definition at line 60 of file PrimaryVertexMonitor.h.

Referenced by cuy.ValElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), bookHistograms(), dirstructure.Directory::calcStats(), dirstructure.Directory::get_fail_rate(), dirstructure.Directory::get_null_rate(), dirstructure.Directory::get_skiped_rate(), dirstructure.Directory::get_success_rate(), dirstructure.Directory::print_report(), pvTracksPlots(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::xDiff[2]

private

Definition at line 54 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::xerr[2]

private

Definition at line 54 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().

MonitorElement* PrimaryVertexMonitor::xerrVsTrks[2]

private

Definition at line 55 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement* PrimaryVertexMonitor::xrec[2]

private

Definition at line 54 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::yDiff[2]

private

Definition at line 54 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::yerr[2]

private

Definition at line 54 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().

MonitorElement * PrimaryVertexMonitor::yerrVsTrks[2]

private

Definition at line 55 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::yrec[2]

private

Definition at line 54 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::zerr[2]

private

Definition at line 54 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::zerrVsTrks[2]

private

Definition at line 55 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().

MonitorElement * PrimaryVertexMonitor::zrec[2]

private

Definition at line 54 of file PrimaryVertexMonitor.h.

Referenced by bookHistograms(), and vertexPlots().