Functions | Variables
makeFocusingCovariance Namespace Reference

Functions

def GetCorrelation (cov_matrix)
 

Variables

bool small_bins = False
 
string beam_config = "opt"
 
list shifts = []
 
list cv = []
 
list uncertainties = ["HornCurrent","WaterThickness","DecayPipeRadius","Horn1XShift","Horn1XTilt"]
 
list uncertainty_titles = ["Horn Current","Water Layer","Decay Pipe Radius","Horn 1 Trans. Offset","Horn 1 Tilt"]
 
string indir = "/dune/app/users/pmadigan/beam_sim/ToleranceStudy/"
 
list files = ["table_ND_numu.csv","table_ND_numubar.csv","table_ND_nue.csv","table_ND_nuebar.csv","table_ND_numubar.csv","table_ND_numu.csv","table_ND_nuebar.csv","table_ND_nue.csv","table_FD_numu.csv","table_FD_numubar.csv","table_FD_nue.csv","table_FD_nuebar.csv","table_FD_numubar.csv","table_FD_numu.csv","table_FD_nuebar.csv","table_FD_nue.csv"]
 
 f_temp = open(indir+afile)
 
 lines = f_temp.readlines()
 
 temp = lines[i].split(",")[7:]
 
int line_iter = -1
 
float scale = 1.0
 
list temp_2 = []
 
int nbins = 7
 
tuple new_shift = (temp[12]/100*cv[12]+temp[13]/100*cv[13])/(cv[12]+cv[13])
 
list temp_uncertainties = ["Target","BeamSigma","BeamOffsetX","TargetDensity"]
 
list sigmas = ["0p5mm","1p7mm","0p5mm","1p8156gcm3"]
 
 unc_iter2 = temp_uncertainties.index(uncertainty)
 
list detectors = ["ND","FD"]
 
list modes = ["neutrino","antineutrino"]
 
list flavors = ["numu","numubar","nue","nuebar"]
 
 varied_file = ROOT.TFile("/dune/data/users/ljf26/fluxfiles/g4lbne/v3r4p3/QGSP_BERT/CP_run15_12388_80GeV_"+uncertainty+sigmas[unc_iter2]+"/"+mode+"/flux/histos_g4lbne_v3r4p3_QGSP_BERT_CP_run15_12388_80GeV_"+uncertainty+sigmas[unc_iter2]+"_"+mode+"_LBNE"+detector+"_fastmc.root")
 
 varied_histo = varied_file.Get(flavor+"_flux")
 
 nominal_file = ROOT.TFile("/dune/data/users/ljf26/fluxfiles/g4lbne/v3r4p2/QGSP_BERT/CP_run15_12388_80GeV/"+mode+"/flux/histos_g4lbne_v3r4p2_QGSP_BERT_CP_run15_12388_80GeV_"+mode+"_LBNE"+detector+"_fastmc.root")
 
 nominal_histo = nominal_file.Get(flavor+"_flux")
 
list varied_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]
 
 baffle_scraping_file = ROOT.TFile("/dune/data/users/ljf26/fluxfiles/g4lbne/v3r4p3/QGSP_BERT/CP_run15_12388_80GeV_BaffleScraping/"+mode+"/flux/histos_g4lbne_v3r4p3_QGSP_BERT_CP_run15_12388_80GeV_BaffleScraping_"+mode+"_LBNE"+detector+"_fastmc.root")
 
 baffle_scraping_histo = baffle_scraping_file.Get(flavor+"_flux")
 
