Replicate histogram by translating R code to Python

Question

The following plot and its code were generated in R. How can I replicate this quality of a histogram in Python code using scipy.stats?

hi