Subversion Repositories svnkaklik

Rev

Rev 945 | 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
 
946 kaklik 93
set output "utlum0.png"
94
plot "./utlumy.txt" using (f($1)):2 with points title "0 mm pozadi", "" using (f($1)):4 with points title "0 mm zarice"
95
 
96
set output "utlum2.png"
97
plot "./utlumy.txt" using (f($1)):7 with points title "2 mm pozadi", "" using (f($1)):9 with points title "2 mm zarice"
98
 
99
 
100
set output "utlum35.png"
101
plot "./utlumy.txt" using (f($1)):12 with points title "3,5 mm pozadi", "" using (f($1)):14 with points title "3,5 mm zarice"
102
 
103
 
945 kaklik 104
print u1
105
print u2
106
print u3