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 f2d91ab7 authored by ulrich_y's avatar ulrich_y
Browse files

Added legacy results for mesa

parent 7984bdd1
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1010 additions and 0 deletions
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/ep2ep-mesa.conf 0 → 100644
+ 47
− 0
View file @ f2d91ab7
## Generated at 14:41 on June 05 2020 by banerjee_p
# git version: muone (633de6a)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 78791 99176 99400 57304 87126 48649 36976 44095 61250 99325 61979 30657 18053 53605 59453 95931 43007 25576 62225 15560 62975 69496 85297 49902 14503 78405 31605 14547 14435 86649 64045 11659 84222 85866 10086 83133 81347 78992 40123 12503 71951 80118 94868 17067 80498 63992 40792 71366 60865 41174 88757 66479 55650 95443 32289 68910 42777 14277 33970 82374 23207 83878 97401 72964 29068 30987 74504 48455 31166 27348 89263 74712 64065 91937 40341 54035 88995 75514 56929 98873 25571 15903 10429 94262 95399 21712 48943 30705 65659 84162 94806 84548 27696 97751 74628 14444 70056 62855 24675 38861 73519 67972 32460 93425 80863 50923 87391 55348 62352 32750 15703 92481 41662 46300 74588 95722 94565 12030 60834 89406 14407 13918 17284 55896 93199 73614 67503 32544 80365 41439 84025 94198 16967 68097 97704 37547 64204 26946 92605 75265 44830 10905 37617 70670 47815 60160 13267 47260 76475 96246 82030 50546 67486 67792 60342 47135 44099 62668 89065 53772 88980 19444 63783 97893 68979 64692 90878 12044 86810 51338 67347 96841 73307 74758 12466 69196 78597 34334 39305 96104 13830 56002 17599 54662 88231 68184 54264 54173 76464 55853 76287 62300 72723 49176 28460 83269 82493 50392 65879 48644 63357 54697 54802 41959 40163 57903 21136 87913 94059 20460 49516 63834 99548 21938 78601 63714 73561 12964 67037 27998 69217 35001 \
# -xi 0.1 0.3 0.5 1.0 \
# --flavour mesa \
# --genprocess mp2mp \
# --output-dir ep2ep_mesa_paper \
# --prog mcmule \
# --stat R15,20000,20,80000,70,6 \
# --stat F,5000,20,10000,70,6 \
# --stat RFco,20000,20,40000,2,6 \
# --stat RR3536536,500000,20,1200000,70,6 \
# --stat 0,5000,20,10000,70,6 \
# --stat RF15,20000,20,80000,70,6 \
# --stat RR1516516,500000,20,1200000,70,6 \
# --stat FF,5000,20,10000,70,6 \
# --stat RF35,20000,20,80000,70,6 \
# --stat R35,20000,20,80000,70,6
# specify the program to run relative to `pwd`
binary=mcmule
# specify the output folder
folder=ep2ep_mesa_paper/
# Specify the variables nenter_ad, itmx_ad, nenter and itmx
# for each piece you want to run.
declare -A STAT=(
["mp2mp0"]="5000\n20\n10000\n70"
["mp2mpF"]="5000\n20\n10000\n70"
["mp2mpR15"]="20000\n20\n80000\n70"
["mp2mpR35"]="20000\n20\n80000\n70"
["mp2mpRco"]="2000000\n1\n1\n1"
["mp2mpFF"]="5000\n20\n10000\n70"
["mp2mpRF35"]="20000\n20\n80000\n70"
["mp2mpRF15"]="20000\n20\n80000\n70"
["mp2mpRR3536"]="500000\n20\n1200000\n70"
["mp2mpRR1516"]="500000\n20\n1200000\n70"
["mp2mpRRc"]="2000000\n1\n1\n1"
["mp2mpA"]="5000\n20\n10000\n70"
["mp2mpAA"]="5000\n20\n10000\n70"
["mp2mpAF"]="5000\n20\n10000\n70"
["mp2mpAR15"]="20000\n20\n80000\n70"
["mp2mpAR35"]="20000\n20\n80000\n70"
["mp2mpNF"]="5000\n20\n10000\n40"
)
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/menu-ep2ep-mesa.menu 0 → 100644
+ 212
− 0
View file @ f2d91ab7
## Generated at 14:41 on June 05 2020 by banerjee_p
# git version: muone (633de6a)
# To re-generate, run python /afs/psi.ch/project/muondecay/.usr/bin/pymule create \
# --seeds 78791 99176 99400 57304 87126 48649 36976 44095 61250 99325 61979 30657 18053 53605 59453 95931 43007 25576 62225 15560 62975 69496 85297 49902 14503 78405 31605 14547 14435 86649 64045 11659 84222 85866 10086 83133 81347 78992 40123 12503 71951 80118 94868 17067 80498 63992 40792 71366 60865 41174 88757 66479 55650 95443 32289 68910 42777 14277 33970 82374 23207 83878 97401 72964 29068 30987 74504 48455 31166 27348 89263 74712 64065 91937 40341 54035 88995 75514 56929 98873 25571 15903 10429 94262 95399 21712 48943 30705 65659 84162 94806 84548 27696 97751 74628 14444 70056 62855 24675 38861 73519 67972 32460 93425 80863 50923 87391 55348 62352 32750 15703 92481 41662 46300 74588 95722 94565 12030 60834 89406 14407 13918 17284 55896 93199 73614 67503 32544 80365 41439 84025 94198 16967 68097 97704 37547 64204 26946 92605 75265 44830 10905 37617 70670 47815 60160 13267 47260 76475 96246 82030 50546 67486 67792 60342 47135 44099 62668 89065 53772 88980 19444 63783 97893 68979 64692 90878 12044 86810 51338 67347 96841 73307 74758 12466 69196 78597 34334 39305 96104 13830 56002 17599 54662 88231 68184 54264 54173 76464 55853 76287 62300 72723 49176 28460 83269 82493 50392 65879 48644 63357 54697 54802 41959 40163 57903 21136 87913 94059 20460 49516 63834 99548 21938 78601 63714 73561 12964 67037 27998 69217 35001 \
# -xi 0.1 0.3 0.5 1.0 \
# --flavour mesa \
# --genprocess mp2mp \
# --output-dir mp2mp_nnlo \
# --prog mcmule \
# --stat R15,20000,20,80000,70,6 \
# --stat F,5000,20,10000,70,6 \
# --stat RFco,20000,20,40000,2,6 \
# --stat RR3536536,500000,20,1200000,70,6 \
# --stat 0,5000,20,10000,70,6 \
# --stat RF15,20000,20,80000,70,6 \
# --stat RR1516516,500000,20,1200000,70,6 \
# --stat FF,5000,20,10000,70,6 \
# --stat RF35,20000,20,80000,70,6 \
# --stat R35,20000,20,80000,70,6
conf ep2ep_mesa_paper/ep2ep-mesa.conf
run 73614 1.000000 mp2mp0 mesa 0
run 93199 1.000000 mp2mp0 mesa 0
run 55896 1.000000 mp2mp0 mesa 0
run 17284 1.000000 mp2mp0 mesa 0
run 50392 0.100000 mp2mpF mesa 0
run 82493 0.100000 mp2mpF mesa 0
run 83269 0.100000 mp2mpF mesa 0
run 28460 0.100000 mp2mpF mesa 0
run 62300 0.300000 mp2mpF mesa 0
run 76287 0.300000 mp2mpF mesa 0
run 55853 0.300000 mp2mpF mesa 0
run 76464 0.300000 mp2mpF mesa 0
run 68184 0.500000 mp2mpF mesa 0
run 88231 0.500000 mp2mpF mesa 0
run 54662 0.500000 mp2mpF mesa 0
run 17599 0.500000 mp2mpF mesa 0
run 96104 1.000000 mp2mpF mesa 0
run 39305 1.000000 mp2mpF mesa 0
run 34334 1.000000 mp2mpF mesa 0
run 78597 1.000000 mp2mpF mesa 0
run 35001 0.100000 mp2mpR15 mesa 0
run 69217 0.100000 mp2mpR15 mesa 0
run 27998 0.100000 mp2mpR15 mesa 0
run 67037 0.100000 mp2mpR15 mesa 0
run 63714 0.300000 mp2mpR15 mesa 0
run 78601 0.300000 mp2mpR15 mesa 0
run 21938 0.300000 mp2mpR15 mesa 0
run 99548 0.300000 mp2mpR15 mesa 0
run 20460 0.500000 mp2mpR15 mesa 0
run 94059 0.500000 mp2mpR15 mesa 0
run 87913 0.500000 mp2mpR15 mesa 0
run 21136 0.500000 mp2mpR15 mesa 0
run 41959 1.000000 mp2mpR15 mesa 0
run 54802 1.000000 mp2mpR15 mesa 0
run 54697 1.000000 mp2mpR15 mesa 0
run 63357 1.000000 mp2mpR15 mesa 0
run 49902 0.100000 mp2mpR35 mesa 0
run 85297 0.100000 mp2mpR35 mesa 0
run 69496 0.100000 mp2mpR35 mesa 0
run 62975 0.100000 mp2mpR35 mesa 0
run 25576 0.300000 mp2mpR35 mesa 0
run 43007 0.300000 mp2mpR35 mesa 0
run 95931 0.300000 mp2mpR35 mesa 0
run 59453 0.300000 mp2mpR35 mesa 0
run 30657 0.500000 mp2mpR35 mesa 0
run 61979 0.500000 mp2mpR35 mesa 0
run 99325 0.500000 mp2mpR35 mesa 0
run 61250 0.500000 mp2mpR35 mesa 0
run 48649 1.000000 mp2mpR35 mesa 0
run 87126 1.000000 mp2mpR35 mesa 0
run 57304 1.000000 mp2mpR35 mesa 0
run 99400 1.000000 mp2mpR35 mesa 0
run 74712 0.100000 mp2mpFF mesa 0
run 89263 0.100000 mp2mpFF mesa 0
run 27348 0.100000 mp2mpFF mesa 0
run 31166 0.100000 mp2mpFF mesa 0
run 30987 0.300000 mp2mpFF mesa 0
run 29068 0.300000 mp2mpFF mesa 0
run 72964 0.300000 mp2mpFF mesa 0
run 97401 0.300000 mp2mpFF mesa 0
run 82374 0.500000 mp2mpFF mesa 0
run 33970 0.500000 mp2mpFF mesa 0
run 14277 0.500000 mp2mpFF mesa 0
run 42777 0.500000 mp2mpFF mesa 0
run 95443 1.000000 mp2mpFF mesa 0
run 55650 1.000000 mp2mpFF mesa 0
run 66479 1.000000 mp2mpFF mesa 0
run 88757 1.000000 mp2mpFF mesa 0
run 14444 0.100000 mp2mpRR1516 mesa 0
run 74628 0.100000 mp2mpRR1516 mesa 0
run 97751 0.100000 mp2mpRR1516 mesa 0
run 27696 0.100000 mp2mpRR1516 mesa 0
run 84162 0.300000 mp2mpRR1516 mesa 0
run 65659 0.300000 mp2mpRR1516 mesa 0
run 30705 0.300000 mp2mpRR1516 mesa 0
run 48943 0.300000 mp2mpRR1516 mesa 0
run 94262 0.500000 mp2mpRR1516 mesa 0
run 10429 0.500000 mp2mpRR1516 mesa 0
run 15903 0.500000 mp2mpRR1516 mesa 0
run 25571 0.500000 mp2mpRR1516 mesa 0
run 75514 1.000000 mp2mpRR1516 mesa 0
run 88995 1.000000 mp2mpRR1516 mesa 0
run 54035 1.000000 mp2mpRR1516 mesa 0
run 40341 1.000000 mp2mpRR1516 mesa 0
run 96246 0.100000 mp2mpRR3536 mesa 0
run 76475 0.100000 mp2mpRR3536 mesa 0
run 47260 0.100000 mp2mpRR3536 mesa 0
run 13267 0.100000 mp2mpRR3536 mesa 0
run 70670 0.300000 mp2mpRR3536 mesa 0
run 37617 0.300000 mp2mpRR3536 mesa 0
run 10905 0.300000 mp2mpRR3536 mesa 0
run 44830 0.300000 mp2mpRR3536 mesa 0
run 26946 0.500000 mp2mpRR3536 mesa 0
run 64204 0.500000 mp2mpRR3536 mesa 0
run 37547 0.500000 mp2mpRR3536 mesa 0
run 97704 0.500000 mp2mpRR3536 mesa 0
run 94198 1.000000 mp2mpRR3536 mesa 0
run 84025 1.000000 mp2mpRR3536 mesa 0
run 41439 1.000000 mp2mpRR3536 mesa 0
run 80365 1.000000 mp2mpRR3536 mesa 0
run 89406 0.100000 mp2mpRF15 mesa 0
run 60834 0.100000 mp2mpRF15 mesa 0
run 12030 0.100000 mp2mpRF15 mesa 0
run 94565 0.100000 mp2mpRF15 mesa 0
run 46300 0.300000 mp2mpRF15 mesa 0
run 41662 0.300000 mp2mpRF15 mesa 0
run 92481 0.300000 mp2mpRF15 mesa 0
run 15703 0.300000 mp2mpRF15 mesa 0
run 55348 0.500000 mp2mpRF15 mesa 0
run 87391 0.500000 mp2mpRF15 mesa 0
run 50923 0.500000 mp2mpRF15 mesa 0
run 80863 0.500000 mp2mpRF15 mesa 0
run 67972 1.000000 mp2mpRF15 mesa 0
run 73519 1.000000 mp2mpRF15 mesa 0
run 38861 1.000000 mp2mpRF15 mesa 0
run 24675 1.000000 mp2mpRF15 mesa 0
run 71366 0.100000 mp2mpRF35 mesa 0
run 40792 0.100000 mp2mpRF35 mesa 0
run 63992 0.100000 mp2mpRF35 mesa 0
run 80498 0.100000 mp2mpRF35 mesa 0
run 80118 0.300000 mp2mpRF35 mesa 0
run 71951 0.300000 mp2mpRF35 mesa 0
run 12503 0.300000 mp2mpRF35 mesa 0
run 40123 0.300000 mp2mpRF35 mesa 0
run 83133 0.500000 mp2mpRF35 mesa 0
run 10086 0.500000 mp2mpRF35 mesa 0
run 85866 0.500000 mp2mpRF35 mesa 0
run 84222 0.500000 mp2mpRF35 mesa 0
run 86649 1.000000 mp2mpRF35 mesa 0
run 14435 1.000000 mp2mpRF35 mesa 0
run 14547 1.000000 mp2mpRF35 mesa 0
run 31605 1.000000 mp2mpRF35 mesa 0
run 47135 1.000000 mp2mpRRc mesa 0
run 47135 1.000000 mp2mpRco mesa 0
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/menu-long.menu 0 → 100644
+ 98
− 0
View file @ f2d91ab7
## Generated at 14:41 on June 05 2020 by banerjee_p
# git version: muone (54e3690d)
image mcmule:muonelegacy ep2ep_mesa_paper/user.f95
conf ep2ep_mesa_paper/ep2ep-mesa.conf
run 14444 0.100000 mp2mpRR1516 mesa 0
run 74628 0.100000 mp2mpRR1516 mesa 0
run 97751 0.100000 mp2mpRR1516 mesa 0
run 27696 0.100000 mp2mpRR1516 mesa 0
run 84162 0.300000 mp2mpRR1516 mesa 0
run 65659 0.300000 mp2mpRR1516 mesa 0
run 30705 0.300000 mp2mpRR1516 mesa 0
run 48943 0.300000 mp2mpRR1516 mesa 0
run 94262 0.500000 mp2mpRR1516 mesa 0
run 10429 0.500000 mp2mpRR1516 mesa 0
run 15903 0.500000 mp2mpRR1516 mesa 0
run 25571 0.500000 mp2mpRR1516 mesa 0
run 75514 1.000000 mp2mpRR1516 mesa 0
run 88995 1.000000 mp2mpRR1516 mesa 0
run 54035 1.000000 mp2mpRR1516 mesa 0
run 40341 1.000000 mp2mpRR1516 mesa 0
run 96246 0.100000 mp2mpRR3536 mesa 0
run 76475 0.100000 mp2mpRR3536 mesa 0
run 47260 0.100000 mp2mpRR3536 mesa 0
run 13267 0.100000 mp2mpRR3536 mesa 0
run 70670 0.300000 mp2mpRR3536 mesa 0
run 37617 0.300000 mp2mpRR3536 mesa 0
run 10905 0.300000 mp2mpRR3536 mesa 0
run 44830 0.300000 mp2mpRR3536 mesa 0
run 26946 0.500000 mp2mpRR3536 mesa 0
run 64204 0.500000 mp2mpRR3536 mesa 0
run 37547 0.500000 mp2mpRR3536 mesa 0
run 97704 0.500000 mp2mpRR3536 mesa 0
run 94198 1.000000 mp2mpRR3536 mesa 0
run 84025 1.000000 mp2mpRR3536 mesa 0
run 41439 1.000000 mp2mpRR3536 mesa 0
run 80365 1.000000 mp2mpRR3536 mesa 0
run 89406 0.100000 mp2mpRF15 mesa 0
run 60834 0.100000 mp2mpRF15 mesa 0
run 12030 0.100000 mp2mpRF15 mesa 0
run 94565 0.100000 mp2mpRF15 mesa 0
run 46300 0.300000 mp2mpRF15 mesa 0
run 41662 0.300000 mp2mpRF15 mesa 0
run 92481 0.300000 mp2mpRF15 mesa 0
run 15703 0.300000 mp2mpRF15 mesa 0
run 55348 0.500000 mp2mpRF15 mesa 0
run 87391 0.500000 mp2mpRF15 mesa 0
run 50923 0.500000 mp2mpRF15 mesa 0
run 80863 0.500000 mp2mpRF15 mesa 0
run 67972 1.000000 mp2mpRF15 mesa 0
run 73519 1.000000 mp2mpRF15 mesa 0
run 38861 1.000000 mp2mpRF15 mesa 0
run 24675 1.000000 mp2mpRF15 mesa 0
run 71366 0.100000 mp2mpRF35 mesa 0
run 40792 0.100000 mp2mpRF35 mesa 0
run 63992 0.100000 mp2mpRF35 mesa 0
run 80498 0.100000 mp2mpRF35 mesa 0
run 80118 0.300000 mp2mpRF35 mesa 0
run 71951 0.300000 mp2mpRF35 mesa 0
run 12503 0.300000 mp2mpRF35 mesa 0
run 40123 0.300000 mp2mpRF35 mesa 0
run 83133 0.500000 mp2mpRF35 mesa 0
run 10086 0.500000 mp2mpRF35 mesa 0
run 85866 0.500000 mp2mpRF35 mesa 0
run 84222 0.500000 mp2mpRF35 mesa 0
run 86649 1.000000 mp2mpRF35 mesa 0
run 14435 1.000000 mp2mpRF35 mesa 0
run 14547 1.000000 mp2mpRF35 mesa 0
run 31605 1.000000 mp2mpRF35 mesa 0
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/menu-short.menu 0 → 100644
+ 102
− 0
View file @ f2d91ab7
## Generated at 14:41 on June 05 2020 by banerjee_p
# git version: muone (54e3690d)
image mcmule:muonelegacy ep2ep_mesa_paper/user.f95
conf ep2ep_mesa_paper/ep2ep-mesa.conf
run 47135 1.000000 mp2mpRRc mesa 0
run 47135 1.000000 mp2mpRco mesa 0
run 73614 1.000000 mp2mp0 mesa 0
run 93199 1.000000 mp2mp0 mesa 0
run 55896 1.000000 mp2mp0 mesa 0
run 17284 1.000000 mp2mp0 mesa 0
run 50392 0.100000 mp2mpF mesa 0
run 82493 0.100000 mp2mpF mesa 0
run 83269 0.100000 mp2mpF mesa 0
run 28460 0.100000 mp2mpF mesa 0
run 62300 0.300000 mp2mpF mesa 0
run 76287 0.300000 mp2mpF mesa 0
run 55853 0.300000 mp2mpF mesa 0
run 76464 0.300000 mp2mpF mesa 0
run 68184 0.500000 mp2mpF mesa 0
run 88231 0.500000 mp2mpF mesa 0
run 54662 0.500000 mp2mpF mesa 0
run 17599 0.500000 mp2mpF mesa 0
run 96104 1.000000 mp2mpF mesa 0
run 39305 1.000000 mp2mpF mesa 0
run 34334 1.000000 mp2mpF mesa 0
run 78597 1.000000 mp2mpF mesa 0
run 35001 0.100000 mp2mpR15 mesa 0
run 69217 0.100000 mp2mpR15 mesa 0
run 27998 0.100000 mp2mpR15 mesa 0
run 67037 0.100000 mp2mpR15 mesa 0
run 63714 0.300000 mp2mpR15 mesa 0
run 78601 0.300000 mp2mpR15 mesa 0
run 21938 0.300000 mp2mpR15 mesa 0
run 99548 0.300000 mp2mpR15 mesa 0
run 20460 0.500000 mp2mpR15 mesa 0
run 94059 0.500000 mp2mpR15 mesa 0
run 87913 0.500000 mp2mpR15 mesa 0
run 21136 0.500000 mp2mpR15 mesa 0
run 41959 1.000000 mp2mpR15 mesa 0
run 54802 1.000000 mp2mpR15 mesa 0
run 54697 1.000000 mp2mpR15 mesa 0
run 63357 1.000000 mp2mpR15 mesa 0
run 49902 0.100000 mp2mpR35 mesa 0
run 85297 0.100000 mp2mpR35 mesa 0
run 69496 0.100000 mp2mpR35 mesa 0
run 62975 0.100000 mp2mpR35 mesa 0
run 25576 0.300000 mp2mpR35 mesa 0
run 43007 0.300000 mp2mpR35 mesa 0
run 95931 0.300000 mp2mpR35 mesa 0
run 59453 0.300000 mp2mpR35 mesa 0
run 30657 0.500000 mp2mpR35 mesa 0
run 61979 0.500000 mp2mpR35 mesa 0
run 99325 0.500000 mp2mpR35 mesa 0
run 61250 0.500000 mp2mpR35 mesa 0
run 48649 1.000000 mp2mpR35 mesa 0
run 87126 1.000000 mp2mpR35 mesa 0
run 57304 1.000000 mp2mpR35 mesa 0
run 99400 1.000000 mp2mpR35 mesa 0
run 74712 0.100000 mp2mpFF mesa 0
run 89263 0.100000 mp2mpFF mesa 0
run 27348 0.100000 mp2mpFF mesa 0
run 31166 0.100000 mp2mpFF mesa 0
run 30987 0.300000 mp2mpFF mesa 0
run 29068 0.300000 mp2mpFF mesa 0
run 72964 0.300000 mp2mpFF mesa 0
run 97401 0.300000 mp2mpFF mesa 0
run 82374 0.500000 mp2mpFF mesa 0
run 33970 0.500000 mp2mpFF mesa 0
run 14277 0.500000 mp2mpFF mesa 0
run 42777 0.500000 mp2mpFF mesa 0
run 95443 1.000000 mp2mpFF mesa 0
run 55650 1.000000 mp2mpFF mesa 0
run 66479 1.000000 mp2mpFF mesa 0
run 88757 1.000000 mp2mpFF mesa 0
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/menu-vp.menu 0 → 100644
+ 142
− 0
View file @ f2d91ab7
image mcmule:muonelegacy ep2ep_mesa_paper/user-vp.f95
conf ep2ep_mesa_paper/ep2ep-mesa.conf
#A
run 88874 1.000000 mp2mpA mesa 0
run 47404 1.000000 mp2mpA mesa 0
run 17404 1.000000 mp2mpA mesa 0
run 93874 1.000000 mp2mpA mesa 1
run 97109 1.000000 mp2mpA mesa 1
run 17109 1.000000 mp2mpA mesa 1
#AA
run 82874 1.000000 mp2mpAA mesa 0
run 82109 1.000000 mp2mpAA mesa 0
run 12109 1.000000 mp2mpAA mesa 0
run 84871 1.000000 mp2mpAA mesa 1
run 84101 1.000000 mp2mpAA mesa 1
run 14101 1.000000 mp2mpAA mesa 1
#NF
run 42404 1.000000 mp2mpNF mesa 0
run 82109 1.000000 mp2mpNF mesa 0
run 12109 1.000000 mp2mpNF mesa 0
run 32474 1.000000 mp2mpNF mesa 1
run 42179 1.000000 mp2mpNF mesa 1
run 12179 1.000000 mp2mpNF mesa 1
run 42404 1.000000 mp2mpNF mesa 2
run 82109 1.000000 mp2mpNF mesa 2
run 12109 1.000000 mp2mpNF mesa 2
#AF
run 65011 0.100000 mp2mpAF mesa 0
run 25811 0.100000 mp2mpAF mesa 0
run 15811 0.100000 mp2mpAF mesa 0
run 51011 0.300000 mp2mpAF mesa 0
run 89611 0.300000 mp2mpAF mesa 0
run 19611 0.300000 mp2mpAF mesa 0
run 67511 0.500000 mp2mpAF mesa 0
run 63511 0.500000 mp2mpAF mesa 0
run 13511 0.500000 mp2mpAF mesa 0
run 94011 1.000000 mp2mpAF mesa 0
run 12911 1.000000 mp2mpAF mesa 0
run 32911 1.000000 mp2mpAF mesa 0
run 15011 0.100000 mp2mpAF mesa 1
run 15811 0.100000 mp2mpAF mesa 1
run 35811 0.100000 mp2mpAF mesa 1
run 11011 0.300000 mp2mpAF mesa 1
run 19611 0.300000 mp2mpAF mesa 1
run 39611 0.300000 mp2mpAF mesa 1
run 17511 0.500000 mp2mpAF mesa 1
run 13511 0.500000 mp2mpAF mesa 1
run 33511 0.500000 mp2mpAF mesa 1
run 14011 1.000000 mp2mpAF mesa 1
run 12911 1.000000 mp2mpAF mesa 1
run 32911 1.000000 mp2mpAF mesa 1
#AR15
run 28026 0.100000 mp2mpAR15 mesa 0
run 14014 0.100000 mp2mpAR15 mesa 0
run 34014 0.100000 mp2mpAR15 mesa 0
run 30193 0.300000 mp2mpAR15 mesa 0
run 72123 0.300000 mp2mpAR15 mesa 0
run 12123 0.300000 mp2mpAR15 mesa 0
run 12430 0.500000 mp2mpAR15 mesa 0
run 21526 0.500000 mp2mpAR15 mesa 0
run 31526 0.500000 mp2mpAR15 mesa 0
run 36226 1.000000 mp2mpAR15 mesa 0
run 90257 1.000000 mp2mpAR15 mesa 0
run 10257 1.000000 mp2mpAR15 mesa 0
run 98026 0.100000 mp2mpAR15 mesa 1
run 94014 0.100000 mp2mpAR15 mesa 1
run 14014 0.100000 mp2mpAR15 mesa 1
run 90193 0.300000 mp2mpAR15 mesa 1
run 92123 0.300000 mp2mpAR15 mesa 1
run 12123 0.300000 mp2mpAR15 mesa 1
run 92430 0.500000 mp2mpAR15 mesa 1
run 91526 0.500000 mp2mpAR15 mesa 1
run 11526 0.500000 mp2mpAR15 mesa 1
run 96226 1.000000 mp2mpAR15 mesa 1
run 90257 1.000000 mp2mpAR15 mesa 1
run 20257 1.000000 mp2mpAR15 mesa 1
#AR35
run 63997 0.100000 mp2mpAR35 mesa 0
run 21847 0.100000 mp2mpAR35 mesa 0
run 11847 0.100000 mp2mpAR35 mesa 0
run 71160 0.300000 mp2mpAR35 mesa 0
run 38322 0.300000 mp2mpAR35 mesa 0
run 18322 0.300000 mp2mpAR35 mesa 0
run 78520 0.500000 mp2mpAR35 mesa 0
run 32026 0.500000 mp2mpAR35 mesa 0
run 12026 0.500000 mp2mpAR35 mesa 0
run 15629 1.000000 mp2mpAR35 mesa 0
run 68351 1.000000 mp2mpAR35 mesa 0
run 18351 1.000000 mp2mpAR35 mesa 0
run 13997 0.100000 mp2mpAR35 mesa 1
run 11847 0.100000 mp2mpAR35 mesa 1
run 41847 0.100000 mp2mpAR35 mesa 1
run 11160 0.300000 mp2mpAR35 mesa 1
run 18322 0.300000 mp2mpAR35 mesa 1
run 48322 0.300000 mp2mpAR35 mesa 1
run 18520 0.500000 mp2mpAR35 mesa 1
run 12026 0.500000 mp2mpAR35 mesa 1
run 52026 0.500000 mp2mpAR35 mesa 1
run 15629 1.000000 mp2mpAR35 mesa 1
run 18351 1.000000 mp2mpAR35 mesa 1
run 28351 1.000000 mp2mpAR35 mesa 1
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/out.tar.bz2 0 → 100644 LFS
+ 3
− 0
View file @ f2d91ab7
File added
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/user-vp.f95 0 → 100644
+ 155
− 0
View file @ f2d91ab7
!!!!!!!!!!!!!!!!!!!!!
MODULE USER
!!!!!!!!!!!!!!!!!!!!!
use functions
implicit none
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
integer, parameter :: nr_q = 10
integer, parameter :: nr_bins =200
real, parameter :: &
min_val(nr_q) = (/ 25., 115., 0.6, 0.6, 25., 25.,25.,25., 4000., 4000. /)
real, parameter :: &
max_val(nr_q) = (/ 45., 155., 1., 1., 45., 45.,45., 45., 14000., 14000. /)
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!! ============================================== !!
!! DO NOT EVEN THINK ABOUT CHANGING ANYTHING HERE !!
!! ============================================== !!
logical :: pass_cut(nr_q)
character (len = 10), 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 ==== !!
mu = Me
musq = mu**2
END SUBROUTINE FIX_MU
SUBROUTINE INITUSER
integer :: nf
print*, "This is a MESA legacy run with the cuts:"
print*," 25 pi/180._prec< thetae < 45 pi/180._prec "
print*, " Ee(out) > 45 MeV "
read*, nf
if(nf==0) then
print*, "filenamesuffix=0 -> nf=nel+nmu+ntau"
nhad = 0._prec
elseif(nf==1) then
print*, "filenamesuffix=1 -> nf=nel+nmu+ntau+nhad"
elseif(nf==2) then
nmu = 0._prec
ntau = 0._prec
nhad = 0._prec
print*, "filenamesuffix=2 -> nf=nel"
endif
write(filenamesuffix,'(I1)') nf
END SUBROUTINE
FUNCTION QUANT(q1,q2,q3,q4,q5,q6,q7)
real (kind=prec), intent(in) :: q1(4),q2(4),q3(4),q4(4), q5(4),q6(4),q7(4)
real (kind=prec) :: q1lab(4),q2lab(4),q3lab(4),q4lab(4),q5lab(4),q6lab(4)
real (kind=prec) :: quant(nr_q)
real (kind=prec) :: thetae, cthetae, qsqE, qsqP, epsE, epsP
pol1 = (/ 0._prec, 0._prec, 0._prec, 0._prec /)
pass_cut = .true.
call fix_mu
q1lab = boost_rf(q2,q1) ! incoming electron
q2lab = boost_rf(q2,q2) ! proton at rest
q3lab = boost_rf(q2,q3) ! outgoing electron
q4lab = boost_rf(q2,q4) ! outgoing proton
q5lab = boost_rf(q2,q5) ! outgoing photon (if present)
q6lab = boost_rf(q2,q6) ! outgoing photon (if present)
cthetae = cos_th(q1lab,q3lab)
thetae = acos(cos_th(q1lab,q3lab))
qsqE = -sq(q1-q3)
qsqP = -sq(q2-q4)
epsE = 1./(1. + 2*(1 + qsqe/(4.*Mproton**2))*(tan(0.5*thetae))**2)
epsP = 1./(1. + 2*(1 + qsqp/(4.*Mproton**2))*(tan(0.5*thetae))**2)
if(thetae < 25*pi/180._prec) pass_cut=.false.
if(thetae > 45*pi/180._prec) pass_cut=.false.
if(q3lab(4) < 45.) pass_cut=.false.
if(abs(qsqE-8000) > 1000) pass_cut(5)=.false.
if(abs(qsqP-8000) > 1000) pass_cut(6)=.false.
if(abs(qsqE-8000) > 2000) pass_cut(7)=.false.
if(abs(qsqP-8000) > 2000) pass_cut(8)=.false.
names(1) = "thetae"
quant(1) = 180*thetae/pi
names(2) = "Ee"
quant(2) = q3lab(4)
names(3) = "epsE"
quant(3) = epse
names(4) = "epsP"
quant(4) = epsp
names(5) = "thetaeE1"
quant(5) = 180*thetae/pi
names(6) = "thetaeP1"
quant(6) = 180*thetae/pi
names(7) = "thetaeE2"
quant(7) = 180*thetae/pi
names(8) = "thetaeP2"
quant(8) = 180*thetae/pi
names(9) = "qsqE"
quant(9) = qsqe
names(10) = "qsqP"
quant(10) = qsqp
END FUNCTION QUANT
!!!!!!!!!!!!!!!!!!!!!!!
END MODULE USER
!!!!!!!!!!!!!!!!!!!!!!!
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/user.f95 0 → 100644
+ 140
− 0
View file @ f2d91ab7
!!!!!!!!!!!!!!!!!!!!!
MODULE USER
!!!!!!!!!!!!!!!!!!!!!
use functions
implicit none
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
integer, parameter :: nr_q = 10
integer, parameter :: nr_bins =200
real, parameter :: &
min_val(nr_q) = (/ 25., 115., 0.6, 0.6, 25., 25.,25.,25., 4000., 4000. /)
real, parameter :: &
max_val(nr_q) = (/ 45., 155., 1., 1., 45., 45.,45., 45., 14000., 14000. /)
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!! ============================================== !!
!! DO NOT EVEN THINK ABOUT CHANGING ANYTHING HERE !!
!! ============================================== !!
logical :: pass_cut(nr_q)
character (len = 10), 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 ==== !!
mu = Me
musq = mu**2
END SUBROUTINE FIX_MU
SUBROUTINE INITUSER
print*, "This is a MESA legacy run with the cuts:"
print*," 25 pi/180._prec< thetae < 45 pi/180._prec "
print*, " Ee(out) > 45 MeV "
END SUBROUTINE
FUNCTION QUANT(q1,q2,q3,q4,q5,q6,q7)
real (kind=prec), intent(in) :: q1(4),q2(4),q3(4),q4(4), q5(4),q6(4),q7(4)
real (kind=prec) :: q1lab(4),q2lab(4),q3lab(4),q4lab(4),q5lab(4),q6lab(4)
real (kind=prec) :: quant(nr_q)
real (kind=prec) :: thetae, cthetae, qsqE, qsqP, epsE, epsP
pol1 = (/ 0._prec, 0._prec, 0._prec, 0._prec /)
pass_cut = .true.
call fix_mu
q1lab = boost_rf(q2,q1) ! incoming electron
q2lab = boost_rf(q2,q2) ! proton at rest
q3lab = boost_rf(q2,q3) ! outgoing electron
q4lab = boost_rf(q2,q4) ! outgoing proton
q5lab = boost_rf(q2,q5) ! outgoing photon (if present)
q6lab = boost_rf(q2,q6) ! outgoing photon (if present)
cthetae = cos_th(q1lab,q3lab)
thetae = acos(cos_th(q1lab,q3lab))
qsqE = -sq(q1-q3)
qsqP = -sq(q2-q4)
epsE = 1./(1. + 2*(1 + qsqe/(4.*Mproton**2))*(tan(0.5*thetae))**2)
epsP = 1./(1. + 2*(1 + qsqp/(4.*Mproton**2))*(tan(0.5*thetae))**2)
if(thetae < 25*pi/180._prec) pass_cut=.false.
if(thetae > 45*pi/180._prec) pass_cut=.false.
if(q3lab(4) < 45.) pass_cut=.false.
if(abs(qsqE-8000) > 1000) pass_cut(5)=.false.
if(abs(qsqP-8000) > 1000) pass_cut(6)=.false.
if(abs(qsqE-8000) > 2000) pass_cut(7)=.false.
if(abs(qsqP-8000) > 2000) pass_cut(8)=.false.
names(1) = "thetae"
quant(1) = 180*thetae/pi
names(2) = "Ee"
quant(2) = q3lab(4)
names(3) = "epsE"
quant(3) = epse
names(4) = "epsP"
quant(4) = epsp
names(5) = "thetaeE1"
quant(5) = 180*thetae/pi
names(6) = "thetaeP1"
quant(6) = 180*thetae/pi
names(7) = "thetaeE2"
quant(7) = 180*thetae/pi
names(8) = "thetaeP2"
quant(8) = 180*thetae/pi
names(9) = "qsqE"
quant(9) = qsqe
names(10) = "qsqP"
quant(10) = qsqp
END FUNCTION QUANT
!!!!!!!!!!!!!!!!!!!!!!!
END MODULE USER
!!!!!!!!!!!!!!!!!!!!!!!
l-p-scattering/mesa-legacy/ep2ep_mesa_paper/workers.tar.xz 0 → 100644 LFS
+ 3
− 0
View file @ f2d91ab7
File added
l-p-scattering/mesa-legacy/legacy.py 0 → 100644
+ 103
− 0
View file @ f2d91ab7
# 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')
###########################################################}}}
#####################################################################}}}
l-p-scattering/mesa-legacy/makefile 0 → 100644
+ 5
− 0
View file @ f2d91ab7
group=l-p-scattering
name=legacy
plots=qsq.pdf thetae.pdf
figures=$(addprefix plots/,$(plots))
include ../../tools/makefile.conf
l-p-scattering/mesa-legacy/plots/.dir 0 → 100644
+ 0
− 0
View file @ f2d91ab7
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