SAVE_PATH, False, colors=COLORS)
#############################################################################
##################### Plot spectra mult / spectra MB ########################
#############################################################################
#Divide by MB
HISTOS_DIV = divide_all_by_first_multovermb(HISTOS)
#Remove MB one
HISTOS_DIVMB = HISTOS_DIV[1:]
ERRS_GR_DIV_TOT = []
ERRS_GR_DIV_WOFD = []
ERRS_GR_DIV_FD = []
for mb, _ in enumerate(HISTOS_DIVMB):
tot_mult_over_MB = calc_systematic_multovermb(ERRS[mb+1], ERRS[0], HISTOS[0].GetNbinsX())
tot_mult_over_MB_FD = calc_systematic_multovermb(ERRS[mb+1], ERRS[0], \
HISTOS[0].GetNbinsX(), True)
tot_mult_over_MB_WOFD = calc_systematic_multovermb(ERRS[mb+1], ERRS[0], \
HISTOS[0].GetNbinsX(), False)
ERRS_GR_DIV_TOT.append(Errors.make_root_asymm(HISTOS_DIVMB[mb], tot_mult_over_MB, \
const_x_err=0.3))
ERRS_GR_DIV_WOFD.append(Errors.make_root_asymm(HISTOS_DIVMB[mb], tot_mult_over_MB_WOFD, \
const_x_err=0.3))
ERRS_GR_DIV_FD.append(Errors.make_root_asymm(HISTOS_DIVMB[mb], tot_mult_over_MB_FD, \
const_x_err=0.3))
ERRS_GR_DIV_TOT[mb].SetName("gr_TotSyst_%d" % mb)
ERRS_GR_DIV_FD[mb].SetName("gr_FDSyst_%d" % mb)
ERRS_GR_DIV_WOFD[mb].SetName("gr_TotSyst_woFD_%d" % mb)
ERRORNB = []
EYHIGH = GRFD.GetEYhigh()
EYLOW = GRFD.GetEYlow()
YVAL = GRFD.GetY()
for i in range(HISTOS[0].GetNbinsX()):
if i >= 5:
continue
ERRORNB.append([0, 0, EYLOW[i + 1], EYHIGH[i + 1], YVAL[i + 1]])
DICTNB["feeddown_NB"] = ERRORNB
ERRFOR19 = Errors(HISTOS_DIVMB[0].GetNbinsX())
NAME19 = "data/errors/LcpKpipp_full/MBvspt_ntrkl/errors_histoSigmaCorr_0_for19.yaml"
ERRFOR19.read(NAME19, DICTNB)
ERRSFOR19.append(ERRFOR19)
TOT_MULT_OVER_MB = calc_systematic_multovermb(ERRS[1], ERRSFOR19[0], \
HISTOS_DIVMB[0].GetNbinsX())
TOT_MULT_OVER_MB_FD = calc_systematic_multovermb(ERRS[1], ERRSFOR19[0], \
HISTOS_DIVMB[0].GetNbinsX(), True)
TOT_MULT_OVER_MB_WOFD = calc_systematic_multovermb(ERRS[1], ERRSFOR19[0], \
HISTOS_DIVMB[0].GetNbinsX(), False)
ERRS_GR_DIV_TOT.append(Errors.make_root_asymm(HISTOS_DIVMB[0], TOT_MULT_OVER_MB, \
const_x_err=0.3))
ERRS_GR_DIV_WOFD.append(Errors.make_root_asymm(HISTOS_DIVMB[0], TOT_MULT_OVER_MB_WOFD, \
const_x_err=0.3))
ERRS_GR_DIV_FD.append(Errors.make_root_asymm(HISTOS_DIVMB[0], TOT_MULT_OVER_MB_FD, \
const_x_err=0.3))
ERRS_GR_DIV_TOT[0].SetName("gr_TotSyst_%d" % 0)
ERRS_GR_DIV_FD[0].SetName("gr_FDSyst_%d" % 0)
ERRS_GR_DIV_WOFD[0].SetName("gr_TotSyst_woFD_%d" % 0)