Code to test the likelihood ratio of a three term fit vs a 2 term fit for our eta distributions of pt slices
just run python likelihood_ratio_implementation.py -b to replicate my results
there are lots of paramaters to play around with in the code (like the number of bins)
Also, the likelihood fitting code is a bit finiky. Check the Chi2/ndf or plot to visualy make sure it is fitting as expected.
Read more about the math here:
https://en.wikipedia.org/wiki/Likelihood-ratio_test
https://en.wikipedia.org/wiki/Chi-squared_distribution
http://maxwell.ucsc.edu/~drip/133/ch4.pdf
warning, lots of plots will be produced! Feel free to comment out the draw or save as parts