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;
}