RctRawToDigi Class Reference

#include <EventFilter/RctRawToDigi/src/RctRawToDigi.cc>

Inheritance diagram for RctRawToDigi:

Public Member Functions
	RctRawToDigi (const edm::ParameterSet &)
	~RctRawToDigi () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	checkHeaders ()
bool	decodeLinkID (const uint32_t inputValue, uint32_t &crateNumber, uint32_t &linkNumber, bool &even)
virtual void	endJob ()
	method called at job end - use to print summary report More...
bool	printAll (const unsigned char *data, const unsigned size)
void	produce (edm::Event &, const edm::EventSetup &) override
void	unpack (const FEDRawData &d, edm::Event &e, RctUnpackCollections *const colls)
	Unpacks the raw data. More...
void	unpackCTP7 (const uint32_t *data, const unsigned block_id, const unsigned size, RctUnpackCollections *const colls)

Private Attributes
std::vector< unsigned >	errorCounters_
	Counts number of errors for each code (index) More...
L1TriggerErrorCollection *	errors_
	pointer to error collection More...
int	fedId_
	RCT FED ID. More...
edm::InputTag	inputLabel_
	FED collection label. More...
unsigned	unpackFailures_
	To count the total number of RCT unpack failures. More...
const bool	verbose_
	If true, then debug print out for each event. More...

Static Private Attributes
static const unsigned	amc13HeaderSize_ = 16
static const unsigned	amc13TrailerSize_ = 8
static const unsigned	ctp7HeaderSize_ = 8
static const unsigned	ctp7TrailerSize_ = 8
static const unsigned	MAX_DATA = 4680
	The maximum number of blocks we will try to unpack before thinking something is wrong. More...
static const unsigned	MAX_ERR_CODE = 6
static const unsigned	MIN_DATA = 900
	The minimum number of blocks we will try to unpack before thinking something is wrong (really this should be 920, to be tested) More...
static const unsigned	sLinkHeaderSize_ = 8
static const unsigned	sLinkTrailerSize_ = 8

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: Produce RCT digis from raw data

Implementation: <Notes on="" implementation>="">

Definition at line 40 of file RctRawToDigi.h.

Constructor & Destructor Documentation

RctRawToDigi()

RctRawToDigi::RctRawToDigi ( const edm::ParameterSet &  iConfig )

explicit

Register Products

Definition at line 34 of file RctRawToDigi.cc.

References fedId_, inputLabel_, and LogDebug.

    : inputLabel_(iConfig.getParameter<edm::InputTag>("inputLabel")),
      fedId_(iConfig.getUntrackedParameter<int>("rctFedId", FEDNumbering::MINRCTFEDID)),
      verbose_(iConfig.getUntrackedParameter<bool>("verbose", false)) {
  LogDebug("RCT") << "RctRawToDigi will unpack FED Id " << fedId_;
  // RCT collections
  produces<L1CaloEmCollection>();
  produces<L1CaloRegionCollection>();

  // Error collection
  consumes<FEDRawDataCollection>(inputLabel_);
}

~RctRawToDigi()

RctRawToDigi::~RctRawToDigi ( )

override

Definition at line 48 of file RctRawToDigi.cc.

                            {
  // do anything here that needs to be done at destruction time
  // (e.g. close files, deallocate resources etc.)
  //delete formatTranslator_;
}

Member Function Documentation

checkHeaders()

void RctRawToDigi::checkHeaders ( )

private

Looks at the firmware version header in the S-Link packet and instantiates relevant format translator. check block headers for consistency

decodeLinkID()

bool RctRawToDigi::decodeLinkID ( const uint32_t  inputValue, uint32_t &  crateNumber, uint32_t &  linkNumber, bool &  even  )

private

endJob()

void RctRawToDigi::endJob ( void  )

privatevirtual

method called at job end - use to print summary report

Definition at line 311 of file RctRawToDigi.cc.

