CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Static Public Member Functions | Static Public Attributes | Private Member Functions | Private Attributes
CSCEventData Class Reference

#include <CSCEventData.h>

Public Member Functions

void add (const CSCStripDigi &, int layer)
 routines to add digis to the data More...
 
void add (const CSCWireDigi &, int layer)
 
void add (const CSCComparatorDigi &, int layer)
 
void add (const CSCComparatorDigi &, const CSCDetId &)
 
void add (const std::vector< CSCALCTDigi > &)
 these go in as vectors, so they get sorted right away More...
 
void add (const std::vector< CSCCLCTDigi > &)
 
void add (const std::vector< CSCCorrelatedLCTDigi > &)
 
void add (const std::vector< CSCShowerDigi > &)
 
void add (const std::vector< GEMPadDigiCluster > &, const GEMDetId &)
 
void addALCTStructures ()
 adds an empty ALCTHeader, trailer, and anode data More...
 
CSCAnodeDataalctData () const
 user must check if nalct > 0 More...
 
CSCALCTHeaderalctHeader () const
 user must check if nalct > 0 More...
 
CSCALCTTraileralctTrailer () const
 user must check if nalct > 0 More...
 
unsigned int calcALCTcrc (std::vector< std::pair< unsigned int, unsigned short * > > &vec)
 
const CSCCFEBDatacfebData (unsigned icfeb) const
 unpacked in long mode: has overflow and error bits decoded More...
 
int chamberType () const
 might not be set in real data More...
 
CSCComparatorDatacomparatorData () const
 user must check if nclct > 0 More...
 
 CSCEventData (int chamberType, uint16_t format_version=2005)
 
 CSCEventData (const uint16_t *buf, uint16_t format_version=2005)
 should make const input soon More...
 
 CSCEventData ()
 
 CSCEventData (const CSCEventData &data)
 
const CSCDMBHeaderdmbHeader () const
 the DAQ motherboard header. A good place for event and chamber info More...
 
CSCDMBHeaderdmbHeader ()
 
const CSCDMBTrailerdmbTrailer () const
 DMB trailer. More...
 
uint16_t getFormatVersion () const
 
bool isALCT (const uint16_t *buf)
 if dealing with ALCT data More...
 
bool isTMB (const uint16_t *buf)
 if dealing with TMB data More...
 
int nalct () const
 the flag for existence of ALCT data More...
 
int nclct () const
 the number of CLCTs More...
 
CSCEventData operator= (const CSCEventData &data)
 
boost::dynamic_bitset pack ()
 returns the packed event data. More...
 
void setEventInformation (int bxnum, int lvl1num)
 
unsigned short size () const
 size of the data buffer used, in bytes More...
 
std::vector< CSCStripDigistripDigis (const CSCDetId &idlayer) const
 returns all the strip digis in the chamber, with the comparator information. More...
 
std::vector< CSCStripDigistripDigis (unsigned idlayer, unsigned icfeb) const
 returns all the strip digis in the chamber's cfeb More...
 
std::vector< std::vector
< CSCStripDigi > > 
stripDigis () const
 deprecated. Use the above methods instead More...
 
CSCTMBDatatmbData () const
 user must check in nclct > 0 More...
 
CSCTMBHeadertmbHeader () const
 user must check if nclct > 0 More...
 
std::vector< CSCWireDigiwireDigis (unsigned ilayer) const
 
std::vector< std::vector
< CSCWireDigi > > 
wireDigis () const
 deprecated. Use the above method instead. More...
 
 ~CSCEventData ()
 

Static Public Member Functions

static void selfTest ()
 
static void setDebug (const bool value)
 

Static Public Attributes

static std::atomic< bool > debug {false}
 

Private Member Functions

void addComparatorInformation (std::vector< CSCStripDigi > &, int layer) const
 adds the comparators to the strip digis More...
 
void checkALCTClasses ()
 makes new ALCT classes, if needed More...
 
void checkTMBClasses ()
 makes new TMB classes, if needed More...
 
void copy (const CSCEventData &)
 
void destroy ()
 
void init ()
 
void unpack_data (const uint16_t *buf)
 

Private Attributes

unsigned short * alctZSErecovered
 Auxiliary bufer to recove the ALCT raw payload from zero suppression. More...
 
int size_
 
CSCALCTHeadertheALCTHeader
 
CSCALCTTrailertheALCTTrailer
 
CSCAnodeDatatheAnodeData
 
CSCCFEBDatatheCFEBData [CSCConstants::MAX_CFEBS_RUN2]
 for up to MAX_CFEB CFEB boards More...
 
int theChamberType
 
CSCDMBHeader theDMBHeader
 
CSCDMBTrailer theDMBTrailer
 
uint16_t theFormatVersion
 Output Format Version (2005, 2013) More...
 
CSCTMBDatatheTMBData
 
int zseEnable
 

Detailed Description

Definition at line 26 of file CSCEventData.h.

Constructor & Destructor Documentation

CSCEventData::CSCEventData ( int  chamberType,
uint16_t  format_version = 2005 
)
explicit

Definition at line 18 of file CSCEventData.cc.

References mps_fire::i, CSCConstants::MAX_CFEBS_RUN2, and theCFEBData.

19  : theDMBHeader(format_version),
20  theALCTHeader(nullptr),
21  theAnodeData(nullptr),
22  theALCTTrailer(nullptr),
23  theTMBData(nullptr),
24  theDMBTrailer(format_version),
26  alctZSErecovered(nullptr),
27  zseEnable(0),
28  theFormatVersion(format_version) {
29  for (unsigned i = 0; i < CSCConstants::MAX_CFEBS_RUN2; ++i) {
30  theCFEBData[i] = nullptr;
31  }
32 }
unsigned short * alctZSErecovered
Auxiliary bufer to recove the ALCT raw payload from zero suppression.
Definition: CSCEventData.h:170
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
CSCDMBTrailer theDMBTrailer
Definition: CSCEventData.h:162
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
int chamberType() const
might not be set in real data
Definition: CSCEventData.h:119
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
CSCEventData::CSCEventData ( const uint16_t *  buf,
uint16_t  format_version = 2005 
)

should make const input soon

Definition at line 34 of file CSCEventData.cc.

References theFormatVersion, and unpack_data().

34  : theFormatVersion(format_version) {
35  theFormatVersion = format_version;
37 }
void unpack_data(const uint16_t *buf)
Definition: CSCEventData.cc:39
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
CSCEventData::CSCEventData ( )
inline

Definition at line 31 of file CSCEventData.h.

31 {}
CSCEventData::~CSCEventData ( )

since we need deep copies, need the Big Three (destructor, copy ctor, op=)

Definition at line 260 of file CSCEventData.cc.

References destroy().

260 { destroy(); }
CSCEventData::CSCEventData ( const CSCEventData data)

Definition at line 258 of file CSCEventData.cc.

References copy().

258 { copy(data); }
void copy(const CSCEventData &)

Member Function Documentation

void CSCEventData::add ( const CSCStripDigi digi,
int  layer 
)

routines to add digis to the data

Definition at line 461 of file CSCEventData.cc.

References CSCCFEBData::add(), CSCDMBHeader::addCFEB(), CSCDMBHeader::format_version(), CSCStripDigi::getADCCounts(), CSCStripDigi::getCFEB(), theCFEBData, theDMBHeader, and theFormatVersion.

Referenced by CSCDigiToRaw::add(), counter.Counter::register(), SequenceTypes.Task::remove(), SequenceTypes.Task::replace(), and selfTest().

