VertexMonitor Class Reference

#include <VertexMonitor.h>

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	initHisto (DQMStore::IBooker &ibooker)
	VertexMonitor (const edm::ParameterSet &, const edm::InputTag &, const edm::InputTag &, std::string pvLabel)
	VertexMonitor (const edm::ParameterSet &, const edm::InputTag &, const edm::InputTag &, std::string pvLabel, edm::ConsumesCollector &iC)
virtual	~VertexMonitor ()

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

Public Attributes
MonitorElement *	Chi2oNDFVsBXlumi
MonitorElement *	Chi2oNDFVsGoodPVtx
MonitorElement *	Chi2ProbVsBXlumi
MonitorElement *	Chi2ProbVsGoodPVtx
edm::ParameterSet	conf_
bool	doAllPlots_
bool	doPlotsVsBXlumi_
bool	doPlotsVsGoodPVtx_
MonitorElement *	FractionOfGoodPVtx
MonitorElement *	FractionOfGoodPVtxVsBXlumi
MonitorElement *	FractionOfGoodPVtxVsGoodPVtx
MonitorElement *	FractionOfGoodPVtxVsPVtx
MonitorElement *	GoodPVtxChi2oNDFVsBXlumi
MonitorElement *	GoodPVtxChi2oNDFVsGoodPVtx
MonitorElement *	GoodPVtxChi2ProbVsBXlumi
MonitorElement *	GoodPVtxChi2ProbVsGoodPVtx
MonitorElement *	GoodPVtxNumberOfTracks
MonitorElement *	GoodPVtxNumberOfTracksVsBXlumi
MonitorElement *	GoodPVtxNumberOfTracksVsGoodPVtx
MonitorElement *	GoodPVtxNumberOfTracksVsGoodPVtxNdof
MonitorElement *	GoodPVtxSumPt
MonitorElement *	GoodPVtxSumPtVsBXlumi
MonitorElement *	GoodPVtxSumPtVsGoodPVtx
std::string	histname
std::string	label_
GetLumi *	lumiDetails_
MonitorElement *	NumberOfBADndofPVtx
MonitorElement *	NumberOfBADndofPVtxVsBXlumi
MonitorElement *	NumberOfBADndofPVtxVsGoodPVtx
MonitorElement *	NumberOfFakePVtx
MonitorElement *	NumberOfFakePVtxVsBXlumi
MonitorElement *	NumberOfFakePVtxVsGoodPVtx
MonitorElement *	NumberOfGoodPVtx
MonitorElement *	NumberOfGoodPVtxVsBXlumi
MonitorElement *	NumberOfPVtx
MonitorElement *	NumberOfPVtxVsBXlumi
MonitorElement *	NumberOfPVtxVsGoodPVtx
edm::InputTag	primaryVertexInputTag_
edm::EDGetTokenT < reco::VertexCollection >	pvToken_
edm::InputTag	selectedPrimaryVertexInputTag_
edm::EDGetTokenT < reco::VertexCollection >	selpvToken_

Detailed Description

Monitoring source for general quantities related to vertex

Definition at line 32 of file VertexMonitor.h.

Constructor & Destructor Documentation

VertexMonitor::VertexMonitor	(	const edm::ParameterSet &	iConfig,
		const edm::InputTag &	primaryVertexInputTag,
		const edm::InputTag &	selectedPrimaryVertexInputTag,
		std::string	pvLabel
	)

Definition at line 19 of file VertexMonitor.cc.

