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 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...
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)
	CSCEventData (unsigned short *buf)
	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...
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 [5]
	for up to 5 CFEB boards More...
int	theChamberType
CSCDMBHeader	theDMBHeader
CSCDMBTrailer	theDMBTrailer
CSCTMBData *	theTMBData
int	zseEnable

Detailed Description

Definition at line 36 of file CSCEventData.h.

Constructor & Destructor Documentation

CSCEventData::CSCEventData ( int  chamberType )

explicit

Definition at line 14 of file CSCEventData.cc.

References i, and theCFEBData.

                                          : 
  theDMBHeader(), 
  theALCTHeader(0), 
  theAnodeData(0),
  theALCTTrailer(0),
  theTMBData(0),
  theDMBTrailer(),
  theChamberType(chamberType),
  alctZSErecovered(0),
  zseEnable(0)
{
  
  for(unsigned i = 0; i < 5; ++i) {
    theCFEBData[i] = 0;
  }
}

CSCEventData::CSCEventData

(

unsigned short *

buf

)

should make const input soon

Definition at line 32 of file CSCEventData.cc.

References unpack_data().

                                              {
  unpack_data(buf);
}

CSCEventData::CSCEventData ( )

inline

Definition at line 41 of file CSCEventData.h.

{}

CSCEventData::~CSCEventData

(

)

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

Definition at line 261 of file CSCEventData.cc.

References destroy().

                            {
  destroy();
}

CSCEventData::CSCEventData

(

const CSCEventData &

data

)

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

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

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)    {
    theCFEBData[cfeb] = new CSCCFEBData(cfeb, sixteenSamples);
    theDMBHeader.addCFEB(cfeb);
  }
  theCFEBData[cfeb]->add(digi, layer);
}

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

Definition at line 468 of file CSCEventData.cc.

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

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

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

Definition at line 474 of file CSCEventData.cc.

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

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

void CSCEventData::add

(

const std::vector< CSCALCTDigi > &

digis

)

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

Definition at line 481 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 487 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 492 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 396 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 385 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 390 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;
}

CSCCFEBData * CSCEventData::cfebData

(

unsigned

icfeb

)

const

unpacked in long mode: has overflow and error bits decoded

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

References theChamberType.

{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(), 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 445 of file CSCEventData.cc.

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

Referenced by add().

{
  if(theTMBData == NULL)    {
    theTMBData = new CSCTMBData();
    theTMBData->tmbHeader()->setEventInformation(theDMBHeader);
    theDMBHeader.addNCLCT();
  }
}

CSCCLCTData * CSCEventData::clctData

(

)

const

user must check if nclct > 0

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

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

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

                                                 {
  init();
  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 < 5; ++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 311 of file CSCEventData.cc.

References alctZSErecovered, 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 < 5; ++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 90 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 91 of file CSCEventData.h.

References theDMBHeader.

{return &theDMBHeader;}

const CSCDMBTrailer* CSCEventData::dmbTrailer ( ) const

inline

DMB trailer.

Definition at line 112 of file CSCEventData.h.

References theDMBTrailer.

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

{return &theDMBTrailer;}

void CSCEventData::init ( void  )

private

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

Definition at line 274 of file CSCEventData.cc.

References alctZSErecovered, 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 < 5; ++icfeb) {
    theCFEBData[icfeb] = 0;
  }
  alctZSErecovered=0;
  zseEnable=0;
}

bool CSCEventData::isALCT

(

const short unsigned int *

buf

)

if dealing with ALCT data

Definition at line 247 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 251 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 84 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 87 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 266 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 511 of file CSCEventData.cc.