list baffle_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]
 
 unc_iter = uncertainties.index(uncertainty)
 
 error = ROOT.TH1D("Focusing_Error","Focusing_Error",nbins,0.5,nbins+0.5)
 
 covariance = ROOT.TH2D("Focusing_Covariance","Focusing_Covariance",nbins,0.5,nbins+0.5,nbins,0.5,nbins+0.5)
 
 c1 = ROOT.TCanvas("c1","",1000,800)
 
 latex0 = ROOT.TLatex( 0.5, 0.93, "Flux Alignment 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.16, 0.124, "#nu_{#mu}" );
 
 latex14 = ROOT.TLatex( 0.225, 0.124, "#bar{#nu}_{#mu}" );
 
 latex15 = ROOT.TLatex( 0.27, 0.124, "#nu_{e}" );
 
 latex16 = ROOT.TLatex( 0.295, 0.124, "#bar{#nu}_{e}" );
 
 latex17 = ROOT.TLatex( 0.34, 0.124, "#nu_{#mu}" );
 
 latex18 = ROOT.TLatex( 0.405, 0.124, "#bar{#nu}_{#mu}" );
 
 latex19 = ROOT.TLatex( 0.45, 0.124, "#nu_{e}" );
 
 latex20 = ROOT.TLatex( 0.475, 0.124, "#bar{#nu}_{e}" );
 
 latex21 = ROOT.TLatex( 0.52, 0.124, "#nu_{#mu}" );
 
 latex22 = ROOT.TLatex( 0.585, 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}" );
 
list line_locations = [19.5,38.5,45.5,52.5,71.5,90.5,97.5,104.5,123.5,142.5,149.5,156.5,175.5,194.5,201.5]
 
 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)
 
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(covariance)
 
string small_bins_string = ""
 
 f = ROOT.TFile("focusing_covariance_"+beam_config+small_bins_string+".root","RECREATE")
 
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 error_bins = [0,2,4,6,8,10]
 
list error_hists = []
 
list total_uncertainties = [0]*(len(error_bins)-1)
 
 temp_hist = ROOT.TH1D(uncertainty_titles[i],uncertainty_titles[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.57,0.5,0.92,0.85)
 

Function Documentation

def makeFocusingCovariance.GetCorrelation (   cov_matrix)

Definition at line 10 of file makeFocusingCovariance.py.

10 def GetCorrelation(cov_matrix):
11  cor_matrix = cov_matrix.Clone("Focusing_Correlation")
12  for i in range(cor_matrix.GetNbinsX()):
13  for j in range(cor_matrix.GetNbinsX()):
14  if cov_matrix.GetBinContent(i+1,i+1)!=0 and cov_matrix.GetBinContent(j+1,j+1)!=0:
15  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)))
16  else:
17  if i==j: cor_matrix.SetBinContent(i+1,j+1,1)
18  else: cor_matrix.SetBinContent(i+1,j+1,0)
19  cor_matrix.SetMinimum(-1.0)
20  cor_matrix.SetMaximum(1.0)
21  return cor_matrix
22 
23 
24 # Parse input files
25 

Variable Documentation

list makeFocusingCovariance.baffle_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 189 of file makeFocusingCovariance.py.

makeFocusingCovariance.baffle_scraping_file = ROOT.TFile("/dune/data/users/ljf26/fluxfiles/g4lbne/v3r4p3/QGSP_BERT/CP_run15_12388_80GeV_BaffleScraping/"+mode+"/flux/histos_g4lbne_v3r4p3_QGSP_BERT_CP_run15_12388_80GeV_BaffleScraping_"+mode+"_LBNE"+detector+"_fastmc.root")

Definition at line 179 of file makeFocusingCovariance.py.

makeFocusingCovariance.baffle_scraping_histo = baffle_scraping_file.Get(flavor+"_flux")

Definition at line 180 of file makeFocusingCovariance.py.

string makeFocusingCovariance.beam_config = "opt"

Definition at line 8 of file makeFocusingCovariance.py.

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

Definition at line 711 of file makeFocusingCovariance.py.

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

Definition at line 313 of file makeFocusingCovariance.py.

list makeFocusingCovariance.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 776 of file makeFocusingCovariance.py.

makeFocusingCovariance.correlation = GetCorrelation(covariance)

Definition at line 714 of file makeFocusingCovariance.py.

makeFocusingCovariance.covariance = ROOT.TH2D("Focusing_Covariance","Focusing_Covariance",nbins,0.5,nbins+0.5,nbins,0.5,nbins+0.5)

