Let R be a root system of type BC in a=Rr of general
positive multiplicity. We introduce certain canonical weight function on
Rr which in the case of symmetric domains corresponds to the
integral kernel of the Berezin transform. We compute its spherical transform
and prove certain Bernstein-Sato type formula. This generalizes earlier work of
Unterberger-Upmeier, van Dijk-Pevsner, Neretin and the author. Associated to
the weight functions there are Heckman-Opdam orthogonal polynomials of Jacobi
type on the compact torus, after a change of variables they form an orthogonal
system on the non-compact space a. We consider their spherical
transform and prove that they are the Macdonald-Koornwinder polynomials
multiplied by the spherical transform of the canonical weight function. For
rank one case this was proved earlier by Koornwinder