# This is a bash script that compile and runs all the MC@NLO codes for
# vector boson pair production. On your system, you need:
#
#    bash shell      AND      gmake
#
# which are rather standard (ask your system manager if they are not
# available).
#
# HOW TO USE THIS SCRIPT:
#  Look for "physical parameters" and "other input parameters" in
#  in this file; they control all the inputs for the MC@NLO codes.
#  After having modified them to suit your needs, execute this
#  file from a bash shell. Notice that the only command in this
#  file is 
#     runMCatNLO
#  which is what you need in order to obtain MC@NLO results. Other
#  commands are available: see at the bottom of this file for a 
#  list of them. In this version, they are all commented out; 
#  uncomment them if you need them.
#
# WHAT THE USER MUST DO PRIOR TO RUNNING
#  The files 
#     mcatnlo_hwdriver.f   mcatnlo_hwlhin.f
#  must be edited in order to insert the 'INCLUDE HERWIGXX.INC' command
#  relevant to the version of Herwig your are going to use. The file(s)
#     mcatnlo_hwanXXX.f
#  contain sample analysis routines, and must be edited for the same reason. 
#  Notice, however, that these analysis routines are provided here to furnish 
#  a ready-to-run package, but they are identical to standard Herwig analysis
#  routines, and should therefore be replaced with your analysis routines.
#  In this case, you will simply set the variable HWUTI (in this file) equal 
#  to the list of object files you need in your routines. 
#  Finally, the variable HERWIGVER below must be set equal to the name
#  of your preferred version of Herwig (matching the one whose common
#  blocks are included in the files above)


#
#!/bin/bash


#
#
# physical parameters
#
#
# CM energy 
ECM=1960
# renormalization scale factor [NLO]
FREN=1
# factorization scale factor [NLO]
FFACT=1
# renormalization scale factor [MC]
FRENMC=1
# factorization scale factor [MC]
FFACTMC=1
# mass of the heavy quark (top for Higgs production)
HVQMASS=173
# W mass 
WMASS=80.42
# W width
WWIDTH=2.124
# Z mass
ZMASS=91.19
# Z width
ZWIDTH=2.495
# Higgs mass
HGGMASS=125
# Higgs width: MC@NLO does not compute the SM width associated with the
# mass set in HGGMASS. The user must set the width by hand
HGGWIDTH=0
# In the computation of the Higgs cross section at the Born level:
# IBORNHGG=1 --> exact M_top dependence, IBORNHGG=2 --> M_top -> infinity
IBORNHGG=1
# If the Higgs, Z or W mass is distributed according to Breit-Wigner, 
# the mass range is
#    M0 - GAMMAX * WIDTH < M < M0 + GAMMAX * WIDTH
# M0 being the nominal mass. In the case of gamma* production, we have
#                MASSINF < Q < MASSSUP
# where Q is the virtuality of the virtual photon
GAMMAX=30
MASSINF=0
MASSSUP=0
# quark and gluon masses
UMASS=0.32
DMASS=0.32
SMASS=0.5
CMASS=1.5
BMASS=5.0
GMASS=0.75
# process number; MC@NLO process codes are negative. A positive process
# code may be used (executing runMC) to run standard Herwig
IPROC=-2850
# incoming left beam
PART1=P
# incoming right beam
PART2=PBAR
# PDF group name
PDFGROUP=MRS
# PDF set number
PDFSET=105
# Lambda_5, used in NLO computations. A negative entry returns the value
# resulting from PDF fit. 
# WARNING: negative entries may lead to errors when using PDFLIB (due to
# bugs in the library): use a positive entry when in doubt
LAMBDAFIVE=-1
# Scheme
SCHEMEOFPDF=MS
# Lambda_5, used by HERWIG. A negative entry returns the HERWIG default value
LAMBDAHERW=-1

#
#
# other input parameters
#
#
# prefix for BASES files; relevant to the integration step
FPREFIX=ww10k
# prefix for event file; relevant to the event generation step
EVPREFIX=ww10k
# prefix for the NLO and MC executables
EXEPREFIX=ww10k
# number of events; set it to 0 to skip the event generation step
NEVENTS=10000
# 0 for weights=+1/-1, 1 for weights whose sum is the total rate
WGTTYPE=1
# seed for random numbers in the generation of events. 0 is default
RNDEVSEED=0
# set BASES=ON to perform integration, =OFF to skip the integration step
BASES=ON
# set PDFLIBRARY=PDFLIB to obtain PDFs from PDFLIB, =THISLIB to obtain
# PDFs from our library of PDFs
PDFLIBRARY=THISLIB
# set HERPDF=DEFAULT to use Herwig default PDFs, HERPDF=EXTPDF to use 
# the same PDFs as used in the NLO; the setting of this parameter is 
# independent of the setting of PDFLIBRARY
HERPDF=EXTPDF
# the variable HWPATH must be set equal to the name of directory 
# which contains the version of Herwig the user wants to link
# to his code
HWPATH="/cdf/home/dwaters/mc22/herwig/"
# prepend this string to prefixes to avoid storage problems
# leave blank to store event and data files in the running directory
SCRTCH="/cdf/home/dwaters/mc22/data/"
# set the following variable equal to the list of object files that
# you need when using Herwig (for analysis purposes, for example)
HWUTI="mcatnlo_to_stdhep.o mcatnlo_hbook.o"
# set the following variable equal to the name of the version of 
# Herwig that you use (version 6.504 in this example)
HERWIGVER="herwig6504.o"
# set the following variable equal to the name of the directory where
# the PDF grid files are stored. Effective only if PDFLIBRARY=THISLIB
PDFPATH="/cdf/home/dwaters/mc22/PDFtemp/"
#
#
#
# NOW LOAD THE SCRIPTS: DO NOT REMOVE THESE LINES
thisdir=`pwd`
. $thisdir/MCatNLO.Script
#
#
#
#

#
#
# HERE, WRITE THE NAME OF THE SHELL FUNCTION THAT YOU NEED TO
# EXECUTE CHOOSING AMONG (ONLY ONE AT A TIME):
#
#    runMCatNLO  runNLO  runMC  compileNLO  compileMC
#
# THEIR MEANINGS ARE DESCRIBED IN WHAT FOLLOWS
#
#
# the following compiles and runs both the NLO and MC codes
# runMCatNLO

# the following compiles and runs the NLO only (thus, the event file 
# is written, but not read by Herwig)
# runNLO

# the following compiles and runs the MC only (thus, the event file must
# be already present, otherwise the program crashes)
# runMC

# the following compiles NLO code
# compileNLO

# the following compiles MC code
# compileMC

runMCatNLO
# runMC