v = std.vector(Double)(2) v[0] = 1.0; v[1] = 0.3; # Fill histogram with the result of kernel density estimation truepdf.project(true_hist) true_hist.Write() poly.project(poly_hist) poly_hist.Write() kde.project(kernel_hist) kernel_hist.Write() # Fill 1D slice in x projection, at the point y=v[1] truepdf.slice(v, 0, true_kpi) true_kpi.Write() poly.slice(v, 0, poly_kpi) poly_kpi.Write() kde.slice(v, 0, kernel_kpi) kernel_kpi.Write() # Fill 1D slice in y projection, at the point x=v[0] truepdf.slice(v, 1, true_pipi) true_pipi.Write() poly.slice(v, 1, poly_pipi) poly_pipi.Write() kde.slice(v, 1, kernel_pipi) kernel_pipi.Write() # Plot 2D and 1D histograms gStyle.SetOptStat(0)