CaloLayer1Packer.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Event.h"
#include "EventFilter/L1TRawToDigi/plugins/PackerFactory.h"

#include "L1Trigger/L1TCalorimeter/interface/CaloTools.h"

#include "CaloLayer1Packer.h"

namespace l1t {
  namespace stage2 {

    // max_iEta_HcalTP = 41; // barrel <= 16, endcap <= 28, hf <= 41
    // there are two TT29’s: one in HE readout in TT28 and another in HF readout in TT30
    // max_iPhi_HcalTP = 72;

    Blocks CaloLayer1Packer::pack(const edm::Event& event, const PackerTokens* toks) {
      edm::Handle<EcalTrigPrimDigiCollection> ecalDigis;
      event.getByToken(static_cast<const CaloLayer1Tokens*>(toks)->getEcalDigiToken(), ecalDigis);
      edm::Handle<HcalTrigPrimDigiCollection> hcalDigis;
      event.getByToken(static_cast<const CaloLayer1Tokens*>(toks)->getHcalDigiToken(), hcalDigis);
      edm::Handle<L1CaloRegionCollection> caloRegions;
      event.getByToken(static_cast<const CaloLayer1Tokens*>(toks)->getCaloRegionToken(), caloRegions);

      std::vector<uint32_t> load;
      load.resize(192, 0u);

      auto ctp7_phi = board();
      uint32_t* ptr = load.data();
      UCTCTP7RawData ctp7Data(ptr);
      makeECalTPGs(ctp7_phi, ctp7Data, ecalDigis.product());
      makeHCalTPGs(ctp7_phi, ctp7Data, hcalDigis.product());
      makeHFTPGs(ctp7_phi, ctp7Data, hcalDigis.product());
      makeRegions(ctp7_phi, ctp7Data, caloRegions.product());

      unsigned bx_per_l1a = 1;
      // CTP7 uses CMS scheme, starting at 0
      // TODO: expected +2, but +1 apparently?
      unsigned calo_bxid = (event.bunchCrossing() + 1) % 3564;

      // a la CTP7Payload::getHeader()
      unsigned blockId = 0;
      unsigned blockSize = 192;
      unsigned capId = 0;
      unsigned blockFlags = ((bx_per_l1a & 0xf) << 16) | (calo_bxid & 0xfff);
      BlockHeader hdr(blockId, blockSize, capId, blockFlags, CTP7);
      Block block(hdr, &*load.begin(), &*load.end());

      Blocks res;
      res.push_back(block);
      return res;
    }

    void CaloLayer1Packer::makeECalTPGs(uint32_t lPhi,
                                        UCTCTP7RawData& ctp7Data,
                                        const EcalTrigPrimDigiCollection* ecalTPGs) {
      UCTCTP7RawData::CaloType cType = UCTCTP7RawData::EBEE;
      for (uint32_t iPhi = 0; iPhi < 4; iPhi++) {  // Loop over all four phi divisions on card
        int cPhi = -1 + lPhi * 4 + iPhi;           // Calorimeter phi index
        if (cPhi == 0)
          cPhi = 72;
        else if (cPhi == -1)
          cPhi = 71;
        else if (cPhi < -1) {
          edm::LogError("CaloLayer1Packer") << "Major error in makeECalTPGs" << std::endl;
          return;
        }
        for (int cEta = -28; cEta <= 28; cEta++) {  // Calorimeter Eta indices (HB/HE for now)
          if (cEta != 0) {                          // Calorimeter eta = 0 is invalid
            bool negativeEta = false;
            if (cEta < 0)
              negativeEta = true;
            uint32_t iEta = abs(cEta);

            int zSide = cEta / ((int)iEta);
            const EcalSubdetector ecalTriggerTower =
                (iEta > 17) ? EcalSubdetector::EcalEndcap : EcalSubdetector::EcalBarrel;
            EcalTrigTowerDetId id(zSide, ecalTriggerTower, iEta, cPhi);
            const auto& tp = ecalTPGs->find(id);
            if (tp != ecalTPGs->end()) {
              ctp7Data.setET(cType, negativeEta, iEta, iPhi, tp->compressedEt());
              ctp7Data.setFB(cType, negativeEta, iEta, iPhi, tp->fineGrain());
            }
          }
        }
      }
    }

