RawToDigiConverter Class Reference

Collection of code to convert TOTEM raw data into digi. More...

#include <RawToDigiConverter.h>

Classes
struct	Record

Public Member Functions
void	printSummaries () const
	Print error summaries. More...
	RawToDigiConverter (const edm::ParameterSet &conf)
void	run (const VFATFrameCollection &coll, const TotemDAQMapping &mapping, const TotemAnalysisMask &mask, edm::DetSetVector< TotemRPDigi > &digi, edm::DetSetVector< TotemVFATStatus > &status)
	Creates RP digi. More...
void	run (const VFATFrameCollection &coll, const TotemDAQMapping &mapping, const TotemAnalysisMask &mask, edm::DetSetVector< CTPPSDiamondDigi > &digi, edm::DetSetVector< TotemVFATStatus > &status)
	Creates Diamond digi. More...
void	run (const VFATFrameCollection &coll, const TotemDAQMapping &mapping, const TotemAnalysisMask &mask, edm::DetSetVector< TotemTimingDigi > &digi, edm::DetSetVector< TotemVFATStatus > &status)
	Creates Totem Timing digi. More...
void	run (const VFATFrameCollection &coll, const TotemDAQMapping &mapping, const TotemAnalysisMask &mask, edmNew::DetSetVector< TotemT2Digi > &digi, edm::DetSetVector< TotemVFATStatus > &status)
	Creates Totem T2 digi. More...

Private Types
enum	TestFlag { tfNoTest, tfWarn, tfErr }

Private Member Functions
void	runCommon (const VFATFrameCollection &input, const TotemDAQMapping &mapping, std::map< TotemFramePosition, Record > &records)
	Common processing for all VFAT based sub-systems. More...

Private Attributes
const double	BC_fraction
const unsigned int	BC_min
const double	EC_fraction
	the minimal required (relative) occupancy of the most frequent counter value to be accepted More...
const unsigned int	EC_min
	the minimal required number of frames to determine the most frequent counter value More...
std::map< TotemFramePosition, std::map< TotemVFATStatus, unsigned int > >	errorSummary
	error summaries More...
const bool	olderTotemT2FileTest
const bool	printErrorSummary
const bool	printUnknownFrameSummary
const unsigned int	testBCMostFrequent
const unsigned int	testCRC
const unsigned int	testECMostFrequent
const unsigned int	testFootprint
	flags for which tests to run More...
const unsigned int	testID
std::map< TotemFramePosition, unsigned int >	unknownSummary
const unsigned char	verbosity

Detailed Description

Collection of code to convert TOTEM raw data into digi.

Definition at line 30 of file RawToDigiConverter.h.

Member Enumeration Documentation

TestFlag

enum RawToDigiConverter::TestFlag

private

Enumerator
tfNoTest
tfWarn
tfErr

Definition at line 77 of file RawToDigiConverter.h.

{ tfNoTest, tfWarn, tfErr };

Constructor & Destructor Documentation

RawToDigiConverter()

RawToDigiConverter::RawToDigiConverter

(

const edm::ParameterSet &

conf

)

Definition at line 26 of file RawToDigiConverter.cc.

    : verbosity(conf.getUntrackedParameter<unsigned int>("verbosity", 0)),
      printErrorSummary(conf.getUntrackedParameter<bool>("printErrorSummary")),
      printUnknownFrameSummary(conf.getUntrackedParameter<bool>("printUnknownFrameSummary")),

      testFootprint(conf.getParameter<unsigned int>("testFootprint")),
      testCRC(conf.getParameter<unsigned int>("testCRC")),
      testID(conf.getParameter<unsigned int>("testID")),
      testECMostFrequent(conf.getParameter<unsigned int>("testECMostFrequent")),
      testBCMostFrequent(conf.getParameter<unsigned int>("testBCMostFrequent")),

      EC_min(conf.getUntrackedParameter<unsigned int>("EC_min", 10)),
      BC_min(conf.getUntrackedParameter<unsigned int>("BC_min", 10)),

      EC_fraction(conf.getUntrackedParameter<double>("EC_fraction", 0.6)),
      BC_fraction(conf.getUntrackedParameter<double>("BC_fraction", 0.6)),

      olderTotemT2FileTest(conf.getParameter<bool>("useOlderT2TestFile")) {}

