CSCEventData.cc

Go to the documentation of this file.

#include "EventFilter/CSCRawToDigi/interface/CSCEventData.h"
#include "EventFilter/CSCRawToDigi/interface/CSCCFEBData.h"
#include "DataFormats/CSCDigi/interface/CSCStripDigi.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "EventFilter/CSCRawToDigi/src/cscPackerCompare.h"
#include <iostream>
#include <iterator>
#include "EventFilter/CSCRawToDigi/src/bitset_append.h"
#include "FWCore/Utilities/interface/Exception.h"


#ifdef LOCAL_UNPACK
bool CSCEventData::debug = false;
#else
std::atomic<bool> CSCEventData::debug{false};
#endif

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

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


CSCEventData::CSCEventData(const uint16_t * buf, uint16_t format_version): theFormatVersion(format_version)
{
  theFormatVersion = format_version;
  unpack_data(buf);
}


void CSCEventData::unpack_data(const uint16_t * buf)
{
  // zero everything
  init();
  const uint16_t * 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
                  const uint16_t * 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] = nullptr;
          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;
}

bool CSCEventData::isALCT(const short unsigned int * buf)
{
  return (((buf[0]&0xFFFF)==0xDB0A)||(((buf[0]&0xF800)==0x6000)&&((buf[1]&0xF800)==0)));
}

bool CSCEventData::isTMB(const short unsigned int * buf)
{
  return ((buf[0]&0xFFF)==0xB0C);
}



CSCEventData::CSCEventData(const CSCEventData & data)
{
  copy(data);
}

CSCEventData::~CSCEventData()
{
  destroy();
}


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


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


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


}


void CSCEventData::destroy()
{
  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;
  */
}


std::vector<CSCStripDigi> CSCEventData::stripDigis(const CSCDetId & idlayer) const
{
  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
{
  //  assert(ilayer > 0 && ilayer <= 6); // off because now idlayer is raw cscdetid
  std::vector<CSCStripDigi> result;
  if (theCFEBData[icfeb] != nullptr)
    {
      std::vector<CSCStripDigi> newDigis = theCFEBData[icfeb]->digis(idlayer);
      result.insert(result.end(), newDigis.begin(), newDigis.end());
    }

  return result;
}


std::vector<CSCWireDigi> CSCEventData::wireDigis(unsigned ilayer) const
{
  if (theAnodeData == nullptr)
    {
      return std::vector<CSCWireDigi>();
    }
  else
    {
      return theAnodeData->wireDigis(ilayer);
    }
}


std::vector < std::vector<CSCStripDigi> > CSCEventData::stripDigis() const
{
  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;
}

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


const CSCCFEBData* CSCEventData::cfebData(unsigned icfeb) const
{
  return theCFEBData[icfeb];
}


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

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


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

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


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

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


void CSCEventData::setEventInformation(int bxnum, int lvl1num)
{
  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); 
    }
}


void CSCEventData::checkALCTClasses()
{
  if (theAnodeData == nullptr)
    {
      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()
{
  int nCFEBs = 5;
  if ((theFormatVersion  == 2013) && ((theChamberType == 1) || (theChamberType == 2)) ) {
      nCFEBs = 7;
  }
  if (theTMBData == nullptr)
    {
      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;
}


void CSCEventData::add(const CSCStripDigi & digi, int layer)
{
  //@@ 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] == nullptr)
    {
      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)
{
  checkALCTClasses();
  theAnodeData->add(digi, layer);
  theALCTHeader->setDAVForChannel(digi.getWireGroup());
  theALCTHeader->setBXNCount(digi.getWireGroupBX());
  
}

void CSCEventData::add(const CSCComparatorDigi & digi, int layer)
{
  checkTMBClasses();
  theTMBData->clctData()->add(digi, layer);
  
}

void CSCEventData::add(const CSCComparatorDigi & digi, const CSCDetId & cid)
{
  checkTMBClasses();
  theTMBData->clctData()->add(digi, cid);
  
}



void CSCEventData::add(const std::vector<CSCALCTDigi> & digis)
{
  checkALCTClasses();
  theALCTHeader->add(digis);
}


void CSCEventData::add(const std::vector<CSCCLCTDigi> & digis)
{
  checkTMBClasses();
  theTMBData->tmbHeader()->add(digis);
}

void CSCEventData::add(const std::vector<CSCCorrelatedLCTDigi> & digis)
{
  checkTMBClasses();
  theTMBData->tmbHeader()->add(digis);
}




std::ostream & operator<<(std::ostream & os, const CSCEventData & evt)
{
  for (int ilayer = 1; ilayer <= 6; ++ilayer)
    {
      std::vector<CSCStripDigi> stripDigis = evt.stripDigis(ilayer);
      //copy(stripDigis.begin(), stripDigis.end(), std::ostream_iterator<CSCStripDigi>(os, "\n"));
      //print your scas here
      std::vector<CSCWireDigi> wireDigis = evt.wireDigis(ilayer);
      //copy(wireDigis.begin(), wireDigis.end(), std::ostream_iterator<CSCWireDigi>(os, "\n"));
    }
  return os;
}

boost::dynamic_bitset<> CSCEventData::pack()
{
  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 != nullptr)
    {
      boost::dynamic_bitset<> alctHeader = theALCTHeader->pack();
      result = bitset_utilities::append(result, alctHeader);
      crcvec.push_back(std::make_pair(theALCTHeader->sizeInWords(), theALCTHeader->data()));
    }
  if (theAnodeData != nullptr)
    {
      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 != nullptr)
    {
      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 != nullptr)
    {
      result  = bitset_utilities::append(result, theTMBData->pack());
    }

  for (int icfeb = 0;  icfeb < MAX_CFEB;  ++icfeb)
    {
      if (theCFEBData[icfeb] != nullptr)
        {
          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;
}

unsigned int CSCEventData::calcALCTcrc(std::vector< std::pair<unsigned int, unsigned short*> > &vec)
{
  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 != nullptr) { 
          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;
}



void CSCEventData::selfTest()
{
  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);


}