#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 13 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 17 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 40 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 123 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 53 of file DigiInvestigatorHistogramMaker.cc.

References _labels, and tablePrinter::labels.

Referenced by MultiplicityInvestigator::MultiplicityInvestigator().

                                                                                                                 {
 
   _labels = labels;
   book(dirname);
 
 }

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

Definition at line 60 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(), dataset::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 168 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);
       }
     }
 
   }
 }