Functions | Variables
makeCovariance Namespace Reference

Functions

def MakeErrorHistogram (universes, asFrac=True)
 
def MakeCovariance (universes, asFrac=True)
 
def MergeHistograms (cvs, universes, name)
 
def GetCorrelation (cov_matrix)
 

Variables

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

Function Documentation

def makeCovariance.GetCorrelation (   cov_matrix)

Definition at line 85 of file makeCovariance.py.

85 def GetCorrelation(cov_matrix):
86  cor_matrix = cov_matrix.Clone()
87  for i in range(cor_matrix.GetNbinsX()):
88  for j in range(cor_matrix.GetNbinsX()):
89  if cov_matrix.GetBinContent(i+1,i+1)!=0 and cov_matrix.GetBinContent(j+1,j+1)!=0:
90  cor_matrix.SetBinContent(i+1,j+1,cov_matrix.GetBinContent(i+1,j+1)/ROOT.TMath.Sqrt(cov_matrix.GetBinContent(i+1,i+1))/ROOT.TMath.Sqrt(cov_matrix.GetBinContent(j+1,j+1)))
91  else:
92  if i==j: cor_matrix.SetBinContent(i+1,j+1,1)
93  else: cor_matrix.SetBinContent(i+1,j+1,0)
94  cor_matrix.SetMinimum(-1.0)
95  cor_matrix.SetMaximum(1.0)
96  return cor_matrix
97 
98 
99 # Open ppfx output file
makeCovariance.GetCorrelation
def GetCorrelation(cov_matrix)
Definition: makeCovariance.py:85
def makeCovariance.MakeCovariance (   universes,
  asFrac = True 
)

Definition at line 41 of file makeCovariance.py.

41 def MakeCovariance(universes,asFrac=True):
42  n_universes = len(universes)
43  avg_hist = universes[0].Clone(universes[0].GetName()+"_avg")
44  avg_hist.Reset()
45  nbins = universes[0].GetNbinsX()
46 
47  for i in range(n_universes):
48  avg_hist.Add(universes[i])
49  avg_hist.Scale(1/float(n_universes))
50 
51  covariance = ROOT.TH2D("Flux_Covariance","Flux_Covariance",nbins,0.5,nbins+0.5,nbins,0.5,nbins+0.5)
52  # Calculate Error = RMS of universes
53  for j in range(avg_hist.GetNbinsX()):
54  for k in range(avg_hist.GetNbinsX()):
55  cov = 0
56  for i in range(n_universes):
57  cov += 1/float(n_universes)*(universes[i].GetBinContent(j+1)-avg_hist.GetBinContent(j+1))*(universes[i].GetBinContent(k+1)-avg_hist.GetBinContent(k+1))
58  if asFrac:
59  if(avg_hist.GetBinContent(j+1)*avg_hist.GetBinContent(k+1) != 0):
60  cov = cov / avg_hist.GetBinContent(j+1) / avg_hist.GetBinContent(k+1)
61  else:
62  cov = 0
63  covariance.SetBinContent(j+1,k+1,cov)
64  return covariance
65 
makeCovariance.MakeCovariance
def MakeCovariance(universes, asFrac=True)
Definition: makeCovariance.py:41
if
if(!yymsg) yymsg
def makeCovariance.MakeErrorHistogram (   universes,
  asFrac = True 
)

Definition at line 17 of file makeCovariance.py.

17 def MakeErrorHistogram(universes,asFrac=True):
18  n_universes = len(universes)
19  avg_hist = universes[0].Clone(universes[0].GetName()+"_avg")
20  error_hist = universes[0].Clone(universes[0].GetName()+"_error")
21  avg_hist.Reset()
22  error_hist.Reset()
23  avg_hist.Reset()
24  error_hist.Reset()
25  for i in range(n_universes):
26  avg_hist.Add(universes[i])
27  avg_hist.Scale(1/float(n_universes))
28 
29  # Calculate Error = RMS of universes
30  for j in range(error_hist.GetNbinsX()+2):
31  error = 0
32  for i in range(n_universes):
33  error += 1/float(n_universes)*ROOT.TMath.Power(universes[i].GetBinContent(j)-avg_hist.GetBinContent(j),2)
34  error_hist.SetBinContent(j,ROOT.TMath.Sqrt(error))
35  if asFrac:
36  error_hist.Divide(avg_hist)
37  return error_hist
38 
39 # Helper function that makes a covariance matrix out of a bunch of universes
40 # by calculating rms of universes in each bin
makeCovariance.MakeErrorHistogram
def MakeErrorHistogram(universes, asFrac=True)
Definition: makeCovariance.py:17
def makeCovariance.MergeHistograms (   cvs,
  universes,
  name 
)

