Go to the source code of this file.

Namespaces
	makeCovariance2

Functions
def	makeCovariance2.MakeErrorHistogram (universes, asFrac=True)

def	makeCovariance2.MakeCovariance (universes, asFrac=True)

def	makeCovariance2.MergeHistograms (cvs, universes, name)

def	makeCovariance2.GetCorrelation (cov_matrix)

Variables
string	makeCovariance2.beam_config = "reference"

list	makeCovariance2.plots = ["numu_neutrino_ND","numubar_neutrino_ND","nue_neutrino_ND","nuebar_neutrino_ND","numu_antineutrino_ND","numubar_antineutrino_ND","nue_antineutrino_ND","nuebar_antineutrino_ND","numu_neutrino_FD","numubar_neutrino_FD","nue_neutrino_FD","nuebar_neutrino_FD","numu_antineutrino_FD","numubar_antineutrino_FD","nue_antineutrino_FD","nuebar_antineutrino_FD"]

list	makeCovariance2.start_indices = [0,19,38,45,52,71,90,97,104,123,142,149,156,175,194,201]

list	makeCovariance2.in_files = []

list	makeCovariance2.cv_histos = []

int	makeCovariance2.n_universes = 100

list	makeCovariance2.error_universes = []

list	makeCovariance2.error_types = ["total","others","targetpcpi","targetpck","targetncpi","targetpcnu","targetmesonInc","targetnuA","tgtabsorption","totabsorption"]

list	makeCovariance2.error_types2 = ["total","others","targetpcpi","targetpcK","targetncpi","targetpcnu","targetmesonInc","targetnucleonA","targetabsorption","absorption"]

list	makeCovariance2.error_titles = ["Total","Other","pC#rightarrow#pi","pC#rightarrowK","nC#rightarrow#pi","pc#rightarrowN","Meson Inc","NucleonA","Target Absorption","Other Absorption"]

list	makeCovariance2.numu_bins = [0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,7 ,8,12,16,20,40,100]

list	makeCovariance2.flavors = ["neutrino_near_numu","neutrino_near_numubar","neutrino_near_nue","neutrino_near_nuebar","antineutrino_near_numu","antineutrino_near_numubar","antineutrino_near_nue","antineutrino_near_nuebar","neutrino_far_numu","neutrino_far_numubar","neutrino_far_nue","neutrino_far_nuebar","antineutrino_far_numu","antineutrino_far_numubar","antineutrino_far_nue","antineutrino_far_nuebar"]

list	makeCovariance2.merged_histograms = []

list	makeCovariance2.merged_universes = []

