l1t::stage2::EtSumUnpacker Class Reference

#include <EtSumUnpacker.h>

Inheritance diagram for l1t::stage2::EtSumUnpacker:

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

Private Attributes
unsigned int	EtSumCopy_

Detailed Description

Definition at line 8 of file EtSumUnpacker.h.

Constructor & Destructor Documentation

l1t::stage2::EtSumUnpacker::EtSumUnpacker

(

)

Definition at line 13 of file EtSumUnpacker.cc.

                                : EtSumCopy_(0)
   {
   }

l1t::stage2::EtSumUnpacker::~EtSumUnpacker ( )

inlineoverride

Definition at line 11 of file EtSumUnpacker.h.

References groupFilesInBlocks::block, coll, and unpack().

{};

Member Function Documentation

void l1t::stage2::EtSumUnpacker::setEtSumCopy ( const unsigned int  copy )

inline

Definition at line 15 of file EtSumUnpacker.h.

References popcon2dropbox::copy(), and EtSumCopy_.

{ EtSumCopy_ = copy; };

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

overridevirtual

Implements l1t::Unpacker.

Definition at line 18 of file EtSumUnpacker.cc.

References coll, DEFINE_L1T_UNPACKER, stringResolutionProvider_cfi::et, EtSumCopy_, l1t::BlockHeader::getID(), l1t::BlockHeader::getSize(), l1t::Block::header(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), createfilelist::int, l1t::EtSum::kMinBiasHFM0, l1t::EtSum::kMinBiasHFM1, l1t::EtSum::kMinBiasHFP0, l1t::EtSum::kMinBiasHFP1, l1t::EtSum::kMissingEt, l1t::EtSum::kMissingEtHF, l1t::EtSum::kMissingHt, l1t::EtSum::kMissingHtHF, l1t::EtSum::kTotalEt, l1t::EtSum::kTotalEtEm, l1t::EtSum::kTotalHt, l1t::EtSum::kTowerCount, LogDebug, RazorAnalyzer::met, l1t::stage2::layer2::mp::nOutputFramePerBX, l1t::CaloTools::p4Demux(), l1t::Block::payload(), l1t::L1Candidate::setHwPhi(), l1t::L1Candidate::setHwPt(), reco::LeafCandidate::setP4(), and l1t::EtSum::setType().

Referenced by ~EtSumUnpacker().

   {

     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);

       // 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);

       
       //  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);

       // 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);

       //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;
   }

Member Data Documentation

unsigned int l1t::stage2::EtSumUnpacker::EtSumCopy_

private

Definition at line 15 of file EtSumUnpacker.h.

Referenced by setEtSumCopy(), and unpack().