Our function does not coincide with the one defined by McCarthy, but coincides with $(-1)$-times of it.
$\rho_r^{s_j}$ should read $\rho_r^{|s_j|}$.
My research interests include arithmetic aspects of differential equations and special functions.