CSCTMBHeader2013.cc

Go to the documentation of this file.

#include "EventFilter/CSCRawToDigi/interface/CSCTMBHeader2013.h"
#include "EventFilter/CSCRawToDigi/interface/CSCDMBHeader.h"
#include "FWCore/Utilities/interface/Exception.h"

CSCTMBHeader2013::CSCTMBHeader2013() {
  bzero(data(), sizeInWords() * 2);
  bits.nHeaderFrames = 42;
  bits.e0bline = 0x6E0B;
  bits.b0cline = 0xDB0C;
  bits.firmRevCode = 0x7a76;
  bits.nTBins = 12;
  bits.nCFEBs = 5;
}

CSCTMBHeader2013::CSCTMBHeader2013(const unsigned short* buf) { memcpy(data(), buf, sizeInWords() * 2); }

void CSCTMBHeader2013::setEventInformation(const CSCDMBHeader& dmbHeader) {
  bits.cscID = dmbHeader.dmbID();
  bits.l1aNumber = dmbHeader.l1a24() & 0xFFF;
  // bits.bxnCount = dmbHeader.bxn12();
}

std::vector<CSCCLCTDigi> CSCTMBHeader2013::CLCTDigis(uint32_t idlayer) {
  std::vector<CSCCLCTDigi> result;
  int halfstrip = bits.clct0_key_low + (bits.clct0_key_high << 7);
  int strip = halfstrip % CSCConstants::NUM_HALF_STRIPS_PER_CFEB;
  int cfeb = halfstrip / CSCConstants::NUM_HALF_STRIPS_PER_CFEB;
  int pattern = bits.clct0_shape;
  int bend = pattern & 0x1;

  //offlineStripNumbering(strip, cfeb, pattern, bend);
  CSCCLCTDigi digi0(
      bits.clct0_valid, bits.clct0_quality, pattern, 1, bend, strip, cfeb, bits.clct_bxn, 1, bits.bxnPreTrigger);

  halfstrip = bits.clct1_key_low + (bits.clct1_key_high << 7);
  strip = halfstrip % CSCConstants::NUM_HALF_STRIPS_PER_CFEB;
  cfeb = halfstrip / CSCConstants::NUM_HALF_STRIPS_PER_CFEB;
  pattern = bits.clct1_shape;
  bend = pattern & 0x1;

  //offlineStripNumbering(strip, cfeb, pattern, bend);
  CSCCLCTDigi digi1(
      bits.clct1_valid, bits.clct1_quality, pattern, 1, bend, strip, cfeb, bits.clct_bxn, 2, bits.bxnPreTrigger);
  result.push_back(digi0);
  result.push_back(digi1);
  return result;
}

std::vector<CSCCorrelatedLCTDigi> CSCTMBHeader2013::CorrelatedLCTDigis(uint32_t idlayer) const {
  std::vector<CSCCorrelatedLCTDigi> result;
  int strip = bits.MPC_Muon0_halfstrip_clct_pattern;  //this goes from 0-159
  //offlineHalfStripNumbering(strip);
  CSCCorrelatedLCTDigi digi(1,
                            bits.MPC_Muon0_vpf_,
                            bits.MPC_Muon0_quality_,
                            bits.MPC_Muon0_wire_,
                            strip,
                            bits.MPC_Muon0_clct_pattern_,
                            bits.MPC_Muon0_bend_,
                            bits.MPC_Muon0_bx_,
                            0,
                            bits.MPC_Muon0_bc0_,
                            bits.MPC_Muon0_SyncErr_,
                            bits.MPC_Muon0_cscid_low | (bits.MPC_Muon0_cscid_bit4 << 3));
  result.push_back(digi);
  strip = bits.MPC_Muon1_halfstrip_clct_pattern;  //this goes from 0-159
  //offlineHalfStripNumbering(strip);
  digi = CSCCorrelatedLCTDigi(2,
                              bits.MPC_Muon1_vpf_,
                              bits.MPC_Muon1_quality_,
                              bits.MPC_Muon1_wire_,
                              strip,
                              bits.MPC_Muon1_clct_pattern_,
                              bits.MPC_Muon1_bend_,
                              bits.MPC_Muon1_bx_,
                              0,
                              bits.MPC_Muon1_bc0_,
                              bits.MPC_Muon1_SyncErr_,
                              bits.MPC_Muon1_cscid_low | (bits.MPC_Muon1_cscid_bit4 << 3));
  result.push_back(digi);
  return result;
}

void CSCTMBHeader2013::addALCT0(const CSCALCTDigi& digi) {
  throw cms::Exception("In CSC TMBHeaderFormat 2007, ALCTs belong in  ALCT header");
}

void CSCTMBHeader2013::addALCT1(const CSCALCTDigi& digi) {
  throw cms::Exception("In CSC TMBHeaderFormat 2007, ALCTs belong in  ALCT header");
}

