Code indexing in gitaly is broken and leads to code not being visible to the user. We work on the issue with highest priority.

Skip to content
Snippets Groups Projects
Commit 7ef7a94e authored by ulrich_y's avatar ulrich_y
Browse files

Added demo calculation for the Michel decay

parent dce20085
No related branches found
No related tags found
No related merge requests found
# vim: foldmethod=marker
from pymule import *
gamma0 = Mmu**5/192/pi**3
## Crazy observable{{{
## Load data{{{
setup(folder='with-cut/out.tar.bz2', flavour='mu-e')
loC = scaleset(mergefks(sigma('m2enn0')), 1/gamma0)
fig, ans = mergefkswithplot([[sigma('m2ennR')],[sigma('m2ennV')]])
nloC = scaleset(ans, alpha/gamma0)
nnloC = scaleset(mergefks(
sigma('m2ennRF'),
sigma('m2ennRR'),
sigma('m2ennFF', folder='../validation/ff.tar.bz2')
), alpha**2/gamma0)
##########################################}}}
## Plot data{{{
fig, _ = kplot(
{'lo': loC['xe'], 'nlo': nloC['xe'], 'nnlo': nnloC['xe']},
labelx=r'$x_e$',
labelsigma=r'$\D\sigma / \D x_e / \sigma_0$',
legendopts={'what':'u'},
legend={
'lo': r'$\sigma_0$',
'nlo': r'$\sigma_1$',
'nnlo': r'$\sigma_2$'
}
)
##########################################}}}
##########################################################################}}}
## Generated at 17:23 on August 09 2019 by banerjee_p
# git version: master (4577359)
# To re-generate, run python tools/create-menu.py \
# --seeds 22968 76093 73757 53126 26356 \
# -xi 0.01 0.03 0.1 0.3 0.0 5.0 \
# --flavour mu-e \
# --piece m2enn \
# --output-dir m2ennRR_cut_nohardcut \
# --prog xs_main \
# --stat RR,20000,50,150000,100
## Generated at 08:44 on August 14 2019 by engel_t
# git version: master (4577359)
# To re-generate, run python tools/create-menu.py \
# --seeds 20449 40056 46724 46674 10257 \
# -xi 0.01 0.03 0.1 0.3 0.5 \
# --flavour mu-e \
# --piece m2enn \
# --output-dir m2ennRV_cut_xicut7 \
# --prog xs_main \
# --stat RV,10000,20,100000,100
## Generated at 17:33 on August 07 2019 by engel_t
# git version: master (4577359)
# To re-generate, run python tools/create-menu.py \
# --seeds 72051 42886 38288 10451 68200 \
# -xi 0.01 0.03 0.1 0.3 0.5 \
# --flavour mu-e \
# --piece m2enn \
# --output-dir mue_cut_nlo \
# --prog xs_main \
# --stat 0,10000,10,20000,20 \
# --stat R,40000,30,80000,100 \
# --stat V,20000,20,35000,70
# specify the program to run relative to `pwd`
binary=xs_main
# specify the output folder
folder=with-cut/
# Specify the variables nenter_ad, itmx_ad, nenter and itmx
# for each piece you want to run.
declare -A STAT=(
["m2ennRR"]="20000\n50\n150000\n100"
["m2ennRV"]="10000\n20\n100000\n100"
["m2enn0"]="10000\n10\n20000\n20"
["m2ennR"]="40000\n30\n80000\n100"
["m2ennV"]="20000\n20\n35000\n70"
)
## Generated at 17:33 on August 07 2019 by engel_t
conf with-cut/m2enn-mu-e.conf
# Do LO calculation
run 72051 1.000000 m2enn0 mu-e 0
run 42886 1.000000 m2enn0 mu-e 0
run 38288 1.000000 m2enn0 mu-e 0
run 10451 1.000000 m2enn0 mu-e 0
run 68200 1.000000 m2enn0 mu-e 0
# Do NLO calculation
run 72051 0.010000 m2ennR mu-e 0
run 42886 0.010000 m2ennR mu-e 0
run 38288 0.010000 m2ennR mu-e 0
run 10451 0.010000 m2ennR mu-e 0
run 68200 0.010000 m2ennR mu-e 0
run 72051 0.030000 m2ennR mu-e 0
run 42886 0.030000 m2ennR mu-e 0
run 38288 0.030000 m2ennR mu-e 0
run 10451 0.030000 m2ennR mu-e 0
run 68200 0.030000 m2ennR mu-e 0
run 72051 0.100000 m2ennR mu-e 0
run 42886 0.100000 m2ennR mu-e 0
run 38288 0.100000 m2ennR mu-e 0
run 10451 0.100000 m2ennR mu-e 0
run 68200 0.100000 m2ennR mu-e 0
run 72051 0.300000 m2ennR mu-e 0
run 42886 0.300000 m2ennR mu-e 0
run 38288 0.300000 m2ennR mu-e 0
run 10451 0.300000 m2ennR mu-e 0
run 68200 0.300000 m2ennR mu-e 0
run 72051 0.500000 m2ennR mu-e 0
run 42886 0.500000 m2ennR mu-e 0
run 38288 0.500000 m2ennR mu-e 0
run 10451 0.500000 m2ennR mu-e 0
run 68200 0.500000 m2ennR mu-e 0
run 72051 0.010000 m2ennV mu-e 0
run 42886 0.010000 m2ennV mu-e 0
run 38288 0.010000 m2ennV mu-e 0
run 10451 0.010000 m2ennV mu-e 0
run 68200 0.010000 m2ennV mu-e 0
run 72051 0.030000 m2ennV mu-e 0
run 42886 0.030000 m2ennV mu-e 0
run 38288 0.030000 m2ennV mu-e 0
run 10451 0.030000 m2ennV mu-e 0
run 68200 0.030000 m2ennV mu-e 0
run 72051 0.100000 m2ennV mu-e 0
run 42886 0.100000 m2ennV mu-e 0
run 38288 0.100000 m2ennV mu-e 0
run 10451 0.100000 m2ennV mu-e 0
run 68200 0.100000 m2ennV mu-e 0
run 72051 0.300000 m2ennV mu-e 0
run 42886 0.300000 m2ennV mu-e 0
run 38288 0.300000 m2ennV mu-e 0
run 10451 0.300000 m2ennV mu-e 0
run 68200 0.300000 m2ennV mu-e 0
run 72051 0.500000 m2ennV mu-e 0
run 42886 0.500000 m2ennV mu-e 0
run 38288 0.500000 m2ennV mu-e 0
run 10451 0.500000 m2ennV mu-e 0
run 68200 0.500000 m2ennV mu-e 0
# Do NNLO
run 22968 0.010000 m2ennRR mu-e 0
run 76093 0.010000 m2ennRR mu-e 0
run 73757 0.010000 m2ennRR mu-e 0
run 53126 0.010000 m2ennRR mu-e 0
run 26356 0.010000 m2ennRR mu-e 0
run 22968 0.030000 m2ennRR mu-e 0
run 76093 0.030000 m2ennRR mu-e 0
run 73757 0.030000 m2ennRR mu-e 0
run 53126 0.030000 m2ennRR mu-e 0
run 26356 0.030000 m2ennRR mu-e 0
run 22968 0.100000 m2ennRR mu-e 0
run 76093 0.100000 m2ennRR mu-e 0
run 73757 0.100000 m2ennRR mu-e 0
run 53126 0.100000 m2ennRR mu-e 0
run 26356 0.100000 m2ennRR mu-e 0
run 22968 0.300000 m2ennRR mu-e 0
run 76093 0.300000 m2ennRR mu-e 0
run 73757 0.300000 m2ennRR mu-e 0
run 53126 0.300000 m2ennRR mu-e 0
run 26356 0.300000 m2ennRR mu-e 0
run 22968 0.500000 m2ennRR mu-e 0
run 76093 0.500000 m2ennRR mu-e 0
run 73757 0.500000 m2ennRR mu-e 0
run 53126 0.500000 m2ennRR mu-e 0
run 26356 0.500000 m2ennRR mu-e 0
run 20449 0.010000 m2ennRV mu-e 0
run 40056 0.010000 m2ennRV mu-e 0
run 46724 0.010000 m2ennRV mu-e 0
run 46674 0.010000 m2ennRV mu-e 0
run 10257 0.010000 m2ennRV mu-e 0
run 20449 0.030000 m2ennRV mu-e 0
run 40056 0.030000 m2ennRV mu-e 0
run 46724 0.030000 m2ennRV mu-e 0
run 46674 0.030000 m2ennRV mu-e 0
run 10257 0.030000 m2ennRV mu-e 0
run 20449 0.100000 m2ennRV mu-e 0
run 40056 0.100000 m2ennRV mu-e 0
run 46724 0.100000 m2ennRV mu-e 0
run 46674 0.100000 m2ennRV mu-e 0
run 10257 0.100000 m2ennRV mu-e 0
run 20449 0.300000 m2ennRV mu-e 0
run 40056 0.300000 m2ennRV mu-e 0
run 46724 0.300000 m2ennRV mu-e 0
run 46674 0.300000 m2ennRV mu-e 0
run 10257 0.300000 m2ennRV mu-e 0
run 20449 0.500000 m2ennRV mu-e 0
run 40056 0.500000 m2ennRV mu-e 0
run 46724 0.500000 m2ennRV mu-e 0
run 46674 0.500000 m2ennRV mu-e 0
run 10257 0.500000 m2ennRV mu-e 0
File added
File added
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment