Not sure what "psi_pq,sigma are the shapelet eigenfunctions" means - aren't psi_pq the eigenfunctions, and sigma the size parameter (variance of the gaussian part)?