-
Notifications
You must be signed in to change notification settings - Fork 0
/
20020921_Polarimetry_1.gp
62 lines (48 loc) · 6.04 KB
/
20020921_Polarimetry_1.gp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
# Copyright 2002-2003 D. C. Hatton
#
# D. C. Hatton can be contacted on [email protected]
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; version 2 of the License.
#
# This program is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program (in the file copying.programs;) if not,
# write to the Free Software Foundation, Inc., 59 Temple Place, Suite
# 330, Boston, MA 02111-1307 USA
# This plots the first data set from 21st September 2002
FIT_LIMIT = 1e-20
c1=1.0
c2=1.0
A=2e-5
f(E, M) = M < 0.5 ? 264199216.0723*c1*(1.0-A)/E**2.0 : M < 1.5 ? 264199216.0723*c2*(1.0+A)/E**2.0 : M < 2.5 ? 264199216.0723*c1*(1.0+A)/E**2.0 : 264199216.0723*c2*(1.0-A)/E**2.0
fit f(x,y) '20020921_Polarimetry_1.dat' using 3:($1+$2):($8/$5):(sqrt(($9/$5)**2.0+($8*$6/$5**2.0)**2.0)) via c1, c2, A
set pointsize 0.01
set style data xyerrorbars
set key below
set size 1, 2
set terminal pslatex
set color
set xlabel "$E_b/\\ut{eV}$"
set ylabel "$f_d/\\ut{Hz}$"
set output '20020921_Polarimetry_1_2.tex'
plot '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 < 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization up, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 >= 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization up, detector 2", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 < 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization down, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 >= 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization down, detector 2", 2000.0*f(x, 2.0) title "Magnetization up, detector 1", 2000.0*f(x, 3.0) title "Magnetization up, detector 2", 2000.0*f(x, 0.0) title "Magnetization down, detector 1", 2000.0*f(x, 1.0) title "Magnetization down, detector 2"
set output '20020921_Polarimetry_1_4.tex'
plot '20020921_Polarimetry_1.dat' using 3:(1 >= 0 && $2 < 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 >= 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 2", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 < 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 >= 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 2", 4000.0*f(x, 2.0) title "Magnetization up, detector 1", 4000.0*f(x, 3.0) title "Magnetization up, detector 2", 4000.0*f(x, 0.0) title "Magnetization down, detector 1", 4000.0*f(x, 1.0) title "Magnetization down, detector 2"
set output '20020921_Polarimetry_1_6.tex'
plot '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 < 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 >= 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 2", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 < 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 >= 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 2", 6000.0*f(x, 2.0) title "Magnetization up, detector 1", 6000.0*f(x, 3.0) title "Magnetization up, detector 2", 6000.0*f(x, 0.0) title "Magnetization down, detector 1", 6000.0*f(x, 1.0) title "Magnetization down, detector 2"
set terminal png
set color
set xlabel "Eb/eV"
set ylabel "fd/Hz"
set output '20020921_Polarimetry_1_2.png'
plot '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 < 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization up, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 >= 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization up, detector 2", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 < 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization down, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 >= 1.5 && $5 < 3000 ? $8 : 1/0):4:9 title "Magnetization down, detector 2", 2000.0*f(x, 2.0) title "Magnetization up, detector 1", 2000.0*f(x, 3.0) title "Magnetization up, detector 2", 2000.0*f(x, 0.0) title "Magnetization down, detector 1", 2000.0*f(x, 1.0) title "Magnetization down, detector 2"
set output '20020921_Polarimetry_1_4.png'
plot '20020921_Polarimetry_1.dat' using 3:(1 >= 0 && $2 < 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 >= 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 2", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 < 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 >= 1.5 && $5 >= 3000 && $5 < 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 2", 4000.0*f(x, 2.0) title "Magnetization up, detector 1", 4000.0*f(x, 3.0) title "Magnetization up, detector 2", 4000.0*f(x, 0.0) title "Magnetization down, detector 1", 4000.0*f(x, 1.0) title "Magnetization down, detector 2"
set output '20020921_Polarimetry_1_6.png'
plot '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 < 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 >= 0 && $2 >= 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization up, detector 2", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 < 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 1", '20020921_Polarimetry_1.dat' using 3:($1 < 0 && $2 >= 1.5 && $5 >= 5000 ? $8 : 1/0):4:9 title "Magnetization down, detector 2", 6000.0*f(x, 2.0) title "Magnetization up, detector 1", 6000.0*f(x, 3.0) title "Magnetization up, detector 2", 6000.0*f(x, 0.0) title "Magnetization down, detector 1", 6000.0*f(x, 1.0) title "Magnetization down, detector 2"