References doPlotsVsBXlumi_, edm::ParameterSet::getParameter(), and lumiDetails_.

    : conf_( iConfig )
    , primaryVertexInputTag_         ( primaryVertexInputTag )
    , selectedPrimaryVertexInputTag_ ( selectedPrimaryVertexInputTag )
    , label_                         ( pvLabel )
    , NumberOfPVtx(NULL)
    , NumberOfPVtxVsBXlumi(NULL)
    , NumberOfPVtxVsGoodPVtx(NULL)
    , NumberOfGoodPVtx(NULL)
    , NumberOfGoodPVtxVsBXlumi(NULL)
    , FractionOfGoodPVtx(NULL)
    , FractionOfGoodPVtxVsBXlumi(NULL)
    , FractionOfGoodPVtxVsGoodPVtx(NULL)
    , FractionOfGoodPVtxVsPVtx(NULL)
    , NumberOfBADndofPVtx(NULL)
    , NumberOfBADndofPVtxVsBXlumi(NULL)
    , NumberOfBADndofPVtxVsGoodPVtx(NULL)
    , GoodPVtxSumPt(NULL)
    , GoodPVtxSumPtVsBXlumi(NULL)
    , GoodPVtxSumPtVsGoodPVtx(NULL)
    , GoodPVtxNumberOfTracks(NULL)
    , GoodPVtxNumberOfTracksVsBXlumi(NULL)
    , GoodPVtxNumberOfTracksVsGoodPVtx(NULL)
    , GoodPVtxNumberOfTracksVsGoodPVtxNdof(NULL)
    , GoodPVtxChi2oNDFVsGoodPVtx(NULL)
    , GoodPVtxChi2oNDFVsBXlumi(NULL)
    , GoodPVtxChi2ProbVsGoodPVtx(NULL)
    , GoodPVtxChi2ProbVsBXlumi(NULL)
    , doAllPlots_       ( conf_.getParameter<bool>("doAllPlots")        )
    , doPlotsVsBXlumi_  ( conf_.getParameter<bool>("doPlotsVsBXlumi")   )
    , doPlotsVsGoodPVtx_( conf_.getParameter<bool>("doPlotsVsGoodPVtx") )

{
   //now do what ever initialization is needed
  if ( doPlotsVsBXlumi_ )
    lumiDetails_ = new GetLumi( iConfig.getParameter<edm::ParameterSet>("BXlumiSetup") );

}

VertexMonitor::VertexMonitor	(	const edm::ParameterSet &	iConfig,
		const edm::InputTag &	primaryVertexInputTag,
		const edm::InputTag &	selectedPrimaryVertexInputTag,
		std::string	pvLabel,
		edm::ConsumesCollector &	iC
	)

Definition at line 58 of file VertexMonitor.cc.

References edm::ConsumesCollector::consumes(), doPlotsVsBXlumi_, edm::ParameterSet::getParameter(), lumiDetails_, primaryVertexInputTag_, pvToken_, selectedPrimaryVertexInputTag_, and selpvToken_.

                                                                                                                                                                                                    : VertexMonitor(iConfig,primaryVertexInputTag,selectedPrimaryVertexInputTag,pvLabel)
{

  if ( doPlotsVsBXlumi_ )
    lumiDetails_ = new GetLumi( iConfig.getParameter<edm::ParameterSet>("BXlumiSetup"), iC );

  pvToken_    = iC.consumes<reco::VertexCollection>(primaryVertexInputTag_);
  selpvToken_ = iC.consumes<reco::VertexCollection>(selectedPrimaryVertexInputTag_); 
}

VertexMonitor::~VertexMonitor ( )

virtual

Definition at line 68 of file VertexMonitor.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
  //  if (lumiDetails_) delete lumiDetails_;
}

Member Function Documentation

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

virtual

Definition at line 85 of file VertexMonitor.cc.

References reco::Vertex::chi2(), gather_cfg::cout, doPlotsVsBXlumi_, doPlotsVsGoodPVtx_, HcalObjRepresent::Fill(), FractionOfGoodPVtx, FractionOfGoodPVtxVsBXlumi, FractionOfGoodPVtxVsGoodPVtx, FractionOfGoodPVtxVsPVtx, edm::Event::getByToken(), GetLumi::getValue(), GoodPVtxChi2oNDFVsBXlumi, GoodPVtxChi2oNDFVsGoodPVtx, GoodPVtxChi2ProbVsBXlumi, GoodPVtxChi2ProbVsGoodPVtx, GoodPVtxNumberOfTracks, GoodPVtxNumberOfTracksVsBXlumi, GoodPVtxNumberOfTracksVsGoodPVtx, GoodPVtxSumPt, GoodPVtxSumPtVsBXlumi, GoodPVtxSumPtVsGoodPVtx, edm::HandleBase::isValid(), lumiDetails_, reco::Vertex::ndof(), reco::Vertex::normalizedChi2(), NumberOfBADndofPVtx, NumberOfBADndofPVtxVsBXlumi, NumberOfBADndofPVtxVsGoodPVtx, NumberOfGoodPVtx, NumberOfGoodPVtxVsBXlumi, NumberOfPVtx, NumberOfPVtxVsBXlumi, NumberOfPVtxVsGoodPVtx, EnergyCorrector::pt, MetAnalyzer::pv(), pvToken_, selpvToken_, reco::Vertex::tracks_begin(), reco::Vertex::tracks_end(), and reco::Vertex::tracksSize().