void CSCTMBHeader2013::addCLCT0(const CSCCLCTDigi& digi) {
  int halfStrip = digi.getKeyStrip();
  int pattern = digi.getPattern();
  //int bend = digi.getBend();
  //hardwareStripNumbering(strip, cfeb, pattern, bend);
  bits.clct0_valid = digi.isValid();
  bits.clct0_quality = digi.getQuality();
  bits.clct0_shape = pattern;
  // first 7 bits of halfstrip
  bits.clct0_key_low = halfStrip & (0x7F);
  // most-significant (8th) bit
  bits.clct0_key_high = (halfStrip >> 7) & (0x1);
  bits.clct_bxn = digi.getBX();
  bits.bxnPreTrigger = digi.getFullBX();
  bits.bxnCount = (digi.getFullBX() + 167) & 0xFFF;
}

void CSCTMBHeader2013::addCLCT1(const CSCCLCTDigi& digi) {
  int halfStrip = digi.getKeyStrip();
  int pattern = digi.getPattern();
  //int bend = digi.getBend();
  //hardwareStripNumbering(strip, cfeb, pattern, bend);
  bits.clct1_valid = digi.isValid();
  bits.clct1_quality = digi.getQuality();
  bits.clct1_shape = pattern;
  // first 7 bits of halfstrip
  bits.clct1_key_low = halfStrip & (0x7F);
  // most-significant (8th) bit
  bits.clct1_key_high = (halfStrip >> 7) & (0x1);
  // There is just one BX field common for CLCT0 and CLCT1 (since both
  // are latched at the same BX); set it in addCLCT0().
  //bits.clct_bxn = digi.getBX();
  bits.bxnPreTrigger = digi.getFullBX();
}

void CSCTMBHeader2013::addCorrelatedLCT0(const CSCCorrelatedLCTDigi& digi) {
  int halfStrip = digi.getStrip();
  //hardwareHalfStripNumbering(halfStrip);

  bits.MPC_Muon0_vpf_ = digi.isValid();
  bits.MPC_Muon0_wire_ = digi.getKeyWG();
  bits.MPC_Muon0_clct_pattern_ = digi.getPattern();
  bits.MPC_Muon0_quality_ = digi.getQuality();
  bits.MPC_Muon0_halfstrip_clct_pattern = halfStrip;
  bits.MPC_Muon0_bend_ = digi.getBend();
  bits.MPC_Muon0_SyncErr_ = digi.getSyncErr();
  bits.MPC_Muon0_bx_ = digi.getBX();
  bits.MPC_Muon0_bc0_ = digi.getBX0();
  bits.MPC_Muon0_cscid_low = digi.getCSCID() & 0x7;
  bits.MPC_Muon0_cscid_bit4 = (digi.getCSCID() >> 3) & 0x1;
}

void CSCTMBHeader2013::addCorrelatedLCT1(const CSCCorrelatedLCTDigi& digi) {
  int halfStrip = digi.getStrip();
  //hardwareHalfStripNumbering(halfStrip);

  bits.MPC_Muon1_vpf_ = digi.isValid();
  bits.MPC_Muon1_wire_ = digi.getKeyWG();
  bits.MPC_Muon1_clct_pattern_ = digi.getPattern();
  bits.MPC_Muon1_quality_ = digi.getQuality();
  bits.MPC_Muon1_halfstrip_clct_pattern = halfStrip;
  bits.MPC_Muon1_bend_ = digi.getBend();
  bits.MPC_Muon1_SyncErr_ = digi.getSyncErr();
  bits.MPC_Muon1_bx_ = digi.getBX();
  bits.MPC_Muon1_bc0_ = digi.getBX0();
  bits.MPC_Muon1_cscid_low = digi.getCSCID() & 0x7;
  bits.MPC_Muon1_cscid_bit4 = (digi.getCSCID() >> 3) & 0x1;
}

void CSCTMBHeader2013::print(std::ostream& os) const {
  os << "...............TMB Header.................."
     << "\n";
  os << std::hex << "BOC LINE " << bits.b0cline << " EOB " << bits.e0bline << "\n";
  os << std::dec << "fifoMode = " << bits.fifoMode << ", nTBins = " << bits.nTBins << "\n";
  os << "boardID = " << bits.boardID << ", cscID = " << bits.cscID << "\n";
  os << "l1aNumber = " << bits.l1aNumber << ", bxnCount = " << bits.bxnCount << "\n";
  os << "trigSourceVect = " << bits.trigSourceVect
     << ", activeCFEBs = " << (bits.activeCFEBs | (bits.activeCFEBs_2 << 5))
     << ", readCFEBs = " << (bits.readCFEBs | (bits.readCFEBs_2 << 5)) << "\n";
  os << "bxnPreTrigger = " << bits.bxnPreTrigger << "\n";
  os << "tmbMatch = " << bits.tmbMatch << " alctOnly = " << bits.alctOnly << " clctOnly = " << bits.clctOnly << "\n";
  os << " bits.clct0_shape = " << bits.clct0_shape << " clct0_quality = " << bits.clct0_quality << "\n";

  os << "..................CLCT....................."
     << "\n";
}