GlobalExtBlkUnpacker.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "EventFilter/L1TRawToDigi/plugins/UnpackerFactory.h"

#include "GTCollections.h"
#include "GlobalExtBlkUnpacker.h"

namespace l1t {
namespace stage2 {
   bool
   GlobalExtBlkUnpacker::unpack(const Block& block, UnpackerCollections *coll)
   {

     LogDebug("L1T") << "Block ID  = " << block.header().getID() << " size = " << block.header().getSize();

     unsigned int wdPerBX = 6; //should this be configured someplace else?
     int nBX = int(ceil(block.header().getSize() / 6.)); // FOR GT Not sure what we have here...put at 6 because of 6 frames. Since there are 12 EGamma objects reported per event (see CMS IN-2013/005)

     // Find the central, first and last BXs
     int firstBX = -(ceil((double)nBX/2.)-1);
     int lastBX;
     if (nBX % 2 == 0) {
       lastBX = ceil((double)nBX/2.);
     } else {
       lastBX = ceil((double)nBX/2.)-1;
     }

     auto res_ = static_cast<GTCollections*>(coll)->getExts();
     res_->setBXRange(firstBX, lastBX);

     LogDebug("L1T") << "nBX = " << nBX << " first BX = " << firstBX << " lastBX = " << lastBX;

     // Loop over multiple BX and then number of EG cands filling collection
     int numBX = 0; //positive int to count BX
     for (int bx=firstBX; bx<=lastBX; bx++){

        // If this is the first block, instantiate GlobalExt so it is there to fill from mult. blocks
       if(block.header().getID()==24) {

      LogDebug("L1T") << "Creating GT External Block for BX =" << bx;
          GlobalExtBlk tExt = GlobalExtBlk();
          res_->push_back(bx,tExt);

       }

       //fetch from collection
       GlobalExtBlk ext = res_->at(bx,0);

       //Determine offset of algorithm bits based on block.ID
       // ID=24  offset = 0;  ID=26  offset=64;  ID=28  offset=128=2*64; ID=30 offset=3*64=192 
       int extOffset = ((block.header().getID()-24)/2)*64;

       for(unsigned int wd=0;  wd<wdPerBX; wd++) {
         uint32_t raw_data = block.payload()[wd+numBX*wdPerBX];
     LogDebug("L1T") << " payload word " << wd << " 0x" << hex << raw_data;

         if(wd < 2) {
           for(unsigned int bt=0; bt<32; bt++) {
         int val = ((raw_data >> bt) & 0x1);
         int extBit = bt+wd*32+extOffset;
         if(val==1) ext.setExternalDecision(extBit,true);
           }
     }  
     
       }

       // Put the object back into place (Must be better way???)
       res_->set(bx,0,ext);

       //ext.print(std::cout);     
       numBX++;
     }

     return true;
   }
}
}

DEFINE_L1T_UNPACKER(l1t::stage2::GlobalExtBlkUnpacker);