References errorCounters_, mps_fire::i, MAX_ERR_CODE, l1tGTMenu_BTagSeeds_cff::os, dqmMemoryStats::total, and verbose_.

                          {
  unsigned total = 0;
  std::ostringstream os;

  for (unsigned i = 0; i <= MAX_ERR_CODE; ++i) {
    total += errorCounters_.at(i);
    os << "Error " << i << " (" << errorCounters_.at(i) << ")";
    if (i < MAX_ERR_CODE) {
      os << ", ";
    }
  }

  if (total > 0 && verbose_) {
    edm::LogError("RCT") << "Encountered " << total << " unpacking errors: " << os.str();
  }
}

printAll()

bool RctRawToDigi::printAll ( const unsigned char *  data, const unsigned  size  )

private

Definition at line 301 of file RctRawToDigi.cc.

References gather_cfg::cout, data, mps_fire::i, and findQualityFiles::size.

                                                                          {
  for (unsigned i = 0; i < size; i++) {
    std::cout << data[i] << " ";
    if (i % 6 == 5)
      std::cout << std::endl;
  }
  return true;
}

produce()

void RctRawToDigi::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 55 of file RctRawToDigi.cc.

References amc13HeaderSize_, amc13TrailerSize_, TauDecayModes::dec, fedId_, l1tstage2emulator_dqm_sourceclient-live_cfg::feds, iEvent, inputLabel_, LogDebug, MIN_DATA, FEDRawData::size(), sLinkHeaderSize_, sLinkTrailerSize_, and unpack().

                                                                        {
  using namespace edm;

  // Instantiate all the collections the unpacker needs; puts them in event when this object goes out of scope.
  std::unique_ptr<RctUnpackCollections> colls(new RctUnpackCollections(iEvent));

  // get raw data collection
  edm::Handle<FEDRawDataCollection> feds;
  iEvent.getByLabel(inputLabel_, feds);

  // if raw data collection is present, check the headers and do the unpacking
  if (feds.isValid()) {
    const FEDRawData& rctRcd = feds->FEDData(fedId_);

    //Check FED size
    LogDebug("RCT") << "Upacking FEDRawData of size " << std::dec << rctRcd.size();

    //check header size
    if (rctRcd.size() < sLinkHeaderSize_ + sLinkTrailerSize_ + amc13HeaderSize_ + amc13TrailerSize_ + MIN_DATA) {
      if (rctRcd.size() > 0) {
        LogError("L1T") << "Cannot unpack: empty/invalid L1T raw data (size = " << rctRcd.size() << ") for ID "
                        << fedId_ << ". Returning empty collections!";
      } else {
        // there's no MC packer for this payload, so totally expected that sometimes it will be absent... no warning issued.
      }
      //continue;
      return;
    }

    // do the unpacking
    unpack(rctRcd, iEvent, colls.get());

  } else {
    LogError("L1T") << "Cannot unpack: no collection found";

    return;
  }
}

unpack()

void RctRawToDigi::unpack ( const FEDRawData &  d, edm::Event &  e, RctUnpackCollections *const  colls  )

private

Unpacks the raw data.

Parameters

invalidDataFlag

- if true, then won't attempt unpack but just output empty collecions.

Definition at line 94 of file RctRawToDigi.cc.

References FEDTrailer::check(), FEDTrailer::crc(), ztail::d, data, FEDTrailer::evtStatus(), FEDTrailer::fragmentLength(), RecoTauValidation_cfi::header, LogDebug, sLinkTrailerSize_, FEDTrailer::ttsBits(), and unpackCTP7().

Referenced by produce().

                                                                                             {
  const unsigned char* data = d.data();  // The 8-bit wide raw-data array.
  // Data offset - starts at 16 as there is a 64-bit S-Link header followed
  // by a 64-bit software-controlled header (for pipeline format version
  // info that is not yet used).

  FEDHeader header(data);

  if (header.check()) {
    LogDebug("L1T") << "Found SLink header:"
                    << " Trigger type " << header.triggerType() << " L1 event ID " << header.lvl1ID() << " BX Number "
                    << header.bxID() << " FED source " << header.sourceID() << " FED version " << header.version();
  } else {
    LogWarning("L1T") << "Did not find a valid SLink header!";
  }

  FEDTrailer trailer(data + (d.size() - sLinkTrailerSize_));

  if (trailer.check()) {
    LogDebug("L1T") << "Found SLink trailer:"
                    << " Length " << trailer.fragmentLength() << " CRC " << trailer.crc() << " Status "
                    << trailer.evtStatus() << " Throttling bits " << trailer.ttsBits();
  } else {
    LogWarning("L1T") << "Did not find a SLink trailer!";
  }

  unpackCTP7(reinterpret_cast<const uint32_t*>(data), 0, sizeof(data), colls);
}