Member Function Documentation

printSummaries()

void RawToDigiConverter::printSummaries

(

)

const

Print error summaries.

Definition at line 469 of file RawToDigiConverter.cc.

References errorSummary, printErrorSummary, printUnknownFrameSummary, and unknownSummary.

Referenced by TotemVFATRawToDigi::endStream().

                                              {
  // print error summary
  if (printErrorSummary) {
    if (!errorSummary.empty()) {
      stringstream ees;
      for (const auto &vit : errorSummary) {
        ees << vit.first << endl;

        for (const auto &it : vit.second)
          ees << "    " << it.first << " : " << it.second << endl;
      }

      LogWarning("Totem") << "RawToDigiConverter: error summary (error signature : number of such events)\n"
                          << endl
                          << ees.rdbuf();
    } else {
      LogInfo("Totem") << "RawToDigiConverter: no errors to be reported.";
    }
  }

  // print summary of unknown frames (found in data but not in the mapping)
  if (printUnknownFrameSummary) {
    if (!unknownSummary.empty()) {
      stringstream ees;
      for (const auto &it : unknownSummary)
        ees << "  " << it.first << " : " << it.second << endl;

      LogWarning("Totem")
          << "RawToDigiConverter: frames found in data, but not in the mapping (frame position : number of events)\n"
          << endl
          << ees.rdbuf();
    } else {
      LogInfo("Totem") << "RawToDigiConverter: no unknown frames to be reported.";
    }
  }
}

run() [1/4]

void RawToDigiConverter::run	(	const VFATFrameCollection &	coll,
		const TotemDAQMapping &	mapping,
		const TotemAnalysisMask &	mask,
		edm::DetSetVector< TotemRPDigi > &	digi,
		edm::DetSetVector< TotemVFATStatus > &	status
	)

Creates RP digi.

Definition at line 190 of file RawToDigiConverter.cc.

References TotemAnalysisMask::analysisMask, TotemRPDetId::chip(), hcalRecHitTable_cff::detId, edm::DetSet< T >::emplace_back(), edm::DetSetVector< T >::find_or_insert(), TotemVFATAnalysisMask::fullMask, input, taus_updatedMVAIds_cff::mapping, TotemVFATAnalysisMask::maskedChannels, hltrates_dqm_sourceclient-live_cfg::offset, AlCaHLTBitMon_ParallelJobs::p, TotemRPDetId::planeId(), edm::DetSet< T >::push_back(), AlCaHarvesting_cff::record, and runCommon().

Referenced by Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), and TotemVFATRawToDigi::run().

                                                                         {
  // structure merging vfat frame data with the mapping
  map<TotemFramePosition, Record> records;

  // common processing - frame validation
  runCommon(input, mapping, records);

  // second loop over data
  for (auto &p : records) {
    Record &record = p.second;

    // calculate ids
    TotemRPDetId chipId(record.info->symbolicID.symbolicID);
    uint8_t chipPosition = chipId.chip();
    TotemRPDetId detId = chipId.planeId();

    // update chipPosition in status
    record.status.setChipPosition(chipPosition);

    // produce digi only for good frames
    if (record.status.isOK()) {
      // find analysis mask (needs a default=no mask, if not in present the mapping)
      TotemVFATAnalysisMask anMa;
      anMa.fullMask = false;

      auto analysisIter = analysisMask.analysisMask.find(record.info->symbolicID);
      if (analysisIter != analysisMask.analysisMask.end()) {
        // if there is some information about masked channels - save it into conversionStatus
        anMa = analysisIter->second;
        if (anMa.fullMask)
          record.status.setFullyMaskedOut();
        else
          record.status.setPartiallyMaskedOut();
      }

      // create the digi
      unsigned short offset = chipPosition * 128;
      const vector<unsigned char> &activeChannels = record.frame->getActiveChannels();

      for (auto ch : activeChannels) {
        // skip masked channels
        if (!anMa.fullMask && anMa.maskedChannels.find(ch) == anMa.maskedChannels.end()) {
          DetSet<TotemRPDigi> &digiDetSet = rpData.find_or_insert(detId);
          digiDetSet.emplace_back(offset + ch);
        }
      }
    }

    // save status
    DetSet<TotemVFATStatus> &statusDetSet = finalStatus.find_or_insert(detId);
    statusDetSet.push_back(record.status);
  }
}

