-
Notifications
You must be signed in to change notification settings - Fork 0
/
DVMXC.jl
129 lines (113 loc) · 4.67 KB
/
DVMXC.jl
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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
using DelimitedFiles
using Plots
using LaTeXStrings
z = [.20,.25,.30,.35,.40,.45,.50,.55,.60,.65,.70,.75,.85]
xexp = [.007,.015,.025,.035,.05,.08,.12,.16,.29]
DIS = readdlm("data/DISFrac.dat")
Had = readdlm("data/HadFrac.dat")
DISa = readdlm("data/DISfractions.txt")
Hada = readdlm("data/HADfractions.txt")
DIS_SIDIS = zeros((9,6))
DIS_rho = zeros((9,6))
DIS_phi = zeros((9,6))
DIS_SIDISa = zeros((9,6))
DIS_rhoa = zeros((9,6))
DIS_phia = zeros((9,6))
Had_SIDIS = zeros((9,6,12))
Had_rho = zeros((9,6,12))
Had_phi = zeros((9,6,12))
Had_SIDISa = zeros((9,6,12))
Had_rhoa = zeros((9,6,12))
Had_phia = zeros((9,6,12))
for i in 1:54
DIS_SIDIS[Int.(DIS[i,1]),Int.(DIS[i,2])] = DIS[i,3]
DIS_rho[Int.(DIS[i,1]),Int.(DIS[i,2])] = DIS[i,4]
DIS_phi[Int.(DIS[i,1]),Int.(DIS[i,2])] = DIS[i,5]
end
for i in 1:45
DIS_SIDISa[Int.(DISa[i,1]+1),Int.(DISa[i,2]+1)] = DISa[i,3]
DIS_rhoa[Int.(DISa[i,1]+1),Int.(DISa[i,2]+1)] = DISa[i,4]
DIS_phia[Int.(DISa[i,1]+1),Int.(DISa[i,2]+1)] = DISa[i,5]
end
for i in 1:648
Had_SIDIS[Int.(Had[i,1]),Int.(Had[i,2]),Int.(Had[i,3])] = Had[i,4]
Had_rho[Int.(Had[i,1]),Int.(Had[i,2]),Int.(Had[i,3])] = Had[i,5]
Had_phi[Int.(Had[i,1]),Int.(Had[i,2]),Int.(Had[i,3])] = Had[i,6]
end
for i in 1:540
Had_SIDISa[Int.(Hada[i,1]+1),Int.(Hada[i,2]+1),Int.(Hada[i,3]+1)] = Hada[i,4]
Had_rhoa[Int.(Hada[i,1]+1),Int.(Hada[i,2]+1),Int.(Hada[i,3]+1)] = Hada[i,5]
Had_phia[Int.(Hada[i,1]+1),Int.(Hada[i,2]+1),Int.(Hada[i,3]+1)] = Hada[i,6]
end
zmid = zeros(12)
for i in 1:12
zmid[i] = (z[i+1]+z[i])/2
end
DIS_s = zeros((9,6))
DIS_r = zeros((9,6))
DIS_p = zeros((9,6))
Had_s = zeros(12,45)
Had_r = zeros(12,45)
Had_p = zeros(12,45)
for i in 1:9
for j in 1:5
global DIS_s, DIS_r, DIS_p
global Had_s, Had_r, Had_p
DIS_s[i,j] = (DIS_SIDIS[i,j] > 0 && DIS_SIDISa[i,j] > 0) ? DIS_SIDIS[i,j]/DIS_SIDISa[i,j] : 0
DIS_r[i,j] = (DIS_rho[i,j] > 0 && DIS_rhoa[i,j] > 0) ? DIS_rho[i,j]/DIS_rhoa[i,j] : 0
DIS_p[i,j] = (DIS_phi[i,j] > 0 && DIS_phia[i,j] > 0) ? DIS_phi[i,j]/DIS_phia[i,j] : 0
for k in 1:12
Had_s[k,i+(j-1)*9] = (Had_SIDIS[i,j,k] > 0 && Had_SIDISa[i,j,k] > 0) ? Had_SIDIS[i,j,k]/Had_SIDISa[i,j,k] : 0
Had_r[k,i+(j-1)*9] = (Had_rho[i,j,k] > 0 && Had_rhoa[i,j,k] > 0) ? Had_rho[i,j,k]/Had_rhoa[i,j,k] : 0
Had_p[k,i+(j-1)*9] = (Had_phi[i,j,k] > 0 && Had_phia[i,j,k] > 0) ? Had_phi[i,j,k]/Had_phia[i,j,k] : 0
end
end
end
plot(xexp,DIS_s,layout=(2,3),seriestype=:scatter,
xscale = :log10,
xlims = (0.01,1),
size = (1600,900),
dpi = 150,
ylims = (0.,2.),
# ylims = (0.75,1.25),
legend=false)
savefig("plots/DVMXC_DIS_SIDIS.png")
plot(xexp,DIS_r,layout=(2,3),seriestype=:scatter,
xscale = :log10,
xlims = (0.01,1),
size = (1600,900),
dpi = 150,
ylims = (0.,2.),
# ylims = (0.75,1.25),
legend=false)
savefig("plots/DVMXC_DIS_rho.png")
plot(xexp,DIS_p,layout=(2,3),seriestype=:scatter,
xscale = :log10,
xlims = (0.01,1),
size = (1600,900),
dpi = 150,
ylims = (0.,2.),
# ylims = (0.75,1.25),
legend=false)
savefig("plots/DVMXC_DIS_phi.png")
plot(zmid,Had_s,layout=(5,9),seriestype=:scatter,
size = (3200,1800),
dpi = 300,
ylims = (0.,2.),
# ylims = (0.75,1.25),
legend=false)
savefig("plots/DVMXC_Had_SIDIS.png")
plot(zmid,Had_r,layout=(5,9),seriestype=:scatter,
size = (3200,1800),
dpi = 300,
ylims = (0.,2.),
# ylims = (0.75,1.25),
legend=false)
savefig("plots/DVMXC_Had_rho.png")
plot(zmid,Had_p,layout=(5,9),seriestype=:scatter,
size = (3200,1800),
dpi = 300,
ylims = (0.,2.),
# ylims = (0.75,1.25),
legend=false)
savefig("plots/DVMXC_Had_phi.png")