list	makeCovariance2.flux_uncertainties = []

	makeCovariance2.c1 = ROOT.TCanvas("c1","",1000,800)

	makeCovariance2.flux_covariance = MakeCovariance(merged_universes[error_types.index("total")])

	makeCovariance2.latex0 = ROOT.TLatex( 0.5, 0.93, "Flux Hadron Production Fractional Covariance Matrix" );

	makeCovariance2.latex1 = ROOT.TLatex( 0.24, 0.035, "Near Detector" );

	makeCovariance2.latex2 = latex1.Clone()

	makeCovariance2.latex3 = latex1.Clone()

	makeCovariance2.latex4 = latex1.Clone()

	makeCovariance2.latex5 = ROOT.TLatex(0.2,0.08,"#nu Mode")

	makeCovariance2.latex6 = latex5.Clone()

	makeCovariance2.latex7 = latex5.Clone()

	makeCovariance2.latex8 = latex5.Clone()

	makeCovariance2.latex9 = ROOT.TLatex(0.085,0.21,"#nu Mode")

	makeCovariance2.latex10 = latex9.Clone()

	makeCovariance2.latex11 = latex9.Clone()

	makeCovariance2.latex12 = latex9.Clone()

	makeCovariance2.latex13 = ROOT.TLatex( 0.175, 0.124, "#nu_{#mu}" );

	makeCovariance2.latex14 = ROOT.TLatex( 0.245, 0.124, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex15 = ROOT.TLatex( 0.285, 0.124, "#nu_{e}" );

	makeCovariance2.latex16 = ROOT.TLatex( 0.31, 0.124, "#bar{#nu}_{e}" );

	makeCovariance2.latex17 = ROOT.TLatex( 0.35, 0.124, "#nu_{#mu}" );

	makeCovariance2.latex18 = ROOT.TLatex( 0.415, 0.124, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex19 = ROOT.TLatex( 0.46, 0.124, "#nu_{e}" );

	makeCovariance2.latex20 = ROOT.TLatex( 0.485, 0.124, "#bar{#nu}_{e}" );

	makeCovariance2.latex21 = ROOT.TLatex( 0.53, 0.124, "#nu_{#mu}" );

	makeCovariance2.latex22 = ROOT.TLatex( 0.595, 0.124, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex23 = ROOT.TLatex( 0.631, 0.124, "#nu_{e}" );

	makeCovariance2.latex24 = ROOT.TLatex( 0.654, 0.124, "#bar{#nu}_{e}" );

	makeCovariance2.latex25 = ROOT.TLatex( 0.70, 0.124, "#nu_{#mu}" );

	makeCovariance2.latex26 = ROOT.TLatex( 0.77, 0.124, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex27 = ROOT.TLatex( 0.81, 0.124, "#nu_{e}" );

	makeCovariance2.latex28 = ROOT.TLatex( 0.835, 0.124, "#bar{#nu}_{e}" );

	makeCovariance2.latex29 = ROOT.TLatex( 0.11,0.18, "#nu_{#mu}" );

	makeCovariance2.latex30 = ROOT.TLatex( 0.11,0.25, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex31 = ROOT.TLatex( 0.11,0.295, "#nu_{e}" );

	makeCovariance2.latex32 = ROOT.TLatex( 0.11,0.32, "#bar{#nu}_{e}" );

	makeCovariance2.latex33 = ROOT.TLatex( 0.11,0.37, "#nu_{#mu}" );

	makeCovariance2.latex34 = ROOT.TLatex( 0.11,0.435, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex35 = ROOT.TLatex( 0.11,0.48, "#nu_{e}" );

	makeCovariance2.latex36 = ROOT.TLatex( 0.11,0.505, "#bar{#nu}_{e}" );

	makeCovariance2.latex37 = ROOT.TLatex( 0.11,0.56, "#nu_{#mu}" );

	makeCovariance2.latex38 = ROOT.TLatex( 0.11,0.62, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex39 = ROOT.TLatex( 0.11,0.67, "#nu_{e}" );

	makeCovariance2.latex40 = ROOT.TLatex( 0.11,0.70, "#bar{#nu}_{e}" );

	makeCovariance2.latex41 = ROOT.TLatex( 0.11,0.75, "#nu_{#mu}" );

	makeCovariance2.latex42 = ROOT.TLatex( 0.11,0.81, "#bar{#nu}_{#mu}" );

	makeCovariance2.latex43 = ROOT.TLatex( 0.11,0.855, "#nu_{e}" );

	makeCovariance2.latex44 = ROOT.TLatex( 0.11,0.88, "#bar{#nu}_{e}" );

	makeCovariance2.line_1 = ROOT.TLine(19.5,0.5,19.5,208)

	makeCovariance2.line_2 = ROOT.TLine(38.5,0.5,38.5,208)

	makeCovariance2.line_3 = ROOT.TLine(45.5,0.5,45.5,208)

	makeCovariance2.line_4 = ROOT.TLine(52.5,0.5,52.5,208)

	makeCovariance2.line_5 = ROOT.TLine(71.5,0.5,71.5,208)

	makeCovariance2.line_6 = ROOT.TLine(90.5,0.5,90.5,208)

	makeCovariance2.line_7 = ROOT.TLine(97.5,0.5,97.5,208)

	makeCovariance2.line_8 = ROOT.TLine(104.5,0.5,104.5,208)

	makeCovariance2.line_9 = ROOT.TLine(123.5,0.5,123.5,208)

	makeCovariance2.line_10 = ROOT.TLine(142.5,0.5,142.5,208)

	makeCovariance2.line_11 = ROOT.TLine(149.5,0.5,149.5,208)

	makeCovariance2.line_12 = ROOT.TLine(156.5,0.5,156.5,208)

	makeCovariance2.line_13 = ROOT.TLine(175.5,0.5,175.5,208)

	makeCovariance2.line_14 = ROOT.TLine(194.5,0.5,194.5,208)

	makeCovariance2.line_15 = ROOT.TLine(201.5,0.5,201.5,208)

	makeCovariance2.lineb_1 = ROOT.TLine(0.5,19.5,208,19.5)

	makeCovariance2.lineb_2 = ROOT.TLine(0.5,38.5,208,38.5)

	makeCovariance2.lineb_3 = ROOT.TLine(0.5,45.5,208,45.5)

	makeCovariance2.lineb_4 = ROOT.TLine(0.5,52.5,208,52.5)

	makeCovariance2.lineb_5 = ROOT.TLine(0.5,71.5,208,71.5)

	makeCovariance2.lineb_6 = ROOT.TLine(0.5,90.5,208,90.5)

	makeCovariance2.lineb_7 = ROOT.TLine(0.5,97.5,208,97.5)

	makeCovariance2.lineb_8 = ROOT.TLine(0.5,104.5,208,104.5)

	makeCovariance2.lineb_9 = ROOT.TLine(0.5,123.5,208,123.5)

	makeCovariance2.lineb_10 = ROOT.TLine(0.5,142.5,208,142.5)

	makeCovariance2.lineb_11 = ROOT.TLine(0.5,149.5,208,149.5)

	makeCovariance2.lineb_12 = ROOT.TLine(0.5,156.5,208,156.5)

	makeCovariance2.lineb_13 = ROOT.TLine(0.5,175.5,208,175.5)

	makeCovariance2.lineb_14 = ROOT.TLine(0.5,194.5,208,194.5)

	makeCovariance2.lineb_15 = ROOT.TLine(0.5,201.5,208,201.5)

	makeCovariance2.f_foc = ROOT.TFile("focusing_covariance_"+beam_config+".root")

	makeCovariance2.covariance_focusing = f_foc.Get("Focusing_Covariance")

	makeCovariance2.total_covariance = flux_covariance.Clone("total_covariance")

	makeCovariance2.latex0b = ROOT.TLatex( 0.5, 0.93, "Fractional Flux Covariance Matrix" );

	makeCovariance2.f_out = ROOT.TFile("total_covariance_DUNE_"+beam_config+".root","RECREATE")

	makeCovariance2.total_error = ROOT.TH1D("Total Error","Total Error",208,0.5,208.5);

	makeCovariance2.linec_1 = ROOT.TLine(19.5,0,19.5,total_error.GetMaximum())

	makeCovariance2.linec_2 = ROOT.TLine(38.5,0,38.5,total_error.GetMaximum())

	makeCovariance2.linec_3 = ROOT.TLine(45.5,0,45.5,total_error.GetMaximum())

	makeCovariance2.linec_4 = ROOT.TLine(52.5,0,52.5,total_error.GetMaximum())

	makeCovariance2.linec_5 = ROOT.TLine(71.5,0,71.5,total_error.GetMaximum())

	makeCovariance2.linec_6 = ROOT.TLine(90.5,0,90.5,total_error.GetMaximum())

	makeCovariance2.linec_7 = ROOT.TLine(97.5,0,97.5,total_error.GetMaximum())

	makeCovariance2.linec_8 = ROOT.TLine(104.5,0,104.5,total_error.GetMaximum())

	makeCovariance2.linec_9 = ROOT.TLine(123.5,0,123.5,total_error.GetMaximum())

	makeCovariance2.linec_10 = ROOT.TLine(142.5,0,142.5,total_error.GetMaximum())

	makeCovariance2.linec_11 = ROOT.TLine(149.5,0,149.5,total_error.GetMaximum())

	makeCovariance2.linec_12 = ROOT.TLine(156.5,0,156.5,total_error.GetMaximum())

	makeCovariance2.linec_13 = ROOT.TLine(175.5,0,175.5,total_error.GetMaximum())

	makeCovariance2.linec_14 = ROOT.TLine(194.5,0,194.5,total_error.GetMaximum())

	makeCovariance2.linec_15 = ROOT.TLine(201.5,0,201.5,total_error.GetMaximum())

int	makeCovariance2.NRGBs = 3

int	makeCovariance2.NCont = 999

list	makeCovariance2.stops = [ 0.00, 0.50, 1.00];

list	makeCovariance2.red = [ 0.00, 1.00, 1.00];

list	makeCovariance2.green = [ 0.00, 1.00, 0.00];

list	makeCovariance2.blue = [ 1.00, 1.00, 0.00];

	makeCovariance2.correlation = GetCorrelation(total_covariance)

list	makeCovariance2.error_bins = [0,2,4,6,8,10]

list	makeCovariance2.error_hists = []

list	makeCovariance2.total_uncertainties = [0]*(len(error_bins)-1)

	makeCovariance2.temp_hist = ROOT.TH1D(error_types[i],error_types[i],len(error_bins)-1,array.array('d',error_bins))

list	makeCovariance2.colors = [2,ROOT.kOrange+3,6,ROOT.kOrange+7,ROOT.kGreen+2,3,ROOT.kAzure+7,ROOT.kBlue-4,7,ROOT.kViolet+1,13,15,17,19]

	makeCovariance2.leg = ROOT.TLegend(0.25,0.6,0.5,0.88)

	makeCovariance2.focusing_unc = ROOT.TH1D("Focusing","Focusing",len(error_bins)-1,array.array('d',error_bins))

	makeCovariance2.total_unc = ROOT.TH1D("Total","Total",len(error_bins)-1,array.array('d',error_bins))

	makeCovariance2.leg2 = ROOT.TLegend(0.25,0.7,0.6,0.88)

Namespaces

Functions

Variables