DigiBXCorrHistogramMaker< T > Class Template Reference

#include <DigiBXCorrHistogramMaker.h>

Public Member Functions
void	beginRun (const unsigned int nrun)
void	book (const char *dirname, const std::map< int, std::string > &labels)
	DigiBXCorrHistogramMaker (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC)
	DigiBXCorrHistogramMaker (edm::ConsumesCollector &&iC, const int ncorbins=1000)
void	fill (const T &he, const std::map< int, int > &ndigi, const edm::Handle< APVCyclePhaseCollection > &phase)
void	fillcorr (const T &he1, const T &he2, const std::map< int, int > &ndigi)
	~DigiBXCorrHistogramMaker ()

Private Attributes
std::map< unsigned int, int >	m_binmax
const bool	m_dbx3Histo
const bool	m_dbx3Histo3D
std::map< int, TProfile * >	m_digicorr
std::string	m_hitname
std::map< int, std::string >	m_labels
const int	m_nbins
int	m_ncorbins
std::map< int, TProfile * >	m_ndigivsbx
std::map< int, TH2F * >	m_ndigivsbx2D
std::map< int, TH2F * >	m_ndigivsbx2Dzoom
std::map< int, TH2F * >	m_ndigivsbx2Dzoom2
std::map< int, TH2F * >	m_ndigivscycle
std::map< int, TProfile2D * >	m_ndigivscycle2dbx
std::map< int, TProfile2D * >	m_ndigivscycledbx
std::map< int, TProfile2D ** >	m_ndigivscycletime
std::map< int, TH2F * >	m_ndigivscyclezoom
std::map< int, TH2F * >	m_ndigivscyclezoom2
std::map< int, TProfile * >	m_ndigivsdbx
std::map< int, TH2F * >	m_ndigivsdbx2D
std::map< int, TH2F * >	m_ndigivsdbx2Dzoom
std::map< int, TH2F * >	m_ndigivsdbx2Dzoom2
std::map< int, TProfile2D * >	m_ndigivsdbx3zoom
std::map< int, TProfile * >	m_ndigivsdbxincycle
std::map< int, TH2F * >	m_ndigivsdbxincycle2D
std::map< int, TProfile2D * >	m_ndigivsdbxincycle3
std::map< int, TH3F * >	m_ndigivsdbxincycle33D
std::map< int, TProfile * >	m_ndigivsdbxzoom
std::map< int, TProfile * >	m_ndigivsdbxzoom2
std::map< int, TProfile * >	m_nmeandigivscycle
std::map< int, std::string >	m_phasepart
RunHistogramManager	m_rhm
const bool	m_runHisto
std::vector< int >	m_scalefact

Detailed Description

template<class T>
class DigiBXCorrHistogramMaker< T >

Definition at line 23 of file DigiBXCorrHistogramMaker.h.

Constructor & Destructor Documentation

DigiBXCorrHistogramMaker() [1/2]

template<class T >

DigiBXCorrHistogramMaker< T >::DigiBXCorrHistogramMaker	(	edm::ConsumesCollector &&	iC,
		const int	ncorbins = 1000
	)

Definition at line 86 of file DigiBXCorrHistogramMaker.h.

    : m_ncorbins(ncorbins),
      m_hitname("digi"),
      m_dbx3Histo(false),
      m_dbx3Histo3D(false),
      m_runHisto(true),
      m_labels(),
      m_binmax(),
      m_phasepart(),
      m_scalefact(),
      m_nbins(200),
      m_rhm(iC) {}

DigiBXCorrHistogramMaker() [2/2]

template<class T >

DigiBXCorrHistogramMaker< T >::DigiBXCorrHistogramMaker	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 100 of file DigiBXCorrHistogramMaker.h.

    : m_ncorbins(iConfig.getUntrackedParameter<int>("corrNbins")),
      m_hitname(iConfig.getUntrackedParameter<std::string>("hitName", "digi")),
      m_dbx3Histo(iConfig.getUntrackedParameter<bool>("dbx3Histo", false)),
      m_dbx3Histo3D(iConfig.getUntrackedParameter<bool>("dbx3Histo3D", false)),
      m_runHisto(iConfig.getUntrackedParameter<bool>("runHisto", true)),
      m_labels(),
      m_scalefact(iConfig.getUntrackedParameter<std::vector<int> >("scaleFactors", std::vector<int>(1, 5))),
      m_nbins(iConfig.getUntrackedParameter<int>("numberOfBins", 200)),
      m_rhm(iC) {
  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_binmax[ps->getParameter<unsigned int>("detSelection")] = ps->getParameter<int>("binMax");
    m_phasepart[ps->getParameter<unsigned int>("detSelection")] =
        ps->getUntrackedParameter<std::string>("phasePartition", "None");
  }
}

