l1t::stage2::emtf::RPCBlockUnpacker Class Reference

Inheritance diagram for l1t::stage2::emtf::RPCBlockUnpacker:

Public Member Functions
virtual int	checkFormat (const Block &block)
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

Detailed Description

Definition at line 13 of file EMTFBlockRPC.cc.

Member Function Documentation

int l1t::stage2::emtf::RPCBlockUnpacker::checkFormat ( const Block &  block )

virtual

Definition at line 34 of file EMTFBlockRPC.cc.

References benchmark_cfg::errors, l1t::stage2::emtf::GetHexBits(), mps_fire::i, jets_cff::payload, and l1t::Block::payload().

Referenced by unpack().

                                                          {                 
    
    auto payload = block.payload();
    int errors = 0;
    
    // Check the number of 16-bit words
    if (payload.size() != 4) { errors += 1; 
      edm::LogError("L1T|EMTF") << "Payload size in 'RPC Data Record' is different than expected"; }
    
    // Check that each word is 16 bits
    for (unsigned int i = 0; i < 4; i++) {
      if (GetHexBits(payload[i], 16, 31) != 0) { errors += 1;
        edm::LogError("L1T|EMTF") << "Payload[" << i << "] has more than 16 bits in 'RPC Data Record'"; }
    }
    
    uint16_t RPCa = payload[0];
    uint16_t RPCb = payload[1];
    uint16_t RPCc = payload[2];
    uint16_t RPCd = payload[3];
    
    // Check Format
    if (GetHexBits(RPCa, 15, 15) != 0) { errors += 1;
      edm::LogError("L1T|EMTF") << "Format identifier bits in RPCa are incorrect"; }
    if (GetHexBits(RPCb, 15, 15) != 0) { errors += 1;
      edm::LogError("L1T|EMTF") << "Format identifier bits in RPCb are incorrect"; }
    if (GetHexBits(RPCc, 15, 15) != 1) { errors += 1;
      edm::LogError("L1T|EMTF") << "Format identifier bits in RPCc are incorrect"; }
    if (GetHexBits(RPCd, 15, 15) != 0) { errors += 1;
      edm::LogError("L1T|EMTF") << "Format identifier bits in RPCd are incorrect"; }

    return errors;
    
      }

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

overridevirtual

Implements l1t::Unpacker.

Definition at line 97 of file EMTFBlockRPC.cc.

