RTIInfo.cpp

Go to the documentation of this file.

#include "Containers/RTIInfo.h"

#include "TMath.h"

#include "fmt/format.h"

namespace NAIA {

TGraph RTIInfo::m_saaEdge = TGraph();

const std::array<float, RTIInfo::nrigbins + 1> RTIInfo::rigbins = {0.,   0.5,  1.0,  1.33, 1.71, 2.15, 2.67,
                                                                   3.29, 3.64, 4.43, 5.37, 6.47, 7.76, 9.26,
                                                                   11.0, 13.0, 15.3, 18.0, 21.1, 24.7, 100000};
const std::array<float, RTIInfo::nbetabins + 1> RTIInfo::betabins = {0.,   0.5, 0.6,  0.7,  0.8,
                                                                     0.85, 0.9, 0.95, 0.99, 0.999};

const std::map<unsigned int, std::pair<float, float>> RTIInfo::qbins = {
    {1, {0.7, 1.3}},   {2, {1.7, 2.3}},   {3, {2.7, 5.3}},   {4, {5.55, 8.45}},
    {5, {9.55, 14.5}}, {6, {14.5, 25.5}}, {7, {25.5, 26.5}}, {8, {26.5, 100}}};

void RTIInfo::UpdateTM() {
  static unsigned long long lastSec = 0;

  if (lastSec != UTCTime) {
    time_t time{UTCTime};
    m_humanTime = *gmtime(&time);
    lastSec = UTCTime;
  }
}

unsigned int RTIInfo::Year() {
  UpdateTM();
  return m_humanTime.tm_year + 1900;
}

unsigned int RTIInfo::Month() {
  UpdateTM();
  return m_humanTime.tm_mon + 1;
}

unsigned int RTIInfo::DayOfMonth() {
  UpdateTM();
  return m_humanTime.tm_mday;
}

unsigned int RTIInfo::DayOfYear() {
  UpdateTM();
  return m_humanTime.tm_yday + 1;
}

unsigned int RTIInfo::Hour() {
  UpdateTM();
  return m_humanTime.tm_hour;
}

unsigned int RTIInfo::Minute() {
  UpdateTM();
  return m_humanTime.tm_min;
}

unsigned int RTIInfo::Second() {
  UpdateTM();
  return m_humanTime.tm_sec;
}

bool RTIInfo::IsInSAA() const {
  if (m_saaEdge.GetN() == 0) {
    m_saaEdge.SetPoint(0, -46.9020531516, 1.07065020944);
    m_saaEdge.SetPoint(1, -42.4324570724, 0.278767621219);
    m_saaEdge.SetPoint(2, -38.7963908651, -0.864383087583);
    m_saaEdge.SetPoint(3, -32.3643804368, -7.86605177782);
    m_saaEdge.SetPoint(4, -32.1564301222, -8.2925091686);
    m_saaEdge.SetPoint(5, -30.4732034929, -10.2380772476);
    m_saaEdge.SetPoint(6, -28.9780074759, -11.1541721933);
    m_saaEdge.SetPoint(7, -24.9729186857, -13.9219707474);
    m_saaEdge.SetPoint(8, -22.3514515118, -15.8035549547);
    m_saaEdge.SetPoint(9, -17.864105537, -17.5329648531);
    m_saaEdge.SetPoint(10, -13.7360395887, -18.4514556443);
    m_saaEdge.SetPoint(11, -9.31195234438, -20.5579164624);
    m_saaEdge.SetPoint(12, -7.18498632474, -21.1949557893);
    m_saaEdge.SetPoint(13, -3.59895538649, -23.8036939586);
    m_saaEdge.SetPoint(14, -2.34504587932, -24.8012082521);
    m_saaEdge.SetPoint(15, -0.819716509312, -27.065064314);
    m_saaEdge.SetPoint(16, -0.0558141297819, -29.4952077252);
    m_saaEdge.SetPoint(17, -1.39465186809, -31.1867845581);
    m_saaEdge.SetPoint(18, -4.48059829104, -33.3607272689);
    m_saaEdge.SetPoint(19, -8.42183566883, -35.9587368092);
    m_saaEdge.SetPoint(20, -9.27876931794, -36.2964186885);
    m_saaEdge.SetPoint(21, -13.8165106415, -37.7311180334);
    m_saaEdge.SetPoint(22, -17.4924677527, -39.0468547418);
    m_saaEdge.SetPoint(23, -18.408230558, -39.407844502);
    m_saaEdge.SetPoint(24, -19.3555119436, -39.7809321206);
    m_saaEdge.SetPoint(25, -26.2402929635, -42.5628157663);
    m_saaEdge.SetPoint(26, -29.3390198221, -43.7512422146);
    m_saaEdge.SetPoint(27, -30.4468807826, -44.1487676525);
    m_saaEdge.SetPoint(28, -36.2608202879, -44.9657315842);
    m_saaEdge.SetPoint(29, -38.6471514882, -45.8362408379);
    m_saaEdge.SetPoint(30, -39.884535768, -46.2670553531);
    m_saaEdge.SetPoint(31, -41.1537099863, -46.7007904009);
    m_saaEdge.SetPoint(32, -46.4231376314, -47.035405079);
    m_saaEdge.SetPoint(33, -50.5517527888, -46.9389068642);
    m_saaEdge.SetPoint(34, -56.1886816966, -45.9179669008);
    m_saaEdge.SetPoint(35, -60.4620401598, -44.5064118787);
    m_saaEdge.SetPoint(36, -65.9797837207, -40.9661101604);
    m_saaEdge.SetPoint(37, -68.6756029734, -39.1105256272);
    m_saaEdge.SetPoint(38, -69.9926612952, -38.1896928803);
    m_saaEdge.SetPoint(39, -72.5874587403, -35.1855408369);
    m_saaEdge.SetPoint(40, -74.9768865609, -32.1834974998);
    m_saaEdge.SetPoint(41, -79.386936716, -24.8704780079);
    m_saaEdge.SetPoint(42, -81.5015374707, -20.6638534546);
    m_saaEdge.SetPoint(43, -83.3018473547, -16.5956399645);
    m_saaEdge.SetPoint(44, -83.4855953538, -15.6706377827);
    m_saaEdge.SetPoint(45, -83.4016563744, -13.3880324461);
    m_saaEdge.SetPoint(46, -79.0157111546, -9.8025691601);
    m_saaEdge.SetPoint(47, -76.0062621477, -8.49383318946);
    m_saaEdge.SetPoint(48, -71.9675022979, -7.9008874195);
    m_saaEdge.SetPoint(49, -69.1989836702, -6.34481990494);
    m_saaEdge.SetPoint(50, -68.320940027, -6.15577355353);
    m_saaEdge.SetPoint(51, -60.8270728968, -4.71604157559);
    m_saaEdge.SetPoint(52, -57.4730300335, -3.41729725291);
    m_saaEdge.SetPoint(53, -53.25104737, -1.64190214253);
    m_saaEdge.SetPoint(54, -51.5811868613, -0.10182123818);
    m_saaEdge.SetPoint(55, -49.9470795728, 1.16581717215);
    m_saaEdge.SetPoint(56, -46.9020531516, 1.07065020944);
  }

  float rtiLong = Phi * TMath::RadToDeg();
  float rtiLat = Theta * TMath::RadToDeg();
  while (rtiLong > 180.0f)
    rtiLong -= 360.0f;

  return m_saaEdge.IsInside(rtiLong, rtiLat);
};

void RTIInfo::Clear() {
  Run = RunTag = FirstEvNo = LastEvNo = good = UTCTime = 0;
  LivetimeFraction = MostProbableHeRig = Theta = Phi = Altitude = Zenith = GalacticLat = GalacticLong = nEvent =
      nError = nTrigger = nHWError = nEventParticle = nTRDHitAvg = nTrClusterAvg = 0.0f;

  DSPError = DSPErrorJLV1 = DSPErrorJINJ = DSPErrorJINF = false;

  MagThetaSphCoo = MagPhiSphCoo = MagAltSphCoo = LShell = 0;

  nTRDTrack = nTRDHits = 0;
  nEcalShower = nEcalHits = 0;
  nTrClusterX = nTrClusterY = nTrClusterL1X = nTrClusterL1Y = 0;
  nTofClusterH = 0;
  nParticle = nGoodParticle = nPartTrk = nPartTrkInn = 0;
  nTrk = nTrkInn = 0;
  nACC = 0;
  TofRaw_nftdc = TofRaw_nstdc = TofRaw_nsumh = TofRaw_nsumsh = 0;
  AntiRaw_nftdc = AntiRaw_ntdct = 0;
  nRoomError = 0;

  // reset arrays
  decltype(MaxStoermerCutoff){0}.swap(MaxStoermerCutoff);
  decltype(MaxIGRFCutoff){0}.swap(MaxIGRFCutoff);
  decltype(MinIGRFCutoff){0}.swap(MinIGRFCutoff);
  decltype(nEventInnerL1L9){0}.swap(nEventInnerL1L9);
  decltype(MeanAlignDiffExtLayer){0}.swap(MeanAlignDiffExtLayer);

  decltype(nTRDLayerHits){0}.swap(nTRDLayerHits);
  decltype(nPartRate){0}.swap(nPartRate);
  decltype(nPartInnRate){0}.swap(nPartInnRate);
  decltype(nPartBetaRate){0}.swap(nPartBetaRate);
  decltype(nPartAboveCutoff){0}.swap(nPartAboveCutoff);
  decltype(nPartBelowCutoff){0}.swap(nPartBelowCutoff);
  decltype(nTrkRate){0}.swap(nTrkRate);
  decltype(nTrkInnRate){0}.swap(nTrkInnRate);
  decltype(TrigPhysBPatt){0}.swap(TrigPhysBPatt);
  decltype(TrigJMembPatt){0}.swap(TrigJMembPatt);
  decltype(TrigRates){0}.swap(TrigRates);
  decltype(TofFlags){0}.swap(TofFlags);
}

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

  fmt::print("{:>15} {:>15} {:>15} {:>20} {:>15} {:>15} {:>15} {:>15} {:>15}\n", "Run", "RunTag", "UTCTime",
             "LivetimeFraction", "Theta", "Phi", "Altitude", "nEvent", "nError");

  fmt::print("{:>15} {:>15x} {:>15} {:>20} {:>15} {:>15} {:>15} {:>15} {:>15}\n", Run, RunTag, UTCTime,
             LivetimeFraction, Theta, Phi, Altitude, nEvent, nError);
  fmt::print("{:-^120}\n", "");
}

} // namespace NAIA