|
| string | makeCovariance.beam_config = "opt" |
| |
| list | makeCovariance.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 | makeCovariance.start_indices = [0,19,38,45,52,71,90,97,104,123,142,149,156,175,194,201] |
| |
| list | makeCovariance.in_files = [] |
| |
| list | makeCovariance.cv_histos = [] |
| |
| int | makeCovariance.n_universes = 100 |
| |
| list | makeCovariance.error_universes = [] |
| |
| list | makeCovariance.error_types = ["total","others","targetpcpi","targetpck","targetncpi","targetpcnu","targetmesonInc","targetnuA","tgtabsorption","totabsorption"] |
| |
| list | makeCovariance.error_types2 = ["total","others","targetpcpi","targetpcK","targetncpi","targetpcnu","targetmesonInc","targetnucleonA","targetabsorption","absorption"] |
| |
| list | makeCovariance.error_titles = ["Total","Other","pC#rightarrow#pi","pC#rightarrowK","nC#rightarrow#pi","pc#rightarrowN","Meson Inc","NucleonA","Target Absorption","Other Absorption"] |
| |
| list | makeCovariance.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 | makeCovariance.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 | makeCovariance.merged_histograms = [] |
| |
| list | makeCovariance.merged_universes = [] |
| |
| list | makeCovariance.flux_uncertainties = [] |
| |
| | makeCovariance.c1 = ROOT.TCanvas("c1","",1000,800) |
| |
| | makeCovariance.flux_covariance = MakeCovariance(merged_universes[error_types.index("total")]) |
| |
| | makeCovariance.latex0 = ROOT.TLatex( 0.5, 0.93, "Flux Hadron Production Fractional Covariance Matrix" ); |
| |
| | makeCovariance.latex1 = ROOT.TLatex( 0.24, 0.035, "Near Detector" ); |
| |
| | makeCovariance.latex2 = latex1.Clone() |
| |
| | makeCovariance.latex3 = latex1.Clone() |
| |
| | makeCovariance.latex4 = latex1.Clone() |
| |
| | makeCovariance.latex5 = ROOT.TLatex(0.2,0.08,"#nu Mode") |
| |
| | makeCovariance.latex6 = latex5.Clone() |
| |
| | makeCovariance.latex7 = latex5.Clone() |
| |
| | makeCovariance.latex8 = latex5.Clone() |
| |
| | makeCovariance.latex9 = ROOT.TLatex(0.085,0.21,"#nu Mode") |
| |
| | makeCovariance.latex10 = latex9.Clone() |
| |
| | makeCovariance.latex11 = latex9.Clone() |
| |
| | makeCovariance.latex12 = latex9.Clone() |
| |
| | makeCovariance.latex13 = ROOT.TLatex( 0.175, 0.124, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex14 = ROOT.TLatex( 0.245, 0.124, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex15 = ROOT.TLatex( 0.285, 0.124, "#nu_{e}" ); |
| |
| | makeCovariance.latex16 = ROOT.TLatex( 0.31, 0.124, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex17 = ROOT.TLatex( 0.35, 0.124, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex18 = ROOT.TLatex( 0.415, 0.124, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex19 = ROOT.TLatex( 0.46, 0.124, "#nu_{e}" ); |
| |
| | makeCovariance.latex20 = ROOT.TLatex( 0.485, 0.124, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex21 = ROOT.TLatex( 0.53, 0.124, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex22 = ROOT.TLatex( 0.595, 0.124, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex23 = ROOT.TLatex( 0.631, 0.124, "#nu_{e}" ); |
| |
| | makeCovariance.latex24 = ROOT.TLatex( 0.654, 0.124, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex25 = ROOT.TLatex( 0.70, 0.124, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex26 = ROOT.TLatex( 0.77, 0.124, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex27 = ROOT.TLatex( 0.81, 0.124, "#nu_{e}" ); |
| |
| | makeCovariance.latex28 = ROOT.TLatex( 0.835, 0.124, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex29 = ROOT.TLatex( 0.11,0.18, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex30 = ROOT.TLatex( 0.11,0.25, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex31 = ROOT.TLatex( 0.11,0.295, "#nu_{e}" ); |
| |
| | makeCovariance.latex32 = ROOT.TLatex( 0.11,0.32, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex33 = ROOT.TLatex( 0.11,0.37, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex34 = ROOT.TLatex( 0.11,0.435, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex35 = ROOT.TLatex( 0.11,0.48, "#nu_{e}" ); |
| |
| | makeCovariance.latex36 = ROOT.TLatex( 0.11,0.505, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex37 = ROOT.TLatex( 0.11,0.56, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex38 = ROOT.TLatex( 0.11,0.62, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex39 = ROOT.TLatex( 0.11,0.67, "#nu_{e}" ); |
| |
| | makeCovariance.latex40 = ROOT.TLatex( 0.11,0.70, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.latex41 = ROOT.TLatex( 0.11,0.75, "#nu_{#mu}" ); |
| |
| | makeCovariance.latex42 = ROOT.TLatex( 0.11,0.81, "#bar{#nu}_{#mu}" ); |
| |
| | makeCovariance.latex43 = ROOT.TLatex( 0.11,0.855, "#nu_{e}" ); |
| |
| | makeCovariance.latex44 = ROOT.TLatex( 0.11,0.88, "#bar{#nu}_{e}" ); |
| |
| | makeCovariance.line_1 = ROOT.TLine(19.5,0.5,19.5,208) |
| |
| | makeCovariance.line_2 = ROOT.TLine(38.5,0.5,38.5,208) |
| |
| | makeCovariance.line_3 = ROOT.TLine(45.5,0.5,45.5,208) |
| |
| | makeCovariance.line_4 = ROOT.TLine(52.5,0.5,52.5,208) |
| |
| | makeCovariance.line_5 = ROOT.TLine(71.5,0.5,71.5,208) |
| |
| | makeCovariance.line_6 = ROOT.TLine(90.5,0.5,90.5,208) |
| |
| | makeCovariance.line_7 = ROOT.TLine(97.5,0.5,97.5,208) |
| |
| | makeCovariance.line_8 = ROOT.TLine(104.5,0.5,104.5,208) |
| |
| | makeCovariance.line_9 = ROOT.TLine(123.5,0.5,123.5,208) |
| |
| | makeCovariance.line_10 = ROOT.TLine(142.5,0.5,142.5,208) |
| |
| | makeCovariance.line_11 = ROOT.TLine(149.5,0.5,149.5,208) |
| |
| | makeCovariance.line_12 = ROOT.TLine(156.5,0.5,156.5,208) |
| |
| | makeCovariance.line_13 = ROOT.TLine(175.5,0.5,175.5,208) |
| |
| | makeCovariance.line_14 = ROOT.TLine(194.5,0.5,194.5,208) |
| |
| | makeCovariance.line_15 = ROOT.TLine(201.5,0.5,201.5,208) |
| |
| | makeCovariance.lineb_1 = ROOT.TLine(0.5,19.5,208,19.5) |
| |
| | makeCovariance.lineb_2 = ROOT.TLine(0.5,38.5,208,38.5) |
| |
| | makeCovariance.lineb_3 = ROOT.TLine(0.5,45.5,208,45.5) |
| |
| | makeCovariance.lineb_4 = ROOT.TLine(0.5,52.5,208,52.5) |
| |
| | makeCovariance.lineb_5 = ROOT.TLine(0.5,71.5,208,71.5) |
| |
| | makeCovariance.lineb_6 = ROOT.TLine(0.5,90.5,208,90.5) |
| |
| | makeCovariance.lineb_7 = ROOT.TLine(0.5,97.5,208,97.5) |
| |
| | makeCovariance.lineb_8 = ROOT.TLine(0.5,104.5,208,104.5) |
| |
| | makeCovariance.lineb_9 = ROOT.TLine(0.5,123.5,208,123.5) |
| |
| | makeCovariance.lineb_10 = ROOT.TLine(0.5,142.5,208,142.5) |
| |
| | makeCovariance.lineb_11 = ROOT.TLine(0.5,149.5,208,149.5) |
| |
| | makeCovariance.lineb_12 = ROOT.TLine(0.5,156.5,208,156.5) |
| |
| | makeCovariance.lineb_13 = ROOT.TLine(0.5,175.5,208,175.5) |
| |
| | makeCovariance.lineb_14 = ROOT.TLine(0.5,194.5,208,194.5) |
| |
| | makeCovariance.lineb_15 = ROOT.TLine(0.5,201.5,208,201.5) |
| |
| | makeCovariance.f_foc = ROOT.TFile("focusing_covariance_"+beam_config+".root") |
| |
| | makeCovariance.covariance_focusing = f_foc.Get("Focusing_Covariance") |
| |
| | makeCovariance.total_covariance = flux_covariance.Clone("total_covariance") |
| |
| | makeCovariance.latex0b = ROOT.TLatex( 0.5, 0.93, "Fractional Flux Covariance Matrix" ); |
| |
| | makeCovariance.f_out = ROOT.TFile("total_covariance_DUNE_"+beam_config+".root","RECREATE") |
| |
| | makeCovariance.total_error = ROOT.TH1D("Total Error","Total Error",208,0.5,208.5); |
| |
| | makeCovariance.linec_1 = ROOT.TLine(19.5,0,19.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_2 = ROOT.TLine(38.5,0,38.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_3 = ROOT.TLine(45.5,0,45.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_4 = ROOT.TLine(52.5,0,52.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_5 = ROOT.TLine(71.5,0,71.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_6 = ROOT.TLine(90.5,0,90.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_7 = ROOT.TLine(97.5,0,97.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_8 = ROOT.TLine(104.5,0,104.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_9 = ROOT.TLine(123.5,0,123.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_10 = ROOT.TLine(142.5,0,142.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_11 = ROOT.TLine(149.5,0,149.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_12 = ROOT.TLine(156.5,0,156.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_13 = ROOT.TLine(175.5,0,175.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_14 = ROOT.TLine(194.5,0,194.5,total_error.GetMaximum()) |
| |
| | makeCovariance.linec_15 = ROOT.TLine(201.5,0,201.5,total_error.GetMaximum()) |
| |
| int | makeCovariance.NRGBs = 3 |
| |
| int | makeCovariance.NCont = 999 |
| |
| list | makeCovariance.stops = [ 0.00, 0.50, 1.00]; |
| |
| list | makeCovariance.red = [ 0.00, 1.00, 1.00]; |
| |
| list | makeCovariance.green = [ 0.00, 1.00, 0.00]; |
| |
| list | makeCovariance.blue = [ 1.00, 1.00, 0.00]; |
| |
| | makeCovariance.correlation = GetCorrelation(total_covariance) |
| |
| list | makeCovariance.error_bins = [0,2,4,6,8,10] |
| |
| list | makeCovariance.error_hists = [] |
| |
| list | makeCovariance.total_uncertainties = [0]*(len(error_bins)-1) |
| |
| | makeCovariance.temp_hist = ROOT.TH1D(error_types[i],error_types[i],len(error_bins)-1,array.array('d',error_bins)) |
| |
| list | makeCovariance.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] |
| |
| | makeCovariance.leg = ROOT.TLegend(0.25,0.6,0.5,0.88) |
| |
| | makeCovariance.focusing_unc = ROOT.TH1D("Focusing","Focusing",len(error_bins)-1,array.array('d',error_bins)) |
| |
| | makeCovariance.total_unc = ROOT.TH1D("Total","Total",len(error_bins)-1,array.array('d',error_bins)) |
| |
| | makeCovariance.leg2 = ROOT.TLegend(0.25,0.7,0.6,0.88) |
| |
1.8.11