{

  double bxlumi = 0.;
  if ( doPlotsVsBXlumi_ )
    bxlumi = lumiDetails_->getValue(iEvent);
  std::cout << "bxlumi : " << bxlumi << std::endl;

  size_t totalNumPV = 0;
  size_t totalNumBADndofPV = 0;
  edm::Handle< reco::VertexCollection > pvHandle;
  iEvent.getByToken(pvToken_, pvHandle );
  if ( pvHandle.isValid() )
    {
      totalNumPV = pvHandle->size();
      std::cout << "totalNumPV : " << totalNumPV << std::endl;
      for (reco::VertexCollection::const_iterator pv = pvHandle->begin();
       pv != pvHandle->end(); ++pv) {
    //--- count pv w/ ndof < 4 
    if (pv->ndof() <  4.) totalNumBADndofPV++;
      }
    } else return;
  NumberOfPVtx        -> Fill( totalNumPV        );
  NumberOfBADndofPVtx -> Fill( totalNumBADndofPV );
  if ( doPlotsVsBXlumi_ ) {
    NumberOfPVtxVsBXlumi        -> Fill( bxlumi, totalNumPV        );
    NumberOfBADndofPVtxVsBXlumi -> Fill( bxlumi, totalNumBADndofPV );
  }
  
  size_t totalNumGoodPV = 0;
  edm::Handle< reco::VertexCollection > selpvHandle;
  iEvent.getByToken(selpvToken_, selpvHandle );
  if ( selpvHandle.isValid() )
    totalNumGoodPV = selpvHandle->size();
  else return;
  std::cout << "totalNumGoodPV: " << totalNumGoodPV << std::endl;
  if ( doPlotsVsGoodPVtx_ ) {
    NumberOfPVtxVsGoodPVtx        -> Fill( totalNumGoodPV, totalNumPV        );
    NumberOfBADndofPVtxVsGoodPVtx -> Fill( totalNumGoodPV, totalNumBADndofPV );
  }

  double fracGoodPV = double(totalNumGoodPV)/double(totalNumPV);
  std::cout << "fracGoodPV: " << fracGoodPV << std::endl;

  NumberOfGoodPVtx    -> Fill( totalNumGoodPV    );
  FractionOfGoodPVtx  -> Fill( fracGoodPV        );
  if ( doPlotsVsBXlumi_ ) {
    NumberOfGoodPVtxVsBXlumi    -> Fill( bxlumi, totalNumGoodPV    );
    FractionOfGoodPVtxVsBXlumi  -> Fill( bxlumi, fracGoodPV        );
  }
  if ( doPlotsVsGoodPVtx_ ) {
    FractionOfGoodPVtxVsGoodPVtx  -> Fill( totalNumGoodPV, fracGoodPV        );
    FractionOfGoodPVtxVsPVtx      -> Fill( totalNumPV,     fracGoodPV        );
  }

  if ( selpvHandle->size() ) {
    double sumpt    = 0;
    size_t ntracks  = 0;
    double chi2ndf  = 0.; 
    double chi2prob = 0.;

    if (!selpvHandle->at(0).isFake()) {
      
      reco::Vertex pv = selpvHandle->at(0);
      
      ntracks  = pv.tracksSize();
      chi2ndf  = pv.normalizedChi2();
      chi2prob = TMath::Prob(pv.chi2(),(int)pv.ndof());
      
      for (reco::Vertex::trackRef_iterator itrk = pv.tracks_begin();
       itrk != pv.tracks_end(); ++itrk) {
    double pt = (**itrk).pt();
    sumpt += pt*pt;
      }
      GoodPVtxSumPt           -> Fill( sumpt   );
      GoodPVtxNumberOfTracks  -> Fill( ntracks );

      if ( doPlotsVsBXlumi_ ) {
    GoodPVtxSumPtVsBXlumi           -> Fill( bxlumi, sumpt    );
    GoodPVtxNumberOfTracksVsBXlumi  -> Fill( bxlumi, ntracks  );
    GoodPVtxChi2oNDFVsBXlumi        -> Fill( bxlumi, chi2ndf  );
    GoodPVtxChi2ProbVsBXlumi        -> Fill( bxlumi, chi2prob );
      }
      if ( doPlotsVsGoodPVtx_ ) {
    GoodPVtxSumPtVsGoodPVtx          -> Fill( totalNumGoodPV, sumpt    );
    GoodPVtxNumberOfTracksVsGoodPVtx -> Fill( totalNumGoodPV, ntracks  );
    GoodPVtxChi2oNDFVsGoodPVtx       -> Fill( totalNumGoodPV, chi2ndf  );
    GoodPVtxChi2ProbVsGoodPVtx       -> Fill( totalNumGoodPV, chi2prob );
      }
    }
  }
}

void VertexMonitor::beginLuminosityBlock ( edm::LuminosityBlock const &  , edm::EventSetup const &    )

virtual

Definition at line 366 of file VertexMonitor.cc.

{
}

void VertexMonitor::endLuminosityBlock ( edm::LuminosityBlock const &  , edm::EventSetup const &    )

virtual

Definition at line 372 of file VertexMonitor.cc.

{
}

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

static

Definition at line 378 of file VertexMonitor.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                          {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

void VertexMonitor::initHisto ( DQMStore::IBooker &  ibooker )

virtual

Definition at line 181 of file VertexMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::bookProfile(), Chi2oNDFVsBXlumi, Chi2ProbVsBXlumi, conf_, doPlotsVsBXlumi_, doPlotsVsGoodPVtx_, FractionOfGoodPVtx, FractionOfGoodPVtxVsBXlumi, FractionOfGoodPVtxVsGoodPVtx, FractionOfGoodPVtxVsPVtx, edm::ParameterSet::getParameter(), MonitorElement::getTH1(), GoodPVtxChi2oNDFVsBXlumi, GoodPVtxChi2oNDFVsGoodPVtx, GoodPVtxChi2ProbVsBXlumi, GoodPVtxChi2ProbVsGoodPVtx, GoodPVtxNumberOfTracks, GoodPVtxNumberOfTracksVsBXlumi, GoodPVtxNumberOfTracksVsGoodPVtx, GoodPVtxSumPt, GoodPVtxSumPtVsBXlumi, GoodPVtxSumPtVsGoodPVtx, histname, label_, NumberOfBADndofPVtx, NumberOfBADndofPVtxVsBXlumi, NumberOfBADndofPVtxVsGoodPVtx, NumberOfGoodPVtx, NumberOfGoodPVtxVsBXlumi, NumberOfPVtx, NumberOfPVtxVsBXlumi, NumberOfPVtxVsGoodPVtx, MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
    // parameters from the configuration
    std::string MEFolderName   = conf_.getParameter<std::string>("PVFolderName"); 

    // get binning from the configuration
    int    GoodPVtxBin   = conf_.getParameter<int>("GoodPVtxBin");
    double GoodPVtxMin   = conf_.getParameter<double>("GoodPVtxMin");
    double GoodPVtxMax   = conf_.getParameter<double>("GoodPVtxMax");

    // book histo
    // ----------------------//
    ibooker.setCurrentFolder(MEFolderName+"/"+label_);

    histname = "NumberOfPVtx_" + label_;
    NumberOfPVtx = ibooker.book1D(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax);
    NumberOfPVtx->setAxisTitle("Number of PV",1);
    NumberOfPVtx->setAxisTitle("Number of Events",2);
    
    histname = "NumberOfGoodPVtx_" + label_;
    NumberOfGoodPVtx = ibooker.book1D(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax);
    NumberOfGoodPVtx->setAxisTitle("Number of Good PV",1);
    NumberOfGoodPVtx->setAxisTitle("Number of Events",2);
    
    histname = "FractionOfGoodPVtx_" + label_;
    FractionOfGoodPVtx = ibooker.book1D(histname,histname, 100,0.,1.);
    FractionOfGoodPVtx->setAxisTitle("fraction of Good PV",1);
    FractionOfGoodPVtx->setAxisTitle("Number of Events",2);
    
    histname = "NumberOfBADndofPVtx_" + label_;
    NumberOfBADndofPVtx = ibooker.book1D(histname,histname,GoodPVtxBin,GoodPVtxMin,GoodPVtxMax);
    NumberOfBADndofPVtx->setAxisTitle("Number of BADndof #PV",1);
    NumberOfBADndofPVtx->setAxisTitle("Number of Events",2);
    
    histname = "GoodPVtxSumPt_" + label_;
    GoodPVtxSumPt = ibooker.book1D(histname,histname,100,0.,500.);
    GoodPVtxSumPt->setAxisTitle("primary vertex #Sum p_{T}^{2} [GeV^{2}/c^{2}]",1);
    GoodPVtxSumPt->setAxisTitle("Number of events",2);

    histname = "GoodPVtxNumberOfTracks_" + label_;
    GoodPVtxNumberOfTracks = ibooker.book1D(histname,histname,100,0.,100.);
    GoodPVtxNumberOfTracks->setAxisTitle("primary vertex number of tracks",1);
    GoodPVtxNumberOfTracks->setAxisTitle("Number of events",2);

    if ( doPlotsVsBXlumi_ ) {
      // get binning from the configuration
      edm::ParameterSet BXlumiParameters = conf_.getParameter<edm::ParameterSet>("BXlumiSetup");
      int    BXlumiBin   = BXlumiParameters.getParameter<int>("BXlumiBin");
      double BXlumiMin   = BXlumiParameters.getParameter<double>("BXlumiMin");
      double BXlumiMax   = BXlumiParameters.getParameter<double>("BXlumiMax");

      ibooker.setCurrentFolder(MEFolderName+"/"+label_+"/PUmonitoring/");

      histname = "NumberOfPVtxVsBXlumi_" + label_;
      NumberOfPVtxVsBXlumi = ibooker.bookProfile(histname,histname, BXlumiBin,BXlumiMin,BXlumiMax,0.,60.,"");
      NumberOfPVtxVsBXlumi->getTH1()->SetCanExtend(TH1::kAllAxes);
      NumberOfPVtxVsBXlumi->setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      NumberOfPVtxVsBXlumi->setAxisTitle("Mean number of PV",2);

      histname = "NumberOfGoodPVtxVsBXlumi_" + label_;
      NumberOfGoodPVtxVsBXlumi = ibooker.bookProfile(histname,histname, BXlumiBin,BXlumiMin,BXlumiMax,0.,60.,"");
      NumberOfGoodPVtxVsBXlumi->getTH1()->SetCanExtend(TH1::kAllAxes);
      NumberOfGoodPVtxVsBXlumi->setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      NumberOfGoodPVtxVsBXlumi->setAxisTitle("Mean number of PV",2);

      histname = "FractionOfGoodPVtxVsBXlumi_" + label_;
      FractionOfGoodPVtxVsBXlumi = ibooker.bookProfile(histname,histname, BXlumiBin,BXlumiMin,BXlumiMax,0.,1.5,"");
      FractionOfGoodPVtxVsBXlumi->getTH1()->SetCanExtend(TH1::kAllAxes);
      FractionOfGoodPVtxVsBXlumi->setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      FractionOfGoodPVtxVsBXlumi->setAxisTitle("Mean number of PV",2);

      histname = "NumberOfBADndofPVtxVsBXlumi_" + label_;
      NumberOfBADndofPVtxVsBXlumi = ibooker.bookProfile(histname,histname, BXlumiBin,BXlumiMin,BXlumiMax,0.,60.,"");
      NumberOfBADndofPVtxVsBXlumi->getTH1()->SetCanExtend(TH1::kAllAxes);
      NumberOfBADndofPVtxVsBXlumi->setAxisTitle("BADndof #PV",1);
      NumberOfBADndofPVtxVsBXlumi->setAxisTitle("Number of Events",2);
    
      histname = "GoodPVtxSumPtVsBXlumi_" + label_;
      GoodPVtxSumPtVsBXlumi = ibooker.bookProfile(histname,histname, BXlumiBin,BXlumiMin,BXlumiMax,0.,500.,"");
      GoodPVtxSumPtVsBXlumi->getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxSumPtVsBXlumi->setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      GoodPVtxSumPtVsBXlumi->setAxisTitle("Mean pv #Sum p_{T}^{2} [GeV^{2}/c]^{2}",2);
      
      histname = "GoodPVtxNumberOfTracksVsBXlumi_" + label_;
      GoodPVtxNumberOfTracksVsBXlumi = ibooker.bookProfile(histname,histname, BXlumiBin,BXlumiMin,BXlumiMax,0.,100.,"");
      GoodPVtxNumberOfTracksVsBXlumi->getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxNumberOfTracksVsBXlumi->setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      GoodPVtxNumberOfTracksVsBXlumi->setAxisTitle("Mean pv number of tracks",2);

      // get binning from the configuration
      double Chi2NDFMin = conf_.getParameter<double>("Chi2NDFMin");
      double Chi2NDFMax = conf_.getParameter<double>("Chi2NDFMax");

      double Chi2ProbMin  = conf_.getParameter<double>("Chi2ProbMin");
      double Chi2ProbMax  = conf_.getParameter<double>("Chi2ProbMax");

      histname = "Chi2oNDFVsBXlumi_" + label_;
      Chi2oNDFVsBXlumi  = ibooker.bookProfile(histname,histname,BXlumiBin, BXlumiMin,BXlumiMax,Chi2NDFMin,Chi2NDFMax,"");
      Chi2oNDFVsBXlumi -> getTH1()->SetCanExtend(TH1::kAllAxes);
      Chi2oNDFVsBXlumi -> setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      Chi2oNDFVsBXlumi -> setAxisTitle("Mean #chi^{2}/ndof",2);

      histname = "Chi2ProbVsBXlumi_" + label_;
      Chi2ProbVsBXlumi  = ibooker.bookProfile(histname,histname,BXlumiBin, BXlumiMin,BXlumiMax,Chi2ProbMin,Chi2ProbMax,"");
      Chi2ProbVsBXlumi -> getTH1()->SetCanExtend(TH1::kAllAxes);
      Chi2ProbVsBXlumi -> setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      Chi2ProbVsBXlumi -> setAxisTitle("Mean #chi^{2}/prob",2);

      histname = "GoodPVtxChi2oNDFVsBXlumi_" + label_;
      GoodPVtxChi2oNDFVsBXlumi  = ibooker.bookProfile(histname,histname,BXlumiBin, BXlumiMin,BXlumiMax,Chi2NDFMin,Chi2NDFMax,"");
      GoodPVtxChi2oNDFVsBXlumi -> getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxChi2oNDFVsBXlumi -> setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      GoodPVtxChi2oNDFVsBXlumi -> setAxisTitle("Mean PV #chi^{2}/ndof",2);

      histname = "GoodPVtxChi2ProbVsBXlumi_" + label_;
      GoodPVtxChi2ProbVsBXlumi  = ibooker.bookProfile(histname,histname,BXlumiBin, BXlumiMin,BXlumiMax,Chi2ProbMin,Chi2ProbMax,"");
      GoodPVtxChi2ProbVsBXlumi -> getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxChi2ProbVsBXlumi -> setAxisTitle("lumi BX [10^{30}Hzcm^{-2}]",1);
      GoodPVtxChi2ProbVsBXlumi -> setAxisTitle("Mean PV #chi^{2}/prob",2);
    }

    if ( doPlotsVsGoodPVtx_ ) {

      ibooker.setCurrentFolder(MEFolderName+"/"+label_+"/PUmonitoring/VsGoodPVtx");

      histname = "NumberOfPVtxVsGoodPVtx_" + label_;
      NumberOfPVtxVsGoodPVtx = ibooker.bookProfile(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax,0.,60.,"");
      NumberOfPVtxVsGoodPVtx->getTH1()->SetCanExtend(TH1::kAllAxes);
      NumberOfPVtxVsGoodPVtx->setAxisTitle("Number of Good PV",1);
      NumberOfPVtxVsGoodPVtx->setAxisTitle("Mean number of PV",2);

      histname = "FractionOfGoodPVtxVsGoodPVtx_" + label_;
      FractionOfGoodPVtxVsGoodPVtx = ibooker.bookProfile(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax,0.,60.,"");
      FractionOfGoodPVtxVsGoodPVtx->getTH1()->SetCanExtend(TH1::kAllAxes);
      FractionOfGoodPVtxVsGoodPVtx->setAxisTitle("Number of Good PV",1);
      FractionOfGoodPVtxVsGoodPVtx->setAxisTitle("Mean fraction of Good PV",2);

      histname = "FractionOfGoodPVtxVsPVtx_" + label_;
      FractionOfGoodPVtxVsPVtx = ibooker.bookProfile(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax,0.,60.,"");
      FractionOfGoodPVtxVsPVtx->getTH1()->SetCanExtend(TH1::kAllAxes);
      FractionOfGoodPVtxVsPVtx->setAxisTitle("Number of Good PV",1);
      FractionOfGoodPVtxVsPVtx->setAxisTitle("Mean number of Good PV",2);

      histname = "NumberOfBADndofPVtxVsGoodPVtx_" + label_;
      NumberOfBADndofPVtxVsGoodPVtx = ibooker.bookProfile(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax,0.,60.,"");
      NumberOfBADndofPVtxVsGoodPVtx->getTH1()->SetCanExtend(TH1::kAllAxes);
      NumberOfBADndofPVtxVsGoodPVtx->setAxisTitle("Number of Good PV",1);
      NumberOfBADndofPVtxVsGoodPVtx->setAxisTitle("Mean Number of BAD PV",2);
    
      histname = "GoodPVtxSumPtVsGoodPVtx_" + label_;
      GoodPVtxSumPtVsGoodPVtx = ibooker.bookProfile(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax,0.,500.,"");
      GoodPVtxSumPtVsGoodPVtx->getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxSumPtVsGoodPVtx->setAxisTitle("Number of Good PV",1);
      GoodPVtxSumPtVsGoodPVtx->setAxisTitle("Mean pv #Sum p_{T}^{2} [GeV^{2}/c]^{2}",2);
      
      histname = "GoodPVtxNumberOfTracksVsGoodPVtx_" + label_;
      GoodPVtxNumberOfTracksVsGoodPVtx = ibooker.bookProfile(histname,histname, GoodPVtxBin,GoodPVtxMin,GoodPVtxMax,0.,100.,"");
      GoodPVtxNumberOfTracksVsGoodPVtx->getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxNumberOfTracksVsGoodPVtx->setAxisTitle("Number of Good PV",1);
      GoodPVtxNumberOfTracksVsGoodPVtx->setAxisTitle("Mean pv number of tracks",2);
      
      // get binning from the configuration
      double Chi2NDFMin = conf_.getParameter<double>("Chi2NDFMin");
      double Chi2NDFMax = conf_.getParameter<double>("Chi2NDFMax");

      double Chi2ProbMin  = conf_.getParameter<double>("Chi2ProbMin");
      double Chi2ProbMax  = conf_.getParameter<double>("Chi2ProbMax");

      histname = "GoodPVtxChi2oNDFVsGoodPVtx_" + label_;
      GoodPVtxChi2oNDFVsGoodPVtx  = ibooker.bookProfile(histname,histname,GoodPVtxBin, GoodPVtxMin,GoodPVtxMax,Chi2NDFMin,Chi2NDFMax,"");
      GoodPVtxChi2oNDFVsGoodPVtx -> getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxChi2oNDFVsGoodPVtx -> setAxisTitle("Number of Good PV",1);
      GoodPVtxChi2oNDFVsGoodPVtx -> setAxisTitle("Mean PV #chi^{2}/ndof",2);

      histname = "GoodPVtxChi2ProbVsGoodPVtx_" + label_;
      GoodPVtxChi2ProbVsGoodPVtx  = ibooker.bookProfile(histname,histname,GoodPVtxBin, GoodPVtxMin,GoodPVtxMax,Chi2ProbMin,Chi2ProbMax,"");
      GoodPVtxChi2ProbVsGoodPVtx -> getTH1()->SetCanExtend(TH1::kAllAxes);
      GoodPVtxChi2ProbVsGoodPVtx -> setAxisTitle("Number of Good PV",1);
      GoodPVtxChi2ProbVsGoodPVtx -> setAxisTitle("Mean PV #chi^{2}/prob",2);

    }
}

