computeNormScales Namespace Reference

Functions
def	GetPercentile (histo, X)
def	GetMedian (histo)
def	GetMinimum (histo)

Variables
list	hierarchies = ["nh","ih"]
list	plotvars = ["cp_75thpercentile","mh_minimum"]
list	tweaks = ["fhc_numu","fhc_numubar","fhc_nue","rhc_numu","rhc_numu","rhc_numubar","rhc_nuebar"]
list	amounts_varied = [0,10,25,50,75,90,100]
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");
list	medians = [[],[],[]]
list	histos = [[],[],[]]
list	amounts_varied_good_files = []
	var = tweak+str(amount_varied)
string	file_prefix = var+"_"
	t_file = ROOT.TFile("/lbne/data/users/lblpwg_tools/FastMC_Data/outputs/ljf26/Sensitivity_Plots/"+file_prefix+"_histos.root");
list	metrics = [[],[],[]]
string	nomfhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_200kA_LBNEFD_fastmc.root"
string	nomrhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/-200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_-200kA_LBNEFD_fastmc.root"
string	varfhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_200kA_LBNEFD_fastmc.root"
string	varrhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/-200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_-200kA_LBNEFD_fastmc.root"
int	power = 120
float	antinufrac = 0.5
	c1 = ROOT.TCanvas(file_prefix);
list	t_graphs = []
int	mgraph = 0
list	line_styles = [1,2,3]
list	line_colors = [1,2,4]
	temp = ROOT.TGraph(len(amounts_varied_good_files),array.array('d',amounts_varied_good_files),array.array('d',medians[i]));
	leg = ROOT.TLegend(0.35,0.2,0.65,0.5);
	m_graph = ROOT.TGraph(len(amounts_varied_good_files),array.array('d',amounts_varied_good_files),array.array('d',metrics[i]));
	ofile = ROOT.TFile(file_prefix+".root","RECREATE")

Function Documentation

def computeNormScales.GetMedian

(

histo

)

Definition at line 20 of file computeNormScales.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 computeNormScales.GetMinimum

(

histo

)

Definition at line 35 of file computeNormScales.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 computeNormScales.GetPercentile	(		histo,
			X
	)

Definition at line 3 of file computeNormScales.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]

Variable Documentation

list computeNormScales.amounts_varied = [0,10,25,50,75,90,100]

Definition at line 52 of file computeNormScales.py.

list computeNormScales.amounts_varied_good_files = []

Definition at line 83 of file computeNormScales.py.

float computeNormScales.antinufrac = 0.5

Definition at line 142 of file computeNormScales.py.

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

Definition at line 59 of file computeNormScales.py.

list computeNormScales.baseline_medians = []

Definition at line 57 of file computeNormScales.py.

computeNormScales.c1 = ROOT.TCanvas(file_prefix);

Definition at line 165 of file computeNormScales.py.

string computeNormScales.file_prefix = var+"_"

Definition at line 88 of file computeNormScales.py.

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

Definition at line 42 of file computeNormScales.py.

list computeNormScales.histos = [[],[],[]]

Definition at line 81 of file computeNormScales.py.

computeNormScales.leg = ROOT.TLegend(0.35,0.2,0.65,0.5);

Definition at line 210 of file computeNormScales.py.

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

Definition at line 171 of file computeNormScales.py.

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

Definition at line 170 of file computeNormScales.py.

computeNormScales.m_graph = ROOT.TGraph(len(amounts_varied_good_files),array.array('d',amounts_varied_good_files),array.array('d',metrics[i]));

Definition at line 221 of file computeNormScales.py.

list computeNormScales.medians = [[],[],[]]

Definition at line 79 of file computeNormScales.py.

list computeNormScales.metrics = [[],[],[]]

Definition at line 132 of file computeNormScales.py.

int computeNormScales.mgraph = 0

Definition at line 168 of file computeNormScales.py.

computeNormScales.mixing_param = plotvar[0:2]

Definition at line 58 of file computeNormScales.py.

string computeNormScales.nomfhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_200kA_LBNEFD_fastmc.root"

Definition at line 135 of file computeNormScales.py.

string computeNormScales.nomrhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/-200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_-200kA_LBNEFD_fastmc.root"

Definition at line 136 of file computeNormScales.py.

computeNormScales.ofile = ROOT.TFile(file_prefix+".root","RECREATE")

Definition at line 246 of file computeNormScales.py.

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

Definition at line 43 of file computeNormScales.py.

int computeNormScales.power = 120

Definition at line 141 of file computeNormScales.py.

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

Definition at line 91 of file computeNormScales.py.

list computeNormScales.t_graphs = []

Definition at line 167 of file computeNormScales.py.

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

Definition at line 175 of file computeNormScales.py.

list computeNormScales.tweaks = ["fhc_numu","fhc_numubar","fhc_nue","rhc_numu","rhc_numu","rhc_numubar","rhc_nuebar"]

Definition at line 45 of file computeNormScales.py.

computeNormScales.var = tweak+str(amount_varied)

Definition at line 87 of file computeNormScales.py.

string computeNormScales.varfhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_200kA_LBNEFD_fastmc.root"

Definition at line 138 of file computeNormScales.py.

string computeNormScales.varrhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r2p4/QGSP_BERT/Nominal/-200kA/flux/histos_g4lbne_v3r2p4_QGSP_BERT_Nominal_-200kA_LBNEFD_fastmc.root"

Definition at line 139 of file computeNormScales.py.