Functions | Variables
plot_xsec_2d Namespace Reference

Functions

def get_genie_universes (basefn)
 
def get_ppfx_universes (fppfx_other)
 
def up_down_counts_1sigma (count_list, cv)
 
def process_multiuniverse (isys, huniv)
 
def process_light_calibration_mue (isys, hidxdown, hidxup)
 
def process_ckv_calibshape (isys, hidx)
 
def process_normalization (isys)
 
def process_focusing (isys, sstart, nsys)
 
def databin (xbin, ybin)
 
def process_data ()
 
def process_mc ()
 
def setMC1D (hIn)
 
def setData1D (gIn)
 
def setUncer1D (hIn, ihist)
 

Variables

string infp = './'
 
list syst_names
 
int Nfinal_fe = 9
 
list final_fe = ["Flux (HP)","GENIE","Light","Calibration","Cherenkov","Calib. Shape","Mu Energy Scale","Normalization","Focusing"]
 
list final_fe_name = ["Flux_HP","GENIE","Light","Calibration","Cherenkov","CalibShape","MuEScale","Normalization","Focusing"]
 
int Nabins2X = 13
 
list leg_angle
 
list col_fe = [2,4,2,4,8,8,kOrange-3,kOrange-3,1]
 
list col_st = [1,1,2,2,1,2,1,2,2]
 
string tit_ke = 'Muon Kinetic Energy (GeV)'
 
string tit_nu = 'Neutrino Energy (GeV)'
 
string tit_fe = 'Fractional Uncertainties'
 
string tit_xs = '#sigma(cm^{2}/nucleon)'
 
 parser = argparse.ArgumentParser()
 
 action
 
 help
 
 args = parser.parse_args()
 
 use_data = args.data
 
list fgenie = []
 
 fmc = TFile(os.path.join(infp, 'nominal_xsec.root'))
 
 fppfx_other = TFile(os.path.join(infp, 'ppfx_others_v2.root'))
 
list hgenie = []
 
list hppfx = []
 
list hother = []
 
 hmc = fmc.Get('hXS2Dsec_nominal')
 
 hdt = fppfx_other.Get('xs_nominal/hXS2Dsec_cv')
 
list lhup = []
 
list lhdw = []
 
list lhfe = []
 
 NbinsX = hdt.GetNbinsX()
 
 NbinsY = hdt.GetNbinsY()
 
float fe_norm = 0.028422526
 
list gdt = [TGraphAsymmErrors() for i in range(Nabins2X)]
 
list count = [0 for i in range(Nabins2X)]
 
list h1MC = [TH1D('h1MC_{}'.format(i),'',20,0.5,2.5) for i in range(Nabins2X)]
 
list h1fe = [[TH1D('h1fe_{}_{}'.format(i,j),'',20,0.5,2.5) for j in range(Nfinal_fe)] for i in range(Nabins2X)]
 
string out_subdir = 'real_data'
 
 dest_path = os.path.join('output', out_subdir)
 
list cxs = []
 
list leg = [TLegend(0.67,0.68,0.87,0.88) for i in range(Nabins2X)]
 
list lat = [TLatex() for i in range(Nabins2X)]
 
int maxY = -1
 
list cfe = []
 
list latfe = [TLatex() for i in range(Nabins2X)]
 
list lfe = []
 
 fileOut = TFile('{}/numucc_inc_{}.root'.format(dest_path,out_subdir), 'recreate')
 

Function Documentation

def plot_xsec_2d.databin (   xbin,
  ybin 
)

Definition at line 203 of file plot_xsec_2d.py.

Referenced by process_data(), and process_mc().

203 def databin(xbin, ybin):
204  answer = False
205  if xbin==1 and ybin<=4: answer = True
206  if (xbin==2 or xbin==3) and ybin<=5: answer = True
207  if xbin==4 and ybin<=7: answer = True
208  if xbin==5 and ybin<=9: answer = True
209  if xbin==6 and ybin<=10: answer = True
210  if xbin==7 and ybin<=12: answer = True
211  if xbin==8 and ybin<=15: answer = True
212  if xbin==9 and ybin<=19: answer = True
213  if (xbin==10 or xbin==11 or xbin==12 or xbin==13) and ybin<=20: answer = True
214 
215  return answer
216 
217 
def databin(xbin, ybin)
def plot_xsec_2d.get_genie_universes (   basefn)

Definition at line 50 of file plot_xsec_2d.py.

References novadaq::HexUtils.format(), Get, makeTrainCVSamples.int, and PandAna.Demos.demo1.range.

