#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/ISyst.h"
#include "CAFAna/Core/SystShifts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "CAFAna/Systs/XSecSystLists.h"
#include "3FlavorAna/Cuts/NueCutsSecondAna.h"
#include "3FlavorAna/Vars/NueVars.h"
#include "3FlavorAna/Vars/NueVarsExtra.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/XSec/GenieMultiverseSyst.h"
#include "CAFAna/XSec/FluxMultiverseSyst.h"
#include "CAFAna/XSec/Flux.h"
#include "CAFAna/XSec/TargetCount.h"
#include "CAFAna/Core/ReweightableSpectrum.h"
#include "CAFAna/Unfold/UnfoldIterative.h"
#include "CAFAna/Analysis/Plots.h"
#include "Utilities/rootlogon.C"
#include "TH1.h"
#include "TH2.h"
#include "THStack.h"
#include "TCanvas.h"
#include "TLegend.h"
#include "TStyle.h"
#include "TFile.h"
#include <vector>
#include <math.h>
#include "TLine.h"
#include "TTree.h"
#include "TF1.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TF2.h"
#include "TCutG.h"
#include "TMath.h"
#include "TRandom.h"
#include "TGraph.h"
#include "TColor.h"
#include "TF3.h"
#include "TH3.h"
#include "TH3F.h"
#include "TMinuit.h"
#include "TMatrix.h"
#include "TGraphErrors.h"
#include "TDecompLU.h"
#include "TDecompSVD.h"
#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
#include <string>
#include <stdint.h>
#include <sstream>
Go to the source code of this file.
Functions | |
std::string | getCaption2D (TH3F *sample, int xbin, int ybin) |
void | PrintCaption (TString str) |
void | fcn (Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag) |
void | fcn2Var (Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag) |
void | MakeWeightAndErrorPlot (TH2F *hWeights, TH2F *hErr, std::string pidName, std::string dataName, std::string outnameWei, std::string outnameErr, std::string ytitle, std::string xtitle, std::string title, std::string ztitle_wei, std::string ztitle_err, float x_lo=0.75, float x_hi=1.0, float y_lo=1.0, float y_hi=10.0) |
void | PlotFitResults (TH1F *hData, std::vector< TH1F * >hMCnom, std::vector< TH1F * >hMCwei, TGraphAsymmErrors *g, std::string pidName, std::string dataName, std::string caption, int xbin, int ybin) |
TGraphAsymmErrors * | PlotSystErrorBand (TH1F *&nom, std::vector< TH1F * > &ups, std::vector< TH1F * > &dns, int col, int errCol, float headroom, bool newaxis) |
void | CalculateCovarianceMatrix (TH1F *nom, std::vector< TH1F * > hSysts, int nPseudo, std::string pidName, int xbin, int ybin) |
void | CalculateTotalCovariance (TH1F *nominal, std::vector< TH1F * > hSysts, std::string pidName, int xbin, int ybin) |
std::vector< TH3F * > | BinByBinTemplateFit (TH3F *data, std::vector< TH3F * > templates, std::vector< TH3F * > systs_hists, std::vector< TH3F * > systs_hists_up, std::vector< TH3F * > systs_hists_down, std::vector< TH3F * > analysis_templates, std::string pidName, std::string varName, std::string dataName) |
std::vector< TH3F * > | BinByBinTemplateFit_TemplateResults (TH3F *data, std::vector< TH3F * > templates, std::vector< TH3F * > systs_hists, std::vector< TH3F * > systs_hists_up, std::vector< TH3F * > systs_hists_down, std::vector< TH3F * > analysis_templates, std::string pidName, std::string varName, std::string dataName) |
Variables | |
const int | pidBins = 21 |
const bool | shouldRebin = false |
double | covMatrix [pidBins][pidBins] |
double | stat_covMatrix [pidBins][pidBins] |
const float | s_b_ratio = 0.40 |
const float | signalRegion = 0.85 |
float | fDA [pidBins] |
float | fSigLike [pidBins] |
float | fNumuLike [pidBins] |
float | fNC [pidBins] |
float | fOther [pidBins] |
const std::vector< int > | kNueCCColorDef |
std::vector<TH3F*> BinByBinTemplateFit | ( | TH3F * | data, |
std::vector< TH3F * > | templates, | ||
std::vector< TH3F * > | systs_hists, | ||
std::vector< TH3F * > | systs_hists_up, | ||
std::vector< TH3F * > | systs_hists_down, | ||
std::vector< TH3F * > | analysis_templates, | ||
std::string | pidName, | ||
std::string | varName, | ||
std::string | dataName | ||
) |
Definition at line 721 of file NueCCIncTemplateFit.h.
References bin, CalculateCovarianceMatrix(), CalculateTotalCovariance(), make_true_q0q3_plots::caption, om::cerr, plot_validation_datamc::Clone(), om::cout, covariance(), cv, allTimeWatchdog::endl, stan::math::fabs(), fcn(), fcn2Var(), fDA, stan::math::fmin(), fNC, fNumuLike, fOther, fSigLike, getCaption2D(), hi(), MECModelEnuComparisons::i, makeTrainCVSamples::int, std::isnan(), kNueCCColorDef, lo(), make_pair(), MakeWeightAndErrorPlot(), nc, ana::nuebar, ana::numubar, fhicl::other, PlotFitResults(), PlotSystErrorBand(), cet::pow(), s_b_ratio, signalRegion, std::sqrt(), stat_covMatrix, string, ana::UniqueName(), submit_syst::x, and test::z.
std::vector<TH3F*> BinByBinTemplateFit_TemplateResults | ( | TH3F * | data, |
std::vector< TH3F * > | templates, | ||
std::vector< TH3F * > | systs_hists, | ||
std::vector< TH3F * > | systs_hists_up, | ||
std::vector< TH3F * > | systs_hists_down, | ||
std::vector< TH3F * > | analysis_templates, | ||
std::string | pidName, | ||
std::string | varName, | ||
std::string | dataName | ||
) |
Definition at line 1283 of file NueCCIncTemplateFit.h.
References bin, CalculateCovarianceMatrix(), CalculateTotalCovariance(), make_true_q0q3_plots::caption, om::cerr, plot_validation_datamc::Clone(), om::cout, covariance(), cv, allTimeWatchdog::endl, stan::math::fabs(), fcn(), fcn2Var(), fDA, stan::math::fmin(), fNC, fNumuLike, fOther, fSigLike, getCaption2D(), hi(), MECModelEnuComparisons::i, makeTrainCVSamples::int, std::isnan(), kNueCCColorDef, lo(), make_pair(), MakeWeightAndErrorPlot(), nc, ana::nuebar, ana::numubar, fhicl::other, PlotFitResults(), PlotSystErrorBand(), cet::pow(), s_b_ratio, signalRegion, std::sqrt(), stat_covMatrix, string, ana::UniqueName(), submit_syst::x, and test::z.
void CalculateCovarianceMatrix | ( | TH1F * | nom, |
std::vector< TH1F * > | hSysts, | ||
int | nPseudo, | ||
std::string | pidName, | ||
int | xbin, | ||
int | ybin | ||
) |
Definition at line 467 of file NueCCIncTemplateFit.h.
References demo5::c1, hadd_many_files::counter, covariance(), covMatrix, MECModelEnuComparisons::i, calib::j, confusionMatrixTree::out, and ana::UniqueName().
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
void CalculateTotalCovariance | ( | TH1F * | nominal, |
std::vector< TH1F * > | hSysts, | ||
std::string | pidName, | ||
int | xbin, | ||
int | ybin | ||
) |
Definition at line 541 of file NueCCIncTemplateFit.h.
References demo5::c1, demo5::c2, covariance(), MECModelEnuComparisons::i, makeTrainCVSamples::int, calib::j, PandAna.Demos.pi0_spectra::labels, confusionMatrixTree::out, std::sqrt(), and ana::UniqueName().
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
void fcn | ( | Int_t & | npar, |
Double_t * | gin, | ||
Double_t & | f, | ||
Double_t * | par, | ||
Int_t | iflag | ||
) |
Definition at line 134 of file NueCCIncTemplateFit.h.
References covMatrix, fDA, fNC, fNumuLike, fOther, fSigLike, MECModelEnuComparisons::i, calib::j, nc, fhicl::other, ana::pidBins, and stat_covMatrix.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
void fcn2Var | ( | Int_t & | npar, |
Double_t * | gin, | ||
Double_t & | f, | ||
Double_t * | par, | ||
Int_t | iflag | ||
) |
Definition at line 195 of file NueCCIncTemplateFit.h.
References covMatrix, fDA, fNC, fNumuLike, fOther, fSigLike, MECModelEnuComparisons::i, calib::j, nc, fhicl::other, ana::pidBins, and stat_covMatrix.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
std::string getCaption2D | ( | TH3F * | sample, |
int | xbin, | ||
int | ybin | ||
) |
Definition at line 103 of file NueCCIncTemplateFit.h.
References make_true_q0q3_plots::caption, and string.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
void MakeWeightAndErrorPlot | ( | TH2F * | hWeights, |
TH2F * | hErr, | ||
std::string | pidName, | ||
std::string | dataName, | ||
std::string | outnameWei, | ||
std::string | outnameErr, | ||
std::string | ytitle, | ||
std::string | xtitle, | ||
std::string | title, | ||
std::string | ztitle_wei, | ||
std::string | ztitle_err, | ||
float | x_lo = 0.75 , |
||
float | x_hi = 1.0 , |
||
float | y_lo = 1.0 , |
||
float | y_hi = 10.0 |
||
) |
Definition at line 251 of file NueCCIncTemplateFit.h.
References demo5::c1, demo5::c2, confusionMatrixTree::out, and ana::UniqueName().
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
void PlotFitResults | ( | TH1F * | hData, |
std::vector< TH1F * > | hMCnom, | ||
std::vector< TH1F * > | hMCwei, | ||
TGraphAsymmErrors * | g, | ||
std::string | pidName, | ||
std::string | dataName, | ||
std::string | caption, | ||
int | xbin, | ||
int | ybin | ||
) |
Definition at line 299 of file NueCCIncTemplateFit.h.
References ana::AutoPlaceLegend(), plot_validation_datamc::Clone(), make_mec_shifts_plots::GetMaximum(), GetYaxis(), kBlue, kGreen, kRed, MECModelEnuComparisons::leg, pad1, pad2, PrintCaption(), Simulation(), and ana::UniqueName().
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
TGraphAsymmErrors* PlotSystErrorBand | ( | TH1F *& | nom, |
std::vector< TH1F * > & | ups, | ||
std::vector< TH1F * > & | dns, | ||
int | col, | ||
int | errCol, | ||
float | headroom, | ||
bool | newaxis | ||
) |
Definition at line 417 of file NueCCIncTemplateFit.h.
References MECModelEnuComparisons::g, hi(), kRed, lo(), std::sqrt(), std::swap(), w, and submit_syst::y.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
void PrintCaption | ( | TString | str | ) |
Definition at line 124 of file NueCCIncTemplateFit.h.
References prelim.
Referenced by PlotFitResults().
Definition at line 78 of file NueCCIncTemplateFit.h.
Referenced by CalculateCovarianceMatrix(), fcn(), and fcn2Var().
float fDA[pidBins] |
Definition at line 85 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), fcn(), and fcn2Var().
float fNC[pidBins] |
Definition at line 88 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), fcn(), and fcn2Var().
float fNumuLike[pidBins] |
Definition at line 87 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), fcn(), and fcn2Var().
float fOther[pidBins] |
Definition at line 89 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), fcn(), and fcn2Var().
float fSigLike[pidBins] |
Definition at line 86 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), fcn(), and fcn2Var().
const std::vector<int> kNueCCColorDef |
Definition at line 91 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
const int pidBins = 21 |
Definition at line 76 of file NueCCIncTemplateFit.h.
const float s_b_ratio = 0.40 |
Definition at line 80 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
const bool shouldRebin = false |
Definition at line 77 of file NueCCIncTemplateFit.h.
const float signalRegion = 0.85 |
Definition at line 83 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), and BinByBinTemplateFit_TemplateResults().
Definition at line 79 of file NueCCIncTemplateFit.h.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), fcn(), and fcn2Var().