def _g_theta(theta,p):
    """
    :param theta: dimensionless electron temperature
    :param p: electron power law slope
    :return: correction term for power law electron distribution
    """
    aa = (6.0 + 15.0 * theta) / (4.0 + 5.0 * theta)
    gamma_m = 1e0 + aa * theta
    gtheta = ((p-1.0)*(1e0+aa*theta)/((p-1.0)*gamma_m - p+2.0))*(gamma_m/(3.0*theta))**(p-1.0)
    return gtheta