run() [2/4]

void RawToDigiConverter::run	(	const VFATFrameCollection &	coll,
		const TotemDAQMapping &	mapping,
		const TotemAnalysisMask &	mask,
		edm::DetSetVector< CTPPSDiamondDigi > &	digi,
		edm::DetSetVector< TotemVFATStatus > &	status
	)

Creates Diamond digi.

Definition at line 248 of file RawToDigiConverter.cc.

References hcalRecHitTable_cff::detId, edm::DetSet< T >::emplace_back(), edm::DetSetVector< T >::find_or_insert(), pps::diamond::vfat::getHptdcErrorFlag(), pps::diamond::vfat::getLeadingEdgeTime(), pps::diamond::vfat::getMultihit(), pps::diamond::vfat::getThresholdVoltage(), pps::diamond::vfat::getTrailingEdgeTime(), taus_updatedMVAIds_cff::mapping, AlCaHLTBitMon_ParallelJobs::p, edm::DetSet< T >::push_back(), AlCaHarvesting_cff::record, runCommon(), and mps_update::status.

Referenced by Types.EventID::cppID(), and Types.LuminosityBlockID::cppID().

                                                                       {
  // structure merging vfat frame data with the mapping
  map<TotemFramePosition, Record> records;

  // common processing - frame validation
  runCommon(coll, mapping, records);

  // second loop over data
  for (auto &p : records) {
    Record &record = p.second;

    // calculate ids
    CTPPSDiamondDetId detId(record.info->symbolicID.symbolicID);

    if (record.status.isOK()) {
      // update Event Counter in status
      record.status.setEC(record.frame->getEC() & 0xFF);

      // create the digi
      DetSet<CTPPSDiamondDigi> &digiDetSet = digi.find_or_insert(detId);
      digiDetSet.emplace_back(pps::diamond::vfat::getLeadingEdgeTime(*record.frame),
                              pps::diamond::vfat::getTrailingEdgeTime(*record.frame),
                              pps::diamond::vfat::getThresholdVoltage(*record.frame),
                              pps::diamond::vfat::getMultihit(*record.frame),
                              pps::diamond::vfat::getHptdcErrorFlag(*record.frame));
    }

    // save status
    DetSet<TotemVFATStatus> &statusDetSet = status.find_or_insert(detId);
    statusDetSet.push_back(record.status);
  }
}

run() [3/4]

void RawToDigiConverter::run	(	const VFATFrameCollection &	coll,
		const TotemDAQMapping &	mapping,
		const TotemAnalysisMask &	mask,
		edm::DetSetVector< TotemTimingDigi > &	digi,
		edm::DetSetVector< TotemVFATStatus > &	status
	)

Creates Totem Timing digi.

Definition at line 285 of file RawToDigiConverter.cc.

