Add code for comparison with 1lep study

This is the code trying to make a comprison with
http://hep.baylor.edu/elog/benwu/1088