L1ExtraDQM::L1ExtraMonElement< CollectionType > Class Template Reference

#include <L1ExtraDQM.h>

Public Types
typedef CollectionType::const_iterator	CIterColl

Public Member Functions
void	bookhistograms (const edm::EventSetup &evSetup, DQMStore::IBooker &ibooker, const std::string &l1ExtraObject, const std::vector< L1GtObject > &l1GtObj, const bool bookPhi=true, const bool bookEta=true)
void	fillCharge (const CollectionType *collType, const bool validColl, const bool isL1Coll, const int bxInEvent)
	fill charge More...
void	fillEtPhiEta (const CollectionType *collType, const bool validColl, const bool bookPhi, const bool bookEta, const bool isL1Coll, const int bxInEvent)
	ET, eta, phi. More...
void	fillEtTotal (const CollectionType *collType, const bool validColl, const bool isL1Coll, const int bxInEvent)
	fill ET total in energy sums More...
void	fillHfBitCounts (const CollectionType *collType, const bool validColl, const int countIndex, const bool isL1Coll, const int bxInEvent)
	fill bit counts in HFRings collections More...
void	fillHfRingEtSums (const CollectionType *collType, const bool validColl, const int countIndex, const bool isL1Coll, const int bxInEvent)
	fill energy sums in HFRings collections More...
void	fillNrObjects (const CollectionType *collType, const bool validColl, const bool isL1Coll, const int bxInEvent)
	number of objects More...
void	fillPtPhiEta (const CollectionType *collType, const bool validColl, const bool bookPhi, const bool bookEta, const bool isL1Coll, const int bxInEvent)
	PT, eta, phi. More...
	L1ExtraMonElement (const L1GetHistLimits::Tokens &, const int)
virtual	~L1ExtraMonElement ()

Private Attributes
int	m_indexCharge
int	m_indexEt
int	m_indexEta
int	m_indexEtTotal
int	m_indexHfBitCounts
int	m_indexHfRingEtSums
int	m_indexNrObjects
	histogram index for each quantity, set during histogram booking More...
int	m_indexPhi
int	m_indexPt
std::vector< MonitorElement * >	m_monElement
L1GetHistLimits::Tokens	m_tokens

Detailed Description

template<class CollectionType>
class L1ExtraDQM::L1ExtraMonElement< CollectionType >

Definition at line 81 of file L1ExtraDQM.h.

Member Typedef Documentation

template<class CollectionType >

typedef CollectionType::const_iterator L1ExtraDQM::L1ExtraMonElement< CollectionType >::CIterColl

Definition at line 90 of file L1ExtraDQM.h.

Constructor & Destructor Documentation

template<class CollectionType >

L1ExtraDQM::L1ExtraMonElement< CollectionType >::L1ExtraMonElement	(	const L1GetHistLimits::Tokens &	tokens,
		const int	nrElements
	)

Definition at line 789 of file L1ExtraDQM.cc.

References L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_monElement.

    : m_tokens(tokens),
      m_indexNrObjects(-1),
      m_indexPt(-1),
      m_indexEt(-1),
      m_indexPhi(-1),
      m_indexEta(-1),
      m_indexEtTotal(-1),
      m_indexCharge(-1),
      m_indexHfBitCounts(-1),
      m_indexHfRingEtSums(-1) {
  m_monElement.reserve(nrElements);
}

template<class CollectionType >

L1ExtraDQM::L1ExtraMonElement< CollectionType >::~L1ExtraMonElement ( )

virtual

Definition at line 806 of file L1ExtraDQM.cc.

                                                                {
  //empty
}

Member Function Documentation

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::bookhistograms	(	const edm::EventSetup &	evSetup,
		DQMStore::IBooker &	ibooker,
		const std::string &	l1ExtraObject,
		const std::vector< L1GtObject > &	l1GtObj,
		const bool	bookPhi = true,
		const bool	bookEta = true
	)

Definition at line 811 of file L1ExtraDQM.cc.

