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 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 39 of file CSCEventData.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 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 34 of file CSCEventData.cc.

References theFormatVersion, and unpack_data().

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

CSCEventData::CSCEventData ( )

inline

Definition at line 44 of file CSCEventData.h.

{}

CSCEventData::~CSCEventData

(

)

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

Definition at line 306 of file CSCEventData.cc.

References destroy().

{
  destroy();
}

CSCEventData::CSCEventData

(

const CSCEventData &

data

)

Definition at line 301 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 565 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 584 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 593 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 600 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 609 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 616 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 622 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 466 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 453 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 459 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 696 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 447 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 138 of file CSCEventData.h.

References theChamberType.

{return theChamberType;}

void CSCEventData::checkALCTClasses ( )

private

makes new ALCT classes, if needed

Definition at line 528 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 544 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 485 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 337 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 364 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 93 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 94 of file CSCEventData.h.

References theDMBHeader.

{return &theDMBHeader;}

const CSCDMBTrailer* CSCEventData::dmbTrailer ( ) const

inline

DMB trailer.

Definition at line 115 of file CSCEventData.h.

References theDMBTrailer.

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

{return &theDMBTrailer;}

uint16_t CSCEventData::getFormatVersion ( ) const

inline

Definition at line 140 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 321 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 289 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 294 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 87 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 90 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 312 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 644 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 729 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 55 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 492 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 52 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 386 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 398 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 425 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 472 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 479 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 41 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 412 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 436 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 186 of file CSCEventData.h.

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

bool CSCEventData::debug = false

static

Definition at line 145 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 180 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 170 of file CSCEventData.h.

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

CSCALCTTrailer* CSCEventData::theALCTTrailer

private

Definition at line 172 of file CSCEventData.h.

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

CSCAnodeData* CSCEventData::theAnodeData

private

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

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

CSCDMBHeader CSCEventData::theDMBHeader

private

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

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

CSCTMBData* CSCEventData::theTMBData

private

Definition at line 173 of file CSCEventData.h.

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

int CSCEventData::zseEnable

private

Definition at line 187 of file CSCEventData.h.

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