print("Plotting B0 tagged dataset...")
    
datasets[0].plotOn(frame, R.RooFit.Binning(nbins), R.RooFit.MarkerColor(R.kRed), R.RooFit.Name("data B0 tag"))
pdfs[0].plotOn(frame, R.RooFit.ProjWData(R.RooArgSet(cond_vars),datasets[0]), R.RooFit.LineColor(R.kRed), R.RooFit.Name("pdf B0 tag"))
    
print("Plotting B0bar tagged dataset...")
datasets[1].plotOn(frame, R.RooFit.Binning(nbins), R.RooFit.MarkerColor(R.kBlue), R.RooFit.Name("data B0 tag"))
pdfs[1].plotOn(frame, R.RooFit.ProjWData(R.RooArgSet(cond_vars),datasets[1]), R.RooFit.LineColor(R.kBlue), R.RooFit.Name("pdf B0 tag"))


import fit_utils_2 as fu

fu.plot_sc(
    pdfs=[reso_ws.pdf("pdf_full_td_btag"), reso_ws.pdf("pdf_full_td_bartag")],
    datasets=[data_set_td_btag, data_set_td_bartag],
    var=reso_ws.var("deltat"),
    cond_vars=[reso_ws.var("deltaterr"), reso_ws.var("mod_mbc"), reso_ws.var("deltae"), reso_ws.var("csobdtprime"), reso_ws.var("r")],
    plotpath="sc.png",
)