Definition at line 66 of file makeCovariance.py.

66 def MergeHistograms(cvs,universes,name):
67  nbins = 0
68  for cv in cvs:
69  nbins += cv.GetNbinsX()
70  new_cv = ROOT.TH1D("Flux_CV","Flux_CV",nbins,0.5,nbins+0.5)
71  new_universes = []
72  for universe in range(len(universes[0])):
73  new_universes.append(ROOT.TH1D("Flux_Universe_"+str(universe)+name,"Flux_Universe_"+str(universe)+name,nbins,0.5,nbins+0.5))
74  global_bin_iter = 1
75  for cv in range(len(cvs)):
76  for local_bin_iter in range(cvs[cv].GetNbinsX()):
77 
78  new_cv.SetBinContent(global_bin_iter,cvs[cv].GetBinContent(local_bin_iter+1))
79  for universe in range(len(universes[0])):
80  new_universes[universe].SetBinContent(global_bin_iter,universes[cv][universe].GetBinContent(local_bin_iter+1))
81  global_bin_iter = global_bin_iter+1
82  return_array = [new_cv,new_universes]
83  return return_array
84 
makeCovariance.MergeHistograms
def MergeHistograms(cvs, universes, name)
Definition: makeCovariance.py:66
str
static QCString str
Definition: fortrancode.cpp:27098

Variable Documentation

string makeCovariance.beam_config = "opt"

Definition at line 4 of file makeCovariance.py.

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

Definition at line 568 of file makeCovariance.py.

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

Definition at line 168 of file makeCovariance.py.

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]

Definition at line 612 of file makeCovariance.py.

makeCovariance.correlation = GetCorrelation(total_covariance)

Definition at line 571 of file makeCovariance.py.

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

Definition at line 477 of file makeCovariance.py.

list makeCovariance.cv_histos = []

Definition at line 114 of file makeCovariance.py.

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

Definition at line 586 of file makeCovariance.py.

list makeCovariance.error_hists = []

Definition at line 589 of file makeCovariance.py.

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

Definition at line 126 of file makeCovariance.py.

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

Definition at line 124 of file makeCovariance.py.

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

Definition at line 125 of file makeCovariance.py.

list makeCovariance.error_universes = []

Definition at line 123 of file makeCovariance.py.

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

Definition at line 476 of file makeCovariance.py.

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

Definition at line 509 of file makeCovariance.py.

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"]

Definition at line 150 of file makeCovariance.py.

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

Definition at line 175 of file makeCovariance.py.

list makeCovariance.flux_uncertainties = []

Definition at line 159 of file makeCovariance.py.

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

Definition at line 628 of file makeCovariance.py.

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

Definition at line 567 of file makeCovariance.py.

list makeCovariance.in_files = []

Definition at line 100 of file makeCovariance.py.

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

Definition at line 188 of file makeCovariance.py.

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

Definition at line 492 of file makeCovariance.py.

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

Definition at line 196 of file makeCovariance.py.

makeCovariance.latex10 = latex9.Clone()

Definition at line 234 of file makeCovariance.py.

makeCovariance.latex11 = latex9.Clone()

Definition at line 237 of file makeCovariance.py.

makeCovariance.latex12 = latex9.Clone()

Definition at line 240 of file makeCovariance.py.

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

Definition at line 244 of file makeCovariance.py.

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

Definition at line 249 of file makeCovariance.py.

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

Definition at line 254 of file makeCovariance.py.

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

Definition at line 259 of file makeCovariance.py.

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

Definition at line 265 of file makeCovariance.py.

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

Definition at line 270 of file makeCovariance.py.

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

Definition at line 275 of file makeCovariance.py.

makeCovariance.latex2 = latex1.Clone()

Definition at line 201 of file makeCovariance.py.

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

Definition at line 280 of file makeCovariance.py.

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

Definition at line 286 of file makeCovariance.py.

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

Definition at line 291 of file makeCovariance.py.

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

Definition at line 296 of file makeCovariance.py.

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

Definition at line 301 of file makeCovariance.py.

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

Definition at line 307 of file makeCovariance.py.

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

Definition at line 312 of file makeCovariance.py.

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

Definition at line 317 of file makeCovariance.py.

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

Definition at line 322 of file makeCovariance.py.

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

Definition at line 328 of file makeCovariance.py.

makeCovariance.latex3 = latex1.Clone()

Definition at line 204 of file makeCovariance.py.

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

Definition at line 334 of file makeCovariance.py.

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

Definition at line 340 of file makeCovariance.py.

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

Definition at line 346 of file makeCovariance.py.

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

Definition at line 353 of file makeCovariance.py.

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

Definition at line 359 of file makeCovariance.py.

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

Definition at line 365 of file makeCovariance.py.

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

Definition at line 371 of file makeCovariance.py.

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

Definition at line 378 of file makeCovariance.py.

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

Definition at line 384 of file makeCovariance.py.

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

Definition at line 390 of file makeCovariance.py.

makeCovariance.latex4 = latex1.Clone()

Definition at line 208 of file makeCovariance.py.

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

Definition at line 396 of file makeCovariance.py.

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

Definition at line 403 of file makeCovariance.py.

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

Definition at line 409 of file makeCovariance.py.

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

Definition at line 415 of file makeCovariance.py.

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

Definition at line 421 of file makeCovariance.py.

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

Definition at line 213 of file makeCovariance.py.

makeCovariance.latex6 = latex5.Clone()

Definition at line 218 of file makeCovariance.py.

makeCovariance.latex7 = latex5.Clone()

Definition at line 221 of file makeCovariance.py.

makeCovariance.latex8 = latex5.Clone()

Definition at line 224 of file makeCovariance.py.

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

Definition at line 228 of file makeCovariance.py.

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

Definition at line 617 of file makeCovariance.py.

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

Definition at line 654 of file makeCovariance.py.

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

Definition at line 428 of file makeCovariance.py.

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

Definition at line 437 of file makeCovariance.py.

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

Definition at line 438 of file makeCovariance.py.

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

Definition at line 439 of file makeCovariance.py.

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

Definition at line 440 of file makeCovariance.py.

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

Definition at line 441 of file makeCovariance.py.

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

Definition at line 442 of file makeCovariance.py.

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

Definition at line 429 of file makeCovariance.py.

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

Definition at line 430 of file makeCovariance.py.

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

Definition at line 431 of file makeCovariance.py.

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

Definition at line 432 of file makeCovariance.py.

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

Definition at line 433 of file makeCovariance.py.

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

Definition at line 434 of file makeCovariance.py.

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

Definition at line 435 of file makeCovariance.py.

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

Definition at line 436 of file makeCovariance.py.

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

Definition at line 443 of file makeCovariance.py.

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

Definition at line 452 of file makeCovariance.py.

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

Definition at line 453 of file makeCovariance.py.

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

Definition at line 454 of file makeCovariance.py.

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

Definition at line 455 of file makeCovariance.py.

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

Definition at line 456 of file makeCovariance.py.

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

Definition at line 457 of file makeCovariance.py.

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

Definition at line 444 of file makeCovariance.py.

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

Definition at line 445 of file makeCovariance.py.

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

Definition at line 446 of file makeCovariance.py.

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

Definition at line 447 of file makeCovariance.py.

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

Definition at line 448 of file makeCovariance.py.

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

Definition at line 449 of file makeCovariance.py.

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

Definition at line 450 of file makeCovariance.py.

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

Definition at line 451 of file makeCovariance.py.

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

Definition at line 526 of file makeCovariance.py.

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

Definition at line 535 of file makeCovariance.py.

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

Definition at line 536 of file makeCovariance.py.

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

Definition at line 537 of file makeCovariance.py.

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

Definition at line 538 of file makeCovariance.py.

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

Definition at line 539 of file makeCovariance.py.

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

Definition at line 540 of file makeCovariance.py.

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

Definition at line 527 of file makeCovariance.py.

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

Definition at line 528 of file makeCovariance.py.

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

Definition at line 529 of file makeCovariance.py.

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

Definition at line 530 of file makeCovariance.py.

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

Definition at line 531 of file makeCovariance.py.

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

Definition at line 532 of file makeCovariance.py.

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

Definition at line 533 of file makeCovariance.py.

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

Definition at line 534 of file makeCovariance.py.

list makeCovariance.merged_histograms = []

Definition at line 157 of file makeCovariance.py.

list makeCovariance.merged_universes = []

Definition at line 158 of file makeCovariance.py.

int makeCovariance.n_universes = 100

Definition at line 122 of file makeCovariance.py.

int makeCovariance.NCont = 999

Definition at line 563 of file makeCovariance.py.

int makeCovariance.NRGBs = 3

Definition at line 562 of file makeCovariance.py.

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]

Definition at line 149 of file makeCovariance.py.

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"]

Definition at line 11 of file makeCovariance.py.

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

Definition at line 566 of file makeCovariance.py.

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

Definition at line 13 of file makeCovariance.py.

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

Definition at line 565 of file makeCovariance.py.

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

Definition at line 593 of file makeCovariance.py.

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

Definition at line 482 of file makeCovariance.py.

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

Definition at line 513 of file makeCovariance.py.

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

Definition at line 629 of file makeCovariance.py.

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

Definition at line 590 of file makeCovariance.py.

Generated by   doxygen 1.8.11