Commit dc46f369 authored by ulrich_y's avatar ulrich_y
Browse files

Added MUSE e115 data

parent ac118f27
## Generated at 21:14 on October 10 2019 by ulrich_y
# git version: muse (0a1d7c2)
# To re-generate, run python tools/create-menu.py \
# --seeds 63869 64438 70353 39853 62257 65233 \
# -xi 0.01 0.03 0.1 0.222623 0.3 1.0 \
# --flavour musee \
# --piece mp2mp \
# --output-dir e115 \
# --prog xs_main \
# --stat FEE35,50000,20,250000,70 \
# --stat RREE,600000,20,2000000,70 \
# --stat FEE15,50000,20,250000,70 \
# --stat REE,20000,20,500000,70,5
conf e115/mp2mp-musee.conf
run 63869 0.010000 mp2mpRF35 musee 0908
run 64438 0.010000 mp2mpRF35 musee 0908
run 70353 0.010000 mp2mpRF35 musee 1008
run 39853 0.010000 mp2mpRF35 musee 1009
run 62257 0.010000 mp2mpRF35 musee 0909
run 65233 0.010000 mp2mpRF35 musee 0907
run 63869 0.030000 mp2mpRF35 musee 0908
run 64438 0.030000 mp2mpRF35 musee 0908
run 70353 0.030000 mp2mpRF35 musee 1008
run 39853 0.030000 mp2mpRF35 musee 1009
run 62257 0.030000 mp2mpRF35 musee 0909
run 65233 0.030000 mp2mpRF35 musee 0907
run 63869 0.100000 mp2mpRF35 musee 0908
run 64438 0.100000 mp2mpRF35 musee 0908
run 70353 0.100000 mp2mpRF35 musee 1008
run 39853 0.100000 mp2mpRF35 musee 1009
run 62257 0.100000 mp2mpRF35 musee 0909
run 65233 0.100000 mp2mpRF35 musee 0909
run 63869 0.222623 mp2mpRF35 musee 0908
run 64438 0.222623 mp2mpRF35 musee 0908
run 70353 0.222623 mp2mpRF35 musee 1008
run 39853 0.222623 mp2mpRF35 musee 1009
run 62257 0.222623 mp2mpRF35 musee 0909
run 65233 0.222623 mp2mpRF35 musee 0909
run 63869 0.300000 mp2mpRF35 musee 0908
run 64438 0.300000 mp2mpRF35 musee 0908
run 70353 0.300000 mp2mpRF35 musee 1008
run 39853 0.300000 mp2mpRF35 musee 1009
run 62257 0.300000 mp2mpRF35 musee 0909
run 65233 0.300000 mp2mpRF35 musee 0909
run 63869 1.000000 mp2mpRF35 musee 0908
run 64438 1.000000 mp2mpRF35 musee 0908
run 70353 1.000000 mp2mpRF35 musee 1008
run 39853 1.000000 mp2mpRF35 musee 1009
run 62257 1.000000 mp2mpRF35 musee 0909
run 65233 1.000000 mp2mpRF35 musee 0909
run 63869 0.010000 mp2mpRR musee 0010
run 64438 0.010000 mp2mpRR musee 0010
run 70353 0.010000 mp2mpRR musee 0010
run 39853 0.010000 mp2mpRR musee 0009
run 62257 0.010000 mp2mpRR musee 0009
run 65233 0.010000 mp2mpRR musee 0009
run 63869 0.030000 mp2mpRR musee 0010
run 64438 0.030000 mp2mpRR musee 0010
run 70353 0.030000 mp2mpRR musee 0010
run 39853 0.030000 mp2mpRR musee 0009
run 62257 0.030000 mp2mpRR musee 0009
run 65233 0.030000 mp2mpRR musee 0009
run 63869 0.100000 mp2mpRR musee 0010
run 64438 0.100000 mp2mpRR musee 0010
run 70353 0.100000 mp2mpRR musee 0010
run 39853 0.100000 mp2mpRR musee 0009
run 62257 0.100000 mp2mpRR musee 0009
run 65233 0.100000 mp2mpRR musee 0009
run 63869 0.222623 mp2mpRR musee 0010
run 64438 0.222623 mp2mpRR musee 0010
run 70353 0.222623 mp2mpRR musee 0010
run 39853 0.222623 mp2mpRR musee 0009
run 62257 0.222623 mp2mpRR musee 0009
run 65233 0.222623 mp2mpRR musee 0009
run 63869 0.300000 mp2mpRR musee 0010
run 64438 0.300000 mp2mpRR musee 0010
run 70353 0.300000 mp2mpRR musee 0010
run 39853 0.300000 mp2mpRR musee 0009
run 62257 0.300000 mp2mpRR musee 0009
run 65233 0.300000 mp2mpRR musee 0009
run 63869 1.000000 mp2mpRR musee 0010
run 64438 1.000000 mp2mpRR musee 0010
run 70353 1.000000 mp2mpRR musee 0010
run 39853 1.000000 mp2mpRR musee 0009
run 62257 1.000000 mp2mpRR musee 0009
run 65233 1.000000 mp2mpRR musee 0009
run 63869 0.010000 mp2mpRF15 musee 0908
run 64438 0.010000 mp2mpRF15 musee 0908
run 70353 0.010000 mp2mpRF15 musee 1008
run 39853 0.010000 mp2mpRF15 musee 1009
run 62257 0.010000 mp2mpRF15 musee 0909
run 65233 0.010000 mp2mpRF15 musee 0907
run 63869 0.030000 mp2mpRF15 musee 0908
run 64438 0.030000 mp2mpRF15 musee 0908
run 70353 0.030000 mp2mpRF15 musee 1008
run 39853 0.030000 mp2mpRF15 musee 1009
run 62257 0.030000 mp2mpRF15 musee 0909
run 65233 0.030000 mp2mpRF15 musee 0907
run 63869 0.100000 mp2mpRF15 musee 0908
run 64438 0.100000 mp2mpRF15 musee 0908
run 70353 0.100000 mp2mpRF15 musee 1008
run 39853 0.100000 mp2mpRF15 musee 1009
run 62257 0.100000 mp2mpRF15 musee 0909
run 65233 0.100000 mp2mpRF15 musee 0909
run 63869 0.222623 mp2mpRF15 musee 0908
run 64438 0.222623 mp2mpRF15 musee 0908
run 70353 0.222623 mp2mpRF15 musee 1008
run 39853 0.222623 mp2mpRF15 musee 1009
run 62257 0.222623 mp2mpRF15 musee 0909
run 65233 0.222623 mp2mpRF15 musee 0909
run 63869 0.300000 mp2mpRF15 musee 0908
run 64438 0.300000 mp2mpRF15 musee 0908
run 70353 0.300000 mp2mpRF15 musee 1008
run 39853 0.300000 mp2mpRF15 musee 1009
run 62257 0.300000 mp2mpRF15 musee 0909
run 65233 0.300000 mp2mpRF15 musee 0909
run 63869 1.000000 mp2mpRF15 musee 0908
run 64438 1.000000 mp2mpRF15 musee 0908
run 70353 1.000000 mp2mpRF15 musee 1008
run 39853 1.000000 mp2mpRF15 musee 1009
run 62257 1.000000 mp2mpRF15 musee 0909
run 65233 1.000000 mp2mpRF15 musee 0909
run 63869 0.010000 mp2mpR musee 0
run 64438 0.010000 mp2mpR musee 0
run 70353 0.010000 mp2mpR musee 0
run 39853 0.010000 mp2mpR musee 0
run 62257 0.010000 mp2mpR musee 0
run 63869 0.030000 mp2mpR musee 0
run 64438 0.030000 mp2mpR musee 0
run 70353 0.030000 mp2mpR musee 0
run 39853 0.030000 mp2mpR musee 0
run 62257 0.030000 mp2mpR musee 0
run 63869 0.100000 mp2mpR musee 0
run 64438 0.100000 mp2mpR musee 0
run 70353 0.100000 mp2mpR musee 0
run 39853 0.100000 mp2mpR musee 0
run 62257 0.100000 mp2mpR musee 0
run 63869 0.222623 mp2mpR musee 0
run 64438 0.222623 mp2mpR musee 0
run 70353 0.222623 mp2mpR musee 0
run 39853 0.222623 mp2mpR musee 0
run 62257 0.222623 mp2mpR musee 0
run 63869 0.300000 mp2mpR musee 0
run 64438 0.300000 mp2mpR musee 0
run 70353 0.300000 mp2mpR musee 0
run 39853 0.300000 mp2mpR musee 0
run 62257 0.300000 mp2mpR musee 0
run 63869 1.000000 mp2mpR musee 0
run 64438 1.000000 mp2mpR musee 0
run 70353 1.000000 mp2mpR musee 0
run 39853 1.000000 mp2mpR musee 0
run 62257 1.000000 mp2mpR musee 0
## Generated at 21:14 on October 10 2019 by ulrich_y
# git version: muse (0a1d7c2)
# To re-generate, run python tools/create-menu.py \
# --seeds 63869 64438 70353 39853 62257 65233 \
# -xi 0.01 0.03 0.1 0.222623 0.3 1.0 \
# --flavour musee \
# --piece em2em \
# --output-dir e115 \
# --prog xs_main \
# --stat FEE35,50000,20,250000,70 \
# --stat RREE,600000,20,2000000,70 \
# --stat FEE15,50000,20,250000,70 \
# --stat REE,20000,20,500000,70,5
# specify the program to run relative to `pwd`
binary=e115/xs_main
# specify the output folder
folder=e115/
# Specify the variables nenter_ad, itmx_ad, nenter and itmx
# for each piece you want to run.
declare -A STAT=(
["mp2mpRF35"]="50000\n20\n250000\n70"
["mp2mpRF15"]="50000\n20\n250000\n70"
["mp2mpRR"]="600000\n20\n2000000\n70"
["mp2mpR"]="20000\n20\n500000\n70"
)
!!!!!!!!!!!!!!!!!!!!!
MODULE USER
!!!!!!!!!!!!!!!!!!!!!
use functions
implicit none
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
integer, parameter :: nr_q = 2
integer, parameter :: nr_bins = 200
real, parameter :: &
min_val(nr_q) = (/ 1./9.*pi , 0. /)
real, parameter :: &
max_val(nr_q) = (/ 5./9.*pi , 200. /)
real(kind=prec) :: hcut,ycut
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!==!
!! ============================================== !!
!! 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
! This is called without arguments once as soon as McMule
! starts and has read all other configuration, i.e. you can
! access which_piece and flavour. Use this to read any
! further information from the user (like cut configuration
! etc). You do not have to print the hashes - this is
! already taken care of - but you are very much invited to
! include information of what it is you are doing
!
! If you are using the cut channel of the menu, you may need to set
! the filenamesuffix variable which is appended to the name of the
! VEGAS file.
! Example for reading a cut:
! integer cut
integer inp
read*, inp
write(filenamesuffix,'(I4)') inp
hcut = 10._prec**(-real(mod(inp, 100)))
ycut = 10._prec**(-real(int(inp/ 100)))
print*,"using hcut=",hcut,"and ycut=",ycut
if (which_piece == 'em2emRREE') then
print*,"Ignoring ycut though"
endif
END SUBROUTINE
FUNCTION QUANT(q1,q2,q3,q4,q5,q6,q7)
real (kind=prec), intent(in) :: q1(4),q2(4),q3(4),q4(4)
real (kind=prec), optional :: q5(4),q6(4),q7(4)
real (kind=prec) :: qq5(4),qq6(4),qq7(4),q1Rest(4),q2Rest(4),q3Rest(4),q4Rest(4)
real (kind=prec) :: quant(nr_q)
real (kind=prec) :: thetal, qsq
qq5 = 0._prec
qq6 = 0._prec
qq7 = 0._prec
if(present(q5)) qq5=q5
if(present(q6)) qq6=q6
if(present(q7)) qq7=q7
! pol1 = (/ 0._prec, 0._prec, 0.85_prec, 0._prec /)
pol1 = (/ 0._prec, 0._prec, 0._prec, 0._prec /)
!! ==== keep the line below in any case ==== !!
call fix_mu
q1Rest = boost_rf(q2,q1)
q2Rest = boost_rf(q2,q2)
q3Rest = boost_rf(q2,q3)
q4Rest = boost_rf(q2,q4)
thetal = acos(cos_th(q1Rest,q3Rest))
qsq = -sq(q2-q4)
set_zero = 1
if (flavour=="muse") then
if (thetal.lt. pi/9.) set_zero=0
if (thetal.gt.5*pi/9.) set_zero=0
elseif (flavour=="musee") then
if (thetal.lt. pi/9.) set_zero=0
if (thetal.gt.5*pi/9.) set_zero=0
elseif(flavour=="mesa") then
if (thetal.lt.5./36.*pi) set_zero=0
if (thetal.gt.1./ 4.*pi) set_zero=0
if (q3rest(4) .lt. 45.) set_zero=0
endif
names(1) = "angle"
quant(1) = thetal
names(2) = "energy"
quant(2) = q3rest(4)
! move variable out of bounds if cut has been applied
quant = quant + 2*max_val*(1-set_zero)
END FUNCTION QUANT
!!!!!!!!!!!!!!!!!!!!!!!
END MODULE USER
!!!!!!!!!!!!!!!!!!!!!!!
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