Begun coding in lisanti et al distribution function - needs debugging