Member Data Documentation

MonitorElement* VertexMonitor::Chi2oNDFVsBXlumi

Definition at line 78 of file VertexMonitor.h.

Referenced by initHisto().

MonitorElement* VertexMonitor::Chi2oNDFVsGoodPVtx

Definition at line 77 of file VertexMonitor.h.

MonitorElement* VertexMonitor::Chi2ProbVsBXlumi

Definition at line 80 of file VertexMonitor.h.

Referenced by initHisto().

MonitorElement* VertexMonitor::Chi2ProbVsGoodPVtx

Definition at line 79 of file VertexMonitor.h.

edm::ParameterSet VertexMonitor::conf_

Definition at line 49 of file VertexMonitor.h.

Referenced by initHisto().

bool VertexMonitor::doAllPlots_

Definition at line 96 of file VertexMonitor.h.

bool VertexMonitor::doPlotsVsBXlumi_

Definition at line 97 of file VertexMonitor.h.

Referenced by analyze(), initHisto(), and VertexMonitor().

bool VertexMonitor::doPlotsVsGoodPVtx_

Definition at line 98 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::FractionOfGoodPVtx

Definition at line 66 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::FractionOfGoodPVtxVsBXlumi

Definition at line 67 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::FractionOfGoodPVtxVsGoodPVtx

