GctUnpackCollections.cc

Go to the documentation of this file.

#include "EventFilter/GctRawToDigi/src/GctUnpackCollections.h"

GctUnpackCollections::GctUnpackCollections(edm::Event& event)
    : m_event(event),
      m_gctFibres(new L1GctFibreCollection()),  // GCT input collections
      m_rctEm(new L1CaloEmCollection()),
      m_rctCalo(new L1CaloRegionCollection()),
      m_gctInternEm(new L1GctInternEmCandCollection()),  // GCT internal collections
      m_gctInternJets(new L1GctInternJetDataCollection()),
      m_gctInternEtSums(new L1GctInternEtSumCollection()),
      m_gctInternHFData(new L1GctInternHFDataCollection()),
      m_gctInternHtMiss(new L1GctInternHtMissCollection()),
      m_gctIsoEm(new L1GctEmCandCollection()),  // GCT output collections
      m_gctNonIsoEm(new L1GctEmCandCollection()),
      m_gctCenJets(new L1GctJetCandCollection()),
      m_gctForJets(new L1GctJetCandCollection()),
      m_gctTauJets(new L1GctJetCandCollection()),
      m_gctHfBitCounts(new L1GctHFBitCountsCollection()),
      m_gctHfRingEtSums(new L1GctHFRingEtSumsCollection()),
      m_gctEtTot(new L1GctEtTotalCollection()),
      m_gctEtHad(new L1GctEtHadCollection()),
      m_gctEtMiss(new L1GctEtMissCollection()),
      m_gctHtMiss(new L1GctHtMissCollection()),
      m_gctJetCounts(new L1GctJetCountsCollection()),  // Deprecated (empty collection still needed by GT)
      m_errors(new L1TriggerErrorCollection())         // Misc
{
  m_gctIsoEm->reserve(4);
  m_gctCenJets->reserve(4);
  m_gctForJets->reserve(4);
  m_gctTauJets->reserve(4);
  // ** DON'T RESERVE SPACE IN VECTORS FOR DEBUG UNPACK ITEMS! **
}

GctUnpackCollections::~GctUnpackCollections() {
  // GCT input collections
  m_event.put(std::move(m_gctFibres));
  m_event.put(std::move(m_rctEm));
  m_event.put(std::move(m_rctCalo));

  // GCT internal collections
  m_event.put(std::move(m_gctInternEm));
  m_event.put(std::move(m_gctInternJets));
  m_event.put(std::move(m_gctInternEtSums));
  m_event.put(std::move(m_gctInternHFData));
  m_event.put(std::move(m_gctInternHtMiss));

  // GCT output collections
  m_event.put(std::move(m_gctIsoEm), "isoEm");
  m_event.put(std::move(m_gctNonIsoEm), "nonIsoEm");
  m_event.put(std::move(m_gctCenJets), "cenJets");
  m_event.put(std::move(m_gctForJets), "forJets");
  m_event.put(std::move(m_gctTauJets), "tauJets");
  m_event.put(std::move(m_gctHfBitCounts));
  m_event.put(std::move(m_gctHfRingEtSums));
  m_event.put(std::move(m_gctEtTot));
  m_event.put(std::move(m_gctEtHad));
  m_event.put(std::move(m_gctEtMiss));
  m_event.put(std::move(m_gctHtMiss));
  m_event.put(std::move(m_gctJetCounts));  // Deprecated (empty collection still needed by GT)

  // Misc
  m_event.put(std::move(m_errors));
}

std::ostream& operator<<(std::ostream& os, const GctUnpackCollections& rhs) {
  // GCT input collections
  os << "Read " << rhs.gctFibres()->size() << " GCT raw fibre data\n"
     << "Read " << rhs.rctEm()->size() << " RCT EM candidates\n"
     << "Read " << rhs.rctCalo()->size()
     << " RCT Calo Regions\n"

     // GCT internal collections
     << "Read " << rhs.gctInternEm()->size() << " GCT intermediate EM candidates\n"
     << "Read " << rhs.gctInternJets()->size() << " GCT intermediate jet candidates\n"
     << "Read " << rhs.gctInternEtSums()->size() << " GCT intermediate et sums\n"
     << "Read " << rhs.gctInternHFData()->size() << " GCT intermediate HF data\n"
     << "Read " << rhs.gctInternHtMiss()->size()
     << " GCT intermediate Missing Ht\n"

     // GCT output collections
     << "Read " << rhs.gctIsoEm()->size() << " GCT iso EM candidates\n"
     << "Read " << rhs.gctNonIsoEm()->size() << " GCT non-iso EM candidates\n"
     << "Read " << rhs.gctCenJets()->size() << " GCT central jet candidates\n"
     << "Read " << rhs.gctForJets()->size() << " GCT forward jet candidates\n"
     << "Read " << rhs.gctTauJets()->size() << " GCT tau jet candidates\n"
     << "Read " << rhs.gctHfBitCounts()->size() << " GCT HF ring bit counts\n"
     << "Read " << rhs.gctHfRingEtSums()->size() << " GCT HF ring et sums\n"
     << "Read " << rhs.gctEtTot()->size() << " GCT total et\n"
     << "Read " << rhs.gctEtHad()->size() << " GCT ht\n"
     << "Read " << rhs.gctEtMiss()->size() << " GCT met\n"
     << "Read " << rhs.gctHtMiss()->size() << " GCT mht";

  // Any point in putting in an m_errors()->size()? Not sure.

  return os;
}