[main.py]: Tested modifications to simulate_cmd.py

- no major changes in this version, just a first commit
  for branch 'binaries'
[myanalysis] Implement Transformation Method

- Previously, simuate_cmd.py used rejection sampling to produce
  mock samples.

- Have now implemented transformation method by calculating cumulative
  mass distribution, and then finding the inverse from spline interpolation.

- To use rejection, need to specify rejection_sampling = 1 at beginning of module.

- Tested transformation method and added additional diagnostic
  plots: both methods result in the same LFs, so both methods are
  working appropriately.

- Transformation method fails when the cumulative F^-1(M) inverse
  function becomes almost vertical. This happens only when dN/dM's
  slope becomes very very close to zero, as in the case of lognormal
  with a small sigma. For all other cases, when F(M) grows gradually,
  no problem. Hence transform method can work for all cases except
  when testing small sigma_c or fitting likelihood with small sigma_c.