BeamSpotHistogramMaker Class Reference

#include <BeamSpotHistogramMaker.h>

Public Member Functions
	BeamSpotHistogramMaker (edm::ConsumesCollector &&iC)
	BeamSpotHistogramMaker (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC)
void	beginRun (const unsigned int nrun)
void	book (const std::string dirname="")
void	fill (const unsigned int orbit, const reco::BeamSpot &bs)
	~BeamSpotHistogramMaker ()

Private Attributes
TFileDirectory *	_currdir
TH1F **	_hbssigmaxrun
TProfile **	_hbssigmaxvsorbrun
TH1F **	_hbssigmayrun
TProfile **	_hbssigmayvsorbrun
TH1F **	_hbssigmazrun
TProfile **	_hbssigmazvsorbrun
TH1F **	_hbsxrun
TProfile **	_hbsxvsorbrun
TH1F **	_hbsyrun
TProfile **	_hbsyvsorbrun
TH1F **	_hbszrun
TProfile **	_hbszvsorbrun
const edm::ParameterSet	_histoParameters
RunHistogramManager	_rhm

Detailed Description

Definition at line 21 of file BeamSpotHistogramMaker.h.

Constructor & Destructor Documentation

BeamSpotHistogramMaker::BeamSpotHistogramMaker

(

edm::ConsumesCollector &&

iC

)

Definition at line 11 of file BeamSpotHistogramMaker.cc.

    : _currdir(nullptr), _histoParameters(), _rhm(iC) {}

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

Definition at line 14 of file BeamSpotHistogramMaker.cc.

    : _currdir(nullptr),
      _histoParameters(iConfig.getUntrackedParameter<edm::ParameterSet>("histoParameters", edm::ParameterSet())),
      _rhm(iC) {}

BeamSpotHistogramMaker::~BeamSpotHistogramMaker

(

)

Definition at line 19 of file BeamSpotHistogramMaker.cc.

References _currdir.

{ delete _currdir; }

Member Function Documentation

void BeamSpotHistogramMaker::beginRun

(

const unsigned int

nrun

)

Definition at line 82 of file BeamSpotHistogramMaker.cc.

References _currdir, _rhm, RunHistogramManager::beginRun(), and TFileService::tFileDirectory().