unpackCTP7()

void RctRawToDigi::unpackCTP7 ( const uint32_t *  data, const unsigned  block_id, const unsigned  size, RctUnpackCollections *const  colls  )

private

Definition at line 123 of file RctRawToDigi.cc.

References cms::cuda::assert(), data, dqmiolumiharvest::j, isotrackApplyRegressor::k, LogDebug, visualization-live-secondInstance_cfg::m, EcalTangentSkim_cfg::o, haddnano::of, HLT_IsoTrack_cff::offset, submitPVResolutionJobs::q, RctUnpackCollections::rctCalo(), RctUnpackCollections::rctEm(), L1CaloEmCand::setBx(), L1CaloRegion::setBx(), findQualityFiles::size, submitPVValidationJobs::t, createJobs::tmp, parallelization::uint, and testProducerWithPsetDescEmpty_cfi::x1.

Referenced by unpack().

                                                                 {
  //offset from 6 link header words
  uint32_t of = 6;
  LogDebug("L1T") << "Block ID  = " << block_id << " size = " << size;

  CTP7Format ctp7Format;
  RctDataDecoder rctDataDecoder;
  uint32_t nBXTemp = 0;
  uint32_t ctp7FWVersion;
  uint32_t L1ID, L1aBCID;
  std::vector<RCTInfo> allCrateRCTInfo[5];

  L1ID = data[1 + of];                          // extract the L1 ID number
  L1aBCID = data[5 + of] & 0x00000FFF;          // extract the BCID number of L1A
  nBXTemp = (data[5 + of] & 0x00FF0000) >> 16;  // extract number of BXs readout per L1A
  ctp7FWVersion = data[4 + of];

  if (nBXTemp != 1 && nBXTemp != 3 && nBXTemp != 5)
    nBXTemp = 1;

  const uint32_t nBX = nBXTemp;

  LogDebug("L1T") << "CTP7 L1ID = " << L1ID << " L1A BCID = " << L1aBCID << " BXs in capture = " << nBX
                  << " CTP7 DAQ FW = " << ctp7FWVersion;

  struct link_data {
    bool even;
    unsigned int crateID;
    unsigned int ctp7LinkNumber;
    std::vector<unsigned int> uint;
  };

  //nBX max 5, nLinks max 36 [nBX][nLinks]
  link_data allLinks[5][36];
  const uint32_t NLinks = ctp7Format.NLINKS;
  assert(NLinks <= 36);

  //change this implementation
  uint32_t iDAQBuffer = 0;

  //Step 1: Grab all data from ctp7 buffer and put into link format
  for (unsigned int iLink = 0; iLink < NLinks; iLink++) {
    iDAQBuffer = of + ctp7Format.EVENT_HEADER_WORDS +
                 iLink * (ctp7Format.CHANNEL_HEADER_WORDS + nBX * ctp7Format.CHANNEL_DATA_WORDS_PER_BX);

    //first decode linkID
    uint32_t linkID = data[iDAQBuffer++];
    uint32_t tmp = data[iDAQBuffer++];
    uint32_t CRCErrCnt = tmp & 0x0000FFFF;
    //uint32_t linkStatus = (tmp & 0xFFFF0000) >> 16;

    uint32_t crateID = 0;
    uint32_t expectedCrateID = 0;
    bool even = false;
    bool expectedEven = false;

    //getExpected Link ID
    rctDataDecoder.getExpectedLinkID(iLink, expectedCrateID, expectedEven);
    //getDecodedLink ID
    rctDataDecoder.decodeLinkID(linkID, crateID, even);

    //Perform a check to see if the link ID is as expected, if not then report an error but continue unpacking
    if (expectedCrateID != crateID || even != expectedEven) {
      LogError("L1T") << "Expected Crate ID " << expectedCrateID << " expectedEven " << expectedEven
                      << "does not match actual Crate ID " << crateID << " even " << even;
    }

    if (CRCErrCnt != 0)
      LogError("L1T") << "WARNING CRC ErrorFound linkID " << linkID << " expected crateID " << expectedCrateID;

    // Loop over multiple BX
    for (uint32_t iBX = 0; iBX < nBX; iBX++) {
      allLinks[iBX][iLink].uint.reserve(6);
      allLinks[iBX][iLink].ctp7LinkNumber = iLink;
      allLinks[iBX][iLink].crateID = expectedCrateID;
      allLinks[iBX][iLink].even = expectedEven;

      //Notice 6 words per BX
      for (unsigned int iWord = 0; iWord < 6; iWord++) {
        allLinks[iBX][iLink].uint.push_back(data[iDAQBuffer + iWord + iBX * 6]);
      }
    }
  }

  //Step 2: Dynamically match links and create RCTInfo Objects
  uint32_t nCratesFound = 0;
  for (unsigned int iCrate = 0; iCrate < 18; iCrate++) {
    bool foundEven = false, foundOdd = false;
    link_data even[5];
    link_data odd[5];

    for (unsigned int iLink = 0; iLink < NLinks; iLink++) {
      if ((allLinks[0][iLink].crateID == iCrate) && (allLinks[0][iLink].even == true)) {
        foundEven = true;
        for (unsigned int iBX = 0; iBX < nBX; iBX++)
          even[iBX] = allLinks[iBX][iLink];
      } else if ((allLinks[0][iLink].crateID == iCrate) && (allLinks[0][iLink].even == false)) {
        foundOdd = true;
        for (unsigned int iBX = 0; iBX < nBX; iBX++)
          odd[iBX] = allLinks[iBX][iLink];
      }

      //if success then create RCTInfoVector and fill output object
      if (foundEven && foundOdd) {
        nCratesFound++;

        //fill rctInfoVector for all BX read out
        for (unsigned int iBX = 0; iBX < nBX; iBX++) {
          //RCTInfoFactory rctInfoFactory;
          std::vector<RCTInfo> rctInfoData;
          rctDataDecoder.decodeLinks(even[iBX].uint, odd[iBX].uint, rctInfoData);
          rctDataDecoder.setRCTInfoCrateID(rctInfoData, iCrate);
          allCrateRCTInfo[iBX].push_back(rctInfoData.at(0));
        }
        break;
      }
    }
  }

  if (nCratesFound != 18)
    LogError("L1T") << "Warning -- only found " << nCratesFound << " valid crates";

  //start assuming 1 BX readout
  int32_t startBX = 0;
  if (nBX == 1)
    startBX = 0;
  else if (nBX == 3)
    startBX = -1;
  else if (nBX == 5)
    startBX = -2;

  //Step 3: Create Collections from RCTInfo Objects
  //Notice, iBX used for grabbing data from array, startBX used for storing in Collection
  for (uint32_t iBX = 0; iBX < nBX; iBX++, startBX++) {
    for (unsigned int iCrate = 0; iCrate < nCratesFound; iCrate++) {
      RCTInfo rctInfo = allCrateRCTInfo[iBX].at(iCrate);
      //Use Crate ID to identify eta/phi of candidate
      for (int j = 0; j < 4; j++) {
        L1CaloEmCand em = L1CaloEmCand(rctInfo.neRank[j], rctInfo.neRegn[j], rctInfo.neCard[j], rctInfo.crateID, false);
        em.setBx(startBX);
        colls->rctEm()->push_back(em);
      }

      for (int j = 0; j < 4; j++) {
        L1CaloEmCand em = L1CaloEmCand(rctInfo.ieRank[j], rctInfo.ieRegn[j], rctInfo.ieCard[j], rctInfo.crateID, true);
        em.setBx(startBX);
        colls->rctEm()->push_back(em);
      }

      for (int j = 0; j < 7; j++) {
        for (int k = 0; k < 2; k++) {
          bool o = (((rctInfo.oBits >> (j * 2 + k)) & 0x1) == 0x1);
          bool t = (((rctInfo.tBits >> (j * 2 + k)) & 0x1) == 0x1);
          bool m = (((rctInfo.mBits >> (j * 2 + k)) & 0x1) == 0x1);
          bool q = (((rctInfo.qBits >> (j * 2 + k)) & 0x1) == 0x1);
          L1CaloRegion rgn = L1CaloRegion(rctInfo.rgnEt[j][k], o, t, m, q, rctInfo.crateID, j, k);
          rgn.setBx(startBX);
          colls->rctCalo()->push_back(rgn);
        }
      }

      for (int k = 0; k < 4; k++) {
        for (int j = 0; j < 2; j++) {
          // 0 1 4 5 2 3 6 7
          uint32_t offset = j * 2 + k % 2 + (k / 2) * 4;
          bool fg = (((rctInfo.hfQBits >> offset) & 0x1) == 0x1);
          L1CaloRegion rgn = L1CaloRegion(rctInfo.hfEt[j][k], fg, rctInfo.crateID, (j * 4 + k));
          rgn.setBx(startBX);
          colls->rctCalo()->push_back(rgn);
        }
      }
    }
  }
}

