ERRORNB = []
EYHIGH = GRFD.GetEYhigh()
EYLOW = GRFD.GetEYlow()
YVAL = GRFD.GetY()
#-1 because D0 has also bin [1-2]
for i in range(histo_.GetNbinsX()-1):
ERRORNB.append([0, 0, EYLOW[i+2], EYHIGH[i+2], YVAL[i+2]])
DICTNB["feeddown_NB"] = ERRORNB
#-1 because D0 has also bin [1-2]
errs = Errors(histo_.GetNbinsX()-1)
errs.read(ERROR_FILESD0[mb], DICTNB)
ERRS_D0.append(errs)
HISTOS_DSOVERD0 = divide_by_eachother(HISTOS_DS, HISTOS_D0, [2.27, 3.89], [2, 4, 6, 8, 12, 24])
for mb, _ in enumerate(HISTOS_DSOVERD0):
tot_Ds_over_D0 = calc_systematic_mesonratio(ERRS_DS[mb], ERRS_D0[mb], \
HISTOS_DSOVERD0[mb].GetNbinsX())
tot_Ds_over_D0_FD = calc_systematic_mesonratio(ERRS_DS[mb], ERRS_D0[mb], \
HISTOS_DSOVERD0[mb].GetNbinsX(), True)
tot_Ds_over_D0_WOFD = calc_systematic_mesonratio(ERRS_DS[mb], ERRS_D0[mb], \
HISTOS_DSOVERD0[mb].GetNbinsX(), False)
ERRS_GR_DIVD0_TOT.append(Errors.make_root_asymm(HISTOS_DSOVERD0[mb], \
tot_Ds_over_D0, const_x_err=0.3))
ERRS_GR_DIVD0_FD.append(Errors.make_root_asymm(HISTOS_DSOVERD0[mb], \
tot_Ds_over_D0_FD, const_x_err=0.3))
ERRS_GR_DIVD0_WOFD.append(Errors.make_root_asymm(HISTOS_DSOVERD0[mb], \
tot_Ds_over_D0_WOFD, const_x_err=0.3))
ERRS_GR_DIVD0_TOT[mb].SetName("gr_TotSyst_%d" % mb)
GRFD = extract_histo_or_error("D0pp", "MBvspt_ntrkl", mb, YEAR_NUMBER, \
"gFcConservative", PATHD0)
ERRORNB = []
EYHIGH = GRFD.GetEYhigh()
EYLOW = GRFD.GetEYlow()
YVAL = GRFD.GetY()
for i in range(histo_.GetNbinsX()):
ERRORNB.append([0, 0, EYLOW[i + 1], EYHIGH[i + 1], YVAL[i + 1]])
DICTNB["feeddown_NB"] = ERRORNB
errs = Errors(histo_.GetNbinsX())
errs.read(ERROR_FILESD0[mb], DICTNB)
ERRS_D0.append(errs)
HISTOS_LCOVERD0 = divide_by_eachother(HISTOS_LC, HISTOS_D0, [6.23, 3.89])
for mb, _ in enumerate(HISTOS_LCOVERD0):
tot_Lc_over_D0 = calc_systematic_mesonratio(ERRS_LC[mb], ERRS_D0[mb], \
HISTOS_LCOVERD0[mb].GetNbinsX())
tot_Lc_over_D0_FD = calc_systematic_mesonratio(ERRS_LC[mb], ERRS_D0[mb], \
HISTOS_LCOVERD0[mb].GetNbinsX(), True)
tot_Lc_over_D0_WOFD = calc_systematic_mesonratio(ERRS_LC[mb], ERRS_D0[mb], \
HISTOS_LCOVERD0[mb].GetNbinsX(), False)
ERRS_GR_DIVD0_TOT.append(Errors.make_root_asymm(HISTOS_LCOVERD0[mb], \
tot_Lc_over_D0, const_x_err=0.3))
ERRS_GR_DIVD0_FD.append(Errors.make_root_asymm(HISTOS_LCOVERD0[mb], \
tot_Lc_over_D0_FD, const_x_err=0.3))
ERRS_GR_DIVD0_WOFD.append(Errors.make_root_asymm(HISTOS_LCOVERD0[mb], \
tot_Lc_over_D0_WOFD, const_x_err=0.3))
ERRS_GR_DIVD0_TOT[mb].SetName("gr_TotSyst_%d" % mb)
GRFD = extract_histo(CASED0, ANA_MB, mb, "gFcConservative", PATH_MIX)
ERRORNB = []
EYHIGH = GRFD.GetEYhigh()
EYLOW = GRFD.GetEYlow()
YVAL = GRFD.GetY()
#-1 because D0 has also bin [1-2]
for i in range(histo_.GetNbinsX() - 1):
ERRORNB.append([0, 0, EYLOW[i + 2], EYHIGH[i + 2], YVAL[i + 2]])
DICTEXTRA["feeddown_NB"] = ERRORNB
#-1 because D0 has also bin [1-2]
errs = Errors(histo_.GetNbinsX() - 1)
errs.read(ERROR_FILESD0[mb], DICTEXTRA)
ERRS_D0.append(errs)
HISTOS_DSOVERD0 = divide_by_eachother(HISTOS_DS, HISTOS_D0, [BRDS, BRD0],
[2, 4, 6, 8, 12, 24])
for mb, _ in enumerate(HISTOS_DSOVERD0):
tot_Ds_over_D0 = calc_systematic_mesonratio(ERRS_DS[mb], ERRS_D0[mb], \
HISTOS_DSOVERD0[mb].GetNbinsX())
tot_Ds_over_D0_FD = calc_systematic_mesonratio(ERRS_DS[mb], ERRS_D0[mb], \
HISTOS_DSOVERD0[mb].GetNbinsX(), True)
tot_Ds_over_D0_WOFD = calc_systematic_mesonratio(ERRS_DS[mb], ERRS_D0[mb], \
HISTOS_DSOVERD0[mb].GetNbinsX(), False)
ERRS_GR_DIVD0_TOT.append(Errors.make_root_asymm(HISTOS_DSOVERD0[mb], \
tot_Ds_over_D0, const_x_err=0.3))
ERRS_GR_DIVD0_FD.append(Errors.make_root_asymm(HISTOS_DSOVERD0[mb], \
tot_Ds_over_D0_FD, const_x_err=0.3))
ERRS_GR_DIVD0_WOFD.append(Errors.make_root_asymm(HISTOS_DSOVERD0[mb], \
tot_Ds_over_D0_WOFD, const_x_err=0.3))
ERRS_GR_DIVD0_TOT[mb].SetName("gr_TotSyst_%d" % mb)
# For now strip the first bin of Lc by hand. Introduce funtion in the future
# This is necessary as Lc spectrum goes down to 1 - 2 whereas D0 starts at 2 - 4
HISTOS_LC_STRIP = []
for mb in range(3):
histo_lc_new = HISTOS_D0[mb].Clone(
f"{HISTOS_LC[mb].GetName()}_lc_over_d0")
histo_lc_new.Reset("ICEMS")
histo_lc_old = HISTOS_LC[mb]
for i in range(1, histo_lc_new.GetNbinsX() + 1):
histo_lc_new.SetBinContent(i,
histo_lc_old.GetBinContent(i + 1))
histo_lc_new.SetBinError(i, histo_lc_old.GetBinError(i + 1))
HISTOS_LC_STRIP.append(histo_lc_new)
HISTOS_LC_OVER_D0 = divide_by_eachother(HISTOS_LC_STRIP, HISTOS_D0,
[1, 1])
SAVE_PATH = f"plot_LcpK0spp_over_D0pp_year_{YEAR_NUMBER}_{hn}.eps"
results(HISTOS_LC_OVER_D0,
None,
"",
LEGEND_TITLES,
"#it{p}_{T} (GeV/#it{c})",
"d#it{N}/(d#it{p}_{T})|_{|y|<0.5} (GeV^{-1} #it{c})",
SAVE_PATH,
False,
colors=COLORS,
y_range=[0, 1],
log_y=False)
##################
# COMPARE TO STD #