Commit d4fcd00b authored by ulrich_y's avatar ulrich_y
Browse files

Add 2D data for MEG

parent c99d3bfe
## Generated at 12:03 on February 07 2020 by ulrich_y
# git version: master (423084d)
# To re-generate, run python tools/create-menu.py \
# --seeds 48239 28709 68031 44313 50699 89382 90658 13285 \
# -xi 0.1 0.3 0.5 \
# --flavour mu-e \
# --piece m2enn \
# --output-dir meg-2d-nnlo \
# --prog xs_main \
# --stat RF,10000,20,100000,100 \
# --stat RR,20000,50,150000,100 \
# --stat 0,20000,20,200000,150 \
# --stat R,20000,20,200000,200 \
# --stat F,20000,20,200000,150
# specify the program to run relative to `pwd`
binary=meg-2d/xs_main
# specify the output folder
folder=meg-2d/
# Specify the variables nenter_ad, itmx_ad, nenter and itmx
# for each piece you want to run.
declare -A STAT=(
["m2enn0"]="20000\n20\n200000\n150"
["m2ennR"]="20000\n20\n200000\n200"
["m2ennF"]="20000\n20\n200000\n150"
["m2ennRF"]="10000\n20\n100000\n80"
["m2ennRR"]="20000\n50\n150000\n80"
)
## Generated at 16:21 on February 05 2020 by gurgone_a
# git version: master (423084d)
# To re-generate, run python tools/create-menu.py \
# --seeds 26351 74892 \
# -xi 0.5 0.25 \
# --flavour mu-e \
# --piece m2enn \
# --output-dir meg-2d \
# --prog xs_main \
# --stat 0,20000,20,200000,150 \
# --stat R,20000,20,200000,200 \
# --stat F,20000,20,200000,150
conf meg-2d/m2enn-mu-e.conf
run 26351 1.000000 m2enn0 mu-e 11
run 26351 1.000000 m2enn0 mu-e 21
run 26351 1.000000 m2enn0 mu-e 12
run 26351 1.000000 m2enn0 mu-e 22
run 26351 1.000000 m2enn0 mu-e 13
run 26351 1.000000 m2enn0 mu-e 23
run 26351 1.000000 m2enn0 mu-e 14
run 26351 1.000000 m2enn0 mu-e 24
run 74892 1.000000 m2enn0 mu-e 11
run 74892 1.000000 m2enn0 mu-e 21
run 74892 1.000000 m2enn0 mu-e 12
run 74892 1.000000 m2enn0 mu-e 22
run 74892 1.000000 m2enn0 mu-e 13
run 74892 1.000000 m2enn0 mu-e 23
run 74892 1.000000 m2enn0 mu-e 14
run 74892 1.000000 m2enn0 mu-e 24
run 26351 0.500000 m2ennR mu-e 11
run 26351 0.500000 m2ennR mu-e 21
run 26351 0.500000 m2ennR mu-e 12
run 26351 0.500000 m2ennR mu-e 22
run 26351 0.500000 m2ennR mu-e 13
run 26351 0.500000 m2ennR mu-e 23
run 26351 0.500000 m2ennR mu-e 14
run 26351 0.500000 m2ennR mu-e 24
run 74892 0.500000 m2ennR mu-e 11
run 74892 0.500000 m2ennR mu-e 21
run 74892 0.500000 m2ennR mu-e 12
run 74892 0.500000 m2ennR mu-e 22
run 74892 0.500000 m2ennR mu-e 13
run 74892 0.500000 m2ennR mu-e 23
run 74892 0.500000 m2ennR mu-e 14
run 74892 0.500000 m2ennR mu-e 24
run 26351 0.250000 m2ennR mu-e 11
run 26351 0.250000 m2ennR mu-e 21
run 26351 0.250000 m2ennR mu-e 12
run 26351 0.250000 m2ennR mu-e 22
run 26351 0.250000 m2ennR mu-e 13
run 26351 0.250000 m2ennR mu-e 23
run 26351 0.250000 m2ennR mu-e 14
run 26351 0.250000 m2ennR mu-e 24
run 74892 0.250000 m2ennR mu-e 11
run 74892 0.250000 m2ennR mu-e 21
run 74892 0.250000 m2ennR mu-e 12
run 74892 0.250000 m2ennR mu-e 22
run 74892 0.250000 m2ennR mu-e 13
run 74892 0.250000 m2ennR mu-e 23
run 74892 0.250000 m2ennR mu-e 14
run 74892 0.250000 m2ennR mu-e 24
run 26351 0.500000 m2ennF mu-e 11
run 26351 0.500000 m2ennF mu-e 21
run 26351 0.500000 m2ennF mu-e 12
run 26351 0.500000 m2ennF mu-e 22
run 26351 0.500000 m2ennF mu-e 13
run 26351 0.500000 m2ennF mu-e 23
run 26351 0.500000 m2ennF mu-e 14
run 26351 0.500000 m2ennF mu-e 24
run 74892 0.500000 m2ennF mu-e 11
run 74892 0.500000 m2ennF mu-e 21
run 74892 0.500000 m2ennF mu-e 12
run 74892 0.500000 m2ennF mu-e 22
run 74892 0.500000 m2ennF mu-e 13
run 74892 0.500000 m2ennF mu-e 23
run 74892 0.500000 m2ennF mu-e 14
run 74892 0.500000 m2ennF mu-e 24
run 26351 0.250000 m2ennF mu-e 11
run 26351 0.250000 m2ennF mu-e 21
run 26351 0.250000 m2ennF mu-e 12
run 26351 0.250000 m2ennF mu-e 22
run 26351 0.250000 m2ennF mu-e 13
run 26351 0.250000 m2ennF mu-e 23
run 26351 0.250000 m2ennF mu-e 14
run 26351 0.250000 m2ennF mu-e 24
run 74892 0.250000 m2ennF mu-e 11
run 74892 0.250000 m2ennF mu-e 21
run 74892 0.250000 m2ennF mu-e 12
run 74892 0.250000 m2ennF mu-e 22
run 74892 0.250000 m2ennF mu-e 13
run 74892 0.250000 m2ennF mu-e 23
run 74892 0.250000 m2ennF mu-e 14
run 74892 0.250000 m2ennF mu-e 24
## Generated at 12:03 on February 07 2020 by ulrich_y
# git version: master (423084d)
# To re-generate, run python tools/create-menu.py \
# --seeds 48239 28709 68031 44313 50699 89382 90658 13285 \
# -xi 0.1 0.3 0.5 \
# --flavour mu-e \
# --piece m2enn \
# --output-dir meg-2d-nnlo \
# --prog xs_main \
# --stat RF,10000,20,100000,100 \
# --stat RR,20000,50,150000,100
conf meg-2d-nnlo/m2enn-mu-e.conf
run 48239 0.100000 m2ennRF mu-e 11
run 28709 0.100000 m2ennRF mu-e 11
run 68031 0.100000 m2ennRF mu-e 11
run 44313 0.100000 m2ennRF mu-e 11
run 50699 0.100000 m2ennRF mu-e 11
run 48239 0.300000 m2ennRF mu-e 11
run 28709 0.300000 m2ennRF mu-e 11
run 68031 0.300000 m2ennRF mu-e 11
run 44313 0.300000 m2ennRF mu-e 11
run 50699 0.300000 m2ennRF mu-e 11
run 48239 0.500000 m2ennRF mu-e 11
run 28709 0.500000 m2ennRF mu-e 11
run 68031 0.500000 m2ennRF mu-e 11
run 44313 0.500000 m2ennRF mu-e 11
run 50699 0.500000 m2ennRF mu-e 11
run 48239 0.100000 m2ennRR mu-e 11
run 28709 0.100000 m2ennRR mu-e 11
run 68031 0.100000 m2ennRR mu-e 11
run 44313 0.100000 m2ennRR mu-e 11
run 50699 0.100000 m2ennRR mu-e 11
run 48239 0.300000 m2ennRR mu-e 11
run 28709 0.300000 m2ennRR mu-e 11
run 68031 0.300000 m2ennRR mu-e 11
run 44313 0.300000 m2ennRR mu-e 11
run 50699 0.300000 m2ennRR mu-e 11
run 48239 0.500000 m2ennRR mu-e 11
run 28709 0.500000 m2ennRR mu-e 11
run 68031 0.500000 m2ennRR mu-e 11
run 44313 0.500000 m2ennRR mu-e 11
run 50699 0.500000 m2ennRR mu-e 11
run 48239 0.100000 m2ennRF mu-e 21
run 28709 0.100000 m2ennRF mu-e 21
run 68031 0.100000 m2ennRF mu-e 21
run 44313 0.100000 m2ennRF mu-e 21
run 50699 0.100000 m2ennRF mu-e 21
run 48239 0.300000 m2ennRF mu-e 21
run 28709 0.300000 m2ennRF mu-e 21
run 68031 0.300000 m2ennRF mu-e 21
run 44313 0.300000 m2ennRF mu-e 21
run 50699 0.300000 m2ennRF mu-e 21
run 48239 0.500000 m2ennRF mu-e 21
run 28709 0.500000 m2ennRF mu-e 21
run 68031 0.500000 m2ennRF mu-e 21
run 44313 0.500000 m2ennRF mu-e 21
run 50699 0.500000 m2ennRF mu-e 21
run 48239 0.100000 m2ennRR mu-e 21
run 28709 0.100000 m2ennRR mu-e 21
run 68031 0.100000 m2ennRR mu-e 21
run 44313 0.100000 m2ennRR mu-e 21
run 50699 0.100000 m2ennRR mu-e 21
run 48239 0.300000 m2ennRR mu-e 21
run 28709 0.300000 m2ennRR mu-e 21
run 68031 0.300000 m2ennRR mu-e 21
run 44313 0.300000 m2ennRR mu-e 21
run 50699 0.300000 m2ennRR mu-e 21
run 48239 0.500000 m2ennRR mu-e 21
run 28709 0.500000 m2ennRR mu-e 21
run 68031 0.500000 m2ennRR mu-e 21
run 44313 0.500000 m2ennRR mu-e 21
run 50699 0.500000 m2ennRR mu-e 21
run 48239 0.100000 m2ennRF mu-e 12
run 28709 0.100000 m2ennRF mu-e 12
run 68031 0.100000 m2ennRF mu-e 12
run 44313 0.100000 m2ennRF mu-e 12
run 50699 0.100000 m2ennRF mu-e 12
run 48239 0.300000 m2ennRF mu-e 12
run 28709 0.300000 m2ennRF mu-e 12
run 68031 0.300000 m2ennRF mu-e 12
run 44313 0.300000 m2ennRF mu-e 12
run 50699 0.300000 m2ennRF mu-e 12
run 48239 0.500000 m2ennRF mu-e 12
run 28709 0.500000 m2ennRF mu-e 12
run 68031 0.500000 m2ennRF mu-e 12
run 44313 0.500000 m2ennRF mu-e 12
run 50699 0.500000 m2ennRF mu-e 12
run 48239 0.100000 m2ennRR mu-e 12
run 28709 0.100000 m2ennRR mu-e 12
run 68031 0.100000 m2ennRR mu-e 12
run 44313 0.100000 m2ennRR mu-e 12
run 50699 0.100000 m2ennRR mu-e 12
run 48239 0.300000 m2ennRR mu-e 12
run 28709 0.300000 m2ennRR mu-e 12
run 68031 0.300000 m2ennRR mu-e 12
run 44313 0.300000 m2ennRR mu-e 12
run 50699 0.300000 m2ennRR mu-e 12
run 48239 0.500000 m2ennRR mu-e 12
run 28709 0.500000 m2ennRR mu-e 12
run 68031 0.500000 m2ennRR mu-e 12
run 44313 0.500000 m2ennRR mu-e 12
run 50699 0.500000 m2ennRR mu-e 12
run 48239 0.100000 m2ennRF mu-e 22
run 28709 0.100000 m2ennRF mu-e 22
run 68031 0.100000 m2ennRF mu-e 22
run 44313 0.100000 m2ennRF mu-e 22
run 50699 0.100000 m2ennRF mu-e 22
run 48239 0.300000 m2ennRF mu-e 22
run 28709 0.300000 m2ennRF mu-e 22
run 68031 0.300000 m2ennRF mu-e 22
run 44313 0.300000 m2ennRF mu-e 22
run 50699 0.300000 m2ennRF mu-e 22
run 48239 0.500000 m2ennRF mu-e 22
run 28709 0.500000 m2ennRF mu-e 22
run 68031 0.500000 m2ennRF mu-e 22
run 44313 0.500000 m2ennRF mu-e 22
run 50699 0.500000 m2ennRF mu-e 22
run 48239 0.100000 m2ennRR mu-e 22
run 28709 0.100000 m2ennRR mu-e 22
run 68031 0.100000 m2ennRR mu-e 22
run 44313 0.100000 m2ennRR mu-e 22
run 50699 0.100000 m2ennRR mu-e 22
run 48239 0.300000 m2ennRR mu-e 22
run 28709 0.300000 m2ennRR mu-e 22
run 68031 0.300000 m2ennRR mu-e 22
run 44313 0.300000 m2ennRR mu-e 22
run 50699 0.300000 m2ennRR mu-e 22
run 48239 0.500000 m2ennRR mu-e 22
run 28709 0.500000 m2ennRR mu-e 22
run 68031 0.500000 m2ennRR mu-e 22
run 44313 0.500000 m2ennRR mu-e 22
run 50699 0.500000 m2ennRR mu-e 22
run 48239 0.100000 m2ennRF mu-e 13
run 28709 0.100000 m2ennRF mu-e 13
run 68031 0.100000 m2ennRF mu-e 13
run 44313 0.100000 m2ennRF mu-e 13
run 50699 0.100000 m2ennRF mu-e 13
run 48239 0.300000 m2ennRF mu-e 13
run 28709 0.300000 m2ennRF mu-e 13
run 68031 0.300000 m2ennRF mu-e 13
run 44313 0.300000 m2ennRF mu-e 13
run 50699 0.300000 m2ennRF mu-e 13
run 48239 0.500000 m2ennRF mu-e 13
run 28709 0.500000 m2ennRF mu-e 13
run 68031 0.500000 m2ennRF mu-e 13
run 44313 0.500000 m2ennRF mu-e 13
run 50699 0.500000 m2ennRF mu-e 13
run 48239 0.100000 m2ennRR mu-e 13
run 28709 0.100000 m2ennRR mu-e 13
run 68031 0.100000 m2ennRR mu-e 13
run 44313 0.100000 m2ennRR mu-e 13
run 50699 0.100000 m2ennRR mu-e 13
run 48239 0.300000 m2ennRR mu-e 13
run 28709 0.300000 m2ennRR mu-e 13
run 68031 0.300000 m2ennRR mu-e 13
run 44313 0.300000 m2ennRR mu-e 13
run 50699 0.300000 m2ennRR mu-e 13
run 48239 0.500000 m2ennRR mu-e 13
run 28709 0.500000 m2ennRR mu-e 13
run 68031 0.500000 m2ennRR mu-e 13
run 44313 0.500000 m2ennRR mu-e 13
run 50699 0.500000 m2ennRR mu-e 13
run 48239 0.100000 m2ennRF mu-e 23
run 28709 0.100000 m2ennRF mu-e 23
run 68031 0.100000 m2ennRF mu-e 23
run 44313 0.100000 m2ennRF mu-e 23
run 50699 0.100000 m2ennRF mu-e 23
run 48239 0.300000 m2ennRF mu-e 23
run 28709 0.300000 m2ennRF mu-e 23
run 68031 0.300000 m2ennRF mu-e 23
run 44313 0.300000 m2ennRF mu-e 23
run 50699 0.300000 m2ennRF mu-e 23
run 48239 0.500000 m2ennRF mu-e 23
run 28709 0.500000 m2ennRF mu-e 23
run 68031 0.500000 m2ennRF mu-e 23
run 44313 0.500000 m2ennRF mu-e 23
run 50699 0.500000 m2ennRF mu-e 23
run 48239 0.100000 m2ennRR mu-e 23
run 28709 0.100000 m2ennRR mu-e 23
run 68031 0.100000 m2ennRR mu-e 23
run 44313 0.100000 m2ennRR mu-e 23
run 50699 0.100000 m2ennRR mu-e 23
run 48239 0.300000 m2ennRR mu-e 23
run 28709 0.300000 m2ennRR mu-e 23
run 68031 0.300000 m2ennRR mu-e 23
run 44313 0.300000 m2ennRR mu-e 23
run 50699 0.300000 m2ennRR mu-e 23
run 48239 0.500000 m2ennRR mu-e 23
run 28709 0.500000 m2ennRR mu-e 23
run 68031 0.500000 m2ennRR mu-e 23
run 44313 0.500000 m2ennRR mu-e 23
run 50699 0.500000 m2ennRR mu-e 23
run 48239 0.100000 m2ennRF mu-e 14
run 28709 0.100000 m2ennRF mu-e 14
run 68031 0.100000 m2ennRF mu-e 14
run 44313 0.100000 m2ennRF mu-e 14
run 50699 0.100000 m2ennRF mu-e 14
run 48239 0.300000 m2ennRF mu-e 14
run 28709 0.300000 m2ennRF mu-e 14
run 68031 0.300000 m2ennRF mu-e 14
run 44313 0.300000 m2ennRF mu-e 14
run 50699 0.300000 m2ennRF mu-e 14
run 48239 0.500000 m2ennRF mu-e 14
run 28709 0.500000 m2ennRF mu-e 14
run 68031 0.500000 m2ennRF mu-e 14
run 44313 0.500000 m2ennRF mu-e 14
run 50699 0.500000 m2ennRF mu-e 14
run 48239 0.100000 m2ennRR mu-e 14
run 28709 0.100000 m2ennRR mu-e 14
run 68031 0.100000 m2ennRR mu-e 14
run 44313 0.100000 m2ennRR mu-e 14
run 50699 0.100000 m2ennRR mu-e 14
run 48239 0.300000 m2ennRR mu-e 14
run 28709 0.300000 m2ennRR mu-e 14
run 68031 0.300000 m2ennRR mu-e 14
run 44313 0.300000 m2ennRR mu-e 14
run 50699 0.300000 m2ennRR mu-e 14
run 48239 0.500000 m2ennRR mu-e 14
run 28709 0.500000 m2ennRR mu-e 14
run 68031 0.500000 m2ennRR mu-e 14
run 44313 0.500000 m2ennRR mu-e 14
run 50699 0.500000 m2ennRR mu-e 14
run 48239 0.100000 m2ennRF mu-e 24
run 28709 0.100000 m2ennRF mu-e 24
run 68031 0.100000 m2ennRF mu-e 24
run 44313 0.100000 m2ennRF mu-e 24
run 50699 0.100000 m2ennRF mu-e 24
run 48239 0.300000 m2ennRF mu-e 24
run 28709 0.300000 m2ennRF mu-e 24
run 68031 0.300000 m2ennRF mu-e 24
run 44313 0.300000 m2ennRF mu-e 24
run 50699 0.300000 m2ennRF mu-e 24
run 48239 0.500000 m2ennRF mu-e 24
run 28709 0.500000 m2ennRF mu-e 24
run 68031 0.500000 m2ennRF mu-e 24
run 44313 0.500000 m2ennRF mu-e 24
run 50699 0.500000 m2ennRF mu-e 24
run 48239 0.100000 m2ennRR mu-e 24
run 28709 0.100000 m2ennRR mu-e 24
run 68031 0.100000 m2ennRR mu-e 24
run 44313 0.100000 m2ennRR mu-e 24
run 50699 0.100000 m2ennRR mu-e 24
run 48239 0.300000 m2ennRR mu-e 24
run 28709 0.300000 m2ennRR mu-e 24
run 68031 0.300000 m2ennRR mu-e 24
run 44313 0.300000 m2ennRR mu-e 24
run 50699 0.300000 m2ennRR mu-e 24
run 48239 0.500000 m2ennRR mu-e 24
run 28709 0.500000 m2ennRR mu-e 24
run 68031 0.500000 m2ennRR mu-e 24
run 44313 0.500000 m2ennRR mu-e 24
run 50699 0.500000 m2ennRR mu-e 24
!!!!!!!!!!!!!!!!!!!!!
MODULE USER
!!!!!!!!!!!!!!!!!!!!!
use functions
implicit none
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
integer, parameter :: nr_q = 5
integer, parameter :: nr_bins = 1000
integer :: whatplot
real, parameter :: &
min_val(nr_q) = (/ 0., 26., 42., 50., -1. /)
real, parameter :: &
max_val(nr_q) = (/ 26., 42., 50., 54., +1. /)
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!! ============================================== !!
!! DO NOT EVEN THINK ABOUT CHANGING ANYTHING HERE !!
!! ============================================== !!
integer :: set_zero(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 !!
!! ----------------------------------------- !!
real (kind=prec) :: bl_div = 1._prec
integer :: bin_flag = 0 !! 0 for standard; +1 for combined;
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 = Mm
musq = mu**2
END SUBROUTINE FIX_MU
SUBROUTINE USERINIT
integer :: p_encrypted
read*, p_encrypted
write(filenamesuffix,'(I2)') p_encrypted
! 1st digit is cos sign (1 = -, 2 = +)
! 2nd digit is the energy range
whatplot = mod(p_encrypted,10)
if(p_encrypted/10 == 1) then
whatplot = - whatplot
endif
END SUBROUTINE
FUNCTION QUANT(q1,q2,q3,q4,q5,q6,q7)
real (kind=prec), intent(in) :: q1(4),q2(4),q3(4)
real (kind=prec), optional :: q4(4),q5(4),q6(4),q7(4)
real (kind=prec) :: qq4(4),qq5(4),qq6(4),qq7(4)
real (kind=prec) :: ga(4), gah(4), gas(4)
real (kind=prec) :: quant(nr_q),ez(4)
qq4 = 0._prec
qq5 = 0._prec
qq6 = 0._prec
qq7 = 0._prec
if(present(q4)) qq4=q4
if(present(q5)) qq5=q5
if(present(q6)) qq6=q6
if(present(q7)) qq7=q7
!! ==== keep the line below in any case ==== !!
set_zero = 1
call fix_mu
! Polarization vector
pol1 = (/ 0._prec, 0._prec, 0.85_prec, 0._prec /)
ez = (/ 0._prec,0._prec, 1._prec,0._prec /)
! Energy range selection
if(abs(whatplot) == 1) then
if(q2(4) > 26.) set_zero = 0
elseif(abs(whatplot) == 2) then
if(q2(4) < 26. .or. q2(4) > 42.) set_zero = 0
elseif(abs(whatplot) == 3) then
if(q2(4) < 42. .or. q2(4) > 50.) set_zero = 0
elseif(abs(whatplot) == 4) then
if(q2(4) < 50.) set_zero = 0
endif
if (whatplot > 0) then
! We are in the cos > 0 region
if (cos_th(q2,ez) < 0) set_zero = 0
else
! We are in the cos < 0 region
if (cos_th(q2,ez) > 0) set_zero = 0
endif
! Observables
names(1) = "Ee1"
quant(1) = q2(4)
names(2) = "Ee2"
quant(2) = q2(4)
names(3) = "Ee3"
quant(3) = q2(4)
names(4) = "Ee4"
quant(4) = q2(4)
names(5) = "cthe"
quant(5) = cos_th(q2,ez)
! move variable out of bounds if cut has been applied
quant = quant + 2*max_val*(1-<