EtSumUnpacker_0x10010057.cc

Go to the documentation of this file.

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

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

#include "L1TObjectCollections.h"

#include "L1TStage2Layer2Constants.h"
#include "EtSumUnpacker_0x10010057.h"

namespace l1t {
namespace stage2 {
  EtSumUnpacker_0x10010057::EtSumUnpacker_0x10010057() : EtSumCopy_(0)
  {
  }
  
   bool
   EtSumUnpacker_0x10010057::unpack(const Block& block, UnpackerCollections *coll)
   {

     using namespace l1t::stage2::layer2;

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

     int nBX = int(ceil(block.header().getSize() / demux::nOutputFramePerBX)); // Since there 6 frames per demux output event
     // expect the first four frames to be the first 4 EtSum 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<L1TObjectCollections*>(coll)->getEtSums(EtSumCopy_);
     res_->setBXRange(firstBX, lastBX);

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

     // Loop over multiple BX and fill EtSums collection
     for (int bx=firstBX; bx<=lastBX; bx++){


       // ET
       int iFrame = (bx-firstBX)*demux::nOutputFramePerBX;

       uint32_t raw_data = block.payload().at(iFrame);

       l1t::EtSum et = l1t::EtSum();
    
       et.setHwPt(raw_data & 0xFFF);
       et.setType(l1t::EtSum::kTotalEt);       
       et.setP4( l1t::CaloTools::p4Demux(&et) );

       LogDebug("L1T") << "ET: pT " << et.hwPt() << " bx " << bx;

       res_->push_back(bx,et);


       // ET EM

       l1t::EtSum etem = l1t::EtSum();
    
       etem.setHwPt( (raw_data >> 12) & 0xFFF);
       etem.setType(l1t::EtSum::kTotalEtEm);       
       etem.setP4( l1t::CaloTools::p4Demux(&etem) );

       LogDebug("L1T") << "ETEM: pT " << etem.hwPt() << " bx " << bx;

       res_->push_back(bx,etem);


       // MBHFPT0

       l1t::EtSum mbp0 = l1t::EtSum();
       mbp0.setHwPt( (raw_data>>28) & 0xf );
       mbp0.setType( l1t::EtSum::kMinBiasHFP0 );

       res_->push_back(bx, mbp0);


       // HT

       raw_data = block.payload()[iFrame+1];

       l1t::EtSum ht = l1t::EtSum();

       ht.setHwPt(raw_data & 0xFFF);
       ht.setType(l1t::EtSum::kTotalHt);       
       ht.setP4( l1t::CaloTools::p4Demux(&ht) );

       LogDebug("L1T") << "HT: pT " << ht.hwPt();

       res_->push_back(bx,ht);

       //MBHFMT0

       l1t::EtSum mbm0 = l1t::EtSum();
       mbm0.setHwPt( (raw_data>>28) & 0xf );
       mbm0.setType( l1t::EtSum::kMinBiasHFM0 );

       res_->push_back(bx, mbm0);


       //  MET (no HF)

       raw_data = block.payload()[iFrame+2];

       l1t::EtSum met = l1t::EtSum();
    
       met.setHwPt(raw_data & 0xFFF);
       met.setHwPhi((raw_data >> 12) & 0xFF);
       met.setType(l1t::EtSum::kMissingEt);       
       met.setP4( l1t::CaloTools::p4Demux(&met) );

       LogDebug("L1T") << "MET: phi " << met.hwPhi() << " pT " << met.hwPt() << " bx " << bx;

       res_->push_back(bx,met);

       // MBHFPT1

       l1t::EtSum mbp1 = l1t::EtSum();
       mbp1.setHwPt( (raw_data>>28) & 0xf );
       mbp1.setType( l1t::EtSum::kMinBiasHFP1 );

       res_->push_back(bx, mbp1);

       // ET ASYMM

       l1t::EtSum etAsym = l1t::EtSum();
       etAsym.setHwPt( (raw_data>>20) & 0xFF );
       etAsym.setType( l1t::EtSum::kAsymEt );

       res_->push_back(bx, etAsym);

       // MHT 

       raw_data = block.payload()[iFrame+3];

       l1t::EtSum mht = l1t::EtSum();
    
       mht.setHwPt(raw_data & 0xFFF);
       mht.setHwPhi((raw_data >> 12) & 0xFF);
       mht.setType(l1t::EtSum::kMissingHt);       
       mht.setP4( l1t::CaloTools::p4Demux(&mht) );

       LogDebug("L1T") << "MHT: phi " << mht.hwPhi() << " pT " << mht.hwPt() << " bx " << bx;

       res_->push_back(bx,mht);

       // MBHFMT1

       l1t::EtSum mbm1 = l1t::EtSum();
       mbm1.setHwPt( (raw_data>>28) & 0xf );
       mbm1.setType( l1t::EtSum::kMinBiasHFM1 );

       res_->push_back(bx, mbm1);

       // HT ASYMM

       l1t::EtSum htAsym = l1t::EtSum();
       htAsym.setHwPt( (raw_data>>20) & 0xFF );
       htAsym.setType( l1t::EtSum::kAsymHt );

       res_->push_back(bx, htAsym);

       
       //  MET (with HF)

       raw_data = block.payload()[iFrame+4];

       l1t::EtSum methf = l1t::EtSum();
    
       methf.setHwPt(raw_data & 0xFFF);
       methf.setHwPhi((raw_data >> 12) & 0xFF);
       methf.setType(l1t::EtSum::kMissingEtHF);       
       methf.setP4( l1t::CaloTools::p4Demux(&methf) );

       LogDebug("L1T") << "METHF: phi " << methf.hwPhi() << " pT " << methf.hwPt() << " bx " << bx;

       res_->push_back(bx,methf);

       // ETHF ASYMM

       l1t::EtSum etHFAsym = l1t::EtSum();
       etHFAsym.setHwPt( (raw_data>>20) & 0xFF );
       etHFAsym.setType( l1t::EtSum::kAsymEtHF );

       res_->push_back(bx, etHFAsym);

       // CENTRALITY LOWER NIB

       uint32_t centLN = ((raw_data >> 28) & 0xF);

       // MHT with HF

       raw_data = block.payload()[iFrame+5];

       l1t::EtSum mhthf = l1t::EtSum();
    
       mhthf.setHwPt(raw_data & 0xFFF);
       mhthf.setHwPhi((raw_data >> 12) & 0xFF);
       mhthf.setType(l1t::EtSum::kMissingHtHF);       
       mhthf.setP4( l1t::CaloTools::p4Demux(&mhthf) );

       LogDebug("L1T") << "MHThf: phi " << mhthf.hwPhi() << " pT " << mhthf.hwPt() << " bx " << bx;

       res_->push_back(bx,mhthf);

       // HTHF ASYMM
       
       l1t::EtSum htHFAsym = l1t::EtSum();
       htHFAsym.setHwPt( (raw_data>>20) & 0xFF );
       htHFAsym.setType( l1t::EtSum::kAsymHtHF );

       res_->push_back(bx, htHFAsym);


       // CENTRALITY UPPER NIB

       uint32_t centUN = ((raw_data >> 28) & 0xF);

       l1t::EtSum cent = l1t::EtSum();
       cent.setHwPt( (centUN << 4) | centLN); 
       cent.setType( l1t::EtSum::kCentrality );
       
       res_->push_back(bx, cent);

       //HI-SUM
       
       raw_data = block.payload()[iFrame+1];

       l1t::EtSum towCount = l1t::EtSum();
       towCount.setHwPt( (raw_data>>12) & 0x1FFF );
       towCount.setType( (l1t::EtSum::kTowerCount) );
       towCount.setP4( l1t::CaloTools::p4Demux(&towCount) );

       res_->push_back(bx, towCount);
    
   
     }

     return true;
   }
}
}

DEFINE_L1T_UNPACKER(l1t::stage2::EtSumUnpacker_0x10010057);