    void CaloLayer1Packer::makeHCalTPGs(uint32_t lPhi,
                                        UCTCTP7RawData& ctp7Data,
                                        const HcalTrigPrimDigiCollection* hcalTPGs) {
      UCTCTP7RawData::CaloType cType = UCTCTP7RawData::HBHE;
      for (uint32_t iPhi = 0; iPhi < 4; iPhi++) {  // Loop over all four phi divisions on card
        int cPhi = -1 + lPhi * 4 + iPhi;           // Calorimeter phi index
        if (cPhi == 0)
          cPhi = 72;
        else if (cPhi == -1)
          cPhi = 71;
        else if (cPhi < -1) {
          edm::LogError("CaloLayer1Packer") << "Major error in makeHCalTPGs" << std::endl;
          return;
        }
        for (int cEta = -28; cEta <= 28; cEta++) {  // Calorimeter Eta indices (HB/HE for now)
          if (cEta != 0) {                          // Calorimeter eta = 0 is invalid
            bool negativeEta = false;
            if (cEta < 0)
              negativeEta = true;
            uint32_t iEta = abs(cEta);

            HcalTrigTowerDetId id(cEta, cPhi);
            const auto tp = hcalTPGs->find(id);

            if (tp != hcalTPGs->end()) {
              uint32_t fg_bits = 0;
              for (int index = 0; index < 6; index++)
                fg_bits |= tp->SOI_fineGrain(index) << index;

              ctp7Data.setET(cType, negativeEta, iEta, iPhi, tp->SOI_compressedEt());
              ctp7Data.setFB(cType, negativeEta, iEta, iPhi, fg_bits);
            }
          }
        }
      }
    }

    void CaloLayer1Packer::makeHFTPGs(uint32_t lPhi,
                                      UCTCTP7RawData& ctp7Data,
                                      const HcalTrigPrimDigiCollection* hcalTPGs) {
      UCTCTP7RawData::CaloType cType = UCTCTP7RawData::HF;
      for (uint32_t side = 0; side <= 1; side++) {
        bool negativeEta = false;
        if (side == 0)
          negativeEta = true;
        for (uint32_t iEta = 30; iEta <= 40; iEta++) {
          for (uint32_t iPhi = 0; iPhi < 2; iPhi++) {
            if (iPhi == 1 && iEta == 40)
              iEta = 41;
            int cPhi = 1 + lPhi * 4 + iPhi * 2;  // Calorimeter phi index: 1, 3, 5, ... 71
            if (iEta == 41)
              cPhi -= 2;                  // Last two HF are 3, 7, 11, ...
            cPhi = (cPhi + 69) % 72 + 1;  // cPhi -= 2 mod 72
            int cEta = iEta;
            if (negativeEta)
              cEta = -iEta;

            HcalTrigTowerDetId id(cEta, cPhi);
            id.setVersion(1);  // To not process these 1x1 HF TPGs with RCT
            const auto tp = hcalTPGs->find(id);
            if (tp != hcalTPGs->end()) {
              ctp7Data.setET(cType, negativeEta, iEta, iPhi, tp->SOI_compressedEt());
              ctp7Data.setFB(cType, negativeEta, iEta, iPhi, ((tp->SOI_fineGrain(1) << 1) | tp->SOI_fineGrain(0)));
            }
          }
        }
      }
    }

    void CaloLayer1Packer::makeRegions(uint32_t lPhi, UCTCTP7RawData& ctp7Data, const L1CaloRegionCollection* regions) {
      for (uint32_t side = 0; side <= 1; side++) {
        bool negativeEta = false;
        if (side == 0)
          negativeEta = true;
        for (uint32_t region = 0; region <= 6; region++) {
          uint32_t lEta = 10 - region;  // GCT eta goes 0-21, 0-3 -HF, 4-10 -B/E, 11-17 +B/E, 18-21 +HF
          if (!negativeEta)
            lEta = region + 11;

          L1CaloRegionDetId id(lEta, lPhi);
          // Can't use find since not an edm::SortedCollection
          // const L1CaloRegion& rtp = *regions->find(id);
          for (const auto& rtp : *regions) {
            if (rtp.id() == id) {
              ctp7Data.setRegionSummary(negativeEta, region, rtp.raw());
              break;
            }
          }
        }
      }
    }

  }  // namespace stage2
}  // namespace l1t

DEFINE_L1T_PACKER(l1t::stage2::CaloLayer1Packer);