Integration wasn't properly normalized for non-SC luminosity functions, I had to change the integral to be in linear space and fix the PDF normalization in Luminosity
Increased the number of luminosity bins (more precise)
Unrelated: added str methods to Evolution.py
To make sure that everything was working properly, I simulated a SC population with FIRESONG.py (histogram) and compared this to the PDF I get from FluxPDF using a SC (blue) or a LG with narrow width (orange) and everything looks good. Plot attached below. This should close https://github.com/icecube/FIRESONG/issues/30 and https://github.com/icecube/FIRESONG/issues/15 as we can just suggest users use Theo's FluxPDF for high densities
A few things:
To make sure that everything was working properly, I simulated a SC population with FIRESONG.py (histogram) and compared this to the PDF I get from FluxPDF using a SC (blue) or a LG with narrow width (orange) and everything looks good. Plot attached below. This should close https://github.com/icecube/FIRESONG/issues/30 and https://github.com/icecube/FIRESONG/issues/15 as we can just suggest users use Theo's FluxPDF for high densities