50 def get_genie_universes(basefn):
51 
52  hl = []
53  fn = os.path.join(infp, basefn)
54  fgenie.append(TFile(fn))
55  # find the genie universe IDs from filename
56  sIdx = int(basefn.split('_')[1])
57  eIdx = int(basefn.split('_')[2])+1
58  for i in range(sIdx,eIdx):
59  hl.append(fgenie[len(fgenie)-1].Get('xs_nominal/hXS2Dsec_genie{}'.format(i)))
60 
61  return hl
62 
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
Definition: AddMC.C:8
std::string format(const int32_t &value, const int &ndigits=8)
Definition: HexUtils.cpp:14
def get_genie_universes(basefn)
Definition: plot_xsec_2d.py:50
def plot_xsec_2d.get_ppfx_universes (   fppfx_other)

Definition at line 63 of file plot_xsec_2d.py.

References novadaq::HexUtils.format(), and PandAna.Demos.demo1.range.

63 def get_ppfx_universes(fppfx_other):
64 
65  hl = []
66  # by default we have 100 ppfx universes
67  for i in range(0,100):
68  hl.append(fppfx_other.Get('xs_nominal/hXS2Dsec_ppfx{}'.format(i)))
69 
70  return hl
71 
72 
std::string format(const int32_t &value, const int &ndigits=8)
Definition: HexUtils.cpp:14
def get_ppfx_universes(fppfx_other)
Definition: plot_xsec_2d.py:63
def plot_xsec_2d.process_ckv_calibshape (   isys,
  hidx 
)

Definition at line 130 of file plot_xsec_2d.py.

References abs(), PandAna.Demos.demo1.range, and SetBinContent().

