newgauss.py 794 Bytes
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
# Versao 0.1 do Gauss usando python
# Usando classes e carregando varios arquivos

from pylab import *
from spectro import *

# Iniciando o espectro de espalhamento vazio
data=spectro_espalhamento(50000,100000,30)


# Definindo os parametros experimentais e de interacao dos ions com a materia
param = dict(dedx=287. , dW2dx=35000 , E0= 100000, Theta_in=0, Theta_out=70, fwhm0=300 )


# Adicionando os sinais correspondentes aos elementos do alvo
# ** Ainda nao considera as massas dos elementos **
  #data.addelements(['distHf.dat','distO.dat'],param,0.9,'gaussiana')


data.addelements2(['profHf.dat','profO.dat'],param,0.9,'gaussiana',0.01)
#pares = array([[1,1],[3,0.5],[1,1]])
#print pares
#x,c = profile(pares, 0.1)
#print x
#print c



# Apresenta o resultado em um grafico
data.plot()