ENH: add score_hessian_factor and hessian

adds analytical hessian but doesn't use it yet (`hessian_1`)