PROPHECY4f: a PROPer description of the Higgs dECaY into 4 fermions
===============================================================================
Prophecy4f, Version 3.0, released on Dec 01, 2019
Prophecy4f can be downloaded from
https://prophecy4f.hepforge.org/
-------------------------------------------------------------------------------
Summary of changes in PROPHECY4f 3.0:
01/12/19:
* new model SESM added
* new model THDM added
* evaluation of one-loop integrals with COLLIER library
Summary of changes in PROPHECY4f 2.0.1:
12/03/14:
* random number generator RANLUX is added to provide seedable random numbers
Summary of changes in PROPHECY4f 2.0:
08/08/11:
* production of unweighted events for leptonic final states
* optional inclusion of a fourth fermion generation
* a bug in the renormalization of the W-boson mass corrected.
This leads to an increase of the partial widths of the WW-mediated
channels by up to 0.5%.
-------------------------------------------------------------------------------
Prophecy4f has been tested under
Operating Systems: LINUX, MacOS
Compilers: GNU Fortran (GCC)
Intel Fortran Compiler
Language: Fortran 77
-------------------------------------------------------------------------------
Authors:
Ansgar Denner Wuerzburg University, Germany
(ansgar.denner@uni-wuerzburg.de)
Stefan Dittmaier Freiburg University, Germany
(stefan.dittmaier@physik.uni-freiburg.de)
Alexander Mueck RWTH Aachen University, Germany
(mueck@physik.rwth-aachen.de)
Former authors:
Axel Bredenstein
Marcus M. Weber
BRIEF DESCRIPTION:
Prophecy4f is a Monte Carlo integrator for H -> WW/ZZ -> 4fermions.
It includes:
- all four-fermion final states
- NLO QCD and electroweak corrections
- all interferences at LO and NLO
- corrections beyond NLO from heavy-Higgs effects in the SM
- alternatively an Improved Born Approximation (IBA) with leading effects of
the corrections in the SM
- production of unweighted events for leptonic final states
- models: Standard Model (SM),
4th fermion generation (SM4) (w or w/o leading 2-loop improvements),
Higgs Singlet Extension of the SM (SESM),
Two-Higgs-Doublet Model (THDM).
The present version does not (yet) include:
- multi-photon final-state radiation
- interface to parton showers
- anomalous HWW/HZZ couplings or SMEFT dimension-6 operators
- production of unweighted events for semi-leptonic and hadronic final states
(at present unweighted events only for Born and IBA for these channels)
PUBLICATIONS:
[1] A. Bredenstein, A. Denner, S. Dittmaier, M.M. Weber,
Precise predictions for the Higgs-boson decay H --> W W / Z Z --> 4 leptons,
Phys. Rev. D 74 (2006) 013004 [arXiv:hep-ph/0604011].
[2] A. Bredenstein, A. Denner, S. Dittmaier, M.M. Weber,
Precision calculations for the Higgs decays H --> Z Z / W W --> 4leptons,
Nucl. Phys. Proc. Suppl. 160 (2006) 131 [arXiv:hep-ph/0607060].
[3] A. Bredenstein, A. Denner, S. Dittmaier, M.M. Weber,
Radiative corrections to the semileptonic and hadronic Higgs-boson decays
H --> W W / Z Z --> 4 fermions,
JHEP 0702 (2007) 080 [arXiv:hep-ph/0611234].
[4] L. Altenkamp, S. Dittmaier, H. Rzehak,
Renormalization schemes for the Two-Higgs-Doublet Model and applications
to h -> WW/ZZ -> 4 fermions,
JHEP 1709 (2017) 134 [arXiv:1704.02645].
[5] L. Altenkamp, S. Dittmaier, H. Rzehak,
Precision calculations for h -> WW/ZZ -> 4 fermions in the Two-Higgs-Doublet
Model with Prophecy4f,
JHEP 1803 (2018) 110 [arXiv:1710.07598].
[6] L. Altenkamp, M. Boggia, S. Dittmaier,
Precision calculations for h -> WW/ZZ -> 4 fermions in a Singlet Extension
of the Standard Model with Prophecy4f,
JHEP 1804 (2018) 062 [arXiv:1801.07291].
[7] A. Denner, S. Dittmaier, J.-N. Lang,
Renormalization of mixing angles,
JHEP 1811 (2018) 104 [arXiv:1808.03466].
INSTALLATION & COMPILATION:
Required library: Collier (download from https://collier.hepforge.org/)
Follow the documentation of Collier for its installation.
Gunzip and untar Prophecy4f-3.0.tar.gz
(it will unpack into the directory ./Prophecy4f-3.0/)
Issue the following make command from the command line
make COLLIERDIR=path FC=compiler
where path has to be the path to the Collier library and FC is an optional
argument that allows to use the Fortran compiler "compiler" instead of the
default one. The executable "Prophecy4f" is generated.
For removing *.o files in the directory obj/ issue "make clean".
EXECUTION
For execution "Prophecy4f" needs an inputfile from the standard input.
The program can be executed using "./Prophecy4f < inputfile"
All output is written to standard output or to a file if a name is
specified in the inputfile.
INPUT
All input should be delivered via the inputfile. Its general format can be
seen from the default inputfile "defaultinput". Only those values that
differ from default have to be specified in the input file.
Do not forget the "d0" after "double precision" quantities.
An inputfile has to be specified via standard input, otherwise
Prophecy4f does not start.
In the following we present the content of the file "defaultinput" with
additional comments added. The specified values correspond to the default.
# global parameters
outputfile='' ! output is written to standard output
*******************************************************************************
* Output will be written to standard output or a given file
* if a name is provided here. The plot data will be written
* to files named plot.* in directory HISTOGRAMS/. The string 'plot'
* is replaced by the file name of the output file if provided.
* Unweighted events are written to the directory UNWEIGHTEDEVENTS/
* (*.lhe files) in the same manner. Unweighted events are also binned
* into distributions written to the directory HISTUNWEIGHTED/.
*******************************************************************************
nevents=10000000 ! nevents: number of weighted events
*******************************************************************************
* Number of weighted events, we recommend to use at least 10^7
* events for the integrated partial decay width, for histograms
* about 5*10^7 should be used.
*******************************************************************************
nunwevents=0 ! nevents: number of weighted events
*******************************************************************************
* Number of unweighted events, which are produced in the Les Houches
* event file format after the generation of "nevents" weighted events
* to find the maximal weights used for unweighting.
* Unweighted events have weight 1 or very rarely weight -1.
* For nunwevents>0 one has to use:
* qsoftcoll=2 (slicing)
* qrecomb=0 (no recombination)
* For nunwevents>0 the two parameters are set accordingly.
*******************************************************************************
model=0 ! model: 0=SM, 1=SESM, 2=THDM, 4=SM4
*******************************************************************************
* The further input that is specific to non-standard models is specified below.
*******************************************************************************
contrib=1 ! SM: contrib: 1=best 2=IBA 3=Born
! SM4, SESM and THDM: 1=best 3=Born
*******************************************************************************
* Specifies whether to calculate the partial decay width including complete
* corrections (NLO corrections as defined by qqcd (see below)
* and some higher order effects), the Improved-Born Approximation in the SM
* (see Ref. [1] for details) or the leading-order result.
*******************************************************************************
qqcd=1 ! qqcd: 0=EW 1=EW+QCD 2=QCD corrections incl.
*******************************************************************************
* Specifies whether to use only EW corrections, both EW and QCD
* corrections or only QCD corrections. For purely leptonic final
* states only EW corrections contribute.
*******************************************************************************
qsoftcoll=1 ! qsoftcoll: 1=subtraction, 2=slicing
*******************************************************************************
* Soft and collinear singularities can be treated with the subtraction
* or the slicing method. For calculations of partial decay widths
* subtraction is the preferred option, while for the production of
* unweighted events slicing has to be used.
*******************************************************************************
channel= e anti-e mu anti-mu ! final state
*******************************************************************************
* Final-state fermions (e,mu,nue,num,dq,uq,sq,cq). If more than one channel is
* specified, by including several "channel" lines in the input file,
* Prophecy4f will calculate the different channels consecutively.
* In addition one can choose the special cases:
* channel= total
* channel= leptonic
* channel= semi-leptonic
* channel= hadronic
* channel= WW
* channel= ZZ
* channel= interference
* Since final-state fermions are treated in the massless limit,
* integrated partial widths usually do not differ between different
* generations of fermions (For example, the integrated partial decay
* width for H -> e anti-e e anti-e is the same as for
* H -> mu anti-mu mu anti-mu). Symmetric final states are an exception,
* here effects of identical particles are taken into account,
* i.e. H -> e anti-e e anti-e is different from H -> e anti-e mu anti-mu.
* Moreover, in distributions (or unweighted events) fermion-mass logarithms do
* show up if no photon recombination is applied, i.e. fermions of different
* generations will in general yield different results.
* Third generation fermions cannot be used as input. However, the partial
* widths including third generation particles like bottom quarks, tau leptons,
* or tau neutrinos do not differ significantly in the massless approximation
* from those into fermions of the first and second generation,
* i.e. use e.g. H -> mu anti-mu sq anti-sq
* to calculate the H -> mu anti-mu bq anti-bq partial width.
* Top quarks in the final state are not supported.
* The old input format of Version 1.0 is also still supported.
*******************************************************************************
qrecomb=1 ! qrecomb: 0=no recomb., 1=photon recomb. invrecomb=5d0 !
recombination condition (for qrecomb=1)
*******************************************************************************
* For qrecomb=0 photons and fermions are not recombined.
* For qrecomb=1 the photon and the fermion with the smallest invariant mass
* are recombined if their invariant mass in GeV is smaller than invrecomb,
* i.e. their 4-momenta are added and attributed to the fermion. (Note that we
* cannot use a proper jet-algorithm for recombination since the lab frame of
* the Higgs decay is not specified). For inclusive partial widths recombination
* does not affect the result.
* Independent of qrecomb we always recombine the two QCD partons with the
* smallest invariant mass in events with gluon emission to form two jets in
* semileptonic decays or four jets in hadronic decays. When producing
* unweighted events for leptonic final states one has to use qrecomb=0, in
* order to create the flexibility to perform the recombination on the event
* files after production.
*******************************************************************************
qrecombcolle=0 ! 1=recombine electrons and photons inside the
! slicing cone, 0=do not recombine
*******************************************************************************
* For qrecombcolle=1 photons and electrons are recombined inside the slicing
* cone around electrons (for qsoftcoll=2). In the slicing approximation
* these recombined electron-photon pairs are strictly collinear.
* The technical slicing cone should always be smaller than the size of the
* "physical recombination cone" in which all photons are recombined.
* The slicing parameters are: deltaE = 4d-4 and deltaTh = 3d-2.
* The option qrecombcolle=1 might be useful to avoid large numbers
* of negative unweighted events for electron final states.
* For qrecombcolle=0 no electron-photon recombination is performed.
*******************************************************************************
mh = 125d0 ! Higgs boson mass in the SM or SM4
alphas = 0.118d0 ! strong coupling constant
gf = 1.1663787d-5 ! Fermi constant
mz = 91.1876d0 ! on-shell Z-boson mass
mw = 80.385d0 ! on-shell W-boson mass
gammaz = 2.4952d0 ! on-shell Z-boson width
! (only to calculate pole mass)
gammaw = 2.085d0 ! on-shell W-boson width
! (only to calculate pole mass)
me = 0.510998928d-3 ! electron mass
mmu = 105.6583715d-3 ! muon mass
mtau = 1.77682d0 ! tau mass
md = 0.100d0 ! d-quark mass
mu = 0.100d0 ! u-quark mass
ms = 0.100d0 ! s-quark mass
mc = 1.51d0 ! c-quark mass
mb = 4.92d0 ! b-quark mass
mt = 172.5d0 ! t-quark mass
*******************************************************************************
* Input parameters:
* The values of the fermion masses are needed but the results are
* practically independent of the specific values in the alpha_GF scheme
* for inclusive quantities. Only if photons and fermions are not recombined,
* logarithms of the fermion masses may appear in distributions.
*
* The values of the on-shell gauge-boson masses and widths as given in the
* input are only used to calculate the pole masses of the gauge bosons.
* For the actual evaluation of the Higgs decay, the gauge-boson widths
* are calculated from the gauge-boson pole masses and the remaining input.
* The required complex gauge-boson masses are built from these calculated
* widths and the pole masses. Upon using calculated W/Z widths, we can assure
* that the effective branching ratios of W and Z bosons add up to 1.
* More details are provided below in the "COMMENTS" section.
*******************************************************************************
randomseed=-1 ! use random numbers of Prophecy4f version 2.0
*******************************************************************************
* Other positive integer values or zero switch to the usage of RANLUX for
* random number generation, where "randomseed" is used as a seed for the
* random number generator to obtain statistically independent samples.
*******************************************************************************
*==============================================================================
* Additional input needed in the SM4:
* A 4th fermion generation of massive fermions can be optionally included
* upon setting model=4.
*==============================================================================
qsm4=1 ! 1,2 NLO,NLO+improvements
! only relevant for model=4
*******************************************************************************
* For qsm4=1 or 2 the full mass dependence of the additional closed
* fermion loops is taken into account at NLO, comprising the
* HWW/HZZ/HZA/HAA vertex corrections as well as all gauge-boson self-energies.
* The qsm4=2 option additionally takes into account the leading corrections
* ~ Gf^2*mf4^4, alpha_s*Gf*mf4^2 to the HVV vertices, which are taken from
* hep-ph/9712330 and hep-ph9602304, respectively.
*******************************************************************************
ml4 = 600d0 ! mass of l in 4th fermion generation
mn4 = 600d0 ! mass of nu in 4th fermion generation
md4 = 600d0 ! mass of d in 4th fermion generation
mu4 = 650d0 ! mass of u in 4th fermion generation
*******************************************************************************
*==============================================================================
* Additional input needed in the SESM (model=1) or THDM (model=2):
*==============================================================================
hboson=h0 ! Higgs boson in the initial state:
! h0 = light Higgs boson in SESM or THDM
! hh = heavy Higgs boson in SESM or THDM
! irrelevant for SM or SM4
*******************************************************************************
* The SESM/THDM input parameters of MSbar type are defined at the
* renormalization scale "mrenbsm1" in the renormalization scheme specified by
* "renscheme" and evolved to the renormalization scale "mrenbsm2" by solving
* the renormalization group equations numerically.
* The parameters at the scale "mrenbsm2" are used in the calculation within
* the scheme "renscheme". In particular, this can be used to evaluate the scale
* variation with respect to the BSM MSbar parameters (alpha_s is not varied).
* To allow for a simple comparison of different renormalization schemes,
* in the output the corresponding input parameters in the other ren. schemes
* are given as well, as resulting from two different parameter conversion
* techniques (see [4-7]).
*******************************************************************************
*******************************************************************************
* SESM input:
*******************************************************************************
renscheme = 5 ! renormalization scheme
! (see [6] arXiv:1801.07291 and [7] arXiv:1808.03466 )
! 0 = alpha MSbar (running l12=lambda12)-> MSbar of [6]
! 1 = alpha MSbar a la FJ (running l12) -> FJ of [6]
! 2 = on-shell scheme (running l12)
! 3 = alpha MSbar (running l1=lambda1) -> MSbar(PRTS) of [7]
! 4 = alpha MSbar a la FJ (running l1) -> MSbar(FJTS) of [7]
! 5 = on-shell scheme (running l1) -> OS of [7]
! 6 = BFM-inspired scheme -> based on Eqs.
! (3.41) and (3.64) of [7]
! 7 = BFM-inspired scheme -> BFMS of [7] based
! on Eqs. (3.41) and (3.68) of [7]
! In general, we recommend to use renscheme = 5.
mrenbsm1 = 125.1d0 ! start renormalization scale for MSbar parameters
! Here, we use mrenbsm1=mh0
mrenbsm2 = 125.1d0 ! target renormalization scale for MSbar parameters
sa = 0.29d0 ! sin(alpha)
!ta = 0.303d0 ! tan(alpha), as alternative to define alpha
mh0 = 125.1d0 ! mass of the light Higgs boson
mhh = 200d0 ! mass of the heavy Higgs boson
l12 = 0.07d0 ! coupling lambda12
*******************************************************************************
* The mixing angle alpha varies within -pi/2 < alpha < pi/2.
*******************************************************************************
*******************************************************************************
* THDM input:
*******************************************************************************
modeltype = 1 ! modeltype: variant of the model as defined in [5]
! 1 = type I: all fermions couple to Higgs doublet phi2 only
! 2 = type II: down-type fermions couple to phi1
! up-type fermions couple to phi2
! 3 = lepton specific: quarks couple to phi2, leptons to phi1
! 4 = flipped: down-type quarks couple to phi1,
! up-type quarks and charged leptons couple to phi2
renscheme = 5 ! renormalization scheme:
! (see [4] arXiv:1704.02645, [5] arXiv:1710.07598,
! [7] arXiv:1808.03466 )
! 0 = alpha/beta MSbar
! -> MSbar(alpha) of [4,5] = MSbar(PRTS) of [7]
! 1 = alpha/beta MSbar a la FJ
! -> FJ(alpha) of [4,5] = MSbar(FJTS) of [7]
! 2 = lambda3/beta MSbar
! -> MSbar(lambda_3) of [4,5]
! 3 = lambda3/beta MSbar a la FJ
! -> FJ(lambda_3) of [4,5]
! 4 = alpha/beta on-shell (nu2) -> OS2 of [7]
! 5 = alpha/beta on-shell (nu1,nu2) -> OS12 of [7]
! 6 = alpha/beta on-shell (nu1) -> OS1 of [7]
! 7 = BFM-inspired scheme
! -> based on Eqs. (3.41) and (3.74) of [7]
! 8 = BFM-inspired scheme -> BFMS of [7]
! based on Eqs. (3.41) and (3.76) of [7]
! In general, we recommend to use renscheme = 5.
mrenbsm1 = 361d0 ! start renormalization scale for MSbar parameters
! Here, we use mrenbsm1=(mh0+mhh+ma0+2mhp)/5
mrenbsm2 = 361d0 ! target renormalization scale for MSbar parameters
sa = -0.355d0! sin(alpha)
!ta = -0.380d0! tan(alpha), as alternative to define alpha
!cba = 0.1d0 ! cos(beta-alpha), as alternative to define alpha
!sgnsba = +1 ! sgn[sin(beta-alpha)], required if cba is input
tb = 2d0 ! tan(beta)
!sb = 0.894d0 ! sin(beta), as alternative to define beta
!cb = 0.447d0 ! cos(beta), as alternative to define beta
mh0 = 125d0 ! mass of the light CP-even Higgs boson
mhh = 300d0 ! mass of the heavy CP-even Higgs boson
ma0 = 460d0 ! mass of the CP-odd Higgs boson
mhp = 460d0 ! mass of the charged Higgs boson
lam5 = -1.9d0 ! coupling lambda5
*******************************************************************************
* The mixing angle beta varies in the range 0 < beta < pi/2, while the
* mixing angle alpha varies within -pi/2 < alpha < pi/2.
*******************************************************************************
OUTPUT:
All output by Prophecy4f is written to standard output or the output file
specified in the input. Information about the Collier library is always
written to standard output and the Collier output directory. In the output,
the input used for the calculation is provided. Along with the full result
and its integration error we separately show the lowest-order result and the
corresponding contributions due to EW and QCD corrections (if non-zero).
HISTOGRAMS:
A few default histograms corresponding to the distributions presented in
Refs. [1-3] are produced in the directory HISTOGRAMS/. They can be modified
in the subroutine "create_histo" in public.F. There, a subroutine called
"histogram" is called. Its first two parameters correspond to the range of
the histogram, the third parameter to the variable of the distribution, and
the number 50 refers to the number of bins.
The output format of the histograms is detailed in the corresponding output
files.
The default histograms are:
The invariant-mass distributions for final-state particles 1/2 or 3/4
for different ranges (75-85 GeV, 85-95 GeV, 50-90 GeV, 60-100 GeV). The
corresponding files read:
outputfile.inv12.5090
outputfile.inv12.7585
outputfile.inv12.60100
outputfile.inv12.8595
outputfile.inv34.5090
outputfile.inv34.7585
outputfile.inv34.60100
outputfile.inv34.8595
The final-state particles 1/2 or 3/4 and the corresponding invariant masses
refer to the first two or the last two particles as listed in the output.
outputfile.cthv2f2:
The cosine of the angle between (k(3)+k(4)) and k(2) in the Higgs rest
frame, see e.g. Fig. 12 in hep-ph/0611234 (however, the particle
numbering is different there).
outputfile.cthv2f3:
The cosine of the angle of k(3) with respect to (k(3)+k(4)) in the (k(3)+k(4))
rest frame, see e.g. Fig. 14 in hep-ph/0604011.
outputfile.phitrf2f3:
Angle between particle 2 and 3 in the transverse plane according to Fig. 15 in
hep-ph/0604011.
outputfile.cthf1f3:
The cosine of the angle between particle 1 and 3 according to Fig. 16 in
hep-ph/0604011.
outputfile.phi:
The distribution of the angle between the decay planes according to Eq. (7.9)
in hep-ph/0604011. The ordering of the momenta corresponds to the ordering of
the final-state particles in the output file. This is only used for fully
leptonic final states.
outputfile.cphihad:
The distribution of the absolute value of the cosine of the angle
between the decay planes (one plane spanned by particles 1 and 2, the
second plane spanned by particles 3 and 4), as defined in Eq. (4.2) of
hep-ph/0611234. This is only used for semi-leptonic final states.
The present version of the program does not provide histograms for
hadronic final states and only a smaller set of distributions for
semi-leptonic final states.
UNWEIGHTED EVENTS:
Unweighted events are written to the directory UNWEIGHTEDEVENTS/ in the
Les Houches event file format (*.lhe). These files contain also the
complete output in their header. As a cross check the unweighted events
are binned into distributions written to the directory HISTUNWEIGHTED/.
These histograms are equivalent to histograms obtained by binning the
events in the *.lhe files accordingly.
COMMENTS:
Electromagnetic coupling constant:
The electromagnetic coupling constant is derived from the Fermi constant.
This procedure takes into account some higher-order effects already at
tree level.
Treatment of gauge-boson resonances:
Gauge-boson resonances are treated using the complex-mass scheme.
(A. Denner, S. Dittmaier, M. Roth, L.H. Wieders,
Nucl.Phys.B724:247-294,2005, hep-ph/0505042)
As input, the program expects the on-shell W/Z masses and widths.
From these it calculates internally the real parts of the complex pole
masses. The imaginary parts of the complex pole masses, i.e. the
vector-boson widths are calculated from the pole masses and the other
input parameters as described in the next paragraph. These pole masses
are then used in propagators and the complex weak mixing angle and other
couplings.
Decay widths of vector bosons:
The decay widths of the W and Z bosons that enter the complex pole masses
are calculated from the pole masses and the other input parameters as
follows: If only LO order results are requested (i.e. for contrib=3) the
LO gauge-boson widths are used. For NLO results and the IBA (contrib=1
or contrib=2) we apply the NLO gauge-boson widths (also for the LO
sub-contribution). This ensures that the effective branching fractions
for the W- and Z-boson decays in both LO and NLO add up to one.
Note that in the Refs. [1-3] we have presented the LO results with NLO
gauge-boson widths. At LO, the difference is, of course, only a
higher-order effect.
Comment on SESM/THDM calculations:
Since the NLO W/Z widths are almost identical to the SM NLO widths,
the latter are used (with a SM Higgs mass set to the mass of the decaying
Higgs boson).
Leading-logarithmic final-state radiation:
Higher-order effects of collinear photon radiation are not yet supported
in this version.
Parton shower
An interface to a parton shower is not provided yet. In fact, naively using a
shower algorithm in combination with the NLO corrections would lead to a
double counting of real radiation.
Runtime of the program
For 10^7 weighted events the program will roughly run 10-90 minutes depending
on the final state, hardware, and compilers. The production of 10^6 unweighted
events will take about 1-2 days.
Total width
The total width for H -> 4 fermions are obtained by summing over all possible
final states:
Gamma(H->4f) = 3*Gamma(nue anti-nue num anti-num)
+ 3*Gamma(e anti-e mu anti-mu)
+ 6*Gamma(nue anti-nue mu anti-mu)
+ 6*Gamma(nue anti-e mu anti-num)
+ 3*Gamma(nue anti-nue nue anti-nue)
+ 3*Gamma(e anti-e e anti-e)
+ 3*Gamma(nue anti-e e anti-nue)
+ Gamma(uq anti-uq cq anti-cq)
+ 3*Gamma(dq anti-dq sq anti-sq)
+ 4*Gamma(uq anti-uq sq anti-sq)
+ 2*Gamma(uq anti-dq sq anti-cq)
+ 2*Gamma(uq anti-uq uq anti-uq)
+ 3*Gamma(dq anti-dq dq anti-dq)
+ 2*Gamma(uq anti-dq dq anti-uq)
+ 6*Gamma(nue anti-nue uq anti-uq)
+ 9*Gamma(nue anti-nue dq anti-dq)
+ 6*Gamma(uq anti-uq e anti-e)
+ 9*Gamma(dq anti-dq e anti-e)
+ 12*Gamma(nue anti-e dq anti-uq).
The above three blocks correspond to the leptonic, hadronic,
and semi-leptonic width, respectively, which can be obtained via
channel=leptonic, channel=hadronic, channel=semi-leptonic. The total
width can be obtained via channel=total.
The individual channels can also be combined to provide the decays
via WW, ZZ and interference as follows:
Gamma(H->W*W*->4f) =
9*Gamma(nue anti-nue mu anti-mu)
+ 12*Gamma(nue anti-e dq anti-uq)
+ 4*Gamma(uq anti-dq sq anti-cq),
Gamma(H->Z*Z*->4f) =
3*Gamma(nue anti-nue num anti-num)
+ 3*Gamma(e anti-e mu anti-mu)
+ 9*Gamma(nue anti-nue mu anti-mu)
+ 3*Gamma(nue anti-nue nue anti-nue)
+ 3*Gamma(e anti-e e anti-e)
+ 6*Gamma(nue anti-nue uq anti-uq)
+ 9*Gamma(nue anti-nue dq anti-dq)
+ 6*Gamma(uq anti-uq e anti-e)
+ 9*Gamma(dq anti-dq e anti-e)
+ Gamma(uq anti-uq cq anti-cq)
+ 3*Gamma(dq anti-dq sq anti-sq)
+ 6*Gamma(uq anti-uq sq anti-sq)
+ 2*Gamma(uq anti-uq uq anti-uq)
+ 3*Gamma(dq anti-dq dq anti-dq),
Gamma(WW/ZZ-interference) =
+ 3*Gamma(nue anti-e e anti-nue)
- 3*Gamma(nue anti-nue mu anti-mu)
- 3*Gamma(nue anti-e mu anti-num)
+ 2*Gamma(uq anti-dq dq anti-uq)
- 2*Gamma(uq anti-du ds anti-us)
- 2*Gamma(uq anti-dq ds anti-uc).
EXAMPLES:
Examples for inputfiles and the resulting output are given in the
directory example-runs/. It contains the subdirectories
example-paper/, example-channels/, example-unweighted/,
example-SESM/, and example-THDM/:
example-paper/
This directory contains input files for the decay modes
H -> ZZ -> e anti-e mu anti-mu, H -> ZZ -> e anti-e e anti-e,
H -> WW -> nue anti-e mu anti-num and H -> WW -> nue anti-e e anti-nue
in the SM with Higgs masses of 140, 170, and 200 GeV for the input
parameters of Ref.[1]. For reference the corresponding output files are
provided in "out.*" and the histograms in the directory
example-paper/HISTOGRAMS/. The results for the WW-mediated channels differ
by up to 0.5% from those given in Table 1 of Ref.[1] due to the bug in
the renormalization of the (complex) W-boson mass. The ZZ-mediated
channels give slightly different results from those in Refs. [1-3] since
the top-mass effects in the Z width calculation are treated in an
improved manner in the present version of Prophecy4f.
example-channels/
This directory contains input and corresponding output files for all final
states in the SM for a Higgs mass of 125 GeV for the default input parameter
set.
example-unweighted/
This directory contains input and corresponding output files for the
production of unweighted events for leptonic final states in the SM, a Higgs
mass of 125 GeV and the default input parameter set.
example-SESM/
This directory contains input files for the decay mode of the light CP-even
Higgs boson h0 -> WW -> num anti-mu e anti-nue of the SESM scenario BHM200 of
Refs. [6,7].
example-THDM/
This directory contains input files for the decay mode of the light CP-even
Higgs boson h0 -> WW -> num anti-mu e anti-nue of the THDM scenario Aa of
Refs. [4,5], which is identical to A1 of Ref. [7].
MPI:
PROPHECY4f supports parallel execution using MPI. To use the parallel
version, one has to compile the program using the preprocessor flag
"-Dmpiuse" and make sure that proper MPI libraries are linked. The program
will produce nevents weighted events in total and nunwevents unweighted events
per core. The parallel version of Prophecy4f has been tested using Intel's
Fortran compiler with Intel's MPI.
Large unweighted event samples:
In unweighting runs with more than 10^6 unweighted events in a single run it is
possible that fewer unweighted events are generated than requested. This is
caused by a 32bit integer overflow and can be solved by using 64bit integers
everywhere. The default size of all integers can be controlled using the
following compiler options for the Fortran compiler:
gfortran: -fdefault-integer-8
ifort: -integer-size 64
These options are the default in the supplied makefile.