changed to expbeta*PTx*n*rho.e0 from expbeta*PTx*PYs*rho.e0