461  {
462  //@@ need special logic here for ME11
463  unsigned cfeb = digi.getCFEB();
464  bool sixteenSamples = false;
465  if (digi.getADCCounts().size() == 16)
466  sixteenSamples = true;
467  if (theCFEBData[cfeb] == nullptr) {
468  bool isDCFEB = false;
469  if (theDMBHeader.format_version() == 2)
470  isDCFEB = true;
471  theCFEBData[cfeb] = new CSCCFEBData(cfeb, sixteenSamples, theFormatVersion, isDCFEB);
472  theDMBHeader.addCFEB(cfeb);
473  }
474  theCFEBData[cfeb]->add(digi, layer);
475 }
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
int getCFEB() const
Get the CFEB number. Counts from 0.
Definition: CSCStripDigi.cc:25
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
std::vector< int > const & getADCCounts() const
Get ADC readings.
Definition: CSCStripDigi.h:47
void add(const CSCStripDigi &, int layer)
Definition: CSCCFEBData.cc:82
constexpr std::array< uint8_t, layerIndexSize > layer
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
void addCFEB(int icfeb)
Definition: CSCDMBHeader.h:23
unsigned format_version() const
Definition: CSCDMBHeader.h:49
void CSCEventData::add ( const CSCWireDigi digi,
int  layer 
)

Definition at line 477 of file CSCEventData.cc.

References CSCAnodeData::add(), checkALCTClasses(), CSCWireDigi::getWireGroup(), CSCWireDigi::getWireGroupBX(), CSCALCTHeader::setBXNCount(), CSCALCTHeader::setDAVForChannel(), theALCTHeader, and theAnodeData.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

477  {
479  theAnodeData->add(digi, layer);
482 }
void checkALCTClasses()
makes new ALCT classes, if needed
void setBXNCount(unsigned int bxn)
Definition: CSCALCTHeader.h:78
constexpr std::array< uint8_t, layerIndexSize > layer
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
void add(const CSCWireDigi &wireDigi, int layer)
Definition: CSCAnodeData.h:24
int getWireGroup() const
default
Definition: CSCWireDigi.h:22
int getWireGroupBX() const
return BX assigned for the wire group (16 upper bits from the wire group number)
Definition: CSCWireDigi.h:24
void setDAVForChannel(int wireGroup)
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
void CSCEventData::add ( const CSCComparatorDigi digi,
int  layer 
)

Definition at line 484 of file CSCEventData.cc.

References CSCComparatorData::add(), checkTMBClasses(), CSCTMBData::comparatorData(), and theTMBData.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

484  {
485  checkTMBClasses();
486  theTMBData->comparatorData()->add(digi, layer);
487 }
void checkTMBClasses()
makes new TMB classes, if needed
constexpr std::array< uint8_t, layerIndexSize > layer
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCComparatorData * comparatorData()
Definition: CSCTMBData.h:41
void add(const CSCComparatorDigi &digi, int layer)
TODO for packing. Doesn&#39;t do flipping yet.
void CSCEventData::add ( const CSCComparatorDigi digi,
const CSCDetId cid 
)

Definition at line 489 of file CSCEventData.cc.

References CSCComparatorData::add(), checkTMBClasses(), CSCTMBData::comparatorData(), and theTMBData.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

489  {
490  checkTMBClasses();
491  theTMBData->comparatorData()->add(digi, cid);
492 }
void checkTMBClasses()
makes new TMB classes, if needed
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCComparatorData * comparatorData()
Definition: CSCTMBData.h:41
void add(const CSCComparatorDigi &digi, int layer)
TODO for packing. Doesn&#39;t do flipping yet.
void CSCEventData::add ( const std::vector< CSCALCTDigi > &  digis)

these go in as vectors, so they get sorted right away

Definition at line 494 of file CSCEventData.cc.

References CSCALCTHeader::add(), checkALCTClasses(), and theALCTHeader.

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

494  {
496  theALCTHeader->add(digis);
497 }
void checkALCTClasses()
makes new ALCT classes, if needed
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
void add(const std::vector< CSCALCTDigi > &digis)
void CSCEventData::add ( const std::vector< CSCCLCTDigi > &  digis)

Definition at line 499 of file CSCEventData.cc.

References CSCTMBHeader::add(), checkTMBClasses(), theTMBData, and CSCTMBData::tmbHeader().

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

499  {
500  checkTMBClasses();
501  theTMBData->tmbHeader()->add(digis);
502 }
void checkTMBClasses()
makes new TMB classes, if needed
void add(const std::vector< CSCCLCTDigi > &digis)
these methods need more brains to figure which one goes first
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCTMBHeader * tmbHeader()
Definition: CSCTMBData.h:40
void CSCEventData::add ( const std::vector< CSCCorrelatedLCTDigi > &  digis)

Definition at line 504 of file CSCEventData.cc.

References CSCTMBHeader::add(), checkTMBClasses(), theTMBData, and CSCTMBData::tmbHeader().

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

504  {
505  checkTMBClasses();
506  theTMBData->tmbHeader()->add(digis);
507 }
void checkTMBClasses()
makes new TMB classes, if needed
void add(const std::vector< CSCCLCTDigi > &digis)
these methods need more brains to figure which one goes first
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCTMBHeader * tmbHeader()
Definition: CSCTMBData.h:40
void CSCEventData::add ( const std::vector< CSCShowerDigi > &  digis)

Definition at line 509 of file CSCEventData.cc.

References checkTMBClasses().

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

509 { checkTMBClasses(); }
void checkTMBClasses()
makes new TMB classes, if needed
void CSCEventData::add ( const std::vector< GEMPadDigiCluster > &  clusters,
const GEMDetId  
)

Definition at line 511 of file CSCEventData.cc.

References checkTMBClasses().

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

511 { checkTMBClasses(); }
void checkTMBClasses()
makes new TMB classes, if needed
void CSCEventData::addALCTStructures ( )

adds an empty ALCTHeader, trailer, and anode data

void CSCEventData::addComparatorInformation ( std::vector< CSCStripDigi > &  ,
int  layer 
) const
private

adds the comparators to the strip digis

CSCAnodeData * CSCEventData::alctData ( ) const

user must check if nalct > 0

Definition at line 381 of file CSCEventData.cc.

References Exception, nalct(), and theAnodeData.

Referenced by cscdqm::EventProcessor::processCSC().

381  {
382  if (nalct() == 0)
383  throw cms::Exception("No ALCT for this chamber");
384  return theAnodeData;
385 }
int nalct() const
the flag for existence of ALCT data
Definition: CSCEventData.h:67
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
CSCALCTHeader * CSCEventData::alctHeader ( ) const

user must check if nalct > 0

Definition at line 369 of file CSCEventData.cc.

References Exception, nalct(), and theALCTHeader.

Referenced by pack(), and cscdqm::EventProcessor::processCSC().

369  {
370  if (nalct() == 0)
371  throw cms::Exception("No ALCT for this chamber");
372  return theALCTHeader;
373 }
int nalct() const
the flag for existence of ALCT data
Definition: CSCEventData.h:67
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
CSCALCTTrailer * CSCEventData::alctTrailer ( ) const

user must check if nalct > 0

Definition at line 375 of file CSCEventData.cc.

References Exception, nalct(), and theALCTTrailer.

Referenced by pack(), and cscdqm::EventProcessor::processCSC().

375  {
376  if (nalct() == 0)
377  throw cms::Exception("No ALCT for this chamber");
378  return theALCTTrailer;
379 }
int nalct() const
the flag for existence of ALCT data
Definition: CSCEventData.h:67
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
unsigned int CSCEventData::calcALCTcrc ( std::vector< std::pair< unsigned int, unsigned short * > > &  vec)

Definition at line 567 of file CSCEventData.cc.