References TotemDAQMapping::TotemTimingPlaneChannelPair::channel, TauDecayModes::dec, hcalRecHitTable_cff::detId, dumpTauVariables_cfi::eventInfo, edm::DetSetVector< T >::find_or_insert(), RawToDigiConverter::Record::frame, VFATFrame::getData(), TotemFramePosition::getGOHId(), TotemFramePosition::getIdxInFiber(), TotemFramePosition::getOptoRxId(), TotemFramePosition::getSubSystemId(), TotemFramePosition::getTOTFEDId(), createfilelist::int, taus_updatedMVAIds_cff::mapping, AlCaHLTBitMon_ParallelJobs::p, TotemDAQMapping::TotemTimingPlaneChannelPair::plane, edm::DetSet< T >::push_back(), AlCaHarvesting_cff::record, runCommon(), and verbosity.

Referenced by Types.EventID::cppID(), and Types.LuminosityBlockID::cppID().

                                                                       {
  // structure merging vfat frame data with the mapping
  map<TotemFramePosition, Record> records;

  // common processing - frame validation
  runCommon(coll, mapping, records);

  // second loop over data
  for (auto &p : records) {
    Record &record = p.second;
    if (!record.status.isOK())
      continue;

    const TotemFramePosition *framepos = &p.first;

    if (((framepos->getIdxInFiber() % 2) == 0) && (framepos->getIdxInFiber() < 14)) {
      //corresponding channel data are always in the neighbouring idx in fiber

      TotemFramePosition frameposdata(framepos->getSubSystemId(),
                                      framepos->getTOTFEDId(),
                                      framepos->getOptoRxId(),
                                      framepos->getGOHId(),
                                      (framepos->getIdxInFiber() + 1));
      TotemFramePosition frameposEvtInfo(
          framepos->getSubSystemId(), framepos->getTOTFEDId(), framepos->getOptoRxId(), framepos->getGOHId(), 0xe);

      auto channelwaveformPtr = records.find(frameposdata);
      auto eventInfoPtr = records.find(frameposEvtInfo);

      if (channelwaveformPtr != records.end() && eventInfoPtr != records.end()) {
        Record &channelwaveform = records[frameposdata];
        Record &eventInfo = records[frameposEvtInfo];

        // Extract all the waveform information from the raw data
        TotemSampicFrame totemSampicFrame((const uint8_t *)record.frame->getData(),
                                          (const uint8_t *)channelwaveform.frame->getData(),
                                          (const uint8_t *)eventInfo.frame->getData());

        if (totemSampicFrame.valid()) {
          // create the digi
          TotemTimingEventInfo eventInfoTmp(totemSampicFrame.getEventHardwareId(),
                                            totemSampicFrame.getL1ATimestamp(),
                                            totemSampicFrame.getBunchNumber(),
                                            totemSampicFrame.getOrbitNumber(),
                                            totemSampicFrame.getEventNumber(),
                                            totemSampicFrame.getChannelMap(),
                                            totemSampicFrame.getL1ALatency(),
                                            totemSampicFrame.getNumberOfSentSamples(),
                                            totemSampicFrame.getOffsetOfSamples(),
                                            totemSampicFrame.getPLLInfo());
          TotemTimingDigi digiTmp(totemSampicFrame.getHardwareId(),
                                  totemSampicFrame.getFPGATimestamp(),
                                  totemSampicFrame.getTimestampA(),
                                  totemSampicFrame.getTimestampB(),
                                  totemSampicFrame.getCellInfo(),
                                  totemSampicFrame.getSamples(),
                                  eventInfoTmp);
          // calculate ids
          TotemTimingDetId detId(record.info->symbolicID.symbolicID);
          const TotemDAQMapping::TotemTimingPlaneChannelPair SWpair =
              mapping.getTimingChannel(totemSampicFrame.getHardwareId());
          // for FW Version > 0 plane and channel are encoded in the dataframe
          if (totemSampicFrame.getFWVersion() == 0)  // Mapping not present in HW, read from SW for FW versions == 0
          {
            if (SWpair.plane == -1 || SWpair.channel == -1) {
              if (verbosity > 0)
                LogWarning("Totem") << "Error in RawToDigiConverter::TotemTiming > "
                                    << "HwId not recognized!  hwId: " << std::hex
                                    << (unsigned int)totemSampicFrame.getHardwareId() << endl;
            } else {
              detId.setPlane(SWpair.plane % 4);
              detId.setChannel(SWpair.channel);
            }
          } else  // Mapping read from HW, checked by SW
          {
            const int HWplane = totemSampicFrame.getDetPlane() % 16;
            const int HWchannel = totemSampicFrame.getDetChannel() % 16;

            if (SWpair.plane == -1 || SWpair.channel == -1) {
              if (verbosity > 0)
                LogWarning("Totem") << "Warning in RawToDigiConverter::TotemTiming > "
                                    << "HwId not recognized!  hwId: " << std::hex
                                    << (unsigned int)totemSampicFrame.getHardwareId()
                                    << "\tUsing plane and ch from HW without check!" << endl;
            } else {
              if (verbosity > 0 && (SWpair.plane != HWplane || SWpair.channel != HWchannel))
                LogWarning("Totem") << "Warning in RawToDigiConverter::TotemTiming > "
                                    << "Hw mapping different from SW mapping. hwId: " << std::hex
                                    << (unsigned int)totemSampicFrame.getHardwareId() << "HW: " << std::dec << HWplane
                                    << ":" << HWchannel << "\tSW " << SWpair.plane << ":" << SWpair.channel
                                    << "\tUsing plane and ch from HW!" << endl;
            }
            detId.setPlane(HWplane % 4);
            detId.setChannel(HWchannel);
          }

          DetSet<TotemTimingDigi> &digiDetSet = digi.find_or_insert(detId);
          digiDetSet.push_back(digiTmp);
        }
      }
    }
  }
}

