Calculate_covariances.C File Reference

#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)

Function Documentation

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().

                                                                          {
  int Nbins = hUniv->GetXaxis()->GetNbins();
  for(int ii=1;ii<=Nbins;ii++){
    double univ_val   = hUniv->GetBinContent(ii);
    double cv_val     = hIncv->GetBinContent(ii);
    double syst_val   = hSystIn->GetBinContent(ii);    
    double shift_val  = syst_val - cv_val;
    shift_val *= shift;    
    univ_val += shift_val;
    hUniv->SetBinContent(ii,univ_val);
  }     
}

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().

                                                                   {
  int Nbins = hUniv->GetXaxis()->GetNbins();
  for(int ii=1;ii<=Nbins;ii++){
    double univ_val   = hUniv->GetBinContent(ii);
    double cv_val     = hIncv->GetBinContent(ii);
    double syst_val   = hSystIn->GetBinContent(ii);    
    double shift_val  = syst_val - cv_val;
    shift_val *= shift;    
    univ_val += shift_val;
    hUniv->SetBinContent(ii,univ_val);
  }     
}

TMatrixD calc_corr_matrix

(

TMatrixD &

mcov

)

Definition at line 352 of file Calculate_covariances.C.

References std::sqrt().

Referenced by MakeCov(), MakeCovASymSyst(), MakeCovSymSyst(), and MakeCovUniverses().

                                         {
  TMatrixD cormx(mcov.GetNrows(),mcov.GetNcols());
  
  for(int ii=0; ii<mcov.GetNrows(); ii++){
    for(int jj=0; jj<mcov.GetNcols(); jj++){
      cormx[ii][jj] = mcov[ii][jj]/(sqrt(mcov[ii][ii])*sqrt(mcov[jj][jj]));
    }
  }    
  return cormx;    
}

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().

                                                              {
  
  unsigned int N =  vhist.size();    
  unsigned int Nbinsx = vhist[0]->GetNbinsX();
  TMatrixD covmx(Nbinsx,Nbinsx);  
  //universes
  for(int iu=0; iu<N; iu++){    
    for(int ii=0; ii< Nbinsx; ii++){
      double xi=vhist[iu]->GetBinContent(ii+1);
      double ximean=hIncv->GetBinContent(ii+1);
      for(int jj=0; jj< Nbinsx; jj++){
        double xj=vhist[iu]->GetBinContent(jj+1);
        double xjmean=hIncv->GetBinContent(jj+1);
        covmx[ii][jj]+=(xi-ximean)*(xj-xjmean) ;
      }
    }   
  }//end of univ
  
  for(int ii=0; ii< Nbinsx; ii++){
    for(int jj=0; jj< Nbinsx; jj++){     
      covmx[ii][jj]/=(N-1.0);
    }
  }  
  return covmx;
}

TMatrixD calc_fer_matrix	(	TMatrixD &	mcov,
		TH1D *	hIncv
	)

Definition at line 362 of file Calculate_covariances.C.

