def METHOD_NAME(obj_pos, psr_dir, T_obj):
    """
    ss_obj_shapiro_delay(obj_pos, psr_dir, T_obj)
    returns Shapiro delay in seconds for a solar system object.
    Inputs:
      obj_pos : position vector from Earth to SS object, with Units
      psr_dir : unit vector in direction of pulsar
      T_obj : mass of object in seconds (GM/c^3)
    """
    # TODO: numpy.sum currently loses units in some cases...
    r = (numpy.sqrt(numpy.sum(obj_pos**2, axis=1))) * obj_pos.unit
    rcostheta = numpy.sum(obj_pos * psr_dir, axis=1)
    # This is the 2nd to last term from Eqn 4.6 in Backer &
    # Hellings, ARAA, 1986 with gamma = 1 (as defined by GR).  We
    # have the opposite sign of the cos(theta) term, since our
    # object position vector is defined from the observatory to
    # the solar system object, rather than from the object to the
    # pulsar (as described after Eqn 4.3 in the paper).
    # See also https://en.wikipedia.org/wiki/Shapiro_time_delay
    # where \Delta t = \frac{2GM}{c^3}\log(1-\vec{R}\cdot\vec{x})
    return -2.0 * T_obj * numpy.log((r - rcostheta) / const.au).value