Classes | Functions
test_predInterp_debug.C File Reference
#include "CAFAna/Analysis/Calcs.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Analysis/Plots.h"
#include "OscLib/IOscCalc.h"
#include "OscLib/OscCalcDumb.h"
#include "CAFAna/nue/Ana2018/joint_fit_2018_tools.h"
#include "3FlavorAna/Systs/JointAna2018Systs.h"
#include "TCanvas.h"
#include <iomanip>

Go to the source code of this file.

Classes

struct  ColumnDef
 

Functions

double Integral (const Spectrum &s, double pot)
 
void test_predInterp_debug (std::string beam="fhc")
 

Function Documentation

double Integral ( const Spectrum s,
double  pot 
)

Definition at line 32 of file test_predInterp_debug.C.

References ana::Spectrum::Integral().

33 {
34  return s.Integral(pot);
35 }
double Integral(double exposure, double *err=0, EExposureType expotype=kPOT) const
Return total number of events scaled to pot.
Definition: Spectrum.cxx:272
#define pot
void test_predInterp_debug ( std::string  beam = "fhc")

Definition at line 37 of file test_predInterp_debug.C.

References POTSpillRate::beam, cols, febshutoff_auto::curr, ana::getAllAna2018Systs(), ana::Flavors::kAll, ana::Current::kBoth, ana::Sign::kBoth, ana::Current::kCC, ana::kCombo, ana::kMAQEGenieReducedSyst2018, ana::Current::kNC, ana::kNueAna2018, ana::Flavors::kNuEToNuE, ana::Flavors::kNuMuToNuE, ana::Flavors::kNuMuToNuMu, ana::Flavors::kNuMuToNuTau, ana::kProportional, plot_validation_datamc::pred, canMan::sign(), and systs.

38 {
39  std::vector<ColumnDef> cols =
40  {{"$\\nu_e$ Signal", Flavors::kNuMuToNuE, Current::kCC, Sign::kBoth},
41  {"Total Beam Bkg", Flavors::kAll, Current::kBoth, Sign::kBoth},
42  {"beam $\\nu_e$ CC", Flavors::kNuEToNuE, Current::kCC, Sign::kBoth},
43  {"NC", Flavors::kAll, Current::kNC, Sign::kBoth},
44  {"$\\nu_\\mu$ CC", Flavors::kNuMuToNuMu, Current::kCC, Sign::kBoth},
45  {"$\\nu_\\tau$ CC", Flavors::kNuMuToNuTau, Current::kCC, Sign::kBoth}
46  };
47 
48  osc::OscCalcDumb dumbosc;
49 
50  auto systs = getAllAna2018Systs(kNueAna2018,false);
51 
53  if(beam == "fhc") extrap = kCombo;
54  if(beam == "rhc") extrap = kProportional;
55 
56  auto *pred = new PredictionSystJoint2018(extrap,&dumbosc,{&kMAQEGenieReducedSyst2018}/*systs*/,beam,false, true,"",false,{&kMAQEGenieReducedSyst2018, &kDISvpCC0pi});
57 
58  pred->DebugPlots(&dumbosc,("debugPlots/debug_signal_"+beam+"_%s.pdf").c_str(),cols[0].flav,cols[0].curr,cols[0].sign);
59  pred->DebugPlots(&dumbosc,("debugPlots/debug_totalbkg_"+beam+"_%s.pdf").c_str(),cols[1].flav,cols[1].curr,cols[1].sign);
60 }
std::vector< SystGroupDef > systs
Definition: syst_header.h:385
(&#39; appearance&#39;)
Definition: IPrediction.h:18
Loads shifted spectra from files.
(&#39;beam &#39;)
Definition: IPrediction.h:15
const NOvARwgtSyst kMAQEGenieReducedSyst2018(genie::rew::GSyst::AsString(genie::rew::EGSyst(rwgt::fReweightMaCCQE))++"_reduced_2018", genie::rew::GSyst::AsString(genie::rew::EGSyst(rwgt::fReweightMaCCQE))++"_reduced_2018", novarwgt::kMAQEGenieReducedSyst2018)
2018 &#39;reduced&#39; M_A^QE shift. See documentation in NOvARwgt
Definition: XSecSysts.h:51
std::vector< const ISyst * > getAllAna2018Systs(const EAnaType2018 ana, const bool smallgenie, const BeamType2018 beam, bool isFit)
Charged-current interactions.
Definition: IPrediction.h:39
const int cols[3]
Simple oscillation probability calculator that has no solar term or mass hierarchy or delta so it&#39;s s...
Definition: OscCalcDumb.h:16
(&#39; survival&#39;)
Definition: IPrediction.h:19
All neutrinos, any flavor.
Definition: IPrediction.h:26
def sign(x)
Definition: canMan.py:197