Commit e5092c5c authored by ulrich_y's avatar ulrich_y
Browse files

Added babar analysis

parent 7fbbb49e
Pipeline #1200 failed with stage
in 10 seconds
......@@ -4,6 +4,8 @@ variables:
before_script:
- date
- export PATH=/afs/psi.ch/sys/psi.x86_64_slp6/Tools/texlive/2017/bin/x86_64-linux/:$PATH
- export PATH=/afs/psi.ch/sys/psi.x86_64_slp6/Programming/psi-python27/4.4.0/bin/:$PATH
- python --version
stages:
- build
......@@ -11,6 +13,8 @@ stages:
compile:
stage: build
script:
- git clone https://"gitlab+deploy-token-2":TToX8HiCR5mV6hLH_UXD@gitlab.psi.ch/mcmule/pymule.git figures/mule/pymule
- export PYTHONPATH=`pwd`/figures/mule/pymule && cd figures/mule && python babar.py
- latexmk -r .latexmkrc manual.tex
- gzip manual.ps
artifacts:
......
## Generated at 11:31 on February 25 2020 by ulrich_y
# git version: redesign (4138c5e)
# To re-generate, run python tools/pymule \
# --seeds 19397 52088 83215 93857 86361 \
# -xi 0.5 0.25 0.125 \
# --flavour tau-e \
# --piece m2enng \
# --output-dir babar-tau-e \
# --prog mcmule \
# --stat 0,1000,10,1000,50 \
# --stat CT,1000,10,1000,50 \
# --stat R,5000,50,10000,100 \
# --stat V,5000,50,10000,100
# specify the program to run relative to `pwd`
binary=mcmule
# specify the output folder
folder=babar-tau-e/
# Specify the variables nenter_ad, itmx_ad, nenter and itmx
# for each piece you want to run.
declare -A STAT=(
["m2enng0"]="1000\n10\n1000\n50"
["m2enngCT"]="1000\n10\n1000\n50"
["m2enngR"]="5000\n50\n10000\n100"
["m2enngV"]="1000\n10\n1000\n50"
)
## Generated at 11:31 on February 25 2020 by ulrich_y
# git version: redesign (4138c5e)
# To re-generate, run python tools/pymule \
# --seeds 19397 52088 83215 93857 86361 \
# -xi 0.5 0.25 0.125 \
# --flavour tau-e \
# --piece m2enng \
# --output-dir babar-tau-e \
# --prog mcmule \
# --stat 0,1000,10,1000,50 \
# --stat CT,1000,10,1000,50 \
# --stat R,5000,50,10000,100 \
# --stat V,5000,50,10000,100
conf babar-tau-e/m2enng-tau-e.conf
#run 19397 1.000000 m2enng0 tau-e 0
#run 52088 1.000000 m2enng0 tau-e 0
#run 83215 1.000000 m2enng0 tau-e 0
#run 93857 1.000000 m2enng0 tau-e 0
#run 86361 1.000000 m2enng0 tau-e 0
#
#
#run 19397 1.000000 m2enngV tau-e 0
#run 52088 1.000000 m2enngV tau-e 0
#run 83215 1.000000 m2enngV tau-e 0
#run 93857 1.000000 m2enngV tau-e 0
#run 86361 1.000000 m2enngV tau-e 0
#
#
#run 19397 0.500000 m2enngCT tau-e 0
#run 52088 0.500000 m2enngCT tau-e 0
#run 83215 0.500000 m2enngCT tau-e 0
#run 93857 0.500000 m2enngCT tau-e 0
#run 86361 0.500000 m2enngCT tau-e 0
run 19397 0.500000 m2enngR tau-e 0
run 52088 0.500000 m2enngR tau-e 0
run 83215 0.500000 m2enngR tau-e 0
run 93857 0.500000 m2enngR tau-e 0
run 86361 0.500000 m2enngR tau-e 0
#run 19397 0.250000 m2enngCT tau-e 0
#run 52088 0.250000 m2enngCT tau-e 0
#run 83215 0.250000 m2enngCT tau-e 0
#run 93857 0.250000 m2enngCT tau-e 0
#run 86361 0.250000 m2enngCT tau-e 0
run 19397 0.250000 m2enngR tau-e 0
run 52088 0.250000 m2enngR tau-e 0
run 83215 0.250000 m2enngR tau-e 0
run 93857 0.250000 m2enngR tau-e 0
run 86361 0.250000 m2enngR tau-e 0
#run 19397 0.125000 m2enngCT tau-e 0
#run 52088 0.125000 m2enngCT tau-e 0
#run 83215 0.125000 m2enngCT tau-e 0
#run 93857 0.125000 m2enngCT tau-e 0
#run 86361 0.125000 m2enngCT tau-e 0
run 19397 0.125000 m2enngR tau-e 0
run 52088 0.125000 m2enngR tau-e 0
run 83215 0.125000 m2enngR tau-e 0
run 93857 0.125000 m2enngR tau-e 0
run 86361 0.125000 m2enngR tau-e 0
import sys
sys.path.append("/home/yannickulrich/Documents/code/mcmule-redesign/tools")
import os
## LaTeX begins in the next line
from pymule import *
# To normalise branching ratios, we need the tau lifetime
lifetime = 1/(1000*(6.582119e-25)/(2.903e-13))
# The folder where McMule has stored the statefiles
setup(folder='babar-tau-e/out.tar.bz2')
# Import LO data and re-scale to branching ratio
LO = scaleset(mergefks(sigma('m2enng0')), GF**2*lifetime*alpha)
# Import NLO corrections from the three pieces
NLO = scaleset(mergefks(
sigma('m2enngR'), # real corrections
sigma('m2enngCT'), # counter term
anyxi=sigma('m2enngV') # virtual corrections
), GF**2*lifetime*alpha**2)
# The branching ratio at NLO = LO + correction
fullNLO = plusnumbers(LO['value'], NLO['value'])
# Print results
print "BR_0 = ", printnumber(LO['value'])
print "dBR = ", printnumber(NLO['value'])
# Produce energy plot
fig1, (ax1, ax2) = kplot(
{'lo': LO['Ee'], 'nlo': NLO['Ee']},
labelx=r"$E_e\,/\,{\rm MeV}$",
labelsigma=r"$\D\mathcal{B}/\D E_e$"
)
ax2.set_ylim(0.8,1.01)
# Produce visible mass plot
fig2, (ax1, ax2) = kplot(
{'lo': LO['minv'], 'nlo': NLO['minv']},
labelx=r"$m_{e\gamma}\,/\,{\rm MeV}$",
labelsigma=r"$\D\mathcal{B}/\D m_{e\gamma}$"
)
ax1.set_yscale('log')
ax1.set_xlim(1000,0) ; ax1.set_ylim(5e-9,1e-3)
ax2.set_ylim(0.8,1.)
# LaTeX ends here
bexp = np.array([1.85e-2, 0.05e-2])
print "Tension to the experimental value: ",np.sqrt(chisq([bexp, fullNLO]))
import pymule.plot,re
pymule.plot.setup_pgf()
fig1.savefig('tau:energy.pgf')
fig2.savefig('tau:minv.pgf')
with open('tau:minv.pgf') as fp:
txt = fp.read()
txt = re.sub(".*\d\d\d\..*\n","", txt)
txt = txt.replace(r"\sffamily\fontsize{10.000000}{12.000000}\selectfont "
+r"\(\displaystyle \textsc{McMule}\)",
r"\(\displaystyle \textsc{McMule}\)")
with open('tau:minv.pgf','w') as fp:
fp.write(txt)
with open('tau:energy.pgf') as fp:
txt = fp.read()
txt = re.sub(".*\d\d\d\..*\n","", txt)
txt = txt.replace(r"\sffamily\fontsize{10.000000}{12.000000}\selectfont "
+r"\(\displaystyle \textsc{McMule}\)",
r"\(\displaystyle \textsc{McMule}\)")
with open('tau:energy.pgf','w') as fp:
fp.write(txt)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment