Detector Simulations¶
Occasionally it is useful to interface to fast detector simulations for phenomenological studies. Two of the publically available programs AcerDet andPGS are written in FORTRAN and use the HEPEVT common block to obtain the event information.
Interfacing to these programs from is tricky as the HepMC event record must be converted to HepEVT and passed to the FORTRAN. We provide an interface to both AcerDet and PGS in the Contrib directory of Herwig, however we warn the user that due to the issues with converting between FORTRAN and C++ these interfaces are not as robust as the rest of the code.
AcerDet¶
In order to run the AcerDet interface ThePEG must be configured to use
HepMC as this is used to convert to the FORTRAN HEPEVT common block, also
for compatibilty with the FORTRAN HepMC should be configured to use units of
GeV and mm.
AcerDet can be obtained from the AcerDet web-page. By default AcerDet is compiled with g77. It will compile using gfortran versions after 4.3 which include support for the -finit-local-zero option which AcerDet requires. Herwig should be configured to use the same FORTRAN compiler as you use to compile AcerDet. In addition there is a minor bug in AcerDet’s ACDHERWIG6 routine which requires the line
N=NHEP
to be added after the declarations at the top of the routine and before
DO 180 I=1,N
The Makefile in the directory Contrib/AcerDetInterface will contain the
correct paths to most of the code, however you will have to change the
location of AcerDet, CERNLIB, and perhaps HepMC to the locations
of your version. In recent versions there may be no phtools library in
cernlib which can be reomved in the makefile and on 64-bit systems problems
with HBOOK. AcerDet does not really need HBOOK so you can comment the
HLIMIT and HOUTPU lines out in AcerDet_files.f in the interface
and all calls to HBOOK, HF1 and HOPERA in the AcerDet code.
A simple example analysis ZAnalysis which plots the reconstructed mass of the Z boson at both hadron and detector level is included.
PGS¶
The PGS program can be obtained from the PGS web-page. The most recent
version, as of writing PGS-120611, uses gfortran however the -fPIC option
needs to be added to the COMPILE-ARGS in the top level makefile and in
the CFLAGS and FFLAGS relevant for your operating system in
src/stdhep-dir/mcfio/arch_mcfio src/stdhep-dir/src/stdhep_Arch.
In order to run the PGS interface ThePEG must be configured to use HepMC
as this is used to convert to the FORTRAN HEPEVT common block, also for
compatibilty with the FORTRAN HepMC should be configured to use units of GeV
and mm.
The Makefile in the directory Contrib/PGSInterface will contain the
correct paths to most of the code, however you will have to change the
location of PGS, CERNLIB, and perhaps HepMC to the locations of
your version.
A simple example analysis ZAnalysis which plots the reconstructed mass of the Z boson at both hadron and detector level is included.
MultiWeight¶
Reading of Les Houches accord files supports reading of multiple alternative weights as described in the Les Houches Event Files section. The MultiWeight analysis provides an example on how to access the multiple weights.
The analysis is loaded like any other analysis
cd /Herwig/Generators
library MultiWeight.so
create Herwig::MultiWeight /Herwig/Analysis/MultiWeight
insert theGenerator:AnalysisHandlers 0 /Herwig/Analysis/MultiWeight
It will output any optional weights according to their ‘internal’ tags. The handler currently handles, apart from the multiple weights: FxFx tags, aMCFast weights and MadGraph5 clustering and scale information.
FxFx merging¶
The FxFx merging module consists of a Les Houches accord file reader that handles the FxFx-related tags, and a shower handler that performs the merging of the NLO events generated with MadGraph/AMC@NLO with the Herwig parton shower. An associated analysis, FxFxAnalysis is provided to allow running multiple weights within Rivet.
Note that the module has only been tested for W+jets and Z+jets events at 7 and 8 TeV. Other processes will be supported in future releases.
The file LHC-FxFx-Rivet-Wjets-15.0.in provides an example on how to load the Les Houches accord files from MadGraph/AMC@NLO.
The following options are necessary, as for MC@NLO events
set KinematicsReconstructor:ReconstructionOption General
set KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
The setting for W+jets and Z+jets events is
set /Herwig/Shower/FxFxHandler:ihrd 3
The highest-multiplicity sample should be set prior to running. For example, if the events consist of Z+0 jets@NLO, Z+1 jets@NLO and Z+2 jets@NLO then the following should be set
set FxFxHandler:njetsmax 2
Options for the merging algorithm should also be set
cd /Herwig/Shower
set FxFxHandler:VetoIsTurnedOff VetoingIsOn
# Is this the highest multiplicity ME in merging?
# merging scale
set FxFxHandler:ETClus 15.0*GeV
# Cone size used in clustering in merging.
set FxFxHandler:RClus 1.0
# Max |eta| for jets in clustering in merging.
set FxFxHandler:EtaClusMax 10
# Default 1.5 factor used to decide if a jet matches a parton
# in merging: if DR(parton,jet)<rclusfactor*rclus the parton
# and jet are said to have been matched.
set FxFxHandler:RClusFactor 1.5
The FxFxAnalysis provides an interface to Rivet just as the default RivetAnalysis
cd /Herwig/Generators
# Stuff for rivet:
create ThePEG::FxFxAnalysis /Herwig/Analysis/FxFxAnalysis FxFxAnalysis.so
# use the optional weights (WARNING: runs rivet multiple times!)
set /Herwig/Analysis/FxFxAnalysis:UseOptWeights No
insert theGenerator:AnalysisHandlers 0 /Herwig/Analysis/FxFxAnalysis
insert /Herwig/Analysis/FxFxAnalysis:Analyses 0 MC_XS
