poly_kpi = TH1F("poly_kpi", "M^{2}(K_{S}#pi) slice", 100, phsp.lowerLimit(0), phsp.upperLimit(0) ) kernel_kpi = TH1F("kernel_kpi", "M^{2}(K_{S}#pi) slice", 100, phsp.lowerLimit(0), phsp.upperLimit(0) ) true_pipi = TH1F("true_pipi", "M^{2}(#pi#pi) slice", 100, phsp.lowerLimit(1), phsp.upperLimit(1) ) poly_pipi = TH1F("poly_pipi", "M^{2}(#pi#pi) slice", 100, phsp.lowerLimit(1), phsp.upperLimit(1) ) kernel_pipi = TH1F("kernel_pipi", "M^{2}(#pi#pi) slice", 100, phsp.lowerLimit(1), phsp.upperLimit(1) ) # Coordinates in the phase space where we will create the 1D slices. from ROOT import std, Double 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()