Skip to content
Snippets Groups Projects
Commit b1717eb3 authored by Orel Gueta's avatar Orel Gueta Committed by Graeme Stewart
Browse files

Forbid the LHEReader from re-opening the LHEF in case it ran out of events. (Herwigpp_i-00-06-19)

parent 3cb176cd
No related branches found
No related tags found
Loading
......@@ -17,6 +17,9 @@ end_private
library Herwigpp_i *.cxx -s=components *.cxx
apply_pattern component_library
# Add the macro definitions for cmake (transparent to CMT)
apply_pattern cmake_add_command command='add_definitions(-DHWVERSION=\"${HERWIGPP_VERSION}\")'
macro Herwigpp_i_linkopts "$(Herwigpp_linkopts) "
macro_append Herwigpp_i_pp_cppflags '-DHWVERSION=\"$(Herwigpp_version)\" '
macro_append Herwigpp_i_pp_cppflags '-DHWREPODIR=\"$(Herwigpp_repodir)\" '
......
......@@ -34,6 +34,8 @@ def pdf_cmds(setname, subname):
cmds = mkpdf_cmds(setname, varname)
cmds += "set /Herwig/Particles/p+:PDF /Herwig/Partons/%s\n" % varname
cmds += "set /Herwig/Particles/pbar-:PDF /Herwig/Partons/%s\n" % varname
cmds += "set /Herwig/Partons/QCDExtractor:FirstPDF /Herwig/Partons/%s\n" % varname
cmds += "set /Herwig/Partons/QCDExtractor:SecondPDF /Herwig/Partons/%s\n" % varname
cmds += "get /Herwig/Particles/p+:PDF\n"
return cmds
......@@ -41,7 +43,13 @@ def pdf_cmds(setname, subname):
## Set a leading order PDF
def lo_pdf_cmds(lo_setname):
"""Set LO PDF as requested by arg, defaulting to current production PDF."""
return pdf_cmds(lo_setname, "LO")
cmds = pdf_cmds(lo_setname, "LO")
cmds += "\n"
cmds += "## Set PDF explicitly for MPI.\n"
cmds += "set /Herwig/Partons/MPIExtractor:FirstPDF /Herwig/Partons/AtlasPDFsetLO\n"
cmds += "set /Herwig/Partons/MPIExtractor:SecondPDF /Herwig/Partons/AtlasPDFsetLO\n"
cmds += "\n"
return cmds
## Set an NLO PDF for the matrix element, plus an LO one for MPI
......@@ -53,6 +61,8 @@ def nlo_pdf_cmds(nlo_setname, lo_setname):
cmds = mkpdf_cmds(lo_setname, "AtlasPDFsetLO")
cmds += "set /Herwig/Shower/ShowerHandler:PDFA /Herwig/Partons/AtlasPDFsetLO\n"
cmds += "set /Herwig/Shower/ShowerHandler:PDFB /Herwig/Partons/AtlasPDFsetLO\n"
cmds += "set /Herwig/Partons/MPIExtractor:FirstPDF /Herwig/Partons/AtlasPDFsetLO\n"
cmds += "set /Herwig/Partons/MPIExtractor:SecondPDF /Herwig/Partons/AtlasPDFsetLO\n"
cmds += "\n"
cmds += pdf_cmds(nlo_setname, "NLO")
cmds += "## Use 2-loop alpha_s\n"
......@@ -236,6 +246,7 @@ insert /Herwig/EventHandlers/LHEHandler:LesHouchesReaders 0 /Herwig/EventHandler
set /Herwig/EventHandlers/LHEReader:MomentumTreatment RescaleEnergy
set /Herwig/EventHandlers/LHEReader:WeightWarnings 0
set /Herwig/EventHandlers/LHEReader:AllowedToReOpen 0
set /Herwig/EventHandlers/LHEHandler:WeightOption VarNegWeight
set /Herwig/EventHandlers/LHEHandler:PartonExtractor /Herwig/Partons/QCDExtractor
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment