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	addALCTStructures ()
	adds an empty ALCTHeader, trailer, and anode data More...
CSCAnodeData *	alctData () const
	user must check if nalct > 0 More...
CSCALCTHeader *	alctHeader () const
	user must check if nalct > 0 More...
CSCALCTTrailer *	alctTrailer () const
	user must check if nalct > 0 More...
unsigned int	calcALCTcrc (std::vector< std::pair< unsigned int, unsigned short * > > &vec)
CSCCFEBData *	cfebData (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...
CSCCLCTData *	clctData () const
	user must check if nclct > 0 More...
	CSCEventData (int chamberType, uint16_t format_version=2005)
	CSCEventData (unsigned short *buf, uint16_t format_version=2005)
	should make const input soon More...
	CSCEventData ()
	CSCEventData (const CSCEventData &data)
const CSCDMBHeader *	dmbHeader () const
	the DAQ motherboard header. A good place for event and chamber info More...
CSCDMBHeader *	dmbHeader ()
const CSCDMBTrailer *	dmbTrailer () const
	DMB trailer. More...
uint16_t	getFormatVersion () const
bool	isALCT (const short unsigned int *buf)
	if dealing with ALCT data More...
bool	isTMB (const short unsigned int *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< CSCStripDigi >	stripDigis (const CSCDetId &idlayer) const
	returns all the strip digis in the chamber, with the comparator information. More...
std::vector< CSCStripDigi >	stripDigis (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...
CSCTMBData *	tmbData () const
	user must check in nclct > 0 More...
CSCTMBHeader *	tmbHeader () const
	user must check if nclct > 0 More...
std::vector< CSCWireDigi >	wireDigis (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 (unsigned short *buf)

Private Attributes
unsigned short *	alctZSErecovered
	Auxiliary bufer to recove the ALCT raw payload from zero suppression. More...
int	size_
CSCALCTHeader *	theALCTHeader
CSCALCTTrailer *	theALCTTrailer
CSCAnodeData *	theAnodeData
CSCCFEBData *	theCFEBData [7]
	for up to MAX_CFEB CFEB boards More...
int	theChamberType
CSCDMBHeader	theDMBHeader
CSCDMBTrailer	theDMBTrailer
uint16_t	theFormatVersion
	Output Format Version (2005, 2013) More...
CSCTMBData *	theTMBData
int	zseEnable

Detailed Description

Definition at line 42 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 i, MAX_CFEB, and theCFEBData.

                                                                   :
    theDMBHeader(format_version),
    theALCTHeader(0),
    theAnodeData(0),
    theALCTTrailer(0),
    theTMBData(0),
    theDMBTrailer(format_version),
    theChamberType(chamberType),
    alctZSErecovered(0),
    zseEnable(0),
    theFormatVersion(format_version)
{

  for (unsigned i = 0; i < MAX_CFEB; ++i)
    {
      theCFEBData[i] = 0;
    }
}

CSCEventData::CSCEventData	(	unsigned short *	buf,
		uint16_t	format_version = 2005
	)

should make const input soon

Definition at line 38 of file CSCEventData.cc.

References theFormatVersion, and unpack_data().

                                                                       : theFormatVersion(format_version)
{
  theFormatVersion = format_version;
  unpack_data(buf);
}

CSCEventData::CSCEventData ( )

inline

Definition at line 47 of file CSCEventData.h.

{}

CSCEventData::~CSCEventData

(

)

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

Definition at line 310 of file CSCEventData.cc.

References destroy().

{
  destroy();
}

CSCEventData::CSCEventData

(

const CSCEventData &

data

)

Definition at line 305 of file CSCEventData.cc.

References copy().

{
  copy(data);
}

Member Function Documentation

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

routines to add digis to the data

Definition at line 569 of file CSCEventData.cc.

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

Referenced by CSCDigiToRaw::add(), and selfTest().

{
  //@@ need special logic here for ME11
  unsigned cfeb = (digi.getStrip()-1)/16;
  bool sixteenSamples = false;
  if (digi.getADCCounts().size()==16) sixteenSamples = true;
  if (theCFEBData[cfeb] == 0)
    {
      bool isDCFEB = false;
      if (theDMBHeader.format_version() == 2) isDCFEB = true;
      theCFEBData[cfeb] = new CSCCFEBData(cfeb, sixteenSamples, theFormatVersion, isDCFEB);
      theDMBHeader.addCFEB(cfeb);
    }
  theCFEBData[cfeb]->add(digi, layer);
  
  
}

void CSCEventData::add	(	const CSCWireDigi &	digi,
		int	layer
	)

Definition at line 588 of file CSCEventData.cc.

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

{
  checkALCTClasses();
  theAnodeData->add(digi, layer);
  theALCTHeader->setDAVForChannel(digi.getWireGroup());
  theALCTHeader->setBXNCount(digi.getWireGroupBX());
  
}

void CSCEventData::add	(	const CSCComparatorDigi &	digi,
		int	layer
	)

Definition at line 597 of file CSCEventData.cc.

References CSCCLCTData::add(), checkTMBClasses(), CSCTMBData::clctData(), and theTMBData.

{
  checkTMBClasses();
  theTMBData->clctData()->add(digi, layer);
  
}

void CSCEventData::add	(	const CSCComparatorDigi &	digi,
		const CSCDetId &	cid
	)

Definition at line 604 of file CSCEventData.cc.

References CSCCLCTData::add(), checkTMBClasses(), CSCTMBData::clctData(), and theTMBData.

{
  checkTMBClasses();
  theTMBData->clctData()->add(digi, cid);
  
}

void CSCEventData::add

(

const std::vector< CSCALCTDigi > &

digis

)

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

Definition at line 613 of file CSCEventData.cc.

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

{
  checkALCTClasses();
  theALCTHeader->add(digis);
}

void CSCEventData::add

(

const std::vector< CSCCLCTDigi > &

digis

)

Definition at line 620 of file CSCEventData.cc.

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

{
  checkTMBClasses();
  theTMBData->tmbHeader()->add(digis);
}

void CSCEventData::add

(

const std::vector< CSCCorrelatedLCTDigi > &

digis

)

Definition at line 626 of file CSCEventData.cc.

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

{
  checkTMBClasses();
  theTMBData->tmbHeader()->add(digis);
}

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 470 of file CSCEventData.cc.

References edm::hlt::Exception, nalct(), and theAnodeData.

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

{
  if (nalct() == 0) throw cms::Exception("No ALCT for this chamber");
  return theAnodeData;
}

CSCALCTHeader * CSCEventData::alctHeader

(

)

const

user must check if nalct > 0

Definition at line 457 of file CSCEventData.cc.

References edm::hlt::Exception, nalct(), and theALCTHeader.

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

{
  if (nalct() == 0) throw cms::Exception("No ALCT for this chamber");
  return theALCTHeader;
}

CSCALCTTrailer * CSCEventData::alctTrailer

(

)

const

user must check if nalct > 0

Definition at line 463 of file CSCEventData.cc.

References edm::hlt::Exception, nalct(), and theALCTTrailer.

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

{
  if (nalct() == 0) throw cms::Exception("No ALCT for this chamber");
  return theALCTTrailer;
}

unsigned int CSCEventData::calcALCTcrc

(

std::vector< std::pair< unsigned int, unsigned short * > > &

vec

)

Definition at line 700 of file CSCEventData.cc.

References i, j, n, NULL, edm::second(), edmStreamStallGrapher::t, and w().

Referenced by pack().

{
  int CRC=0;
//  int size=0;

  for (unsigned int n = 0; n < vec.size(); n++)
    {
//      size += vec[n].first;
      for (uint16_t j=0, w=0; j<vec[n].first; j++ )
        {
     
          if (vec[n].second != NULL) { 
          w = vec[n].second[j] & 0xffff;
          for (uint32_t i=15, t=0, ncrc=0; i<16; i--)
            {
              t = ((w >> i) & 1) ^ ((CRC >> 21) & 1);
              ncrc = (CRC << 1) & 0x3ffffc;
              ncrc |= (t ^ (CRC & 1)) << 1;
              ncrc |= t;
              CRC = ncrc;
            }

      }

        }
    }

  //  std::cout << "ALCT CRC vector size: " << size <<  ", crc: 0x" << std::hex << CRC<< std::endl;
  return CRC;
}

CSCCFEBData * CSCEventData::cfebData

(

unsigned

icfeb

)

const

unpacked in long mode: has overflow and error bits decoded

Definition at line 451 of file CSCEventData.cc.

References theCFEBData.

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

{
  return theCFEBData[icfeb];
}

int CSCEventData::chamberType ( ) const

inline

might not be set in real data

Definition at line 141 of file CSCEventData.h.

References theChamberType.

{return theChamberType;}

void CSCEventData::checkALCTClasses ( )

private

makes new ALCT classes, if needed

Definition at line 532 of file CSCEventData.cc.

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

Referenced by add().

{
  if (theAnodeData == NULL)
    {
      assert(theChamberType>0);
      theALCTHeader = new CSCALCTHeader(theChamberType);
      theALCTHeader->setEventInformation(theDMBHeader);
      theAnodeData = new CSCAnodeData(*theALCTHeader);
      int size = theALCTHeader->sizeInWords() + theAnodeData->sizeInWords() + CSCALCTTrailer::sizeInWords();
      theALCTTrailer = new CSCALCTTrailer(size, theALCTHeader->alctFirmwareVersion());
      // set data available flag
      theDMBHeader.addNALCT();
    }
}

void CSCEventData::checkTMBClasses ( )

private

makes new TMB classes, if needed

Definition at line 548 of file CSCEventData.cc.

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

Referenced by add().

{
  int nCFEBs = 5;
  if ((theFormatVersion  == 2013) && ((theChamberType == 1) || (theChamberType == 2)) ) {
      nCFEBs = 7;
  }
  if (theTMBData == NULL)
    {
      if (theFormatVersion  == 2013) { // Set to TMB format for Post-LS1 data
        theTMBData = new CSCTMBData(2013, 0x7a76, nCFEBs);
      } else {
    theTMBData = new CSCTMBData(2007, 0x50c3);
      }
      theTMBData->tmbHeader()->setEventInformation(theDMBHeader);
      theDMBHeader.addNCLCT();
    }
      theTMBData->tmbHeader()->setNCFEBs(nCFEBs);
   // std::cout << "nCFEBs: " << theTMBData->tmbHeader()->NCFEBs() << std::endl;
}

CSCCLCTData * CSCEventData::clctData

(

)

const

user must check if nclct > 0

Definition at line 489 of file CSCEventData.cc.

References CSCTMBData::clctData(), edm::hlt::Exception, nclct(), NULL, and tmbData().

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

{
  if ((nclct() == 0)||(tmbData()==NULL)) throw cms::Exception("No CLCT data for this chamber");
  return tmbData()->clctData();
}

void CSCEventData::copy ( const CSCEventData &  data )

private

Definition at line 341 of file CSCEventData.cc.

References init(), MAX_CFEB, NULL, size_, theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theChamberType, theDMBHeader, theDMBTrailer, theFormatVersion, and theTMBData.

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

{
  init();
  theFormatVersion = data.theFormatVersion;
  theDMBHeader  = data.theDMBHeader;
  theDMBTrailer = data.theDMBTrailer;
  if (data.theALCTHeader != NULL)
    theALCTHeader  = new CSCALCTHeader(*(data.theALCTHeader));
  if (data.theAnodeData != NULL)
    theAnodeData   = new CSCAnodeData(*(data.theAnodeData));
  if (data.theALCTTrailer != NULL)
    theALCTTrailer = new CSCALCTTrailer(*(data.theALCTTrailer));
  if (data.theTMBData != NULL)
    theTMBData     = new CSCTMBData(*(data.theTMBData));
  for (int icfeb = 0; icfeb < MAX_CFEB; ++icfeb)
    {
      theCFEBData[icfeb] = 0;
      if (data.theCFEBData[icfeb] != NULL)
        theCFEBData[icfeb] = new CSCCFEBData(*(data.theCFEBData[icfeb]));
    }
  size_  = data.size_;
  theChamberType = data.theChamberType;


}

void CSCEventData::destroy ( )

private

Definition at line 368 of file CSCEventData.cc.

References alctZSErecovered, MAX_CFEB, theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theTMBData, and zseEnable.

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

{
  if (zseEnable)
    {
      delete [] alctZSErecovered;
    }
  delete theALCTHeader;
  delete theAnodeData;
  delete theALCTTrailer;
  delete theTMBData;
  for (int icfeb = 0; icfeb < MAX_CFEB; ++icfeb)
    {
      delete theCFEBData[icfeb];
    }
  /*
    std::cout << "Before delete alctZSErecovered " << std::endl;
    delete [] alctZSErecovered;
    std::cout << "After delete alctZSErecovered " << std::endl;
  */
}

const CSCDMBHeader* CSCEventData::dmbHeader ( ) const

inline

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

Definition at line 96 of file CSCEventData.h.

References theDMBHeader.

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

{return &theDMBHeader;}

CSCDMBHeader* CSCEventData::dmbHeader ( )

inline

Definition at line 97 of file CSCEventData.h.

References theDMBHeader.

{return &theDMBHeader;}

const CSCDMBTrailer* CSCEventData::dmbTrailer ( ) const

inline

DMB trailer.

Definition at line 118 of file CSCEventData.h.

References theDMBTrailer.

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

{return &theDMBTrailer;}

uint16_t CSCEventData::getFormatVersion ( ) const

inline

Definition at line 143 of file CSCEventData.h.

References theFormatVersion.

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

{ return theFormatVersion; }

void CSCEventData::init ( void  )

private

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

Definition at line 325 of file CSCEventData.cc.

References alctZSErecovered, MAX_CFEB, theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theTMBData, and zseEnable.

Referenced by copy(), and unpack_data().

{
  //dataPresent = 0;
  theALCTHeader = 0;
  theAnodeData = 0;
  theALCTTrailer = 0;
  theTMBData = 0;
  for (int icfeb = 0; icfeb < MAX_CFEB; ++icfeb)
    {
      theCFEBData[icfeb] = 0;
    }
  alctZSErecovered=0;
  zseEnable=0;
}

bool CSCEventData::isALCT

(

const short unsigned int *

buf

)

if dealing with ALCT data

Definition at line 293 of file CSCEventData.cc.

Referenced by unpack_data().

{
  return (((buf[0]&0xFFFF)==0xDB0A)||(((buf[0]&0xF800)==0x6000)&&((buf[1]&0xF800)==0)));
}

bool CSCEventData::isTMB

(

const short unsigned int *

buf

)

if dealing with TMB data

Definition at line 298 of file CSCEventData.cc.

Referenced by unpack_data().

{
  return ((buf[0]&0xFFF)==0xB0C);
}

int CSCEventData::nalct ( ) const

inline

the flag for existence of ALCT data

Definition at line 90 of file CSCEventData.h.

References CSCDMBHeader::nalct(), and theDMBHeader.

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

{return theDMBHeader.nalct();}

int CSCEventData::nclct ( ) const

inline

the number of CLCTs

Definition at line 93 of file CSCEventData.h.

References CSCDMBHeader::nclct(), and theDMBHeader.

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

{return theDMBHeader.nclct();}

CSCEventData CSCEventData::operator=

(

const CSCEventData &

data

)

Definition at line 316 of file CSCEventData.cc.

References copy(), and destroy().

{
  // check for self-assignment before destructing
  if (&data != this) destroy();
  copy(data);
  return *this;
}

boost::dynamic_bitset CSCEventData::pack

(

)

returns the packed event data.

Definition at line 648 of file CSCEventData.cc.

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

Referenced by pyrootRender.interactiveRender::draw().

{
  boost::dynamic_bitset<> result = bitset_utilities::ushortToBitset( theDMBHeader.sizeInWords()*16,
                                   theDMBHeader.data());

  // Container for CRC calculations
  std::vector<std::pair<unsigned int, unsigned short*> > crcvec;

  if (theALCTHeader != NULL)
    {
      boost::dynamic_bitset<> alctHeader = theALCTHeader->pack();
      result = bitset_utilities::append(result, alctHeader);
      crcvec.push_back(std::make_pair(theALCTHeader->sizeInWords(), theALCTHeader->data()));
    }
  if (theAnodeData != NULL)
    {
      boost::dynamic_bitset<> anodeData = bitset_utilities::ushortToBitset (theAnodeData->sizeInWords()*16,
                                          theAnodeData->data());
      result = bitset_utilities::append(result, anodeData);
      crcvec.push_back(std::make_pair(theAnodeData->sizeInWords(), theAnodeData->data()));

    }
  if (theALCTTrailer != NULL)
    {
      unsigned int crc = calcALCTcrc(crcvec);
      theALCTTrailer->setCRC(crc);
      boost::dynamic_bitset<> alctTrailer =bitset_utilities::ushortToBitset(theALCTTrailer->sizeInWords()*16,
                                           theALCTTrailer->data());
      result = bitset_utilities::append(result, alctTrailer);

    }
  if (theTMBData != NULL)
    {
      result  = bitset_utilities::append(result, theTMBData->pack());
    }

  for (int icfeb = 0;  icfeb < MAX_CFEB;  ++icfeb)
    {
      if (theCFEBData[icfeb] != NULL)
        {
          boost::dynamic_bitset<> cfebData = bitset_utilities::ushortToBitset(theCFEBData[icfeb]->sizeInWords()*16,
                                             theCFEBData[icfeb]->data());
          result = bitset_utilities::append(result, cfebData);
        }
    }

  boost::dynamic_bitset<> dmbTrailer = bitset_utilities::ushortToBitset( theDMBTrailer.sizeInWords()*16,
                                       theDMBTrailer.data());
  result = bitset_utilities::append(result, dmbTrailer);
  return result;
}

void CSCEventData::selfTest ( )

static

Definition at line 733 of file CSCEventData.cc.

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

{
  CSCEventData chamberData(5);
  CSCDetId detId(1, 3, 2, 1, 3);
  std::vector<CSCCLCTDigi> clctDigis;
  // Both CLCTs are read-out at the same (pre-trigger) bx, so the last-but-one
  // arguments in both digis must be the same.
  clctDigis.push_back(CSCCLCTDigi(1, 1, 4, 1, 0, 30, 3, 2, 1)); // valid for 2007
  clctDigis.push_back(CSCCLCTDigi(1, 1, 2, 1, 1, 31, 1, 2, 2));

  // BX of LCT (8th argument) is 1-bit word (the least-significant bit
  // of ALCT's bx).
  std::vector<CSCCorrelatedLCTDigi> corrDigis;
  corrDigis.push_back(CSCCorrelatedLCTDigi(1, 1, 2, 10, 98, 5, 0, 1, 0, 0, 0, 0));
  corrDigis.push_back(CSCCorrelatedLCTDigi(2, 1, 2, 20, 15, 9, 1, 0, 0, 0, 0, 0));

  chamberData.add(clctDigis);
  chamberData.add(corrDigis);

  CSCWireDigi wireDigi(10, 6);
  CSCComparatorDigi comparatorDigi(30, 1, 6);
  chamberData.add(wireDigi, 3);
  chamberData.add(comparatorDigi, 3);

  CSCEventData newData = cscPackAndUnpack(chamberData);

  std::vector<CSCCLCTDigi> clcts = newData.tmbHeader()->CLCTDigis(detId.rawId());
  assert(cscPackerCompare(clcts[0],clctDigis[0]));
  assert(cscPackerCompare(clcts[1],clctDigis[1]));

  std::vector<CSCCorrelatedLCTDigi> lcts = newData.tmbHeader()->CorrelatedLCTDigis(detId.rawId());
  assert(cscPackerCompare(lcts[0], corrDigis[0]));
  assert(cscPackerCompare(lcts[1], corrDigis[1]));

  // test strip digis
  CSCDetId me1adet1(1, 1, 1, 4, 1);
  CSCDetId me1bdet1(1, 1, 4, 4, 6);
  CSCDetId me1adet2(2, 1, 1, 4, 2);
  CSCDetId me1bdet2(2, 1, 4, 4, 5);

  std::vector<int> sca(16, 600);
  std::vector<unsigned short> overflow(16, 0), overlap(16, 0), errorfl(16,0);
  CSCStripDigi me1a(5, sca, overflow, overlap, errorfl);
  CSCStripDigi me1b(8, sca, overflow, overlap, errorfl);

  CSCEventData forward(1);
  CSCEventData backward(1);

  forward.add(me1a, me1adet1.layer());
  forward.add(me1b, me1bdet1.layer());
  backward.add(me1a, me1adet2.layer());
  backward.add(me1b, me1adet2.layer());
  std::vector<CSCStripDigi> me1afs = forward.stripDigis(me1adet1);
  std::vector<CSCStripDigi> me1bfs = forward.stripDigis(me1bdet1);
  std::vector<CSCStripDigi> me1abs = backward.stripDigis(me1adet2);
  std::vector<CSCStripDigi> me1bbs = backward.stripDigis(me1bdet2);
  //FIXME The current code works under the assumption that ME11 and ME1A
  // go into separate EventData.  They need to be combined.
  assert(me1afs.size() == 16);
  assert(me1bfs.size() == 16);
  assert(me1abs.size() == 16);
  assert(me1bbs.size() == 16);

  assert(me1afs[4].getStrip() == 5);
  assert(me1bfs[7].getStrip() == 8);
  assert(me1abs[4].getStrip() == 5);
  assert(me1bbs[7].getStrip() == 8);
  assert(me1afs[4].pedestal() == 600);
  assert(me1bfs[7].pedestal() == 600);
  assert(me1abs[4].pedestal() == 600);
  assert(me1bbs[7].pedestal() == 600);


}

static void CSCEventData::setDebug ( const bool  value )

inlinestatic

turns on/off debug flag for this class

Definition at line 58 of file CSCEventData.h.

References debug, and relativeConstraints::value.

Referenced by CSCDCCUnpacker::CSCDCCUnpacker().

{debug = value;}

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 496 of file CSCEventData.cc.

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

Referenced by CSCDDUEventData::add().

{
  theDMBHeader.setBXN(bxnum);
  theDMBHeader.setL1A(lvl1num);
  theDMBHeader.setL1A24(lvl1num);
  if (theALCTHeader)
    {
      theALCTHeader->setEventInformation(theDMBHeader);
    }
  if (theTMBData)
    {
      theTMBData->tmbHeader()->setEventInformation(theDMBHeader);

      assert(theChamberType>0);

      theTMBData->tmbHeader()->setNCFEBs(5);

      // Set number of CFEBs to 7 for Post-LS1 ME11
      if ((theFormatVersion  == 2013) && ((theChamberType == 1) || (theChamberType == 2)) ) {
        theTMBData->tmbHeader()->setNCFEBs(7);
      }
/*
      // Set number of CFEBs to 4 for ME13 chambers
      if (theChamberType == 4) { 
    theTMBData->tmbHeader()->setNCFEBs(4);
      }
*/
    
    }
  for (unsigned cfeb=0; cfeb< 7; cfeb++)
    {
      if (theCFEBData[cfeb]) theCFEBData[cfeb]->setL1A(lvl1num); 
    }
}

unsigned short CSCEventData::size ( void  ) const

inline

size of the data buffer used, in bytes

Definition at line 55 of file CSCEventData.h.

References size_.

Referenced by checkALCTClasses().

{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 390 of file CSCEventData.cc.

References MAX_CFEB, query::result, and stripDigis().

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

{
  std::vector<CSCStripDigi> result;
  for (unsigned icfeb = 0; icfeb < MAX_CFEB; ++icfeb)
    {
      std::vector<CSCStripDigi> newDigis = stripDigis(idlayer, icfeb);
      result.insert(result.end(), newDigis.begin(), newDigis.end());
    }
  return result;
}

std::vector< CSCStripDigi > CSCEventData::stripDigis	(	unsigned	idlayer,
		unsigned	icfeb
	)		const

returns all the strip digis in the chamber's cfeb

Definition at line 402 of file CSCEventData.cc.

References CSCCFEBData::digis(), NULL, query::result, and theCFEBData.

{
  //  assert(ilayer > 0 && ilayer <= 6); // off because now idlayer is raw cscdetid
  std::vector<CSCStripDigi> result;
  if (theCFEBData[icfeb] != NULL)
    {
      std::vector<CSCStripDigi> newDigis = theCFEBData[icfeb]->digis(idlayer);
      result.insert(result.end(), newDigis.begin(), newDigis.end());
    }

  return result;
}

std::vector< std::vector< CSCStripDigi > > CSCEventData::stripDigis

(

)

const

deprecated. Use the above methods instead

Definition at line 429 of file CSCEventData.cc.

References query::result.

Referenced by stripDigis().

{
  std::vector < std::vector<CSCStripDigi> > result;
  for (int layer = 1; layer <= 6; ++layer)
    {
      std::vector<CSCStripDigi> digis = stripDigis(layer);
      result.push_back(digis);
    }
  return result;
}

CSCTMBData * CSCEventData::tmbData

(

)

const

user must check in nclct > 0

Definition at line 476 of file CSCEventData.cc.

References edm::hlt::Exception, nclct(), and theTMBData.

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

{
  if (nclct() == 0) throw cms::Exception("No CLCT for this chamber");
  return theTMBData;
}

CSCTMBHeader * CSCEventData::tmbHeader

(

)

const

user must check if nclct > 0

Definition at line 483 of file CSCEventData.cc.

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

Referenced by selfTest().

{
  if ((nclct() == 0)||(tmbData()==NULL)) throw cms::Exception("No CLCT header for this chamber");
  return tmbData()->tmbHeader();
}

void CSCEventData::unpack_data ( unsigned short *  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 45 of file CSCEventData.cc.

References CSCALCTHeader::ALCTDigis(), alctZSErecovered, CSCDMBTrailer::cfeb_endtimeout(), CSCDMBTrailer::cfeb_starttimeout(), CSCDMBHeader::cfebAvailable(), CSCDMBHeader::check(), CSCALCTHeader::check(), CSCALCTHeader::data(), debug, CSCDMBHeader::format_version(), i, init(), isALCT(), isTMB(), j, prof2calltree::l, LogTrace, MAX_CFEB, 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().

{
  // zero everything
  init();
  unsigned short * pos = buf;
  if (debug)
    {
      LogTrace ("CSCEventData|CSCRawToDigi") << "The event data ";
      for (int i = 0; i < 16; ++i)
        {
          LogTrace ("CSCEventData|CSCRawToDigi") << std::hex << pos[i ] << " ";
        }
    }

  theDMBHeader = CSCDMBHeader(pos, theFormatVersion);
  if (!(theDMBHeader.check()))
    {
      LogTrace ("CSCEventData|CSCRawToDigi")  << "Bad DMB Header??? " << " first four words: ";
      for (int i = 0; i < 4; ++i)
        {
          LogTrace ("CSCEventData|CSCRawToDigi") << std::hex << pos[i ] << " ";
        }
    }


  if (debug)
    {
      LogTrace ("CSCEventData|CSCRawToDigi") << "nalct = " << nalct();
      LogTrace ("CSCEventData|CSCRawToDigi") << "nclct = " << nclct();
    }

  if (debug)
    {
      LogTrace ("CSCEventData|CSCRawToDigi") << "size in words of DMBHeader" << theDMBHeader.sizeInWords();
      LogTrace ("CSCEventData|CSCRawToDigi") << "sizeof(DMBHeader)" << sizeof(theDMBHeader);
    }

  pos += theDMBHeader.sizeInWords();

  if (nalct() ==1)
    {
      if (isALCT(pos))  //checking for ALCTData
        {
          theALCTHeader = new CSCALCTHeader( pos );
          if (!theALCTHeader->check())
            {
              LogTrace ("CSCEventData|CSCRawToDigi") <<"+++WARNING: Corrupt ALCT data - won't attempt to decode";
            }
          else
            {
              //dataPresent|=0x40;
              pos += theALCTHeader->sizeInWords(); //size of the header
              //fill ALCT Digis
              theALCTHeader->ALCTDigis();

              //theAnodeData = new CSCAnodeData(*theALCTHeader, pos);


              /*
              std::cout << " ****The ALCT information from CSCEventData.cc (begin)**** " << std::endl; ///to_rm
              std::cout << " alctHeader2007().size: " << theALCTHeader->alctHeader2007().sizeInWords() << std::endl; ///to_rm
              std::cout << " ALCT Header Content: " << std::endl; ///to_rm
              for(int k=0; k<theALCTHeader->sizeInWords(); k+=4){
                 std::cout << std::hex << theALCTHeader->data()[k+3]
                           << " " << theALCTHeader->data()[k+2]
                           << " " << theALCTHeader->data()[k+1]
                           << " " << theALCTHeader->data()[k] << std::dec << std::endl;
                 }
               */
              //std::cout << " ALCT Size: " << theAnodeData->sizeInWords() << std::endl;
              // int zseEnable = 0;
              zseEnable = (theALCTHeader->data()[5] & 0x1000) >> 12;
              //std::cout << " ZSE Bit: " <<  zseEnable << std::endl; /// to_rm
              int sizeInWord_ZSE =0;

              //alctZSErecovered = new unsigned short [theAnodeData->sizeInWords()];

              if (zseEnable)
                {
                  int  nWGs_per_layer = ( (theALCTHeader->data()[6]&0x0007) + 1 ) * 16 ;
                  int nWG_round_up   = int(nWGs_per_layer/12)+(nWGs_per_layer%3?1:0);
                  //std::cout << " Words per layer: " << nWG_round_up << std::endl; ///to_rm
                  unsigned short * posZSE = pos;
                  std::vector<unsigned short> alctZSErecoveredVector;
                  alctZSErecoveredVector.clear();

                  //alctZSErecovered = new unsigned short [theAnodeData->sizeInWords()];
                  //delete [] alctZSErecovered;
                  //std::cout << " ALCT Buffer with ZSE: " << std::endl; ///to_rm
                  //unsigned short * posZSEdebug = pos; ///to_rm

                  /*
                  while (*posZSEdebug != 0xDE0D){
                        unsigned short d = *posZSEdebug;
                        unsigned short c = *(posZSEdebug+1);
                        unsigned short b = *(posZSEdebug+2);
                        unsigned short a = *(posZSEdebug+3);
                        posZSEdebug+=4;
                        std::cout << std::hex << a << " " << b << " " << c << " " << d << std::dec << std::endl;
                  }
                  */

                  int alctZSErecoveredPos=0;
                  while (*posZSE != 0xDE0D)
                    {
                      if ( (*posZSE == 0x1000) && (*posZSE != 0x3000))
                        {
                          for (int j=0; j<nWG_round_up; j++)
                            {
                              alctZSErecoveredVector.push_back(0x0000);
                            }
                          alctZSErecoveredPos+=nWG_round_up;
                        }
                      else
                        {
                          alctZSErecoveredVector.push_back(*posZSE);
                          ++alctZSErecoveredPos;
                        }
                      posZSE++;
                      sizeInWord_ZSE++;
                    }

                  alctZSErecovered = new unsigned short [alctZSErecoveredVector.size()];

                  for (int l=0; l<(int)alctZSErecoveredVector.size(); l++)
                    {
                      alctZSErecovered[l]=alctZSErecoveredVector[l];
                    }

                  unsigned short *posRecovered = alctZSErecovered;
                  theAnodeData = new CSCAnodeData(*theALCTHeader, posRecovered);

                  /*
                  std::cout << " The ALCT payload recovered: " << std::endl;
                  for(int k=0; k<theAnodeData->sizeInWords(); k+=4){
                     std::cout << std::hex << alctZSErecovered[k+3] << " "
                                           << alctZSErecovered[k+2] << " "
                                           << alctZSErecovered[k+1] << " "
                                           << alctZSErecovered[k] << std::dec << std::endl;
                  }
                  */
                  //delete [] alctZSErecovered;
                  //std::cout << " ALCT SizeZSE : " << sizeInWord_ZSE << std::endl; ///to_rm
                  //std::cout << " ALCT SizeZSE Recovered: " << alctZSErecoveredPos << std::endl; ///to_rm
                  //std::cout << " ALCT Size Expected: " << theAnodeData->sizeInWords() << std::endl; ///to_rm
                  pos +=sizeInWord_ZSE;
                }
              else
                {
                  //pos +=sizeInWord_ZSE;
                  theAnodeData = new CSCAnodeData(*theALCTHeader, pos);
                  pos += theAnodeData->sizeInWords(); // size of the data is determined during unpacking
                }
              //std::cout << " ****The ALCT information from CSCEventData.cc (end)**** " << std::endl; ///to_rm
              theALCTTrailer = new CSCALCTTrailer( pos );
              pos += theALCTTrailer->sizeInWords();
            }
        }
      else
        {
          LogTrace ("CSCEventData|CSCRawToDigi") << "Error:nalct reported but no ALCT data found!!!";
        }
    }

  if (nclct() ==1)
    {
      if (isTMB(pos))
        {
          //dataPresent|=0x20;
          theTMBData = new CSCTMBData(pos);  //fill all TMB data
          pos += theTMBData->size();
        }
      else
        {
          LogTrace ("CSCEventData|CSCRawToDigi") << "Error:nclct reported but no TMB data found!!!";
        }
    }

  //now let's try to find and unpack the DMBTrailer
  bool dmbTrailerReached= false;
  for (int i=0; i<12000; ++i)  //8000 max for cfeb + 1980ALCT + 287 TMB
    {
      dmbTrailerReached =
        (*(i+pos) & 0xF000) == 0xF000 && (*(i+pos+1) & 0xF000) == 0xF000
        && (*(i+pos+2) & 0xF000) == 0xF000 && (*(i+pos+3) & 0xF000) == 0xF000
        && (*(i+pos+4) & 0xF000) == 0xE000 && (*(i+pos+5) & 0xF000) == 0xE000
        && (*(i+pos+6) & 0xF000) == 0xE000 && (*(i+pos+7) & 0xF000) == 0xE000;
      if (dmbTrailerReached)
        {
          // theDMBTrailer = *( (CSCDMBTrailer *) (pos+i) );
          theDMBTrailer = CSCDMBTrailer(pos+i, theFormatVersion);
          break;
        }
    }
  if (dmbTrailerReached)
    {
      for (int icfeb = 0; icfeb < MAX_CFEB; ++icfeb)
        {
          theCFEBData[icfeb] = 0;
          int cfeb_available = theDMBHeader.cfebAvailable(icfeb);
          unsigned int cfebTimeout = theDMBTrailer.cfeb_starttimeout() | theDMBTrailer.cfeb_endtimeout();
          //cfeb_available cannot be trusted - need additional verification!
          if ( cfeb_available==1 )
            {
              if ((cfebTimeout >> icfeb) & 1)
                {
                  if (debug) LogTrace ("CSCEventData|CSCRawToDigi") << "CFEB Timed out! ";
                }
              else
                {
                  //dataPresent|=(0x1>>icfeb);
                  // Fill CFEB data and convert it into cathode digis

                  // Check if we have here DCFEB  using DMB format version field (new ME11 with DCFEBs - 0x2, other chamber types 0x1)
                  bool isDCFEB = false;
                  if (theDMBHeader.format_version() == 2) isDCFEB = true;

                  theCFEBData[icfeb] = new CSCCFEBData(icfeb, pos, theFormatVersion, isDCFEB);
                  pos += theCFEBData[icfeb]->sizeInWords();
                }
            }
        }
      pos += theDMBTrailer.sizeInWords();
      size_ = pos-buf;
    }
  else
    {
      LogTrace ("CSCEventData|CSCRawToDigi") << "Critical Error: DMB Trailer was not found!!! ";
    }

  // std::cout << "CSC format: " << theFormatVersion << " " << getFormatVersion() << std::endl;
}

std::vector< CSCWireDigi > CSCEventData::wireDigis

(

unsigned

ilayer

)

const

Definition at line 416 of file CSCEventData.cc.

References theAnodeData, and CSCAnodeData::wireDigis().

Referenced by operator<<().

{
  if (theAnodeData == 0)
    {
      return std::vector<CSCWireDigi>();
    }
  else
    {
      return theAnodeData->wireDigis(ilayer);
    }
}

std::vector< std::vector< CSCWireDigi > > CSCEventData::wireDigis

(

)

const

deprecated. Use the above method instead.

Definition at line 440 of file CSCEventData.cc.

References query::result.

{
  std::vector < std::vector<CSCWireDigi> > result;
  for (int layer = 1; layer <= 6; ++layer)
    {
      result.push_back(wireDigis(layer));
    }
  return result;
}

Member Data Documentation

unsigned short* CSCEventData::alctZSErecovered

private

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

Definition at line 193 of file CSCEventData.h.

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

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

static

Definition at line 151 of file CSCEventData.h.

Referenced by rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), setDebug(), unpack_data(), and pkg.AbstractPkg::write().

int CSCEventData::size_

private

Definition at line 187 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 177 of file CSCEventData.h.

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

CSCALCTTrailer* CSCEventData::theALCTTrailer

private

Definition at line 179 of file CSCEventData.h.

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

CSCAnodeData* CSCEventData::theAnodeData

private

Definition at line 178 of file CSCEventData.h.

Referenced by add(), alctData(), checkALCTClasses(), copy(), destroy(), init(), pack(), unpack_data(), and wireDigis().

CSCCFEBData* CSCEventData::theCFEBData[7]

private

for up to MAX_CFEB CFEB boards

Definition at line 183 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 190 of file CSCEventData.h.

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

CSCDMBHeader CSCEventData::theDMBHeader

private

Definition at line 172 of file CSCEventData.h.

Referenced by add(), checkALCTClasses(), checkTMBClasses(), copy(), dmbHeader(), nalct(), nclct(), pack(), setEventInformation(), and unpack_data().

CSCDMBTrailer CSCEventData::theDMBTrailer

private

Definition at line 185 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 197 of file CSCEventData.h.

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

CSCTMBData* CSCEventData::theTMBData

private

Definition at line 180 of file CSCEventData.h.

Referenced by add(), checkTMBClasses(), copy(), destroy(), init(), pack(), setEventInformation(), tmbData(), and unpack_data().

int CSCEventData::zseEnable

private

Definition at line 194 of file CSCEventData.h.

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