References mps_fire::i, dqmiolumiharvest::j, dqmiodumpmetadata::n, edm::second(), submitPVValidationJobs::t, and w.

Referenced by pack().

567  {
568  int CRC = 0;
569  // int size=0;
570 
571  for (unsigned int n = 0; n < vec.size(); n++) {
572  // size += vec[n].first;
573  for (uint16_t j = 0, w = 0; j < vec[n].first; j++) {
574  if (vec[n].second != nullptr) {
575  w = vec[n].second[j] & 0xffff;
576  for (uint32_t i = 15, t = 0, ncrc = 0; i < 16; i--) {
577  t = ((w >> i) & 1) ^ ((CRC >> 21) & 1);
578  ncrc = (CRC << 1) & 0x3ffffc;
579  ncrc |= (t ^ (CRC & 1)) << 1;
580  ncrc |= t;
581  CRC = ncrc;
582  }
583  }
584  }
585  }
586 
587  // std::cout << "ALCT CRC vector size: " << size << ", crc: 0x" << std::hex << CRC<< std::endl;
588  return CRC;
589 }
const double w
Definition: UKUtility.cc:23
U second(std::pair< T, U > const &p)
const CSCCFEBData * CSCEventData::cfebData ( unsigned  icfeb) const

unpacked in long mode: has overflow and error bits decoded

Definition at line 367 of file CSCEventData.cc.

References theCFEBData.

Referenced by pack(), and cscdqm::EventProcessor::processCSC().

367 { return theCFEBData[icfeb]; }
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
int CSCEventData::chamberType ( ) const
inline

might not be set in real data

Definition at line 119 of file CSCEventData.h.

References theChamberType.

119 { return theChamberType; }
void CSCEventData::checkALCTClasses ( )
private

makes new ALCT classes, if needed

Definition at line 430 of file CSCEventData.cc.

References CSCDMBHeader::addNALCT(), CSCALCTHeader::alctFirmwareVersion(), cms::cuda::assert(), CSCALCTHeader::setEventInformation(), size(), CSCAnodeData::sizeInWords(), CSCALCTTrailer::sizeInWords(), CSCALCTHeader::sizeInWords(), theALCTHeader, theALCTTrailer, theAnodeData, theChamberType, and theDMBHeader.

Referenced by add().

430  {
431  if (theAnodeData == nullptr) {
432  assert(theChamberType > 0);
438  // set data available flag
440  }
441 }
unsigned short int sizeInWords() const
the amount of the input binary buffer read, in 16-bit words
Definition: CSCAnodeData.h:18
assert(be >=bs)
void setEventInformation(const CSCDMBHeader &)
static int sizeInWords()
in 16-bit frames
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
int sizeInWords()
in 16-bit words
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
unsigned short int alctFirmwareVersion() const
unsigned short size() const
size of the data buffer used, in bytes
Definition: CSCEventData.h:39
void addNALCT()
Definition: CSCDMBHeader.h:26
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
void CSCEventData::checkTMBClasses ( )
private

makes new TMB classes, if needed

Definition at line 443 of file CSCEventData.cc.

References CSCDMBHeader::addNCLCT(), CSCConstants::MAX_CFEBS_RUN1, CSCConstants::MAX_CFEBS_RUN2, CSCTMBHeader::setEventInformation(), CSCTMBHeader::setNCFEBs(), theChamberType, theDMBHeader, theFormatVersion, theTMBData, and CSCTMBData::tmbHeader().

Referenced by add().

443  {
444  int nCFEBs = CSCConstants::MAX_CFEBS_RUN1;
445  if ((theFormatVersion == 2013) && ((theChamberType == 1) || (theChamberType == 2))) {
447  }
448  if (theTMBData == nullptr) {
449  if (theFormatVersion == 2013) { // Set to TMB format for Post-LS1 data
450  theTMBData = new CSCTMBData(2013, 0x7a76, nCFEBs);
451  } else {
452  theTMBData = new CSCTMBData(2007, 0x50c3);
453  }
456  }
457  theTMBData->tmbHeader()->setNCFEBs(nCFEBs);
458  // std::cout << "nCFEBs: " << theTMBData->tmbHeader()->NCFEBs() << std::endl;
459 }
void setEventInformation(const CSCDMBHeader &dmbHeader)
fills fields like bxn and l1a
Definition: CSCTMBHeader.h:38
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
void addNCLCT()
Definition: CSCDMBHeader.h:24
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
void setNCFEBs(uint16_t ncfebs)
Definition: CSCTMBHeader.h:74
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCTMBHeader * tmbHeader()
Definition: CSCTMBData.h:40
CSCComparatorData * CSCEventData::comparatorData ( ) const

user must check if nclct > 0

Definition at line 399 of file CSCEventData.cc.

References CSCTMBData::comparatorData(), Exception, nclct(), and tmbData().

Referenced by cscdqm::EventProcessor::processCSC().

399  {
400  if ((nclct() == 0) || (tmbData() == nullptr))
401  throw cms::Exception("No CLCT data for this chamber");
402  return tmbData()->comparatorData();
403 }
int nclct() const
the number of CLCTs
Definition: CSCEventData.h:70
CSCTMBData * tmbData() const
user must check in nclct &gt; 0
CSCComparatorData * comparatorData()
Definition: CSCTMBData.h:41
void CSCEventData::copy ( const CSCEventData data)
private

Definition at line 283 of file CSCEventData.cc.

References init(), CSCConstants::MAX_CFEBS_RUN2, size_, theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theChamberType, theDMBHeader, theDMBTrailer, theFormatVersion, and theTMBData.

Referenced by SequenceTypes.Task::copyAndAdd(), CSCEventData(), and operator=().

283  {
284  init();
286  theDMBHeader = data.theDMBHeader;
288  if (data.theALCTHeader != nullptr)
290  if (data.theAnodeData != nullptr)
291  theAnodeData = new CSCAnodeData(*(data.theAnodeData));
292  if (data.theALCTTrailer != nullptr)
294  if (data.theTMBData != nullptr)
295  theTMBData = new CSCTMBData(*(data.theTMBData));
296  for (int icfeb = 0; icfeb < CSCConstants::MAX_CFEBS_RUN2; ++icfeb) {
297  theCFEBData[icfeb] = nullptr;
298  if (data.theCFEBData[icfeb] != nullptr)
299  theCFEBData[icfeb] = new CSCCFEBData(*(data.theCFEBData[icfeb]));
300  }
301  size_ = data.size_;
303 }
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
CSCDMBTrailer theDMBTrailer
Definition: CSCEventData.h:162
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
void CSCEventData::destroy ( )
private

Definition at line 305 of file CSCEventData.cc.

References alctZSErecovered, CSCConstants::MAX_CFEBS_RUN2, theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theTMBData, and zseEnable.

Referenced by operator=(), and ~CSCEventData().

305  {
306  if (zseEnable) {
307  delete[] alctZSErecovered;
308  }
309  delete theALCTHeader;
310  delete theAnodeData;
311  delete theALCTTrailer;
312  delete theTMBData;
313  for (int icfeb = 0; icfeb < CSCConstants::MAX_CFEBS_RUN2; ++icfeb) {
314  delete theCFEBData[icfeb];
315  }
316  /*
317  std::cout << "Before delete alctZSErecovered " << std::endl;
318  delete [] alctZSErecovered;
319  std::cout << "After delete alctZSErecovered " << std::endl;
320  */
321 }
unsigned short * alctZSErecovered
Auxiliary bufer to recove the ALCT raw payload from zero suppression.
Definition: CSCEventData.h:170
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
const CSCDMBHeader* CSCEventData::dmbHeader ( ) const
inline

the DAQ motherboard header. A good place for event and chamber info