Referenced by MakeCovSymSyst().

                                                     {
  TMatrixD fermx(mcov.GetNrows(),mcov.GetNcols());
  for(int ii=0; ii<mcov.GetNrows(); ii++){
    for(int jj=0; jj<mcov.GetNcols(); jj++){
      if(hIncv->GetBinContent(ii+1)!=0 && hIncv->GetBinContent(jj+1)!=0){
        fermx[ii][jj] = mcov[ii][jj] / ( (hIncv->GetBinContent(ii+1)) * (hIncv->GetBinContent(jj+1)) );
      }
    }
  }
  return fermx;
}

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().

                            {
  std::cout<<"Usage: Calculate_covariances [file name In (2D)] [file name In (1D)] [file name Out] [mode] [comp]"<<std::endl;    
  std::cout<<"mode = total, shape or notppfx"<<std::endl;    
  std::cout<<"comp = all"<<std::endl;
  std::cout<<"       GENIE"<<std::endl;
  std::cout<<"       PPFX"<<std::endl;
  for(int ii=0;ii<Nsymsyst;ii++){
    std::cout<<"       "<<symsyst[ii]<<std::endl;
  }
  for(int ii=0;ii<Nasymsyst;ii++){
    std::cout<<"       "<<asymsyst[ii]<<std::endl;
  }
  exit (1);
}

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.

                                                                                                                                     {
  
  //mode:
  std::vector<int> ioption = starting(fnameIn, fnameonly1DIn,fnameOut, cmode,ccomp);
  int imode = ioption[0];
  int icomp = ioption[1];
  std::cout<<"=> check: (imode,icomp) = ("<<imode<<", "<<icomp<<")"<<std::endl;
  bool is_shape = false;
  if(imode==1)is_shape = true;
  
  TFile* fOut = new TFile(fnameOut,"recreate");
  
  //Declarations:
  TFile* file1D = new TFile(fnameonly1DIn,"read"); 
  TFile* file2D = new TFile(fnameIn,"read");   
  //ENu
  std::vector<TH1D*> h1cv_ENu, h1cv_ENu_input;
  std::vector<TH1D*> h1ppfx_ENu, h1genie_ENu, h1ppfx_ENu_input, h1genie_ENu_input;
  std::vector<TH1D*> h1bfoc_ENu, h1dsyst_ENu, h1bfoc_ENu_input, h1dsyst_ENu_input;
  //Q2
  std::vector<TH1D*> h1cv_Q2, h1cv_Q2_input;
  std::vector<TH1D*> h1ppfx_Q2, h1genie_Q2, h1ppfx_Q2_input, h1genie_Q2_input;
  std::vector<TH1D*> h1bfoc_Q2, h1dsyst_Q2, h1bfoc_Q2_input, h1dsyst_Q2_input;
  //MuKin
  std::vector<TH2D*> h2cv_MuKin_input;
  std::vector<TH1D*> h1cv_MuKin;
  std::vector<TH2D*> h2ppfx_MuKin_input, h2genie_MuKin_input;   
  std::vector<TH1D*> h1ppfx_MuKin, h1genie_MuKin;
  std::vector<TH2D*> h2bfoc_MuKin_input, h2dsyst_MuKin_input;
  std::vector<TH1D*> h1bfoc_MuKin, h1dsyst_MuKin;
  
  //Getting
  //CV
  double integral_cv_ENu, integral_cv_Q2, integral_cv_MuKin;
  double integral_syst_ENu, integral_syst_Q2, integral_syst_MuKin;
  double integral_aux;
  for(unsigned i=0;i<Ncv;i++){
    h1cv_ENu_input.push_back(   (TH1D*)file1D->Get(Form("%s/xsec_%s_ENu"   ,cv[i].c_str(),cv[i].c_str())) );
    h1cv_Q2_input.push_back(    (TH1D*)file1D->Get(Form("%s/xsec_%s_Q2"    ,cv[i].c_str(),cv[i].c_str())) );
    h2cv_MuKin_input.push_back( (TH2D*)file2D->Get(Form("%s/xsec_%s_2DProj",cv[i].c_str(),cv[i].c_str())) );
    if(i==0){
      h1cv_ENu.push_back(         GetWanted1DFromENu(   h1cv_ENu_input[i],   integral_cv_ENu,  Form("henu_%s",  cv[i].c_str())) );
      h1cv_Q2.push_back(          GetWanted1DFromQ2(    h1cv_Q2_input[i] ,   integral_cv_Q2,   Form("hq2_%s",   cv[i].c_str())) );      
      h1cv_MuKin.push_back(       GetWanted1DFromMuKin( h2cv_MuKin_input[i], integral_cv_MuKin,Form("hmukin_%s",cv[i].c_str())) );   
    }   
    else{ //we currently do not use the cv for "nom_good_seventh" and "neutron_nom"
      h1cv_ENu.push_back(         GetWanted1DFromENu(   h1cv_ENu_input[i],   integral_aux,  Form("henu_%s",  cv[i].c_str())) );
      h1cv_Q2.push_back(          GetWanted1DFromQ2(    h1cv_Q2_input[i] ,   integral_aux,  Form("hq2_%s",   cv[i].c_str())) );      
      h1cv_MuKin.push_back(       GetWanted1DFromMuKin( h2cv_MuKin_input[i], integral_aux,  Form("hmukin_%s",cv[i].c_str())) );
    }
  }
  
  //PPFX
  for(unsigned i=0;i<Nppfx;i++){
    h1ppfx_ENu_input.push_back(   (TH1D*)file1D->Get(Form("ppfx%d/xsec_ppfx%d_ENu"    ,i,i)));
    h1ppfx_Q2_input.push_back(    (TH1D*)file1D->Get(Form("ppfx%d/xsec_ppfx%d_Q2"     ,i,i)));
    h2ppfx_MuKin_input.push_back( (TH2D*)file2D->Get(Form("ppfx%d/xsec_ppfx%d_2DProj" ,i,i)));
    h1ppfx_ENu.push_back(         GetWanted1DFromENu(  h1ppfx_ENu_input[i],  integral_syst_ENu,   Form("henu_ppfx%d",  i)) );
    h1ppfx_Q2.push_back(          GetWanted1DFromQ2(   h1ppfx_Q2_input[i],   integral_syst_Q2,    Form("hq2_ppfx%d" ,  i)) );
    h1ppfx_MuKin.push_back(       GetWanted1DFromMuKin(h2ppfx_MuKin_input[i],integral_syst_MuKin, Form("hmukin_ppfx%d",i)) );
    if(is_shape){
      scaling(h1ppfx_ENu[i]  ,integral_cv_ENu  / integral_syst_ENu);
      scaling(h1ppfx_Q2[i]   ,integral_cv_Q2   / integral_syst_Q2);
      scaling(h1ppfx_MuKin[i],integral_cv_MuKin/ integral_syst_MuKin);
    }
  }
  h1ppfx_ENu_input.clear();
  h1ppfx_Q2_input.clear();
  h2ppfx_MuKin_input.clear();

  //GENIE
  for(unsigned i=0;i<Ngenie;i++){
    h1genie_ENu_input.push_back(  (TH1D*)file1D->Get(Form("genie%d/xsec_genie%d_ENu"   ,i,i)));
    h1genie_Q2_input.push_back(   (TH1D*)file1D->Get(Form("genie%d/xsec_genie%d_Q2"    ,i,i)));
    h2genie_MuKin_input.push_back((TH2D*)file2D->Get(Form("genie%d/xsec_genie%d_2DProj",i,i)));
    h1genie_ENu.push_back(         GetWanted1DFromENu(  h1genie_ENu_input[i],  integral_syst_ENu,   Form("henu_genie%d",  i)) );
    h1genie_Q2.push_back(          GetWanted1DFromQ2(   h1genie_Q2_input[i],   integral_syst_Q2,    Form("hq2_genie%d",   i)) );
    h1genie_MuKin.push_back(       GetWanted1DFromMuKin(h2genie_MuKin_input[i],integral_syst_MuKin, Form("hmukin_genie%d",i)) );
    if(is_shape){
      scaling(h1genie_ENu[i]  ,integral_cv_ENu  / integral_syst_ENu);
      scaling(h1genie_Q2[i]   ,integral_cv_Q2   / integral_syst_Q2);
      scaling(h1genie_MuKin[i],integral_cv_MuKin/ integral_syst_MuKin);
    }
  }
  h1genie_ENu_input.clear();
  h1genie_Q2_input.clear();
  h2genie_MuKin_input.clear();
  
  //BEAM FOC and DET SYST
  for(unsigned i=0;i<Nbfoc;i++){
    h1bfoc_ENu_input.push_back(   (TH1D*)file1D->Get(Form("%s/xsec_%s_ENu"   ,bfoc[i].c_str(),bfoc[i].c_str())) );
    h1bfoc_Q2_input.push_back(    (TH1D*)file1D->Get(Form("%s/xsec_%s_Q2"    ,bfoc[i].c_str(),bfoc[i].c_str())) );
    h2bfoc_MuKin_input.push_back( (TH2D*)file2D->Get(Form("%s/xsec_%s_2DProj",bfoc[i].c_str(),bfoc[i].c_str())) );
  }
  for(unsigned i=0;i<Ndsyst;i++){
    h1dsyst_ENu_input.push_back(   (TH1D*)file1D->Get(Form("%s/xsec_%s_ENu"   ,dsyst[i].c_str(),dsyst[i].c_str())) );
    h1dsyst_Q2_input.push_back(    (TH1D*)file1D->Get(Form("%s/xsec_%s_Q2"    ,dsyst[i].c_str(),dsyst[i].c_str())) );
    h2dsyst_MuKin_input.push_back( (TH2D*)file2D->Get(Form("%s/xsec_%s_2DProj",dsyst[i].c_str(),dsyst[i].c_str())) );
  }

  //Getting systematics histograms:
  //Symmetricals:
  std::vector<TH1D*> h1symsyst_ENu_input, h1symsyst_Q2_input;
  std::vector<TH2D*> h2symsyst_MuKin_input;
  std::vector<TH1D*> h1symsyst_ENu,       h1symsyst_Q2,       h1symsyst_MuKin;
  //Asymmetricals:
  //one side syst cal shape
  std::vector<TH1D*> h1asymsyst_ENu_input, h1asymsyst_Q2_input;
  std::vector<TH2D*> h2asymsyst_MuKin_input;
  std::vector<TH1D*> h1asymsyst_ENu,       h1asymsyst_Q2,       h1asymsyst_MuKin;
  for(unsigned i=0;i<Nsymsyst;i++){
    //Foc: 11 systematics
    if(i<11){   
      int ifoc = 2*i;
      h1symsyst_ENu_input.push_back(   GetSymSyst(h1bfoc_ENu_input[ifoc]  ,h1bfoc_ENu_input[ifoc+1]  ,h1cv_ENu_input[0]  ,symsyst[i], "enu")   );
      h1symsyst_Q2_input.push_back(    GetSymSyst(h1bfoc_Q2_input[ifoc]   ,h1bfoc_Q2_input[ifoc+1]   ,h1cv_Q2_input[0]   ,symsyst[i], "q2")    );
      h2symsyst_MuKin_input.push_back( GetSymSyst(h2bfoc_MuKin_input[ifoc],h2bfoc_MuKin_input[ifoc+1],h2cv_MuKin_input[0],symsyst[i], "mukin") );
    }
    //Det sym syst: 5 systematics
    else{
      int idsyst = 2*(i-11);
      h1symsyst_ENu_input.push_back(   GetSymSyst(h1dsyst_ENu_input[idsyst]  ,h1dsyst_ENu_input[idsyst+1]  ,h1cv_ENu_input[0]  ,symsyst[i], "enu")   );
      h1symsyst_Q2_input.push_back(    GetSymSyst(h1dsyst_Q2_input[idsyst]   ,h1dsyst_Q2_input[idsyst+1]   ,h1cv_Q2_input[0]   ,symsyst[i], "q2")    );
      h2symsyst_MuKin_input.push_back( GetSymSyst(h2dsyst_MuKin_input[idsyst],h2dsyst_MuKin_input[idsyst+1],h2cv_MuKin_input[0],symsyst[i], "mukin") );
    }
    h1symsyst_ENu.push_back(   GetWanted1DFromENu(   h1symsyst_ENu_input[i],   integral_syst_ENu,   Form("henu_%s"   ,symsyst[i].c_str())) );
    h1symsyst_Q2.push_back(    GetWanted1DFromQ2(    h1symsyst_Q2_input[i],    integral_syst_Q2,    Form("hq2_%s"    ,symsyst[i].c_str())) );
    h1symsyst_MuKin.push_back( GetWanted1DFromMuKin( h2symsyst_MuKin_input[i], integral_syst_MuKin, Form("hmukin_%s" ,symsyst[i].c_str())) );
    if(is_shape){
      scaling(h1symsyst_ENu[i]  ,integral_cv_ENu  / integral_syst_ENu);
      scaling(h1symsyst_Q2[i]   ,integral_cv_Q2   / integral_syst_Q2);
      scaling(h1symsyst_MuKin[i],integral_cv_MuKin/ integral_syst_MuKin);
    }
  }
  //Det sym syst: 2 systematics
  for(unsigned i=0;i<Nasymsyst;i++){
    h1asymsyst_ENu_input.push_back(  h1dsyst_ENu_input[10+i]);
    h1asymsyst_Q2_input.push_back(   h1dsyst_Q2_input[10+i]);
    h2asymsyst_MuKin_input.push_back(h2dsyst_MuKin_input[10+i]);
    h1asymsyst_ENu.push_back(   GetWanted1DFromENu(   h1asymsyst_ENu_input[i],  integral_syst_ENu,   Form("henu_%s"   ,asymsyst[i].c_str())) );
    h1asymsyst_Q2.push_back(    GetWanted1DFromQ2(    h1asymsyst_Q2_input[i],   integral_syst_Q2,    Form("hq2_%s"    ,asymsyst[i].c_str())) );
    h1asymsyst_MuKin.push_back( GetWanted1DFromMuKin( h2asymsyst_MuKin_input[i],integral_syst_MuKin, Form("hmukin_%s" ,asymsyst[i].c_str())) );
    if(is_shape){
      scaling(h1asymsyst_ENu[i]  ,integral_cv_ENu  / integral_syst_ENu);
      scaling(h1asymsyst_Q2[i]   ,integral_cv_Q2   / integral_syst_Q2);
      scaling(h1asymsyst_MuKin[i],integral_cv_MuKin/ integral_syst_MuKin);
    }
  }
   
  //Making universes:
  fOut->mkdir("check");
  //Total
  if(icomp==0){
    MakeCov("enu"  , fOut, h1cv_ENu[0]  , h1genie_ENu  , h1ppfx_ENu  , h1symsyst_ENu  , h1asymsyst_ENu  , imode);
    MakeCov("q2"   , fOut, h1cv_Q2[0]   , h1genie_Q2   , h1ppfx_Q2   , h1symsyst_Q2   , h1asymsyst_Q2   , imode);
    MakeCov("mukin", fOut, h1cv_MuKin[0], h1genie_MuKin, h1ppfx_MuKin, h1symsyst_MuKin, h1asymsyst_MuKin, imode);
  }
  //Only GENIE and PPFX
  if(icomp==1){
    MakeCovUniverses("enu_genie"  , fOut, h1cv_ENu[0]  , h1genie_ENu  , imode);
    MakeCovUniverses("q2_gene"    , fOut, h1cv_Q2[0]   , h1genie_Q2   , imode);
    MakeCovUniverses("mukin_genie", fOut, h1cv_MuKin[0], h1genie_MuKin, imode);
  }
  if(icomp==2){
    MakeCovUniverses("enu_ppfx"   , fOut, h1cv_ENu[0]  , h1ppfx_ENu   , imode);
    MakeCovUniverses("q2_ppfx"    , fOut, h1cv_Q2[0]   , h1ppfx_Q2    , imode);
    MakeCovUniverses("mukin_ppfx" , fOut, h1cv_MuKin[0], h1ppfx_MuKin , imode);
  }
  //Only symmetric systematics:
  int local_idx = -1;
  if(icomp>=3 && icomp<=18){
    local_idx = icomp -3;
    MakeCovSymSyst("enu",local_idx,fOut,h1cv_ENu[0]  ,h1symsyst_ENu[local_idx]  ,imode);
    MakeCovSymSyst("q2",local_idx,fOut,h1cv_Q2[0]   ,h1symsyst_Q2[local_idx]   ,imode);
    MakeCovSymSyst("mukin",local_idx,fOut,h1cv_MuKin[0],h1symsyst_MuKin[local_idx],imode);
  }
  //Only asymmetric systematics:
  if(icomp>=19 && icomp<=20){
    local_idx = icomp -19;
    MakeCovASymSyst("enu",local_idx,fOut,h1cv_ENu[0]  ,h1asymsyst_ENu[local_idx]  ,imode);
    MakeCovASymSyst("q2",local_idx,fOut,h1cv_Q2[0]   ,h1asymsyst_Q2[local_idx]   ,imode);
    MakeCovASymSyst("mukin",local_idx,fOut,h1cv_MuKin[0],h1asymsyst_MuKin[local_idx],imode);
  }
  
 
  fOut->Close();
  std::cout<<"Ending "<<std::endl;
}

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().

                                                                                                  {  
  TH1D* houtput = (TH1D*)hIncv->Clone(Form("hSymSist%s_%s",SystName.c_str(),VarName.c_str() ));
  int Nbins = hIncv->GetXaxis()->GetNbins();
  for(int ii=1;ii<=Nbins;ii++){
    double pp = hInshiftp->GetBinContent(ii);
    double nn = hInshiftn->GetBinContent(ii);
    double cv = hIncv->GetBinContent(ii);
    double diff_pp = std::abs(pp-cv);
    double diff_nn = std::abs(nn-cv);
    if(diff_pp>diff_nn)houtput->SetBinContent(ii,pp);
    else               houtput->SetBinContent(ii,2.*cv-nn);
  }
  return houtput;
}

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.

                                                                                                  {  
  TH2D* houtput = (TH2D*)hIncv->Clone(Form("hSymSist%s_%s",SystName.c_str(),VarName.c_str() ));
  int NbinsX = hIncv->GetXaxis()->GetNbins();
  int NbinsY = hIncv->GetYaxis()->GetNbins();
  for(int ii=1;ii<=NbinsX;ii++){
    for(int jj=1;jj<=NbinsY;jj++){
      double pp = hInshiftp->GetBinContent(ii,jj);
      double nn = hInshiftn->GetBinContent(ii,jj);
      double cv = hIncv->GetBinContent(ii,jj);
      double diff_pp = std::abs(pp-cv);
      double diff_nn = std::abs(nn-cv);
      if(diff_pp>diff_nn)houtput->SetBinContent(ii,jj,pp);
      else               houtput->SetBinContent(ii,jj,2.*cv-nn);
    }
  }
  return houtput;
}

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* h1Wanted_ENu = new TH1D(name.c_str(),"",Nwanted_ENuBins,0,Nwanted_ENuBins);
  int counter = 1;
  int Nbinsx  = h1In->GetXaxis()->GetNbins();
  IntegralVal = 0;
  for(int ii=1;ii<=Nbinsx;ii++){
    if(ii<=3 || ii==Nbinsx)continue;
    IntegralVal += ( h1In->GetBinContent(ii) ) * ( h1In->GetXaxis()->GetBinWidth(ii) );
    h1Wanted_ENu->SetBinContent(counter,h1In->GetBinContent(ii));
    h1Wanted_ENu->SetBinError( counter,0);
    counter++;
  }
  return h1Wanted_ENu;
}

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* h1Wanted_MuKin = new TH1D(name.c_str(),"",Nwanted_MuKinBins,0,Nwanted_MuKinBins);
  int counter = 1;
  int NbinX = h2In->GetXaxis()->GetNbins();
  int NbinY = h2In->GetYaxis()->GetNbins();
  IntegralVal = 0;
  for(int ii=1;ii<=NbinX;ii++){
    if(ii==1)continue;
    for(int jj=1;jj<=NbinY;jj++){     
      if(jj==1 || jj==NbinY)continue;
      if( (jj-1) > MuTperCosBins[ii-2] )continue;
      IntegralVal += ( h2In->GetBinContent(ii,jj) ) * ( h2In->GetXaxis()->GetBinWidth(ii) ) * ( h2In->GetYaxis()->GetBinWidth(jj) );
      h1Wanted_MuKin->SetBinContent(counter,h2In->GetBinContent(ii,jj));
      h1Wanted_MuKin->SetBinError( counter,0);
      counter++;
    }
  }
  return h1Wanted_MuKin;
}

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().

                                                                        {
  TH1D* h1Wanted_Q2 = new TH1D(name.c_str(),"",Nwanted_Q2Bins,0,Nwanted_Q2Bins);
  int counter = 1;
  int Nbinsx  = h1In->GetXaxis()->GetNbins();
  IntegralVal = 0;
  for(int ii=1;ii<=Nbinsx;ii++){
    if(ii==Nbinsx)continue;
    IntegralVal += ( h1In->GetBinContent(ii) ) * ( h1In->GetXaxis()->GetBinWidth(ii) );
    h1Wanted_Q2->SetBinContent(counter,h1In->GetBinContent(ii));
    h1Wanted_Q2->SetBinError( counter,0);
    counter++;
  }
  return h1Wanted_Q2;
}