130 def process_ckv_calibshape(isys, hidx):
131 
132  for i in range(1, NbinsX+1):
133  for j in range(1, NbinsY+1):
134  dt_cv = hdt.GetBinContent(i,j)
135 
136  lhup[isys].SetBinContent(i,j,0)
137  lhdw[isys].SetBinContent(i,j,0)
138  lhfe[isys].SetBinContent(i,j,0)
139 
140  sigma = dt_cv - hother[hidx].GetBinContent(i,j)
141 
142  # if sigma > 0:
143  # lhdw[isys].SetBinContent(i,j,abs(sigma))
144  # if sigma <= 0:
145  # lhup[isys].SetBinContent(i,j,abs(sigma))
146  lhdw[isys].SetBinContent(i,j,abs(sigma))
147  lhup[isys].SetBinContent(i,j,abs(sigma))
148 
149  if dt_cv > 0:
150  lhfe[isys].SetBinContent(i,j,sigma/dt_cv)
151  else:
152  lhfe[isys].SetBinContent(i,j,0)
153 
154 
void abs(TH1 *hist)
def process_ckv_calibshape(isys, hidx)
cout<< "--"<< endl;for(Int_t iP=1;iP<=hyz->GetNbinsX();iP++){for(Int_t iC=1;iC<=hyz->GetNbinsY();iC++){if(hyv->GetBinContent(iP, iC)>-999){goal_hyv-> SetBinContent(iP, iC,-(dy[iP-1][iC-1]))
def plot_xsec_2d.process_data ( )

Definition at line 218 of file plot_xsec_2d.py.

References databin(), and PandAna.Demos.demo1.range.

219 
220  for i in range(1, NbinsX+1):
221  for j in range(1, NbinsY+1):
222  kin = hdt.GetYaxis().GetBinCenter(j)
223  binw = hdt.GetYaxis().GetBinWidth(j)
224  cont = hdt.GetBinContent(i,j)
225  errUp = 0
226  errDw = 0
227 
228  for p in range(Nfinal_fe):
229  up = lhup[p].GetBinContent(i,j)
230  dw = lhdw[p].GetBinContent(i,j)
231  errUp += up**2
232  errDw += dw**2
233 
234  errUp = math.sqrt(errUp)
235  errDw = math.sqrt(errDw)
236 
237  if databin(i,j):
238  count[i-1] += 1
239  gdt[i-1].SetPoint(count[i-1], kin, cont)
240  gdt[i-1].SetPointError(count[i-1],binw/2, binw/2, errDw, errUp)
241 
242 
def process_data()
def databin(xbin, ybin)
def plot_xsec_2d.process_focusing (   isys,
  sstart,
  nsys 
)

Definition at line 172 of file plot_xsec_2d.py.

References abs(), cet::sqlite.max(), print, PandAna.Demos.demo1.range, SetBinContent(), and sum.

172 def process_focusing(isys, sstart, nsys):
173 
174  # debug information
175  print('Focusing systematics:')
176  for k in range(nsys):
177  print(syst_names[sstart+k*2], syst_names[sstart+k*2+1])
178 
179  for i in range(1, NbinsX+1):
180  for j in range(1, NbinsY+1):
181  dt_cv = hdt.GetBinContent(i,j)
182 
183  lhup[isys].SetBinContent(i,j,0)
184  lhdw[isys].SetBinContent(i,j,0)
185  lhfe[isys].SetBinContent(i,j,0)
186 
187  larger_sigma = []
188  for k in range(nsys):
189  sigma_dw = dt_cv - hother[sstart+k*2].GetBinContent(i,j)
190  sigma_up = dt_cv - hother[sstart+k*2+1].GetBinContent(i,j)
191  larger_sigma.append(max(abs(sigma_dw), abs(sigma_up)))
192 
193  qsum = math.sqrt(sum([x**2 for x in larger_sigma]))
194  lhup[isys].SetBinContent(i,j,qsum)
195  lhdw[isys].SetBinContent(i,j,qsum)
196 
197  if dt_cv > 0:
198  lhfe[isys].SetBinContent(i,j,qsum/dt_cv)
199  else:
200  lhfe[isys].SetBinContent(i,j,0)
201 
202 
def process_focusing(isys, sstart, nsys)
void abs(TH1 *hist)
bool print
cout<< "--"<< endl;for(Int_t iP=1;iP<=hyz->GetNbinsX();iP++){for(Int_t iC=1;iC<=hyz->GetNbinsY();iC++){if(hyv->GetBinContent(iP, iC)>-999){goal_hyv-> SetBinContent(iP, iC,-(dy[iP-1][iC-1]))
Double_t sum
Definition: plot.C:31
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:68
def plot_xsec_2d.process_light_calibration_mue (   isys,
  hidxdown,
  hidxup 
)

Definition at line 102 of file plot_xsec_2d.py.

References abs(), cet::sqlite.max(), PandAna.Demos.demo1.range, and SetBinContent().

102 def process_light_calibration_mue(isys, hidxdown, hidxup):
103 
104  for i in range(1, NbinsX+1):
105  for j in range(1, NbinsY+1):
106  dt_cv = hdt.GetBinContent(i,j)
107 
108  lhup[isys].SetBinContent(i,j,0)
109  lhdw[isys].SetBinContent(i,j,0)
110  lhfe[isys].SetBinContent(i,j,0)
111 
112  sigma_dw = dt_cv - hother[hidxdown].GetBinContent(i,j)
113  sigma_up = dt_cv - hother[hidxup].GetBinContent(i,j)
114  sigma = max(abs(sigma_dw), abs(sigma_up))
115 
116  if sigma_dw < 0 and sigma_up < 0:
117  lhup[isys].SetBinContent(i,j,sigma)
118  elif sigma_dw > 0 and sigma_up > 0:
119  lhdw[isys].SetBinContent(i,j,sigma)
120  else:
121  lhup[isys].SetBinContent(i,j,abs(sigma_up))
122  lhdw[isys].SetBinContent(i,j,abs(sigma_dw))
123 
124  if dt_cv > 0:
125  lhfe[isys].SetBinContent(i,j,sigma/dt_cv)
126  else:
127  lhfe[isys].SetBinContent(i,j,0)
128 
129 
def process_light_calibration_mue(isys, hidxdown, hidxup)
void abs(TH1 *hist)
cout<< "--"<< endl;for(Int_t iP=1;iP<=hyz->GetNbinsX();iP++){for(Int_t iC=1;iC<=hyz->GetNbinsY();iC++){if(hyv->GetBinContent(iP, iC)>-999){goal_hyv-> SetBinContent(iP, iC,-(dy[iP-1][iC-1]))
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:68
def plot_xsec_2d.process_mc ( )

Definition at line 243 of file plot_xsec_2d.py.

References databin(), PandAna.Demos.demo1.range, and SetBinContent().

244 
245  for i in range(1, NbinsX+1):
246  for j in range(1, NbinsY+1):
247  if databin(i,j):
248  binw = hdt.GetYaxis().GetBinWidth(j)
249  cont = hmc.GetBinContent(i,j)
250  h1MC[i-1].SetBinContent(j,cont)
251  for p in range(Nfinal_fe):
252  ferr = lhfe[p].GetBinContent(i,j)
253  h1fe[i-1][p].SetBinContent(j,ferr)
254 
def databin(xbin, ybin)
cout<< "--"<< endl;for(Int_t iP=1;iP<=hyz->GetNbinsX();iP++){for(Int_t iC=1;iC<=hyz->GetNbinsY();iC++){if(hyv->GetBinContent(iP, iC)>-999){goal_hyv-> SetBinContent(iP, iC,-(dy[iP-1][iC-1]))
def process_mc()
def plot_xsec_2d.process_multiuniverse (   isys,
  huniv 
)

Definition at line 81 of file plot_xsec_2d.py.

References cet::sqlite.max(), PandAna.Demos.demo1.range, SetBinContent(), and up_down_counts_1sigma().

81 def process_multiuniverse(isys, huniv):
82 
83  for i in range(1, NbinsX+1):
84  for j in range(1, NbinsY+1):
85  dt_cv = hdt.GetBinContent(i,j)
86 
87  if dt_cv <= 0:
88  lhup[isys].SetBinContent(i,j,0)
89  lhdw[isys].SetBinContent(i,j,0)
90  lhfe[isys].SetBinContent(i,j,0)
91  else:
92  # accumulate bin content from all ppfx universes
93  binc_all_uni = [h.GetBinContent(i,j) for h in huniv]
94  cup, cdown = up_down_counts_1sigma(binc_all_uni, dt_cv)
95  lhup[isys].SetBinContent(i,j,cup)
96  lhdw[isys].SetBinContent(i,j,cdown)
97  # calculate conservative fractional uncertainty
98  conserv_error = max(cup,cdown)
99  lhfe[isys].SetBinContent(i,j,conserv_error/dt_cv)
100 
101 
cout<< "--"<< endl;for(Int_t iP=1;iP<=hyz->GetNbinsX();iP++){for(Int_t iC=1;iC<=hyz->GetNbinsY();iC++){if(hyv->GetBinContent(iP, iC)>-999){goal_hyv-> SetBinContent(iP, iC,-(dy[iP-1][iC-1]))
def process_multiuniverse(isys, huniv)
Definition: plot_xsec_2d.py:81
T max(sqlite3 *const db, std::string const &table_name, std::string const &column_name)
Definition: statistics.h:68
def up_down_counts_1sigma(count_list, cv)
Definition: plot_xsec_2d.py:73
def plot_xsec_2d.process_normalization (   isys)

Definition at line 155 of file plot_xsec_2d.py.

References PandAna.Demos.demo1.range, and SetBinContent().

156 
157  for i in range(1, NbinsX+1):
158  for j in range(1, NbinsY+1):
159  dt_cv = hdt.GetBinContent(i,j)
160 
161  sigma = dt_cv * fe_norm
162  # lhup[isys].SetBinContent(i,j,dt_cv+sigma)
163  # lhdw[isys].SetBinContent(i,j,dt_cv-sigma)
164  lhup[isys].SetBinContent(i,j,sigma)
165  lhdw[isys].SetBinContent(i,j,sigma)
166  if dt_cv > 0:
167  lhfe[isys].SetBinContent(i,j,fe_norm)
168  else:
169  lhfe[isys].SetBinContent(i,j,0)
170 
171 
def process_normalization(isys)
cout<< "--"<< endl;for(Int_t iP=1;iP<=hyz->GetNbinsX();iP++){for(Int_t iC=1;iC<=hyz->GetNbinsY();iC++){if(hyv->GetBinContent(iP, iC)>-999){goal_hyv-> SetBinContent(iP, iC,-(dy[iP-1][iC-1]))
def plot_xsec_2d.setData1D (   gIn)

Definition at line 260 of file plot_xsec_2d.py.

260 def setData1D(gIn):
261  gIn.SetMarkerColor(kBlack)
262  gIn.SetMarkerStyle(20)
263  gIn.SetLineColor(kBlack)
264  gIn.SetLineWidth(3)
265 
def setData1D(gIn)
def plot_xsec_2d.setMC1D (   hIn)

Definition at line 255 of file plot_xsec_2d.py.

255 def setMC1D(hIn):
256  hIn.SetLineColor(2)
257  hIn.SetLineWidth(2)
258  hIn.SetStats(0)
259 
def setMC1D(hIn)
def plot_xsec_2d.setUncer1D (   hIn,
  ihist 
)

Definition at line 266 of file plot_xsec_2d.py.

266 def setUncer1D(hIn, ihist):
267  hIn.SetLineColor(col_fe[ihist])
268  hIn.SetLineWidth(3)
269  hIn.SetLineStyle(col_st[ihist])
270  hIn.SetStats(0)
271 
272 
def setUncer1D(hIn, ihist)
def plot_xsec_2d.up_down_counts_1sigma (   count_list,
  cv 
)

Definition at line 73 of file plot_xsec_2d.py.

References makeTrainCVSamples.int.

Referenced by process_multiuniverse().

73 def up_down_counts_1sigma(count_list, cv):
74  sorted_list = sorted(count_list)
75  cv_pos = bisect(sorted_list, cv)
76  up_pos = int((len(sorted_list) - cv_pos)*.68 + cv_pos)
77  down_pos = int(cv_pos - cv_pos*.68)
78  return sorted_list[up_pos]-cv, cv-sorted_list[down_pos]
79 
80 
def up_down_counts_1sigma(count_list, cv)
Definition: plot_xsec_2d.py:73

Variable Documentation

plot_xsec_2d.action

Definition at line 282 of file plot_xsec_2d.py.

plot_xsec_2d.args = parser.parse_args()

Definition at line 283 of file plot_xsec_2d.py.

list plot_xsec_2d.cfe = []

Definition at line 417 of file plot_xsec_2d.py.

list plot_xsec_2d.col_fe = [2,4,2,4,8,8,kOrange-3,kOrange-3,1]

Definition at line 41 of file plot_xsec_2d.py.

list plot_xsec_2d.col_st = [1,1,2,2,1,2,1,2,2]

Definition at line 42 of file plot_xsec_2d.py.

list plot_xsec_2d.count = [0 for i in range(Nabins2X)]

Definition at line 375 of file plot_xsec_2d.py.

list plot_xsec_2d.cxs = []

Definition at line 390 of file plot_xsec_2d.py.

plot_xsec_2d.dest_path = os.path.join('output', out_subdir)

Definition at line 385 of file plot_xsec_2d.py.

float plot_xsec_2d.fe_norm = 0.028422526

Definition at line 366 of file plot_xsec_2d.py.

list plot_xsec_2d.fgenie = []

Definition at line 289 of file plot_xsec_2d.py.

plot_xsec_2d.fileOut = TFile('{}/numucc_inc_{}.root'.format(dest_path,out_subdir), 'recreate')

Definition at line 447 of file plot_xsec_2d.py.

Referenced by make_nominal_xs(), and ND_DataMC().

list plot_xsec_2d.final_fe = ["Flux (HP)","GENIE","Light","Calibration","Cherenkov","Calib. Shape","Mu Energy Scale","Normalization","Focusing"]

Definition at line 22 of file plot_xsec_2d.py.

list plot_xsec_2d.final_fe_name = ["Flux_HP","GENIE","Light","Calibration","Cherenkov","CalibShape","MuEScale","Normalization","Focusing"]

Definition at line 23 of file plot_xsec_2d.py.

plot_xsec_2d.fmc = TFile(os.path.join(infp, 'nominal_xsec.root'))

Definition at line 290 of file plot_xsec_2d.py.

plot_xsec_2d.fppfx_other = TFile(os.path.join(infp, 'ppfx_others_v2.root'))

Definition at line 292 of file plot_xsec_2d.py.

list plot_xsec_2d.gdt = [TGraphAsymmErrors() for i in range(Nabins2X)]

Definition at line 374 of file plot_xsec_2d.py.

list plot_xsec_2d.h1fe = [[TH1D('h1fe_{}_{}'.format(i,j),'',20,0.5,2.5) for j in range(Nfinal_fe)] for i in range(Nabins2X)]

Definition at line 380 of file plot_xsec_2d.py.

list plot_xsec_2d.h1MC = [TH1D('h1MC_{}'.format(i),'',20,0.5,2.5) for i in range(Nabins2X)]

Definition at line 379 of file plot_xsec_2d.py.

plot_xsec_2d.hdt = fppfx_other.Get('xs_nominal/hXS2Dsec_cv')

Definition at line 304 of file plot_xsec_2d.py.

plot_xsec_2d.help

Definition at line 282 of file plot_xsec_2d.py.

list plot_xsec_2d.hgenie = []

Definition at line 297 of file plot_xsec_2d.py.

plot_xsec_2d.hmc = fmc.Get('hXS2Dsec_nominal')

Definition at line 302 of file plot_xsec_2d.py.

list plot_xsec_2d.hother = []

Definition at line 299 of file plot_xsec_2d.py.

list plot_xsec_2d.hppfx = []

Definition at line 298 of file plot_xsec_2d.py.

string plot_xsec_2d.infp = './'

Definition at line 13 of file plot_xsec_2d.py.

list plot_xsec_2d.lat = [TLatex() for i in range(Nabins2X)]

Definition at line 392 of file plot_xsec_2d.py.

list plot_xsec_2d.latfe = [TLatex() for i in range(Nabins2X)]

Definition at line 418 of file plot_xsec_2d.py.

list plot_xsec_2d.leg = [TLegend(0.67,0.68,0.87,0.88) for i in range(Nabins2X)]

Definition at line 391 of file plot_xsec_2d.py.

list plot_xsec_2d.leg_angle
Initial value:
1 = [
2  "0.50<Cos#theta_{#mu}<0.56",
3  "0.56<Cos#theta_{#mu}<0.62",
4  "0.62<Cos#theta_{#mu}<0.68",
5  "0.68<Cos#theta_{#mu}<0.74",
6  "0.74<Cos#theta_{#mu}<0.80",
7  "0.80<Cos#theta_{#mu}<0.85",
8  "0.85<Cos#theta_{#mu}<0.88",
9  "0.88<Cos#theta_{#mu}<0.91",
10  "0.91<Cos#theta_{#mu}<0.94",
11  "0.94<Cos#theta_{#mu}<0.96",
12  "0.96<Cos#theta_{#mu}<0.98",
13  "0.98<Cos#theta_{#mu}<0.99",
14  "0.99<Cos#theta_{#mu}<1.00"
15 ]

Definition at line 26 of file plot_xsec_2d.py.

plot_xsec_2d.lfe = []

Definition at line 419 of file plot_xsec_2d.py.

list plot_xsec_2d.lhdw = []

Definition at line 325 of file plot_xsec_2d.py.

list plot_xsec_2d.lhfe = []

Definition at line 326 of file plot_xsec_2d.py.

list plot_xsec_2d.lhup = []

Definition at line 324 of file plot_xsec_2d.py.

plot_xsec_2d.maxY = -1

Definition at line 393 of file plot_xsec_2d.py.

int plot_xsec_2d.Nabins2X = 13

Definition at line 25 of file plot_xsec_2d.py.

plot_xsec_2d.NbinsX = hdt.GetNbinsX()
plot_xsec_2d.NbinsY = hdt.GetNbinsY()
int plot_xsec_2d.Nfinal_fe = 9

Definition at line 21 of file plot_xsec_2d.py.

string plot_xsec_2d.out_subdir = 'real_data'

Definition at line 384 of file plot_xsec_2d.py.

plot_xsec_2d.parser = argparse.ArgumentParser()

Definition at line 281 of file plot_xsec_2d.py.

list plot_xsec_2d.syst_names
Initial value:
1 = [
2  'lightdown','lightup','ckv','calibneg','calibpos','calibshape','MuES_onlyMuKEDw','MuES_onlyMuKEUp',
3  '2kA_Dw','2kA_Up','02mmBeamSpotSize_Dw','02mmBeamSpotSize_Up','1mmBeamShiftX_Dw','1mmBeamShiftX_Up',
4  '1mmBeamShiftY_Dw','1mmBeamShiftY_Up','3mmHorn1X_Dw','3mmHorn1X_Up','3mmHorn1Y_Dw','3mmHorn1Y_Up',
5  '3mmHorn2X_Dw','3mmHorn2X_Up','3mmHorn2Y_Dw','3mmHorn2Y_Up','7mmTargetZ_Dw','7mmTargetZ_Up',
6  'MagneticFieldinDecayPipe_Dw','MagneticFieldinDecayPipe_Up','1mmHornWater_Dw','1mmHornWater_Up'
7 ]

Definition at line 14 of file plot_xsec_2d.py.

string plot_xsec_2d.tit_fe = 'Fractional Uncertainties'

Definition at line 46 of file plot_xsec_2d.py.

string plot_xsec_2d.tit_ke = 'Muon Kinetic Energy (GeV)'

Definition at line 44 of file plot_xsec_2d.py.

string plot_xsec_2d.tit_nu = 'Neutrino Energy (GeV)'

Definition at line 45 of file plot_xsec_2d.py.

string plot_xsec_2d.tit_xs = '#sigma(cm^{2}/nucleon)'

Definition at line 47 of file plot_xsec_2d.py.

plot_xsec_2d.use_data = args.data

Definition at line 284 of file plot_xsec_2d.py.