Working version of constrained PMF