l1t::stage2::CaloTowerPacker Class Reference

Inheritance diagram for l1t::stage2::CaloTowerPacker:

Public Member Functions
virtual Blocks	pack (const edm::Event &, const PackerTokens *) override

Detailed Description

Definition at line 10 of file CaloTowerPacker.cc.

Member Function Documentation

Blocks l1t::stage2::CaloTowerPacker::pack ( const edm::Event &  event, const PackerTokens *  toks  )

overridevirtual

Implements l1t::Packer.

Definition at line 22 of file CaloTowerPacker.cc.

References l1t::CaloTools::caloTowerHash(), eta, l1t::CaloTower::hwEtRatio(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), i, svgfig::load(), min(), phi(), and HLT_25ns14e33_v1_cff::towers.

Referenced by pyrootRender.interactiveRender::draw().

   {

      edm::Handle<CaloTowerBxCollection> towers;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getCaloTowerToken(), towers);

      Blocks res;

      for (int i = towers->getFirstBX(); i <= towers->getLastBX(); ++i) {

        for (int phi = 1; phi <=72; phi=phi+2) { // Two phi values per link

          unsigned int id = 2*phi - 2; // Block IDs start at zero and span even numbers up to 142
          std::vector<uint32_t> load;

          for (int eta = 1; eta <=41; eta++) { // This is abs(eta) since +/- eta are interleaved in time

            // Get four towers +/- eta and phi and phi+1 to all be packed in this loop
            l1t::CaloTower t1 = towers->at(i,l1t::CaloTools::caloTowerHash(eta,phi));
            l1t::CaloTower t2 = towers->at(i,l1t::CaloTools::caloTowerHash(eta,phi+1));
            l1t::CaloTower t3 = towers->at(i,l1t::CaloTools::caloTowerHash(-1*eta,phi));
            l1t::CaloTower t4 = towers->at(i,l1t::CaloTools::caloTowerHash(-1*eta,phi+1));

            // Merge phi and phi+1 into one block (phi is LSW, phi+1 is MSW)
            uint32_t word1 = \
	      std::min(t1.hwPt(), 0x1FF) |
          (t1.hwEtRatio() & 0x7) << 9 |
          (t1.hwQual() & 0xF) << 12;

        word1 = word1 |
          std::min(t2.hwPt(), 0x1FF) << 16 |
          (t2.hwEtRatio() & 0x7) << 25 |
          (t2.hwQual() & 0xF) << 28;

            load.push_back(word1);

            // Do it all again for -eta

            uint32_t word2 = \
              std::min(t3.hwPt(), 0x1FF) |
              (t3.hwEtRatio() & 0x7) << 9 |
              (t3.hwQual() & 0xF) << 12;

            word2 = word2 |
              std::min(t4.hwPt(), 0x1FF) << 16 |
              (t4.hwEtRatio() & 0x7) << 25 |
              (t4.hwQual() & 0xF) << 28;

            load.push_back(word2);

          }

          res.push_back(Block(id, load));

        }
      }
      return res;
   }