Functions
def	myfunc (x, par)

def	GetPercentile (histo, X)

def	GetMedian (histo)

def	GetMinimum (histo)

def	GetFractionAtMoreThanXSigma (histo, X)

Variables
bool	use_realistic_powers = True

bool	use_fixed_scale = True

list	tempgraphs = []

list	amounts_varied = [0,50,90,110,200,500]

list	modes = ["fhc","rhc"]

list	nus = ["numu","numubar","nuenuebar"]

list	vars = ["0to0p5GeVPlus10per","0p5to1GeVPlus10per","1to2GeVPlus10per","2to3GeVPlus10per","3to4GeVPlus10per","4to5GeVPlus10per","5to6GeVPlus10per","6to7GeVPlus10per","7to8GeVPlus10per","8to9GeVPlus10per","9to10GeVPlus10per","10to15GeVPlus10per","15to20GeVPlus10per","20to120GeVPlus10per"]

list	bin_centers = [0.25,0.75,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,12.5,17.5,70.0]

	bin_centers_array = array.array('d',bin_centers)

list	hierarchies = ["nh","ih"]

list	plotvars = ["cp_75thpercentile","mh_minimum"]

	outfile = ROOT.TFile("/lbne/data/users/ljf26/fluxfiles/sensitivity_comps_"+str(amount_varied)+"_improved.root","RECREATE");

list	medians = [[],[],[]]

list	baseline_medians = []

	mixing_param = plotvar[0:2]

	baseline_file = ROOT.TFile("/lbne/data/users/lblpwg_tools/FastMC_Data/outputs/ljf26/Sensitivity_Plots/ProtonP120GeV_"+hierarchy+"_"+mixing_param+"_histos.root");

	var = vars[i]+str(amount_varied)

string	file_prefix = mode+nu+var+"_"

	t_file = ROOT.TFile("/lbne/data/users/lblpwg_tools/FastMC_Data/outputs/ljf26/Sensitivity_Plots/"+file_prefix+"_histos.root");

	c1 = ROOT.TCanvas(file_prefix);

list	t_graphs = []

list	line_styles = [1,2,3]

list	line_colors = [1,2,4]

	temp = ROOT.TGraph(len(bin_centers_array),bin_centers_array,array.array('d',medians[i]));

string	nu_string = "#nu_{#mu}"

	leg = ROOT.TLegend(0.45,0.3,0.85,0.55);

list	sensitivity_changes = []

list	flux_changes = []

	variation = amounts_varied[variation_iter]

	temp_file = ROOT.TFile("/lbne/data/users/ljf26/fluxfiles/sensitivity_comps_"+str(variation)+"_improved.root");

	mygraph = temp_file.Get(plotvar+"_"+mode+"_"+nu+"_"+hierarchy)

list	temps = []

list	lines = []

list	tempfuncs = []

int	xlow = -20

int	xhigh = 1000

float	maximum = 0.6

float	minimum = -0.6

	myrange = ROOT.TMath.Max(ROOT.TMath.Abs(maximum),ROOT.TMath.Abs(minimum));

	title = ROOT.gPad.GetPrimitive("title");

int	iter_100 = -1

int	iter_110 = -1

tuple	slope = (sensitivity_changes[bin_iter][iter_110]-sensitivity_changes[bin_iter][iter_100])/(flux_changes[bin_iter][iter_110]-flux_changes[bin_iter][iter_100]);

	y_int = sensitivity_changes[bin_iter][iter_110]-slope*flux_changes[bin_iter][iter_110];

	ylow_line = slope*xlow+y_int

	yhigh_line = slope*xhigh+y_int

	xlow_line = xlow

	xhigh_line = xhigh

string	scale_string = "fixedScale"

Function Documentation

def drawSensitivitySummaries_simplevars_improved.GetFractionAtMoreThanXSigma	(	histo,
		X
	)