run() [4/4]

void RawToDigiConverter::run	(	const VFATFrameCollection &	coll,
		const TotemDAQMapping &	mapping,
		const TotemAnalysisMask &	mask,
		edmNew::DetSetVector< TotemT2Digi > &	digi,
		edm::DetSetVector< TotemVFATStatus > &	status
	)

Creates Totem T2 digi.

Definition at line 393 of file RawToDigiConverter.cc.

References totem::nt2::vfat::channelMarker(), hcalRecHitTable_cff::detId, edmNew::DetSetVector< T >::FastFiller, createfilelist::int, totem::nt2::vfat::leadingEdgeTime(), taus_updatedMVAIds_cff::mapping, totem::nt2::vfat::newChannelId(), totem::nt2::vfat::num_channels_per_payload, olderTotemT2FileTest, AlCaHLTBitMon_ParallelJobs::p, edm::DetSet< T >::push_back(), AlCaHarvesting_cff::record, runCommon(), mps_update::status, totem::nt2::vfat::statusMarker(), totem::nt2::vfat::trailingEdgeTime(), verbosity, and y.

Referenced by Types.EventID::cppID(), and Types.LuminosityBlockID::cppID().

                                                                       {
  // structure merging vfat frame data with the mapping
  map<TotemFramePosition, Record> records;

  // common processing - frame validation
  runCommon(coll, mapping, records);

  int allT2 = 0;
  int goodT2 = 0;
  int foundT2 = 0;
  const int T2shiftOld = (olderTotemT2FileTest ? 8 : 0);  //Run on TOTEM T2 test file (ver 2.1) or final T2 data ver 2.3

  // second loop over data
  for (auto &p : records) {
    Record &record = p.second;

    allT2++;
    // calculate ids
    TotemT2DetId detId(record.info->symbolicID.symbolicID);

    if (record.status.isOK()) {
      // update Event Counter in status
      record.status.setEC(record.frame->getEC() & 0xFF);
      goodT2++;
      if (verbosity > 2) {
        LogWarning("Totem") << "RawToDigiConverter: VFAT frame number " << allT2
                            << " is OK , mapping HW_ID (decimal) is: " << (record.info->hwID)
                            << ", T2DetId arm/plane/channel = " << (detId) << endl;
        LogWarning("Totem") << "HW_id_16b CH0 (dec), LE CH0, TE CH0, marker CH0, HW_id_16b CH1 (dec), LE CH1,"
                            << " TE CH1, marker CH1 = ";
        for (size_t y = 0; y < 2; y++) {
          LogWarning("Totem") << ((unsigned int)totem::nt2::vfat::newChannelId(*record.frame, y)) << "/"
                              << ((unsigned int)totem::nt2::vfat::leadingEdgeTime(*record.frame, y)) << "/"
                              << ((unsigned int)totem::nt2::vfat::trailingEdgeTime(*record.frame, y)) << "/"
                              << ((unsigned int)totem::nt2::vfat::channelMarker(*record.frame, y)) << "/";
        }
      }
      for (size_t frame_id = 0; frame_id < totem::nt2::vfat::num_channels_per_payload; ++frame_id) {
        if (const uint16_t hw_id = totem::nt2::vfat::newChannelId(*record.frame, frame_id) >> T2shiftOld;
            hw_id == record.info->hwID) {  // only unpack the payload associated to this hardware ID
          // create the digi
          edmNew::DetSetVector<TotemT2Digi>::FastFiller(digi, detId)
              .emplace_back(hw_id,
                            totem::nt2::vfat::channelMarker(*record.frame, frame_id),
                            totem::nt2::vfat::leadingEdgeTime(*record.frame, frame_id),
                            totem::nt2::vfat::trailingEdgeTime(*record.frame, frame_id),
                            totem::nt2::vfat::statusMarker(*record.frame));
          foundT2++;
        } else {
          if (verbosity > 2)
            LogWarning("Totem") << "HW_ID comparison fail (CH#/Channel HwID/Mapping HwID): " << ((int)frame_id) << "/"
                                << ((unsigned int)hw_id) << "/" << (record.info->hwID) << endl;
        }
      }
    } else {
      if (verbosity > 1)
        LogWarning("Totem") << "Bad T2 record, is missing/IDmismatch/footprintError"
                            << "/CRCerror/ECprogressBad/BCprogressBad: " << record.status.isMissing() << "/"
                            << record.status.isIDMismatch() << "/" << record.status.isFootprintError() << "/"
                            << record.status.isCRCError() << "/" << record.status.isECProgressError() << "/"
                            << record.status.isBCProgressError() << "/" << endl;
    }

    // save status
    DetSet<TotemVFATStatus> &statusDetSet = status.find_or_insert(detId);
    statusDetSet.push_back(record.status);
  }
  if (verbosity > 1)
    LogWarning("Totem") << "RawToDigiConverter:: VFAT frames per event, total/good/matched the xml mapping"
                        << " (T2Digi created): " << allT2 << "/" << goodT2 << "/" << foundT2 << endl;
}

runCommon()

void RawToDigiConverter::runCommon ( const VFATFrameCollection &  input, const TotemDAQMapping &  mapping, std::map< TotemFramePosition, Record > &  records  )

private

Common processing for all VFAT based sub-systems.

Definition at line 45 of file RawToDigiConverter.cc.

References CounterChecker::Analyze(), BC_fraction, BC_min, CounterChecker::BCChecker, TauDecayModes::dec, EC_fraction, EC_min, CounterChecker::ECChecker, errorSummary, CounterChecker::Fill(), input, VFATFrameCollection::Iterator::IsEnd(), visualization-live-secondInstance_cfg::m, taus_updatedMVAIds_cff::mapping, AlCaHLTBitMon_ParallelJobs::p, printErrorSummary, AlCaHarvesting_cff::record, TotemVFATStatus::setMissing(), testBCMostFrequent, testCRC, testECMostFrequent, testFootprint, testID, tfErr, tfNoTest, unknownSummary, and verbosity.

Referenced by run().

                                                                                               {
  // EC and BC checks (wrt. the most frequent value), BC checks per subsystem
  CounterChecker ECChecker(CounterChecker::ECChecker, "EC", EC_min, EC_fraction, verbosity);
  CounterChecker BCChecker(CounterChecker::BCChecker, "BC", BC_min, BC_fraction, verbosity);

  // initialise structure merging vfat frame data with the mapping
  for (auto &p : mapping.VFATMapping) {
    TotemVFATStatus st;
    st.setMissing(true);
    records[p.first] = {&p.second, nullptr, st};
  }

  // event error message buffer
  stringstream ees;

  // associate data frames with records
  for (VFATFrameCollection::Iterator fr(&input); !fr.IsEnd(); fr.Next()) {
    // frame error message buffer
    stringstream fes;

    bool problemsPresent = false;
    bool stopProcessing = false;

    // skip data frames not listed in the DAQ mapping
    auto records_it = records.find(fr.Position());
    if (records_it == records.end()) {
      unknownSummary[fr.Position()]++;
      continue;
    }

    // update record
    Record &record = records_it->second;
    record.frame = fr.Data();
    record.status.setMissing(false);
    record.status.setNumberOfClustersSpecified(record.frame->isNumberOfClustersPresent());
    record.status.setNumberOfClusters(record.frame->getNumberOfClusters());

    // check for T2 payload bits
    int rawT2 = fr.Position().getRawPosition();
    bool isT2Frame = (rawT2 >> 18);

    // check footprint
    if (((!isT2Frame) && testFootprint != tfNoTest && !record.frame->checkFootprint()) ||
        (isT2Frame && testFootprint != tfNoTest && !record.frame->checkFootprintT2())) {
      problemsPresent = true;

      if (verbosity > 0)
        fes << "    invalid footprint" << endl;

      if (testFootprint == tfErr) {
        record.status.setFootprintError();
        stopProcessing = true;
      }
    }

    // check CRC
    if (((!isT2Frame) && (testCRC != tfNoTest && !record.frame->checkCRC())) ||
        (isT2Frame && testCRC != tfNoTest && !record.frame->checkCRCT2())) {
      problemsPresent = true;

      if (verbosity > 0)
        fes << "    CRC failure" << endl;

      if (testCRC == tfErr) {
        record.status.setCRCError();
        stopProcessing = true;
      }
    }

    // check the id mismatch
    if (testID != tfNoTest && record.frame->isIDPresent() &&
        (record.frame->getChipID() & 0xFFF) != (record.info->hwID & 0xFFF)) {
      if (verbosity > 0)
        fes << "    ID mismatch (data: 0x" << hex << record.frame->getChipID() << ", mapping: 0x" << record.info->hwID
            << dec << ", symbId: " << record.info->symbolicID.symbolicID << ")" << endl;

      if (testID == tfErr) {
        record.status.setIDMismatch();
        stopProcessing = true;
      }
    }

    // if there were errors, put the information to ees buffer
    if (verbosity > 0 && problemsPresent) {
      string message = (stopProcessing) ? "(and will be dropped)" : "(but will be used though)";
      if (verbosity > 2) {
        if (isT2Frame && verbosity > 3)
          record.frame->PrintT2();
        ees << "  Frame at " << fr.Position() << " seems corrupted " << message << ":" << endl;
        ees << fes.rdbuf();
      } else
        ees << "  Frame at " << fr.Position() << " seems corrupted " << message << "." << endl;
    }

    // if there were serious errors, do not process this frame
    if (stopProcessing)
      continue;

    // fill EC and BC values to the statistics
    if (fr.Data()->isECPresent())
      ECChecker.Fill(fr.Data()->getEC(), fr.Position());

    if (fr.Data()->isBCPresent())
      BCChecker.Fill(fr.Data()->getBC(), fr.Position());
  }

  // analyze EC and BC statistics
  if (testECMostFrequent != tfNoTest)
    ECChecker.Analyze(records, (testECMostFrequent == tfErr), ees);

  if (testBCMostFrequent != tfNoTest)
    BCChecker.Analyze(records, (testBCMostFrequent == tfErr), ees);

  // add error message for missing frames
  if (verbosity > 1) {
    for (const auto &p : records) {
      if (p.second.status.isMissing())
        ees << "Frame for VFAT " << p.first << " is not present in the data." << endl;
    }
  }

  // print error message
  if (verbosity > 0 && !ees.rdbuf()->str().empty()) {
    if (verbosity > 1)
      LogWarning("Totem") << "Error in RawToDigiConverter::runCommon > "
                          << "event contains the following problems:" << endl
                          << ees.rdbuf() << endl;
    else
      LogWarning("Totem") << "Error in RawToDigiConverter::runCommon > "
                          << "event contains problems." << endl;
  }

  // increase error counters
  if (printErrorSummary) {
    for (const auto &it : records) {
      if (!it.second.status.isOK()) {
        auto &m = errorSummary[it.first];
        m[it.second.status]++;
      }
    }
  }
}