Member Data Documentation

amc13HeaderSize_

const unsigned RctRawToDigi::amc13HeaderSize_ = 16

staticprivate

Definition at line 75 of file RctRawToDigi.h.

Referenced by produce().

amc13TrailerSize_

const unsigned RctRawToDigi::amc13TrailerSize_ = 8

staticprivate

Definition at line 78 of file RctRawToDigi.h.

Referenced by produce().

ctp7HeaderSize_

const unsigned RctRawToDigi::ctp7HeaderSize_ = 8

staticprivate

Definition at line 81 of file RctRawToDigi.h.

ctp7TrailerSize_

const unsigned RctRawToDigi::ctp7TrailerSize_ = 8

staticprivate

Definition at line 84 of file RctRawToDigi.h.

errorCounters_

std::vector<unsigned> RctRawToDigi::errorCounters_

private

Counts number of errors for each code (index)

Definition at line 104 of file RctRawToDigi.h.

Referenced by endJob().

errors_

L1TriggerErrorCollection* RctRawToDigi::errors_

private

pointer to error collection

Definition at line 103 of file RctRawToDigi.h.

fedId_

int RctRawToDigi::fedId_

private

RCT FED ID.

Definition at line 94 of file RctRawToDigi.h.

Referenced by produce(), and RctRawToDigi().

