#include "TFile.h"
#include "TH1D.h"
#include "TH2D.h"
#include <TRandom.h>
#include <TMatrixD.h>
#include <TRandom3.h>
#include <vector>
#include <iostream>
#include <string>
Go to the source code of this file.
Functions | |
TH1D * | GetWanted1DFromENu (TH1D *h1In, double &IntegralVal, std::string name) |
TH1D * | GetWanted1DFromQ2 (TH1D *h1In, double &IntegralVal, std::string name) |
TH1D * | GetWanted1DFromMuKin (TH2D *h2In, double &IntegralVal, std::string name) |
TMatrixD | calc_cov_matrix (std::vector< TH1D * > &vhist, TH1D *hIncv) |
TMatrixD | calc_corr_matrix (TMatrixD &mcov) |
TMatrixD | calc_fer_matrix (TMatrixD &mcov, TH1D *hIncv) |
TH1D * | GetSymSyst (TH1D *hInshiftp, TH1D *hInshiftn, TH1D *hIncv, std::string SystName, std::string VarName) |
TH2D * | GetSymSyst (TH2D *hInshiftp, TH2D *hInshiftn, TH2D *hIncv, std::string SystName, std::string VarName) |
void | AddShift (TH1D *hUniv, double &shift, TH1D *hIncv, TH1D *hSystIn) |
void | AddOneSideShift (TH1D *hUniv, double &shift, TH1D *hIncv, TH1D *hSystIn) |
void | scaling (TH1D *hIn, const double shape_scale) |
void | MakeCov (std::string histName, TFile *fIn, TH1D *hIncv, std::vector< TH1D * > &hIngenie, std::vector< TH1D * > &hInppfx, std::vector< TH1D * > &hInsym, std::vector< TH1D * > &hInasym, int thismode) |
void | MakeCovUniverses (std::string VarName, TFile *fIn, TH1D *hIncv, std::vector< TH1D * > &hInuniv, int thismode) |
void | MakeCovSymSyst (std::string varName, int iVar, TFile *fIn, TH1D *hIncv, TH1D *hInsym, int thismode) |
void | MakeCovASymSyst (std::string varName, int iVar, TFile *fIn, TH1D *hIncv, TH1D *hInasym, int thismode) |
void | Calculate_covariances () |
std::vector< int > | starting (const char *fnameIn, const char *fnameonly1DIn, const char *fnameOut, const char *cmode, const char *ccomp) |
void | Calculate_covariances (const char *fnameIn, const char *fnameonly1DIn, const char *fnameOut, const char *cmode, const char *ccomp) |
int | main (int argc, const char *argv[]) |
Variables | |
const int | Nppfx = 100 |
const int | Ngenie = 1000 |
const int | Nuniv = Ngenie*Nppfx |
const int | Nbfoc = 22 |
const std::string | bfoc [Nbfoc] |
const int | Ndsyst = 12 |
const std::string | dsyst [Ndsyst] |
const int | Ncv = 3 |
const std::string | cv [Ncv] = {"cv","nom_good_seventh","neutron_nom"} |
const int | Nwanted_MuKinBins = 172 |
const int | NMuCosBins = 13 |
const int | MuTperCosBins [NMuCosBins] = {6,7,8,9,9,9,13,14,17,20,20,20,20} |
const int | Nwanted_ENuBins = 18 |
const int | Nwanted_Q2Bins = 15 |
const int | Nsymsyst = 16 |
const std::string | symsyst [Nsymsyst] |
const int | basesymsyst_rdm [Nsymsyst] |
const int | Nasymsyst = 2 |
const std::string | asymsyst [Nasymsyst] = {"CalibShape","Cherenkov"} |
const int | baseasymsyst_rdm [Nasymsyst] = {26000000,27000000} |
TRandom3 * | r3 = new TRandom3(0) |
void AddOneSideShift | ( | TH1D * | hUniv, |
double & | shift, | ||
TH1D * | hIncv, | ||
TH1D * | hSystIn | ||
) |
Definition at line 419 of file Calculate_covariances.C.
References plot_xsec_1d::Nbins.
Referenced by MakeCov(), and MakeCovASymSyst().
void AddShift | ( | TH1D * | hUniv, |
double & | shift, | ||
TH1D * | hIncv, | ||
TH1D * | hSystIn | ||
) |
Definition at line 407 of file Calculate_covariances.C.
References plot_xsec_1d::Nbins.
Referenced by ana::SystMaker::GetNominal(), MakeCov(), and MakeCovSymSyst().
Definition at line 352 of file Calculate_covariances.C.
References std::sqrt().
Referenced by MakeCov(), MakeCovASymSyst(), MakeCovSymSyst(), and MakeCovUniverses().
TMatrixD calc_cov_matrix | ( | std::vector< TH1D * > & | vhist, |
TH1D * | hIncv | ||
) |
Definition at line 327 of file Calculate_covariances.C.
Referenced by MakeCov(), MakeCovASymSyst(), MakeCovSymSyst(), and MakeCovUniverses().
void Calculate_covariances | ( | ) |
Definition at line 616 of file Calculate_covariances.C.
References asymsyst, om::cout, allTimeWatchdog::endl, exit(), Nasymsyst, Nsymsyst, and symsyst.
Referenced by main(), and starting().
void Calculate_covariances | ( | const char * | fnameIn, |
const char * | fnameonly1DIn, | ||
const char * | fnameOut, | ||
const char * | cmode, | ||
const char * | ccomp | ||
) |
Definition at line 91 of file Calculate_covariances.C.
References asymsyst, bfoc, om::cout, cv, dsyst, allTimeWatchdog::endl, makeBrightnessMap::fOut, GetSymSyst(), GetWanted1DFromENu(), GetWanted1DFromMuKin(), GetWanted1DFromQ2(), MECModelEnuComparisons::i, MakeCov(), MakeCovASymSyst(), MakeCovSymSyst(), MakeCovUniverses(), Nasymsyst, Nbfoc, Ncv, Ndsyst, Ngenie, Nppfx, Nsymsyst, scaling(), starting(), and symsyst.
TH1D * GetSymSyst | ( | TH1D * | hInshiftp, |
TH1D * | hInshiftn, | ||
TH1D * | hIncv, | ||
std::string | SystName, | ||
std::string | VarName | ||
) |
Definition at line 374 of file Calculate_covariances.C.
References std::abs(), cv, and plot_xsec_1d::Nbins.
Referenced by Calculate_covariances().
TH2D * GetSymSyst | ( | TH2D * | hInshiftp, |
TH2D * | hInshiftn, | ||
TH2D * | hIncv, | ||
std::string | SystName, | ||
std::string | VarName | ||
) |
Definition at line 389 of file Calculate_covariances.C.
References std::abs(), cv, plot_xsec_2d::NbinsX, and plot_xsec_2d::NbinsY.
TH1D * GetWanted1DFromENu | ( | TH1D * | h1In, |
double & | IntegralVal, | ||
std::string | name | ||
) |
Definition at line 299 of file Calculate_covariances.C.
References hadd_many_files::counter, and Nwanted_ENuBins.
Referenced by Calculate_covariances().
TH1D * GetWanted1DFromMuKin | ( | TH2D * | h2In, |
double & | IntegralVal, | ||
std::string | name | ||
) |
Definition at line 279 of file Calculate_covariances.C.
References hadd_many_files::counter, MuTperCosBins, and Nwanted_MuKinBins.
Referenced by Calculate_covariances().
TH1D * GetWanted1DFromQ2 | ( | TH1D * | h1In, |
double & | IntegralVal, | ||
std::string | name | ||
) |
Definition at line 313 of file Calculate_covariances.C.
References hadd_many_files::counter, and Nwanted_Q2Bins.
Referenced by Calculate_covariances().
void MakeCov | ( | std::string | histName, |
TFile * | fIn, | ||
TH1D * | hIncv, | ||
std::vector< TH1D * > & | hIngenie, | ||
std::vector< TH1D * > & | hInppfx, | ||
std::vector< TH1D * > & | hInsym, | ||
std::vector< TH1D * > & | hInasym, | ||
int | thismode | ||
) |
Definition at line 545 of file Calculate_covariances.C.
References std::abs(), AddOneSideShift(), AddShift(), baseasymsyst_rdm, basesymsyst_rdm, calc_corr_matrix(), calc_cov_matrix(), om::cout, allTimeWatchdog::endl, Nasymsyst, Ngenie, Nppfx, Nsymsyst, r3, plot_timing_data::swap, and Write().
Referenced by Calculate_covariances().
void MakeCovASymSyst | ( | std::string | varName, |
int | iVar, | ||
TFile * | fIn, | ||
TH1D * | hIncv, | ||
TH1D * | hInasym, | ||
int | thismode | ||
) |
Definition at line 510 of file Calculate_covariances.C.
References std::abs(), AddOneSideShift(), asymsyst, baseasymsyst_rdm, calc_corr_matrix(), calc_cov_matrix(), om::cout, allTimeWatchdog::endl, Nuniv, and r3.
Referenced by Calculate_covariances().
void MakeCovSymSyst | ( | std::string | varName, |
int | iVar, | ||
TFile * | fIn, | ||
TH1D * | hIncv, | ||
TH1D * | hInsym, | ||
int | thismode | ||
) |
Definition at line 467 of file Calculate_covariances.C.
References AddShift(), basesymsyst_rdm, calc_corr_matrix(), calc_cov_matrix(), calc_fer_matrix(), om::cout, allTimeWatchdog::endl, Nuniv, r3, and symsyst.
Referenced by Calculate_covariances().
void MakeCovUniverses | ( | std::string | VarName, |
TFile * | fIn, | ||
TH1D * | hIncv, | ||
std::vector< TH1D * > & | hInuniv, | ||
int | thismode | ||
) |
Definition at line 435 of file Calculate_covariances.C.
References calc_corr_matrix(), calc_cov_matrix(), om::cout, allTimeWatchdog::endl, and makeTrainCVSamples::int.
Referenced by Calculate_covariances().
void scaling | ( | TH1D * | hIn, |
const double | shape_scale | ||
) |
Definition at line 431 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), genie::BaryonResonanceDecayer::EvolveDeltaDecayWidth(), LoadGenerator(), print_nus17_fd_cut_tables(), and print_nus17_fd_cut_tables2D().
std::vector< int > starting | ( | const char * | fnameIn, |
const char * | fnameonly1DIn, | ||
const char * | fnameOut, | ||
const char * | cmode, | ||
const char * | ccomp | ||
) |
Definition at line 631 of file Calculate_covariances.C.
References asymsyst, Calculate_covariances(), om::cout, allTimeWatchdog::endl, Nasymsyst, Nsymsyst, string, and symsyst.
Referenced by Calculate_covariances().
const std::string asymsyst[Nasymsyst] = {"CalibShape","Cherenkov"} |
Definition at line 63 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), MakeCovASymSyst(), and starting().
Definition at line 64 of file Calculate_covariances.C.
Referenced by MakeCov(), and MakeCovASymSyst().
Definition at line 56 of file Calculate_covariances.C.
Referenced by MakeCov(), and MakeCovSymSyst().
const std::string bfoc[Nbfoc] |
Definition at line 19 of file Calculate_covariances.C.
Referenced by Calculate_covariances().
const std::string cv[Ncv] = {"cv","nom_good_seventh","neutron_nom"} |
Definition at line 39 of file Calculate_covariances.C.
Referenced by BinByBinTemplateFit(), BinByBinTemplateFit_TemplateResults(), Calculate_covariances(), ana::NumuCCIncPionTemplateFitter::CalculateOneSigmaShift(), caf::Proxy< caf::SRFluxWeights >::CheckEquals(), flat::Flat< caf::SRFluxWeights >::Clear(), flat::Flat< caf::SRFluxWeights >::Fill(), ana::FluxMultiverseSyst::FindSigmaBoundaries(), ana::Multiverse::GetNSigmaShift(), ana::GetNSigmaShift(), GetSymSyst(), MakePlots(), caf::Proxy< caf::SRFluxWeights >::operator=(), ppfx_syst_pca_fn(), jmshower::JMClusterMerge::produce(), ana::TrivialCrossSectionAnalysis::RelativeUncertainty(), fxwgt::FluxWeights::SetAdditionalCentralValue(), fxwgt::FluxWeights::SetBeamFocusingCentralValue(), g4rwgt::G4WeightTable::SetCV(), fxwgt::FluxWeights::SetHadronProductionCentralValue(), and ana::FluxMultiverseSyst::ToAreaNormalizedTH1().
const std::string dsyst[Ndsyst] |
Definition at line 31 of file Calculate_covariances.C.
Referenced by Calculate_covariances().
const int MuTperCosBins[NMuCosBins] = {6,7,8,9,9,9,13,14,17,20,20,20,20} |
Definition at line 44 of file Calculate_covariances.C.
Referenced by GetWanted1DFromMuKin().
const int Nasymsyst = 2 |
Definition at line 62 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), MakeCov(), and starting().
const int Nbfoc = 22 |
Definition at line 18 of file Calculate_covariances.C.
Referenced by Calculate_covariances().
const int Ncv = 3 |
Definition at line 38 of file Calculate_covariances.C.
Referenced by Calculate_covariances().
const int Ndsyst = 12 |
Definition at line 30 of file Calculate_covariances.C.
Referenced by Calculate_covariances().
const int Ngenie = 1000 |
Definition at line 16 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), and MakeCov().
const int NMuCosBins = 13 |
Definition at line 43 of file Calculate_covariances.C.
const int Nppfx = 100 |
Definition at line 15 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), and MakeCov().
const int Nsymsyst = 16 |
Definition at line 50 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), MakeCov(), and starting().
Definition at line 17 of file Calculate_covariances.C.
Referenced by ana::DeriveFluxCVAndUniverses(), MakeCovASymSyst(), and MakeCovSymSyst().
const int Nwanted_ENuBins = 18 |
Definition at line 46 of file Calculate_covariances.C.
Referenced by GetWanted1DFromENu().
const int Nwanted_MuKinBins = 172 |
Definition at line 42 of file Calculate_covariances.C.
Referenced by GetWanted1DFromMuKin().
const int Nwanted_Q2Bins = 15 |
Definition at line 47 of file Calculate_covariances.C.
Referenced by GetWanted1DFromQ2().
TRandom3* r3 = new TRandom3(0) |
Definition at line 85 of file Calculate_covariances.C.
Referenced by MakeCov(), MakeCovASymSyst(), MakeCovSymSyst(), and genie::StrumiaVissaniIBDPXSec::MAterm().
const std::string symsyst[Nsymsyst] |
Definition at line 51 of file Calculate_covariances.C.
Referenced by Calculate_covariances(), MakeCovSymSyst(), and starting().