Dear David,
Thank you so much for your response!
By <a^2> * \pi^2/(48 \cosh(\pi R)):
<a^2> is the average of the Hecke eigenvalues^2. For example, for the first cusp form, which has spectral parameter R=9.534, it would be the average of: (1)^2, (-1.068)^2, (-0.456)^2, (0.141)^2, etc. I think this can be derived by considering the triple integral (Eisenstein * cusp^2) and using the Rankin-Selberg unfolding trick on the Eisenstein, and then taking the residue as s approaches 1. If you know of a better way to get this numerically (I am only interested in SL2Z, which I think is level 1), I would be very interested!
If it is not too much trouble, I would be extremely grateful if you could send me the Petersson norms based on the formula above using your internal data. The first e.g. 30 or so cusp forms would be fine, if that is easy for you to compute. Thank you so much for your help!
Best
Nathan