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(dathist   != nullptr) delete dathist;
    if(dathist2d != nullptr) delete dathist2d;
    if(_hPDF1D   != nullptr) delete _hPDF1D;
    if(_hPDF2D   != nullptr) delete _hPDF2D;
  }
 
  std::string XVarStr = "";
  std::string YVarStr = "";
 
  std::string SampleTitle = "";
 
  int nDimensions = M3::_BAD_INT_;
 
  TH1D *dathist = nullptr;
  TH2D *dathist2d = nullptr;
 
  TH1D* _hPDF1D = nullptr;
  TH2D* _hPDF2D = nullptr;
 
  void InitialiseHistograms() {
    TString histname1d = (XVarStr).c_str();
    TString histname2d = (XVarStr + "_" + YVarStr).c_str();
    TString histtitle = SampleTitle;
 
    //The binning here is arbitrary, now we get info from cfg so the
    //set1DBinning and set2Dbinning calls below will make the binning
    //to be what we actually want
    _hPDF1D   = new TH1D("h" + histname1d + SampleTitle, histtitle, 1, 0, 1);
    dathist   = new TH1D("d" + histname1d + SampleTitle, histtitle, 1, 0, 1);
    _hPDF2D   = new TH2D("h" + histname2d + SampleTitle, histtitle, 1, 0, 1, 1, 0, 1);
    dathist2d = new TH2D("d" + histname2d + SampleTitle, histtitle, 1, 0, 1, 1, 0, 1);
 
    _hPDF1D->GetXaxis()->SetTitle(XVarStr.c_str());
    dathist->GetXaxis()->SetTitle(XVarStr.c_str());
    _hPDF2D->GetXaxis()->SetTitle(XVarStr.c_str());
    dathist2d->GetXaxis()->SetTitle(XVarStr.c_str());
  }
};
 
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;
};