#include <DigiLumiCorrHistogramMaker.h>

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

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

void	book (const std::string dirname, edm::ConsumesCollector &&iC)

void	book (const std::string dirname, edm::ConsumesCollector &iC)

	DigiLumiCorrHistogramMaker (edm::ConsumesCollector &&iC)

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

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

	~DigiLumiCorrHistogramMaker ()

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

std::map< unsigned int, RunHistogramManager * >	m_fhm

std::string	m_hitname

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

edm::EDGetTokenT< LumiDetails >	m_lumiProducerToken

const double	m_maxlumi

const int	m_nbins

std::map< unsigned int, TH2F * >	m_nmultvslumi

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

std::map< unsigned int, TProfile2D ** >	m_nmultvslumivsbxprofrun

bool	m_runHisto

const int	m_scalefact

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

Detailed Description

Definition at line 21 of file DigiLumiCorrHistogramMaker.h.

Constructor & Destructor Documentation

◆ DigiLumiCorrHistogramMaker() [1/2]

DigiLumiCorrHistogramMaker::DigiLumiCorrHistogramMaker ( edm::ConsumesCollector && iC )

Definition at line 14 of file DigiLumiCorrHistogramMaker.cc.

     : m_lumiProducerToken(iC.consumes<LumiDetails, edm::InLumi>(edm::InputTag("lumiProducer"))),
       m_fhm(),
       m_runHisto(false),
       m_hitname(),
       m_nbins(500),
       m_scalefact(),
       m_maxlumi(10.),
       m_binmax(),
       m_labels(),
       m_nmultvslumi(),
       m_nmultvslumiprof(),
       m_subdirs() {}

◆ DigiLumiCorrHistogramMaker() [2/2]

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

Definition at line 28 of file DigiLumiCorrHistogramMaker.cc.

