VertexHistogramMaker Class Reference

#include <VertexHistogramMaker.h>

Public Member Functions
void	beginRun (const edm::Run &iRun)
void	book (const std::string dirname="")
void	fill (const edm::Event &iEvent, const reco::VertexCollection &vertices, const double weight=1.)
	VertexHistogramMaker (edm::ConsumesCollector &&iC)
	VertexHistogramMaker (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC)
	~VertexHistogramMaker ()

Private Member Functions
void	fill (const unsigned int orbit, const int bx, const float bxlumi, const reco::VertexCollection &vertices, const double weight=1.)

Private Attributes
const bool	m_bsConstrained
TFileDirectory *	m_currdir
RunHistogramManager	m_fhm
TH1F *	m_haveweight
TProfile *	m_haveweightvsvtxz
TProfile *	m_haveweightvsvtxzchk
const edm::ParameterSet	m_histoParameters
TH1F *	m_hndof
TH2F *	m_hndofvstracks
TProfile *	m_hndofvsvtxz
TH1F *	m_hnheavytracks
TH1F *	m_hntracks
TProfile *	m_hntracksvsvtxz
TH1F *	m_hntruevtx
TProfile *	m_hntruevtxvslumi
TH2D *	m_hntruevtxvslumi2D
TH1F *	m_hnvtx
TProfile **	m_hnvtxvsbxrun
TProfile2D **	m_hnvtxvsbxvslumirun
TProfile2D **	m_hnvtxvsbxvsorbrun
TProfile **	m_hnvtxvsorbrun
TH2F **	m_hnvtxvsorbrun2D
TH1F *	m_hsqsumptsq
TH1F *	m_hsqsumptsqheavy
TH1F *	m_hvtxx
TH1F **	m_hvtxxrun
TH2F **	m_hvtxxvsbx2drun
TProfile **	m_hvtxxvsbxrun
TProfile **	m_hvtxxvsorbrun
TH1F *	m_hvtxy
TH1F **	m_hvtxyrun
TH2F **	m_hvtxyvsbx2drun
TProfile **	m_hvtxyvsbxrun
TProfile **	m_hvtxyvsorbrun
TH1F *	m_hvtxz
TH1F **	m_hvtxzrun
TH2F **	m_hvtxzvsbx2drun
TProfile **	m_hvtxzvsbxrun
TProfile **	m_hvtxzvsorbrun
TH1F *	m_hweights
edm::EDGetTokenT< LumiDetails >	m_lumiDetailsToken
const unsigned int	m_maxLS
RunHistogramManager	m_rhm
const bool	m_runHisto
const bool	m_runHisto2D
const bool	m_runHistoBXProfile
const bool	m_runHistoBXProfile2D
const bool	m_runHistoProfile
const bool	m_trueOnly
const double	m_weightThreshold

Detailed Description

Definition at line 23 of file VertexHistogramMaker.h.

Constructor & Destructor Documentation

VertexHistogramMaker::VertexHistogramMaker

(

edm::ConsumesCollector &&

iC

)

Definition at line 15 of file VertexHistogramMaker.cc.

    : m_currdir(nullptr),
      m_maxLS(100),
      m_weightThreshold(0.5),
      m_trueOnly(true),
      m_runHisto(true),
      m_runHistoProfile(true),
      m_runHistoBXProfile(true),
      m_runHistoBXProfile2D(false),
      m_runHisto2D(false),
      m_bsConstrained(false),
      m_histoParameters(),
      m_rhm(iC),
      m_fhm(iC) {}

VertexHistogramMaker::VertexHistogramMaker	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 30 of file VertexHistogramMaker.cc.

    : m_currdir(nullptr),
      m_maxLS(iConfig.getParameter<unsigned int>("maxLSBeforeRebin")),
      m_weightThreshold(iConfig.getUntrackedParameter<double>("weightThreshold", 0.5)),
      m_trueOnly(iConfig.getUntrackedParameter<bool>("trueOnly", true)),
      m_runHisto(iConfig.getUntrackedParameter<bool>("runHisto", true)),
      m_runHistoProfile(iConfig.getUntrackedParameter<bool>("runHistoProfile", true)),
      m_runHistoBXProfile(iConfig.getUntrackedParameter<bool>("runHistoBXProfile", true)),
      m_runHistoBXProfile2D(iConfig.getUntrackedParameter<bool>("runHistoBXProfile2D", false)),
      m_runHisto2D(iConfig.getUntrackedParameter<bool>("runHisto2D", false)),
      m_bsConstrained(iConfig.getParameter<bool>("bsConstrained")),
      m_histoParameters(iConfig.getUntrackedParameter<edm::ParameterSet>("histoParameters", edm::ParameterSet())),
      m_lumiDetailsToken(iC.consumes<LumiDetails, edm::InLumi>(edm::InputTag(std::string("lumiProducer")))),
      m_rhm(iC, false),
      m_fhm(iC, true) {}

