BMTFUnpackerInputs.cc

Go to the documentation of this file.

#include "EventFilter/L1TRawToDigi/plugins/UnpackerFactory.h"

#include "BMTFUnpackerInputs.h"

#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambPhContainer.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambThContainer.h"

namespace l1t {
  namespace stage2 {
    void numWheelSectorTrTag_bmtf(int& wheelNo, int& tagSegID, int linkNo, int amcNo) {
      if (linkNo >= 0 && linkNo < 6)
        wheelNo = -2;
      else if (linkNo >= 8 && linkNo < 14)
        wheelNo = -1;
      else if (linkNo >= 16 && linkNo < 22)
        wheelNo = 0;
      else if (linkNo >= 22 && linkNo < 28)
        wheelNo = 1;
      else if ((linkNo >= 28 && linkNo < 30) || (linkNo >= 32 && linkNo < 36))
        wheelNo = 2;

      if (linkNo % 2 == 0)
        tagSegID = 0;
      else
        tagSegID = 1;
    }

    bool checkQual_bmtf(const unsigned int& value, const bool& isNewFw) {
      if (isNewFw)
        return (value == 7);
      else
        return (value == 0);
    }

    bool unpacking_bmtf(const Block& block, UnpackerCollections* coll, qualityHits& linkAndQual_, const bool& isNewFw) {
      unsigned int ownLinks[] = {4, 5, 12, 13, 20, 21, 22, 23, 28, 29};
      bool ownFlag(false);

      //Checks if the given block coresponds to 1 of the OWN links
      for (int i = 0; i < 10; i++) {
        if (block.header().getID() / 2 == ownLinks[i])
          ownFlag = true;
      }
      if (!ownFlag)  //if not returns that the "job here done"
        return true;

      //Get header ID and payload from the given Block
      unsigned int blockId = block.header().getID();
      LogDebug("L1T") << "Block ID: " << blockId << " size: " << block.header().getSize();

      //Make output CMSSW collections
      L1MuDTChambPhContainer* resPhi = static_cast<BMTFCollections*>(coll)->getInMuonsPh();
      L1MuDTChambThContainer* resThe = static_cast<BMTFCollections*>(coll)->getInMuonsTh();

      //Get input phi & eta Containers
      L1MuDTChambPhContainer::Phi_Container phiData = *(resPhi->getContainer());
      L1MuDTChambThContainer::The_Container theData = *(resThe->getContainer());

      //ZeroSuppresion Handler
      BxBlocks bxBlocks;
      bool ZS_enabled =
          (bool)((block.header().getFlags() >> 1) & 0x01);  //getFlags() returns first 8-bits from the amc header
      if (ZS_enabled)
        bxBlocks =
            block.getBxBlocks((unsigned int)6, true);  //it returnes 7-32bit bxBlocks originated from the amc13 Block
      else
        bxBlocks =
            block.getBxBlocks((unsigned int)6, false);  //it returnes 6-32bit bxBlocks originated from the amc13 Block

      for (const auto& ibx : bxBlocks)  //Bx iteration
      {
        int bxNum = ibx.header().getBx();
        uint32_t inputWords[ibx.getSize()];  //array of 6 uint32_t payload-words (size of the payload in the BxBlock)

        //Note
        /*In the non-ZS fashion, the expression "block.header().getSize()/nBX" was 6 in any case
        the reason is that the size is 6 or 30, and these numbers are divided by 1 or 5 respectively.*/

        //Fill the above uint32_t array
        for (unsigned int iw = 0; iw < ibx.getSize(); iw++)
          inputWords[iw] = (ibx.payload())[iw];

        int wheel, sector, trTag;                                                         //Container information
        numWheelSectorTrTag_bmtf(wheel, trTag, blockId / 2, block.amc().getAMCNumber());  //this returns wheel & tsTag
        sector = block.amc().getBoardID() - 1;

        //Check if the sector is "out of range" - (trys then to use AMC13 information?)
        if (sector < 0 || sector > 11) {
          edm::LogInfo("l1t:stage2::BMTFUnpackerInputs::unpack")
              << "Sector found out of range so it will be calculated by the slot number";
          if (block.amc().getAMCNumber() % 2 != 0)
            sector = block.amc().getAMCNumber() / 2;
          else
            sector = 6 + (block.amc().getAMCNumber() / 2 - 1);
        }

        int mbPhi[4], mbPhiB[4], mbQual[4], mbBxC[4], mbRPC[4];  //Container information
        //mbPhiB[2] = 0;

        for (int iw = 0; iw < 4; iw++)  // 4 phi (32-bit) words
        {
          if (((inputWords[iw] & 0xfffffff) == 0) || (inputWords[iw] == 0x505050bc))
            continue;
          else if ((inputWords[iw] != 0x505050bc) && (inputWords[iw + 2] == 0x505050bc))
            continue;

          if (((inputWords[iw] >> 11) & 0x1) == 1)
            mbPhi[iw] = (inputWords[iw] & 0x7FF) - 2048;
          else
            mbPhi[iw] = (inputWords[iw] & 0xFFF);

          if (((inputWords[iw] >> 21) & 0x1) == 1)
            mbPhiB[iw] = ((inputWords[iw] >> 12) & 0x1FF) - 512;
          else
            mbPhiB[iw] = (inputWords[iw] >> 12) & 0x3FF;

          mbQual[iw] = (inputWords[iw] >> 22) & 0x7;
          mbRPC[iw] = (inputWords[iw] >> 26) & 0x1;
          mbBxC[iw] = (inputWords[iw] >> 30) & 0x3;

          //if (mbQual[iw] == 0)
          if (checkQual_bmtf(mbQual[iw], isNewFw))
            continue;

          phiData.push_back(L1MuDTChambPhDigi(
              bxNum, wheel, sector, iw + 1, mbPhi[iw], mbPhiB[iw], mbQual[iw], trTag, mbBxC[iw], mbRPC[iw]));

        }  //4 phi words

        int etaHits[3][7];  //Container information
        bool zeroFlag[3];
        for (int i = 0; i < 3; i++)  // 3 eta (7-bit) words
        {
          zeroFlag[i] = false;
          for (int j = 0; j < 7; j++) {
            etaHits[i][j] = (inputWords[4] >> (i * 7 + j)) & 0x1;
            if (etaHits[i][j] != 0)
              zeroFlag[i] = true;
          }
        }  //3 eta words

        if (trTag == 1) {
          for (int i = 0; i < 3; i++) {
            if (zeroFlag[i])
              theData.push_back(L1MuDTChambThDigi(bxNum, wheel, sector, i + 1, etaHits[i], linkAndQual_.hits[i]));
          }

        } else {
          /*
        qualityHits temp;
        temp.linkNo = blockId/2;
        std::copy(&etaHits[0][0], &etaHits[0][0]+3*7,&temp.hits[0][0]);
        linkAndQual_[blockId/2] = temp; 
          */
          linkAndQual_.linkNo = blockId / 2;
          std::copy(&etaHits[0][0], &etaHits[0][0] + 3 * 7, &linkAndQual_.hits[0][0]);
        }

      }  //iBxBlock

      //Fill Containers
      resThe->setContainer(theData);
      resPhi->setContainer(phiData);

      return true;
    }

    bool BMTFUnpackerInputsOldQual::unpack(const Block& block, UnpackerCollections* coll) {
      return unpacking_bmtf(block, coll, linkAndQual_, false);
    }  //unpack old quality

    bool BMTFUnpackerInputsNewQual::unpack(const Block& block, UnpackerCollections* coll) {
      return unpacking_bmtf(block, coll, linkAndQual_, true);
    }  //unpack new quality
  }    // namespace stage2
}  // namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage2::BMTFUnpackerInputsOldQual);
DEFINE_L1T_UNPACKER(l1t::stage2::BMTFUnpackerInputsNewQual);