inputLabel_

edm::InputTag RctRawToDigi::inputLabel_

private

FED collection label.

Definition at line 93 of file RctRawToDigi.h.

Referenced by produce(), and RctRawToDigi().

MAX_DATA

const unsigned RctRawToDigi::MAX_DATA = 4680

staticprivate

The maximum number of blocks we will try to unpack before thinking something is wrong.

Definition at line 87 of file RctRawToDigi.h.

MAX_ERR_CODE

const unsigned RctRawToDigi::MAX_ERR_CODE = 6

staticprivate

Definition at line 102 of file RctRawToDigi.h.

Referenced by endJob().

MIN_DATA

const unsigned RctRawToDigi::MIN_DATA = 900

staticprivate

The minimum number of blocks we will try to unpack before thinking something is wrong (really this should be 920, to be tested)

Definition at line 90 of file RctRawToDigi.h.

Referenced by produce().

sLinkHeaderSize_

const unsigned RctRawToDigi::sLinkHeaderSize_ = 8

staticprivate

Definition at line 69 of file RctRawToDigi.h.

Referenced by produce().

sLinkTrailerSize_

const unsigned RctRawToDigi::sLinkTrailerSize_ = 8

staticprivate

Definition at line 72 of file RctRawToDigi.h.

Referenced by produce(), and unpack().

unpackFailures_

unsigned RctRawToDigi::unpackFailures_

private

To count the total number of RCT unpack failures.

Definition at line 105 of file RctRawToDigi.h.

verbose_

const bool RctRawToDigi::verbose_

private

If true, then debug print out for each event.

Definition at line 96 of file RctRawToDigi.h.

Referenced by endJob().