def do_endcap_lowr9_lownv_fits():
'''
End Cap Low R9 Low NV
'''
systematics = 0.2
#______________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EE_sixie_LowNV_R9Low_0_R9High_0.94_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Endcaps_LowNV_LowR9',
title = 'Endcaps NVtx < 19 R_{9}^{#gamma} < 0.94, MC Truth',
systematics = systematics,
)
fitter.run()
fitters.append(fitter)
eefitter = fitter
#______________________________________________________________________________
## Fix S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EE_sixie_LowNV_R9Low_0_R9High_0.94_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Endcaps_LowNV_LowR9_SfromTB',
title = 'Endcaps NVtx < 19 R_{9}^{#gamma} < 0.94, MC Truth, S from TB',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(eefitter.N.getVal())
fitter.C.setVal(eefitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix N to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EE_sixie_LowNV_R9Low_0_R9High_0.94_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Endcaps_LowNV_LowR9_NfromTB',
title = 'Endcaps NVtx < 19 R_{9}^{#gamma} < 0.94, MC Truth, N from TB',
systematics = systematics,
)
fitter.S.setVal(eefitter.S.getVal())
fitter.N.setVal(_noise_from_tb/_mean_cosh_eta_endcaps)
fitter.C.setVal(eefitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
def do_barrel_highr9_fits():
'''
Barrel High R9
'''
systematics = 0.2
#__________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EB_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Barrel_HighR9',
title = 'Barrel R_{9}^{#gamma} > 0.94, MC Truth',
systematics = systematics,
)
fitter.run()
fitters.append(fitter)
ebfitter = fitter
#__________________________________________________________________________
## Fix S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EB_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Barrel_HighR9_SfromTB',
title = 'Barrel R_{9}^{#gamma} > 0.94, MC Truth, S from TB',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(ebfitter.N.getVal())
fitter.C.setVal(ebfitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#__________________________________________________________________________
## Fix N to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EB_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Barrel_HighR9_NfromTB',
title = 'Barrel R_{9}^{#gamma} > 0.94, MC Truth, N from TB',
systematics = systematics,
)
fitter.S.setVal(ebfitter.S.getVal())
fitter.N.setVal(_noise_from_tb/_mean_cosh_eta_barrel)
fitter.C.setVal(ebfitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
def do_barrel_lownv_fits():
'''
Barrel Low NV
'''
systematics = 0.2
#__________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Barrel_LowNV',
title = 'Barrel NVtx < 19, MC Truth',
systematics = systematics,
)
fitter.run()
fitters.append(fitter)
ebfitter = fitter
#__________________________________________________________________________
## S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Barrel_LowNV_SfromTB',
title = 'Barrel NVtx < 19, MC Truth, S from TB',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(ebfitter.N.getVal())
fitter.C.setVal(ebfitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#__________________________________________________________________________
## N from TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_restruth_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCTruth_Barrel_LowNV_NfromTB',
title = 'Barrel NVtx < 19, MC Truth, N from TB',
systematics = systematics,
)
fitter.S.setVal(ebfitter.S.getVal())
fitter.N.setVal(_noise_from_tb/_mean_cosh_eta_barrel)
fitter.C.setVal(ebfitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
def do_barrel_lownv_fits():
'''
Barrel Low NV
'''
systematics = 0.2
#______________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EB_sixie_LowNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Barrel_LowNV',
title = 'Barrel NVtx < 18 , Data Fit',
systematics = systematics,
)
fitter.run()
fitters.append(fitter)
ebfitter = fitter
#______________________________________________________________________________
## Fix S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EB_sixie_LowNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Barrel_LowNV_SfromTB',
title = 'Barrel NVtx < 18 , Data Fit, S from TB',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(ebfitter.N.getVal())
fitter.C.setVal(ebfitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EB_sixie_LowNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Barrel_LowNV_NfromMC',
title = 'Barrel NVtx < 18 , Data Fit, N from MC',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(_noise_from_mc_barrel_lownv)
fitter.C.setVal(ebfitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix S to TB and N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EB_sixie_LowNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Barrel_LowNV_SfromTB_NfromMC',
title = 'Barrel NVtx < 18 , Data Fit, S from TB, N from MC',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(_noise_from_mc_barrel_lownv)
fitter.C.setVal(ebfitter.C.getVal())
fitter.S.setConstant()
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
def do_endcap_highr9_fits():
'''
Endcaps High R9
'''
systematics = 0.2
#______________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EE_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Endcaps_HighR9',
title = 'Endcaps R_{9}^{#gamma} > 0.94, Data Fit',
systematics = systematics,
)
fitter.run()
fitters.append(fitter)
eefitter = fitter
#______________________________________________________________________________
## Fix S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EE_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Endcaps_HighR9_SfromTB',
title = 'Endcaps R_{9}^{#gamma} > 0.94, Data Fit, S from TB',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(eefitter.N.getVal())
fitter.C.setVal(eefitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EE_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Endcaps_HighR9_NfromMC',
title = 'Endcaps R_{9}^{#gamma} > 0.94, Data Fit, N from MC',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(_noise_from_mc_endcaps_highr9)
fitter.C.setVal(eefitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix S to TB and N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resdata_EE_sixie_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_DataFit_Endcaps_HighR9_SfromTB_NfromMC',
title = 'Endcaps R_{9}^{#gamma} > 0.94, Data Fit, S from TB, N from MC',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(_noise_from_mc_endcaps_highr9)
fitter.C.setVal(eefitter.C.getVal())
fitter.S.setConstant()
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
def do_barrel_highr9_lownv_fits():
'''
Barrel High R9 Low NV
'''
systematics = 0.2
#__________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Barrel_LowNV_HighR9',
title = 'Barrel NVtx < 19 R_{9}^{#gamma} > 0.94, MC Fit',
systematics = systematics,
yrange = (-1, 6),
)
fitter.run()
fitters.append(fitter)
ebfitter = fitter
#__________________________________________________________________________
## Fix S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Barrel_LowNV_HighR9_SfromTB',
title = 'Barrel NVtx < 19 R_{9}^{#gamma} > 0.94, MC Fit, S from TB',
systematics = systematics,
yrange = (-1 ,6),
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(ebfitter.N.getVal())
fitter.C.setVal(ebfitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#__________________________________________________________________________
## Fix N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Barrel_LowNV_HighR9_NfromMC',
title = 'Barrel NVtx < 19 R_{9}^{#gamma} > 0.94, MC Fit, N from MC',
systematics = systematics,
yrange = (-1 , 6),
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(_noise_from_mc_barrel_highr9_lownv)
fitter.C.setVal(ebfitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
#__________________________________________________________________________
## Fix S to TB and N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EB_sixie_LowNV_R9Low_0.94_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Barrel_LowNV_HighR9_SfromTB_NfromMC',
title = 'Barrel NVtx < 19 R_{9}^{#gamma} > 0.94, MC Fit, S from TB, N from MC',
systematics = systematics,
yrange = (-1 , 6),
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_barrel)
fitter.N.setVal(_noise_from_mc_barrel_highr9_lownv)
fitter.C.setVal(ebfitter.C.getVal())
fitter.S.setConstant()
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
def do_endcap_highnv_fits():
'''
End Cap High NV
'''
systematics = 0.2
#______________________________________________________________________________
## Float everything
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EE_sixie_HighNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Endcaps_HighNV',
title = 'Endcaps NVtx > 18 , MC Fit',
systematics = systematics,
)
fitter.run()
fitters.append(fitter)
eefitter = fitter
#______________________________________________________________________________
## Fix S to TB
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EE_sixie_HighNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Endcaps_HighNV_SfromTB',
title = 'Endcaps NVtx > 18 , MC Fit, S from TB',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(eefitter.N.getVal())
fitter.C.setVal(eefitter.C.getVal())
fitter.S.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EE_sixie_HighNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Endcaps_HighNV_NfromMC',
title = 'Endcaps NVtx > 18 , MC Fit, N from MC',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(_noise_from_mc_endcaps_highnv)
fitter.C.setVal(eefitter.C.getVal())
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)
#______________________________________________________________________________
## Fix S to TB and N to MC
fitter = Fitter(
filename = _filename,
graphname = 'regressions_resmc_EE_sixie_HighNV_R9Low_0_R9High_999_sixie',
name = 'PhotonResolutionVsEt_MCFit_Endcaps_HighNV_SfromTB_NfromMC',
title = 'Endcaps NVtx > 18 , MC Fit, S from TB, N from MC',
systematics = systematics,
)
fitter.S.setVal(_stochastic_from_tb/_mean_sqrt_cosh_eta_endcaps)
fitter.N.setVal(_noise_from_mc_endcaps_highnv)
fitter.C.setVal(eefitter.C.getVal())
fitter.S.setConstant()
fitter.N.setConstant()
fitter.run()
fitters.append(fitter)