Stephan/gpthist
1
0
Fork 0
mirror of https://codeberg.org/ET-Kiel/gpthist.git synced 2026-05-01 13:34:22 +02:00
Code Issues Projects Releases Packages Wiki Activity

Compare commits

base: Stephan:d1ff49ab3e10767603af2c08e54d3480faeb5c64
Branches Tags
Stephan:master
Stephan:main
..
compare: Stephan:b9fca1e33fe9d7cea504a91217039eef10197a43
Branches Tags
Stephan:master
Stephan:main

3 commits
d1ff49ab3e ... b9fca1e33f

Author SHA1 Message Date
Stephan I. Böttcher
b9fca1e33f Merge branch 'master' of codeberg.org:ET-Kiel/gpthist 2024-03-01 15:03:27 +01:00
Stephan I. Böttcher
df6ccab1a4 bg models w/ gauss domain 2024-03-01 15:01:38 +01:00
Stephan I. Böttcher
7b4acd9da1 landau: add bg decay 2023-12-16 23:36:29 +01:00
1 changed files with 6 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

12
landau.gpt
View file

@ -2,8 +2,8 @@
landau(l)=sqrt(exp(-l-exp(-l))/2/pi)
mips(a,e,s, x)=a*landau((x-e)/s)
pk(a,e,s, x) = a*exp(-((x-e)/s)**2/2)
compton(a,e,s, x) = a*(1+erf(-(x-e)/s))/2
bg(a,e,s, x) = a*exp((e-x)/s)
compton(a,e,s, x) = a*(1-erf((x-e)/s))/2
bg(a,e,s,ss, x) = exp((e-x)/s)*pk(a,e,ss, x)
keV75 = 74.969 # X-ray
keV73 = 72.805 # X-ray
@ -16,13 +16,13 @@ keV858 = 857.6427 # Comptom edge
Bi_X1(a,e,s ,x) = pk(a*35.7,e*keV75,s, x) + pk(a*23.7*1.11, e*keV73,s, x)
Bi_X2(a,e,s ,x) = pk(a*4.31*1.02,e*keV84,s, x) + pk(a*8.27,e*keV85,s, x) + pk(a*3.02/1.10,e*keV87,s, x)
Bi_C1(a,e,s,b, x) = bg(a,e*keV393,b, x)*(1+erf(-(x-e*keV393)/s))
Bi_C2(a,e,s,b, x) = bg(a,e*keV858,b, x)*(1+erf(-(x-e*keV858)/s))
Bi_C1(a,e,s,b,bb, x) = bg(a,e*keV393,b,bb, x)*(1-erf((x-e*keV393)/s))
Bi_C2(a,e,s,b,bb, x) = bg(a,e*keV858,b,bb, x)*(1-erf((x-e*keV858)/s))
ax21 = (75./85.)**3.5
Bi_spec(x) = Bi_X1(ax1,e,s ,x) + Bi_X2(ax1*ax21,e,s ,x) \
+ Bi_C1(ac1,e,sc1,bc1, x) + Bi_C2(ac2,e,sc2,bc2, x) \
+ bg(ba,e*keV858,bs, x)
+ Bi_C1(ac1,e,sc1,bc1,bbc1, x) + Bi_C2(ac2,e,sc2,bc2,bbc2, x) \
+ bg(ba,e*keV858,bs,bbs, x)
Bi_G3(a,e,s, x) = pk(6.87*a, 1.770228*e, s, x)
Bi_G12(a,e,s, x) = pk(72.9*a, 1.633354*e, s, x)