References edm::ParameterSet::getUntrackedParameter(), m_binmax, m_labels, and AlCaHLTBitMon_QueryRunRegistry::string.

     : m_lumiProducerToken(iC.consumes<LumiDetails, edm::InLumi>(iConfig.getParameter<edm::InputTag>("lumiProducer"))),
       m_fhm(),
       m_runHisto(iConfig.getUntrackedParameter<bool>("runHisto", false)),
       m_hitname(iConfig.getUntrackedParameter<std::string>("hitName", "digi")),
       m_nbins(iConfig.getUntrackedParameter<int>("numberOfBins", 500)),
       m_scalefact(iConfig.getUntrackedParameter<int>("scaleFactor", 5)),
       m_maxlumi(iConfig.getUntrackedParameter<double>("maxLumi", 10.)),
       m_labels(),
       m_nmultvslumi(),
       m_nmultvslumiprof(),
       m_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) {
     m_labels[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<std::string>("detLabel");
     m_binmax[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<int>("binMax");
   }
 }

◆ ~DigiLumiCorrHistogramMaker()

DigiLumiCorrHistogramMaker::~DigiLumiCorrHistogramMaker ( )

Definition at line 49 of file DigiLumiCorrHistogramMaker.cc.

References mps_fire::i, m_fhm, m_labels, m_subdirs, and AlCaHLTBitMon_QueryRunRegistry::string.

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

Member Function Documentation

◆ beginRun()

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

Definition at line 122 of file DigiLumiCorrHistogramMaker.cc.

References mps_fire::i, m_fhm, m_labels, m_nmultvslumivsbxprofrun, m_runHisto, m_subdirs, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MultiplicityInvestigator::beginRun().

                                                             {
   edm::Service<TFileService> tfserv;
 
   for (std::map<unsigned int, std::string>::const_iterator lab = m_labels.begin(); lab != m_labels.end(); ++lab) {
     const int i = lab->first;
     const std::string slab = lab->second;
     m_fhm[i]->beginRun(iRun, *m_subdirs[i]);
     if (m_runHisto) {
       (*m_nmultvslumivsbxprofrun[i])->GetXaxis()->SetTitle("BX");
       (*m_nmultvslumivsbxprofrun[i])->GetYaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
     }
   }
 }

◆ book() [1/3]

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

Definition at line 59 of file DigiLumiCorrHistogramMaker.cc.

References BTVHLTOfflineSource_cfi::dirname, SummaryClient_cfi::labels, and m_labels.

Referenced by MultiplicityInvestigator::MultiplicityInvestigator().

                                                                  {
   m_labels = labels;
   book(dirname, iC);
 }

◆ book() [2/3]

void DigiLumiCorrHistogramMaker::book	(	const std::string	dirname,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 29 of file DigiLumiCorrHistogramMaker.h.

References book(), and BTVHLTOfflineSource_cfi::dirname.

Referenced by book().

29 { book(dirname, iC); }

DigiLumiCorrHistogramMaker::book

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

Definition: DigiLumiCorrHistogramMaker.cc:59

BTVHLTOfflineSource_cfi.dirname

dirname

Definition: BTVHLTOfflineSource_cfi.py:7

◆ book() [3/3]

void DigiLumiCorrHistogramMaker::book	(	const std::string	dirname,
		edm::ConsumesCollector &	iC
	)

Definition at line 66 of file DigiLumiCorrHistogramMaker.cc.

References BTVHLTOfflineSource_cfi::dirname, mps_fire::i, m_binmax, m_fhm, m_hitname, m_labels, m_maxlumi, m_nbins, m_nmultvslumi, m_nmultvslumiprof, m_nmultvslumivsbxprofrun, m_runHisto, m_scalefact, m_subdirs, 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: " << m_nbins;
   edm::LogInfo("ScaleFactors") << "y-axis range scale factor: " << m_scalefact;
   edm::LogInfo("MaxLumi") << "max lumi value: " << m_maxlumi;
   edm::LogInfo("BinMaxValue") << "Setting bin max values";
 
   for (std::map<unsigned int, std::string>::const_iterator lab = m_labels.begin(); lab != m_labels.end(); lab++) {
     const unsigned int i = lab->first;
     const std::string slab = lab->second;
 
     if (m_binmax.find(i) == m_binmax.end()) {
       edm::LogVerbatim("NotConfiguredBinMax")
           << "Bin max for " << lab->second << " not configured: " << trnumb.nstrips(i) << " used";
       m_binmax[i] = trnumb.nstrips(i);
     }
 
     edm::LogVerbatim("BinMaxValue") << "Bin max for " << lab->second << " is " << m_binmax[i];
   }
 
   for (std::map<unsigned int, std::string>::const_iterator lab = m_labels.begin(); lab != m_labels.end(); ++lab) {
     const int i = lab->first;
     const std::string slab = lab->second;
 
     char name[200];
     char title[500];
 
     m_subdirs[i] = new TFileDirectory(subev.mkdir(slab));
     m_fhm[i] = new RunHistogramManager(iC, true);
 
     if (m_subdirs[i]) {
       sprintf(name, "n%sdigivslumi", slab.c_str());
       sprintf(title, "%s %s multiplicity vs BX lumi", slab.c_str(), m_hitname.c_str());
       m_nmultvslumi[i] = m_subdirs[i]->make<TH2F>(
           name, title, 250, 0., m_maxlumi, m_nbins, 0., (1 + m_binmax[i] / (m_scalefact * m_nbins)) * m_nbins);
       m_nmultvslumi[i]->GetXaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
       m_nmultvslumi[i]->GetYaxis()->SetTitle("Number of Hits");
       sprintf(name, "n%sdigivslumiprof", slab.c_str());
       m_nmultvslumiprof[i] = m_subdirs[i]->make<TProfile>(name, title, 250, 0., m_maxlumi);
       m_nmultvslumiprof[i]->GetXaxis()->SetTitle("BX lumi [10^{30}cm^{-2}s^{-1}]");
       m_nmultvslumiprof[i]->GetYaxis()->SetTitle("Number of Hits");
 
       if (m_runHisto) {
         edm::LogInfo("RunHistos") << "Pseudo-booking run histos " << slab.c_str();
         sprintf(name, "n%sdigivslumivsbxprofrun", slab.c_str());
         sprintf(title, "%s %s multiplicity vs BX lumi vs BX", slab.c_str(), m_hitname.c_str());
         m_nmultvslumivsbxprofrun[i] = m_fhm[i]->makeTProfile2D(name, title, 3564, -0.5, 3563.5, 250, 0., m_maxlumi);
       }
     }
   }
 }

◆ fill()

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

Definition at line 136 of file DigiLumiCorrHistogramMaker.cc.

References HcalObjRepresent::Fill(), mps_fire::i, iEvent, LumiDetails::isValid(), edm::HandleBase::isValid(), LumiDetails::kOCC1, LumiDetails::lumiValue(), m_labels, m_lumiProducerToken, m_nmultvslumi, m_nmultvslumiprof, and m_nmultvslumivsbxprofrun.

Referenced by MultiplicityInvestigator::analyze().

                                                                                                     {
   edm::Handle<LumiDetails> ld;
   iEvent.getLuminosityBlock().getByToken(m_lumiProducerToken, ld);
 
   if (ld.isValid()) {
     if (ld->isValid()) {
       float bxlumi = ld->lumiValue(LumiDetails::kOCC1, iEvent.bunchCrossing()) * 6.37;
 
       for (std::map<unsigned int, int>::const_iterator digi = ndigi.begin(); digi != ndigi.end(); digi++) {
         if (m_labels.find(digi->first) != m_labels.end()) {
           const unsigned int i = digi->first;
           m_nmultvslumi[i]->Fill(bxlumi, digi->second);
           m_nmultvslumiprof[i]->Fill(bxlumi, digi->second);
 
           if (m_nmultvslumivsbxprofrun[i] && *m_nmultvslumivsbxprofrun[i])
             (*m_nmultvslumivsbxprofrun[i])->Fill(iEvent.bunchCrossing() % 3564, bxlumi, digi->second);
         }
       }
     }
   }
 }