ENH: Adding smoothing of hit and total images (-l option) and quadratic constraint on prior