BeamSpotOnline.cc

Go to the documentation of this file.

/*
 *   File: DataFormats/Scalers/src/BeamSpotOnline.cc   (W.Badgett)
 */

#include "DataFormats/Scalers/interface/BeamSpotOnline.h"
#include "DataFormats/Scalers/interface/ScalersRaw.h"
#include <cstdio>
#include <ostream>

BeamSpotOnline::BeamSpotOnline() : 
   trigType_(0),
   eventID_(0),
   sourceID_(0),
   bunchNumber_(0),
   version_(0),
   collectionTime_(0,0),
   x_((float)0.0),
   y_((float)0.0),
   z_((float)0.0),
   dxdz_((float)0.0),
   dydz_((float)0.0),
   err_x_((float)0.0),
   err_y_((float)0.0),
   err_z_((float)0.0),
   err_dxdz_((float)0.0),
   err_dydz_((float)0.0),
   width_x_((float)0.0),
   width_y_((float)0.0),
   sigma_z_((float)0.0),
   err_width_x_((float)0.0),
   err_width_y_((float)0.0),
   err_sigma_z_((float)0.0)
{ 
}

BeamSpotOnline::BeamSpotOnline(const unsigned char * rawData)
{ 
  BeamSpotOnline();

  struct ScalersEventRecordRaw_v4 const* raw 
    = reinterpret_cast<struct ScalersEventRecordRaw_v4 const*>(rawData);
  trigType_     = ( raw->header >> 56 ) &        0xFULL;
  eventID_      = ( raw->header >> 32 ) & 0x00FFFFFFULL;
  sourceID_     = ( raw->header >>  8 ) & 0x00000FFFULL;
  bunchNumber_  = ( raw->header >> 20 ) &      0xFFFULL;

  version_ = raw->version;
  if ( version_ >= 4 )
  {
    collectionTime_.set_tv_sec(static_cast<long>(raw->beamSpotOnline.collectionTime_sec));
    collectionTime_.set_tv_nsec(raw->beamSpotOnline.collectionTime_nsec);
    x_           = raw->beamSpotOnline.x;
    y_           = raw->beamSpotOnline.y;
    z_           = raw->beamSpotOnline.z;
    dxdz_        = raw->beamSpotOnline.dxdz;
    dydz_        = raw->beamSpotOnline.dydz;
    err_x_       = raw->beamSpotOnline.err_x;
    err_y_       = raw->beamSpotOnline.err_y;
    err_z_       = raw->beamSpotOnline.err_z;
    err_dxdz_    = raw->beamSpotOnline.err_dxdz;
    err_dydz_    = raw->beamSpotOnline.err_dydz;
    width_x_     = raw->beamSpotOnline.width_x;
    width_y_     = raw->beamSpotOnline.width_y;
    sigma_z_     = raw->beamSpotOnline.sigma_z;
    err_width_x_ = raw->beamSpotOnline.err_width_x;
    err_width_y_ = raw->beamSpotOnline.err_width_y;
    err_sigma_z_ = raw->beamSpotOnline.err_sigma_z;
  }
}

BeamSpotOnline::~BeamSpotOnline() { } 


std::ostream& operator<<(std::ostream& s, const BeamSpotOnline& c) 
{
  char zeit[128];
  constexpr size_t kLineBufferSize = 157;
  char line[kLineBufferSize];
  struct tm * hora;

  s << "BeamSpotOnline    Version: " << c.version() << 
    "   SourceID: "<< c.sourceID() << std::endl;

  timespec ts = c.collectionTime();
  hora = gmtime(&ts.tv_sec);
  strftime(zeit, sizeof(zeit), "%Y.%m.%d %H:%M:%S", hora);
  snprintf(line, kLineBufferSize,
           " CollectionTime: %s.%9.9d", zeit, 
           (int)ts.tv_nsec);
  s << line << std::endl;

  snprintf(line, kLineBufferSize,
           " TrigType: %d   EventID: %d    BunchNumber: %d", 
           c.trigType(), c.eventID(), c.bunchNumber());
  s << line << std::endl;

  snprintf(line, kLineBufferSize,
           "    x: %e +/- %e   width: %e +/- %e",
           c.x(), c.err_x(), c.width_x(), c.err_width_x());
  s << line << std::endl;

  snprintf(line, kLineBufferSize,
           "    y: %e +/- %e   width: %e +/- %e",
           c.y(), c.err_y(), c.width_y(), c.err_width_y());
  s << line << std::endl;

  snprintf(line, kLineBufferSize,
           "    z: %e +/- %e   sigma: %e +/- %e",
           c.z(), c.err_z(), c.sigma_z(), c.err_sigma_z());
  s << line << std::endl;

  snprintf(line, kLineBufferSize,
           " dxdy: %e +/- %e    dydz: %e +/- %e",
           c.dxdz(), c.err_dxdz(), c.dydz(), c.err_dydz());
  s << line << std::endl;
  return s;
}