Member Data Documentation

BC_fraction

const double RawToDigiConverter::BC_fraction

private

Definition at line 90 of file RawToDigiConverter.h.

Referenced by runCommon().

BC_min

const unsigned int RawToDigiConverter::BC_min

private

Definition at line 87 of file RawToDigiConverter.h.

Referenced by runCommon().

EC_fraction

const double RawToDigiConverter::EC_fraction

private

the minimal required (relative) occupancy of the most frequent counter value to be accepted

Definition at line 90 of file RawToDigiConverter.h.

Referenced by runCommon().

EC_min

const unsigned int RawToDigiConverter::EC_min

private

the minimal required number of frames to determine the most frequent counter value

Definition at line 87 of file RawToDigiConverter.h.

Referenced by runCommon().

errorSummary

std::map<TotemFramePosition, std::map<TotemVFATStatus, unsigned int> > RawToDigiConverter::errorSummary

private

error summaries

Definition at line 96 of file RawToDigiConverter.h.

Referenced by printSummaries(), and runCommon().

olderTotemT2FileTest

const bool RawToDigiConverter::olderTotemT2FileTest

private

Definition at line 93 of file RawToDigiConverter.h.

Referenced by run().

printErrorSummary

const bool RawToDigiConverter::printErrorSummary

private

