単一エネルギーガンマ線平均自由行程

私は支配的な相互作用が光電吸収である単一エネルギーガンマビームに関するコードを書いています、mu = 2 cm -1、50000個の乱数を生成し、相互作用の深さをサンプリングする必要があります（私がやったかどうかはわかりませんか否か）。 平均自由行程= mu-1であることは知っていますが、シミュレーションとmuから平均自由行程を見つけて比較する必要があります。

import random
import matplotlib.pyplot as plt
import numpy as np
mu=(2)
random.seed=()
data = np.random.randn(50000)*10
bins = np.arange(data.min(), data.max()+1e-8, 0.1)
meanfreepath = 1/mu
print(meanfreepath)
plt.hist(data, bins=bins)
plt.show()