Skip to content
Snippets Groups Projects
Commit a6750134 authored by Carlos Marcelo Rubin de Celis's avatar Carlos Marcelo Rubin de Celis
Browse files

Upload New File

parent 6a194489
No related branches found
No related tags found
No related merge requests found
%% Cell type:code id: tags:
``` python
import math
#Funcion para calcular una integral tipo gaussiana uniforme.
#a = Altura del pico
#b = Centro de la campana
#sigma = Desviación estándar
#interv = intervalos de división para integración
#li = límite inferior de la curva
#ls = límite superior de la curva
def _gauss_integral(a, b, sigma, interv, li, ls):
incr = (ls-li)/interv
res = 0
for i in range(interv):
r = li+i*incr
r_1 = li+(i+1)*incr
b_m = math.exp(-(r-b)*(r-b)/(2*sigma*sigma))
b_M = math.exp(-(r_1-b)*(r_1-b)/(2*sigma*sigma))
res += (b_M+b_m)*incr/2
return(a*res)
#Variables
inter= 1000000
```
%% Cell type:markdown id: tags:
Centelleador LYSO - Fuente Cobalto 60
%% Cell type:code id: tags:
``` python
#PICO 1 1173 KeV
height = 4118.21
center = 1169.99
sigma = 131.73
li = 1075
ls = 1250
print("Para el PICO 1 1173 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Para el PICO 1 1173 KeV: 670087.9485019445
%% Cell type:code id: tags:
``` python
#PICO 2 1333 KeV
height = 2516.36715-2200
center = 1326.59876
sigma = 128.014016
li = 1260
ls = 1380
print("Para el PICO 2 1333 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Para el PICO 2 1333 KeV: 36573.54376706954
%% Cell type:markdown id: tags:
Centelleador BGO - Fuente Cobalto 60
%% Cell type:code id: tags:
``` python
#Tipo de Centelleador : BGO
#Primer fotopico 1173 KeV
height = 5530.68837-2500
center = 1171.05450
sigma = 85.0379959
li = 1050
ls = 1250
print("Primer fotopico 1173 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Primer fotopico 1173 KeV: 481992.1546793344
%% Cell type:code id: tags:
``` python
#Tipo de Centelleador : BGO
#Segundo Fotopico 1333 KeV
height = 3942.70339-1500
center = 1323.43289
sigma = 66.7403072
li = 1280
ls = 1400
print("Segundo Fotopico 1333 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Segundo Fotopico 1333 KeV: 252038.5177349729
%% Cell type:markdown id: tags:
Centelleador CSI - Fuente Cobalto 60
%% Cell type:code id: tags:
``` python
#Tipo de Centelleador : CSI
#PICO 1173 KeV
height = 1406.14
center = 1179.80
sigma = 66.8607
li = 1140
ls = 1240
print("PICO 1173 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
PICO 1173 KeV: 127306.56774417317
%% Cell type:code id: tags:
``` python
#Tipo de Centelleador : CSI
#PICO 1333 KeV
height = 864.25
center = 1326.22
sigma = 48.70
li = 1260
ls = 1400
print("PICO 1173 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
PICO 1173 KeV: 89481.92612266203
%% Cell type:markdown id: tags:
Centelleador BGO - Fuente Cesio 137
%% Cell type:code id: tags:
``` python
#Tipo de Centelleador : BGO
#Fotopico 662. KeV
height = 14365.0902-2000
center = 658.936978
sigma = 37.2170682
li = 550
ls = 800
print("Fotopico 662. KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Fotopico 662. KeV: 1151470.976408768
%% Cell type:markdown id: tags:
Centelleador LYSO - Fuente Cesio 137
%% Cell type:code id: tags:
``` python
#Tipo de Centelleador : LYSO
height = 3066.55380
center = 656.387332
sigma = 41.0754300
li = 550
ls = 800
print("Tipo de Centelleador : LYSO : "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Tipo de Centelleador : LYSO : 314145.5297374663
%% Cell type:markdown id: tags:
Centelleador CSI - fuente cesio 137
%% Cell type:code id: tags:
``` python
# Tipo de Centelleador : CSI
height = 273.31
center = 654.55
sigma = 27.2305
li = 550
ls = 800
print("Tipo de Centelleador : CSI: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Tipo de Centelleador : CSI: 18654.09886979009
%% Cell type:markdown id: tags:
Centelleador Tipo NaI
%% Cell type:code id: tags:
``` python
#Fuente CESIO 137 fotopico GAMMA 662 KeV
height = 19529.4580-2500
center = 662.255849
sigma = 28.3688220
li = 600
ls = 750
print("Fuente CESIO 137 fotopico GAMMA: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Fuente CESIO 137 fotopico GAMMA: 1192692.9032746102
%% Cell type:code id: tags:
``` python
#Fuente CO 60
#Fotopico 1173 KeV
height = 12861.6378-4000
center = 1165.65474
sigma = 59.9694536
li = 1100
ls = 1250
print("Fotopico 1173 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Fotopico 1173 KeV: 1043568.7670718948
%% Cell type:code id: tags:
``` python
#Fuente CO 60
#Fotopico 1333 KeV
height = 11155.3780
center = 1331.19634
sigma = 59.6862446
li = 1250
ls = 1450
print("Fotopico 1333 KeV: "+str(_gauss_integral(height, center, sigma, inter, li, ls)))
```
%% Output
Fotopico 1333 KeV: 1485177.0230626967
%% Cell type:code id: tags:
``` python
```
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment