"git@gitlab.psi.ch:mcmule/user-library.git" did not exist on "87955c9e4bca9c45754f282223f1a86ad52e9ec8"
Newer
Older
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
# vim: foldmethod=marker
from pymule import *
setup(folder='ep2ep_mesa_paper/out.tar.bz2')
## Load{{{
lo = scaleset(mergefks(sigma('mp2mp0')), alpha**2*conv)
fig, nloNoVP = mergefkswithplot([
[sigma('mp2mpF')],
[sigma('mp2mpR15'), sigma('mp2mpR35')]
], scale=alpha**3*conv, xlim=[-3,0])
fig.savefig('plots/xi-nlo.pdf')
fig, nnloNoVP = mergefkswithplot([
[sigma('mp2mpFF')],
[sigma('mp2mpRF15'), sigma('mp2mpRF35')],
[sigma('mp2mpRR1516'), sigma('mp2mpRR3536')],
], scale=alpha**4*conv, xlim=[-3, 0])
fig.savefig('plots/xi-nnlo.pdf')
nlo = scaleset(mergefks(
sigma('mp2mpF'),
sigma('mp2mpR15'), sigma('mp2mpR35'),
anyxiVP=sigma('mp2mpA', obs='1')
), alpha**3*conv)
nnlo = scaleset(mergefks(
sigma('mp2mpFF'), sigma('mp2mpAF', obs='1'),
sigma('mp2mpRF15'), sigma('mp2mpRF35'),
sigma('mp2mpAR15', obs='1'),sigma('mp2mpAR35', obs='1'),
sigma('mp2mpRR1516'), sigma('mp2mpRR3536'),
anyxiF=sigma('mp2mpAA', obs='1'),
anyxiNF=sigma('mp2mpNF', obs='1')
), alpha**4*conv)
nnloNF = scaleset(mergefks(
sigma('mp2mpFF'),
sigma('mp2mpRF15'), sigma('mp2mpRF35'),
sigma('mp2mpRR1516'), sigma('mp2mpRR3536'),
anyxiNF=sigma('mp2mpNF', obs='2')
), alpha**4*conv)
#####################################################################}}}
## Print numbers{{{
print("run time %f days" % ((lo['time'] + nlo['time'] + nnlo['time']) / 3600. / 24.))
print("\\sigma^(0) = " + printnumber(lo['value']))
print("\\sigma^(1) = +" + printnumber(nlo['value']))
print("\\sigma^(2) = " + printnumber(nnlo['value']))
#####################################################################}}}
## Make pictures{{{
### $\theta_e${{{
fig, (ax1, ax2, ax3) = kplot(
{
'lo': mergebins( lo['thetae'],4),
'nlo': mergebins( nlo['thetae'],4),
'nnlo': mergebins(nnlo['thetae'],4),
'nlo2': mergebins(nloNoVP ['thetae'],4),
'nnlo2':mergebins(nnloNoVP['thetae'],4)
},
labelx="$\\theta_e\,/\,{\\rm deg}$",
labelsigma="$\\D\\sigma/\\D\\theta_e\ /\ {\\rm\\upmu b}$",
legend={
'lo': '$\\sigma^{(0)}$',
'nlo': '$\\sigma^{(1)}$',
'nnlo': '$\\sigma^{(2)}$',
'nlo2': '$\\sigma^{(1)}_\\text{no VP}$',
'nnlo2': '$\\sigma^{(2)}_\\text{no VP}$'
},
legendopts={'what': 'u', 'loc': 'upper right'}
)
fig.savefig('plots/thetae.pdf')
###########################################################}}}
### $Q_e^2$ and $Q_p^2${{{
fig, (ax1, ax2, ax3) = kplot(
{
'lo': mergebins( lo['qsqP'][10:189],4),
'nlo': mergebins( nlo['qsqP'][10:189],4),
'nnlo': mergebins(nnlo['qsqP'][10:189],4),
'nlo2': mergebins( nlo['qsqE'][10:189],4),
'nnlo2':mergebins(nnlo['qsqE'][10:189],4),
},
labelx="$|t|\,/\,{\\rm MeV}^2$",
labelsigma="$\\D\\sigma/\\D|t|\ /\ {\\rm\\upmu b}$",
legend={
'lo': '$\\sigma^{(0)}$',
'nlo': '$\\sigma^{(1)}$', 'nlo2': '$\\sigma^{(1)}(|t|)$',
'nnlo': '$\\sigma^{(2)}$', 'nnlo2': '$\\sigma^{(2)}(|t_e|)$'
},
legendopts={'what': 'u', 'loc': 'upper right'}
)
ax2.set_ylim(-0.07,0.07)
ax3.set_ylim(-0.0042,0.0042)
updateaxis(ax3, fig,2)
fig.savefig('plots/qsq.pdf')
###########################################################}}}
#####################################################################}}}