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 FarDetectorCoreInfo {
  FarDetectorCoreInfo(){}
  FarDetectorCoreInfo(FarDetectorCoreInfo const &other) = delete;
  FarDetectorCoreInfo(FarDetectorCoreInfo &&other) = default;
  FarDetectorCoreInfo& operator=(FarDetectorCoreInfo const &other) = delete;
  FarDetectorCoreInfo& operator=(FarDetectorCoreInfo &&other) = delete;
 
  ~FarDetectorCoreInfo(){}
 
  const int* Target = 0; 
  const int* nupdg  = 0;
  const int* nupdgUnosc = 0;
 
  //THe x_var and y_vars that you're binning in
  const double* x_var = &M3::Unity_D;
  const double* y_var = &M3::Unity_D;
  const double* rw_etru = &M3::_BAD_DOUBLE_;
  const double* rw_truecz = &M3::_BAD_DOUBLE_;
 
  std::vector<const double*> xsec_norm_pointers;
  std::vector<const M3::float_t*> xsec_spline_pointers;
  M3::float_t xsec_w = 1.;
  const M3::float_t* osc_w_pointer = &M3::Unity;
  std::vector<const M3::float_t*> total_weight_pointers;
 
  //M3::float_t total_w  = M3::_BAD_INT_;
 
  int NomXBin = -1;
  int NomYBin = -1;
 
  bool isNC = false;
 
  const double* mode = &M3::Unity_D;
};