Definition at line 73 of file CSCEventData.h.

References theDMBHeader.

Referenced by CSCDDUEventData::add(), CSCDigiToRaw::findEventData(), and cscdqm::EventProcessor::processCSC().

73 { return &theDMBHeader; }
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
CSCDMBHeader* CSCEventData::dmbHeader ( )
inline

Definition at line 74 of file CSCEventData.h.

References theDMBHeader.

74 { return &theDMBHeader; }
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
const CSCDMBTrailer* CSCEventData::dmbTrailer ( ) const
inline

DMB trailer.

Definition at line 95 of file CSCEventData.h.

References theDMBTrailer.

Referenced by pack(), and cscdqm::EventProcessor::processCSC().

95 { return &theDMBTrailer; }
CSCDMBTrailer theDMBTrailer
Definition: CSCEventData.h:162
uint16_t CSCEventData::getFormatVersion ( ) const
inline

Definition at line 121 of file CSCEventData.h.

References theFormatVersion.

Referenced by cscdqm::EventProcessor::processCSC().

121 { return theFormatVersion; }
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
void CSCEventData::init ( void  )
private

helpers for ctors, dtor, and op= zeroes all pointers

Definition at line 270 of file CSCEventData.cc.

References alctZSErecovered, CSCConstants::MAX_CFEBS_RUN2, theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theTMBData, and zseEnable.

Referenced by copy(), and unpack_data().

270  {
271  //dataPresent = 0;
272  theALCTHeader = nullptr;
273  theAnodeData = nullptr;
274  theALCTTrailer = nullptr;
275  theTMBData = nullptr;
276  for (int icfeb = 0; icfeb < CSCConstants::MAX_CFEBS_RUN2; ++icfeb) {
277  theCFEBData[icfeb] = nullptr;
278  }
279  alctZSErecovered = nullptr;
280  zseEnable = 0;
281 }
unsigned short * alctZSErecovered
Auxiliary bufer to recove the ALCT raw payload from zero suppression.
Definition: CSCEventData.h:170
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
bool CSCEventData::isALCT ( const uint16_t *  buf)

if dealing with ALCT data

Definition at line 252 of file CSCEventData.cc.

Referenced by unpack_data().

252  {
253  return (((buf[0] & 0xFFFF) == 0xDB0A) || (((buf[0] & 0xF800) == 0x6000) && ((buf[1] & 0xF800) == 0)));
254 }
bool CSCEventData::isTMB ( const uint16_t *  buf)

if dealing with TMB data

Definition at line 256 of file CSCEventData.cc.

Referenced by unpack_data().

256 { return ((buf[0] & 0xFFF) == 0xB0C); }
int CSCEventData::nalct ( ) const
inline

the flag for existence of ALCT data

Definition at line 67 of file CSCEventData.h.

References CSCDMBHeader::nalct(), and theDMBHeader.

Referenced by alctData(), alctHeader(), alctTrailer(), cscdqm::EventProcessor::processCSC(), and unpack_data().

67 { return theDMBHeader.nalct(); }
unsigned nalct() const
Definition: CSCDMBHeader.h:42
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
int CSCEventData::nclct ( ) const
inline

the number of CLCTs

Definition at line 70 of file CSCEventData.h.

References CSCDMBHeader::nclct(), and theDMBHeader.

Referenced by comparatorData(), cscdqm::EventProcessor::processCSC(), tmbData(), tmbHeader(), and unpack_data().

70 { return theDMBHeader.nclct(); }
unsigned nclct() const
Definition: CSCDMBHeader.h:43
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
CSCEventData CSCEventData::operator= ( const CSCEventData data)

Definition at line 262 of file CSCEventData.cc.

References copy(), and destroy().

262  {
263  // check for self-assignment before destructing
264  if (&data != this)
265  destroy();
266  copy(data);
267  return *this;
268 }
void copy(const CSCEventData &)
boost::dynamic_bitset CSCEventData::pack ( )

returns the packed event data.

Definition at line 524 of file CSCEventData.cc.

References alctHeader(), alctTrailer(), bitset_utilities::append(), calcALCTcrc(), cfebData(), CSCAnodeData::data(), CSCDMBHeader::data(), CSCDMBTrailer::data(), CSCALCTTrailer::data(), data, CSCALCTHeader::data(), dmbTrailer(), CSCConstants::MAX_CFEBS_RUN2, CSCTMBData::pack(), CSCALCTHeader::pack(), mps_fire::result, CSCALCTTrailer::setCRC(), CSCAnodeData::sizeInWords(), CSCDMBHeader::sizeInWords(), CSCDMBTrailer::sizeInWords(), CSCALCTTrailer::sizeInWords(), CSCALCTHeader::sizeInWords(), theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theDMBHeader, theDMBTrailer, theTMBData, and bitset_utilities::ushortToBitset().

524  {
525  boost::dynamic_bitset<> result =
527 
528  // Container for CRC calculations
529  std::vector<std::pair<unsigned int, unsigned short*> > crcvec;
530 
531  if (theALCTHeader != nullptr) {
532  boost::dynamic_bitset<> alctHeader = theALCTHeader->pack();
533  result = bitset_utilities::append(result, alctHeader);
534  crcvec.push_back(std::make_pair(theALCTHeader->sizeInWords(), theALCTHeader->data()));
535  }
536  if (theAnodeData != nullptr) {
537  boost::dynamic_bitset<> anodeData =
539  result = bitset_utilities::append(result, anodeData);
540  crcvec.push_back(std::make_pair(theAnodeData->sizeInWords(), theAnodeData->data()));
541  }
542  if (theALCTTrailer != nullptr) {
543  unsigned int crc = calcALCTcrc(crcvec);
544  theALCTTrailer->setCRC(crc);
545  boost::dynamic_bitset<> alctTrailer =
547  result = bitset_utilities::append(result, alctTrailer);
548  }
549  if (theTMBData != nullptr) {
550  result = bitset_utilities::append(result, theTMBData->pack());
551  }
552 
553  for (int icfeb = 0; icfeb < CSCConstants::MAX_CFEBS_RUN2; ++icfeb) {
554  if (theCFEBData[icfeb] != nullptr) {
555  boost::dynamic_bitset<> cfebData =
556  bitset_utilities::ushortToBitset(theCFEBData[icfeb]->sizeInWords() * 16, theCFEBData[icfeb]->data());
557  result = bitset_utilities::append(result, cfebData);
558  }
559  }
560 
561  boost::dynamic_bitset<> dmbTrailer =
563  result = bitset_utilities::append(result, dmbTrailer);
564  return result;
565 }
unsigned short int sizeInWords() const
the amount of the input binary buffer read, in 16-bit words
Definition: CSCAnodeData.h:18
unsigned int calcALCTcrc(std::vector< std::pair< unsigned int, unsigned short * > > &vec)
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
boost::dynamic_bitset pack()
not const because it sets size int TMBTrailer
Definition: CSCTMBData.cc:454
const CSCDMBTrailer * dmbTrailer() const
DMB trailer.
Definition: CSCEventData.h:95
boost::dynamic_bitset append(const boost::dynamic_bitset<> &bs1, const boost::dynamic_bitset<> &bs2)
this method takes two bitsets bs1 and bs2 and returns result of bs2 appended to the end of bs1 ...
unsigned sizeInWords() const
Definition: CSCDMBTrailer.h:61
CSCALCTHeader * alctHeader() const
user must check if nalct &gt; 0
tuple result
Definition: mps_fire.py:311
static int sizeInWords()
in 16-bit frames
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
unsigned short * data()
Definition: CSCDMBTrailer.h:58
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTTrailer * alctTrailer() const
user must check if nalct &gt; 0
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
unsigned sizeInWords() const
Definition: CSCDMBHeader.h:48
boost::dynamic_bitset pack()
CSCDMBTrailer theDMBTrailer
Definition: CSCEventData.h:162
unsigned short * data()
int sizeInWords()
in 16-bit words
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
unsigned short * data()
Definition: CSCDMBHeader.h:53
const CSCCFEBData * cfebData(unsigned icfeb) const
unpacked in long mode: has overflow and error bits decoded
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:79
boost::dynamic_bitset ushortToBitset(const unsigned int numberOfBits, unsigned short *buf)
this method takes numberOfBits bits from unsigned short * array and returns them in the bitset obj...
unsigned short * data()
Definition: CSCAnodeData.h:16
void setCRC(unsigned int crc)
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
unsigned short int * data()
void CSCEventData::selfTest ( )
static

Definition at line 591 of file CSCEventData.cc.

References add(), cms::cuda::assert(), CSCTMBHeader::CLCTDigis(), CSCTMBHeader::CorrelatedLCTDigis(), cscPackAndUnpack(), cscPackerCompare(), CSCDetId::layer(), muon::overlap(), EcalCondDBWriter_cfi::pedestal, DetId::rawId(), stripDigis(), and tmbHeader().

591  {
592  CSCEventData chamberData(5);
593  CSCDetId detId(1, 3, 2, 1, 3);
594  std::vector<CSCCLCTDigi> clctDigis;
595  // Both CLCTs are read-out at the same (pre-trigger) bx, so the last-but-one
596  // arguments in both digis must be the same.
597  clctDigis.push_back(CSCCLCTDigi(1, 1, 4, 1, 0, 30, 3, 2, 1)); // valid for 2007
598  clctDigis.push_back(CSCCLCTDigi(1, 1, 2, 1, 1, 31, 1, 2, 2));
599 
600  // BX of LCT (8th argument) is 1-bit word (the least-significant bit
601  // of ALCT's bx).
602  std::vector<CSCCorrelatedLCTDigi> corrDigis;
603  corrDigis.push_back(CSCCorrelatedLCTDigi(1, 1, 2, 10, 98, 5, 0, 1, 0, 0, 0, 0));
604  corrDigis.push_back(CSCCorrelatedLCTDigi(2, 1, 2, 20, 15, 9, 1, 0, 0, 0, 0, 0));
605 
606  chamberData.add(clctDigis);
607  chamberData.add(corrDigis);
608 
609  CSCWireDigi wireDigi(10, 6);
610  CSCComparatorDigi comparatorDigi(30, 1, 6);
611  chamberData.add(wireDigi, 3);
612  chamberData.add(comparatorDigi, 3);
613 
614  CSCEventData newData = cscPackAndUnpack(chamberData);
615 
616  std::vector<CSCCLCTDigi> clcts = newData.tmbHeader()->CLCTDigis(detId.rawId());
617  assert(cscPackerCompare(clcts[0], clctDigis[0]));
618  assert(cscPackerCompare(clcts[1], clctDigis[1]));
619 
620  std::vector<CSCCorrelatedLCTDigi> lcts = newData.tmbHeader()->CorrelatedLCTDigis(detId.rawId());
621  assert(cscPackerCompare(lcts[0], corrDigis[0]));
622  assert(cscPackerCompare(lcts[1], corrDigis[1]));
623 
624  // test strip digis
625  CSCDetId me1adet1(1, 1, 1, 4, 1);
626  CSCDetId me1bdet1(1, 1, 4, 4, 6);
627  CSCDetId me1adet2(2, 1, 1, 4, 2);
628  CSCDetId me1bdet2(2, 1, 4, 4, 5);
629 
630  std::vector<int> sca(16, 600);
631  std::vector<unsigned short> overflow(16, 0), overlap(16, 0), errorfl(16, 0);
632  CSCStripDigi me1a(5, sca, overflow, overlap, errorfl);
633  CSCStripDigi me1b(8, sca, overflow, overlap, errorfl);
634 
635  CSCEventData forward(1);
636  CSCEventData backward(1);
637 
638  forward.add(me1a, me1adet1.layer());
639  forward.add(me1b, me1bdet1.layer());
640  backward.add(me1a, me1adet2.layer());
641  backward.add(me1b, me1adet2.layer());
642  std::vector<CSCStripDigi> me1afs = forward.stripDigis(me1adet1);
643  std::vector<CSCStripDigi> me1bfs = forward.stripDigis(me1bdet1);
644  std::vector<CSCStripDigi> me1abs = backward.stripDigis(me1adet2);
645  std::vector<CSCStripDigi> me1bbs = backward.stripDigis(me1bdet2);
646  //FIXME The current code works under the assumption that ME11 and ME1A
647  // go into separate EventData. They need to be combined.
648  assert(me1afs.size() == 16);
649  assert(me1bfs.size() == 16);
650  assert(me1abs.size() == 16);
651  assert(me1bbs.size() == 16);
652 
653  assert(me1afs[4].getStrip() == 5);
654  assert(me1bfs[7].getStrip() == 8);
655  assert(me1abs[4].getStrip() == 5);
656  assert(me1bbs[7].getStrip() == 8);
657  assert(me1afs[4].pedestal() == 600);
658  assert(me1bfs[7].pedestal() == 600);
659  assert(me1abs[4].pedestal() == 600);
660  assert(me1bbs[7].pedestal() == 600);
661 }
std::vector< CSCCorrelatedLCTDigi > CorrelatedLCTDigis(uint32_t idlayer) const
returns CorrelatedLCT digis
Definition: CSCTMBHeader.h:97
assert(be >=bs)
bool overlap(const reco::Muon &muon1, const reco::Muon &muon2, double pullX=1.0, double pullY=1.0, bool checkAdjacentChambers=false)
std::vector< CSCCLCTDigi > CLCTDigis(uint32_t idlayer)
returns CLCT digis
Definition: CSCTMBHeader.h:94
bool cscPackerCompare(const T &t1, const T &t2)
CSCTMBHeader * tmbHeader() const
user must check if nclct &gt; 0
T cscPackAndUnpack(T &t)
static void CSCEventData::setDebug ( const bool  value)
inlinestatic

turns on/off debug flag for this class

Definition at line 42 of file CSCEventData.h.

References debug, and relativeConstraints::value.

Referenced by CSCDCCUnpacker::CSCDCCUnpacker().

42 { debug = value; }
static std::atomic< bool > debug
Definition: CSCEventData.h:128
void CSCEventData::setEventInformation ( int  bxnum,
int  lvl1num 
)

this will fill the DMB header, and change all related fields in the DMBTrailer, ALCTHeader, and TMBHeader

Definition at line 405 of file CSCEventData.cc.

References cms::cuda::assert(), CSCConstants::MAX_CFEBS_RUN1, CSCConstants::MAX_CFEBS_RUN2, CSCDMBHeader::setBXN(), CSCALCTHeader::setEventInformation(), CSCTMBHeader::setEventInformation(), CSCDMBHeader::setL1A(), CSCCFEBData::setL1A(), CSCDMBHeader::setL1A24(), CSCTMBHeader::setNCFEBs(), theALCTHeader, theCFEBData, theChamberType, theDMBHeader, theFormatVersion, theTMBData, and CSCTMBData::tmbHeader().

Referenced by CSCDDUEventData::add().

405  {
406  theDMBHeader.setBXN(bxnum);
407  theDMBHeader.setL1A(lvl1num);
408  theDMBHeader.setL1A24(lvl1num);
409  if (theALCTHeader) {
411  }
412  if (theTMBData) {
414 
415  assert(theChamberType > 0);
416 
418 
419  // Set number of CFEBs to 7 for Post-LS1 ME11
420  if ((theFormatVersion == 2013) && ((theChamberType == 1) || (theChamberType == 2))) {
422  }
423  }
424  for (unsigned cfeb = 0; cfeb < CSCConstants::MAX_CFEBS_RUN2; cfeb++) {
425  if (theCFEBData[cfeb])
426  theCFEBData[cfeb]->setL1A(lvl1num);
427  }
428 }
void setEventInformation(const CSCDMBHeader &dmbHeader)
fills fields like bxn and l1a
Definition: CSCTMBHeader.h:38
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
void setL1A(unsigned l1a)
Definition: CSCCFEBData.cc:159
void setL1A(int l1a)
Definition: CSCDMBHeader.h:28
assert(be >=bs)
void setEventInformation(const CSCDMBHeader &)
void setBXN(int bxn)
Definition: CSCDMBHeader.h:27
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
void setNCFEBs(uint16_t ncfebs)
Definition: CSCTMBHeader.h:74
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
void setL1A24(int l1a)
Definition: CSCDMBHeader.h:29
CSCTMBHeader * tmbHeader()
Definition: CSCTMBData.h:40
unsigned short CSCEventData::size ( void  ) const
inline

size of the data buffer used, in bytes

Definition at line 39 of file CSCEventData.h.

References size_.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), and checkALCTClasses().

39 { return size_; }
std::vector< CSCStripDigi > CSCEventData::stripDigis ( const CSCDetId idlayer) const

returns all the strip digis in the chamber, with the comparator information.

Definition at line 323 of file CSCEventData.cc.

References CSCConstants::MAX_CFEBS_RUN2, mps_fire::result, and stripDigis().

Referenced by operator<<(), and selfTest().

323  {
324  std::vector<CSCStripDigi> result;
325  for (unsigned icfeb = 0; icfeb < CSCConstants::MAX_CFEBS_RUN2; ++icfeb) {
326  std::vector<CSCStripDigi> newDigis = stripDigis(idlayer, icfeb);
327  result.insert(result.end(), newDigis.begin(), newDigis.end());
328  }
329  return result;
330 }
std::vector< std::vector< CSCStripDigi > > stripDigis() const
deprecated. Use the above methods instead
tuple result
Definition: mps_fire.py:311
std::vector< CSCStripDigi > CSCEventData::stripDigis ( unsigned  idlayer,
unsigned  icfeb 
) const

returns all the strip digis in the chamber's cfeb

Definition at line 332 of file CSCEventData.cc.

References CSCCFEBData::digis(), mps_fire::result, and theCFEBData.

332  {
333  std::vector<CSCStripDigi> result;
334  if (theCFEBData[icfeb] != nullptr) {
335  std::vector<CSCStripDigi> newDigis = theCFEBData[icfeb]->digis(idlayer);
336  result.insert(result.end(), newDigis.begin(), newDigis.end());
337  }
338 
339  return result;
340 }
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
tuple result
Definition: mps_fire.py:311
void digis(uint32_t idlayer, std::vector< CSCStripDigi > &result) const
faster way to get to digis
Definition: CSCCFEBData.cc:202
std::vector< std::vector< CSCStripDigi > > CSCEventData::stripDigis ( ) const

deprecated. Use the above methods instead

Definition at line 350 of file CSCEventData.cc.

References phase1PixelTopology::layer, CSCDetId::maxLayerId(), CSCDetId::minLayerId(), and mps_fire::result.

Referenced by stripDigis().

350  {
351  std::vector<std::vector<CSCStripDigi> > result;
353  std::vector<CSCStripDigi> digis = stripDigis(layer);
354  result.push_back(digis);
355  }
356  return result;
357 }
std::vector< std::vector< CSCStripDigi > > stripDigis() const
deprecated. Use the above methods instead
constexpr std::array< uint8_t, layerIndexSize > layer
tuple result
Definition: mps_fire.py:311
static int minLayerId()
Definition: CSCDetId.h:242
static int maxLayerId()
Definition: CSCDetId.h:243
CSCTMBData * CSCEventData::tmbData ( ) const

user must check in nclct > 0

Definition at line 387 of file CSCEventData.cc.

References Exception, nclct(), and theTMBData.

Referenced by comparatorData(), cscdqm::EventProcessor::processCSC(), and tmbHeader().

387  {
388  if (nclct() == 0)
389  throw cms::Exception("No CLCT for this chamber");
390  return theTMBData;
391 }
int nclct() const
the number of CLCTs
Definition: CSCEventData.h:70
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCTMBHeader * CSCEventData::tmbHeader ( ) const

user must check if nclct > 0

Definition at line 393 of file CSCEventData.cc.

References Exception, nclct(), tmbData(), and CSCTMBData::tmbHeader().

Referenced by selfTest().

393  {
394  if ((nclct() == 0) || (tmbData() == nullptr))
395  throw cms::Exception("No CLCT header for this chamber");
396  return tmbData()->tmbHeader();
397 }
int nclct() const
the number of CLCTs
Definition: CSCEventData.h:70
CSCTMBData * tmbData() const
user must check in nclct &gt; 0
CSCTMBHeader * tmbHeader()
Definition: CSCTMBData.h:40
void CSCEventData::unpack_data ( const uint16_t *  buf)
private

The size of the ALCT payload is determined here

Check if Zero Suppression ALCT Enabled

Aauxilary variables neede to recover zero suppression Calculate the number of wire groups per layer

Calculate the number of words in the layer

unsigned short * posZSEtmpALCT = pos; This is just to dump the actual ALCT payload ** begin ** For debuggin purposes

to_rm (8 lines)

This is just to dump the actual ALCT payload ** end **

Actual word counting and recovering the original ALCT payload

Convert the recovered vector into the array

This is to check the content of the recovered ALCT payload to_rm (7 lines)

Definition at line 39 of file CSCEventData.cc.

References CSCALCTHeader::ALCTDigis(), alctZSErecovered, visDQMUpload::buf, CSCDMBTrailer::cfeb_endtimeout(), CSCDMBTrailer::cfeb_starttimeout(), CSCDMBHeader::cfebAvailable(), CSCDMBHeader::check(), CSCALCTHeader::check(), CSCALCTHeader::data(), debug, CSCDMBHeader::format_version(), mps_fire::i, init(), isALCT(), isTMB(), dqmiolumiharvest::j, cmsLHEtoEOSManager::l, LogTrace, CSCConstants::MAX_CFEBS_RUN2, nalct(), nclct(), CSCTMBData::size(), size_, CSCAnodeData::sizeInWords(), CSCCFEBData::sizeInWords(), CSCDMBHeader::sizeInWords(), CSCDMBTrailer::sizeInWords(), CSCALCTTrailer::sizeInWords(), CSCALCTHeader::sizeInWords(), theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theDMBHeader, theDMBTrailer, theFormatVersion, theTMBData, and zseEnable.

Referenced by CSCEventData().

39  {
40  // zero everything
41  init();
42  const uint16_t* pos = buf;
43  if (debug) {
44  LogTrace("CSCEventData|CSCRawToDigi") << "The event data ";
45  for (int i = 0; i < 16; ++i) {
46  LogTrace("CSCEventData|CSCRawToDigi") << std::hex << pos[i] << " ";
47  }
48  }
49 
51  if (!(theDMBHeader.check())) {
52  LogTrace("CSCEventData|CSCRawToDigi") << "Bad DMB Header??? "
53  << " first four words: ";
54  for (int i = 0; i < 4; ++i) {
55  LogTrace("CSCEventData|CSCRawToDigi") << std::hex << pos[i] << " ";
56  }
57  }
58 
59  if (debug) {
60  LogTrace("CSCEventData|CSCRawToDigi") << "nalct = " << nalct();
61  LogTrace("CSCEventData|CSCRawToDigi") << "nclct = " << nclct();
62  }
63 
64  if (debug) {
65  LogTrace("CSCEventData|CSCRawToDigi") << "size in words of DMBHeader" << theDMBHeader.sizeInWords();
66  LogTrace("CSCEventData|CSCRawToDigi") << "sizeof(DMBHeader)" << sizeof(theDMBHeader);
67  }
68 
69  pos += theDMBHeader.sizeInWords();
70 
71  if (nalct() == 1) {
72  if (isALCT(pos)) //checking for ALCTData
73  {
74  theALCTHeader = new CSCALCTHeader(pos);
75  if (!theALCTHeader->check()) {
76  LogTrace("CSCEventData|CSCRawToDigi") << "+++WARNING: Corrupt ALCT data - won't attempt to decode";
77  } else {
78  //dataPresent|=0x40;
79  pos += theALCTHeader->sizeInWords(); //size of the header
80  //fill ALCT Digis
82 
83  //theAnodeData = new CSCAnodeData(*theALCTHeader, pos);
84 
86  /*
87  std::cout << " ****The ALCT information from CSCEventData.cc (begin)**** " << std::endl; ///to_rm
88  std::cout << " alctHeader2007().size: " << theALCTHeader->alctHeader2007().sizeInWords() << std::endl; ///to_rm
89  std::cout << " ALCT Header Content: " << std::endl; ///to_rm
91  for(int k=0; k<theALCTHeader->sizeInWords(); k+=4){
92  std::cout << std::hex << theALCTHeader->data()[k+3]
93  << " " << theALCTHeader->data()[k+2]
94  << " " << theALCTHeader->data()[k+1]
95  << " " << theALCTHeader->data()[k] << std::dec << std::endl;
96  }
97  */
98  //std::cout << " ALCT Size: " << theAnodeData->sizeInWords() << std::endl;
100  // int zseEnable = 0;
101  zseEnable = (theALCTHeader->data()[5] & 0x1000) >> 12;
102  //std::cout << " ZSE Bit: " << zseEnable << std::endl; /// to_rm
103  int sizeInWord_ZSE = 0;
104 
105  //alctZSErecovered = new unsigned short [theAnodeData->sizeInWords()];
106 
107  if (zseEnable) {
110  int nWGs_per_layer = ((theALCTHeader->data()[6] & 0x0007) + 1) * 16;
112  int nWG_round_up = int(nWGs_per_layer / 12) + (nWGs_per_layer % 3 ? 1 : 0);
113  //std::cout << " Words per layer: " << nWG_round_up << std::endl; ///to_rm
114  const uint16_t* posZSE = pos;
115  std::vector<unsigned short> alctZSErecoveredVector;
116  alctZSErecoveredVector.clear();
117 
118  //alctZSErecovered = new unsigned short [theAnodeData->sizeInWords()];
119  //delete [] alctZSErecovered;
120  //std::cout << " ALCT Buffer with ZSE: " << std::endl; ///to_rm
124  //unsigned short * posZSEdebug = pos; ///to_rm
125 
127  /*
128  while (*posZSEdebug != 0xDE0D){
129  unsigned short d = *posZSEdebug;
130  unsigned short c = *(posZSEdebug+1);
131  unsigned short b = *(posZSEdebug+2);
132  unsigned short a = *(posZSEdebug+3);
133  posZSEdebug+=4;
134  std::cout << std::hex << a << " " << b << " " << c << " " << d << std::dec << std::endl;
135  }
136  */
138 
140  int alctZSErecoveredPos = 0;
141  while (*posZSE != 0xDE0D) {
142  if ((*posZSE == 0x1000) && (*posZSE != 0x3000)) {
143  for (int j = 0; j < nWG_round_up; j++) {
144  alctZSErecoveredVector.push_back(0x0000);
145  }
146  alctZSErecoveredPos += nWG_round_up;
147  } else {
148  alctZSErecoveredVector.push_back(*posZSE);
149  ++alctZSErecoveredPos;
150  }
151  posZSE++;
152  sizeInWord_ZSE++;
153  }
154 
155  alctZSErecovered = new unsigned short[alctZSErecoveredVector.size()];
156 
158  for (int l = 0; l < (int)alctZSErecoveredVector.size(); l++) {
159  alctZSErecovered[l] = alctZSErecoveredVector[l];
160  }
161 
162  unsigned short* posRecovered = alctZSErecovered;
163  theAnodeData = new CSCAnodeData(*theALCTHeader, posRecovered);
164 
167  /*
168  std::cout << " The ALCT payload recovered: " << std::endl;
169  for(int k=0; k<theAnodeData->sizeInWords(); k+=4){
170  std::cout << std::hex << alctZSErecovered[k+3] << " "
171  << alctZSErecovered[k+2] << " "
172  << alctZSErecovered[k+1] << " "
173  << alctZSErecovered[k] << std::dec << std::endl;
174  }
175  */
176  //delete [] alctZSErecovered;
177  //std::cout << " ALCT SizeZSE : " << sizeInWord_ZSE << std::endl; ///to_rm
178  //std::cout << " ALCT SizeZSE Recovered: " << alctZSErecoveredPos << std::endl; ///to_rm
179  //std::cout << " ALCT Size Expected: " << theAnodeData->sizeInWords() << std::endl; ///to_rm
180  pos += sizeInWord_ZSE;
181  } else {
182  //pos +=sizeInWord_ZSE;
184  pos += theAnodeData->sizeInWords(); // size of the data is determined during unpacking
185  }
186  //std::cout << " ****The ALCT information from CSCEventData.cc (end)**** " << std::endl; ///to_rm
187  theALCTTrailer = new CSCALCTTrailer(pos);
188  pos += theALCTTrailer->sizeInWords();
189  }
190  } else {
191  LogTrace("CSCEventData|CSCRawToDigi") << "Error:nalct reported but no ALCT data found!!!";
192  }
193  }
194 
195  if (nclct() == 1) {
196  if (isTMB(pos)) {
197  //dataPresent|=0x20;
198  theTMBData = new CSCTMBData(pos); //fill all TMB data
199  pos += theTMBData->size();
200  } else {
201  LogTrace("CSCEventData|CSCRawToDigi") << "Error:nclct reported but no TMB data found!!!";
202  }
203  }
204 
205  //now let's try to find and unpack the DMBTrailer
206  bool dmbTrailerReached = false;
207  for (int i = 0; i < 12000; ++i) //8000 max for cfeb + 1980ALCT + 287 TMB
208  {
209  dmbTrailerReached = (*(i + pos) & 0xF000) == 0xF000 && (*(i + pos + 1) & 0xF000) == 0xF000 &&
210  (*(i + pos + 2) & 0xF000) == 0xF000 && (*(i + pos + 3) & 0xF000) == 0xF000 &&
211  (*(i + pos + 4) & 0xF000) == 0xE000 && (*(i + pos + 5) & 0xF000) == 0xE000 &&
212  (*(i + pos + 6) & 0xF000) == 0xE000 && (*(i + pos + 7) & 0xF000) == 0xE000;
213  if (dmbTrailerReached) {
214  // theDMBTrailer = *( (CSCDMBTrailer *) (pos+i) );
216  break;
217  }
218  }
219  if (dmbTrailerReached) {
220  for (int icfeb = 0; icfeb < CSCConstants::MAX_CFEBS_RUN2; ++icfeb) {
221  theCFEBData[icfeb] = nullptr;
222  int cfeb_available = theDMBHeader.cfebAvailable(icfeb);
223  unsigned int cfebTimeout = theDMBTrailer.cfeb_starttimeout() | theDMBTrailer.cfeb_endtimeout();
224  //cfeb_available cannot be trusted - need additional verification!
225  if (cfeb_available == 1) {
226  if ((cfebTimeout >> icfeb) & 1) {
227  if (debug)
228  LogTrace("CSCEventData|CSCRawToDigi") << "CFEB Timed out! ";
229  } else {
230  //dataPresent|=(0x1>>icfeb);
231  // Fill CFEB data and convert it into cathode digis
232 
233  // Check if we have here DCFEB using DMB format version field (new ME11 with DCFEBs - 0x2, other chamber types 0x1)
234  bool isDCFEB = false;
235  if (theDMBHeader.format_version() == 2)
236  isDCFEB = true;
237 
238  theCFEBData[icfeb] = new CSCCFEBData(icfeb, pos, theFormatVersion, isDCFEB);
239  pos += theCFEBData[icfeb]->sizeInWords();
240  }
241  }
242  }
243  pos += theDMBTrailer.sizeInWords();
244  size_ = pos - buf;
245  } else {
246  LogTrace("CSCEventData|CSCRawToDigi") << "Critical Error: DMB Trailer was not found!!! ";
247  }
248 
249  // std::cout << "CSC format: " << theFormatVersion << " " << getFormatVersion() << std::endl;
250 }
bool check() const
Definition: CSCDMBHeader.h:51
unsigned short int sizeInWords() const
the amount of the input binary buffer read, in 16-bit words
Definition: CSCAnodeData.h:18
unsigned short * alctZSErecovered
Auxiliary bufer to recove the ALCT raw payload from zero suppression.
Definition: CSCEventData.h:170
int nclct() const
the number of CLCTs
Definition: CSCEventData.h:70
uint16_t theFormatVersion
Output Format Version (2005, 2013)
Definition: CSCEventData.h:174
CSCCFEBData * theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
for up to MAX_CFEB CFEB boards
Definition: CSCEventData.h:160
unsigned cfeb_endtimeout() const
Definition: CSCDMBTrailer.h:33
const unsigned short size() const
Definition: CSCTMBData.h:38
unsigned sizeInWords() const
Definition: CSCDMBTrailer.h:61
bool check() const
#define LogTrace(id)
int nalct() const
the flag for existence of ALCT data
Definition: CSCEventData.h:67
static int sizeInWords()
in 16-bit frames
bool isALCT(const uint16_t *buf)
if dealing with ALCT data
CSCDMBHeader theDMBHeader
Definition: CSCEventData.h:149
CSCTMBData * theTMBData
Definition: CSCEventData.h:157
CSCALCTHeader * theALCTHeader
Definition: CSCEventData.h:154
unsigned cfeb_starttimeout() const
Definition: CSCDMBTrailer.h:37
unsigned sizeInWords() const
Definition: CSCDMBHeader.h:48
CSCDMBTrailer theDMBTrailer
Definition: CSCEventData.h:162
int sizeInWords()
in 16-bit words
std::vector< CSCALCTDigi > ALCTDigis() const
CSCALCTTrailer * theALCTTrailer
Definition: CSCEventData.h:156
unsigned format_version() const
Definition: CSCDMBHeader.h:49
bool isTMB(const uint16_t *buf)
if dealing with TMB data
unsigned sizeInWords() const
Definition: CSCCFEBData.h:46
static std::atomic< bool > debug
Definition: CSCEventData.h:128
bool cfebAvailable(unsigned icfeb)
Definition: CSCDMBHeader.h:21
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
unsigned short int * data()
std::vector< CSCWireDigi > CSCEventData::wireDigis ( unsigned  ilayer) const

