FarDetectorCoreInfoStruct.h

Go to the documentation of this file.

#pragma once
#include <vector>
#include <string>
 
#include "Samples/SampleStructs.h"
 
struct OscChannelInfo {
  std::string flavourName;
  std::string flavourName_Latex;
 
  int InitPDG;
  int FinalPDG;
 
  double ChannelIndex;
};
 
inline int GetOscChannel(const std::vector<OscChannelInfo>& OscChannel, const int InitFlav, const int FinalFlav) {
  for (size_t i = 0; i < OscChannel.size(); ++i) {
    if (InitFlav == OscChannel[i].InitPDG && FinalFlav == OscChannel[i].FinalPDG) {
      return static_cast<int>(OscChannel[i].ChannelIndex);
    }
  }
 
  MACH3LOG_ERROR("Didn't find Osc channel for InitFlav = {}, FinalFlav = {}", InitFlav, FinalFlav);
  throw MaCh3Exception(__FILE__, __LINE__);
}
 
struct SampleInfo {
  SampleInfo() = default;
 
  ~SampleInfo() {
    if(DataHist != nullptr) delete DataHist;
    if(MCHist   != nullptr) delete MCHist;
    if(W2Hist   != nullptr) delete W2Hist;
  }
 
  std::vector<std::string> VarStr;
 
  std::string SampleTitle = "";
 
  int nDimensions = M3::_BAD_INT_;
 
  std::vector<std::vector<std::string>> mc_files;
  std::vector<std::string> spline_files;
 
  std::vector<OscChannelInfo> OscChannels;
 
  TH1 *DataHist = nullptr;
  TH1* MCHist = nullptr;
  TH1* W2Hist = nullptr;
};
 
struct EventInfo {
  EventInfo(){}
  EventInfo(EventInfo const &other) = delete;
  EventInfo(EventInfo &&other) = default;
  EventInfo& operator=(EventInfo const &other) = delete;
  EventInfo& operator=(EventInfo &&other) = delete;
  ~EventInfo(){}
 
  std::vector<const M3::float_t*> total_weight_pointers;
 
  std::vector<const double*> KinVar;
  std::vector<int> NomBin;
 
  int NominalSample = M3::_BAD_INT_;
 
  int nupdg  = M3::_BAD_INT_;
  int nupdgUnosc = M3::_BAD_INT_;
  bool isNC = false;
 
  double enu_true = M3::_BAD_DOUBLE_;
  double coszenith_true = M3::_BAD_DOUBLE_;
};