code: def _jmomentsurfaceMCIntegrand(vz,vR,vT,R,z,df,sigmaR1,gamma,sigmaz1,mvT,n,m,o):
      
      return df(R,vR*sigmaR1,vT*sigmaR1*gamma,z,vz*sigmaz1,use_physical=False,
                func=(lambda x,y,z: x**n*y**m*z**o))\
                *numpy.exp(vR**2./2.+(vT-mvT)**2./2.+vz**2./2.)
docs: Internal function that is the integrand for the vmomentsurface mass integration