#include <DigiInvestigatorHistogramMaker.h>

Public Member Functions
void	beginRun (const edm::Run &iRun)

void	book (const std::string dirname, const std::map< unsigned int, std::string > &labels)

void	book (const std::string dirname)

	DigiInvestigatorHistogramMaker (edm::ConsumesCollector &&iC)

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

void	fill (const edm::Event &iEvent, const std::map< unsigned int, int > &ndigi)

	~DigiInvestigatorHistogramMaker ()

Private Attributes
std::map< unsigned int, int >	_binmax

RunHistogramManager	_fhm

const bool	_fillHisto

std::string	_hitname

std::map< unsigned int, std::string >	_labels

const int	_nbins

std::map< unsigned int, TH1F * >	_nmult

std::map< unsigned int, TProfile ** >	_nmultvsbxfill

std::map< unsigned int, TProfile ** >	_nmultvsbxrun

std::map< unsigned int, TProfile ** >	_nmultvsorbrun

RunHistogramManager	_rhm

const bool	_runHisto

int	_scalefact

std::map< unsigned int, TFileDirectory * >	_subdirs

const unsigned int	m_LSfrac

const unsigned int	m_maxLS

Detailed Description

Definition at line 18 of file DigiInvestigatorHistogramMaker.h.

Constructor & Destructor Documentation

DigiInvestigatorHistogramMaker::DigiInvestigatorHistogramMaker ( edm::ConsumesCollector && iC )

Definition at line 12 of file DigiInvestigatorHistogramMaker.cc.

     : _hitname(),
       _nbins(500),
       m_maxLS(100),
       m_LSfrac(4),
       _scalefact(),
       _runHisto(true),
       _fillHisto(false),
       _binmax(),
       _labels(),
       _rhm(iC),
       _fhm(iC, true),
       _nmultvsorbrun(),
       _nmultvsbxrun(),
       _nmultvsbxfill(),
       _nmult() {}

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

Definition at line 29 of file DigiInvestigatorHistogramMaker.cc.