Definition at line 304 of file makeFocusingCovariance.py.

makeFocusingCovariance.cv = []

Definition at line 27 of file makeFocusingCovariance.py.

list makeFocusingCovariance.detectors = ["ND","FD"]

Definition at line 124 of file makeFocusingCovariance.py.

makeFocusingCovariance.error = ROOT.TH1D("Focusing_Error","Focusing_Error",nbins,0.5,nbins+0.5)

Definition at line 303 of file makeFocusingCovariance.py.

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

Definition at line 743 of file makeFocusingCovariance.py.

list makeFocusingCovariance.error_hists = []

Definition at line 746 of file makeFocusingCovariance.py.

makeFocusingCovariance.f = ROOT.TFile("focusing_covariance_"+beam_config+small_bins_string+".root","RECREATE")

Definition at line 728 of file makeFocusingCovariance.py.

makeFocusingCovariance.f_temp = open(indir+afile)

Definition at line 45 of file makeFocusingCovariance.py.

list makeFocusingCovariance.files = ["table_ND_numu.csv","table_ND_numubar.csv","table_ND_nue.csv","table_ND_nuebar.csv","table_ND_numubar.csv","table_ND_numu.csv","table_ND_nuebar.csv","table_ND_nue.csv","table_FD_numu.csv","table_FD_numubar.csv","table_FD_nue.csv","table_FD_nuebar.csv","table_FD_numubar.csv","table_FD_numu.csv","table_FD_nuebar.csv","table_FD_nue.csv"]

Definition at line 39 of file makeFocusingCovariance.py.

list makeFocusingCovariance.flavors = ["numu","numubar","nue","nuebar"]

Definition at line 126 of file makeFocusingCovariance.py.

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

Definition at line 710 of file makeFocusingCovariance.py.

string makeFocusingCovariance.indir = "/dune/app/users/pmadigan/beam_sim/ToleranceStudy/"

Definition at line 36 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex0 = ROOT.TLatex( 0.5, 0.93, "Flux Alignment Fractional Covariance Matrix" );

Definition at line 329 of file makeFocusingCovariance.py.

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

Definition at line 337 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex10 = latex9.Clone()

Definition at line 375 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex11 = latex9.Clone()

Definition at line 378 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex12 = latex9.Clone()

Definition at line 381 of file makeFocusingCovariance.py.

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

Definition at line 385 of file makeFocusingCovariance.py.

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

Definition at line 390 of file makeFocusingCovariance.py.

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

Definition at line 395 of file makeFocusingCovariance.py.

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

Definition at line 400 of file makeFocusingCovariance.py.

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

Definition at line 406 of file makeFocusingCovariance.py.

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

Definition at line 411 of file makeFocusingCovariance.py.

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

Definition at line 416 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex2 = latex1.Clone()

Definition at line 342 of file makeFocusingCovariance.py.

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

Definition at line 421 of file makeFocusingCovariance.py.

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

Definition at line 427 of file makeFocusingCovariance.py.

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

Definition at line 432 of file makeFocusingCovariance.py.

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

Definition at line 437 of file makeFocusingCovariance.py.

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

Definition at line 442 of file makeFocusingCovariance.py.

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

Definition at line 448 of file makeFocusingCovariance.py.

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

Definition at line 453 of file makeFocusingCovariance.py.

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

Definition at line 458 of file makeFocusingCovariance.py.

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

Definition at line 463 of file makeFocusingCovariance.py.

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

Definition at line 469 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex3 = latex1.Clone()

Definition at line 345 of file makeFocusingCovariance.py.

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

Definition at line 475 of file makeFocusingCovariance.py.

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

Definition at line 481 of file makeFocusingCovariance.py.

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

Definition at line 487 of file makeFocusingCovariance.py.

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

Definition at line 494 of file makeFocusingCovariance.py.

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

Definition at line 500 of file makeFocusingCovariance.py.

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

Definition at line 506 of file makeFocusingCovariance.py.

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