Definition at line 74 of file RawToDigiConverter.h.

Referenced by printSummaries(), and runCommon().

printUnknownFrameSummary

const bool RawToDigiConverter::printUnknownFrameSummary

private

Definition at line 75 of file RawToDigiConverter.h.

Referenced by printSummaries().

testBCMostFrequent

const unsigned int RawToDigiConverter::testBCMostFrequent

private

Definition at line 84 of file RawToDigiConverter.h.

Referenced by runCommon().

testCRC

const unsigned int RawToDigiConverter::testCRC

private

Definition at line 81 of file RawToDigiConverter.h.

Referenced by runCommon().

testECMostFrequent

const unsigned int RawToDigiConverter::testECMostFrequent

private

Definition at line 83 of file RawToDigiConverter.h.

Referenced by runCommon().

testFootprint

const unsigned int RawToDigiConverter::testFootprint

private

flags for which tests to run

Definition at line 80 of file RawToDigiConverter.h.

Referenced by runCommon().

testID

const unsigned int RawToDigiConverter::testID

private

Definition at line 82 of file RawToDigiConverter.h.

Referenced by runCommon().

unknownSummary

std::map<TotemFramePosition, unsigned int> RawToDigiConverter::unknownSummary

private

Definition at line 97 of file RawToDigiConverter.h.

Referenced by printSummaries(), and runCommon().

verbosity

const unsigned char RawToDigiConverter::verbosity

private

Definition at line 72 of file RawToDigiConverter.h.

Referenced by run(), and runCommon().