LumiScalers.cc

Go to the documentation of this file.

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

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

LumiScalers::LumiScalers()
    : trigType_(0),
      eventID_(0),
      sourceID_(0),
      bunchNumber_(0),
      version_(0),
      normalization_(0.0),
      deadTimeNormalization_(0.0),
      lumiFill_(0.0),
      lumiRun_(0.0),
      liveLumiFill_(0.0),
      liveLumiRun_(0.0),
      instantLumi_(0.0),
      instantLumiErr_(0.0),
      instantLumiQlty_(0),
      lumiETFill_(0.0),
      lumiETRun_(0.0),
      liveLumiETFill_(0.0),
      liveLumiETRun_(0.0),
      instantETLumi_(0.0),
      instantETLumiErr_(0.0),
      instantETLumiQlty_(0),
      lumiOccFill_(nOcc),
      lumiOccRun_(nOcc),
      liveLumiOccFill_(nOcc),
      liveLumiOccRun_(nOcc),
      instantOccLumi_(nOcc),
      instantOccLumiErr_(nOcc),
      instantOccLumiQlty_(nOcc),
      lumiNoise_(nOcc),
      sectionNumber_(0),
      startOrbit_(0),
      numOrbits_(0),
      pileup_(0.0),
      pileupRMS_(0.0),
      bunchLumi_(0.0),
      spare_(0.0) {}

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

  struct ScalersEventRecordRaw_v1 const* raw = reinterpret_cast<struct ScalersEventRecordRaw_v1 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;

  struct LumiScalersRaw_v1 const* lumi = nullptr;

  if (version_ >= 1) {
    if (version_ <= 2) {
      lumi = &(raw->lumi);
    } else {
      struct ScalersEventRecordRaw_v3 const* raw3 = reinterpret_cast<struct ScalersEventRecordRaw_v3 const*>(rawData);
      lumi = &(raw3->lumi);
    }
    collectionTime_.set_tv_sec(static_cast<long>(lumi->collectionTime_sec));
    collectionTime_.set_tv_nsec(lumi->collectionTime_nsec);
    deadTimeNormalization_ = lumi->DeadtimeNormalization;
    normalization_ = lumi->Normalization;
    lumiFill_ = lumi->LumiFill;
    lumiRun_ = lumi->LumiRun;
    liveLumiFill_ = lumi->LiveLumiFill;
    liveLumiRun_ = lumi->LiveLumiRun;
    instantLumi_ = lumi->InstantLumi;
    instantLumiErr_ = lumi->InstantLumiErr;
    instantLumiQlty_ = lumi->InstantLumiQlty;
    lumiETFill_ = lumi->LumiETFill;
    lumiETRun_ = lumi->LumiETRun;
    liveLumiETFill_ = lumi->LiveLumiETFill;
    liveLumiETRun_ = lumi->LiveLumiETRun;
    instantETLumi_ = lumi->InstantETLumi;
    instantETLumiErr_ = lumi->InstantETLumiErr;
    instantETLumiQlty_ = lumi->InstantETLumiQlty;
    for (int i = 0; i < ScalersRaw::N_LUMI_OCC_v1; i++) {
      lumiOccFill_.push_back(lumi->LumiOccFill[i]);
      lumiOccRun_.push_back(lumi->LumiOccRun[i]);
      liveLumiOccFill_.push_back(lumi->LiveLumiOccFill[i]);
      liveLumiOccRun_.push_back(lumi->LiveLumiOccRun[i]);
      instantOccLumi_.push_back(lumi->InstantOccLumi[i]);
      instantOccLumiErr_.push_back(lumi->InstantOccLumiErr[i]);
      instantOccLumiQlty_.push_back(lumi->InstantOccLumiQlty[i]);
      lumiNoise_.push_back(lumi->lumiNoise[i]);
    }
    sectionNumber_ = lumi->sectionNumber;
    startOrbit_ = lumi->startOrbit;
    numOrbits_ = lumi->numOrbits;

    if (version_ >= 7) {
      struct ScalersEventRecordRaw_v6 const* raw6 = (struct ScalersEventRecordRaw_v6 const*)rawData;
      float const* fspare = reinterpret_cast<float const*>(raw6->spare);
      pileup_ = fspare[ScalersRaw::I_SPARE_PILEUP_v7];
      pileupRMS_ = fspare[ScalersRaw::I_SPARE_PILEUPRMS_v7];
      if (version_ >= 8) {
        bunchLumi_ = fspare[ScalersRaw::I_SPARE_BUNCHLUMI_v8];
        spare_ = fspare[ScalersRaw::I_SPARE_SPARE_v8];
      } else {
        bunchLumi_ = (float)0.0;
        spare_ = (float)0.0;
      }
    } else {
      pileup_ = (float)0.0;
      pileupRMS_ = (float)0.0;
      bunchLumi_ = (float)0.0;
      spare_ = (float)0.0;
    }
  }
}

