#include "CAFAna/Analysis/Calcs.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Core/LoadFromFile.h"
#include "CAFAna/Experiment/SingleSampleExperiment.h"
#include "3FlavorAna/Prediction/PredictionSystJoint2018.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"
#include "CAFAna/Prediction/PredictionExtrap.h"
#include "3FlavorAna/Prediction/PredictionAddRock.h"
#include "CAFAna/Prediction/PredictionCombinePeriods.h"
#include "CAFAna/Analysis/CountingExperimentStatistics.h"
#include "3FlavorAna/Systs/JointAna2018Systs.h"
#include "CAFAna/Systs/XSecSysts.h"
#include "CAFAna/Vars/FitVars.h"
#include "CAFAna/Analysis/Style.h"
#include "3FlavorAna/Plotting/NuePlotStyle.h"
#include "Utilities/func/MathUtil.h"
#include "TCanvas.h"
#include "TBox.h"
#include "TLatex.h"
#include "TColor.h"
#include "TGraph.h"
#include "TVectorD.h"
#include "TF1.h"
#include "TLegend.h"
#include "TLine.h"
#include "TMarker.h"
#include "TStyle.h"
#include "TSystem.h"
#include "TGaxis.h"
#include "OscLib/IOscCalc.h"
#include "TFile.h"
Go to the source code of this file.
Classes | |
class | FitSinSqTheta23LowerOctant |
class | FitSinSqTheta23UpperOctant |
Functions | |
TFile * | GetGaussianSurface (TString par, TString hier) |
void | CylindricalSmoothing (TH2 *h) |
std::string | GetNuePredPath (std::string beam, bool isFakeData, bool NERSC) |
const IPrediction * | GetNuePrediction2018 (std::string decomp, osc::IOscCalc *calc, bool corrSysts, std::string beam, bool isFakeData, bool GetFromUPS=false, bool minimizeMemory=false, bool NERSC=false) |
std::pair< TH1D *, double > | GetNueCosmics2018 (std::string beam, bool GetFromUPS=false, bool NERSC=false) |
Spectrum * | GetNueData2018 (std::string beam) |
std::vector< const IPrediction * > | GetNumuPredictions2018 (const int nq=4, bool useSysts=true, std::string beam="fhc", bool GetFromUPS=false, ENu2018ExtrapType numuExtrap=kNuMu, bool minimizeMemory=false, bool NERSC=false) |
std::vector< std::pair< TH1D *, double > > | GetNumuCosmics2018 (const int nq=4, std::string beam="fhc", bool GetFromUPS=false, bool NERSC=false) |
std::vector< Spectrum * > | GetNumuData2018 (const int nq=4, std::string beam="fhc") |
Spectrum * | GetFakeData (const IPrediction *pred, osc::IOscCalc *calc, const double pot, TH1D *cosmics=0) |
Spectrum * | GetMockData (const IPrediction *pred, osc::IOscCalc *calc, const double pot, TH1D *cosmics=0) |
void | SetFakeCalc (osc::IOscCalcAdjustable *calc, double ssth23=-5, double dmsq32=-5, double dCP_Pi=-5) |
std::vector< const ISyst * > | GetJointFitSystematicList (bool corrSysts, bool nueExclusive=false, bool numuExclusive=false, bool isFHC=true, bool isRHC=true, bool intersection=true) |
std::vector< std::pair< const ISyst *, const ISyst * > > | GetCorrelations (bool isNue, bool isFHC) |
void | DrawSliceCanvas (TString slicename, const double ymax, const double xmin=0, const double xmax=2., bool overlay=false) |
TLegend * | SliceLegend (std::vector< std::pair< TGraph *, TString > > &graphs, bool isDelta, double low=0.7, bool isDmsq=false, bool overlay=false) |
void | DrawContourCanvas (TString surfName, double minx=0, double maxx=2, double miny=0, double maxy=1) |
TLegend * | ContourLegend (int hie, bool fillContour, bool fccorr, Int_t kFillColor1, Int_t kFillColor2, Int_t kFillColor3, Int_t kDarkColor, bool bestFit) |
TLegend * | FCOverlayContourLegend (int hie, bool fillContour, Int_t kFillColor1, Int_t kFillColor2, Int_t kFillColor3, Int_t kDarkColor) |
Variables | |
const IConstrainedFitVar * | kFitSinSqTheta23LowerOctant = new FitSinSqTheta23LowerOctant() |
const IConstrainedFitVar * | kFitSinSqTheta23UpperOctant = new FitSinSqTheta23UpperOctant() |
TLegend* ContourLegend | ( | int | hie, |
bool | fillContour, | ||
bool | fccorr, | ||
Int_t | kFillColor1, | ||
Int_t | kFillColor2, | ||
Int_t | kFillColor3, | ||
Int_t | kDarkColor, | ||
bool | bestFit | ||
) |
Definition at line 643 of file joint_fit_2018_tools.h.
References APDHVSetting::dummy, ana::k2SigmaConfidenceStyle, ana::kBlessedLabelFontSize, and MECModelEnuComparisons::leg.
void CylindricalSmoothing | ( | TH2 * | h | ) |
Definition at line 63 of file joint_fit_2018_tools.h.
References MECModelEnuComparisons::i, calib::j, maxY, minY, APDHVSetting::temp, x1, submit_syst::x2, submit_syst::y, y1, and submit_syst::y2.
Referenced by joint_fit_2018_contours(), joint_fit_2019_contours(), nova_official_data_release(), and plot_joint_fit_2020_contours().
void DrawContourCanvas | ( | TString | surfName, |
double | minx = 0 , |
||
double | maxx = 2 , |
||
double | miny = 0 , |
||
double | maxy = 1 |
||
) |
Definition at line 606 of file joint_fit_2018_tools.h.
References file_size_ana::axes, demo5::c1, ana::kBlessedLabelFontSize, ana::kBlessedTitleFontSize, maxy, plotROC::title, ana::UniqueName(), and ana::XAxisDeltaCPLabels().
void DrawSliceCanvas | ( | TString | slicename, |
const double | ymax, | ||
const double | xmin = 0 , |
||
const double | xmax = 2. , |
||
bool | overlay = false |
||
) |
Definition at line 536 of file joint_fit_2018_tools.h.
References file_size_ana::axes, demo5::c1, ana::kBlessedLabelFontSize, ana::kBlessedTitleFontSize, allInOneTrainingPlots::overlay, plotROC::title, ana::UniqueName(), ana::XAxisDeltaCPLabels(), xmax, and make_mec_shifts_plots::xmin.
TLegend* FCOverlayContourLegend | ( | int | hie, |
bool | fillContour, | ||
Int_t | kFillColor1, | ||
Int_t | kFillColor2, | ||
Int_t | kFillColor3, | ||
Int_t | kDarkColor | ||
) |
Definition at line 696 of file joint_fit_2018_tools.h.
References APDHVSetting::dummy, ana::k2SigmaConfidenceStyle, ana::kBlessedLabelFontSize, and MECModelEnuComparisons::leg.
Definition at line 489 of file joint_fit_2018_tools.h.
References om::cout, ana::GetJointFitSystematicList(), and runNovaSAM::ret.
Spectrum* GetFakeData | ( | const IPrediction * | pred, |
osc::IOscCalc * | calc, | ||
const double | pot, | ||
TH1D * | cosmics = 0 |
||
) |
Definition at line 424 of file joint_fit_2018_tools.h.
References make_associated_cosmic_defs::cosmics, ana::Spectrum::FakeData(), ana::IPrediction::Predict(), and APDHVSetting::temp.
TFile* GetGaussianSurface | ( | TString | par, |
TString | hier | ||
) |
Definition at line 45 of file joint_fit_2018_tools.h.
References om::cout, allTimeWatchdog::endl, runNovaSAM::ret, and cvnie::subdir.
Referenced by nova_official_data_release().
std::vector<const ISyst*> GetJointFitSystematicList | ( | bool | corrSysts, |
bool | nueExclusive = false , |
||
bool | numuExclusive = false , |
||
bool | isFHC = true , |
||
bool | isRHC = true , |
||
bool | intersection = true |
||
) |
Definition at line 452 of file joint_fit_2018_tools.h.
References ana::AddJointAna2018NotCorrelSysts(), ana::AddNonLoadable2018Systs(), om::cout, allTimeWatchdog::endl, ana::getAllAna2018Systs(), isFHC, make_mec_shifts_plots::isRHC, ana::k2ndClassCurrs, ana::kAna2018NormFHC, ana::kAna2018NormRHC, ana::kBoth, ana::kFHC, ana::kJointAna2018, ana::kNueAna2018, ana::kNumuAna2018, ana::kRadCorrNue, ana::kRadCorrNuebar, ana::kRHC, and runNovaSAM::ret.
Spectrum* GetMockData | ( | const IPrediction * | pred, |
osc::IOscCalc * | calc, | ||
const double | pot, | ||
TH1D * | cosmics = 0 |
||
) |
Definition at line 432 of file joint_fit_2018_tools.h.
References make_associated_cosmic_defs::cosmics, ana::Spectrum::MockData(), ana::IPrediction::Predict(), and APDHVSetting::temp.
std::pair<TH1D*, double > GetNueCosmics2018 | ( | std::string | beam, |
bool | GetFromUPS = false , |
||
bool | NERSC = false |
||
) |
Definition at line 240 of file joint_fit_2018_tools.h.
References om::cout, allTimeWatchdog::endl, plot_validation_datamc::fname, cet::getenv(), make_syst_table_plots::h, ana::Spectrum::Integral(), ana::kAna2018FHCLivetime, ana::kAna2018RHCLivetime, ana::kLivetime, livetime, outtime(), runNovaSAM::release, std::sqrt(), string, and ana::Spectrum::ToTH1().
Referenced by bin_composition_pie_chart(), demoFitContours(), demoFitSlices(), fake_future_data(), getHists_FNEX(), joint_fit_2018_contours(), joint_fit_2018_slices(), jointsensitivity(), make_fc_mass_and_oct_nersc_2018(), make_fc_mh_nersc_2018(), make_fc_oct_nersc_2018(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), MakeBkgEMPlot(), MakeCAFSensitivities_for_FNEX(), monoprob(), nuebar_signif(), plot_datamcpred(), sensitivity2018(), and syst_table_fit().
Spectrum* GetNueData2018 | ( | std::string | beam | ) |
Definition at line 281 of file joint_fit_2018_tools.h.
References om::cout, allTimeWatchdog::endl, exit(), ana::Spectrum::Integral(), intime(), ana::kAna2018FHCLivetime, ana::kAna2018FHCPOT, ana::kAna2018RHCLivetime, ana::kAna2018RHCPOT, livetime, ana::Spectrum::Livetime(), POT, ana::Spectrum::POT(), and runNovaSAM::release.
Referenced by joint_fit_2018_contours(), joint_fit_2018_slices(), MakeCAFSensitivities_for_FNEX(), and nuebar_signif().
const IPrediction* GetNuePrediction2018 | ( | std::string | decomp, |
osc::IOscCalc * | calc, | ||
bool | corrSysts, | ||
std::string | beam, | ||
bool | isFakeData, | ||
bool | GetFromUPS = false , |
||
bool | minimizeMemory = false , |
||
bool | NERSC = false |
||
) |
Definition at line 136 of file joint_fit_2018_tools.h.
References POTSpillRate::beam, om::cout, dirName, allTimeWatchdog::endl, exit(), ana::getAllAna2018Systs(), cet::getenv(), GetNuePredPath(), ana::Spectrum::Integral(), ana::kAna2018FHCPOT, ana::kAna2018RHCPOT, ana::kCombo, ana::kFake, ana::kFHC, ana::kNoExtrap, ana::kNueAna2018, ana::kProportional, ana::kRHC, pot, ana::PredictionInterp::Predict(), runNovaSAM::release, rock(), and string.
Referenced by bicount_2018(), bin_composition_pie_chart(), demoFitContours(), demoFitSlices(), fake_future_data(), GetAllPredictions(), getHists_FNEX(), joint_fit_2018_contours(), joint_fit_2018_slices(), jointsensitivity(), make_fc_mass_and_oct_nersc_2018(), make_fc_mh_nersc_2018(), make_fc_oct_nersc_2018(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), MakeCAFSensitivities_for_FNEX(), monoprob(), nuebar_signif(), plot_datamcpred(), sensitivity2018(), syst_table_fit(), systematics_extrap_comp_from_pred_interp(), and systematics_summary_from_pred_interp().
std::string GetNuePredPath | ( | std::string | beam, |
bool | isFakeData, | ||
bool | NERSC | ||
) |
Definition at line 111 of file joint_fit_2018_tools.h.
References dir, cet::getenv(), isFHC, runNovaSAM::ret, and string.
Referenced by GetNuePrediction2018().
std::vector<std::pair <TH1D*, double> > GetNumuCosmics2018 | ( | const int | nq = 4 , |
std::string | beam = "fhc" , |
||
bool | GetFromUPS = false , |
||
bool | NERSC = false |
||
) |
Definition at line 364 of file joint_fit_2018_tools.h.
References POTSpillRate::beam, om::cerr, om::cout, dir, allTimeWatchdog::endl, exit(), shutoffs::filename, cet::getenv(), make_syst_table_plots::h, MECModelEnuComparisons::i, std::sqrt(), string, art::to_string(), and ana::UniqueName().
Referenced by demoFitContours(), demoFitSlices(), fake_future_data(), joint_fit_2018_contours(), joint_fit_2018_slices(), jointsensitivity(), make_fc_mass_and_oct_nersc_2018(), make_fc_mh_nersc_2018(), make_fc_oct_nersc_2018(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), MakeCAFSensitivities_for_FNEX(), sensitivity2018(), and syst_table_fit().
std::vector<Spectrum * > GetNumuData2018 | ( | const int | nq = 4 , |
std::string | beam = "fhc" |
||
) |
Definition at line 404 of file joint_fit_2018_tools.h.
References POTSpillRate::beam, om::cerr, exit(), MakeMiniprodValidationCuts::f, file, MECModelEnuComparisons::i, ana::Spectrum::LoadFrom(), string, and art::to_string().
Referenced by joint_fit_2018_contours(), joint_fit_2018_slices(), and MakeCAFSensitivities_for_FNEX().
std::vector<const IPrediction*> GetNumuPredictions2018 | ( | const int | nq = 4 , |
bool | useSysts = true , |
||
std::string | beam = "fhc" , |
||
bool | GetFromUPS = false , |
||
ENu2018ExtrapType | numuExtrap = kNuMu , |
||
bool | minimizeMemory = false , |
||
bool | NERSC = false |
||
) |
Definition at line 310 of file joint_fit_2018_tools.h.
References POTSpillRate::beam, calc, om::cout, ana::DefaultOscCalc(), shutoffs::filename, plot_validation_datamc::fname, ana::getAllAna2018Systs(), cet::getenv(), ana::kNumuAna2018, and string.
Referenced by demoFitContours(), demoFitSlices(), fake_future_data(), GetAllPredictions(), joint_fit_2018_contours(), joint_fit_2018_slices(), jointsensitivity(), make_fc_mass_and_oct_nersc_2018(), make_fc_mh_nersc_2018(), make_fc_oct_nersc_2018(), make_fc_slices_nersc_2018(), make_fc_slices_nersc_2018_stats(), make_fc_surfaces_nersc_2018(), make_fc_surfaces_nersc_2018_stats(), MakeCAFSensitivities_for_FNEX(), sensitivity2018(), and syst_table_fit().
void SetFakeCalc | ( | osc::IOscCalcAdjustable * | calc, |
double | ssth23 = -5 , |
||
double | dmsq32 = -5 , |
||
double | dCP_Pi = -5 |
||
) |
Definition at line 442 of file joint_fit_2018_tools.h.
References dmsq32, ana::kFitDeltaInPiUnits, ana::kFitDmSq32, ana::kFitSinSqTheta23, ana::FitVarWithPrior::SetValue(), ana::ConstrainedFitVarWithPrior::SetValue(), ana::FitDmSq32::SetValue(), and ssth23.
TLegend* SliceLegend | ( | std::vector< std::pair< TGraph *, TString > > & | graphs, |
bool | isDelta, | ||
double | low = 0.7 , |
||
bool | isDmsq = false , |
||
bool | overlay = false |
||
) |
Definition at line 571 of file joint_fit_2018_tools.h.
References MECModelEnuComparisons::graphs, ana::kBlessedLabelFontSize, MECModelEnuComparisons::leg, allInOneTrainingPlots::overlay, and submit_syst::str.
const IConstrainedFitVar* kFitSinSqTheta23LowerOctant = new FitSinSqTheta23LowerOctant() |
Definition at line 514 of file joint_fit_2018_tools.h.
const IConstrainedFitVar* kFitSinSqTheta23UpperOctant = new FitSinSqTheta23UpperOctant() |
Definition at line 528 of file joint_fit_2018_tools.h.