~DigiBXCorrHistogramMaker()

template<class T>

DigiBXCorrHistogramMaker< T >::~DigiBXCorrHistogramMaker ( )

inline

Definition at line 28 of file DigiBXCorrHistogramMaker.h.

{};

Member Function Documentation

beginRun()

template<class T >

void DigiBXCorrHistogramMaker< T >::beginRun

(

const unsigned int

nrun

)

Definition at line 354 of file DigiBXCorrHistogramMaker.h.

                                                                  {
  m_rhm.beginRun(nrun);
 
  for (std::map<int, std::string>::const_iterator lab = m_labels.begin(); lab != m_labels.end(); lab++) {
    const int i = lab->first;
    if (m_runHisto) {
      if (m_ndigivscycletime[i]) {
        (*m_ndigivscycletime[i])->GetXaxis()->SetTitle("Event 1 BX mod(70)");
        (*m_ndigivscycletime[i])->GetYaxis()->SetTitle("time [Orb#]");
        (*m_ndigivscycletime[i])->SetCanExtend(TH1::kAllAxes);
      }
    }
  }
}

Referenced by MultiplicityTimeCorrelations::beginRun().

book()

template<class T >

void DigiBXCorrHistogramMaker< T >::book	(	const char *	dirname,
		const std::map< int, std::string > &	labels
	)