Definition at line 49 of file drawSensitivitySummaries_simplevars_improved.py.

 def GetFractionAtMoreThanXSigma(histo, X):
     n = histo.GetNbinsX();
 
     greaterThanX = 0.0;
     total = 0.0;
 
     for i in range(0,n):
         total = total + float(histo.GetBinWidth(i+1));
         if(histo.GetBinContent(i+1)>=X):
             greaterThanX = greaterThanX + float(histo.GetBinWidth(i+1));
 
 
     return greaterThanX/total
 

def drawSensitivitySummaries_simplevars_improved.GetMedian ( histo )

Definition at line 30 of file drawSensitivitySummaries_simplevars_improved.py.

 def GetMedian(histo):
     n = histo.GetNbinsX();
     x = []
     y = []
     
     for i in range(0,n):
         x.append(histo.GetBinCenter(i+1));
         y.append(histo.GetBinContent(i+1));
         y_array = array.array('d',y)
     
     return ROOT.TMath.Median(n,y_array);
 

def drawSensitivitySummaries_simplevars_improved.GetMinimum ( histo )

Definition at line 42 of file drawSensitivitySummaries_simplevars_improved.py.

 def GetMinimum(histo):
     min = 99999
     for i in range(0,histo.GetNbinsX()):
         if histo.GetBinContent(i+1) < min:
             min = histo.GetBinContent(i+1)
     return min
 

def drawSensitivitySummaries_simplevars_improved.GetPercentile	(	histo,
		X
	)

Definition at line 13 of file drawSensitivitySummaries_simplevars_improved.py.

 def GetPercentile(histo,X):
     n = histo.GetNbinsX();
     x = []
     y = []
     
     for i in range(0,n):
         x.append(histo.GetBinCenter(i+1));
         y.append(histo.GetBinContent(i+1));
     y_array = array.array('d',y)
 
     percentiles = array.array('d',[0.0]);
     probs = array.array('d',[1-float(X)/100.0])
 
     ROOT.TMath.Quantiles(n,1,y_array,percentiles,probs,False);
     return percentiles[0]
 
 

def drawSensitivitySummaries_simplevars_improved.myfunc	(	x,
		par
	)

Definition at line 10 of file drawSensitivitySummaries_simplevars_improved.py.

 def myfunc(x, par):
     return tempgraphs[int(par[0])].Eval(x[0]);
 

Variable Documentation

list drawSensitivitySummaries_simplevars_improved.amounts_varied = [0,50,90,110,200,500]

Definition at line 63 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.baseline_file = ROOT.TFile("/lbne/data/users/lblpwg_tools/FastMC_Data/outputs/ljf26/Sensitivity_Plots/ProtonP120GeV_"+hierarchy+"_"+mixing_param+"_histos.root");

Definition at line 88 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.baseline_medians = []

Definition at line 86 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.bin_centers = [0.25,0.75,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,12.5,17.5,70.0]

Definition at line 69 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.bin_centers_array = array.array('d',bin_centers)

Definition at line 71 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.c1 = ROOT.TCanvas(file_prefix);

Definition at line 144 of file drawSensitivitySummaries_simplevars_improved.py.

string drawSensitivitySummaries_simplevars_improved.file_prefix = mode+nu+var+"_"

Definition at line 116 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.flux_changes = []

Definition at line 235 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.hierarchies = ["nh","ih"]

Definition at line 72 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.iter_100 = -1

Definition at line 301 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.iter_110 = -1

Definition at line 302 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.leg = ROOT.TLegend(0.45,0.3,0.85,0.55);

Definition at line 191 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.line_colors = [1,2,4]

Definition at line 149 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.line_styles = [1,2,3]

Definition at line 148 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.lines = []

Definition at line 255 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.maximum = 0.6

Definition at line 265 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.medians = [[],[],[]]

Definition at line 84 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.minimum = -0.6

Definition at line 266 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.mixing_param = plotvar[0:2]

Definition at line 87 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.modes = ["fhc","rhc"]

Definition at line 65 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.mygraph = temp_file.Get(plotvar+"_"+mode+"_"+nu+"_"+hierarchy)

