#include <DigiPileupCorrHistogramMaker.h>

Public Member Functions
void	beginRun (const unsigned int nrun)

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

void	book (const std::string dirname)

	DigiPileupCorrHistogramMaker (edm::ConsumesCollector &&iC)

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

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

	~DigiPileupCorrHistogramMaker ()

Private Attributes
const bool	m_2dhisto

std::map< unsigned int, int >	m_binmax

std::string	m_hitname

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

const int	m_nbins

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

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

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

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

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

const int	m_ootBX

edm::EDGetTokenT< std::vector< PileupSummaryInfo > >	m_pileupcollectionToken

const int	m_scalefact

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

bool	m_useVisibleVertices

Detailed Description

Definition at line 19 of file DigiPileupCorrHistogramMaker.h.

Constructor & Destructor Documentation

DigiPileupCorrHistogramMaker::DigiPileupCorrHistogramMaker ( edm::ConsumesCollector && iC )

Definition at line 15 of file DigiPileupCorrHistogramMaker.cc.

     : m_pileupcollectionToken(iC.consumes<std::vector<PileupSummaryInfo> >(edm::InputTag("addPileupInfo"))),
       m_useVisibleVertices(false),
       m_hitname(),
       m_nbins(500),
       m_scalefact(),
       m_binmax(),
       m_labels(),
       m_nmultvsmclumi(),
       m_nmultvsmclumiprof(),
       m_nmultvsmcnvtx(),
       m_nmultvsmcnvtxprof(),
       m_nmultvsmcnvtxprof2d(),
       m_subdirs(),
       m_2dhisto(),
       m_ootBX(-1) {}

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

Definition at line 32 of file DigiPileupCorrHistogramMaker.cc.

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

     : m_pileupcollectionToken(
           iC.consumes<std::vector<PileupSummaryInfo> >(iConfig.getParameter<edm::InputTag>("pileupSummaryCollection"))),
       m_useVisibleVertices(iConfig.getParameter<bool>("useVisibleVertices")),
       m_hitname(iConfig.getUntrackedParameter<std::string>("hitName", "digi")),
       m_nbins(iConfig.getUntrackedParameter<int>("numberOfBins", 500)),
       m_scalefact(iConfig.getUntrackedParameter<int>("scaleFactor", 5)),
       m_labels(),
       m_nmultvsmclumi(),
       m_nmultvsmclumiprof(),
       m_nmultvsmcnvtx(),
       m_nmultvsmcnvtxprof(),
       m_nmultvsmcnvtxprof2d(),
       m_subdirs(),
       m_2dhisto(iConfig.getUntrackedParameter<bool>("wanted2DHisto", false)),
       m_ootBX(iConfig.getUntrackedParameter<int>("ootBX", -1)) {
   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");
   }
 }

DigiPileupCorrHistogramMaker::~DigiPileupCorrHistogramMaker ( )

Definition at line 58 of file DigiPileupCorrHistogramMaker.cc.

References mps_fire::i, 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];
   }
 }

Member Function Documentation

void DigiPileupCorrHistogramMaker::beginRun ( const unsigned int nrun )

Definition at line 137 of file DigiPileupCorrHistogramMaker.cc.

137 {}

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

Definition at line 67 of file DigiPileupCorrHistogramMaker.cc.

References SummaryClient_cfi::labels, and m_labels.

Referenced by MultiplicityInvestigator::MultiplicityInvestigator().

                                                                                                               {
   m_labels = labels;
   book(dirname);
 }

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

Definition at line 72 of file DigiPileupCorrHistogramMaker.cc.

