real function one_temp_pro(x) *-------------------------------------------- * Function to fit the proton mt spectra with one temperature * written by Daniel Cebra 14-Jan-01 * Following a formula provided by Jenn Klay * Modified 19-Apr-02 to work for deuteron spectra *----------------------------------------------- implicit none real T1,N1,A1,E,m,mt,P,gamma,pi integer i common /pawpar/ par(2) T1 = par(1) N1 = par(2) m = 0.93826 pi = 3.1415926 mt = x + m c E = mt c P = Sqrt(E**2-m**2) A1 = N1 / (2.0 * pi * T1 * (m*m + 2.0*m*T1 + 2.0*T1*T1)) one_temp_pro= A1 * mt * exp(-x/T1) end