Subversion Repositories svnkaklik

Rev

Rev 776 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log

Rev Author Line No. Line
767 kaklik 1 
set terminal png
2 
 
3 
set key off
4 
set grid
5 
 
774 kaklik 6 
f(x)=k*x+q
945 kaklik 7 
g(x)=p/x+t
767 kaklik 8 
 
945 kaklik 9 
#energeticka kalibrace
10 
 
776 kaklik 11 
fit f(x) "kalibrace.txt" using 1:2 via k,q
774 kaklik 12 
set output "kalibrace.png"
13 
set xlabel "Číslo kanálu"
14 
set ylabel "Energie [keV]"
15 
plot "kalibrace.txt" using 1:2 with points, f(x)
16 
 
945 kaklik 17 
#set xrange [0:2000]
18 
 
19 
# energeticke rozliseni
776 kaklik 20 
fit g(x) "kalibrace.txt" using ($2):(($3)/($2)) via p,t
945 kaklik 21 
set output "kalibrace_sigma.png"
776 kaklik 22 
set xlabel "Energie [keV]"
23 
set ylabel "S [-]"
24 
plot "kalibrace.txt" using 2:($3/$2) with points, g(x)
25 
 
774 kaklik 26 
set xrange [0:1800]
767 kaklik 27 
 
774 kaklik 28 
set output "pozadi.png"
29 
set xlabel "Energie [keV]"
30 
set ylabel "Intenzita [Imp]"
945 kaklik 31 
plot "./data2/pozadi.txt" using (f($1)):2 with points
767 kaklik 32 
 
945 kaklik 33 
 
34 
set key top
35 
 
776 kaklik 36 
set xrange [0:300]
774 kaklik 37 
set output "Am241.png"
945 kaklik 38 
plot "./data2/Am.txt" using (f($1)):2 with points title "Am"
774 kaklik 39 
 
776 kaklik 40 
set xrange [0:800]
945 kaklik 41 
set yrange [0:8000]
774 kaklik 42 
set output "Ba133.png"
945 kaklik 43 
plot "./data2/Ba.txt" using (f($1)):2 with points title "Ba"
774 kaklik 44 
 
776 kaklik 45 
set xrange [0:1800]
945 kaklik 46 
set yrange [0:500]
774 kaklik 47 
set output "Co60.png"
945 kaklik 48 
plot "./data2/Co.txt" using (f($1)):2 with points title "Co"
774 kaklik 49 
 
945 kaklik 50 
set yrange [0:3000]
774 kaklik 51 
set output "Cs137.png"
945 kaklik 52 
plot "./data2/Cs.txt" using (f($1)):2 with points title "Cs"
774 kaklik 53 
 
945 kaklik 54 
set key off
55 
set yrange [0:600]
774 kaklik 56 
set output "neznamy.png"
945 kaklik 57 
plot "./data2/neznamy.txt" using (f($1)):2 with points title "Uknown"
774 kaklik 58 
 
776 kaklik 59 
set key under
945 kaklik 60 
#set autoscale x
61 
set autoscale y
62 
set xrange [0:4]
63 
#set logscale y
774 kaklik 64 
 
945 kaklik 65 
c1(x)=I1*exp(u1*x)
66 
c2(x)=I2*exp(u2*x)
67 
c3(x)=I3*exp(u3*x)
774 kaklik 68 
 
945 kaklik 69 
I1=10000
70 
I2=10000
71 
I3=10000
776 kaklik 72 
 
945 kaklik 73 
fit [0:2] c1(x) "./stineni.txt" using ($1):2 via I1,u1
74 
fit [0:2] c2(x) "./stineni.txt" using ($1):3 via I2,u2
75 
fit [0:2] c3(x) "./stineni.txt" using ($1):4 via I3,u3
76 
 
77 
set output "stineni.png"
78 
plot "./stineni.txt" using ($1):2:5 with errorbars title "1332,492 keV", c1(x), "" using ($1):3:6 with errorbars title "1173,228 keV",c2(x)
79 
 
80 
 
81 
set xrange [0:4]
82 
set yrange [150000:200000]
83 
 
84 
set output "stineni2.png"
85 
plot "./stineni.txt" using ($1):4:7 with errorbars title "661,657 keV", c3(x)
86 
 
87 
set autoscale y
88 
#set autoscale x
89 
set xrange [100:1500]
90 
set output "utlum.png"
91 
plot "./utlumy.txt" using (f($1)):5 with points title "0 mm", "" using (f($1)):10 with points title "2 mm", "" using (f($1)):15 with points title "3,5 mm"
92 
 
93 
print u1
94 
print u2
95 
print u3