// param[0]*inv_p_gev + param[1]*ty^2 inv_p_gev + param[2]*tx^2 inv_p_gev +
// param[3]*tx inv_p_gev pol_qop_gev + param[4]*inv_p_gev^3 + param[5]*tx^3
// pol_qop_gev + param[6]*tx^2 inv_p_gev^2 + param[7]*tx inv_p_gev^2 pol_qop_gev
// + param[8]*inv_p_gev^4 + param[9]*ty^2 tx^2 inv_p_gev + param[10]*ty^2 tx
// inv_p_gev pol_qop_gev + param[11]*ty^2 inv_p_gev^3 + param[12]*tx^4 inv_p_gev
static constexpr std::array<float, 13> momentumWindowParamsRef{
    4013.5392952527636f, 6713.931835651405f,  3698.1175175030185f,
    -4205.958156564871f, 1563.16527675433f,   5436.545339478851f,
    7829.470965204402f,  -5946.6465946298f,   1408.4004460016097f,
    103314.80145944108f, -93751.44471183348f, 26807.566464613057f,
    -25400.206891250582f};