hgcal::econd::HGCalModuleTreeReader Class Reference

#include <HGCalModuleTreeReader.h>

Inheritance diagram for hgcal::econd::HGCalModuleTreeReader:

Classes
struct	HGCModuleTreeEvent
	Input tree collections. More...

Public Member Functions
	HGCalModuleTreeReader (const EmulatorParameters &, const std::string &tree_name, const std::vector< std::string > &filenames)
ECONDInput	next () override
	Fetch the next ECON-D event. More...
Public Member Functions inherited from hgcal::econd::Emulator
	Emulator (const EmulatorParameters &params)
virtual	~Emulator ()=default

Private Attributes
ECONDInputColl	data_
ECONDInputColl::const_iterator	it_data_

Additional Inherited Members
Protected Attributes inherited from hgcal::econd::Emulator
const EmulatorParameters	params_

Detailed Description

Read out a the relevant raw data produced by a module to memory and returns ECON-D frames on request

Note

The format is as agreed with system tests convenors so that it can be used in integration/beam tests

Definition at line 18 of file HGCalModuleTreeReader.h.

Constructor & Destructor Documentation

HGCalModuleTreeReader()

HGCalModuleTreeReader::HGCalModuleTreeReader ( const EmulatorParameters &  params, const std::string &  tree_name, const std::vector< std::string > &  filenames  )

explicit

Parameters

[in]	tree_name	Name of the TB events tree
[in]	filenames	List of filenames to loop on
[in]	num_channels	Channels multiplicity

Definition at line 20 of file HGCalModuleTreeReader.cc.

References submitPVResolutionJobs::count, data_, edmPickEvents::event, corrVsCorr::filename, amptDefault_cfi::frame, mps_fire::i, it_data_, crabWrapper::key, hgcal::econd::EmulatorParameters::num_channels_per_erx, hgcal::econd::Emulator::params_, and siStripShotFilter_cfi::ZeroSuppressed.

    : Emulator(params) {
  TChain chain(tree_name.data());
  for (const auto& filename : filenames)
    chain.Add(filename.c_str());

  HGCModuleTreeEvent event;
  chain.SetBranchAddress("event", &event.event);
  chain.SetBranchAddress("chip", &event.chip);
  chain.SetBranchAddress("half", &event.half);
  chain.SetBranchAddress("daqdata", &event.daqdata);
  chain.SetBranchAddress("bxcounter", &event.bxcounter);
  chain.SetBranchAddress("eventcounter", &event.eventcounter);
  chain.SetBranchAddress("orbitcounter", &event.orbitcounter);
  chain.SetBranchAddress("trigtime", &event.trigtime);
  chain.SetBranchAddress("trigwidth", &event.trigwidth);

  for (long long i = 0; i < chain.GetEntries(); ++i) {
    chain.GetEntry(i);

    // check if event already exists
    EventId key{(uint32_t)event.eventcounter, (uint32_t)event.bxcounter, (uint32_t)event.orbitcounter};
    if (data_.count(key) == 0)
      data_[key] = ERxInput{};

    // check if chip already exists
    ERxId_t erxKey{(uint8_t)event.chip, (uint8_t)event.half};
    if (data_[key].count(erxKey) == 0)
      data_[key][erxKey] = ERxData{};

    // daqdata: header, CM, 37 ch, CRC32, idle
    if (const auto nwords = event.daqdata->size(); nwords != 41)
      throw cms::Exception("HGCalModuleTreeReader")
          << "Invalid number of words retrieved for event {" << event.eventcounter << ":" << event.bxcounter << ":"
          << event.orbitcounter << "}: should be 41, got " << nwords << ".";

    // 1st word is the header: discard
    // 2nd word are the common mode words

    const uint32_t cmword(event.daqdata->at(1));
    if (((cmword >> 20) & 0xfff) != 0)
      throw cms::Exception("HGCalModuleTreeReader")
          << "Consistency check failed for common mode word: " << ((cmword >> 20) & 0xfff) << " != 0.";

    data_[key][erxKey].cm1 = cmword & 0x3ff;
    data_[key][erxKey].cm0 = (cmword >> 10) & 0x3ff;

    // next 37 words are channel data
    for (size_t i = 2; i < 2 + params_.num_channels_per_erx; i++) {
      HGCROCChannelDataFrame<uint32_t> frame(0, event.daqdata->at(i));
      const auto tctp = static_cast<ToTStatus>(frame.tctp());
      data_[key][erxKey].tctp.push_back(tctp);
      data_[key][erxKey].adcm.push_back(frame.adcm1());
      data_[key][erxKey].adc.push_back(tctp == ToTStatus::ZeroSuppressed ? frame.adc() : 0);
      data_[key][erxKey].tot.push_back(tctp == ToTStatus::ZeroSuppressed ? frame.rawtot() : 0);
      data_[key][erxKey].toa.push_back(frame.toa());
    }

    // copy CRC32
    data_[key][erxKey].crc32 = event.daqdata->at(39);

    // we could assert the idle word from #40 if needed

    // copy metadata
    data_[key][erxKey].meta.push_back(event.trigtime);
    data_[key][erxKey].meta.push_back(event.trigwidth);
  }

  edm::LogInfo("HGCalModuleTreeReader") << "read " << data_.size() << " events.";

  it_data_ = data_.begin();
}

Member Function Documentation

next()

ECONDInput HGCalModuleTreeReader::next ( )

overridevirtual

Fetch the next ECON-D event.

Implements hgcal::econd::Emulator.

Definition at line 96 of file HGCalModuleTreeReader.cc.

References data_, and it_data_.

                                       {
  if (it_data_ == data_.end())
    throw cms::Exception("HGCalModuleTreeReader") << "Insufficient number of events were retrieved from input tree to "
                                                     "proceed with the generation of emulated events.";

  ++it_data_;
  return ECONDInput{it_data_->first, it_data_->second};
}

Member Data Documentation

data_

ECONDInputColl hgcal::econd::HGCalModuleTreeReader::data_

private

Definition at line 36 of file HGCalModuleTreeReader.h.

Referenced by HGCalModuleTreeReader(), and next().

it_data_

ECONDInputColl::const_iterator hgcal::econd::HGCalModuleTreeReader::it_data_

private

Definition at line 37 of file HGCalModuleTreeReader.h.

Referenced by HGCalModuleTreeReader(), and next().