TauPacker.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Event.h"

#include "EventFilter/L1TRawToDigi/interface/Packer.h"

#include "CaloTokens.h"

#include "L1TStage2Layer2Constants.h"

namespace l1t {
   namespace stage2 {
      class TauPacker : public Packer {
         public:
        TauPacker(int b1, int b2) : b1_(b1), b2_(b2) {}
            virtual Blocks pack(const edm::Event&, const PackerTokens*) override;
        int b1_, b2_;
      };

      class GTTauPacker : public TauPacker {
         public:
             GTTauPacker() : TauPacker(16,18) {}
      };
      class CaloTauPacker : public TauPacker {
         public:
         CaloTauPacker() : TauPacker(17,19) {}
      };

   }
}

// Implementation

namespace l1t {
namespace stage2 {
   Blocks
   TauPacker::pack(const edm::Event& event, const PackerTokens* toks)
   {
      edm::Handle<TauBxCollection> taus;
      event.getByToken(static_cast<const CommonTokens*>(toks)->getTauToken(), taus);

      std::vector<uint32_t> load1, load2;

      int nBx=0;
      for (int i = taus->getFirstBX(); i <= taus->getLastBX(); ++i) {

    for (auto j = taus->begin(i); j != taus->end(i); ++j) {

        uint32_t packed_eta = abs(j->hwEta()) & 0x7F;
        if (j->hwEta() < 0){
           packed_eta = (128 - packed_eta) | 1<<7;
        }

            uint32_t word = \
	      std::min(j->hwPt(), 0x1FF) |
          packed_eta << 9 |
          (j->hwPhi() & 0xFF) << 17 |
          (j->hwIso() & 0x3) << 25 |//two bits for tau isolation
          (j->hwQual() & 0x7) << 27;

        // Depending on the number, push onto appropriate link
        if (load1.size() - nBx*l1t::stage2::layer2::demux::nTauPerLink < l1t::stage2::layer2::demux::nTauPerLink) load1.push_back(word);
        else load2.push_back(word);

        }
    
    // push zeroes if not enough taus must do this for each BX and each link                                     
    while (load1.size() - nBx*l1t::stage2::layer2::demux::nOutputFramePerBX <l1t::stage2::layer2::demux::nOutputFramePerBX) load1.push_back(0);
    while (load2.size() - nBx*l1t::stage2::layer2::demux::nOutputFramePerBX <l1t::stage2::layer2::demux::nOutputFramePerBX) load2.push_back(0);
    
    nBx++;
      }
      
     

      return {Block(b1_, load1), Block(b2_, load2)};

   }
}
}

DEFINE_L1T_PACKER(l1t::stage2::GTTauPacker);
DEFINE_L1T_PACKER(l1t::stage2::CaloTauPacker);