Definition at line 248 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.myrange = ROOT.TMath.Max(ROOT.TMath.Abs(maximum),ROOT.TMath.Abs(minimum));

Definition at line 270 of file drawSensitivitySummaries_simplevars_improved.py.

string drawSensitivitySummaries_simplevars_improved.nu_string = "#nu_{#mu}"

Definition at line 155 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.nus = ["numu","numubar","nuenuebar"]

Definition at line 66 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.outfile = ROOT.TFile("/lbne/data/users/ljf26/fluxfiles/sensitivity_comps_"+str(amount_varied)+"_improved.root","RECREATE");

Definition at line 78 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.plotvars = ["cp_75thpercentile","mh_minimum"]

Definition at line 74 of file drawSensitivitySummaries_simplevars_improved.py.

string drawSensitivitySummaries_simplevars_improved.scale_string = "fixedScale"

Definition at line 335 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.sensitivity_changes = []

Definition at line 234 of file drawSensitivitySummaries_simplevars_improved.py.

tuple drawSensitivitySummaries_simplevars_improved.slope = (sensitivity_changes[bin_iter][iter_110]-sensitivity_changes[bin_iter][iter_100])/(flux_changes[bin_iter][iter_110]-flux_changes[bin_iter][iter_100]);

Definition at line 311 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.t_file = ROOT.TFile("/lbne/data/users/lblpwg_tools/FastMC_Data/outputs/ljf26/Sensitivity_Plots/"+file_prefix+"_histos.root");

Definition at line 118 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.t_graphs = []

Definition at line 146 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.temp = ROOT.TGraph(len(bin_centers_array),bin_centers_array,array.array('d',medians[i]));

Definition at line 153 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.temp_file = ROOT.TFile("/lbne/data/users/ljf26/fluxfiles/sensitivity_comps_"+str(variation)+"_improved.root");

Definition at line 247 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.tempfuncs = []

Definition at line 256 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.tempgraphs = []

Definition at line 7 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.temps = []

Definition at line 254 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.title = ROOT.gPad.GetPrimitive("title");

Definition at line 291 of file drawSensitivitySummaries_simplevars_improved.py.

bool drawSensitivitySummaries_simplevars_improved.use_fixed_scale = True

Definition at line 5 of file drawSensitivitySummaries_simplevars_improved.py.

bool drawSensitivitySummaries_simplevars_improved.use_realistic_powers = True

Definition at line 3 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.var = vars[i]+str(amount_varied)

Definition at line 113 of file drawSensitivitySummaries_simplevars_improved.py.

drawSensitivitySummaries_simplevars_improved.variation = amounts_varied[variation_iter]

Definition at line 241 of file drawSensitivitySummaries_simplevars_improved.py.

list drawSensitivitySummaries_simplevars_improved.vars = ["0to0p5GeVPlus10per","0p5to1GeVPlus10per","1to2GeVPlus10per","2to3GeVPlus10per","3to4GeVPlus10per","4to5GeVPlus10per","5to6GeVPlus10per","6to7GeVPlus10per","7to8GeVPlus10per","8to9GeVPlus10per","9to10GeVPlus10per","10to15GeVPlus10per","15to20GeVPlus10per","20to120GeVPlus10per"]

Definition at line 68 of file drawSensitivitySummaries_simplevars_improved.py.

int drawSensitivitySummaries_simplevars_improved.xhigh = 1000

Definition at line 262 of file drawSensitivitySummaries_simplevars_improved.py.

tuple drawSensitivitySummaries_simplevars_improved.xhigh_line = xhigh

Definition at line 316 of file drawSensitivitySummaries_simplevars_improved.py.

int drawSensitivitySummaries_simplevars_improved.xlow = -20

Definition at line 261 of file drawSensitivitySummaries_simplevars_improved.py.

tuple drawSensitivitySummaries_simplevars_improved.xlow_line = xlow

Definition at line 315 of file drawSensitivitySummaries_simplevars_improved.py.