References mps_fire::i, m_2dhisto, m_binmax, m_hitname, m_labels, m_nbins, m_nmultvsmclumi, m_nmultvsmclumiprof, m_nmultvsmcnvtx, m_nmultvsmcnvtxprof, m_nmultvsmcnvtxprof2d, m_scalefact, m_subdirs, TFileDirectory::mkdir(), TFileService::mkdir(), Skims_PA_cff::name, SiStripTKNumbers::nstrips(), AlCaHLTBitMon_QueryRunRegistry::string, and overlapproblemtsosanalyzer_cfi::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("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));
 
     if (m_subdirs[i]) {
       sprintf(name, "n%sdigivsmclumi", slab.c_str());
       sprintf(title, "%s %s multiplicity ave pileup interactions", slab.c_str(), m_hitname.c_str());
       m_nmultvsmclumi[i] = m_subdirs[i]->make<TH2F>(
           name, title, 200, 0., 50., m_nbins, 0., m_binmax[i] / (m_scalefact * m_nbins) * m_nbins);
       m_nmultvsmclumi[i]->GetXaxis()->SetTitle("Average Pileup Interactions");
       m_nmultvsmclumi[i]->GetYaxis()->SetTitle("Number of Hits");
       sprintf(name, "n%sdigivsmclumiprof", slab.c_str());
       m_nmultvsmclumiprof[i] = m_subdirs[i]->make<TProfile>(name, title, 200, 0., 50.);
       m_nmultvsmclumiprof[i]->GetXaxis()->SetTitle("Average Pileup Interactions");
       m_nmultvsmclumiprof[i]->GetYaxis()->SetTitle("Number of Hits");
 
       sprintf(name, "n%sdigivsmcnvtx", slab.c_str());
       sprintf(title, "%s %s multiplicity vs pileup interactions", slab.c_str(), m_hitname.c_str());
       m_nmultvsmcnvtx[i] = m_subdirs[i]->make<TH2F>(
           name, title, 60, -0.5, 59.5, m_nbins, 0., m_binmax[i] / (m_scalefact * m_nbins) * m_nbins);
       m_nmultvsmcnvtx[i]->GetXaxis()->SetTitle("Pileup Interactions");
       m_nmultvsmcnvtx[i]->GetYaxis()->SetTitle("Number of Hits");
       sprintf(name, "n%sdigivsmcnvtxprof", slab.c_str());
       m_nmultvsmcnvtxprof[i] = m_subdirs[i]->make<TProfile>(name, title, 60, -0.5, 59.5);
       m_nmultvsmcnvtxprof[i]->GetXaxis()->SetTitle("Pileup Interactions");
       m_nmultvsmcnvtxprof[i]->GetYaxis()->SetTitle("Number of Hits");
       if (m_2dhisto) {
         sprintf(name, "n%sdigivsmcnvtxprof2d", slab.c_str());
         m_nmultvsmcnvtxprof2d[i] = m_subdirs[i]->make<TProfile2D>(name, title, 60, -0.5, 59.5, 60, -0.5, 59.5);
         m_nmultvsmcnvtxprof2d[i]->GetXaxis()->SetTitle("Pileup Interactions");
         m_nmultvsmcnvtxprof2d[i]->GetYaxis()->SetTitle("Pileup Interactions OOT");
         m_nmultvsmcnvtxprof2d[i]->GetZaxis()->SetTitle("Number of Hits");
       }
     }
   }
 }

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

Definition at line 139 of file DigiPileupCorrHistogramMaker.cc.

References edm::Event::getByToken(), mps_fire::i, m_2dhisto, m_labels, m_nmultvsmclumi, m_nmultvsmclumiprof, m_nmultvsmcnvtx, m_nmultvsmcnvtxprof, m_nmultvsmcnvtxprof2d, m_ootBX, m_pileupcollectionToken, and m_useVisibleVertices.

Referenced by MultiplicityInvestigator::analyze().

                                                                                                       {
   edm::Handle<std::vector<PileupSummaryInfo> > pileupinfos;
   iEvent.getByToken(m_pileupcollectionToken, pileupinfos);
 
   // look for the intime PileupSummaryInfo
 
   std::vector<PileupSummaryInfo>::const_iterator pileupinfoInTime = pileupinfos->end();
   std::vector<PileupSummaryInfo>::const_iterator pileupinfoMinusOne = pileupinfos->end();
 
   for (std::vector<PileupSummaryInfo>::const_iterator pileupinfo = pileupinfos->begin();
        pileupinfo != pileupinfos->end();
        ++pileupinfo) {
     if (pileupinfo->getBunchCrossing() == 0)
       pileupinfoInTime = pileupinfo;
     if (pileupinfo->getBunchCrossing() == m_ootBX)
       pileupinfoMinusOne = pileupinfo;
   }
 
   if (pileupinfoInTime == pileupinfos->end()) {
     edm::LogError("NoInTimePileUpInfo") << "Cannot find the in-time pileup info ";
 
   } else {
     int npileup = pileupinfoInTime->getPU_NumInteractions();
 
     if (m_useVisibleVertices)
       npileup = pileupinfoInTime->getPU_zpositions().size();
 
     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_nmultvsmcnvtx[i]->Fill(npileup, digi->second);
         m_nmultvsmcnvtxprof[i]->Fill(npileup, digi->second);
         m_nmultvsmclumi[i]->Fill(pileupinfoInTime->getTrueNumInteractions(), digi->second);
         m_nmultvsmclumiprof[i]->Fill(pileupinfoInTime->getTrueNumInteractions(), digi->second);
         if (m_2dhisto) {
           if (pileupinfoMinusOne != pileupinfos->end()) {
             int npileupminusone = pileupinfoMinusOne->getPU_NumInteractions();
             if (m_useVisibleVertices)
               npileupminusone = pileupinfoMinusOne->getPU_zpositions().size();
             m_nmultvsmcnvtxprof2d[i]->Fill(npileup, npileupminusone, digi->second);
           }
         }
       }
     }
   }
 }