Commit c125fa99 authored by ulrich_y's avatar ulrich_y
Browse files

Added legacy results for muse

parent f2d91ab7
# vim: foldmethod=marker
from pymule import *
from pymule.plot import twopanel, threepanel
from pymule.colours import orderscheme
import pymule.plot
setup(folder='mp2mp_muse_paper/out.tar.bz2')
## Load{{{
lo = scaleset(mergefks(sigma('mp2mp0')), alpha**2*conv)
nloNoVP = scaleset(mergefks(
sigma('mp2mpF'),
sigma('mp2mpR15'), sigma('mp2mpR35')
), alpha**3*conv)
nnloNoVP = scaleset(mergefks(
sigma('mp2mpFF'),
sigma('mp2mpRF15'), sigma('mp2mpRF35'),
sigma('mp2mpRR1516'), sigma('mp2mpRR3536'),
), alpha**4*conv)
nlo = scaleset(mergefks(
sigma('mp2mpF'),
sigma('mp2mpR15'), sigma('mp2mpR35'),
anyxiVP=sigma('mp2mpA')
), alpha**3*conv)
nnlo = scaleset(mergefks(
sigma('mp2mpFF'), sigma('mp2mpAF'),
sigma('mp2mpRF15'), sigma('mp2mpRF35'),
sigma('mp2mpAR15'),sigma('mp2mpAR35'),
sigma('mp2mpRR1516'), sigma('mp2mpRR3536'),
anyxiF=sigma('mp2mpAA'),
anyxiNF=sigma('mp2mpNF')
), alpha**4*conv)
#####################################################################}}}
## Make pictures{{{
### $\theta_e$ w/o & w/ VP{{{
fig, (ax1, ax2, ax3) = kplot(
{
'lo': lo['thetal'],
'nlo': nlo['thetal'],
'nnlo': nnlo['thetal'],
'nlo2': nloNoVP['thetal'],
'nnlo2': nnloNoVP['thetal']
},
labelx="$\\theta_e\,/\,{\\rm deg}$",
labelsigma="$\\D\\sigma/\\D\\theta_\\mu\ /\ {\\rm\\upmu b}$",
legend={
'lo': '$\\sigma^{(0)}$',
'nlo': '$\\sigma^{(1)}$',
'nnlo': '$\\sigma^{(2)}$'
},
legendopts={'what': 'u', 'loc': 'upper right'}
)
ylim(0.99995-1, 1.0002-1)
pymule.plot.mpl_axes_aligner.yaxes(ax2, ax3,y1=0.)
fig.savefig('plots/thetamu.pdf')
###########################################################}}}
### bands{{{
def sanitise(x):
x[x[:,1]==0,1] = -10
return x
bands = ['EmuB1', 'EmuB2', 'EmuB3', 'EmuB4']
bandxs = [addplots(lo[k], addplots(nlo[k], nnlo[k])) for k in bands]
bandk1 = [sanitise(divideplots(nlo[k], lo[k])) for k in bands]
bandk2 = [sanitise(divideplots(nnlo[k], addplots(lo[k], nlo[k]))) for k in bands]
fig, (ax1, ax2, ax3) = threepanel(
labelx="$E_\\mu^{\\rm kin}\,/\,{\\rm MeV}$",
upleft=bandxs, labupleft="$\\D\\sigma/\\D E_\\mu^{\\rm kin}\,/\,{\\rm \\upmu b}$",
middleleft=bandk1, labmiddleleft="$\\delta K^{(1)}$",
downleft=bandk2, labdownleft="$\\delta K^{(2)}$",
colupleft =[orderscheme['nnlo'], 'paleturquoise', 'olivedrab', 'gold'],
colmiddleleft=[orderscheme['nnlo'], 'paleturquoise', 'olivedrab', 'gold'],
coldownleft =[orderscheme['nnlo'], 'paleturquoise', 'olivedrab', 'gold']
)
ax1.set_yscale('log')
ax2.set_ylim(0.982-1,1.016-1)
ax3.set_ylim(0.999-1,1.018-1)
thbounds = np.array([
[
0.3875697, # about 22 degrees -> 126 MeV
0.7708983 # about 44 degrees -> 117 MeV
],[
0.8054385, # about 46 degrees -> 116 MeV
1.0939359 # about 63 degrees -> 107 MeV
],[
1.124739 , # about 64 degrees -> 106 MeV
1.403278 # about 80 degrees -> 97 MeV
],[
1.435036 , # about 82 degrees -> 96 MeV
1.739436 # about 99 degrees -> 87 MeV
]
])
ax2.legend([
"$%4.1f^\\circ \\le \\theta_\\mu \\le %4.1f^\\circ$" % tuple(i)
for i in thbounds/pi * 180
])
mulify(fig)
fig.savefig('plots/bands.pdf')
###########################################################}}}
#####################################################################}}}
group=l-p-scattering
name=legacy
plots=thetamu.pdf bands.pdf
figures=$(addprefix plots/,$(plots))
include ../../tools/makefile.conf
## Generated at 14:19 on June 20 2020 by signer_a
# git version: not on git (not on git)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 39021 33569 12524 12257 19633 88944 73213 13385 44677 10486 65428 67096 92983 54709 43317 74995 16302 57157 52644 84477 94053 20924 21591 49750 43465 86307 22057 62008 84698 64901 44884 98484 62615 32553 38429 34230 68618 46984 27074 80887 48712 96101 74711 68670 45961 45897 37275 65828 96181 45509 89895 20190 22516 92687 78920 26602 86643 72981 47361 26589 38999 20912 13953 88274 18146 51598 79255 55207 80447 80776 42039 93027 40161 56723 69626 \
# -xi 0.5 0.25 0.125 \
# --flavour muse210 \
# --genprocess mp2mp \
# --output-dir mp2mp_muse_paper \
# --prog mcmule \
# --stat R15,20000,20,250000,100,3 \
# --stat F,2000,20,20000,100,3 \
# --stat RR3536,10000,20,100000,100,3 \
# --stat 0,10000,20,500000,100,3 \
# --stat RF15,10000,20,50000,50,3 \
# --stat RR1516,10000,20,100000,100,3 \
# --stat FF,2000,20,20000,100,3 \
# --stat RF35,10000,20,50000,50,3 \
# --stat R35,20000,20,250000,100,3
image mcmule:muonelegacy mp2mp_muse_paper/user.f95
conf mp2mp_muse_paper/mp2mp-muse210.conf
run 69626 0.500000 mp2mpR15 muse210 0
run 56723 0.500000 mp2mpR15 muse210 0
run 40161 0.500000 mp2mpR15 muse210 0
run 93027 0.250000 mp2mpR15 muse210 0
run 42039 0.250000 mp2mpR15 muse210 0
run 80776 0.250000 mp2mpR15 muse210 0
run 80447 0.125000 mp2mpR15 muse210 0
run 55207 0.125000 mp2mpR15 muse210 0
run 79255 0.125000 mp2mpR15 muse210 0
run 51598 0.500000 mp2mpF muse210 0
run 18146 0.500000 mp2mpF muse210 0
run 88274 0.500000 mp2mpF muse210 0
run 13953 0.250000 mp2mpF muse210 0
run 20912 0.250000 mp2mpF muse210 0
run 38999 0.250000 mp2mpF muse210 0
run 26589 0.125000 mp2mpF muse210 0
run 47361 0.125000 mp2mpF muse210 0
run 72981 0.125000 mp2mpF muse210 0
run 65828 1.000000 mp2mp0 muse210 0
run 37275 1.000000 mp2mp0 muse210 0
run 45897 1.000000 mp2mp0 muse210 0
run 22057 0.500000 mp2mpFF muse210 0
run 86307 0.500000 mp2mpFF muse210 0
run 43465 0.500000 mp2mpFF muse210 0
run 49750 0.250000 mp2mpFF muse210 0
run 21591 0.250000 mp2mpFF muse210 0
run 20924 0.250000 mp2mpFF muse210 0
run 94053 0.125000 mp2mpFF muse210 0
run 84477 0.125000 mp2mpFF muse210 0
run 52644 0.125000 mp2mpFF muse210 0
run 44677 0.500000 mp2mpR35 muse210 0
run 13385 0.500000 mp2mpR35 muse210 0
run 73213 0.500000 mp2mpR35 muse210 0
run 88944 0.250000 mp2mpR35 muse210 0
run 19633 0.250000 mp2mpR35 muse210 0
run 12257 0.250000 mp2mpR35 muse210 0
run 12524 0.125000 mp2mpR35 muse210 0
run 33569 0.125000 mp2mpR35 muse210 0
run 39021 0.125000 mp2mpR35 muse210 0
## Generated at 14:19 on June 20 2020 by signer_a
# git version: not on git (not on git)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 39021 33569 12524 12257 19633 88944 73213 13385 44677 10486 65428 67096 92983 54709 43317 74995 16302 57157 52644 84477 94053 20924 21591 49750 43465 86307 22057 62008 84698 64901 44884 98484 62615 32553 38429 34230 68618 46984 27074 80887 48712 96101 74711 68670 45961 45897 37275 65828 96181 45509 89895 20190 22516 92687 78920 26602 86643 72981 47361 26589 38999 20912 13953 88274 18146 51598 79255 55207 80447 80776 42039 93027 40161 56723 69626 \
# -xi 0.5 0.25 0.125 \
# --flavour muse210 \
# --genprocess mp2mp \
# --output-dir mp2mp_muse_paper \
# --prog mcmule \
# --stat R15,20000,20,250000,100,3 \
# --stat F,2000,20,20000,100,3 \
# --stat RR3536,10000,20,100000,100,3 \
# --stat 0,10000,20,500000,100,3 \
# --stat RF15,10000,20,50000,50,3 \
# --stat RR1516,10000,20,100000,100,3 \
# --stat FF,2000,20,20000,100,3 \
# --stat RF35,10000,20,50000,50,3 \
# --stat R35,20000,20,250000,100,3
image mcmule:muonelegacy mp2mp_muse_paper/user.f95
conf mp2mp_muse_paper/mp2mp-muse210.conf
run 86643 0.500000 mp2mpRR3536 muse210 0
run 26602 0.500000 mp2mpRR3536 muse210 0
run 78920 0.500000 mp2mpRR3536 muse210 0
run 92687 0.250000 mp2mpRR3536 muse210 0
run 22516 0.250000 mp2mpRR3536 muse210 0
run 20190 0.250000 mp2mpRR3536 muse210 0
run 89895 0.125000 mp2mpRR3536 muse210 0
run 45509 0.125000 mp2mpRR3536 muse210 0
run 96181 0.125000 mp2mpRR3536 muse210 0
run 45961 0.500000 mp2mpRF15 muse210 0
run 68670 0.500000 mp2mpRF15 muse210 0
run 74711 0.500000 mp2mpRF15 muse210 0
run 96101 0.250000 mp2mpRF15 muse210 0
run 48712 0.250000 mp2mpRF15 muse210 0
run 80887 0.250000 mp2mpRF15 muse210 0
run 27074 0.125000 mp2mpRF15 muse210 0
run 46984 0.125000 mp2mpRF15 muse210 0
run 68618 0.125000 mp2mpRF15 muse210 0
run 34230 0.500000 mp2mpRR1516 muse210 0
run 38429 0.500000 mp2mpRR1516 muse210 0
run 32553 0.500000 mp2mpRR1516 muse210 0
run 62615 0.250000 mp2mpRR1516 muse210 0
run 98484 0.250000 mp2mpRR1516 muse210 0
run 44884 0.250000 mp2mpRR1516 muse210 0
run 64901 0.125000 mp2mpRR1516 muse210 0
run 84698 0.125000 mp2mpRR1516 muse210 0
run 62008 0.125000 mp2mpRR1516 muse210 0
run 57157 0.500000 mp2mpRF35 muse210 0
run 16302 0.500000 mp2mpRF35 muse210 0
run 74995 0.500000 mp2mpRF35 muse210 0
run 43317 0.250000 mp2mpRF35 muse210 0
run 54709 0.250000 mp2mpRF35 muse210 0
run 92983 0.250000 mp2mpRF35 muse210 0
run 67096 0.125000 mp2mpRF35 muse210 0
run 65428 0.125000 mp2mpRF35 muse210 0
run 10486 0.125000 mp2mpRF35 muse210 0
## Generated at 14:19 on June 20 2020 by signer_a
# git version: not on git (not on git)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 39021 33569 12524 12257 19633 88944 73213 13385 44677 10486 65428 67096 92983 54709 43317 74995 16302 57157 52644 84477 94053 20924 21591 49750 43465 86307 22057 62008 84698 64901 44884 98484 62615 32553 38429 34230 68618 46984 27074 80887 48712 96101 74711 68670 45961 45897 37275 65828 96181 45509 89895 20190 22516 92687 78920 26602 86643 72981 47361 26589 38999 20912 13953 88274 18146 51598 79255 55207 80447 80776 42039 93027 40161 56723 69626 \
# -xi 0.5 0.25 0.125 \
# --flavour muse210 \
# --genprocess mp2mp \
# --output-dir mp2mp_muse_paper \
# --prog mcmule \
# --stat R15,20000,20,250000,100,3 \
# --stat F,2000,20,20000,100,3 \
# --stat RR3536,10000,20,100000,100,3 \
# --stat 0,10000,20,500000,100,3 \
# --stat RF15,10000,20,50000,50,3 \
# --stat RR1516,10000,20,100000,100,3 \
# --stat FF,2000,20,20000,100,3 \
# --stat RF35,10000,20,50000,50,3 \
# --stat R35,20000,20,250000,100,3
conf mp2mp_muse_paper/mp2mp-muse210.conf
run 69626 0.500000 mp2mpR15 muse210 0
run 56723 0.500000 mp2mpR15 muse210 0
run 40161 0.500000 mp2mpR15 muse210 0
run 93027 0.250000 mp2mpR15 muse210 0
run 42039 0.250000 mp2mpR15 muse210 0
run 80776 0.250000 mp2mpR15 muse210 0
run 80447 0.125000 mp2mpR15 muse210 0
run 55207 0.125000 mp2mpR15 muse210 0
run 79255 0.125000 mp2mpR15 muse210 0
run 51598 0.500000 mp2mpF muse210 0
run 18146 0.500000 mp2mpF muse210 0
run 88274 0.500000 mp2mpF muse210 0
run 13953 0.250000 mp2mpF muse210 0
run 20912 0.250000 mp2mpF muse210 0
run 38999 0.250000 mp2mpF muse210 0
run 26589 0.125000 mp2mpF muse210 0
run 47361 0.125000 mp2mpF muse210 0
run 72981 0.125000 mp2mpF muse210 0
run 86643 0.500000 mp2mpRR3536 muse210 0
run 26602 0.500000 mp2mpRR3536 muse210 0
run 78920 0.500000 mp2mpRR3536 muse210 0
run 92687 0.250000 mp2mpRR3536 muse210 0
run 22516 0.250000 mp2mpRR3536 muse210 0
run 20190 0.250000 mp2mpRR3536 muse210 0
run 89895 0.125000 mp2mpRR3536 muse210 0
run 45509 0.125000 mp2mpRR3536 muse210 0
run 96181 0.125000 mp2mpRR3536 muse210 0
run 65828 1.000000 mp2mp0 muse210 0
run 37275 1.000000 mp2mp0 muse210 0
run 45897 1.000000 mp2mp0 muse210 0
run 45961 0.500000 mp2mpRF15 muse210 0
run 68670 0.500000 mp2mpRF15 muse210 0
run 74711 0.500000 mp2mpRF15 muse210 0
run 96101 0.250000 mp2mpRF15 muse210 0
run 48712 0.250000 mp2mpRF15 muse210 0
run 80887 0.250000 mp2mpRF15 muse210 0
run 27074 0.125000 mp2mpRF15 muse210 0
run 46984 0.125000 mp2mpRF15 muse210 0
run 68618 0.125000 mp2mpRF15 muse210 0
run 34230 0.500000 mp2mpRR1516 muse210 0
run 38429 0.500000 mp2mpRR1516 muse210 0
run 32553 0.500000 mp2mpRR1516 muse210 0
run 62615 0.250000 mp2mpRR1516 muse210 0
run 98484 0.250000 mp2mpRR1516 muse210 0
run 44884 0.250000 mp2mpRR1516 muse210 0
run 64901 0.125000 mp2mpRR1516 muse210 0
run 84698 0.125000 mp2mpRR1516 muse210 0
run 62008 0.125000 mp2mpRR1516 muse210 0
run 22057 0.500000 mp2mpFF muse210 0
run 86307 0.500000 mp2mpFF muse210 0
run 43465 0.500000 mp2mpFF muse210 0
run 49750 0.250000 mp2mpFF muse210 0
run 21591 0.250000 mp2mpFF muse210 0
run 20924 0.250000 mp2mpFF muse210 0
run 94053 0.125000 mp2mpFF muse210 0
run 84477 0.125000 mp2mpFF muse210 0
run 52644 0.125000 mp2mpFF muse210 0
run 57157 0.500000 mp2mpRF35 muse210 0
run 16302 0.500000 mp2mpRF35 muse210 0
run 74995 0.500000 mp2mpRF35 muse210 0
run 43317 0.250000 mp2mpRF35 muse210 0
run 54709 0.250000 mp2mpRF35 muse210 0
run 92983 0.250000 mp2mpRF35 muse210 0
run 67096 0.125000 mp2mpRF35 muse210 0
run 65428 0.125000 mp2mpRF35 muse210 0
run 10486 0.125000 mp2mpRF35 muse210 0
run 44677 0.500000 mp2mpR35 muse210 0
run 13385 0.500000 mp2mpR35 muse210 0
run 73213 0.500000 mp2mpR35 muse210 0
run 88944 0.250000 mp2mpR35 muse210 0
run 19633 0.250000 mp2mpR35 muse210 0
run 12257 0.250000 mp2mpR35 muse210 0
run 12524 0.125000 mp2mpR35 muse210 0
run 33569 0.125000 mp2mpR35 muse210 0
run 39021 0.125000 mp2mpR35 muse210 0
## Generated at 14:19 on June 20 2020 by signer_a
# git version: not on git (not on git)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 39021 33569 12524 12257 19633 88944 73213 13385 44677 10486 65428 67096 92983 54709 43317 74995 16302 57157 52644 84477 94053 20924 21591 49750 43465 86307 22057 62008 84698 64901 44884 98484 62615 32553 38429 34230 68618 46984 27074 80887 48712 96101 74711 68670 45961 45897 37275 65828 96181 45509 89895 20190 22516 92687 78920 26602 86643 72981 47361 26589 38999 20912 13953 88274 18146 51598 79255 55207 80447 80776 42039 93027 40161 56723 69626 \
# -xi 0.5 0.25 0.125 \
# --flavour muse210 \
# --genprocess mp2mp \
# --output-dir mp2mp_muse_paper \
# --prog mcmule \
# --stat R15,20000,20,250000,100,3 \
# --stat F,2000,20,20000,100,3 \
# --stat RR3536,10000,20,100000,100,3 \
# --stat 0,10000,20,500000,100,3 \
# --stat RF15,10000,20,50000,50,3 \
# --stat RR1516,10000,20,100000,100,3 \
# --stat FF,2000,20,20000,100,3 \
# --stat RF35,10000,20,50000,50,3 \
# --stat R35,20000,20,250000,100,3
image mcmule:muonelegacy mp2mp_muse_paper/user.f95
conf mp2mp_muse_paper/mp2mp-muse210.conf
#A
run 88874 1.000000 mp2mpA muse210 0
run 47404 1.000000 mp2mpA muse210 0
run 17404 1.000000 mp2mpA muse210 0
run 27404 1.000000 mp2mpA muse210 0
#AA
run 82874 1.000000 mp2mpAA muse210 0
run 82109 1.000000 mp2mpAA muse210 0
run 12109 1.000000 mp2mpAA muse210 0
run 32109 1.000000 mp2mpAA muse210 0
#AF
run 65011 0.100000 mp2mpAF muse210 0
run 25811 0.100000 mp2mpAF muse210 0
run 15811 0.100000 mp2mpAF muse210 0
run 51011 0.300000 mp2mpAF muse210 0
run 89611 0.300000 mp2mpAF muse210 0
run 19611 0.300000 mp2mpAF muse210 0
run 67511 0.500000 mp2mpAF muse210 0
run 63511 0.500000 mp2mpAF muse210 0
run 13511 0.500000 mp2mpAF muse210 0
run 94011 1.000000 mp2mpAF muse210 0
run 12911 1.000000 mp2mpAF muse210 0
run 32911 1.000000 mp2mpAF muse210 0
#AR15
run 28026 0.100000 mp2mpAR15 muse210 0
run 14014 0.100000 mp2mpAR15 muse210 0
run 34014 0.100000 mp2mpAR15 muse210 0
run 30193 0.300000 mp2mpAR15 muse210 0
run 72123 0.300000 mp2mpAR15 muse210 0
run 12123 0.300000 mp2mpAR15 muse210 0
run 12430 0.500000 mp2mpAR15 muse210 0
run 21526 0.500000 mp2mpAR15 muse210 0
run 31526 0.500000 mp2mpAR15 muse210 0
run 36226 1.000000 mp2mpAR15 muse210 0
run 90257 1.000000 mp2mpAR15 muse210 0
run 10257 1.000000 mp2mpAR15 muse210 0
#AR35
run 63997 0.100000 mp2mpAR35 muse210 0
run 21847 0.100000 mp2mpAR35 muse210 0
run 11847 0.100000 mp2mpAR35 muse210 0
run 71160 0.300000 mp2mpAR35 muse210 0
run 38322 0.300000 mp2mpAR35 muse210 0
run 18322 0.300000 mp2mpAR35 muse210 0
run 78520 0.500000 mp2mpAR35 muse210 0
run 32026 0.500000 mp2mpAR35 muse210 0
run 12026 0.500000 mp2mpAR35 muse210 0
run 15629 1.000000 mp2mpAR35 muse210 0
run 68351 1.000000 mp2mpAR35 muse210 0
run 18351 1.000000 mp2mpAR35 muse210 0
## Generated at 14:19 on June 20 2020 by signer_a
# git version: not on git (not on git)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 39021 33569 12524 12257 19633 88944 73213 13385 44677 10486 65428 67096 92983 54709 43317 74995 16302 57157 52644 84477 94053 20924 21591 49750 43465 86307 22057 62008 84698 64901 44884 98484 62615 32553 38429 34230 68618 46984 27074 80887 48712 96101 74711 68670 45961 45897 37275 65828 96181 45509 89895 20190 22516 92687 78920 26602 86643 72981 47361 26589 38999 20912 13953 88274 18146 51598 79255 55207 80447 80776 42039 93027 40161 56723 69626 \
# -xi 0.5 0.25 0.125 \
# --flavour muse210 \
# --genprocess mp2mp \
# --output-dir mp2mp_muse_paper \
# --prog mcmule \
# --stat R15,20000,20,250000,100,3 \
# --stat F,2000,20,20000,100,3 \
# --stat RR3536,10000,20,100000,100,3 \
# --stat 0,10000,20,500000,100,3 \
# --stat RF15,10000,20,50000,50,3 \
# --stat RR1516,10000,20,100000,100,3 \
# --stat FF,2000,20,20000,100,3 \
# --stat RF35,10000,20,50000,50,3 \
# --stat R35,20000,20,250000,100,3
# specify the program to run relative to `pwd`
binary=mcmule
# specify the output folder
folder=mp2mp_muse_paper/
# Specify the variables nenter_ad, itmx_ad, nenter and itmx
# for each piece you want to run.
declare -A STAT=(
["mp2mpR15"]="20000\n20\n250000\n100"
["mp2mpF"]="2000\n20\n20000\n100"
["mp2mpRR3536"]="10000\n20\n100000\n100"
["mp2mp0"]="10000\n20\n500000\n100"
["mp2mpRF15"]="10000\n20\n50000\n50"
["mp2mpRR1516"]="10000\n20\n100000\n100"
["mp2mpFF"]="2000\n20\n20000\n100"
["mp2mpRF35"]="10000\n20\n50000\n50"
["mp2mpR35"]="20000\n20\n250000\n100"
["mp2mpA"]="1000\n20\n2000\n50"
["mp2mpAA"]="1000\n20\n2000\n50"
["mp2mpAF"]="1000\n20\n2000\n50"
["mp2mpAR15"]="20000\n20\n80000\n70"
["mp2mpAR35"]="20000\n20\n80000\n70"
)
!!!!!!!!!!!!!!!!!!!!!
MODULE USER
!!!!!!!!!!!!!!!!!!!!!
use functions
implicit none
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
integer, parameter :: nr_q = 6
integer, parameter :: nr_bins = 120
real, parameter :: &
min_val(nr_q) = (/ 20., 10., 10., 10., 10., 10. /)
real, parameter :: &
max_val(nr_q) = (/ 100., 130.,130., 130., 130., 130. /)
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!! ============================================== !!
!! DO NOT EVEN THINK ABOUT CHANGING ANYTHING HERE !!
!! ============================================== !!
logical :: pass_cut(nr_q)
character (len = 6), dimension(nr_q) :: names
character(len=10) :: filenamesuffix
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!! ----------------------------------------- !!
!! There are two versions of binning !!
!! One for computing d \sigma/ d Q !!
!! One for computing Q d \sigma/ d Q !!
!! choose by setting the variable bin_kind !!
!! ----------------------------------------- !!
integer :: bin_kind = 0 !! 0 for d \sig/dQ; +1 for Q d \sig/dQ;
contains
SUBROUTINE FIX_MU
!! ==== Specify the scale mu AND musq==mu**2 ==== !!