Definition at line 121 of file DigiBXCorrHistogramMaker.h.

                                                                                                {
  m_labels = labels;
 
  edm::Service<TFileService> tfserv;
  TFileDirectory subev = tfserv->mkdir(dirname);
 
  SiStripTKNumbers trnumb;
 
  edm::LogInfo("NumberOfBins") << "Number of Bins: " << m_nbins;
  edm::LogInfo("ScaleFactors") << "x-axis range scale factors: ";
  for (std::vector<int>::const_iterator sf = m_scalefact.begin(); sf != m_scalefact.end(); ++sf) {
    edm::LogVerbatim("ScaleFactors") << *sf;
  }
  edm::LogInfo("BinMaxValue") << "Setting bin max values";
 
  for (std::map<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;
    const unsigned int ui = i;
 
    if (m_binmax.find(ui) == m_binmax.end()) {
      edm::LogVerbatim("NotConfiguredBinMax")
          << "Bin max for " << lab->second << " not configured: " << trnumb.nstrips(i) << " used";
      m_binmax[ui] = trnumb.nstrips(i);
    }
 
    edm::LogVerbatim("BinMaxValue") << "Bin max for " << lab->second << " is " << m_binmax[ui];
  }
 
  edm::LogInfo("PhasePartitions") << "Partitions for APV Cycle Phase";
 
  for (std::map<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;
    const unsigned int ui = i;
    edm::LogVerbatim("PhasePartitions") << "Partition for " << lab->second << " is "
                                        << ((m_phasepart.find(ui) != m_phasepart.end()) ? m_phasepart[ui]
                                                                                        : "not found");
  }
 
  for (std::map<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;
    const unsigned int ui = i;
 
    char name[200];
    char title[500];
 
    // vs DBX
 
    if (!m_scalefact.empty()) {
      sprintf(title, "%s %s multiplicity vs BX separation", slab.c_str(), m_hitname.c_str());
      sprintf(name, "n%sdigivsdbx2D", slab.c_str());
      m_ndigivsdbx2D[i] = subev.make<TH2F>(
          name, title, 1000, -0.5, 500000 - 0.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[0] * m_nbins)) * m_nbins);
      sprintf(name, "n%sdigivsdbx2Dzoom2", slab.c_str());
      m_ndigivsdbx2Dzoom2[i] = subev.make<TH2F>(
          name, title, 1000, -0.5, 50000 - 0.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[0] * m_nbins)) * m_nbins);
      sprintf(name, "n%sdigivsdbx2Dzoom", slab.c_str());
      m_ndigivsdbx2Dzoom[i] = subev.make<TH2F>(
          name, title, 1000, -0.5, 999.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[0] * m_nbins)) * m_nbins);
 
      m_ndigivsdbx2D[i]->GetXaxis()->SetTitle("#DeltaBX");
      m_ndigivsdbx2D[i]->GetYaxis()->SetTitle("Number of Hits");
      m_ndigivsdbx2Dzoom2[i]->GetXaxis()->SetTitle("#DeltaBX");
      m_ndigivsdbx2Dzoom2[i]->GetYaxis()->SetTitle("Number of Hits");
      m_ndigivsdbx2Dzoom[i]->GetXaxis()->SetTitle("#DeltaBX");
      m_ndigivsdbx2Dzoom[i]->GetYaxis()->SetTitle("Number of Hits");
    }
 
    sprintf(title, "%s %s multiplicity vs BX separation", slab.c_str(), m_hitname.c_str());
    sprintf(name, "n%sdigivsdbx", slab.c_str());
    m_ndigivsdbx[i] = subev.make<TProfile>(name, title, 1000, -0.5, 500000. - 0.5);
    sprintf(name, "n%sdigivsdbxzoom2", slab.c_str());
    m_ndigivsdbxzoom2[i] = subev.make<TProfile>(name, title, 1000, -0.5, 50000. - 0.5);
    sprintf(name, "n%sdigivsdbxzoom", slab.c_str());
    m_ndigivsdbxzoom[i] = subev.make<TProfile>(name, title, 1000, -0.5, 999.5);
    m_ndigivsdbx[i]->GetXaxis()->SetTitle("#DeltaBX");
    m_ndigivsdbx[i]->GetYaxis()->SetTitle("Number of Hits");
    m_ndigivsdbxzoom2[i]->GetXaxis()->SetTitle("#DeltaBX");
    m_ndigivsdbxzoom2[i]->GetYaxis()->SetTitle("Number of Hits");
    m_ndigivsdbxzoom[i]->GetXaxis()->SetTitle("#DeltaBX");
    m_ndigivsdbxzoom[i]->GetYaxis()->SetTitle("Number of Hits");
 
    sprintf(name, "n%sdigivsdbx3zoom", slab.c_str());
    sprintf(title, "%s %s multiplicity vs Triplets BX separation", slab.c_str(), m_hitname.c_str());
    m_ndigivsdbx3zoom[i] = subev.make<TProfile2D>(name, title, 100, -0.5, 999.5, 100, -0.5, 999.5);
    m_ndigivsdbx3zoom[i]->GetXaxis()->SetTitle("#DeltaBX(n,n-1)");
    m_ndigivsdbx3zoom[i]->GetYaxis()->SetTitle("#DeltaBX(n,n-2)");
 
    sprintf(name, "%sdigicorr", slab.c_str());
    sprintf(title, "%s %s DBX correlation", slab.c_str(), m_hitname.c_str());
    m_digicorr[i] = subev.make<TProfile>(name, title, m_ncorbins, -0.5, m_ncorbins - 0.5);
    m_digicorr[i]->GetXaxis()->SetTitle("#DeltaBX");
    m_digicorr[i]->GetYaxis()->SetTitle("Number of Hits");
 
    // vs DBX w.r.t. cycle
 
    if (!m_scalefact.empty()) {
      if (m_phasepart.find(ui) != m_phasepart.end() && m_phasepart[ui] != "None") {
        sprintf(name, "n%sdigivsdbxincycle", slab.c_str());
        sprintf(title, "%s %s multiplicity vs BX separation w.r.t. cycle", slab.c_str(), m_hitname.c_str());
        m_ndigivsdbxincycle[i] = subev.make<TProfile>(name, title, 1000, -0.5, 999.5);
        m_ndigivsdbxincycle[i]->GetXaxis()->SetTitle("#DeltaBX w.r.t. cycle");
        m_ndigivsdbxincycle[i]->GetYaxis()->SetTitle("Number of Hits");
 
        sprintf(name, "n%sdigivsdbxincycle2D", slab.c_str());
        sprintf(title, "%s %s multiplicity vs BX separation w.r.t. cycle", slab.c_str(), m_hitname.c_str());
        m_ndigivsdbxincycle2D[i] = subev.make<TH2F>(
            name, title, 1000, -0.5, 999.5, m_nbins, 0., (1 + m_binmax[ui] / (m_scalefact[0] * m_nbins)) * m_nbins);
        m_ndigivsdbxincycle2D[i]->GetXaxis()->SetTitle("#DeltaBX w.r.t. cycle");
        m_ndigivsdbxincycle2D[i]->GetYaxis()->SetTitle("Number of Hits");
 
        if (m_dbx3Histo) {
          sprintf(name, "n%sdigivsdbxincycle3", slab.c_str());
          sprintf(title, "%s %s multiplicity vs Triplets BX separation w.r.t. cycle", slab.c_str(), m_hitname.c_str());
          m_ndigivsdbxincycle3[i] = subev.make<TProfile2D>(name, title, 2000, -0.5, 1999.5, 30, -0.5, 2099.5);
          m_ndigivsdbxincycle3[i]->GetXaxis()->SetTitle("#DeltaBX(n,n-1)");
          m_ndigivsdbxincycle3[i]->GetYaxis()->SetTitle("#DeltaBX(n,n-2)-#DeltaBX(n,n-1)");
 
          if (m_dbx3Histo3D) {
            sprintf(name, "n%sdigivsdbxincycle33D", slab.c_str());
            sprintf(
                title, "%s %s multiplicity vs Triplets BX separation w.r.t. cycle", slab.c_str(), m_hitname.c_str());
            m_ndigivsdbxincycle33D[i] = subev.make<TH3F>(name,
                                                         title,
                                                         2000,
                                                         -0.5,
                                                         1999.5,
                                                         30,
                                                         -0.5,
                                                         2099.5,
                                                         50,
                                                         0.,
                                                         (1 + m_binmax[ui] / (m_scalefact[0] * 50)) * 50);
            m_ndigivsdbxincycle33D[i]->GetXaxis()->SetTitle("#DeltaBX(n,n-1)");
            m_ndigivsdbxincycle33D[i]->GetYaxis()->SetTitle("#DeltaBX(n,n-2)-#DeltaBX(n,n-1)");
          }
        }
      }
    }
 
    // vs absolute BX mod 70
 
    if (m_phasepart.find(ui) != m_phasepart.end() && m_phasepart[ui] != "None") {
      sprintf(title, "%s Mean %s multiplicity vs BX mod(70)", slab.c_str(), m_hitname.c_str());
      sprintf(name, "n%smeandigivscycle", slab.c_str());
      m_nmeandigivscycle[i] = subev.make<TProfile>(name, title, 70, -0.5, 69.5);
      m_nmeandigivscycle[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
      m_nmeandigivscycle[i]->GetYaxis()->SetTitle("Mean number of Hits");
 
      sprintf(title, "%s %s multiplicity vs BX mod(70)", slab.c_str(), m_hitname.c_str());
 
      if (!m_scalefact.empty()) {
        sprintf(name, "n%sdigivscycle", slab.c_str());
        m_ndigivscycle[i] = subev.make<TH2F>(
            name, title, 70, -0.5, 69.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[0] * m_nbins)) * m_nbins);
        m_ndigivscycle[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
        m_ndigivscycle[i]->GetYaxis()->SetTitle("Number of Hits");
      }
      if (m_scalefact.size() >= 2) {
        sprintf(name, "n%sdigivscyclezoom", slab.c_str());
        m_ndigivscyclezoom[i] = subev.make<TH2F>(
            name, title, 70, -0.5, 69.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[1] * m_nbins)) * m_nbins);
        m_ndigivscyclezoom[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
        m_ndigivscyclezoom[i]->GetYaxis()->SetTitle("Number of Hits");
      }
      if (m_scalefact.size() >= 3) {
        sprintf(name, "n%sdigivscyclezoom2", slab.c_str());
        m_ndigivscyclezoom2[i] = subev.make<TH2F>(
            name, title, 70, -0.5, 69.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[2] * m_nbins)) * m_nbins);
        m_ndigivscyclezoom2[i]->GetXaxis()->SetTitle("absolute BX mod(70)");
        m_ndigivscyclezoom2[i]->GetYaxis()->SetTitle("Number of Hits");
      }
 
      sprintf(name, "n%sdigivscycledbx", slab.c_str());
      sprintf(title, "%s %s multiplicity vs BX mod(70) and DBX", slab.c_str(), m_hitname.c_str());
      m_ndigivscycledbx[i] = subev.make<TProfile2D>(name, title, 70, -0.5, 69.5, 1000, -0.5, 999.5);
      m_ndigivscycledbx[i]->GetXaxis()->SetTitle("Event 1 BX mod(70)");
      m_ndigivscycledbx[i]->GetYaxis()->SetTitle("#DeltaBX event 2-1");
 
      sprintf(name, "n%sdigivscycle2dbx", slab.c_str());
      sprintf(title, "%s %s multiplicity vs BX mod(70) and DBX", slab.c_str(), m_hitname.c_str());
      m_ndigivscycle2dbx[i] = subev.make<TProfile2D>(name, title, 70, -0.5, 69.5, 1000, -0.5, 999.5);
      m_ndigivscycle2dbx[i]->GetXaxis()->SetTitle("Event 2 BX mod(70)");
      m_ndigivscycle2dbx[i]->GetYaxis()->SetTitle("#DeltaBX event 2-1");
    }
 
    // Multiplicity in cycle vs time is booked also if the phase is not corrected
 
    if (m_runHisto) {
      sprintf(name, "n%sdigivscycletime", slab.c_str());
      sprintf(title, "%s %s multiplicity vs BX mod(70) and Orbit", slab.c_str(), m_hitname.c_str());
      m_ndigivscycletime[i] = m_rhm.makeTProfile2D(name, title, 70, -0.5, 69.5, 90, 0., 90 * 262144);
      //      m_ndigivscycletime[i]->GetXaxis()->SetTitle("Event 1 BX mod(70)"); m_ndigivscycletime[i]->GetYaxis()->SetTitle("time [Orb#]");
      //      m_ndigivscycletime[i]->SetCanExtend(TH1::kYaxis);
    }
 
    // vs BX number
 
    sprintf(title, "%s %s mean multiplicity vs BX", slab.c_str(), m_hitname.c_str());
    sprintf(name, "n%sdigivsbx", slab.c_str());
    m_ndigivsbx[i] = subev.make<TProfile>(name, title, 3564, -0.5, 3563.5);
    m_ndigivsbx[i]->GetXaxis()->SetTitle("BX#");
    m_ndigivsbx[i]->GetYaxis()->SetTitle("Mean Number of Hits");
 
    sprintf(title, "%s %s multiplicity vs BX", slab.c_str(), m_hitname.c_str());
 
    if (!m_scalefact.empty()) {
      sprintf(name, "n%sdigivsbx2D", slab.c_str());
      m_ndigivsbx2D[i] = subev.make<TH2F>(
          name, title, 3564, -0.5, 3563.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[0] * m_nbins)) * m_nbins);
      m_ndigivsbx2D[i]->GetXaxis()->SetTitle("BX#");
      m_ndigivsbx2D[i]->GetYaxis()->SetTitle("Number of Hits");
    }
    if (m_scalefact.size() >= 2) {
      sprintf(name, "n%sdigivsbx2Dzoom", slab.c_str());
      m_ndigivsbx2Dzoom[i] = subev.make<TH2F>(
          name, title, 3564, -0.5, 3563.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[1] * m_nbins)) * m_nbins);
      m_ndigivsbx2Dzoom[i]->GetXaxis()->SetTitle("BX#");
      m_ndigivsbx2Dzoom[i]->GetYaxis()->SetTitle("Number of Hits");
    }
    if (m_scalefact.size() >= 3) {
      sprintf(name, "n%sdigivsbx2Dzoom2", slab.c_str());
      m_ndigivsbx2Dzoom2[i] = subev.make<TH2F>(
          name, title, 3564, -0.5, 3563.5, m_nbins, 0, (1 + m_binmax[ui] / (m_scalefact[2] * m_nbins)) * m_nbins);
      m_ndigivsbx2Dzoom2[i]->GetXaxis()->SetTitle("BX#");
      m_ndigivsbx2Dzoom2[i]->GetYaxis()->SetTitle("Number of Hits");
    }
  }
}

