FileInfo.cpp

Go to the documentation of this file.

#include "Containers/FileInfo.h"

#include "fmt/format.h"

namespace NAIA {
void FileInfo::Clear() {
  FileName = "";
  Run = NEvents = 0;
  UTCTime = EventNo = {0, 0};

  BadRunReason = "";
  BadRunTag = 0;
}

void FileInfo::Dump() const {
  fmt::print("{:=^120} \n", " FileInfo ");

  fmt::print("File: {}\n", FileName);

  fmt::print("{:>15} {:>15} {:^24} {:^24}\n", "Run", "NEvents", "UTCTime", "EventNo");

  fmt::print("{:>15} {:>15} [{:>10}, {:>10}] [{:>10}, {:>10}]\n", Run, NEvents, UTCTime.first, UTCTime.second,
             EventNo.first, EventNo.second);

  fmt::print("{:-^120}\n", "");
}

unsigned int MCFileInfo::GetNGen() const {
  int ngen1 = EventNo.second;
  if (EventNo.first > 1000)
    ngen1 -= EventNo.first;
  int ngen2 = DatacardNGen;
  int ngen = 0;
  // if file and datacard agree at percent, use datacard
  if (fabs(ngen2 - ngen1) < 0.01 * ngen2)
    ngen = ngen2;
  // otherwise use last-first
  else
    ngen = ngen1;
  // if number of generated is too low return 0
  if (ngen < 1000)
    return 0;
  return ngen;
}

double MCFileInfo::GetRMin() const {
  if (Charge == 0)
    return 0;
  // by default take it from datacard
  double rmin = fabs(DatacardMomentumRange.first / Charge);
  return rmin;
}

double MCFileInfo::GetRMax() const {
  if (Charge == 0)
    return 0;
  // by default take it from datacard
  double rmin = fabs(DatacardMomentumRange.second / Charge);
  return rmin;
}

void MCFileInfo::Clear() {
  Charge = DatacardNGen = 0;
  Mass = 0;
  DatacardPID = 0;

  EventNo = {0, 0};
  MomentumRange = DatacardMomentumRange = {0, 0};

  Focus = SimFocus::NONE;
}

void MCFileInfo::Dump() const {
  fmt::print("{:=^120} \n", " MCFileInfo ");

  std::string focusName;
  switch (Focus) {
  case SimFocus::L1:
    focusName = "L1";
    break;
  case SimFocus::L19:
    focusName = "L19";
    break;
  case SimFocus::TB:
    focusName = "TB";
    break;
  default:
    break;
  }

  fmt::print("Focus: {}\n", focusName);
  fmt::print("Charge: {}    Mass: {}\n", Charge, Mass);

  fmt::print("{:>15} {:^24} {:^24}\n", "DatacardPID", "MomentumRange", "DatacardMomentumRange");
  fmt::print("{:>15} [{:>10}, {:>10}] [{:>10},{:>10}]\n", DatacardPID, MomentumRange.first, MomentumRange.second,
             DatacardMomentumRange.first, DatacardMomentumRange.second);

  fmt::print("{:-^120}\n", "");
}

} // namespace NAIA