l1t::stage2::MPUnpacker_0x1001000b Class Reference

#include <MPUnpacker_0x1001000b.h>

Inheritance diagram for l1t::stage2::MPUnpacker_0x1001000b:

Public Types
enum	{   BLK_TOT_POS =123, BLK_X_POS =121, BLK_Y_POS =127, BLK_TOT_NEG =125,   BLK_X_NEG =131, BLK_Y_NEG =129 }

Public Member Functions
bool	unpack (const Block &block, UnpackerCollections *coll) override
Public Member Functions inherited from l1t::Unpacker
unsigned int	getAlgoVersion ()
void	setAlgoVersion (const unsigned int version)
	Unpacker ()
virtual	~Unpacker ()=default

Private Member Functions
int	etaSign (int blkId)

Detailed Description

Definition at line 8 of file MPUnpacker_0x1001000b.h.

Member Enumeration Documentation

anonymous enum

Enumerator
BLK_TOT_POS
BLK_X_POS
BLK_Y_POS
BLK_TOT_NEG
BLK_X_NEG
BLK_Y_NEG

Definition at line 10 of file MPUnpacker_0x1001000b.h.

{ BLK_TOT_POS=123, BLK_X_POS=121, BLK_Y_POS=127, BLK_TOT_NEG=125, BLK_X_NEG=131, BLK_Y_NEG=129};

Member Function Documentation

int l1t::stage2::MPUnpacker_0x1001000b::etaSign ( int  blkId )

private

Definition at line 12 of file MPUnpacker_0x1001000b.cc.

References BLK_TOT_NEG, BLK_X_NEG, and BLK_Y_NEG.

Referenced by unpack().

                                              {
     if ((blkId == BLK_TOT_NEG) ||
         (blkId == BLK_X_NEG) ||
     (blkId == BLK_Y_NEG)) {
       return -1;
     }
     return 1;
   }

bool l1t::stage2::MPUnpacker_0x1001000b::unpack ( const Block &  block, UnpackerCollections *  coll  )

overridevirtual

Implements l1t::Unpacker.

Definition at line 22 of file MPUnpacker_0x1001000b.cc.