Referenced by MultiplicityTimeCorrelations::MultiplicityTimeCorrelations().

fill()

template<class T>

void DigiBXCorrHistogramMaker< T >::fill	(	const T &	he,
		const std::map< int, int > &	ndigi,
		const edm::Handle< APVCyclePhaseCollection > &	phase
	)

Definition at line 370 of file DigiBXCorrHistogramMaker.h.

                                                                                        {
  for (std::map<int, int>::const_iterator digi = ndigi.begin(); digi != ndigi.end(); digi++) {
    if (m_labels.find(digi->first) != m_labels.end()) {
      const int i = digi->first;
      const unsigned int ui = i;
 
      int thephase = APVCyclePhaseCollection::invalid;
      if (m_phasepart.find(ui) != m_phasepart.end() && m_phasepart[ui] != "None") {
        if (!phase.failedToGet() && phase.isValid()) {
          thephase = phase->getPhase(m_phasepart[ui]);
        }
      }
 
      long long tbx = he.absoluteBX();
      if (thephase != APVCyclePhaseCollection::nopartition && thephase != APVCyclePhaseCollection::multiphase &&
          thephase != APVCyclePhaseCollection::invalid) {
        tbx -= thephase;
 
        if (m_nmeandigivscycle.find(i) != m_nmeandigivscycle.end())
          m_nmeandigivscycle[i]->Fill(tbx % 70, digi->second);
 
        if (m_ndigivscycle.find(i) != m_ndigivscycle.end())
          m_ndigivscycle[i]->Fill(tbx % 70, digi->second);
        if (m_ndigivscyclezoom.find(i) != m_ndigivscyclezoom.end())
          m_ndigivscyclezoom[i]->Fill(tbx % 70, digi->second);
        if (m_ndigivscyclezoom2.find(i) != m_ndigivscyclezoom2.end())
          m_ndigivscyclezoom2[i]->Fill(tbx % 70, digi->second);
      }
 
      if (m_runHisto) {
        if (m_ndigivscycletime.find(i) != m_ndigivscycletime.end()) {
          if (m_ndigivscycletime[i] != nullptr && (*m_ndigivscycletime[i]) != nullptr)
            (*m_ndigivscycletime[i])->Fill(tbx % 70, (int)he._orbit, digi->second);
        }
      }
 
      m_ndigivsbx[i]->Fill(he.bx() % 3564, digi->second);
      if (m_ndigivsbx2D.find(i) != m_ndigivsbx2D.end())
        m_ndigivsbx2D[i]->Fill(he.bx() % 3564, digi->second);
      if (m_ndigivsbx2Dzoom.find(i) != m_ndigivsbx2Dzoom.end())
        m_ndigivsbx2Dzoom[i]->Fill(he.bx() % 3564, digi->second);
      if (m_ndigivsbx2Dzoom2.find(i) != m_ndigivsbx2Dzoom2.end())
        m_ndigivsbx2Dzoom2[i]->Fill(he.bx() % 3564, digi->second);
 
      if (he.depth() > 0) {
        long long dbx = he.deltaBX();
 
        m_ndigivsdbx[i]->Fill(dbx, digi->second);
        m_ndigivsdbxzoom[i]->Fill(dbx, digi->second);
        m_ndigivsdbxzoom2[i]->Fill(dbx, digi->second);
 
        if (m_ndigivsdbx2D.find(i) != m_ndigivsdbx2D.end())
          m_ndigivsdbx2D[i]->Fill(dbx, digi->second);
        if (m_ndigivsdbx2Dzoom.find(i) != m_ndigivsdbx2Dzoom.end())
          m_ndigivsdbx2Dzoom[i]->Fill(dbx, digi->second);
        if (m_ndigivsdbx2Dzoom2.find(i) != m_ndigivsdbx2Dzoom2.end())
          m_ndigivsdbx2Dzoom2[i]->Fill(dbx, digi->second);
 
        long long prevtbx = he.absoluteBX(1);
        if (thephase != APVCyclePhaseCollection::nopartition && thephase != APVCyclePhaseCollection::multiphase &&
            thephase != APVCyclePhaseCollection::invalid) {
          long long dbxincycle = he.deltaBXinCycle(thephase);
          if (m_ndigivsdbxincycle2D.find(i) != m_ndigivsdbxincycle2D.end())
            m_ndigivsdbxincycle2D[i]->Fill(dbxincycle, digi->second);
          if (m_ndigivsdbxincycle.find(i) != m_ndigivsdbxincycle.end())
            m_ndigivsdbxincycle[i]->Fill(dbxincycle, digi->second);
 
          prevtbx -= thephase;
          if (m_ndigivscycledbx.find(i) != m_ndigivscycledbx.end())
            m_ndigivscycledbx[i]->Fill(prevtbx % 70, dbx, digi->second);
          if (m_ndigivscycle2dbx.find(i) != m_ndigivscycle2dbx.end())
            m_ndigivscycle2dbx[i]->Fill(tbx % 70, dbx, digi->second);
        }
 
        if (he.depth() > 1) {
          long long dbx2 = he.deltaBX(2);
          m_ndigivsdbx3zoom[i]->Fill(dbx2, dbx, digi->second);
 
          if (thephase != APVCyclePhaseCollection::nopartition && thephase != APVCyclePhaseCollection::multiphase &&
              thephase != APVCyclePhaseCollection::invalid) {
            long long dbxincycle = he.deltaBXinCycle(thephase);
            long long dbxincycle2 = he.deltaBXinCycle(2, thephase);
            if (m_dbx3Histo) {
              if (m_ndigivsdbxincycle3.find(i) != m_ndigivsdbxincycle3.end())
                m_ndigivsdbxincycle3[i]->Fill(dbxincycle, dbxincycle2 - dbxincycle, digi->second);
              if (m_dbx3Histo3D) {
                if (m_ndigivsdbxincycle33D.find(i) != m_ndigivsdbxincycle33D.end())
                  m_ndigivsdbxincycle33D[i]->Fill(dbxincycle, dbxincycle2 - dbxincycle, digi->second);
              }
            }
          }
        }
      }
 
    } else {
      edm::LogWarning("MissingKey") << " Key " << digi->first << " is missing ";
    }
  }
}

