#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 (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC)

	VertexHistogramMaker (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() [1/2]

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() [2/2]

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::~VertexHistogramMaker ( )

Definition at line 46 of file VertexHistogramMaker.cc.

46 { delete m_currdir; }

References m_currdir.

Member Function Documentation

◆ beginRun()

void VertexHistogramMaker::beginRun ( const edm::Run & iRun )

Definition at line 257 of file VertexHistogramMaker.cc.

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

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().

◆ book()

void VertexHistogramMaker::book ( const std::string dirname = "" )

Definition at line 48 of file VertexHistogramMaker.cc.

                                                        {
   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.));
       }
     }
   }
 }

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().

◆ fill() [1/2]

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

Definition at line 440 of file VertexHistogramMaker.cc.

                                                                                                                  {
   // 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);
 }

References iEvent, LumiDetails::isValid(), edm::HandleBase::isValid(), LumiDetails::kOCC1, LumiDetails::lumiValue(), m_lumiDetailsToken, pwdgSkimBPark_cfi::vertices, and mps_merge::weight.

Referenced by AnotherPrimaryVertexAnalyzer::analyze().

◆ fill() [2/2]

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.

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

References l1GtPatternGenerator_cfi::bx, 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(), pwdgSkimBPark_cfi::vertices, extraflags_cff::vtx, and mps_merge::weight.

Member Data Documentation

◆ m_bsConstrained

const bool VertexHistogramMaker::m_bsConstrained

private

Definition at line 50 of file VertexHistogramMaker.h.

Referenced by fill().

◆ m_currdir

TFileDirectory* VertexHistogramMaker::m_currdir

private

Definition at line 41 of file VertexHistogramMaker.h.

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

◆ m_fhm

RunHistogramManager VertexHistogramMaker::m_fhm

private

Definition at line 55 of file VertexHistogramMaker.h.

Referenced by beginRun(), and book().

◆ m_haveweight

TH1F* VertexHistogramMaker::m_haveweight

private

Definition at line 65 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_haveweightvsvtxz

TProfile* VertexHistogramMaker::m_haveweightvsvtxz

private

Definition at line 69 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_haveweightvsvtxzchk

TProfile* VertexHistogramMaker::m_haveweightvsvtxzchk

private

Definition at line 70 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_histoParameters

const edm::ParameterSet VertexHistogramMaker::m_histoParameters

private

Definition at line 51 of file VertexHistogramMaker.h.

Referenced by book().

◆ m_hndof

TH1F* VertexHistogramMaker::m_hndof

private

Definition at line 64 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hndofvstracks

TH2F* VertexHistogramMaker::m_hndofvstracks

private

Definition at line 66 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hndofvsvtxz

TProfile* VertexHistogramMaker::m_hndofvsvtxz

private

Definition at line 67 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnheavytracks

TH1F* VertexHistogramMaker::m_hnheavytracks

private

Definition at line 63 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hntracks

TH1F* VertexHistogramMaker::m_hntracks

private

Definition at line 60 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hntracksvsvtxz

TProfile* VertexHistogramMaker::m_hntracksvsvtxz

private

Definition at line 68 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hntruevtx

TH1F* VertexHistogramMaker::m_hntruevtx

private

Definition at line 57 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hntruevtxvslumi

TProfile* VertexHistogramMaker::m_hntruevtxvslumi

private

Definition at line 58 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hntruevtxvslumi2D

TH2D* VertexHistogramMaker::m_hntruevtxvslumi2D

private

Definition at line 59 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnvtx

TH1F* VertexHistogramMaker::m_hnvtx

private

Definition at line 56 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnvtxvsbxrun

TProfile** VertexHistogramMaker::m_hnvtxvsbxrun

private

Definition at line 88 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnvtxvsbxvslumirun

TProfile2D** VertexHistogramMaker::m_hnvtxvsbxvslumirun

private

Definition at line 90 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnvtxvsbxvsorbrun

TProfile2D** VertexHistogramMaker::m_hnvtxvsbxvsorbrun

private

Definition at line 82 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnvtxvsorbrun

TProfile** VertexHistogramMaker::m_hnvtxvsorbrun

private

Definition at line 81 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hnvtxvsorbrun2D

TH2F** VertexHistogramMaker::m_hnvtxvsorbrun2D

private

Definition at line 83 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hsqsumptsq

TH1F* VertexHistogramMaker::m_hsqsumptsq

private