Definition at line 342 of file CSCEventData.cc.

References theAnodeData, and CSCAnodeData::wireDigis().

Referenced by operator<<().

342  {
343  if (theAnodeData == nullptr) {
344  return std::vector<CSCWireDigi>();
345  } else {
346  return theAnodeData->wireDigis(ilayer);
347  }
348 }
std::vector< CSCWireDigi > wireDigis(int layer) const
input layer is from 1 to 6
Definition: CSCAnodeData.h:21
CSCAnodeData * theAnodeData
Definition: CSCEventData.h:155
std::vector< std::vector< CSCWireDigi > > CSCEventData::wireDigis ( ) const

deprecated. Use the above method instead.

Definition at line 359 of file CSCEventData.cc.

References phase1PixelTopology::layer, CSCDetId::maxLayerId(), CSCDetId::minLayerId(), and mps_fire::result.

359  {
360  std::vector<std::vector<CSCWireDigi> > result;
362  result.push_back(wireDigis(layer));
363  }
364  return result;
365 }
std::vector< std::vector< CSCWireDigi > > wireDigis() const
deprecated. Use the above method instead.
constexpr std::array< uint8_t, layerIndexSize > layer
tuple result
Definition: mps_fire.py:311
static int minLayerId()
Definition: CSCDetId.h:242
static int maxLayerId()
Definition: CSCDetId.h:243

Member Data Documentation

unsigned short* CSCEventData::alctZSErecovered
private

Auxiliary bufer to recove the ALCT raw payload from zero suppression.

Definition at line 170 of file CSCEventData.h.

Referenced by destroy(), init(), and unpack_data().

std::atomic< bool > CSCEventData::debug {false}
static
int CSCEventData::size_
private

Definition at line 164 of file CSCEventData.h.

Referenced by copy(), size(), and unpack_data().

CSCALCTHeader* CSCEventData::theALCTHeader
private

these may or may not be present. I decided to make them dynamic because most CSC chambers don't have LCTs, therefore don't have data, except for DMB headers and trailers.

Definition at line 154 of file CSCEventData.h.

Referenced by add(), alctHeader(), checkALCTClasses(), copy(), destroy(), init(), pack(), setEventInformation(), and unpack_data().

CSCALCTTrailer* CSCEventData::theALCTTrailer
private

Definition at line 156 of file CSCEventData.h.

Referenced by alctTrailer(), checkALCTClasses(), copy(), destroy(), init(), pack(), and unpack_data().

CSCAnodeData* CSCEventData::theAnodeData
private
CSCCFEBData* CSCEventData::theCFEBData[CSCConstants::MAX_CFEBS_RUN2]
private

for up to MAX_CFEB CFEB boards

Definition at line 160 of file CSCEventData.h.

Referenced by add(), cfebData(), copy(), CSCEventData(), destroy(), init(), pack(), setEventInformation(), stripDigis(), and unpack_data().

int CSCEventData::theChamberType
private

this won't be filled when real data is read it. It's only used when packing simulated data, so we know how many wire and strip channels to make.

Definition at line 167 of file CSCEventData.h.

Referenced by chamberType(), checkALCTClasses(), checkTMBClasses(), copy(), and setEventInformation().

CSCDMBHeader CSCEventData::theDMBHeader
private
CSCDMBTrailer CSCEventData::theDMBTrailer
private

Definition at line 162 of file CSCEventData.h.

Referenced by copy(), dmbTrailer(), pack(), and unpack_data().

uint16_t CSCEventData::theFormatVersion
private

Output Format Version (2005, 2013)

Definition at line 174 of file CSCEventData.h.

Referenced by add(), checkTMBClasses(), copy(), CSCEventData(), getFormatVersion(), setEventInformation(), and unpack_data().

CSCTMBData* CSCEventData::theTMBData
private
int CSCEventData::zseEnable
private

Definition at line 171 of file CSCEventData.h.

Referenced by destroy(), init(), and unpack_data().