References _binmax, _labels, edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : _hitname(iConfig.getUntrackedParameter<std::string>("hitName", "digi")),
       _nbins(iConfig.getUntrackedParameter<int>("numberOfBins", 500)),
       m_maxLS(iConfig.getUntrackedParameter<unsigned int>("maxLSBeforeRebin", 100)),
       m_LSfrac(iConfig.getUntrackedParameter<unsigned int>("startingLSFraction", 4)),
       _scalefact(iConfig.getUntrackedParameter<int>("scaleFactor", 5)),
       _runHisto(iConfig.getUntrackedParameter<bool>("runHisto", true)),
       _fillHisto(iConfig.getUntrackedParameter<bool>("fillHisto", false)),
       _labels(),
       _rhm(iC),
       _fhm(iC, true),
       _nmultvsorbrun(),
       _nmultvsbxrun(),
       _nmultvsbxfill(),
       _nmult(),
       _subdirs() {
   std::vector<edm::ParameterSet> wantedsubds(iConfig.getUntrackedParameter<std::vector<edm::ParameterSet> >(
       "wantedSubDets", std::vector<edm::ParameterSet>()));
 
   for (std::vector<edm::ParameterSet>::iterator ps = wantedsubds.begin(); ps != wantedsubds.end(); ++ps) {
     _labels[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<std::string>("detLabel");
     _binmax[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<int>("binMax");
   }
 }

DigiInvestigatorHistogramMaker::~DigiInvestigatorHistogramMaker ( )

Definition at line 55 of file DigiInvestigatorHistogramMaker.cc.

References _labels, _subdirs, mps_fire::i, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                 {
   for (std::map<unsigned int, std::string>::const_iterator lab = _labels.begin(); lab != _labels.end(); lab++) {
     const unsigned int i = lab->first;
     const std::string slab = lab->second;
 
     delete _subdirs[i];
   }
 }

Member Function Documentation

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

Definition at line 128 of file DigiInvestigatorHistogramMaker.cc.

References _fhm, _fillHisto, _labels, _nmultvsbxfill, _nmultvsbxrun, _nmultvsorbrun, _rhm, _runHisto, RunHistogramManager::beginRun(), mps_fire::i, AlCaHLTBitMon_QueryRunRegistry::string, and TFileService::tFileDirectory().

Referenced by MultiplicityInvestigator::beginRun().

                                                                 {
   //  char runname[100];
   //  sprintf(runname,"run_%d",nrun);
 
   edm::Service<TFileService> tfserv;
 
   //  currdir = &(*tfserv);
   //  _rhm.beginRun(nrun,*currdir);
 
   _rhm.beginRun(iRun, tfserv->tFileDirectory());
   _fhm.beginRun(iRun, tfserv->tFileDirectory());
 
   for (std::map<unsigned int, std::string>::const_iterator lab = _labels.begin(); lab != _labels.end(); ++lab) {
     const int i = lab->first;
     const std::string slab = lab->second;
 
     //    char name[200];
     //    char title[500];
 
     //    TFileDirectory subd =_subdirs[i]->mkdir(runname);
 
     //    sprintf(name,"n%sdigivsorbrun",slab.c_str());
     //    sprintf(title,"%s %s mean multiplicity vs orbit",slab.c_str(),_hitname.c_str());
     //    _nmultvsorbrun[i] = subd.make<TProfile>(name,title,_norbbin,0.5,11223*_norbbin+0.5);
     if (_runHisto) {
       if (*_nmultvsorbrun[i]) {
         (*_nmultvsorbrun[i])->GetXaxis()->SetTitle("time [orbit#]");
         (*_nmultvsorbrun[i])->GetYaxis()->SetTitle("Hits");
         (*_nmultvsorbrun[i])->SetCanExtend(TH1::kXaxis);
       }
       if (*_nmultvsbxrun[i]) {
         (*_nmultvsbxrun[i])->GetXaxis()->SetTitle("BX#");
         (*_nmultvsbxrun[i])->GetYaxis()->SetTitle("Mean Number of Hits");
       }
     }
     if (_fillHisto) {
       if (*_nmultvsbxfill[i]) {
         (*_nmultvsbxfill[i])->GetXaxis()->SetTitle("BX#");
         (*_nmultvsbxfill[i])->GetYaxis()->SetTitle("Mean Number of Hits");
       }
     }
   }
 }

void DigiInvestigatorHistogramMaker::book	(	const std::string	dirname,
		const std::map< unsigned int, std::string > &	labels
	)

Definition at line 64 of file DigiInvestigatorHistogramMaker.cc.

References _labels, and HLT_FULL_cff::labels.

Referenced by MultiplicityInvestigator::MultiplicityInvestigator().

                                                                                          {
   _labels = labels;
   book(dirname);
 }

void DigiInvestigatorHistogramMaker::book ( const std::string dirname )

Definition at line 70 of file DigiInvestigatorHistogramMaker.cc.

References _binmax, _fhm, _fillHisto, _hitname, _labels, _nbins, _nmult, _nmultvsbxfill, _nmultvsbxrun, _nmultvsorbrun, _rhm, _runHisto, _scalefact, _subdirs, mps_fire::i, m_LSfrac, m_maxLS, RunHistogramManager::makeTProfile(), TFileDirectory::mkdir(), TFileService::mkdir(), mergeVDriftHistosByStation::name, SiStripTKNumbers::nstrips(), AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

                                                                  {
   edm::Service<TFileService> tfserv;
   TFileDirectory subev = tfserv->mkdir(dirname);
 
   SiStripTKNumbers trnumb;
 
   edm::LogInfo("NumberOfBins") << "Number of Bins: " << _nbins;
   edm::LogInfo("NumberOfMaxLS") << "Max number of LS before rebinning: " << m_maxLS;
   edm::LogInfo("StartingLSFrac") << "Fraction of LS in one bin before rebinning: " << m_LSfrac;
   edm::LogInfo("ScaleFactors") << "x-axis range scale factor: " << _scalefact;
   edm::LogInfo("BinMaxValue") << "Setting bin max values";
 
   for (std::map<unsigned int, std::string>::const_iterator lab = _labels.begin(); lab != _labels.end(); lab++) {
     const unsigned int i = lab->first;
     const std::string slab = lab->second;
 
     if (_binmax.find(i) == _binmax.end()) {
       edm::LogVerbatim("NotConfiguredBinMax")
           << "Bin max for " << lab->second << " not configured: " << trnumb.nstrips(i) << " used";
       _binmax[i] = trnumb.nstrips(i);
     }
 
     edm::LogVerbatim("BinMaxValue") << "Bin max for " << lab->second << " is " << _binmax[i];
   }
 
   for (std::map<unsigned int, std::string>::const_iterator lab = _labels.begin(); lab != _labels.end(); ++lab) {
     const int i = lab->first;
     const std::string slab = lab->second;
 
     char name[200];
     char title[500];
 
     _subdirs[i] = new TFileDirectory(subev.mkdir(slab));
 
     if (_subdirs[i]) {
       sprintf(name, "n%sdigi", slab.c_str());
       sprintf(title, "%s %s multiplicity", slab.c_str(), _hitname.c_str());
       _nmult[i] = _subdirs[i]->make<TH1F>(name, title, _nbins, 0., (1 + _binmax[i] / (_scalefact * _nbins)) * _nbins);
       _nmult[i]->GetXaxis()->SetTitle("Number of Hits");
       _nmult[i]->GetYaxis()->SetTitle("Events");
 
       if (_runHisto) {
         sprintf(name, "n%sdigivsorbrun", slab.c_str());
         sprintf(title, "%s %s mean multiplicity vs orbit", slab.c_str(), _hitname.c_str());
         _nmultvsorbrun[i] = _rhm.makeTProfile(name, title, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
         sprintf(name, "n%sdigivsbxrun", slab.c_str());
         sprintf(title, "%s %s mean multiplicity vs BX", slab.c_str(), _hitname.c_str());
         _nmultvsbxrun[i] = _rhm.makeTProfile(name, title, 3564, -0.5, 3563.5);
       }
       if (_fillHisto) {
         sprintf(name, "n%sdigivsbxfill", slab.c_str());
         sprintf(title, "%s %s mean multiplicity vs BX", slab.c_str(), _hitname.c_str());
         _nmultvsbxfill[i] = _fhm.makeTProfile(name, title, 3564, -0.5, 3563.5);
       }
     }
   }
 }

void DigiInvestigatorHistogramMaker::fill	(	const edm::Event &	iEvent,
		const std::map< unsigned int, int > &	ndigi
	)

Definition at line 172 of file DigiInvestigatorHistogramMaker.cc.

References _fillHisto, _labels, _nmult, _nmultvsbxfill, _nmultvsbxrun, _nmultvsorbrun, _runHisto, edm::EventBase::bunchCrossing(), HcalObjRepresent::Fill(), mps_fire::i, and edm::EventBase::orbitNumber().

Referenced by MultiplicityInvestigator::analyze().

                                                                                                         {
   for (std::map<unsigned int, int>::const_iterator digi = ndigi.begin(); digi != ndigi.end(); digi++) {
     if (_labels.find(digi->first) != _labels.end()) {
       const unsigned int i = digi->first;
 
       _nmult[i]->Fill(digi->second);
       if (_runHisto) {
         if (_nmultvsorbrun[i] && *_nmultvsorbrun[i])
           (*_nmultvsorbrun[i])->Fill(iEvent.orbitNumber(), digi->second);
         if (_nmultvsbxrun[i] && *_nmultvsbxrun[i])
           (*_nmultvsbxrun[i])->Fill(iEvent.bunchCrossing() % 3564, digi->second);
       }
       if (_fillHisto) {
         if (_nmultvsbxfill[i] && *_nmultvsbxfill[i])
           (*_nmultvsbxfill[i])->Fill(iEvent.bunchCrossing() % 3564, digi->second);
       }
     }
   }
 }