sim_semicontinuous_data function

Generates data from a two part distribution with a point mass at zero and heterogeneous treatment effects