rewrite the functions in MeV and improve the plotting (logscale,lineshape)