LumiScalers::~LumiScalers() {}

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

  s << "LumiScalers    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,
           " SectionNumber: %10d   StartOrbit: %10d  NumOrbits: %10d",
           c.sectionNumber(),
           c.startOrbit(),
           c.numOrbits());
  s << line << std::endl;

  snprintf(line,
           kLineBufferSize,
           " Normalization: %e  DeadTimeNormalization: %e",
           c.normalization(),
           c.deadTimeNormalization());
  s << line << std::endl;

  // Integrated Luminosity

  snprintf(line, kLineBufferSize, " LumiFill:            %e   LumiRun:            %e", c.lumiFill(), c.lumiRun());
  s << line << std::endl;
  snprintf(
      line, kLineBufferSize, " LiveLumiFill:        %e   LiveLumiRun:        %e", c.liveLumiFill(), c.liveLumiRun());
  s << line << std::endl;

  snprintf(line, kLineBufferSize, " LumiETFill:          %e   LumiETRun:          %e", c.lumiFill(), c.lumiRun());
  s << line << std::endl;

  snprintf(
      line, kLineBufferSize, " LiveLumiETFill:      %e   LiveLumETiRun:      %e", c.liveLumiETFill(), c.liveLumiETRun());
  s << line << std::endl;

  int length = c.instantOccLumi().size();
  for (int i = 0; i < length; i++) {
    snprintf(line,
             kLineBufferSize,
             " LumiOccFill[%d]:      %e   LumiOccRun[%d]:      %e",
             i,
             c.lumiOccFill()[i],
             i,
             c.lumiOccRun()[i]);
    s << line << std::endl;

    snprintf(line,
             kLineBufferSize,
             " LiveLumiOccFill[%d]:  %e   LiveLumiOccRun[%d]:  %e",
             i,
             c.liveLumiOccFill()[i],
             i,
             c.liveLumiOccRun()[i]);
    s << line << std::endl;
  }

  // Instantaneous Luminosity

  snprintf(line,
           kLineBufferSize,
           " InstantLumi:       %e  Err: %e  Qlty: %d",
           c.instantLumi(),
           c.instantLumiErr(),
           c.instantLumiQlty());
  s << line << std::endl;

  snprintf(line,
           kLineBufferSize,
           " InstantETLumi:     %e  Err: %e  Qlty: %d",
           c.instantETLumi(),
           c.instantETLumiErr(),
           c.instantETLumiQlty());
  s << line << std::endl;

  for (int i = 0; i < length; i++) {
    snprintf(line,
             kLineBufferSize,
             " InstantOccLumi[%d]: %e  Err: %e  Qlty: %d",
             i,
             c.instantOccLumi()[i],
             c.instantOccLumiErr()[i],
             c.instantOccLumiQlty()[i]);
    s << line << std::endl;
    snprintf(line, kLineBufferSize, "      LumiNoise[%d]: %e", i, c.lumiNoise()[i]);
    s << line << std::endl;
  }

  snprintf(line, kLineBufferSize, " Pileup: %f       PileupRMS: %f", c.pileup(), c.pileupRMS());
  s << line << std::endl;

  snprintf(line, kLineBufferSize, " BunchLumi: %f    Spare: %f", c.bunchLumi(), c.spare());
  s << line << std::endl;

  return s;
}