Definition at line 512 of file makeFocusingCovariance.py.

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

Definition at line 519 of file makeFocusingCovariance.py.

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

Definition at line 525 of file makeFocusingCovariance.py.

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

Definition at line 531 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex4 = latex1.Clone()

Definition at line 349 of file makeFocusingCovariance.py.

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

Definition at line 537 of file makeFocusingCovariance.py.

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

Definition at line 544 of file makeFocusingCovariance.py.

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

Definition at line 550 of file makeFocusingCovariance.py.

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

Definition at line 556 of file makeFocusingCovariance.py.

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

Definition at line 562 of file makeFocusingCovariance.py.

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

Definition at line 354 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex6 = latex5.Clone()

Definition at line 359 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex7 = latex5.Clone()

Definition at line 362 of file makeFocusingCovariance.py.

makeFocusingCovariance.latex8 = latex5.Clone()

Definition at line 365 of file makeFocusingCovariance.py.

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

Definition at line 369 of file makeFocusingCovariance.py.

makeFocusingCovariance.leg = ROOT.TLegend(0.57,0.5,0.92,0.85)

Definition at line 788 of file makeFocusingCovariance.py.

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

Definition at line 657 of file makeFocusingCovariance.py.

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

Definition at line 666 of file makeFocusingCovariance.py.

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

Definition at line 667 of file makeFocusingCovariance.py.

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

Definition at line 668 of file makeFocusingCovariance.py.

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

Definition at line 669 of file makeFocusingCovariance.py.

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

Definition at line 670 of file makeFocusingCovariance.py.

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

Definition at line 671 of file makeFocusingCovariance.py.

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

Definition at line 658 of file makeFocusingCovariance.py.

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

Definition at line 659 of file makeFocusingCovariance.py.

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

Definition at line 660 of file makeFocusingCovariance.py.

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

Definition at line 661 of file makeFocusingCovariance.py.

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

Definition at line 662 of file makeFocusingCovariance.py.

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

Definition at line 663 of file makeFocusingCovariance.py.

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

Definition at line 664 of file makeFocusingCovariance.py.

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

Definition at line 665 of file makeFocusingCovariance.py.

makeFocusingCovariance.line_iter = -1

Definition at line 56 of file makeFocusingCovariance.py.

list makeFocusingCovariance.line_locations = [19.5,38.5,45.5,52.5,71.5,90.5,97.5,104.5,123.5,142.5,149.5,156.5,175.5,194.5,201.5]

Definition at line 653 of file makeFocusingCovariance.py.

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

Definition at line 672 of file makeFocusingCovariance.py.

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

Definition at line 681 of file makeFocusingCovariance.py.

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

Definition at line 682 of file makeFocusingCovariance.py.

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

Definition at line 683 of file makeFocusingCovariance.py.

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

Definition at line 684 of file makeFocusingCovariance.py.

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

Definition at line 685 of file makeFocusingCovariance.py.

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

Definition at line 686 of file makeFocusingCovariance.py.

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

Definition at line 673 of file makeFocusingCovariance.py.

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

Definition at line 674 of file makeFocusingCovariance.py.

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

Definition at line 675 of file makeFocusingCovariance.py.

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

Definition at line 676 of file makeFocusingCovariance.py.

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

Definition at line 677 of file makeFocusingCovariance.py.

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

Definition at line 678 of file makeFocusingCovariance.py.

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

Definition at line 679 of file makeFocusingCovariance.py.

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

Definition at line 680 of file makeFocusingCovariance.py.

makeFocusingCovariance.lines = f_temp.readlines()

Definition at line 46 of file makeFocusingCovariance.py.

list makeFocusingCovariance.modes = ["neutrino","antineutrino"]

Definition at line 125 of file makeFocusingCovariance.py.

int makeFocusingCovariance.nbins = 7

Definition at line 79 of file makeFocusingCovariance.py.

int makeFocusingCovariance.NCont = 999

Definition at line 706 of file makeFocusingCovariance.py.