References l1t::Block::amc(), BLK_TOT_NEG, BLK_TOT_POS, BLK_X_NEG, BLK_X_POS, BLK_Y_NEG, BLK_Y_POS, l1t::CaloTools::caloEta(), coll, DEFINE_L1T_UNPACKER, ModuleToSequenceAssign::EGamma, stringResolutionProvider_cfi::et, etaSign(), l1t::BlockHeader::getID(), l1t::BlockHeader::getSize(), l1t::Block::header(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::L1Candidate::hwQual(), metsig::jet, nanoDQM_cff::Jet, l1t::EtSum::kTotalEt, l1t::EtSum::kTotalEtx, l1t::EtSum::kTotalEty, l1t::EtSum::kTotalHt, l1t::EtSum::kTotalHtx, l1t::EtSum::kTotalHty, LogDebug, l1t::CaloTools::p4MP(), l1t::Block::payload(), l1t::L1Candidate::setHwEta(), l1t::L1Candidate::setHwIso(), l1t::L1Candidate::setHwPhi(), l1t::L1Candidate::setHwPt(), l1t::L1Candidate::setHwQual(), reco::LeafCandidate::setP4(), l1t::EtSum::setType(), metsig::tau, and globals_cff::x1.

   {

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

     // check this is the correct MP
     unsigned int amc  = block.amc().getAMCNumber();
     unsigned int bxid = block.amc().getBX();
     LogDebug("L1T") << "Unpacking AMC " << amc << " for BX " << bxid;

     auto res1_ = static_cast<CaloCollections*>(coll)->getMPJets();
     auto res2_ = static_cast<CaloCollections*>(coll)->getMPEtSums();
     auto res3_ = static_cast<CaloCollections*>(coll)->getMPEGammas();
     auto res4_ = static_cast<CaloCollections*>(coll)->getMPTaus();
     
     res1_->setBXRange(0,0);
     res2_->setBXRange(0,0);
     res3_->setBXRange(0,0);
     res4_->setBXRange(0,0);

     // Initialise frame indices for each data type
     int unsigned fet  = 0;
     int unsigned fht  = 2;
     int unsigned feg  = 3;
     int unsigned ftau = 5;
     int unsigned fjet = 7;



     //      ===== Jets and Sums =====

     // ET / MET(x) / MET (y)

     uint32_t raw_data = block.payload()[fet];

     l1t::EtSum et = l1t::EtSum();

     switch(block.header().getID()){
     case BLK_TOT_POS:
       et.setType(l1t::EtSum::kTotalEt);
       et.setHwPt( ( uint32_t(raw_data & 0xFFFFFF) << 16 ) >> 16 );
       break;
     case BLK_X_POS:
       et.setType(l1t::EtSum::kTotalEtx);
       et.setHwPt( ( uint32_t(raw_data) ) );
       break;
     case BLK_Y_POS:
       et.setType(l1t::EtSum::kTotalEty);
       et.setHwPt( ( uint32_t(raw_data) ) );
       break;
     case BLK_TOT_NEG:
       et.setType(l1t::EtSum::kTotalEt);
       et.setHwPt( ( uint32_t(raw_data & 0xFFFFFF) << 16 ) >> 16 );
       break;
     case BLK_X_NEG:
       et.setType(l1t::EtSum::kTotalEtx);
       et.setHwPt( ( uint32_t(raw_data) ) );
       break;
     case BLK_Y_NEG:
       et.setType(l1t::EtSum::kTotalEty); 
       et.setHwPt( ( uint32_t(raw_data) ) );
       break;
     default: 
       break;
     }

     LogDebug("L1T") << "ET/METx/METy: pT " << et.hwPt();

     res2_->push_back(0,et);


     // HT / MHT(x)/ MHT (y)

     raw_data = block.payload()[fht];

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

     //ht.setHwPt(raw_data & 0xFFFFF);
     ht.setHwPt( ( uint32_t(raw_data & 0xFFFFFF) << 16 ) >> 16 );
     switch(block.header().getID()){
     case BLK_TOT_POS:  ht.setType(l1t::EtSum::kTotalHt);  break;
     case BLK_X_POS:  ht.setType(l1t::EtSum::kTotalHtx); break;
     case BLK_Y_POS:  ht.setType(l1t::EtSum::kTotalHty); break;
     case BLK_TOT_NEG:  ht.setType(l1t::EtSum::kTotalHt);  break;
     case BLK_X_NEG:  ht.setType(l1t::EtSum::kTotalHtx); break;
     case BLK_Y_NEG: ht.setType(l1t::EtSum::kTotalHty); break;
     default: break;
     }

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

     res2_->push_back(0,ht);

     // Two jets
     for (unsigned nJet=0; nJet < 2; nJet++){

       raw_data = block.payload()[fjet+nJet];

       if (raw_data == 0)
            continue;

       l1t::Jet jet = l1t::Jet();

       int etasign = etaSign(block.header().getID());

       LogDebug("L1") << "block ID=" << block.header().getID() << " etasign=" << etasign;

       int mpEta = etasign*(raw_data & 0x3F);
       jet.setHwEta(CaloTools::caloEta(mpEta));
       jet.setHwPhi((raw_data >> 6) & 0x7F);
       jet.setHwPt((raw_data >> 13) & 0xFFFF);
       jet.setHwQual((raw_data>>29) & 0x1 );
         
       LogDebug("L1T") << "Jet: eta " << jet.hwEta() << " phi " << jet.hwPhi() << " pT " << jet.hwPt() << " qual " << jet.hwQual();

       jet.setP4( l1t::CaloTools::p4MP(&jet) );
       res1_->push_back(0,jet);

     }

     //      ===== EGammas =====
     
     // Two EGammas


     for (unsigned nEG=0; nEG < 2; nEG++){
       
       raw_data = block.payload()[feg+nEG];
       
       if (raw_data == 0)
     continue;
       
       l1t::EGamma eg = l1t::EGamma();

       int etasign = etaSign(block.header().getID());

       LogDebug("L1") << "block ID=" << block.header().getID() << " etasign=" << etasign;

       eg.setHwEta(etasign*((raw_data >> 3) & 0x3F));
       eg.setHwPhi((raw_data >> 9) & 0x7F);
       eg.setHwPt((raw_data >> 20) & 0xFFF);
       eg.setHwIso((raw_data>>1) & 0x3);       
       eg.setHwQual((raw_data>>16) & 0xf );

       
       LogDebug("L1T") << "Egamma: eta " << eg.hwEta() << " phi " << eg.hwPhi() << " pT " << eg.hwPt() << " qual " << eg.hwQual();
       
       eg.setP4( l1t::CaloTools::p4MP(&eg) );
       res3_->push_back(0,eg);
     }

     
     //      ===== Taus =====
     
     // Two taus

     for (unsigned nTau=0; nTau < 2; nTau++){
       
       raw_data = block.payload()[ftau+nTau];
       
       if (raw_data == 0)
     continue;
       
       l1t::Tau tau = l1t::Tau();

       int etasign = etaSign(block.header().getID());
       
       LogDebug("L1") << "block ID=" << block.header().getID() << " etasign=" << etasign;
       
       tau.setHwEta(etasign*((raw_data >> 3) & 0x3F));
       tau.setHwPhi((raw_data >> 9) & 0x7F);
       tau.setHwPt((raw_data >> 20) & 0xFFF);
       tau.setHwIso((raw_data>>1) & 0x3);       
       tau.setHwQual((raw_data>>16) & 0xf );
       
       LogDebug("L1T") << "Tau: eta " << tau.hwEta() << " phi " << tau.hwPhi() << " pT " << tau.hwPt() << " qual " << tau.hwQual();
       
       tau.setP4( l1t::CaloTools::p4MP(&tau) );
       res4_->push_back(0,tau);
     }



     return true;
   }