EGammaPacker.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 EGammaPacker : public Packer {
         public:
            virtual Blocks pack(const edm::Event&, const PackerTokens*) override;
      };
   }
}

// Implementation

namespace l1t {
namespace stage2 {
   Blocks
   EGammaPacker::pack(const edm::Event& event, const PackerTokens* toks)
   {
      edm::Handle<EGammaBxCollection> egs;
      event.getByToken(static_cast<const CaloTokens*>(toks)->getEGammaToken(), egs);

      std::vector<uint32_t> load1, load2;

      for (int i = egs->getFirstBX(); i <= egs->getLastBX(); ++i) {
    
    for (auto j = egs->begin(i); j != egs->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 |
        (j->hwQual() & 0x7) << 27;
    
      if (load1.size() < l1t::stage2::layer2::demux::nEGPerLink) load1.push_back(word);
      else load2.push_back(word);
    
    }
      }
      
      // push zeroes if jets are missing                                      
      while (load1.size()<l1t::stage2::layer2::demux::nOutputFramePerBX) load1.push_back(0);
      
      while (load2.size()<l1t::stage2::layer2::demux::nOutputFramePerBX) load2.push_back(0);
  
      return {Block(9, load1), Block(11, load2)};

   }

}
}

DEFINE_L1T_PACKER(l1t::stage2::EGammaPacker);