References L1GetHistLimits::L1HistLimits::binThresholds, dqm::implementation::IBooker::book1D(), CenJet, filterCSVwithJSON::copy, ForJet, HfBitCounts, HfRingEtSums, IsoEG, l1extra::L1HFRings::kNumRings, L1GetHistLimits::l1HistLimits(), LogDebug, LogTrace, L1GetHistLimits::L1HistLimits::lowerBinValue, Mu, NoIsoEG, L1GetHistLimits::L1HistLimits::nrBins, AlCaHLTBitMon_QueryRunRegistry::string, TauJet, cond::impl::to_string(), and L1GetHistLimits::L1HistLimits::upperBinValue.

                                                                                       {
  // FIXME
  L1GtObject gtObj = l1GtObj.at(0);

  //
  std::string histName;
  std::string histTitle;
  std::string xAxisTitle;
  std::string yAxisTitle;

  std::string quantity = "";

  int indexHistogram = -1;

  if (gtObj == HfBitCounts) {
    L1GetHistLimits l1GetHistLimits(m_tokens, evSetup);
    const L1GetHistLimits::L1HistLimits& histLimits = l1GetHistLimits.l1HistLimits(gtObj, quantity);

    const int histNrBins = histLimits.nrBins;
    const double histMinValue = histLimits.lowerBinValue;
    const double histMaxValue = histLimits.upperBinValue;

    indexHistogram++;
    m_indexHfBitCounts = indexHistogram;

    for (int iCount = 0; iCount < l1extra::L1HFRings::kNumRings; ++iCount) {
      histName = l1ExtraObject + "_Count_" + std::to_string(iCount);
      histTitle = l1ExtraObject + ": count " + std::to_string(iCount);
      xAxisTitle = l1ExtraObject;
      yAxisTitle = "Entries";

      m_monElement.push_back(ibooker.book1D(histName, histTitle, histNrBins, histMinValue, histMaxValue));
      m_monElement[m_indexHfBitCounts + iCount]->setAxisTitle(xAxisTitle, 1);
      m_monElement[m_indexHfBitCounts + iCount]->setAxisTitle(yAxisTitle, 2);
    }

    return;
  }

  // number of objects per event
  if ((gtObj == Mu) || (gtObj == IsoEG) || (gtObj == NoIsoEG) || (gtObj == CenJet) || (gtObj == ForJet) ||
      (gtObj == TauJet)) {
    quantity = "NrObjects";

    L1GetHistLimits l1GetHistLimits(m_tokens, evSetup);
    const L1GetHistLimits::L1HistLimits& histLimits = l1GetHistLimits.l1HistLimits(gtObj, quantity);

    const int histNrBins = histLimits.nrBins;
    const double histMinValue = histLimits.lowerBinValue;
    const double histMaxValue = histLimits.upperBinValue;

    histName = l1ExtraObject + "_NrObjectsPerEvent";
    histTitle = l1ExtraObject + ": number of objects per event";
    xAxisTitle = "Nr_" + l1ExtraObject;
    yAxisTitle = "Entries";

    m_monElement.push_back(ibooker.book1D(histName, histTitle, histNrBins, histMinValue, histMaxValue));
    indexHistogram++;

    m_monElement[indexHistogram]->setAxisTitle(xAxisTitle, 1);
    m_monElement[indexHistogram]->setAxisTitle(yAxisTitle, 2);
    m_indexNrObjects = indexHistogram;
  }

  // transverse momentum (energy)  PT (ET) [GeV]

  quantity = "ET";
  std::string quantityLongName = " transverse energy ";

  if (gtObj == Mu) {
    quantity = "PT";
    quantityLongName = " transverse momentum ";
  }

  L1GetHistLimits l1GetHistLimits(m_tokens, evSetup);
  const L1GetHistLimits::L1HistLimits& histLimits = l1GetHistLimits.l1HistLimits(gtObj, quantity);

  const int histNrBinsET = histLimits.nrBins;
  const double histMinValueET = histLimits.lowerBinValue;
  const double histMaxValueET = histLimits.upperBinValue;
  const std::vector<float>& binThresholdsET = histLimits.binThresholds;

  float* binThresholdsETf;
  size_t sizeBinThresholdsET = binThresholdsET.size();
  binThresholdsETf = new float[sizeBinThresholdsET];
  copy(binThresholdsET.begin(), binThresholdsET.end(), binThresholdsETf);

  LogDebug("L1ExtraDQM") << "\n PT/ET histogram for " << l1ExtraObject << "\n histNrBinsET = " << histNrBinsET
                         << "\n histMinValueET = " << histMinValueET << "\n histMaxValueET = " << histMaxValueET
                         << "\n Last bin value represents the upper limit of the histogram" << std::endl;
  for (size_t iBin = 0; iBin < sizeBinThresholdsET; ++iBin) {
    LogTrace("L1ExtraDQM") << "Bin " << iBin << ": " << quantity << " = " << binThresholdsETf[iBin] << " GeV"
                           << std::endl;
  }

  histName = l1ExtraObject + "_" + quantity;
  histTitle = l1ExtraObject + ":" + quantityLongName + quantity + " [GeV]";
  xAxisTitle = l1ExtraObject + "_" + quantity + " [GeV]";
  yAxisTitle = "Entries";

  if (gtObj == HfRingEtSums) {
    indexHistogram++;
    m_indexHfRingEtSums = indexHistogram;

    for (int iCount = 0; iCount < l1extra::L1HFRings::kNumRings; ++iCount) {
      histName = l1ExtraObject + "_Count_" + std::to_string(iCount);
      histTitle = l1ExtraObject + ": count " + std::to_string(iCount);
      xAxisTitle = l1ExtraObject;
      yAxisTitle = "Entries";

      m_monElement.push_back(ibooker.book1D(histName, histTitle, histNrBinsET, binThresholdsETf));

      m_monElement[m_indexHfRingEtSums + iCount]->setAxisTitle(xAxisTitle, 1);
      m_monElement[m_indexHfRingEtSums + iCount]->setAxisTitle(yAxisTitle, 2);
    }

  } else {
    m_monElement.push_back(ibooker.book1D(histName, histTitle, histNrBinsET, binThresholdsETf));
    indexHistogram++;

    m_monElement[indexHistogram]->setAxisTitle(xAxisTitle, 1);
    m_monElement[indexHistogram]->setAxisTitle(yAxisTitle, 2);
    m_indexPt = indexHistogram;
    m_indexEt = indexHistogram;
    m_indexEtTotal = indexHistogram;
  }

  delete[] binThresholdsETf;

  //

  if (bookPhi) {
    quantity = "phi";

    // get limits and binning from L1Extra
    L1GetHistLimits l1GetHistLimits(m_tokens, evSetup);
    const L1GetHistLimits::L1HistLimits& histLimits = l1GetHistLimits.l1HistLimits(gtObj, quantity);

    const int histNrBinsPhi = histLimits.nrBins;
    const double histMinValuePhi = histLimits.lowerBinValue;
    const double histMaxValuePhi = histLimits.upperBinValue;
    const std::vector<float>& binThresholdsPhi = histLimits.binThresholds;

    float* binThresholdsPhif;
    size_t sizeBinThresholdsPhi = binThresholdsPhi.size();
    binThresholdsPhif = new float[sizeBinThresholdsPhi];
    copy(binThresholdsPhi.begin(), binThresholdsPhi.end(), binThresholdsPhif);

    LogDebug("L1ExtraDQM") << "\n phi histogram for " << l1ExtraObject << "\n histNrBinsPhi = " << histNrBinsPhi
                           << "\n histMinValuePhi = " << histMinValuePhi << "\n histMaxValuePhi = " << histMaxValuePhi
                           << "\n Last bin value represents the upper limit of the histogram" << std::endl;
    for (size_t iBin = 0; iBin < sizeBinThresholdsPhi; ++iBin) {
      LogTrace("L1ExtraDQM") << "Bin " << iBin << ": phi = " << binThresholdsPhif[iBin] << " deg" << std::endl;
    }

    histName = l1ExtraObject + "_phi";
    histTitle = l1ExtraObject + ": phi distribution ";
    xAxisTitle = l1ExtraObject + "_phi [deg]";
    yAxisTitle = "Entries";

    m_monElement.push_back(ibooker.book1D(histName, histTitle, histNrBinsPhi, histMinValuePhi, histMaxValuePhi));
    indexHistogram++;

    m_monElement[indexHistogram]->setAxisTitle(xAxisTitle, 1);
    m_monElement[indexHistogram]->setAxisTitle(yAxisTitle, 2);
    m_indexPhi = indexHistogram;

    delete[] binThresholdsPhif;
  }

  //

  if (bookEta) {
    quantity = "eta";

    // get limits and binning from L1Extra
    L1GetHistLimits l1GetHistLimits(m_tokens, evSetup);
    const L1GetHistLimits::L1HistLimits& histLimits = l1GetHistLimits.l1HistLimits(gtObj, quantity);

    const int histNrBinsEta = histLimits.nrBins;
    const double histMinValueEta = histLimits.lowerBinValue;
    const double histMaxValueEta = histLimits.upperBinValue;
    const std::vector<float>& binThresholdsEta = histLimits.binThresholds;

    //
    float* binThresholdsEtaf;
    size_t sizeBinThresholdsEta = binThresholdsEta.size();
    binThresholdsEtaf = new float[sizeBinThresholdsEta];
    copy(binThresholdsEta.begin(), binThresholdsEta.end(), binThresholdsEtaf);

    LogDebug("L1ExtraDQM") << "\n eta histogram for " << l1ExtraObject << "\n histNrBinsEta = " << histNrBinsEta
                           << "\n histMinValueEta = " << histMinValueEta << "\n histMaxValueEta = " << histMaxValueEta
                           << "\n Last bin value represents the upper limit of the histogram" << std::endl;
    for (size_t iBin = 0; iBin < sizeBinThresholdsEta; ++iBin) {
      LogTrace("L1ExtraDQM") << "Bin " << iBin << ": eta = " << binThresholdsEtaf[iBin] << std::endl;
    }

    histName = l1ExtraObject + "_eta";
    histTitle = l1ExtraObject + ": eta distribution ";
    xAxisTitle = l1ExtraObject + "_eta";
    yAxisTitle = "Entries";

    m_monElement.push_back(ibooker.book1D(histName, histTitle, histNrBinsEta, binThresholdsEtaf));
    indexHistogram++;

    m_monElement[indexHistogram]->setAxisTitle(xAxisTitle, 1);
    m_monElement[indexHistogram]->setAxisTitle(yAxisTitle, 2);
    m_indexEta = indexHistogram;

    delete[] binThresholdsEtaf;
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillCharge	(	const CollectionType *	collType,
		const bool	validColl,
		const bool	isL1Coll,
		const int	bxInEvent
	)

fill charge

Definition at line 1119 of file L1ExtraDQM.cc.

                                                                                    {
  if (validColl) {
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (isL1Coll && (iterColl->bx() != bxInEvent)) {
        continue;
      }

      m_monElement[m_indexCharge]->Fill(iterColl->charge());
    }
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillEtPhiEta	(	const CollectionType *	collType,
		const bool	validColl,
		const bool	bookPhi,
		const bool	bookEta,
		const bool	isL1Coll,
		const int	bxInEvent
	)

ET, eta, phi.

Definition at line 1076 of file L1ExtraDQM.cc.

References alignCSCRings::e, and rad2deg().

                                                                                      {
  if (validColl) {
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (isL1Coll && (iterColl->bx() != bxInEvent)) {
        continue;
      }

      m_monElement[m_indexEt]->Fill(iterColl->et());

      if (bookPhi) {
        // add a very small quantity to get off the bin edge
        m_monElement[m_indexPhi]->Fill(rad2deg(iterColl->phi()) + 1.e-6);
      }

      if (bookEta) {
        m_monElement[m_indexEta]->Fill(iterColl->eta());
      }
    }
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillEtTotal	(	const CollectionType *	collType,
		const bool	validColl,
		const bool	isL1Coll,
		const int	bxInEvent
	)

fill ET total in energy sums

Definition at line 1103 of file L1ExtraDQM.cc.

                                                                                     {
  if (validColl) {
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (isL1Coll && (iterColl->bx() != bxInEvent)) {
        continue;
      }

      m_monElement[m_indexEtTotal]->Fill(iterColl->etTotal());
    }
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillHfBitCounts	(	const CollectionType *	collType,
		const bool	validColl,
		const int	countIndex,
		const bool	isL1Coll,
		const int	bxInEvent
	)

fill bit counts in HFRings collections

Definition at line 1135 of file L1ExtraDQM.cc.

                                                                                         {
  if (validColl) {
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (isL1Coll && (iterColl->bx() != bxInEvent)) {
        continue;
      }

      m_monElement[m_indexHfBitCounts + countIndex]->Fill(
          iterColl->hfBitCount((l1extra::L1HFRings::HFRingLabels)countIndex));
    }
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillHfRingEtSums	(	const CollectionType *	collType,
		const bool	validColl,
		const int	countIndex,
		const bool	isL1Coll,
		const int	bxInEvent
	)

fill energy sums in HFRings collections

Definition at line 1153 of file L1ExtraDQM.cc.

                                                                                          {
  if (validColl) {
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (isL1Coll && (iterColl->bx() != bxInEvent)) {
        continue;
      }

      m_monElement[m_indexHfRingEtSums + countIndex]->Fill(
          iterColl->hfEtSum((l1extra::L1HFRings::HFRingLabels)countIndex));
    }
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillNrObjects	(	const CollectionType *	collType,
		const bool	validColl,
		const bool	isL1Coll,
		const int	bxInEvent
	)

number of objects

Definition at line 1030 of file L1ExtraDQM.cc.

                                                                                       {
  if (validColl && isL1Coll) {
    size_t collSize = 0;
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (iterColl->bx() == bxInEvent) {
        collSize++;
      }
    }
    m_monElement[m_indexNrObjects]->Fill(collSize);
  } else {
    size_t collSize = collType->size();
    m_monElement[m_indexNrObjects]->Fill(collSize);
  }
}

template<class CollectionType >

void L1ExtraDQM::L1ExtraMonElement< CollectionType >::fillPtPhiEta	(	const CollectionType *	collType,
		const bool	validColl,
		const bool	bookPhi,
		const bool	bookEta,
		const bool	isL1Coll,
		const int	bxInEvent
	)

PT, eta, phi.

Definition at line 1049 of file L1ExtraDQM.cc.

References alignCSCRings::e, and rad2deg().

                                                                                      {
  if (validColl) {
    for (CIterColl iterColl = collType->begin(); iterColl != collType->end(); ++iterColl) {
      if (isL1Coll && (iterColl->bx() != bxInEvent)) {
        continue;
      }

      m_monElement[m_indexPt]->Fill(iterColl->pt());

      if (bookPhi) {
        // add a very small quantity to get off the bin edge
        m_monElement[m_indexPhi]->Fill(rad2deg(iterColl->phi()) + 1.e-6);
      }

      if (bookEta) {
        m_monElement[m_indexEta]->Fill(iterColl->eta());
      }
    }
  }
}

Member Data Documentation

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexCharge

private

Definition at line 149 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexEt

private

Definition at line 145 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexEta

private

Definition at line 147 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexEtTotal

private

Definition at line 148 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexHfBitCounts

private

Definition at line 150 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexHfRingEtSums

private

Definition at line 151 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexNrObjects

private

histogram index for each quantity, set during histogram booking

Definition at line 143 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexPhi

private

Definition at line 146 of file L1ExtraDQM.h.

template<class CollectionType >

int L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_indexPt

private

Definition at line 144 of file L1ExtraDQM.h.

template<class CollectionType >

std::vector<MonitorElement*> L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_monElement

private

Definition at line 140 of file L1ExtraDQM.h.

Referenced by L1ExtraDQM::L1ExtraMonElement< CollectionType >::L1ExtraMonElement().

template<class CollectionType >

L1GetHistLimits::Tokens L1ExtraDQM::L1ExtraMonElement< CollectionType >::m_tokens

private

Definition at line 139 of file L1ExtraDQM.h.