Referenced by MultiplicityTimeCorrelations::analyze().

fillcorr()

template<class T>

void DigiBXCorrHistogramMaker< T >::fillcorr	(	const T &	he1,
		const T &	he2,
		const std::map< int, int > &	ndigi
	)

Definition at line 473 of file DigiBXCorrHistogramMaker.h.

                                                                                                    {
  for (std::map<int, int>::const_iterator digi = ndigi.begin(); digi != ndigi.end(); digi++) {
    if (m_labels.find(digi->first) != m_labels.end()) {
      const int i = digi->first;
 
      long long dbx = he2.deltaBX(he1);
      m_digicorr[i]->Fill(dbx, digi->second);
 
    } else {
      edm::LogWarning("MissingKey") << " Key " << digi->first << " is missing ";
    }
  }
}

Member Data Documentation

m_binmax

template<class T>

std::map<unsigned int, int> DigiBXCorrHistogramMaker< T >::m_binmax

private

Definition at line 43 of file DigiBXCorrHistogramMaker.h.

Referenced by DigiBXCorrHistogramMaker< EventWithHistory >::DigiBXCorrHistogramMaker().

m_dbx3Histo

template<class T>

const bool DigiBXCorrHistogramMaker< T >::m_dbx3Histo

private

Definition at line 38 of file DigiBXCorrHistogramMaker.h.

