#!/usr/bin/env python
"""

Plotter den ikke-relativistiske Breit-Wigner fordelingen.

Vi antar m_0 = 91.19 GeV/c^2 og at \Gamma/c^2= 2.49 GeV/c^2
som er massen og vidden til Z-bosonet.

@author Are Raklev
"""

from numpy import *
from matplotlib.pyplot import *

# Definer Breit-Wigner fordelingen
def BW( m ):
    m0 = 91.19 # GeV/c^2
    Gamma = 2.49 # GeV/c^2
    return Gamma/2./pi/((m-m0)**2+(Gamma/2.)**2)

# Lager x-verdier
m = linspace( 60, 120, 1e3 )

# Lager foerst en figur
figure()
plot( m, BW( m ))

# Setter tittel
title('Breit-Wigner fordeling for Z-massen')

# Tekst langs x-aksen
xlabel('$m$ [GeV/c$^2$]')

# Tekst langs y-aksen
ylabel('$|\psi (m)|^2$ [c$^2$/GeV]')

# Lagre som .eps fil
savefig('bw.eps')

# Viser plottet
show()