Definition at line 61 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hsqsumptsqheavy

TH1F* VertexHistogramMaker::m_hsqsumptsqheavy

private

Definition at line 62 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxx

TH1F* VertexHistogramMaker::m_hvtxx

private

Definition at line 72 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxxrun

TH1F** VertexHistogramMaker::m_hvtxxrun

private

Definition at line 75 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxxvsbx2drun

TH2F** VertexHistogramMaker::m_hvtxxvsbx2drun

private

Definition at line 92 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxxvsbxrun

TProfile** VertexHistogramMaker::m_hvtxxvsbxrun

private

Definition at line 85 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxxvsorbrun

TProfile** VertexHistogramMaker::m_hvtxxvsorbrun

private

Definition at line 78 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxy

TH1F* VertexHistogramMaker::m_hvtxy

private

Definition at line 73 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxyrun

TH1F** VertexHistogramMaker::m_hvtxyrun

private

Definition at line 76 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxyvsbx2drun

TH2F** VertexHistogramMaker::m_hvtxyvsbx2drun

private

Definition at line 93 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxyvsbxrun

TProfile** VertexHistogramMaker::m_hvtxyvsbxrun

private

Definition at line 86 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxyvsorbrun

TProfile** VertexHistogramMaker::m_hvtxyvsorbrun

private

Definition at line 79 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxz

TH1F* VertexHistogramMaker::m_hvtxz

private

Definition at line 74 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxzrun

TH1F** VertexHistogramMaker::m_hvtxzrun

private

Definition at line 77 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxzvsbx2drun

TH2F** VertexHistogramMaker::m_hvtxzvsbx2drun

private

Definition at line 94 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxzvsbxrun

TProfile** VertexHistogramMaker::m_hvtxzvsbxrun

private

Definition at line 87 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hvtxzvsorbrun

TProfile** VertexHistogramMaker::m_hvtxzvsorbrun

private

Definition at line 80 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_hweights

TH1F* VertexHistogramMaker::m_hweights

private

Definition at line 71 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

◆ m_lumiDetailsToken

edm::EDGetTokenT<LumiDetails> VertexHistogramMaker::m_lumiDetailsToken

private

Definition at line 52 of file VertexHistogramMaker.h.

Referenced by fill().

◆ m_maxLS

const unsigned int VertexHistogramMaker::m_maxLS

private

Definition at line 42 of file VertexHistogramMaker.h.

Referenced by book().

◆ m_rhm

RunHistogramManager VertexHistogramMaker::m_rhm

private

Definition at line 54 of file VertexHistogramMaker.h.

Referenced by beginRun(), and book().

◆ m_runHisto

const bool VertexHistogramMaker::m_runHisto

private

Definition at line 45 of file VertexHistogramMaker.h.

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

◆ m_runHisto2D

const bool VertexHistogramMaker::m_runHisto2D

private

Definition at line 49 of file VertexHistogramMaker.h.

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

◆ m_runHistoBXProfile

const bool VertexHistogramMaker::m_runHistoBXProfile

private

Definition at line 47 of file VertexHistogramMaker.h.

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

◆ m_runHistoBXProfile2D

const bool VertexHistogramMaker::m_runHistoBXProfile2D

private

Definition at line 48 of file VertexHistogramMaker.h.

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

◆ m_runHistoProfile

const bool VertexHistogramMaker::m_runHistoProfile

private

Definition at line 46 of file VertexHistogramMaker.h.

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

◆ m_trueOnly

const bool VertexHistogramMaker::m_trueOnly

private

Definition at line 44 of file VertexHistogramMaker.h.

Referenced by fill().

◆ m_weightThreshold

const double VertexHistogramMaker::m_weightThreshold

private

Definition at line 43 of file VertexHistogramMaker.h.

Referenced by book(), and fill().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ VertexHistogramMaker() [1/2]

◆ VertexHistogramMaker() [2/2]

◆ ~VertexHistogramMaker()

Member Function Documentation

◆ beginRun()

◆ book()

◆ fill() [1/2]

◆ fill() [2/2]

Member Data Documentation

◆ m_bsConstrained

◆ 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_lumiDetailsToken

◆ m_maxLS

◆ m_rhm

◆ m_runHisto

◆ m_runHisto2D

◆ m_runHistoBXProfile

◆ m_runHistoBXProfile2D

◆ m_runHistoProfile

◆ m_trueOnly

◆ m_weightThreshold