Neutrina (Irvine-Michigan-Brookhaven experiment)¶
In [21]:
import numpy as np
import matplotlib.pyplot as plt
from math import factorial
In [22]:
events=np.arange(0,9,1)
n_intervals=[1042,860,307,78,15,3,0,0,1]
In [23]:
#funkce počítající předpověď Poissonova rozdělení
def Poisson(k,nu):
return nu**k*np.exp(-nu)/factorial(k)
#vektorizace funkce Poisson aby mohla pracovat najednou s polem
vPoisson=np.vectorize(Poisson)
In [24]:
nu=np.sum(events*n_intervals)/np.sum(n_intervals)
print(f"střední počet detekovaných neutrin v jednom intervalu {nu}")
střední počet detekovaných neutrin v jednom intervalu 0.7766695576756288
In [31]:
plt.bar(events,n_intervals,width=1,label='measured data')
plt.bar(events,np.sum(n_intervals)*vPoisson(events,nu),width=0.5,label='Possion prediction')
plt.legend()
plt.xlabel("number of events")
plt.ylabel("number of intervals")
plt.show()
In [ ]: