Distribution Fitting Python SciPy

I want to fit the following distributions to the Factor data (column C) in the csv file and generate plots as attached.
Normal, LogNormal (2 parameter), Weibull, Johnson SB
Provide 4 plots, parameters of the distributions and the source code.