Referenced by AnotherBeamSpotAnalyzer::beginRun().

                                                             {
  char runname[100];
  sprintf(runname, "run_%d", nrun);

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

  _rhm.beginRun(nrun, *currdir);

  (*_hbsxrun)->GetXaxis()->SetTitle("X [cm]");
  (*_hbsxrun)->GetYaxis()->SetTitle("Events");
  (*_hbsyrun)->GetXaxis()->SetTitle("Y [cm]");
  (*_hbsyrun)->GetYaxis()->SetTitle("Events");
  (*_hbszrun)->GetXaxis()->SetTitle("Z [cm]");
  (*_hbszrun)->GetYaxis()->SetTitle("Events");
  (*_hbssigmaxrun)->GetXaxis()->SetTitle("sigmaX [cm]");
  (*_hbssigmaxrun)->GetYaxis()->SetTitle("Events");
  (*_hbssigmayrun)->GetXaxis()->SetTitle("sigmaY [cm]");
  (*_hbssigmayrun)->GetYaxis()->SetTitle("Events");
  (*_hbssigmazrun)->GetXaxis()->SetTitle("sigmaZ [cm]");
  (*_hbssigmazrun)->GetYaxis()->SetTitle("Events");

  (*_hbsxvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
  (*_hbsxvsorbrun)->GetYaxis()->SetTitle("X [cm]");
  (*_hbsxvsorbrun)->SetCanExtend(TH1::kAllAxes);
  (*_hbsyvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
  (*_hbsyvsorbrun)->GetYaxis()->SetTitle("Y [cm]");
  (*_hbsyvsorbrun)->SetCanExtend(TH1::kAllAxes);
  (*_hbszvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
  (*_hbszvsorbrun)->GetYaxis()->SetTitle("Z [cm]");
  (*_hbszvsorbrun)->SetCanExtend(TH1::kAllAxes);
  (*_hbssigmaxvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
  (*_hbssigmaxvsorbrun)->GetYaxis()->SetTitle("sigmaX [cm]");
  (*_hbssigmaxvsorbrun)->SetCanExtend(TH1::kAllAxes);
  (*_hbssigmayvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
  (*_hbssigmayvsorbrun)->GetYaxis()->SetTitle("sigmaY [cm]");
  (*_hbssigmayvsorbrun)->SetCanExtend(TH1::kAllAxes);
  (*_hbssigmazvsorbrun)->GetXaxis()->SetTitle("time [orbit#]");
  (*_hbssigmazvsorbrun)->GetYaxis()->SetTitle("sigmaZ [cm]");
  (*_hbssigmazvsorbrun)->SetCanExtend(TH1::kAllAxes);
}

void BeamSpotHistogramMaker::book

(

const std::string

dirname = ""

)

Definition at line 21 of file BeamSpotHistogramMaker.cc.

References _currdir, _hbssigmaxrun, _hbssigmaxvsorbrun, _hbssigmayrun, _hbssigmayvsorbrun, _hbssigmazrun, _hbssigmazvsorbrun, _hbsxrun, _hbsxvsorbrun, _hbsyrun, _hbsyvsorbrun, _hbszrun, _hbszvsorbrun, _histoParameters, _rhm, BTVHLTOfflineSource_cfi::dirname, edm::ParameterSet::getUntrackedParameter(), RunHistogramManager::makeTH1F(), RunHistogramManager::makeTProfile(), TFileService::mkdir(), and TFileService::tFileDirectory().

Referenced by AnotherBeamSpotAnalyzer::AnotherBeamSpotAnalyzer().

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

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

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

  _hbsxrun = _rhm.makeTH1F("bsxrun",
                           "BeamSpot X position",
                           _histoParameters.getUntrackedParameter<unsigned int>("nBinX", 200),
                           _histoParameters.getUntrackedParameter<double>("xMin", -1.),
                           _histoParameters.getUntrackedParameter<double>("xMax", 1.));

  _hbsyrun = _rhm.makeTH1F("bsyrun",
                           "BeamSpot Y position",
                           _histoParameters.getUntrackedParameter<unsigned int>("nBinY", 200),
                           _histoParameters.getUntrackedParameter<double>("yMin", -1.),
                           _histoParameters.getUntrackedParameter<double>("yMax", 1.));

  _hbszrun = _rhm.makeTH1F("bszrun",
                           "BeamSpot Z position",
                           _histoParameters.getUntrackedParameter<unsigned int>("nBinZ", 200),
                           _histoParameters.getUntrackedParameter<double>("zMin", -1.),
                           _histoParameters.getUntrackedParameter<double>("zMax", 1.));

  _hbssigmaxrun = _rhm.makeTH1F("bssigmaxrun",
                                "BeamSpot sigmaX",
                                _histoParameters.getUntrackedParameter<unsigned int>("nBinSigmaX", 200),
                                _histoParameters.getUntrackedParameter<double>("sigmaXMin", 0.),
                                _histoParameters.getUntrackedParameter<double>("sigmaXMax", 0.025));

  _hbssigmayrun = _rhm.makeTH1F("bssigmayrun",
                                "BeamSpot sigmaY",
                                _histoParameters.getUntrackedParameter<unsigned int>("nBinSigmaY", 200),
                                _histoParameters.getUntrackedParameter<double>("sigmaYMin", 0.),
                                _histoParameters.getUntrackedParameter<double>("sigmaYMax", 0.025));

  _hbssigmazrun = _rhm.makeTH1F("bssigmazrun",
                                "BeamSpot sigmaZ",
                                _histoParameters.getUntrackedParameter<unsigned int>("nBinSigmaZ", 200),
                                _histoParameters.getUntrackedParameter<double>("sigmaZMin", 0.),
                                _histoParameters.getUntrackedParameter<double>("sigmaZMax", 15.));

  _hbsxvsorbrun =
      _rhm.makeTProfile("bsxvsorbrun", "BeamSpot X position vs orbit number", 1600, 0.5, 1600. * 16384 + 0.5);
  _hbsyvsorbrun =
      _rhm.makeTProfile("bsyvsorbrun", "BeamSpot Y position vs orbit number", 1600, 0.5, 1600. * 16384 + 0.5);
  _hbszvsorbrun =
      _rhm.makeTProfile("bszvsorbrun", "BeamSpot Z position vs orbit number", 1600, 0.5, 1600. * 16384 + 0.5);
  _hbssigmaxvsorbrun =
      _rhm.makeTProfile("bssigmaxvsorbrun", "BeamSpot sigmaX vs orbit number", 1600, 0.5, 1600. * 16384 + 0.5);
  _hbssigmayvsorbrun =
      _rhm.makeTProfile("bssigmayvsorbrun", "BeamSpot sigmaY vs orbit number", 1600, 0.5, 1600. * 16384 + 0.5);
  _hbssigmazvsorbrun =
      _rhm.makeTProfile("bssigmazvsorbrun", "BeamSpot sigmaZ vs orbit number", 1600, 0.5, 1600. * 16384 + 0.5);
}

void BeamSpotHistogramMaker::fill	(	const unsigned int	orbit,
		const reco::BeamSpot &	bs
	)

Definition at line 127 of file BeamSpotHistogramMaker.cc.

References _hbssigmaxrun, _hbssigmaxvsorbrun, _hbssigmayrun, _hbssigmayvsorbrun, _hbssigmazrun, _hbssigmazvsorbrun, _hbsxrun, _hbsxvsorbrun, _hbsyrun, _hbsyvsorbrun, _hbszrun, _hbszvsorbrun, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::sigmaZ(), reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by AnotherBeamSpotAnalyzer::analyze().

                                                                                  {
  if (_hbsxrun && *_hbsxrun)
    (*_hbsxrun)->Fill(bs.x0());
  if (_hbsxvsorbrun && *_hbsxvsorbrun)
    (*_hbsxvsorbrun)->Fill(orbit, bs.x0());

  if (_hbsyrun && *_hbsyrun)
    (*_hbsyrun)->Fill(bs.y0());
  if (_hbsyvsorbrun && *_hbsyvsorbrun)
    (*_hbsyvsorbrun)->Fill(orbit, bs.y0());

  if (_hbszrun && *_hbszrun)
    (*_hbszrun)->Fill(bs.z0());
  if (_hbszvsorbrun && *_hbszvsorbrun)
    (*_hbszvsorbrun)->Fill(orbit, bs.z0());

  if (_hbssigmaxrun && *_hbssigmaxrun)
    (*_hbssigmaxrun)->Fill(bs.BeamWidthX());
  if (_hbssigmayrun && *_hbssigmayrun)
    (*_hbssigmayrun)->Fill(bs.BeamWidthY());
  if (_hbssigmazrun && *_hbssigmazrun)
    (*_hbssigmazrun)->Fill(bs.sigmaZ());
  if (_hbssigmaxvsorbrun && *_hbssigmaxvsorbrun)
    (*_hbssigmaxvsorbrun)->Fill(orbit, bs.BeamWidthX());
  if (_hbssigmayvsorbrun && *_hbssigmayvsorbrun)
    (*_hbssigmayvsorbrun)->Fill(orbit, bs.BeamWidthY());
  if (_hbssigmazvsorbrun && *_hbssigmazvsorbrun)
    (*_hbssigmazvsorbrun)->Fill(orbit, bs.sigmaZ());
}

Member Data Documentation

TFileDirectory* BeamSpotHistogramMaker::_currdir

private

Definition at line 33 of file BeamSpotHistogramMaker.h.

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

TH1F** BeamSpotHistogramMaker::_hbssigmaxrun

private

Definition at line 40 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TProfile** BeamSpotHistogramMaker::_hbssigmaxvsorbrun

private

Definition at line 46 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TH1F** BeamSpotHistogramMaker::_hbssigmayrun

private

Definition at line 41 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TProfile** BeamSpotHistogramMaker::_hbssigmayvsorbrun

private

Definition at line 47 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TH1F** BeamSpotHistogramMaker::_hbssigmazrun

private

Definition at line 42 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TProfile** BeamSpotHistogramMaker::_hbssigmazvsorbrun

private

Definition at line 48 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TH1F** BeamSpotHistogramMaker::_hbsxrun

private

Definition at line 37 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TProfile** BeamSpotHistogramMaker::_hbsxvsorbrun

private

Definition at line 43 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TH1F** BeamSpotHistogramMaker::_hbsyrun

private

Definition at line 38 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TProfile** BeamSpotHistogramMaker::_hbsyvsorbrun

private

Definition at line 44 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TH1F** BeamSpotHistogramMaker::_hbszrun

private

Definition at line 39 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

TProfile** BeamSpotHistogramMaker::_hbszvsorbrun

private

Definition at line 45 of file BeamSpotHistogramMaker.h.

Referenced by book(), and fill().

const edm::ParameterSet BeamSpotHistogramMaker::_histoParameters

private

Definition at line 34 of file BeamSpotHistogramMaker.h.

Referenced by book().

RunHistogramManager BeamSpotHistogramMaker::_rhm

private

Definition at line 36 of file BeamSpotHistogramMaker.h.

Referenced by beginRun(), and book().