LTCDigi.cc

Go to the documentation of this file.

#include "DataFormats/LTCDigi/interface/LTCDigi.h"


LTCDigi::LTCDigi(const unsigned char *data)
{
  // six 64 bit words
  cms_uint64_t const *ld = reinterpret_cast<cms_uint64_t const*>(data);

  trigType_   = (ld[0]>>56)&        0xFULL; // 4 bits

  eventID_     = (ld[0]>>32)&0x00FFFFFFULL; // 24 bits
  runNumber_   = (ld[2]>>32)&0xFFFFFFFFULL; // 32 bits
  eventNumber_ = (ld[2])    &0xFFFFFFFFULL; // 32 bits
  
  sourceID_      = (ld[0]>> 8)&0x00000FFFULL; // 12 bits
  // this should always be 815?
  //assert(sourceID_ == 815);

  bunchNumber_   = (ld[0]>>20)&     0xFFFULL; // 12 bits
  orbitNumber_   = (ld[1]>>32)&0xFFFFFFFFULL; // 32 bits

  versionNumber_ = (ld[1]>>24)&0xFFULL;       // 8 bits
  
  daqPartition_  = (ld[1]    )&0xFULL;        // 4 bits


  trigInputStat_ = (ld[3]    )&0xFFFFFFFFULL; // 32 bits

  trigInhibitNumber_ = (ld[3]>>32)&0xFFFFFFFFULL; // 32 bits

  bstGpsTime_    = ld[4]; // 64 bits

}
//static
cms_uint32_t LTCDigi::GetEventNumberFromBuffer(const unsigned char *data) 
{
  // six 64 bit words
  cms_uint64_t const *ld = reinterpret_cast<cms_uint64_t const*>(data);
  cms_uint32_t eventNumber = (ld[2])    &0xFFFFFFFFULL; // 32 bits
  return eventNumber;
}
//static
cms_uint32_t LTCDigi::GetRunNumberFromBuffer(const unsigned char *data) 
{
  // six 64 bit words
  cms_uint64_t const *ld = reinterpret_cast<cms_uint64_t const*>(data);
  cms_uint32_t runNumber   = (ld[2]>>32)&0xFFFFFFFFULL; // 32 bits
  return runNumber;
}


std::ostream & 
operator<<(std::ostream & stream, const LTCDigi & myDigi)
{
   stream << "----------------------------------------"<< std::endl;
   stream << "Dumping LTC digi. " << std::endl;
   stream << "Source ID: " << myDigi.sourceID() << std::endl;
   stream << "Run, event: " << myDigi.runNumber()
      << ", " << myDigi.eventNumber () << std::endl;
   stream << "N_Inhibit:" << myDigi.triggerInhibitNumber() << std::endl;
   stream << LTCDigi::utcTime(myDigi.bstGpsTime()) << std::endl;
   stream << LTCDigi::locTime(myDigi.bstGpsTime()) << std::endl;
   stream << "Partition: " << myDigi.daqPartition() << std::endl;
   stream << "Bunch #:   " << myDigi.bunchNumber()  << std::endl;
   stream << "Orbit #:   " << myDigi.orbitNumber()  << std::endl;

   // Trigger information
   stream << "Trigger Bits(0-5):" ;
   for (int i = 0; i < 6; ++i ) {
      if ( myDigi.HasTriggered(i) )
     stream << "1";
      else 
     stream << "0";
      stream << " ";
   }
   stream << std::endl;

   //
   stream << "Ram trigger: " << myDigi.ramTrigger() << std::endl;
   stream << "VME trigger: " << myDigi.vmeTrigger() << std::endl;

   stream << "++++++++++++++++++++++++++++++++++++++++"<< std::endl;
   
   stream << "Raw Data" << std::endl;
   stream << "Trigger Input status: 0x" 
      << std::hex << myDigi.triggerInputStatus() << std::endl;
   stream << "GPS time:             0x"
      << std::hex << myDigi.bstGpsTime() << std::endl;

   stream << "----------------------------------------"<< std::endl;
   stream << std::dec << std::endl;

   return stream;
}


std::string LTCDigi::utcTime(cms_uint64_t t) //const
{
  // note that gmtime isn't reentrant and we don't own the data
  time_t tsmall = t/1000000;
  tm *utct = gmtime(&tsmall);
  std::string tstr("UTC: ");
  tstr += asctime(utct);
  tstr.replace(tstr.find("\n"),tstr.size(), "");
  return tstr;
}

std::string LTCDigi::locTime(cms_uint64_t t) //const
{
  time_t tsmall = t/1000000;
  std::string a("LOC: ");
  a+= std::string(ctime(&tsmall));
  a.replace(a.find("\n"),a.size(), "");
  return a;
}