EMTFBlockME.cc

Go to the documentation of this file.

// Code to unpack the "ME Data Record"

#include "EventFilter/L1TRawToDigi/plugins/UnpackerFactory.h"

#include "EMTFCollections.h"
#include "EMTFUnpackerTools.h"

// This is the "header" - no EMTFBlockME.h file is needed
namespace l1t {
  namespace stage2 {
    namespace emtf {
      
      class MEBlockUnpacker : public Unpacker { // "MEBlockUnpacker" inherits from "Unpacker"
      public:
    virtual int  checkFormat(const Block& block);
    virtual bool unpack(const Block& block, UnpackerCollections *coll) override; // Apparently it's always good to use override in C++
    // virtual bool packBlock(const Block& block, UnpackerCollections *coll) override;
      };
      
      // class MEBlockPacker : public Packer { // "MEBlockPacker" inherits from "Packer"
      // public:
      //    virtual bool unpack(const Block& block, UnpackerCollections *coll) override; // Apparently it's always good to use override in C++
      // };
      
    }
  }
}

namespace l1t {
  namespace stage2 {
    namespace emtf {

      int MEBlockUnpacker::checkFormat(const Block& block) {
    
    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 'ME 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 'ME Data Record'"; }
    }
    
    uint16_t MEa = payload[0];
    uint16_t MEb = payload[1];
    uint16_t MEc = payload[2];
    uint16_t MEd = payload[3];

    //Check Format
    if (GetHexBits(MEa, 15, 15) != 1) { errors += 1; 
      edm::LogError("L1T|EMTF") << "Format identifier bits in MEa are incorrect"; }
    if (GetHexBits(MEb, 15, 15) != 1) { errors += 1; 
      edm::LogError("L1T|EMTF") << "Format identifier bits in MEb are incorrect"; }
    if (GetHexBits(MEc, 15, 15) != 0) { errors += 1; 
      edm::LogError("L1T|EMTF") << "Format identifier bits in MEc are incorrect"; }
    if (GetHexBits(MEd, 15, 15) != 0) { errors += 1; 
      edm::LogError("L1T|EMTF") << "Format identifier bits in MEd are incorrect"; }

    return errors;

      }