VertexHistogramMaker::~VertexHistogramMaker

(

)

Definition at line 46 of file VertexHistogramMaker.cc.

References m_currdir.

{ delete m_currdir; }

Member Function Documentation

void VertexHistogramMaker::beginRun

(

const edm::Run &

iRun

)

Definition at line 257 of file VertexHistogramMaker.cc.

References RunHistogramManager::beginRun(), m_currdir, m_fhm, m_rhm, m_runHisto, m_runHisto2D, m_runHistoBXProfile, m_runHistoBXProfile2D, m_runHistoProfile, and TFileService::tFileDirectory().

Referenced by AnotherPrimaryVertexAnalyzer::beginRun().

                                                      {
  TFileDirectory* currdir = m_currdir;
  if (currdir == nullptr) {
    edm::Service<TFileService> tfserv;
    currdir = &(tfserv->tFileDirectory());
  }

  m_rhm.beginRun(iRun, *currdir);
  m_fhm.beginRun(iRun, *currdir);

  if (m_runHisto) {
    (*m_hvtxxrun)->GetXaxis()->SetTitle("X [cm]");
    (*m_hvtxxrun)->GetYaxis()->SetTitle("Vertices");
    (*m_hvtxyrun)->GetXaxis()->SetTitle("Y [cm]");
    (*m_hvtxyrun)->GetYaxis()->SetTitle("Vertices");
    (*m_hvtxzrun)->GetXaxis()->SetTitle("Z [cm]");
    (*m_hvtxzrun)->GetYaxis()->SetTitle("Vertices");

    if (m_runHistoProfile) {
      (*m_hvtxxvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
      (*m_hvtxxvsorbrun)->GetYaxis()->SetTitle("X [cm]");
      (*m_hvtxxvsorbrun)->SetCanExtend(TH1::kAllAxes);
      (*m_hvtxyvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
      (*m_hvtxyvsorbrun)->GetYaxis()->SetTitle("Y [cm]");
      (*m_hvtxyvsorbrun)->SetCanExtend(TH1::kAllAxes);
      (*m_hvtxzvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
      (*m_hvtxzvsorbrun)->GetYaxis()->SetTitle("Z [cm]");
      (*m_hvtxzvsorbrun)->SetCanExtend(TH1::kAllAxes);
      (*m_hnvtxvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
      (*m_hnvtxvsorbrun)->GetYaxis()->SetTitle("Nvertices");
      (*m_hnvtxvsorbrun)->SetCanExtend(TH1::kAllAxes);
    }

    if (m_runHistoBXProfile) {
      (*m_hvtxxvsbxrun)->GetXaxis()->SetTitle("BX");
      (*m_hvtxxvsbxrun)->GetYaxis()->SetTitle("X [cm]");
      (*m_hvtxyvsbxrun)->GetXaxis()->SetTitle("BX");
      (*m_hvtxyvsbxrun)->GetYaxis()->SetTitle("Y [cm]");
      (*m_hvtxzvsbxrun)->GetXaxis()->SetTitle("BX");
      (*m_hvtxzvsbxrun)->GetYaxis()->SetTitle("Z [cm]");
      (*m_hnvtxvsbxrun)->GetXaxis()->SetTitle("BX");
      (*m_hnvtxvsbxrun)->GetYaxis()->SetTitle("Nvertices");
      if (m_runHistoBXProfile2D) {
        (*m_hnvtxvsbxvslumirun)->GetXaxis()->SetTitle("BX");
        (*m_hnvtxvsbxvslumirun)->GetYaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
      }
      if (m_runHisto2D) {
        (*m_hvtxxvsbx2drun)->GetXaxis()->SetTitle("BX");
        (*m_hvtxxvsbx2drun)->GetYaxis()->SetTitle("X [cm]");
        (*m_hvtxyvsbx2drun)->GetXaxis()->SetTitle("BX");
        (*m_hvtxyvsbx2drun)->GetYaxis()->SetTitle("Y [cm]");
        (*m_hvtxzvsbx2drun)->GetXaxis()->SetTitle("BX");
        (*m_hvtxzvsbx2drun)->GetYaxis()->SetTitle("Z [cm]");
      }
    }

    if (m_runHisto2D) {
      (*m_hnvtxvsbxvsorbrun)->GetXaxis()->SetTitle("BX#");
      (*m_hnvtxvsbxvsorbrun)->GetYaxis()->SetTitle("time [orbit#]");
      (*m_hnvtxvsbxvsorbrun)->SetCanExtend(TH1::kAllAxes);
      (*m_hnvtxvsorbrun2D)->GetXaxis()->SetTitle("time [orbit#]");
      (*m_hnvtxvsorbrun2D)->GetYaxis()->SetTitle("Nvertices");
      (*m_hnvtxvsorbrun2D)->SetCanExtend(TH1::kAllAxes);
    }
  }
}

void VertexHistogramMaker::book

(

const std::string

dirname = ""

)

Definition at line 48 of file VertexHistogramMaker.cc.

References BTVHLTOfflineSource_cfi::dirname, edm::ParameterSet::getUntrackedParameter(), m_currdir, m_fhm, m_haveweight, m_haveweightvsvtxz, m_haveweightvsvtxzchk, m_histoParameters, m_hndof, m_hndofvstracks, m_hndofvsvtxz, m_hnheavytracks, m_hntracks, m_hntracksvsvtxz, m_hntruevtx, m_hntruevtxvslumi, m_hntruevtxvslumi2D, m_hnvtx, m_hnvtxvsbxrun, m_hnvtxvsbxvslumirun, m_hnvtxvsbxvsorbrun, m_hnvtxvsorbrun, m_hnvtxvsorbrun2D, m_hsqsumptsq, m_hsqsumptsqheavy, m_hvtxx, m_hvtxxrun, m_hvtxxvsbx2drun, m_hvtxxvsbxrun, m_hvtxxvsorbrun, m_hvtxy, m_hvtxyrun, m_hvtxyvsbx2drun, m_hvtxyvsbxrun, m_hvtxyvsorbrun, m_hvtxz, m_hvtxzrun, m_hvtxzvsbx2drun, m_hvtxzvsbxrun, m_hvtxzvsorbrun, m_hweights, m_maxLS, m_rhm, m_runHisto, m_runHisto2D, m_runHistoBXProfile, m_runHistoBXProfile2D, m_runHistoProfile, m_weightThreshold, TFileDirectory::make(), RunHistogramManager::makeTH1F(), RunHistogramManager::makeTH2F(), RunHistogramManager::makeTProfile(), RunHistogramManager::makeTProfile2D(), TFileService::mkdir(), and TFileService::tFileDirectory().

Referenced by AnotherPrimaryVertexAnalyzer::AnotherPrimaryVertexAnalyzer().

                                                       {
  edm::Service<TFileService> tfserv;
  TFileDirectory* currdir = &(tfserv->tFileDirectory());

  if (!dirname.empty()) {
    currdir = new TFileDirectory(tfserv->mkdir(dirname));
    m_currdir = currdir;
  }

  edm::LogInfo("HistogramBooking") << "Vertex histogram booking in directory " << dirname;

  m_hnvtx = currdir->make<TH1F>("nvtx", "Number of Vertices", 60, -0.5, 59.5);
  m_hnvtx->GetXaxis()->SetTitle("vertices");
  m_hnvtx->GetYaxis()->SetTitle("Events");

  m_hntruevtx = currdir->make<TH1F>("ntruevtx", "Number of True Vertices", 60, -0.5, 59.5);
  m_hntruevtx->GetXaxis()->SetTitle("vertices");
  m_hntruevtx->GetYaxis()->SetTitle("Events");

  m_hntruevtxvslumi = currdir->make<TProfile>("ntruevtxvslumi", "Number of True Vertices vs BX lumi", 250, 0., 10.);
  m_hntruevtxvslumi->GetXaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
  m_hntruevtxvslumi->GetYaxis()->SetTitle("Vertices");

  m_hntruevtxvslumi2D =
      currdir->make<TH2D>("ntruevtxvslumi2D", "Number of True Vertices vs BX lumi", 250, 0., 10., 100, -0.5, 99.5);
  m_hntruevtxvslumi2D->GetXaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
  m_hntruevtxvslumi2D->GetYaxis()->SetTitle("Vertices");

  m_hntracks = currdir->make<TH1F>("ntracks", "Number of Tracks", 300, -0.5, 299.5);
  m_hntracks->GetXaxis()->SetTitle("tracks");
  m_hntracks->GetYaxis()->SetTitle("Vertices");

  m_hsqsumptsq = currdir->make<TH1F>("sqsumptsq", "sqrt(sum pt**2)", 1000, 0., 1000.);
  m_hsqsumptsq->GetXaxis()->SetTitle("sqrt(#Sigma pt^{2}) (GeV)");
  m_hsqsumptsq->GetYaxis()->SetTitle("Vertices");

  char htitle[300];
  sprintf(htitle, "sqrt(sum pt**2) of Tracks weight > %f", m_weightThreshold);
  m_hsqsumptsqheavy = currdir->make<TH1F>("sqsumptsqheavy", htitle, 1000, 0., 1000.);
  m_hsqsumptsqheavy->GetXaxis()->SetTitle("sqrt(#Sigma pt^{2}) (GeV)");
  m_hsqsumptsqheavy->GetYaxis()->SetTitle("Vertices");

  sprintf(htitle, "Number of Tracks weight > %f", m_weightThreshold);
  m_hnheavytracks = currdir->make<TH1F>("nheavytracks", htitle, 200, -0.5, 199.5);
  m_hnheavytracks->GetXaxis()->SetTitle("tracks");
  m_hnheavytracks->GetYaxis()->SetTitle("Vertices");

  m_hndof = currdir->make<TH1F>("ndof", "Number of degree of freedom", 250, -0.5, 499.5);
  m_hndof->GetXaxis()->SetTitle("ndof");
  m_hndof->GetYaxis()->SetTitle("Vertices");

  m_hndofvstracks = currdir->make<TH2F>("ndofvstracks", "Ndof vs Ntracks", 300, -0.5, 299.5, 250, -0.5, 499.5);
  m_hndofvstracks->GetXaxis()->SetTitle("tracks");
  m_hndofvstracks->GetYaxis()->SetTitle("ndof");

  m_hndofvsvtxz = currdir->make<TProfile>("ndofvsvtxz",
                                          "Ndof vs Vertex Z position",
                                          200,
                                          m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                          m_histoParameters.getUntrackedParameter<double>("zMax", 20.));
  m_hndofvsvtxz->GetXaxis()->SetTitle("Z [cm]");
  m_hndofvsvtxz->GetYaxis()->SetTitle("ndof");

  m_hntracksvsvtxz = currdir->make<TProfile>("ntracksvsvtxz",
                                             "Ntracks vs Vertex Z position",
                                             200,
                                             m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                             m_histoParameters.getUntrackedParameter<double>("zMax", 20.));
  m_hntracksvsvtxz->GetXaxis()->SetTitle("Z [cm]");
  m_hntracksvsvtxz->GetYaxis()->SetTitle("tracks");

  m_haveweightvsvtxz = currdir->make<TProfile>("aveweightvsvtxz",
                                               "Average weight vs Vertex Z position",
                                               200,
                                               m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                               m_histoParameters.getUntrackedParameter<double>("zMax", 20.));
  m_haveweightvsvtxz->GetXaxis()->SetTitle("Z [cm]");
  m_haveweightvsvtxz->GetYaxis()->SetTitle("Average weight");

  m_haveweightvsvtxzchk = currdir->make<TProfile>("aveweightvsvtxzchk",
                                                  "Average weight vs Vertex Z position (check)",
                                                  200,
                                                  m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                                  m_histoParameters.getUntrackedParameter<double>("zMax", 20.));
  m_haveweightvsvtxzchk->GetXaxis()->SetTitle("Z [cm]");
  m_haveweightvsvtxzchk->GetYaxis()->SetTitle("Average weight");

  m_hweights = currdir->make<TH1F>("weights", "Tracks weights", 51, 0., 1.02);
  m_hweights->GetXaxis()->SetTitle("weights");
  m_hweights->GetYaxis()->SetTitle("Tracks");

  m_haveweight = currdir->make<TH1F>("aveweight", "Tracks average weights sum", 51, 0., 1.02);
  m_haveweight->GetXaxis()->SetTitle("Average weight");
  m_haveweight->GetYaxis()->SetTitle("Vertices");

  m_hvtxx = currdir->make<TH1F>("vtxx",
                                "Vertex X position",
                                m_histoParameters.getUntrackedParameter<unsigned int>("nBinX", 200),
                                m_histoParameters.getUntrackedParameter<double>("xMin", -1.),
                                m_histoParameters.getUntrackedParameter<double>("xMax", 1.));
  m_hvtxx->GetXaxis()->SetTitle("X [cm]");
  m_hvtxx->GetYaxis()->SetTitle("Vertices");

  m_hvtxy = currdir->make<TH1F>("vtxy",
                                "Vertex Y position",
                                m_histoParameters.getUntrackedParameter<unsigned int>("nBinY", 200),
                                m_histoParameters.getUntrackedParameter<double>("yMin", -1.),
                                m_histoParameters.getUntrackedParameter<double>("yMax", 1.));
  m_hvtxy->GetXaxis()->SetTitle("Y [cm]");
  m_hvtxy->GetYaxis()->SetTitle("Vertices");

  m_hvtxz = currdir->make<TH1F>("vtxz",
                                "Vertex Z position",
                                m_histoParameters.getUntrackedParameter<unsigned int>("nBinZ", 200),
                                m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                m_histoParameters.getUntrackedParameter<double>("zMax", 20.));
  m_hvtxz->GetXaxis()->SetTitle("Z [cm]");
  m_hvtxz->GetYaxis()->SetTitle("Vertices");

  if (m_runHisto) {
    m_hvtxxrun = m_rhm.makeTH1F("vtxxrun",
                                "Vertex X position",
                                m_histoParameters.getUntrackedParameter<unsigned int>("nBinX", 200),
                                m_histoParameters.getUntrackedParameter<double>("xMin", -1.),
                                m_histoParameters.getUntrackedParameter<double>("xMax", 1.));

    m_hvtxyrun = m_rhm.makeTH1F("vtxyrun",
                                "Vertex Y position",
                                m_histoParameters.getUntrackedParameter<unsigned int>("nBinY", 200),
                                m_histoParameters.getUntrackedParameter<double>("yMin", -1.),
                                m_histoParameters.getUntrackedParameter<double>("yMax", 1.));

    m_hvtxzrun = m_rhm.makeTH1F("vtxzrun",
                                "Vertex Z position",
                                m_histoParameters.getUntrackedParameter<unsigned int>("nBinZ", 200),
                                m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                m_histoParameters.getUntrackedParameter<double>("zMax", 20.));

    if (m_runHistoProfile) {
      m_hvtxxvsorbrun = m_rhm.makeTProfile(
          "vtxxvsorbrun", "Vertex X position vs orbit number", 4 * m_maxLS, 0.5, m_maxLS * 262144 + 0.5);
      m_hvtxyvsorbrun = m_rhm.makeTProfile(
          "vtxyvsorbrun", "Vertex Y position vs orbit number", 4 * m_maxLS, 0.5, m_maxLS * 262144 + 0.5);
      m_hvtxzvsorbrun = m_rhm.makeTProfile(
          "vtxzvsorbrun", "Vertex Z position vs orbit number", 4 * m_maxLS, 0.5, m_maxLS * 262144 + 0.5);
      m_hnvtxvsorbrun = m_rhm.makeTProfile(
          "nvtxvsorbrun", "Number of true vertices vs orbit number", m_maxLS, 0.5, m_maxLS * 262144 + 0.5);
    }

    if (m_runHisto2D) {
      m_hnvtxvsbxvsorbrun = m_rhm.makeTProfile2D("nvtxvsbxvsorbrun",
                                                 "Number of true vertices vs BX vs orbit number",
                                                 3564,
                                                 -0.5,
                                                 3563.5,
                                                 m_maxLS,
                                                 0.5,
                                                 m_maxLS * 262144 + 0.5);
      m_hnvtxvsorbrun2D = m_rhm.makeTH2F("nvtxvsorbrun2D",
                                         "Number of true vertices vs orbit number",
                                         m_maxLS,
                                         0.5,
                                         m_maxLS * 262144 + 0.5,
                                         60,
                                         -0.5,
                                         59.5);
    }

    if (m_runHistoBXProfile) {
      m_hvtxxvsbxrun = m_fhm.makeTProfile("vtxxvsbxrun", "Vertex X position vs BX number", 3564, -0.5, 3563.5);
      m_hvtxyvsbxrun = m_fhm.makeTProfile("vtxyvsbxrun", "Vertex Y position vs BX number", 3564, -0.5, 3563.5);
      m_hvtxzvsbxrun = m_fhm.makeTProfile("vtxzvsbxrun", "Vertex Z position vs BX number", 3564, -0.5, 3563.5);

      m_hnvtxvsbxrun = m_rhm.makeTProfile("nvtxvsbxrun", "Number of true vertices vs BX number", 3564, -0.5, 3563.5);

      if (m_runHistoBXProfile2D) {
        m_hnvtxvsbxvslumirun = m_fhm.makeTProfile2D(
            "nvtxvsbxvslumirun", "Number of vertices vs BX and BX lumi", 3564, -0.5, 3563.5, 250, 0., 10.);
      }
      if (m_runHisto2D) {
        m_hvtxxvsbx2drun = m_fhm.makeTH2F("vtxxvsbx2drun",
                                          "Vertex X position vs BX number",
                                          3564,
                                          -0.5,
                                          3563.5,
                                          m_histoParameters.getUntrackedParameter<unsigned int>("nBinX", 200),
                                          m_histoParameters.getUntrackedParameter<double>("xMin", -1.),
                                          m_histoParameters.getUntrackedParameter<double>("xMax", 1.));
        m_hvtxyvsbx2drun = m_fhm.makeTH2F("vtxyvsbx2drun",
                                          "Vertex Y position vs BX number",
                                          3564,
                                          -0.5,
                                          3563.5,
                                          m_histoParameters.getUntrackedParameter<unsigned int>("nBinY", 200),
                                          m_histoParameters.getUntrackedParameter<double>("yMin", -1.),
                                          m_histoParameters.getUntrackedParameter<double>("yMax", 1.));
        m_hvtxzvsbx2drun = m_fhm.makeTH2F("vtxzvsbx2drun",
                                          "Vertex Z position vs BX number",
                                          3564,
                                          -0.5,
                                          3563.5,
                                          m_histoParameters.getUntrackedParameter<unsigned int>("nBinZ", 200),
                                          m_histoParameters.getUntrackedParameter<double>("zMin", -20.),
                                          m_histoParameters.getUntrackedParameter<double>("zMax", 20.));
      }
    }
  }
}

void VertexHistogramMaker::fill	(	const edm::Event &	iEvent,
		const reco::VertexCollection &	vertices,
		const double	weight = 1.
	)

Definition at line 440 of file VertexHistogramMaker.cc.

References edm::EventBase::bunchCrossing(), edm::LuminosityBlock::getByToken(), edm::Event::getLuminosityBlock(), LumiDetails::isValid(), edm::HandleBase::isValid(), LumiDetails::kOCC1, LumiDetails::lumiValue(), m_lumiDetailsToken, edm::EventBase::orbitNumber(), pwdgSkimBPark_cfi::vertices, and mps_merge::weight.

Referenced by AnotherPrimaryVertexAnalyzer::analyze().

                                                                                                                 {
  // get luminosity

  edm::Handle<LumiDetails> ld;
  iEvent.getLuminosityBlock().getByToken(m_lumiDetailsToken, ld);

  float bxlumi = -1.;

  if (ld.isValid()) {
    if (ld->isValid()) {
      bxlumi = ld->lumiValue(LumiDetails::kOCC1, iEvent.bunchCrossing()) * 6.37;
    }
  }

  fill(iEvent.orbitNumber(), iEvent.bunchCrossing(), bxlumi, vertices, weight);
}

void VertexHistogramMaker::fill ( const unsigned int  orbit, const int  bx, const float  bxlumi, const reco::VertexCollection &  vertices, const double  weight = 1.  )

private

Definition at line 324 of file VertexHistogramMaker.cc.

References m_bsConstrained, m_haveweight, m_haveweightvsvtxz, m_haveweightvsvtxzchk, m_hndof, m_hndofvstracks, m_hndofvsvtxz, m_hnheavytracks, m_hntracks, m_hntracksvsvtxz, m_hntruevtx, m_hntruevtxvslumi, m_hntruevtxvslumi2D, m_hnvtx, m_hnvtxvsbxrun, m_hnvtxvsbxvslumirun, m_hnvtxvsbxvsorbrun, m_hnvtxvsorbrun, m_hnvtxvsorbrun2D, m_hsqsumptsq, m_hsqsumptsqheavy, m_hvtxx, m_hvtxxrun, m_hvtxxvsbx2drun, m_hvtxxvsbxrun, m_hvtxxvsorbrun, m_hvtxy, m_hvtxyrun, m_hvtxyvsbx2drun, m_hvtxyvsbxrun, m_hvtxyvsorbrun, m_hvtxz, m_hvtxzrun, m_hvtxzvsbx2drun, m_hvtxzvsbxrun, m_hvtxzvsorbrun, m_hweights, m_runHisto, m_runHisto2D, m_runHistoBXProfile, m_runHistoBXProfile2D, m_runHistoProfile, m_trueOnly, m_weightThreshold, mathSSE::sqrt(), badGlobalMuonTaggersAOD_cff::vtx, and mps_merge::weight.

                                                     {
  m_hnvtx->Fill(vertices.size(), weight);

  int ntruevtx = 0;
  for (reco::VertexCollection::const_iterator vtx = vertices.begin(); vtx != vertices.end(); ++vtx) {
    if (!vtx->isFake())
      ntruevtx++;

    if (!(m_trueOnly && vtx->isFake())) {
      double aveweight =
          m_bsConstrained ? vtx->ndof() / (2. * vtx->tracksSize()) : (vtx->ndof() + 3) / (2. * vtx->tracksSize());

      m_hntracks->Fill(vtx->tracksSize(), weight);
      m_hndof->Fill(vtx->ndof(), weight);
      m_haveweight->Fill(aveweight, weight);
      m_hndofvstracks->Fill(vtx->tracksSize(), vtx->ndof(), weight);
      m_hndofvsvtxz->Fill(vtx->z(), vtx->ndof(), weight);
      m_hntracksvsvtxz->Fill(vtx->z(), vtx->tracksSize(), weight);
      m_haveweightvsvtxz->Fill(vtx->z(), aveweight, weight);

      m_hvtxx->Fill(vtx->x(), weight);
      m_hvtxy->Fill(vtx->y(), weight);
      m_hvtxz->Fill(vtx->z(), weight);

      if (m_runHisto) {
        if (m_hvtxxrun && *m_hvtxxrun)
          (*m_hvtxxrun)->Fill(vtx->x(), weight);
        if (m_hvtxyrun && *m_hvtxyrun)
          (*m_hvtxyrun)->Fill(vtx->y(), weight);
        if (m_hvtxzrun && *m_hvtxzrun)
          (*m_hvtxzrun)->Fill(vtx->z(), weight);
        if (m_runHistoProfile) {
          if (m_hvtxxvsorbrun && *m_hvtxxvsorbrun)
            (*m_hvtxxvsorbrun)->Fill(orbit, vtx->x(), weight);
          if (m_hvtxyvsorbrun && *m_hvtxyvsorbrun)
            (*m_hvtxyvsorbrun)->Fill(orbit, vtx->y(), weight);
          if (m_hvtxzvsorbrun && *m_hvtxzvsorbrun)
            (*m_hvtxzvsorbrun)->Fill(orbit, vtx->z(), weight);
        }
        if (m_runHistoBXProfile) {
          if (m_hvtxxvsbxrun && *m_hvtxxvsbxrun)
            (*m_hvtxxvsbxrun)->Fill(bx % 3564, vtx->x(), weight);
          if (m_hvtxyvsbxrun && *m_hvtxyvsbxrun)
            (*m_hvtxyvsbxrun)->Fill(bx % 3564, vtx->y(), weight);
          if (m_hvtxzvsbxrun && *m_hvtxzvsbxrun)
            (*m_hvtxzvsbxrun)->Fill(bx % 3564, vtx->z(), weight);
          if (m_runHisto2D) {
            if (m_hvtxxvsbx2drun && *m_hvtxxvsbx2drun)
              (*m_hvtxxvsbx2drun)->Fill(bx % 3564, vtx->x(), weight);
            if (m_hvtxyvsbx2drun && *m_hvtxyvsbx2drun)
              (*m_hvtxyvsbx2drun)->Fill(bx % 3564, vtx->y(), weight);
            if (m_hvtxzvsbx2drun && *m_hvtxzvsbx2drun)
              (*m_hvtxzvsbx2drun)->Fill(bx % 3564, vtx->z(), weight);
          }
        }
      }

      int nheavytracks = 0;
      double sumpt2 = 0.;
      double sumpt2heavy = 0.;

      for (reco::Vertex::trackRef_iterator trk = vtx->tracks_begin(); trk != vtx->tracks_end(); ++trk) {
        sumpt2 += (*trk)->pt() * (*trk)->pt();

        if (vtx->trackWeight(*trk) > m_weightThreshold) {
          nheavytracks++;
          sumpt2heavy += (*trk)->pt() * (*trk)->pt();
        }

        m_hweights->Fill(vtx->trackWeight(*trk), weight);
        m_haveweightvsvtxzchk->Fill(vtx->z(), vtx->trackWeight(*trk), weight);
      }

      m_hnheavytracks->Fill(nheavytracks, weight);
      m_hsqsumptsq->Fill(sqrt(sumpt2), weight);
      m_hsqsumptsqheavy->Fill(sqrt(sumpt2heavy), weight);
    }
  }

  m_hntruevtx->Fill(ntruevtx, weight);

  if (bxlumi >= 0.) {
    m_hntruevtxvslumi->Fill(bxlumi, ntruevtx, weight);
    m_hntruevtxvslumi2D->Fill(bxlumi, ntruevtx, weight);
  }

  if (m_runHisto) {
    if (m_runHistoProfile) {
      if (m_hnvtxvsorbrun && *m_hnvtxvsorbrun)
        (*m_hnvtxvsorbrun)->Fill(orbit, ntruevtx, weight);
    }
    if (m_runHistoBXProfile) {
      if (m_hnvtxvsbxrun && *m_hnvtxvsbxrun)
        (*m_hnvtxvsbxrun)->Fill(bx % 3564, ntruevtx, weight);
      if (m_runHistoBXProfile2D) {
        if (m_hnvtxvsbxvslumirun && *m_hnvtxvsbxvslumirun && bxlumi >= 0.)
          (*m_hnvtxvsbxvslumirun)->Fill(bx % 3564, bxlumi, ntruevtx, weight);
      }
    }
    if (m_runHisto2D) {
      if (m_hnvtxvsbxvsorbrun && *m_hnvtxvsbxvsorbrun)
        (*m_hnvtxvsbxvsorbrun)->Fill(bx % 3564, orbit, ntruevtx, weight);
      if (m_hnvtxvsorbrun2D && *m_hnvtxvsorbrun2D) {
        if (ntruevtx < (*m_hnvtxvsorbrun2D)->GetYaxis()->GetXmax() &&
            ntruevtx > (*m_hnvtxvsorbrun2D)->GetYaxis()->GetXmin()) {
          (*m_hnvtxvsorbrun2D)->Fill(orbit, ntruevtx, weight);
        }
      }
    }
  }
}

Member Data Documentation

const bool VertexHistogramMaker::m_bsConstrained

private

Definition at line 50 of file VertexHistogramMaker.h.

Referenced by fill().

TFileDirectory* VertexHistogramMaker::m_currdir

private

Definition at line 41 of file VertexHistogramMaker.h.

Referenced by beginRun(), book(), and ~VertexHistogramMaker().

RunHistogramManager VertexHistogramMaker::m_fhm

private

Definition at line 55 of file VertexHistogramMaker.h.

Referenced by beginRun(), and book().

TH1F* VertexHistogramMaker::m_haveweight

private

Definition at line 65 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile* VertexHistogramMaker::m_haveweightvsvtxz

private

Definition at line 69 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile* VertexHistogramMaker::m_haveweightvsvtxzchk

private

Definition at line 70 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

const edm::ParameterSet VertexHistogramMaker::m_histoParameters

private

Definition at line 51 of file VertexHistogramMaker.h.

Referenced by book().

TH1F* VertexHistogramMaker::m_hndof

private

Definition at line 64 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH2F* VertexHistogramMaker::m_hndofvstracks

private

Definition at line 66 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile* VertexHistogramMaker::m_hndofvsvtxz

private

Definition at line 67 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hnheavytracks

private

Definition at line 63 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hntracks

private

Definition at line 60 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile* VertexHistogramMaker::m_hntracksvsvtxz

private

Definition at line 68 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hntruevtx

private

Definition at line 57 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile* VertexHistogramMaker::m_hntruevtxvslumi

private

Definition at line 58 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH2D* VertexHistogramMaker::m_hntruevtxvslumi2D

private

Definition at line 59 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hnvtx

private

Definition at line 56 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hnvtxvsbxrun

private

Definition at line 88 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile2D** VertexHistogramMaker::m_hnvtxvsbxvslumirun

private

Definition at line 90 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile2D** VertexHistogramMaker::m_hnvtxvsbxvsorbrun

private

Definition at line 82 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hnvtxvsorbrun

private

Definition at line 81 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH2F** VertexHistogramMaker::m_hnvtxvsorbrun2D

private

Definition at line 83 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hsqsumptsq

private

Definition at line 61 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hsqsumptsqheavy

private

Definition at line 62 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hvtxx

private

Definition at line 72 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F** VertexHistogramMaker::m_hvtxxrun

private

Definition at line 75 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH2F** VertexHistogramMaker::m_hvtxxvsbx2drun

private

Definition at line 92 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hvtxxvsbxrun

private

Definition at line 85 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hvtxxvsorbrun

private

Definition at line 78 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hvtxy

private

Definition at line 73 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F** VertexHistogramMaker::m_hvtxyrun

private

Definition at line 76 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH2F** VertexHistogramMaker::m_hvtxyvsbx2drun

private

Definition at line 93 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hvtxyvsbxrun

private

Definition at line 86 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hvtxyvsorbrun

private

Definition at line 79 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hvtxz

private

Definition at line 74 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F** VertexHistogramMaker::m_hvtxzrun

private

Definition at line 77 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH2F** VertexHistogramMaker::m_hvtxzvsbx2drun

private

Definition at line 94 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hvtxzvsbxrun

private

Definition at line 87 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TProfile** VertexHistogramMaker::m_hvtxzvsorbrun

private

Definition at line 80 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

TH1F* VertexHistogramMaker::m_hweights

private

Definition at line 71 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

edm::EDGetTokenT<LumiDetails> VertexHistogramMaker::m_lumiDetailsToken

private

Definition at line 52 of file VertexHistogramMaker.h.

Referenced by fill().

const unsigned int VertexHistogramMaker::m_maxLS

private

Definition at line 42 of file VertexHistogramMaker.h.

Referenced by book().

RunHistogramManager VertexHistogramMaker::m_rhm

private

Definition at line 54 of file VertexHistogramMaker.h.

Referenced by beginRun(), and book().

const bool VertexHistogramMaker::m_runHisto

private

Definition at line 45 of file VertexHistogramMaker.h.

Referenced by beginRun(), book(), and fill().

const bool VertexHistogramMaker::m_runHisto2D

private

Definition at line 49 of file VertexHistogramMaker.h.

Referenced by beginRun(), book(), and fill().

const bool VertexHistogramMaker::m_runHistoBXProfile

private

Definition at line 47 of file VertexHistogramMaker.h.

Referenced by beginRun(), book(), and fill().

const bool VertexHistogramMaker::m_runHistoBXProfile2D

private

Definition at line 48 of file VertexHistogramMaker.h.

Referenced by beginRun(), book(), and fill().

const bool VertexHistogramMaker::m_runHistoProfile

private

Definition at line 46 of file VertexHistogramMaker.h.

Referenced by beginRun(), book(), and fill().

const bool VertexHistogramMaker::m_trueOnly

private

Definition at line 44 of file VertexHistogramMaker.h.

Referenced by fill().

const double VertexHistogramMaker::m_weightThreshold

private

Definition at line 43 of file VertexHistogramMaker.h.

Referenced by book(), and fill().