exampleSensitivityCalculation.py
Go to the documentation of this file.
1 import OptimizationUtils
2 
3 
4 nominal_fhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r3p2/QGSP_BERT/Nominal/200kA/flux/histos_g4lbne_v3r3p2_QGSP_BERT_Nominal_200kA_LBNEFD_fastmc.root"
5 nominal_rhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r3p2/QGSP_BERT/Nominal/-200kA/flux/histos_g4lbne_v3r3p2_QGSP_BERT_Nominal_-200kA_LBNEFD_fastmc.root"
6 nominal_sensitivity = OptimizationUtils.ComputeMean75PercentSensitivityAndError(nominal_fhcfile,nominal_rhcfile, 120.0, 120.0,0.5, "cp")
7 print "Nominal Sensitivity: ",nominal_sensitivity[1][0]
8 
9 varied_fhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r3p6/QGSP_BERT/CP_run5_9116_80GeV/200kA/flux/histos_g4lbne_v3r3p6_QGSP_BERT_CP_run5_9116_80GeV_200kA_LBNEFD_fastmc.root"
10 varied_rhcfile = "/lbne/data/users/ljf26/fluxfiles/g4lbne/v3r3p6/QGSP_BERT/CP_run5_9116_80GeV/-200kA/flux/histos_g4lbne_v3r3p6_QGSP_BERT_CP_run5_9116_80GeV_-200kA_LBNEFD_fastmc.root"
11 
12 varied_sensitivity = OptimizationUtils.ComputeMean75PercentSensitivityAndError(varied_fhcfile,varied_rhcfile, 80.0, 80.0,0.5, "cp")
13 print "Varied Sensitivity: ",varied_sensitivity[1][0]
14 
15 
OptimizationUtils.ComputeMean75PercentSensitivityAndError
def ComputeMean75PercentSensitivityAndError(varfhcfile, varrhcfile, fhcvarenergy=120.0, rhcvarenergy=120.0, antinufrac=0.5, quantity="cp")
Definition: OptimizationUtils.py:81
Generated by   doxygen 1.8.11