EventSummary.cpp

Go to the documentation of this file.

#include "Utility/Logging/Logging.h"
#include <bitset>

#include "Containers/EventSummary.h"

namespace NAIA {
const std::string EventSummary::BranchName{"EventSummaryData"};

float EventSummaryData::TriggerWeight() const {
  float weight = 1.0f;

  if (!IsPhysicsTrigger()) {
    switch (JMembPatt & 0x802) {
    case 0:
      spdlog::error(
          "EventSummaryData::TriggerWeight: Trigger is non-physics ({}) but no unbiased but set in JMembPatt ({})",
          std::bitset<8>(PhysBPatt).to_string(), std::bitset<16>(JMembPatt).to_string());
      throw std::runtime_error("Invalid trigger");
    case 0x2:
      weight = 100.0f;
      break;
    case 0x800:
      weight = 1000.0f;
      break;
    case 0x802:
      weight = 100.0f;
      break;
    default:
      break;
    }
  }

  return weight;
}

void EventSummaryData::RestorePhysBPatt(bool force) {
  // same as in root.C:14634

  if (force || PhysBPatt == 0) {
    if ((JMembPatt & 0x0001) == 0x0001)
      PhysBPatt |= 0x0001; // unbiased 3/4 TOF-CP
    if ((JMembPatt & 0x0090) == 0x0090)
      PhysBPatt |= 0x0002; // 4/4 TOF-CT & ACC0
    if ((JMembPatt & 0x0300) == 0x0300)
      PhysBPatt |= 0x0004; // 4/4 TOF-BZ & ACC1
    if ((JMembPatt & 0x0020) == 0x0020)
      PhysBPatt |= 0x0008; // 4/4 FTZ
    if ((JMembPatt & 0x0410) == 0x0410)
      PhysBPatt |= 0x0010; // 4/4 TOF-CT & ECALF&
    if ((JMembPatt & 0x1000) == 0x1000)
      PhysBPatt |= 0x0020; // ECALA&
    if ((JMembPatt & 0x0800) == 0x0800)
      PhysBPatt |= 0x0040; // unbiased ECALA||
  }
}

void EventSummaryData::Clear() {
  NAntiCluster = 0;
  NTofCluster = 0;
  NTrRecHit = 0;
  NTrTrack = 0;
  NRichRing = 0;
  NEcalShower = 0;
  NParticle = 0;

  HighestInnerChargeParticle = -1;
  HighestEcalEnergyParticle = -1;
  HighestBetaHChargeParticle = -1;

  JMembPatt = std::numeric_limits<unsigned short>::max();
  PhysBPatt = std::numeric_limits<unsigned short>::max();
  NAcc = std::numeric_limits<unsigned short>::max();
}

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

  fmt::print("{:>15} {:>15} {:>15} {:>15} {:>15} {:>15} {:>15} {:>15} {:>10} {:>16}\n", "NAntiCluster", "NTofCluster",
             "NTrRecHit", "NTrTrack", "NRichRing", "NEcalShower", "NParticle", "NAcc", "PhysBPatt", "JMembPatt");

  fmt::print("{:>15} {:>15} {:>15} {:>15} {:>15} {:>15} {:>15} {:>15} {:>10} {:>16}\n", NAntiCluster, NTofCluster,
             NTrRecHit, NTrTrack, NRichRing, NEcalShower, NParticle, NAcc, std::bitset<8>(PhysBPatt).to_string(),
             std::bitset<16>(JMembPatt).to_string());
  fmt::print("{:-^120}\n", "");
}
} // namespace NAIA