Subversion Repositories svnkaklik

Rev

Rev 776 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log

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