tuple makeFocusingCovariance.new_shift = (temp[12]/100*cv[12]+temp[13]/100*cv[13])/(cv[12]+cv[13])

Definition at line 93 of file makeFocusingCovariance.py.

makeFocusingCovariance.nominal_file = ROOT.TFile("/dune/data/users/ljf26/fluxfiles/g4lbne/v3r4p2/QGSP_BERT/CP_run15_12388_80GeV/"+mode+"/flux/histos_g4lbne_v3r4p2_QGSP_BERT_CP_run15_12388_80GeV_"+mode+"_LBNE"+detector+"_fastmc.root")

Definition at line 136 of file makeFocusingCovariance.py.

makeFocusingCovariance.nominal_histo = nominal_file.Get(flavor+"_flux")

Definition at line 137 of file makeFocusingCovariance.py.

int makeFocusingCovariance.NRGBs = 3

Definition at line 705 of file makeFocusingCovariance.py.

list makeFocusingCovariance.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 734 of file makeFocusingCovariance.py.

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

Definition at line 709 of file makeFocusingCovariance.py.

float makeFocusingCovariance.scale = 1.0

Definition at line 72 of file makeFocusingCovariance.py.

list makeFocusingCovariance.shifts = []

Definition at line 26 of file makeFocusingCovariance.py.

list makeFocusingCovariance.sigmas = ["0p5mm","1p7mm","0p5mm","1p8156gcm3"]

Definition at line 117 of file makeFocusingCovariance.py.

bool makeFocusingCovariance.small_bins = False

Definition at line 5 of file makeFocusingCovariance.py.

makeFocusingCovariance.small_bins_string = ""

Definition at line 723 of file makeFocusingCovariance.py.

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

Definition at line 736 of file makeFocusingCovariance.py.

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

Definition at line 708 of file makeFocusingCovariance.py.

makeFocusingCovariance.temp = lines[i].split(",")[7:]

Definition at line 50 of file makeFocusingCovariance.py.

makeFocusingCovariance.temp_2 = []

Definition at line 77 of file makeFocusingCovariance.py.

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

Definition at line 750 of file makeFocusingCovariance.py.

list makeFocusingCovariance.temp_uncertainties = ["Target","BeamSigma","BeamOffsetX","TargetDensity"]

Definition at line 114 of file makeFocusingCovariance.py.

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

Definition at line 747 of file makeFocusingCovariance.py.

makeFocusingCovariance.unc_iter = uncertainties.index(uncertainty)

Definition at line 216 of file makeFocusingCovariance.py.

makeFocusingCovariance.unc_iter2 = temp_uncertainties.index(uncertainty)

Definition at line 122 of file makeFocusingCovariance.py.

list makeFocusingCovariance.uncertainties = ["HornCurrent","WaterThickness","DecayPipeRadius","Horn1XShift","Horn1XTilt"]

Definition at line 32 of file makeFocusingCovariance.py.

list makeFocusingCovariance.uncertainty_titles = ["Horn Current","Water Layer","Decay Pipe Radius","Horn 1 Trans. Offset","Horn 1 Tilt"]

Definition at line 34 of file makeFocusingCovariance.py.

list makeFocusingCovariance.varied_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 141 of file makeFocusingCovariance.py.

makeFocusingCovariance.varied_file = ROOT.TFile("/dune/data/users/ljf26/fluxfiles/g4lbne/v3r4p3/QGSP_BERT/CP_run15_12388_80GeV_"+uncertainty+sigmas[unc_iter2]+"/"+mode+"/flux/histos_g4lbne_v3r4p3_QGSP_BERT_CP_run15_12388_80GeV_"+uncertainty+sigmas[unc_iter2]+"_"+mode+"_LBNE"+detector+"_fastmc.root")

Definition at line 131 of file makeFocusingCovariance.py.

makeFocusingCovariance.varied_histo = varied_file.Get(flavor+"_flux")

Definition at line 132 of file makeFocusingCovariance.py.