References alctHeader(), alctTrailer(), bitset_utilities::append(), cfebData(), CSCAnodeData::data(), CSCDMBTrailer::data(), CSCDMBHeader::data(), CSCALCTTrailer::data(), data, dmbTrailer(), NULL, CSCTMBData::pack(), CSCALCTHeader::pack(), query::result, CSCAnodeData::sizeInWords(), CSCDMBHeader::sizeInWords(), CSCALCTTrailer::sizeInWords(), CSCDMBTrailer::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());

  if(theALCTHeader != NULL)     {
    boost::dynamic_bitset<> alctHeader = theALCTHeader->pack();
    result = bitset_utilities::append(result, alctHeader);
  }
  if(theAnodeData != NULL) {
    boost::dynamic_bitset<> anodeData = bitset_utilities::ushortToBitset (theAnodeData->sizeInWords()*16,
                                      theAnodeData->data());
    result = bitset_utilities::append(result, anodeData);
  }
  if(theALCTTrailer != NULL)  {
    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 < 5;  ++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 548 of file CSCEventData.cc.

References add(), CSCTMBHeader::CLCTDigis(), CSCTMBHeader::CorrelatedLCTDigis(), cscPackAndUnpack(), cscPackerCompare(), CSCDetId::layer(), muon::overlap(), 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 52 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 418 of file CSCEventData.cc.

References CSCDMBHeader::setBXN(), CSCTMBHeader::setEventInformation(), CSCALCTHeader::setEventInformation(), CSCDMBHeader::setL1A(), theALCTHeader, theDMBHeader, theTMBData, and CSCTMBData::tmbHeader().

Referenced by CSCDDUEventData::add().

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

unsigned short CSCEventData::size ( void  ) const

inline

size of the data buffer used, in bytes

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

References query::result, and stripDigis().

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

                                                                               {
  std::vector<CSCStripDigi> result;
  for(unsigned icfeb = 0; icfeb < 5; ++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 340 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 362 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 401 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 407 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 37 of file CSCEventData.cc.

References CSCALCTHeader::ALCTDigis(), alctZSErecovered, CSCDMBTrailer::cfeb_endtimeout, CSCDMBTrailer::cfeb_starttimeout, CSCDMBHeader::cfebAvailable(), CSCDMBHeader::check(), CSCALCTHeader::check(), CSCALCTHeader::data(), debug, i, init(), isALCT(), isTMB(), j, ConfigFiles::l, LogTrace, nalct(), nclct(), CSCTMBData::size(), size_, CSCAnodeData::sizeInWords(), CSCDMBHeader::sizeInWords(), CSCCFEBData::sizeInWords(), CSCALCTTrailer::sizeInWords(), CSCDMBTrailer::sizeInWords(), CSCALCTHeader::sizeInWords(), theALCTHeader, theALCTTrailer, theAnodeData, theCFEBData, theDMBHeader, theDMBTrailer, 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);
  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) );
      break;
    }
  }
  if (dmbTrailerReached) {
    for(int icfeb = 0; icfeb < 5; ++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
      theCFEBData[icfeb] = new CSCCFEBData(icfeb, pos);
      pos += theCFEBData[icfeb]->sizeInWords();
    }
      }
    }   
    pos += theDMBTrailer.sizeInWords();
    size_ = pos-buf;
  }
  else {
    LogTrace ("CSCEventData|CSCRawToDigi") << "Critical Error: DMB Trailer was not found!!! ";
  }
}

std::vector< CSCWireDigi > CSCEventData::wireDigis

(

unsigned

ilayer

)

const

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

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

bool CSCEventData::debug = false

static

Definition at line 137 of file CSCEventData.h.

Referenced by pkg.AbstractPkg::generate(), pkg.AbstractPkg::get_kwds(), setDebug(), unpack_data(), and pkg.AbstractPkg::write().

int CSCEventData::size_

private

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

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

CSCALCTTrailer* CSCEventData::theALCTTrailer

private

Definition at line 164 of file CSCEventData.h.

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

CSCAnodeData* CSCEventData::theAnodeData

private

Definition at line 163 of file CSCEventData.h.

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

CSCCFEBData* CSCEventData::theCFEBData[5]

private

for up to 5 CFEB boards

Definition at line 168 of file CSCEventData.h.

Referenced by add(), cfebData(), copy(), CSCEventData(), destroy(), init(), pack(), 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 175 of file CSCEventData.h.

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

CSCDMBHeader CSCEventData::theDMBHeader

private

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

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

CSCTMBData* CSCEventData::theTMBData

private

Definition at line 165 of file CSCEventData.h.

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

int CSCEventData::zseEnable

private

Definition at line 179 of file CSCEventData.h.

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