m_dbx3Histo3D

template<class T>

const bool DigiBXCorrHistogramMaker< T >::m_dbx3Histo3D

private

Definition at line 39 of file DigiBXCorrHistogramMaker.h.

m_digicorr

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_digicorr

private

Definition at line 82 of file DigiBXCorrHistogramMaker.h.

m_hitname

template<class T>

std::string DigiBXCorrHistogramMaker< T >::m_hitname

private

Definition at line 37 of file DigiBXCorrHistogramMaker.h.

m_labels

template<class T>

std::map<int, std::string> DigiBXCorrHistogramMaker< T >::m_labels

private

Definition at line 42 of file DigiBXCorrHistogramMaker.h.

m_nbins

template<class T>

const int DigiBXCorrHistogramMaker< T >::m_nbins

private

Definition at line 46 of file DigiBXCorrHistogramMaker.h.

m_ncorbins

template<class T>

int DigiBXCorrHistogramMaker< T >::m_ncorbins

private

Definition at line 36 of file DigiBXCorrHistogramMaker.h.

m_ndigivsbx

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx

private

Definition at line 63 of file DigiBXCorrHistogramMaker.h.

m_ndigivsbx2D

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx2D

private

Definition at line 64 of file DigiBXCorrHistogramMaker.h.

m_ndigivsbx2Dzoom

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx2Dzoom