References l1t::emtf::RPC::add_format_error(), l1t::EMTFHit::BX(), l1t::EMTFHit::Chamber(), checkFormat(), coll, l1t::stage2::emtf::convert_RPC_location(), l1t::EMTFHit::CreateCPPFDigi(), DEFINE_L1T_UNPACKER, l1t::EMTFHit::Endcap(), l1t::emtf::RPC::Frame(), l1t::stage2::emtf::GetHexBits(), l1t::stage2::emtf::ImportRPC(), l1t::emtf::RPC::Link(), l1t::EMTFHit::Neighbor(), jets_cff::payload, l1t::Block::payload(), l1t::EMTFHit::Phi_fp(), l1t::EMTFHit::Ring(), l1t::EMTFHit::Sector(), l1t::EMTFHit::Sector_RPC(), l1t::emtf::RPC::set_bc0(), l1t::EMTFHit::set_bt_segment(), l1t::EMTFHit::set_chamber(), l1t::emtf::RPC::set_frame(), l1t::EMTFHit::set_fs_segment(), l1t::emtf::RPC::set_link(), l1t::EMTFHit::set_neighbor(), l1t::EMTFHit::set_pc_segment(), l1t::emtf::RPC::set_phi(), l1t::EMTFHit::set_ring(), l1t::emtf::RPC::set_rpc_bxn(), l1t::EMTFHit::set_sector(), l1t::EMTFHit::set_sector_RPC(), l1t::EMTFHit::set_station(), l1t::EMTFHit::set_stub_num(), l1t::EMTFHit::set_subsector(), l1t::EMTFHit::set_subsector_RPC(), l1t::emtf::RPC::set_tbin(), l1t::emtf::RPC::set_theta(), l1t::emtf::RPC::set_vp(), l1t::emtf::RPC::set_word(), l1t::EMTFHit::Station(), l1t::EMTFHit::Stub_num(), l1t::EMTFHit::Subsector_RPC(), l1t::EMTFHit::Theta_fp(), and l1t::emtf::RPC::Word().

                                                                                 {
    
    // std::cout << "Inside EMTFBlockRPC.cc: unpack" << std::endl;
    
    // Get the payload for this block, made up of 16-bit words (0xffff)
    // Format defined in MTF7Payload::getBlock() in src/Block.cc
    // payload[0] = bits 0-15, payload[1] = 16-31, payload[3] = 32-47, etc.
    auto payload = block.payload();
    
    // Check Format of Payload
    l1t::emtf::RPC RPC_;    
    for (int err = 0; err < checkFormat(block); err++) RPC_.add_format_error();
    
    // Assign payload to 16-bit words
    uint16_t RPCa = payload[0];
    uint16_t RPCb = payload[1];
    uint16_t RPCc = payload[2];
    uint16_t RPCd = payload[3];
    
    // res is a pointer to a collection of EMTFDaqOut class objects
    // There is one EMTFDaqOut for each MTF7 (60 deg. sector) in the event
    EMTFDaqOutCollection* res;
    res = static_cast<EMTFCollections*>(coll)->getEMTFDaqOuts();
    int iOut = res->size() - 1;

        EMTFHitCollection* res_hit;
        res_hit = static_cast<EMTFCollections*>(coll)->getEMTFHits();
        EMTFHit Hit_;

        CPPFDigiCollection* res_CPPF;
        res_CPPF = static_cast<EMTFCollections*>(coll)->getEMTFCPPFs();

    // Unpack the RPC Data Record
    
    RPC_.set_phi     ( GetHexBits(RPCa,  0, 10) );

    RPC_.set_theta   ( GetHexBits(RPCb,  0,  4) );
    RPC_.set_word    ( GetHexBits(RPCb,  8,  9) );
    RPC_.set_frame   ( GetHexBits(RPCb, 10, 11) );
    RPC_.set_link    ( GetHexBits(RPCb, 12, 14) ); // Link index (0 - 6); link number runs 1 - 7
    
    RPC_.set_rpc_bxn ( GetHexBits(RPCc,  0, 11) );
    RPC_.set_bc0     ( GetHexBits(RPCc, 14, 14) );
    
    RPC_.set_tbin    ( GetHexBits(RPCd,  0,  2) );
    RPC_.set_vp      ( GetHexBits(RPCd,  3,  3) );
    
    // RPC_.set_dataword            ( uint64_t dataword);


    // Convert specially-encoded RPC quantities
    int _station, _ring, _sector, _subsector, _neighbor, _segment;
    convert_RPC_location( _station, _ring, _sector, _subsector, _neighbor, _segment,
                  (res->at(iOut)).PtrEventHeader()->Sector(), RPC_.Frame(), RPC_.Word(), RPC_.Link() );

    Hit_.set_station       ( _station   );
    Hit_.set_ring          ( _ring      );
    Hit_.set_sector        ( _sector    );
    Hit_.set_subsector     ( _subsector );
    Hit_.set_sector_RPC    ( _subsector < 5 ? _sector : (_sector % 6) + 1);  // Rotate by 20 deg to match RPC convention in CMSSW
    Hit_.set_subsector_RPC ( ((_subsector + 1) % 6) + 1 );  // Rotate by 2 to match RPC convention in CMSSW (RPCDetId.h) 
    Hit_.set_chamber       ( (Hit_.Sector_RPC() - 1)*6 + Hit_.Subsector_RPC() );
    Hit_.set_neighbor      ( _neighbor  );
    Hit_.set_pc_segment    ( _segment   );
    Hit_.set_fs_segment    ( _segment   );
    Hit_.set_bt_segment    ( _segment   );

    // Fill the EMTFHit
    ImportRPC( Hit_, RPC_, (res->at(iOut)).PtrEventHeader()->Endcap(), (res->at(iOut)).PtrEventHeader()->Sector() );

    // Set the stub number for this hit
    // Each chamber can send up to 2 stubs per BX
    // Also count stubs in corresponding CSC chamber; RPC hit counting is on top of LCT counting
    Hit_.set_stub_num(0);
    // See if matching hit is already in event record
    bool exact_duplicate = false;
    for (auto const & iHit : *res_hit) {
      
      if ( Hit_.BX()      == iHit.BX()      && 
           Hit_.Endcap()  == iHit.Endcap()  &&
           Hit_.Station() == iHit.Station() &&
           Hit_.Chamber() == iHit.Chamber() ) {

        if ( (iHit.Is_CSC() == 1 && iHit.Ring() == 2) ||
         (iHit.Is_RPC() == 1) ) { // RPC rings 2 and 3 both map to CSC ring 2
          if ( Hit_.Neighbor() == iHit.Neighbor() ) {
        Hit_.set_stub_num( Hit_.Stub_num() + 1);
        if ( iHit.Is_RPC() == 1                 &&
             iHit.Ring()     == Hit_.Ring()     &&
             iHit.Theta_fp() == Hit_.Theta_fp() &&
             iHit.Phi_fp()   == Hit_.Phi_fp()   ) {
          exact_duplicate = true;
        }
          }
        }
      }
    } // End loop: for (auto const & iHit : *res_hit)

        if (exact_duplicate) edm::LogWarning("L1T|EMTF") << "EMTF unpacked duplicate CPPF digis: BX " << Hit_.BX()
                                                         << ", endcap " << Hit_.Endcap() << ", station " << Hit_.Station()
                                                         << ", sector " << Hit_.Sector() << ", neighbor " << Hit_.Neighbor()
                                                         << ", ring " << Hit_.Ring() << ", chamber " << Hit_.Chamber()
                                                         << ", theta " << Hit_.Theta_fp() / 4 << ", phi " << Hit_.Phi_fp() / 4 << std::endl;


    (res->at(iOut)).push_RPC(RPC_);
    if (!exact_duplicate)
      res_hit->push_back(Hit_);
    if (!exact_duplicate)
      res_CPPF->push_back( Hit_.CreateCPPFDigi() );
    
    // Finished with unpacking one RPC Data Record
    return true;
    
      } // End bool RPCBlockUnpacker::unpack