int main	(	int	argc,
		const char *	argv[]
	)

Definition at line 689 of file Calculate_covariances.C.

References Calculate_covariances().

                                      {
  if(argc==1)Calculate_covariances();
  Calculate_covariances(argv[1],argv[2],argv[3],argv[4],argv[5]);
  return 0;
}

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().

                                                                                                                                                                                 {
  std::cout<<"=> Making covariances"<<std::endl;
  //General:
  std::vector<TH1D*> vout; 
  int univ_id = -1;
  double this_shift = 0;
  
  for(int ii=0;ii<Ngenie;ii++){  
    for(int jj=0;jj<Nppfx;jj++){   
      univ_id++;
      if(univ_id%5000==0)std::cout<<"=> Shifting var: "<<univ_id<<" from "<<(Nppfx*Ngenie)<<std::endl;      
      vout.push_back( (TH1D*)hIncv->Clone(Form("h%s_u%d",histName.c_str(),univ_id)) );
      //genie
      vout[univ_id]->Add(hIngenie[ii]);
      vout[univ_id]->Add(hIncv,-1);
      //ppfx
      if(thismode == 0 || thismode == 1){
        if(ii==0 && jj==0)std::cout<<"=> Doing ppfx here"<<std::endl;
        vout[univ_id]->Add(hInppfx[jj]);
        vout[univ_id]->Add(hIncv,-1);           
      }
      //sym syst
      for(unsigned kk=0;kk<Nsymsyst;kk++){
        r3->SetSeed(basesymsyst_rdm[kk]+univ_id);
        this_shift = r3->Gaus(0.0,1.0);
        AddShift(vout[univ_id],this_shift,hIncv, hInsym[kk]);
      }
      //one side syst cal shape:
      for(unsigned kk=0;kk<Nasymsyst;kk++){
        r3->SetSeed(baseasymsyst_rdm[kk]+univ_id);
        this_shift = std::abs(r3->Gaus(0.0,1.0));
        AddOneSideShift(vout[univ_id],this_shift,hIncv, hInasym[kk]);
      }
    }
  } //end of new universes
  
  std::cout<<"We generated  "<<vout.size()<< " of new universes for "<<histName<<std::endl;
  TMatrixD cov = calc_cov_matrix(vout,hIncv);
  TMatrixD cor = calc_corr_matrix(cov);
  int nofbins  = hIncv->GetXaxis()->GetNbins();
  
  TH2D*  hcov = new TH2D(Form("hcov_%s",histName.c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  TH2D*  hcor = new TH2D(Form("hcor_%s",histName.c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  for(int ii=1;ii<=nofbins;ii++){
    for(int jj=1;jj<=nofbins;jj++){
      hcov->SetBinContent(ii,jj,cov[ii-1][jj-1]);
      hcor->SetBinContent(ii,jj,cor[ii-1][jj-1]);
    }
  }
  std::cout<<"Storing results for "<<histName<<std::endl;
  fIn->cd(); 
  hIncv->Write(Form("hcv_%s",histName.c_str()));
  hcov->Write( Form("cov_%s",histName.c_str()));
  hcor->Write( Form("cor_%s",histName.c_str()));

  //storing sym systematics
  fIn->mkdir(("check_"+histName).c_str()); 
  fIn->cd(("check_"+histName).c_str()); 
  for(int ii=0;ii<hInsym.size();ii++ )hInsym[ii]->Write();
  for(int ii=0;ii<hInasym.size();ii++)hInasym[ii]->Write();
  
  //Clearing vectors:
  hIngenie.clear();
  hInppfx.clear();
  hInsym.clear();
  hInasym.clear();
  vout.clear();
  std::vector<TH1D*>().swap(vout);
  vout.shrink_to_fit();
}

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().

                                                                                                       {
  if(thismode==2)return; //this calculation is not for "not ppfx mode"
  std::cout<<"=> Making covariances asymmetric systematics: "<< varName<<", "<<asymsyst[iSyst] <<", mode: "<< thismode <<std::endl;

  //General:
  std::vector<TH1D*> vout; 
  int univ_id = -1;
  double this_shift = 0;
  for(int ii=0;ii<Nuniv;ii++){  
    univ_id++;
    vout.push_back( (TH1D*)hIncv->Clone(Form("h%s_%s_u%d",varName.c_str(),asymsyst[iSyst].c_str(),univ_id)) );    
    r3->SetSeed(baseasymsyst_rdm[iSyst]+univ_id);
    this_shift = std::abs(r3->Gaus(0.0,1.0));
    AddOneSideShift(vout[ii],this_shift,hIncv, hInasym);
  } //end of new universes
  
  std::cout<<"We generated  "<<vout.size()<< " of new universes for "<<varName<<", "<<asymsyst[iSyst]<<std::endl;
  TMatrixD cov = calc_cov_matrix(vout,hIncv);
  TMatrixD cor = calc_corr_matrix(cov);
  int nofbins  = hIncv->GetXaxis()->GetNbins();
  
  TH2D*  hcov = new TH2D(Form("hcov_%s_%s",varName.c_str(),asymsyst[iSyst].c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  TH2D*  hcor = new TH2D(Form("hcor_%s_%s",varName.c_str(),asymsyst[iSyst].c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  for(int ii=1;ii<=nofbins;ii++){
    for(int jj=1;jj<=nofbins;jj++){
      hcov->SetBinContent(ii,jj,cov[ii-1][jj-1]);
      hcor->SetBinContent(ii,jj,cor[ii-1][jj-1]);
    }
  }
  std::cout<<"Storing results for "<<asymsyst[iSyst]<<std::endl;
  fIn->cd(); 
  hcov->Write( Form("cov_%s_%s",varName.c_str(),asymsyst[iSyst].c_str()));
  hcor->Write( Form("cor_%s_%s",varName.c_str(),asymsyst[iSyst].c_str()));
}

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().

                                                                                                       {
  if(thismode==2)return; //this calculation is not for the "not ppfx mode"
  std::cout<<"=> Making covariances symmetric systematics: "<< varName<<", "<<symsyst[iSyst]<<", mode: "<< thismode <<std::endl;

  //General:
  std::vector<TH1D*> vout; 
  int univ_id = -1;
  double this_shift = 0;

  for(int ii=0;ii<Nuniv;ii++){  
    univ_id++;
    vout.push_back( (TH1D*)hIncv->Clone(Form("h%s_%s_u%d",varName.c_str(),symsyst[iSyst].c_str(),univ_id)) );    
    r3->SetSeed(basesymsyst_rdm[iSyst]+univ_id);
    this_shift = r3->Gaus(0.0,1.0);
    AddShift(vout[ii],this_shift,hIncv, hInsym);
  } //end of new universes
  
  std::cout<<"We generated  "<<vout.size()<< " of new universes for "<<varName<<", "<<symsyst[iSyst]<<std::endl;
  TMatrixD cov = calc_cov_matrix(vout,hIncv);
  TMatrixD cor = calc_corr_matrix(cov);
  TMatrixD fer = calc_fer_matrix(cov,hIncv);
  int nofbins  = hIncv->GetXaxis()->GetNbins();
  
  TH2D*  hcov = new TH2D(Form("hcov_%s_%s",varName.c_str(),symsyst[iSyst].c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  TH2D*  hcor = new TH2D(Form("hcor_%s_%s",varName.c_str(),symsyst[iSyst].c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  TH2D*  hfer = new TH2D(Form("hfer_%s_%s",varName.c_str(),symsyst[iSyst].c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  for(int ii=1;ii<=nofbins;ii++){
    for(int jj=1;jj<=nofbins;jj++){
      hcov->SetBinContent(ii,jj,cov[ii-1][jj-1]);
      hcor->SetBinContent(ii,jj,cor[ii-1][jj-1]);
      hfer->SetBinContent(ii,jj,fer[ii-1][jj-1]);
    }
  }
  std::cout<<"Storing results"<<std::endl;
  fIn->cd(); 
  hcov->Write( Form("cov_%s_%s",varName.c_str(),symsyst[iSyst].c_str()));
  hcor->Write( Form("cor_%s_%s",varName.c_str(),symsyst[iSyst].c_str()));
  fIn->cd("check");
  hInsym->Write();
  hIncv->Write();
  hfer->Write();
}

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().

                                                                                                           {
  if(thismode==2)return; //this calculation is not for the "not ppfx mode"
  std::cout<<"=> Making covariances universes: "<< histName<<", mode: "<< thismode <<std::endl;

  std::vector<TH1D*> vout; 
  int univ_id = -1;
  int NNuniv = int(hInuniv.size());
  for(int ii=0;ii<NNuniv;ii++){  
    univ_id++;
    vout.push_back( (TH1D*)hIncv->Clone(Form("h%s_u%d",histName.c_str(),univ_id)) );
    vout[ii]->Add(hInuniv[ii]);
    vout[ii]->Add(hIncv,-1);
  } //end of new universes
  
  std::cout<<"We generated  "<<vout.size()<< " of new universes for "<<histName<<std::endl;
  TMatrixD cov = calc_cov_matrix(vout,hIncv);
  TMatrixD cor = calc_corr_matrix(cov);
  int nofbins  = hIncv->GetXaxis()->GetNbins();
  
  TH2D* hcov = new TH2D(Form("hcov_%s",histName.c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  TH2D* hcor = new TH2D(Form("hcor_%s",histName.c_str()),"",nofbins,0,nofbins,nofbins,0,nofbins);
  for(int ii=1;ii<=nofbins;ii++){
    for(int jj=1;jj<=nofbins;jj++){
      hcov->SetBinContent(ii,jj,cov[ii-1][jj-1]);
      hcor->SetBinContent(ii,jj,cor[ii-1][jj-1]);
    }
  }
  std::cout<<"Storing results for "<<histName<<std::endl;
  fIn->cd(); 
  hcov->Write( Form("cov_%s",histName.c_str()));
  hcor->Write( Form("cor_%s",histName.c_str()));
}

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().

                                                 {
  hIn->Scale(shape_scale);
}

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().

                                                                                                                                  {
  
  std::vector<int> voption(2,-1);
  std::cout<<"Starting Calculate_covariances()"<<std::endl;
  std::cout<<"Input:"<<std::endl;
  std::cout<<"file name In (2D): "<< fnameIn      <<std::endl;
  std::cout<<"file name In (1D): "<< fnameonly1DIn<<std::endl;
  std::cout<<"file name Out    : "<< fnameOut     <<std::endl;
  std::cout<<"mode             : "<< cmode        <<std::endl;
  std::cout<<"component        : "<< ccomp        <<std::endl;
  
  //mode
  if(std::string(cmode)=="total"){
    std::cout<<"Calculating total covariance"<<std::endl;
    voption[0] = 0;
  }
  else if(std::string(cmode)=="shape"){
    std::cout<<"Calculating shape covariance"<<std::endl;
    voption[0] = 1;
  }
  else if(std::string(cmode)=="notppfx"){
    std::cout<<"Calculating covariance without ppfx"<<std::endl;
    voption[0] = 2;
  }
  else{
    Calculate_covariances();
  }
  //component
  if(std::string(ccomp)=="all"){
    std::cout<<"Calculating all components"<<std::endl;
    voption[1] = 0;
  }
  if(std::string(ccomp)=="GENIE"){
    std::cout<<"Calculating "<<ccomp<<" systematic"<<std::endl;
    voption[1] = 1;
  }
  if(std::string(ccomp)=="PPFX"){
    std::cout<<"Calculating "<<ccomp<<" systematic"<<std::endl;
    voption[1] = 2;
  }
  for(int ii=0;ii<Nsymsyst;ii++){
    if(std::string(ccomp) == symsyst[ii]){
      std::cout<<"Calculating "<<ccomp<<" systematic"<<std::endl;
      voption[1] = 3 + ii;
    }
  }
  for(int ii=0;ii<Nasymsyst;ii++){
    if(std::string(ccomp) == asymsyst[ii]){
      std::cout<<"Calculating "<<ccomp<<" systematic"<<std::endl;
      voption[1] = 19 + ii;
    }
  }
  if(voption[1] == -1)Calculate_covariances();
  
  return voption;
}

Variable Documentation

const std::string asymsyst[Nasymsyst] = {"CalibShape","Cherenkov"}

Definition at line 63 of file Calculate_covariances.C.

Referenced by Calculate_covariances(), MakeCovASymSyst(), and starting().

const int baseasymsyst_rdm[Nasymsyst] = {26000000,27000000}

Definition at line 64 of file Calculate_covariances.C.

Referenced by MakeCov(), and MakeCovASymSyst().

const int basesymsyst_rdm[Nsymsyst]

Initial value:

= {10000000,11000000,12000000,13000000,14000000,
                                       15000000,16000000,17000000,18000000,19000000,
                                       20000000,21000000,22000000,23000000,24000000,
                                       25000000}

Definition at line 56 of file Calculate_covariances.C.

Referenced by MakeCov(), and MakeCovSymSyst().

const std::string bfoc[Nbfoc]

Initial value:

= {"2kA_Up","2kA_Dw",
                                 "02mmBeamSpotSize_Up","02mmBeamSpotSize_Dw",
                                 "1mmBeamShiftX_Up","1mmBeamShiftX_Dw",
                                 "1mmBeamShiftY_Up","1mmBeamShiftY_Dw",
                                 "3mmHorn1X_Up","3mmHorn1X_Dw",
                                 "3mmHorn1Y_Up","3mmHorn1Y_Dw",
                                 "3mmHorn2X_Up","3mmHorn2X_Dw",
                                 "3mmHorn2Y_Up","3mmHorn2Y_Dw",
                                 "7mmTargetZ_Up","7mmTargetZ_Dw",
                                 "MagneticFieldinDecayPipe_Up","MagneticFieldinDecayPipe_Dw",
                                 "1mmHornWater_Up","1mmHornWater_Dw"}

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]

Initial value:

= {"angleup","angledw",
                                   "calibpos","calibneg",
                                   "lightup","lightdown",
                                   "muoneup","muonedw",
                                   "neutron_Up","neutron_Dw",
                                   "calibshift",
                                   "ckv"}

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().

const int Nuniv = Ngenie*Nppfx

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]

Initial value:

= {"FocCurrMag","FocBeamSpotSize","FocBeamShiftX","FocBeamShiftY",
                                       "FocHorn1X", "FocHorn1Y", "FocHorn2X", "FocHorn2Y", 
                                       "FocTargetZ", "FocBFieldDPIP", "FocHornWaterIC",
                                       "MuAngle", "CalibShift", "light", "MuE","Neu"}

Definition at line 51 of file Calculate_covariances.C.

Referenced by Calculate_covariances(), MakeCovSymSyst(), and starting().