private

Definition at line 65 of file DigiBXCorrHistogramMaker.h.

m_ndigivsbx2Dzoom2

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsbx2Dzoom2

private

Definition at line 66 of file DigiBXCorrHistogramMaker.h.

m_ndigivscycle

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivscycle

private

Definition at line 59 of file DigiBXCorrHistogramMaker.h.

m_ndigivscycle2dbx

template<class T>

std::map<int, TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivscycle2dbx

private

Definition at line 70 of file DigiBXCorrHistogramMaker.h.

m_ndigivscycledbx

template<class T>

std::map<int, TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivscycledbx

private

Definition at line 68 of file DigiBXCorrHistogramMaker.h.

m_ndigivscycletime

template<class T>

std::map<int, TProfile2D**> DigiBXCorrHistogramMaker< T >::m_ndigivscycletime

private

Definition at line 72 of file DigiBXCorrHistogramMaker.h.

m_ndigivscyclezoom

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivscyclezoom

private

Definition at line 60 of file DigiBXCorrHistogramMaker.h.

m_ndigivscyclezoom2

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivscyclezoom2

private

Definition at line 61 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbx

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx

private

Definition at line 50 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbx2D

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx2D

private

Definition at line 74 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbx2Dzoom

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx2Dzoom

private

