real function one_temp_deu(x)
*--------------------------------------------
*  Function to fit the pion mt spectra with two temperatures
*  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,T2,N1,N2,A1,A2,E,m,mt,P,gamma,pi
	integer i

        common /pawpar/ par(4)

        T1  = par(1)	
        N1  = par(2)	
        T2  = par(3)	
        N2  = par(4)	

	m=1.88
	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))
	A2 = N2 / (2.0 * pi * T2 * (m*m + 2.0*m*T2 + 2.0*T2*T2))

	one_temp_deu= A1 * mt * exp(-x/T1) 

      end