Definition at line 68 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::FractionOfGoodPVtxVsPVtx

Definition at line 69 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxChi2oNDFVsBXlumi

Definition at line 92 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxChi2oNDFVsGoodPVtx

Definition at line 91 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxChi2ProbVsBXlumi

Definition at line 94 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxChi2ProbVsGoodPVtx

Definition at line 93 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxNumberOfTracks

Definition at line 86 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxNumberOfTracksVsBXlumi

Definition at line 87 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxNumberOfTracksVsGoodPVtx

Definition at line 88 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxNumberOfTracksVsGoodPVtxNdof

Definition at line 89 of file VertexMonitor.h.

MonitorElement* VertexMonitor::GoodPVtxSumPt

Definition at line 82 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxSumPtVsBXlumi

Definition at line 83 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::GoodPVtxSumPtVsGoodPVtx

Definition at line 84 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

std::string VertexMonitor::histname

Definition at line 100 of file VertexMonitor.h.

Referenced by initHisto().

std::string VertexMonitor::label_

Definition at line 53 of file VertexMonitor.h.

Referenced by Mixins._Labelable::_findDependencies(), Mixins._TypedParameterizable::dumpPython(), initHisto(), and Mixins._Labelable::setLabel().

GetLumi* VertexMonitor::lumiDetails_

Definition at line 59 of file VertexMonitor.h.

Referenced by analyze(), and VertexMonitor().

MonitorElement* VertexMonitor::NumberOfBADndofPVtx

Definition at line 73 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfBADndofPVtxVsBXlumi

Definition at line 74 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfBADndofPVtxVsGoodPVtx

Definition at line 75 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfFakePVtx

Definition at line 70 of file VertexMonitor.h.

MonitorElement* VertexMonitor::NumberOfFakePVtxVsBXlumi

Definition at line 71 of file VertexMonitor.h.

MonitorElement* VertexMonitor::NumberOfFakePVtxVsGoodPVtx

Definition at line 72 of file VertexMonitor.h.

MonitorElement* VertexMonitor::NumberOfGoodPVtx

Definition at line 64 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfGoodPVtxVsBXlumi

Definition at line 65 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfPVtx

Definition at line 61 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfPVtxVsBXlumi

Definition at line 62 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

MonitorElement* VertexMonitor::NumberOfPVtxVsGoodPVtx

Definition at line 63 of file VertexMonitor.h.

Referenced by analyze(), and initHisto().

edm::InputTag VertexMonitor::primaryVertexInputTag_

Definition at line 51 of file VertexMonitor.h.

Referenced by VertexMonitor().

edm::EDGetTokenT<reco::VertexCollection> VertexMonitor::pvToken_

Definition at line 55 of file VertexMonitor.h.

Referenced by analyze(), and VertexMonitor().

edm::InputTag VertexMonitor::selectedPrimaryVertexInputTag_

Definition at line 52 of file VertexMonitor.h.

Referenced by VertexMonitor().

edm::EDGetTokenT<reco::VertexCollection> VertexMonitor::selpvToken_

Definition at line 56 of file VertexMonitor.h.

Referenced by analyze(), and VertexMonitor().