Definition at line 76 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbx2Dzoom2

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx2Dzoom2

private

Definition at line 75 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbx3zoom

template<class T>

std::map<int, TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbx3zoom

private

Definition at line 78 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbxincycle

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle

private

Definition at line 54 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbxincycle2D

template<class T>

std::map<int, TH2F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle2D

private

Definition at line 55 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbxincycle3

template<class T>

std::map<int, TProfile2D*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle3

private

Definition at line 79 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbxincycle33D

template<class T>

std::map<int, TH3F*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxincycle33D

private

Definition at line 80 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbxzoom

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxzoom

private

Definition at line 52 of file DigiBXCorrHistogramMaker.h.

m_ndigivsdbxzoom2

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_ndigivsdbxzoom2

private

Definition at line 51 of file DigiBXCorrHistogramMaker.h.

m_nmeandigivscycle

template<class T>

std::map<int, TProfile*> DigiBXCorrHistogramMaker< T >::m_nmeandigivscycle

private

Definition at line 57 of file DigiBXCorrHistogramMaker.h.

m_phasepart

template<class T>

std::map<int, std::string> DigiBXCorrHistogramMaker< T >::m_phasepart

private

Definition at line 44 of file DigiBXCorrHistogramMaker.h.

Referenced by DigiBXCorrHistogramMaker< EventWithHistory >::DigiBXCorrHistogramMaker().

m_rhm

template<class T>

RunHistogramManager DigiBXCorrHistogramMaker< T >::m_rhm

private

Definition at line 48 of file DigiBXCorrHistogramMaker.h.

m_runHisto

template<class T>

const bool DigiBXCorrHistogramMaker< T >::m_runHisto

private

Definition at line 40 of file DigiBXCorrHistogramMaker.h.

m_scalefact

template<class T>

std::vector<int> DigiBXCorrHistogramMaker< T >::m_scalefact

private

Definition at line 45 of file DigiBXCorrHistogramMaker.h.