      // Converts station, CSC_ID, sector, subsector, and neighbor from the ME output
      std::vector<int> convert_ME_location(int _station, int _csc_ID, int _sector) {
    int new_sector = _sector;
    int new_csc_ID = _csc_ID; // Before FW update on 05.05.16, shift by +1 from 0,1,2... convention to 1,2,3...
    if      (_station == 0) { int arr[] = {       1, new_csc_ID, new_sector,  1, 0}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (_station == 1) { int arr[] = {       1, new_csc_ID, new_sector,  2, 0}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (_station <= 4) { int arr[] = {_station, new_csc_ID, new_sector, -1, 0}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (_station == 5) new_sector = (_sector != 1) ? _sector-1 : 6;
    else { int arr[] = {_station, _csc_ID, _sector, -99, -99}; std::vector<int> vec(arr, arr+5); return vec; }
    
    if      (new_csc_ID == 1) { int arr[] = {1, 3, new_sector,  2, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 2) { int arr[] = {1, 6, new_sector,  2, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 3) { int arr[] = {1, 9, new_sector,  2, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 4) { int arr[] = {2, 3, new_sector, -1, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 5) { int arr[] = {2, 9, new_sector, -1, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 6) { int arr[] = {3, 3, new_sector, -1, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 7) { int arr[] = {3, 9, new_sector, -1, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 8) { int arr[] = {4, 3, new_sector, -1, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else if (new_csc_ID == 9) { int arr[] = {4, 9, new_sector, -1, 1}; std::vector<int> vec(arr, arr+5); return vec; }
    else                   { int arr[] = {_station, _csc_ID, _sector, -99, -99}; std::vector<int> vec(arr, arr+5); return vec; }
      }


      bool MEBlockUnpacker::unpack(const Block& block, UnpackerCollections *coll) {

    // std::cout << "Inside EMTFBlockME.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();

    // Assign payload to 16-bit words
        uint16_t MEa = payload[0];
        uint16_t MEb = payload[1];
        uint16_t MEc = payload[2];
        uint16_t MEd = payload[3];

    // Check Format of Payload
    l1t::emtf::ME ME_;
    for (int err = 0; err < checkFormat(block); err++) ME_.add_format_error();

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

    CSCCorrelatedLCTDigiCollection* res_LCT;
    res_LCT = static_cast<EMTFCollections*>(coll)->getEMTFLCTs();

    // Unpack the ME Data Record

    ME_.set_clct_pattern ( GetHexBits(MEa,  0,  3) );
    ME_.set_quality      ( GetHexBits(MEa,  4,  7) );
    ME_.set_wire         ( GetHexBits(MEa,  8, 14) );

    ME_.set_strip        ( GetHexBits(MEb,  0,  7) );
    ME_.set_csc_ID       ( GetHexBits(MEb,  8, 11) );
    ME_.set_lr           ( GetHexBits(MEb, 12, 12) );
    ME_.set_bxe          ( GetHexBits(MEb, 13, 13) );
    ME_.set_bc0          ( GetHexBits(MEb, 14, 14) );

    ME_.set_me_bxn       ( GetHexBits(MEc,  0, 11) );
    ME_.set_nit          ( GetHexBits(MEc, 12, 12) );
    ME_.set_cik          ( GetHexBits(MEc, 13, 13) );
    ME_.set_afff         ( GetHexBits(MEc, 14, 14) );

    ME_.set_tbin         ( GetHexBits(MEd,  0,  2) );
    ME_.set_vp           ( GetHexBits(MEd,  3,  3) );
    ME_.set_station      ( GetHexBits(MEd,  4,  6) );
    ME_.set_af           ( GetHexBits(MEd,  7,  7) );
    ME_.set_epc          ( GetHexBits(MEd,  8, 11) );
    ME_.set_sm           ( GetHexBits(MEd, 12, 12) );
    ME_.set_se           ( GetHexBits(MEd, 13, 13) );
    ME_.set_afef         ( GetHexBits(MEd, 14, 14) );

    // ME_.set_dataword     ( uint64_t dataword);

    // Convert specially-encoded ME quantities
    std::vector<int> conv_vals = convert_ME_location( ME_.Station(), ME_.CSC_ID(), 
                              (res->at(iOut)).PtrEventHeader()->Sector() );
    Hit_.set_station   ( conv_vals.at(0) );
    Hit_.set_csc_ID    ( conv_vals.at(1) );
    Hit_.set_sector    ( conv_vals.at(2) );
    Hit_.set_subsector ( conv_vals.at(3) );
    Hit_.set_neighbor  ( conv_vals.at(4) );

    // Fill the EMTFHit
    ImportME( Hit_, ME_, (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
    ME_.set_stub_num(0);
    Hit_.set_stub_num(0);
    // See if matching hit is already in event record (from neighboring sector)
    bool duplicate_hit_exists = false;
    for (unsigned int iHit = 0; iHit < res_hit->size(); iHit++) {

      if ( res_hit->at(iHit).Is_CSC() == 1                     && 
           Hit_.BX()         == res_hit->at(iHit).BX()         && 
           Hit_.Endcap()     == res_hit->at(iHit).Endcap()     &&
           Hit_.Station()    == res_hit->at(iHit).Station()    &&
           Hit_.Chamber()    == res_hit->at(iHit).Chamber()    &&
           (Hit_.Ring() % 3) == (res_hit->at(iHit).Ring() % 3) ) { // ME1/1a and ME1/1b (rings "4" and 1) are the same chamber
        
        if ( Hit_.Neighbor() == res_hit->at(iHit).Neighbor() ) {
          ME_.set_stub_num( ME_.Stub_num() + 1 );
          Hit_.set_stub_num( Hit_.Stub_num() + 1); }
        else if ( Hit_.Ring()  == res_hit->at(iHit).Ring()  && 
              Hit_.Strip() == res_hit->at(iHit).Strip() && 
              Hit_.Wire()  == res_hit->at(iHit).Wire()  )
          duplicate_hit_exists = true;
      }
    } // End loop: for (uint iHit = 0; iHit < res_hit->size(); iHit++)

    (res->at(iOut)).push_ME(ME_);
    res_hit->push_back(Hit_);
    if (not duplicate_hit_exists) // Don't write duplicate LCTs from adjacent sectors
      res_LCT->insertDigi( Hit_.CSC_DetId(), Hit_.CSC_LCTDigi() );

    // Finished with unpacking one ME Data Record
    return true;
    
      } // End bool MEBlockUnpacker::unpack

      // bool MEBlockPacker::pack(const Block& block, UnpackerCollections *coll) {
      //    std::cout << "Inside MEBlockPacker::pack" << std::endl;
      //    return true;
      // } // End bool MEBlockPacker::pack

    } // End namespace emtf
  } // End namespace stage2
} // End namespace l1t

DEFINE_L1T_UNPACKER(l1t::stage2::emtf::MEBlockUnpacker);
// DEFINE_L1T_PACKER(l1t::stage2::MEBlockPacker);