Commit 64c1cfbd authored by Jürgen Reuter's avatar Jürgen Reuter
Browse files

Manual update.

parent da4c1c8f
......@@ -18,6 +18,7 @@
<li><a href="/trac/wiki">Wiki Page</a></li>
<li><a href="/news.html">News</a></li>
<li><a href="http://whizard.hepforge.org/tutorials">Tutorials</a></li>
<li><a href="https://whizard.hepforge.org/manual/manual014.html#sec288">Delphes Fast Simulation</a></li>
<li><a href="http://whizard.hepforge.org/talks.html">WHIZARD talks</a></li>
<li><a href="/CHANGES">ChangeLog</a></li>
</ul>
......@@ -27,6 +28,7 @@
<ul>
<li><a href="https://launchpad.net/whizard">Launchpad Support Page</a></li>
<li><a href="https://phab.hepforge.org/diffusion/WHIZARDSVN">Subversion Repository</a></li>
<li><a href="https://gitlab.tp.nt.uni-siegen.de/whizard/public">Public Git Repository</a></li>
<li><a href="https://answers.launchpad.net/whizard">Support Questions</a></li>
<li><a href="https://bugs.launchpad.net/whizard">Bug Tracker</a></li>
</ul>
......
No preview for this file type
......@@ -125,53 +125,52 @@ Hadronization</a>
</li><li><a href="manual012.html#sec270">11.7  Rescanning and recalculating events</a>
</li><li><a href="manual012.html#sec271">11.8  Negative weight events</a>
</li></ul>
</li><li><a href="manual013.html">Chapter 12  User Code Plug-Ins</a>
</li><li><a href="manual013.html">Chapter 12  Internal Data Visualization</a>
<ul>
<li><a href="manual013.html#sec273">12.1  The plug-in mechanism</a>
</li><li><a href="manual013.html#sec274">12.2  Data Types Used for Communication</a>
</li><li><a href="manual013.html#sec275">12.3  User-defined Observables and Functions</a>
</li><li><a href="manual013.html#sec283">12.4  User Code and Static Executables</a>
<li><a href="manual013.html#sec273">12.1  GAMELAN</a>
</li><li><a href="manual013.html#sec275">12.2  Histogram Display</a>
</li><li><a href="manual013.html#sec276">12.3  Plot Display</a>
</li><li><a href="manual013.html#sec277">12.4  Graphs</a>
</li><li><a href="manual013.html#sec278">12.5  Drawing options</a>
</li></ul>
</li><li><a href="manual014.html">Chapter 13  Data Visualization</a>
</li><li><a href="manual014.html">Chapter 13  Fast Detector Simulation and External Analysis</a>
<ul>
<li><a href="manual014.html#sec285">13.1  GAMELAN</a>
</li><li><a href="manual014.html#sec287">13.2  Histogram Display</a>
</li><li><a href="manual014.html#sec288">13.3  Plot Display</a>
</li><li><a href="manual014.html#sec289">13.4  Graphs</a>
</li><li><a href="manual014.html#sec290">13.5  Drawing options</a>
<li><a href="manual014.html#sec286">13.1  Interfacing ROOT</a>
</li><li><a href="manual014.html#sec287">13.2  Interfacing RIVET</a>
</li><li><a href="manual014.html#sec288">13.3  Fast Detector Simulation with DELPHES</a>
</li></ul>
</li><li><a href="manual015.html">Chapter 14  User Interfaces for WHIZARD</a>
<ul>
<li><a href="manual015.html#sec298">14.1  Command Line and <span style="font-family:monospace">SINDARIN</span> Input Files</a>
</li><li><a href="manual015.html#sec299">14.2  WHISH – The <span style="font-family:monospace">WHIZARD</span> Shell/Interactive mode</a>
</li><li><a href="manual015.html#sec300">14.3  Graphical user interface</a>
</li><li><a href="manual015.html#sec301">14.4  WHIZARD as a library</a>
<li><a href="manual015.html#sec290">14.1  Command Line and <span style="font-family:monospace">SINDARIN</span> Input Files</a>
</li><li><a href="manual015.html#sec291">14.2  WHISH – The <span style="font-family:monospace">WHIZARD</span> Shell/Interactive mode</a>
</li><li><a href="manual015.html#sec292">14.3  Graphical user interface</a>
</li><li><a href="manual015.html#sec293">14.4  WHIZARD as a library</a>
</li></ul>
</li><li><a href="manual016.html">Chapter 15  Examples</a>
<ul>
<li><a href="manual016.html#sec303">15.1  <span style="font-style:italic">Z</span> lineshape at LEP I</a>
</li><li><a href="manual016.html#sec304">15.2  <span style="font-style:italic">W</span> pairs at LEP II</a>
</li><li><a href="manual016.html#sec305">15.3  Higgs search at LEP II</a>
</li><li><a href="manual016.html#sec306">15.4  Deep Inelastic Scattering at HERA</a>
</li><li><a href="manual016.html#sec307">15.5  <span style="font-style:italic">W</span> endpoint at LHC</a>
</li><li><a href="manual016.html#sec308">15.6  SUSY Cascades at LHC</a>
</li><li><a href="manual016.html#sec309">15.7  Polarized <span style="font-style:italic">WW</span> at ILC</a>
<li><a href="manual016.html#sec295">15.1  <span style="font-style:italic">Z</span> lineshape at LEP I</a>
</li><li><a href="manual016.html#sec296">15.2  <span style="font-style:italic">W</span> pairs at LEP II</a>
</li><li><a href="manual016.html#sec297">15.3  Higgs search at LEP II</a>
</li><li><a href="manual016.html#sec298">15.4  Deep Inelastic Scattering at HERA</a>
</li><li><a href="manual016.html#sec299">15.5  <span style="font-style:italic">W</span> endpoint at LHC</a>
</li><li><a href="manual016.html#sec300">15.6  SUSY Cascades at LHC</a>
</li><li><a href="manual016.html#sec301">15.7  Polarized <span style="font-style:italic">WW</span> at ILC</a>
</li></ul>
</li><li><a href="manual017.html">Chapter 16  Technical details – Advanced Spells</a>
<ul>
<li><a href="manual017.html#sec311">16.1  Efficiency and tuning</a>
<li><a href="manual017.html#sec303">16.1  Efficiency and tuning</a>
</li></ul>
</li><li><a href="manual018.html">Chapter 17  New External Physics Models</a>
<ul>
<li><a href="manual018.html#sec313">17.1  New physics models via <span style="font-family:monospace">SARAH</span></a>
</li><li><a href="manual018.html#sec319">17.2  New physics models via <span style="font-family:monospace">FeynRules</span></a>
</li><li><a href="manual018.html#sec335">17.3  New physics models via the <span style="font-family:monospace">UFO</span> file format</a>
<li><a href="manual018.html#sec305">17.1  New physics models via <span style="font-family:monospace">SARAH</span></a>
</li><li><a href="manual018.html#sec311">17.2  New physics models via <span style="font-family:monospace">FeynRules</span></a>
</li><li><a href="manual018.html#sec327">17.3  New physics models via the <span style="font-family:monospace">UFO</span> file format</a>
</li></ul>
</li><li><a href="manual019.html">Appendix A  <span style="font-family:monospace">SINDARIN</span> Reference</a>
<ul>
<li><a href="manual019.html#sec337">A.1  Commands and Operators</a>
</li><li><a href="manual019.html#sec338">A.2  Variables</a>
</li><li><a href="manual019.html#sec341">Acknowledgements</a>
<li><a href="manual019.html#sec329">A.1  Commands and Operators</a>
</li><li><a href="manual019.html#sec330">A.2  Variables</a>
</li><li><a href="manual019.html#sec333">Acknowledgements</a>
</li></ul>
</li><li><a href="manual020.html">References</a>
</li></ul>
......
......@@ -342,75 +342,68 @@ Hadronization</a>
</li><li class="li-toc"><a href="#sec270">11.7  Rescanning and recalculating events</a>
</li><li class="li-toc"><a href="#sec271">11.8  Negative weight events</a>
</li></ul>
</li><li class="li-toc"><a href="#sec272">Chapter 12  User Code Plug-Ins</a>
</li><li class="li-toc"><a href="#sec272">Chapter 12  Internal Data Visualization</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec273">12.1  The plug-in mechanism</a>
</li><li class="li-toc"><a href="#sec274">12.2  Data Types Used for Communication</a>
</li><li class="li-toc"><a href="#sec275">12.3  User-defined Observables and Functions</a>
<a href="#sec273">12.1  GAMELAN</a>
<ul class="ftoc3"><li class="li-toc">
<a href="#sec276">12.3.1  Cut function</a>
</li><li class="li-toc"><a href="#sec277">12.3.2  Event-shape function</a>
</li><li class="li-toc"><a href="#sec278">12.3.3  Observable</a>
</li><li class="li-toc"><a href="#sec279">12.3.4  Examples</a>
<a href="#sec274">12.1.1  User-specific changes</a>
</li></ul>
</li><li class="li-toc"><a href="#sec283">12.4  User Code and Static Executables</a>
</li><li class="li-toc"><a href="#sec275">12.2  Histogram Display</a>
</li><li class="li-toc"><a href="#sec276">12.3  Plot Display</a>
</li><li class="li-toc"><a href="#sec277">12.4  Graphs</a>
</li><li class="li-toc"><a href="#sec278">12.5  Drawing options</a>
</li></ul>
</li><li class="li-toc"><a href="#sec284">Chapter 13  Data Visualization</a>
</li><li class="li-toc"><a href="#sec285">Chapter 13  Fast Detector Simulation and External Analysis</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec285">13.1  GAMELAN</a>
<ul class="ftoc3"><li class="li-toc">
<a href="#sec286">13.1.1  User-specific changes</a>
</li></ul>
</li><li class="li-toc"><a href="#sec287">13.2  Histogram Display</a>
</li><li class="li-toc"><a href="#sec288">13.3  Plot Display</a>
</li><li class="li-toc"><a href="#sec289">13.4  Graphs</a>
</li><li class="li-toc"><a href="#sec290">13.5  Drawing options</a>
<a href="#sec286">13.1  Interfacing ROOT</a>
</li><li class="li-toc"><a href="#sec287">13.2  Interfacing RIVET</a>
</li><li class="li-toc"><a href="#sec288">13.3  Fast Detector Simulation with DELPHES</a>
</li></ul>
</li><li class="li-toc"><a href="#sec297">Chapter 14  User Interfaces for WHIZARD</a>
</li><li class="li-toc"><a href="#sec289">Chapter 14  User Interfaces for WHIZARD</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec298">14.1  Command Line and <span style="font-family:monospace">SINDARIN</span> Input Files</a>
</li><li class="li-toc"><a href="#sec299">14.2  WHISH – The <span style="font-family:monospace">WHIZARD</span> Shell/Interactive mode</a>
</li><li class="li-toc"><a href="#sec300">14.3  Graphical user interface</a>
</li><li class="li-toc"><a href="#sec301">14.4  WHIZARD as a library</a>
<a href="#sec290">14.1  Command Line and <span style="font-family:monospace">SINDARIN</span> Input Files</a>
</li><li class="li-toc"><a href="#sec291">14.2  WHISH – The <span style="font-family:monospace">WHIZARD</span> Shell/Interactive mode</a>
</li><li class="li-toc"><a href="#sec292">14.3  Graphical user interface</a>
</li><li class="li-toc"><a href="#sec293">14.4  WHIZARD as a library</a>
</li></ul>
</li><li class="li-toc"><a href="#sec302">Chapter 15  Examples</a>
</li><li class="li-toc"><a href="#sec294">Chapter 15  Examples</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec303">15.1  <span style="font-style:italic">Z</span> lineshape at LEP I</a>
</li><li class="li-toc"><a href="#sec304">15.2  <span style="font-style:italic">W</span> pairs at LEP II</a>
</li><li class="li-toc"><a href="#sec305">15.3  Higgs search at LEP II</a>
</li><li class="li-toc"><a href="#sec306">15.4  Deep Inelastic Scattering at HERA</a>
</li><li class="li-toc"><a href="#sec307">15.5  <span style="font-style:italic">W</span> endpoint at LHC</a>
</li><li class="li-toc"><a href="#sec308">15.6  SUSY Cascades at LHC</a>
</li><li class="li-toc"><a href="#sec309">15.7  Polarized <span style="font-style:italic">WW</span> at ILC</a>
<a href="#sec295">15.1  <span style="font-style:italic">Z</span> lineshape at LEP I</a>
</li><li class="li-toc"><a href="#sec296">15.2  <span style="font-style:italic">W</span> pairs at LEP II</a>
</li><li class="li-toc"><a href="#sec297">15.3  Higgs search at LEP II</a>
</li><li class="li-toc"><a href="#sec298">15.4  Deep Inelastic Scattering at HERA</a>
</li><li class="li-toc"><a href="#sec299">15.5  <span style="font-style:italic">W</span> endpoint at LHC</a>
</li><li class="li-toc"><a href="#sec300">15.6  SUSY Cascades at LHC</a>
</li><li class="li-toc"><a href="#sec301">15.7  Polarized <span style="font-style:italic">WW</span> at ILC</a>
</li></ul>
</li><li class="li-toc"><a href="#sec310">Chapter 16  Technical details – Advanced Spells</a>
</li><li class="li-toc"><a href="#sec302">Chapter 16  Technical details – Advanced Spells</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec311">16.1  Efficiency and tuning</a>
<a href="#sec303">16.1  Efficiency and tuning</a>
</li></ul>
</li><li class="li-toc"><a href="#sec312">Chapter 17  New External Physics Models</a>
</li><li class="li-toc"><a href="#sec304">Chapter 17  New External Physics Models</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec313">17.1  New physics models via <span style="font-family:monospace">SARAH</span></a>
<a href="#sec305">17.1  New physics models via <span style="font-family:monospace">SARAH</span></a>
<ul class="ftoc3"><li class="li-toc">
<a href="#sec314">17.1.1  <span style="font-family:monospace">WHIZARD</span>/<span style="font-family:monospace">O’Mega</span> model files from <span style="font-family:monospace">SARAH</span></a>
</li><li class="li-toc"><a href="#sec317">17.1.2  Linking <span style="font-family:monospace">SPheno</span> and <span style="font-family:monospace">WHIZARD</span></a>
</li><li class="li-toc"><a href="#sec318">17.1.3  BSM Toolbox</a>
<a href="#sec306">17.1.1  <span style="font-family:monospace">WHIZARD</span>/<span style="font-family:monospace">O’Mega</span> model files from <span style="font-family:monospace">SARAH</span></a>
</li><li class="li-toc"><a href="#sec309">17.1.2  Linking <span style="font-family:monospace">SPheno</span> and <span style="font-family:monospace">WHIZARD</span></a>
</li><li class="li-toc"><a href="#sec310">17.1.3  BSM Toolbox</a>
</li></ul>
</li><li class="li-toc"><a href="#sec319">17.2  New physics models via <span style="font-family:monospace">FeynRules</span></a>
</li><li class="li-toc"><a href="#sec311">17.2  New physics models via <span style="font-family:monospace">FeynRules</span></a>
<ul class="ftoc3"><li class="li-toc">
<a href="#sec320">17.2.1  Installation and Usage of the <span style="font-family:monospace">WHIZARD</span>-<span style="font-family:monospace">FeynRules</span> interface</a>
</li><li class="li-toc"><a href="#sec326">17.2.2  Options of the <span style="font-family:monospace">WHIZARD</span>-<span style="font-family:monospace">FeynRules</span> interface</a>
</li><li class="li-toc"><a href="#sec327">17.2.3  Validation of the interface</a>
</li><li class="li-toc"><a href="#sec328">17.2.4  Examples for the <span style="font-family:monospace">WHIZARD</span>-/<span style="font-family:monospace">FeynRules</span> interface</a>
<a href="#sec312">17.2.1  Installation and Usage of the <span style="font-family:monospace">WHIZARD</span>-<span style="font-family:monospace">FeynRules</span> interface</a>
</li><li class="li-toc"><a href="#sec318">17.2.2  Options of the <span style="font-family:monospace">WHIZARD</span>-<span style="font-family:monospace">FeynRules</span> interface</a>
</li><li class="li-toc"><a href="#sec319">17.2.3  Validation of the interface</a>
</li><li class="li-toc"><a href="#sec320">17.2.4  Examples for the <span style="font-family:monospace">WHIZARD</span>-/<span style="font-family:monospace">FeynRules</span> interface</a>
</li></ul>
</li><li class="li-toc"><a href="#sec335">17.3  New physics models via the <span style="font-family:monospace">UFO</span> file format</a>
</li><li class="li-toc"><a href="#sec327">17.3  New physics models via the <span style="font-family:monospace">UFO</span> file format</a>
</li></ul>
</li><li class="li-toc"><a href="#sec336">Appendix A  <span style="font-family:monospace">SINDARIN</span> Reference</a>
</li><li class="li-toc"><a href="#sec328">Appendix A  <span style="font-family:monospace">SINDARIN</span> Reference</a>
<ul class="ftoc2"><li class="li-toc">
<a href="#sec337">A.1  Commands and Operators</a>
</li><li class="li-toc"><a href="#sec338">A.2  Variables</a>
<a href="#sec329">A.1  Commands and Operators</a>
</li><li class="li-toc"><a href="#sec330">A.2  Variables</a>
<ul class="ftoc3"><li class="li-toc">
<a href="#sec339">A.2.1  Rebuild Variables</a>
</li><li class="li-toc"><a href="#sec340">A.2.2  Standard Variables</a>
<a href="#sec331">A.2.1  Rebuild Variables</a>
</li><li class="li-toc"><a href="#sec332">A.2.2  Standard Variables</a>
</li></ul>
</li></ul>
</li></ul>
......@@ -495,10 +488,7 @@ for their extension. After that, the next chapter discusses parton
showering, matching and hadronization as well as options for event
normalizations and supported event formats. Also weighted event
generation is explained along the lines with options for negative
weights. Then, in Chap. <a href="#chap%3Auser">12</a>, options for user to plug-in
self-written code into the <span style="font-family:monospace">WHIZARD</span> framework are detailed, e.g. for
observables, selections and cut functions, or for spectra and
structure functions. Also, static executables are discussed.</p><p>Chap. <a href="#chap%3Avisualization">13</a> is a stand-alone documentation of
weights.</p><p>Chap. <a href="#chap%3Avisualization">12</a> is a stand-alone documentation of
GAMELAN, the interal graphics support for the visualization of data
and analysis. The next chapter, Chap. <a href="#chap%3Auserint">14</a> details user
interfaces: how to use more options of the <span style="font-family:monospace">WHIZARD</span> command on the
......@@ -823,9 +813,8 @@ especially the the early version of the <span style="font-family:monospace">gfor
some bugs. <span style="font-family:monospace">gfortran</span> 6.5.0 has a severe regression and cannot
be used.
</li><li class="li-itemize"><span style="font-family:monospace">nagfor</span> (NAG). You will need version 6.2 or higher.
</li><li class="li-itemize"><span style="font-family:monospace">ifort</span> (Intel). You will need version 17.0.4 or
higher. Version 19.0.0 unfortunately has a severe regression and
cannot be used.
</li><li class="li-itemize"><span style="font-family:monospace">ifort</span> (Intel). You will need version 19.0.2 or
higher
</li></ul>
<!--TOC subsection id="sec14" 2.2.5  LHAPDF-->
<h3 id="sec14" class="subsection">2.2.5  LHAPDF</h3><!--SEC END --><p>
......@@ -984,7 +973,7 @@ respectively. After configuration, do <span style="font-family:monospace">make</
<span style="font-family:monospace">make check</span> (which might sometimes fail for non-standard values
of momentum and length), and finally <span style="font-family:monospace">make install</span>.</p><p>The latest version of <span style="font-family:monospace">HepMC</span> (2.6.10) as well as the new relase
series use <span style="font-family:monospace">cmake</span> for their build process. For more
information, confer the <span style="font-family:monospace">HepMC</span> webpage.</p><p>A <span style="font-family:monospace">WHIZARD</span> configuration for HepMC looks like this:
information, confer the <span style="font-family:monospace">HepMC</span> webpage.</p><p>A <span style="font-family:monospace">WHIZARD</span> configuration for <span style="font-family:monospace">HepMC</span> looks like this:
<span style="font-size:small">
</span></p><pre class="verbatim"><span style="font-size:small"> configure: --------------------------------------------------------------
configure: --- HepMC ---
......@@ -1001,7 +990,31 @@ would not find it, set the environment variable <span style="font-family:monospa
the correct installation path when configuring <span style="font-family:monospace">WHIZARD</span>. Furthermore,
the environment variable <span style="font-family:monospace">CXXFLAGS</span> allows you to set specific
<span style="font-family:monospace">C/C++</span> preprocessor flags, e.g. non-standard include paths for
header files.</p>
header files.</p><p>A typical configuration of <span style="font-family:monospace">HepMC3</span> will look like this:
<span style="font-size:small">
</span></p><pre class="verbatim"><span style="font-size:small">configure: --------------------------------------------------------------
configure: --- ROOT ---
configure:
checking for root-config... /usr/local/bin/root-config
checking for root... /usr/local/bin/root
checking for rootcint... /usr/local/bin/rootcint
checking for dlopen in -ldl... (cached) yes
configure: --------------------------------------------------------------
configure: --- HepMC ---
configure:
checking for HepMC3-config... /usr/local/bin/HepMC3-config
checking if HepMC3 is built with ROOT interface... yes
checking if HepMC3 is functional... yes
checking for HepMC3... yes
checking the HepMC3 version... 3.02.01
configure: --------------------------------------------------------------
</span></pre><p><span style="font-size:small">
</span>
As can be seen, <span style="font-family:monospace">WHIZARD</span> will check for the <span style="font-family:monospace">ROOT</span> environment as
well as whether <span style="font-family:monospace">HepMC3</span> has been built with support for the
<span style="font-family:monospace">ROOT</span> and <span style="font-family:monospace">RootTree</span> writer classes. This is an easy option to
use <span style="font-family:monospace">WHIZARD</span> to write out <span style="font-family:monospace">ROOT</span> events. For more information see
Sec. <a href="#sec%3Aroot">13.1</a>. </p>
<!--TOC subsection id="sec17" 2.2.8  PYTHIA8-->
<h3 id="sec17" class="subsection">2.2.8  PYTHIA8</h3><!--SEC END --><p>
<a id="sec:pythia8"></a></p><p><em>NOTE: This is at the moment not yet supported, but merely a stub
......@@ -2819,7 +2832,7 @@ useful for pretty-printing single histograms or plots. Example:
} = hist1 &amp; hist2
</pre> </td></tr>
</table><p>
See Sec. <a href="#sec%3Agraphs">13.4</a>.</p>
See Sec. <a href="#sec%3Agraphs">12.4</a>.</p>
<!--TOC subsubsection id="sec88" write_analysis-->
<h4 id="sec88" class="subsubsection">write_analysis</h4><!--SEC END --><blockquote class="quote"><div class="flushleft"><span style="font-family:monospace">
write_analysis (</span>⟨<span style="font-family:monospace"><span style="font-style:italic">analysis-objects</span></span>⟩<span style="font-family:monospace">)
......@@ -5164,9 +5177,198 @@ the <span style="font-style:italic">e</span> → <span style="font-style:italic"
photon-induced backgrounds in lepton collider physics. The original
concept is that of the Weizsäcker-Williams
approximation [<a href="#vonWeizsacker%3A1934sx">21</a>, <a href="#Williams%3A1934ad">22</a>, <a href="#Budnev%3A1974de">23</a>]. This
is a single-beam structure
function that can be applied to both beams, or also to one beam
only. Examples are:
is a single-beam structure function that can be applied to both beams,
or also to one beam only. Usually, there are some simplifications
being made in the derivation. The formula which is implemented here
and seems to be the best for the QCD background for low-<span style="font-style:italic">p</span><sub><span style="font-style:italic">T</span></sub> hadrons,
corresponds to Eq. (6.17) of Ref. [<a href="#Budnev%3A1974de">23</a>]. As this
reference already found, this leads to an "overshooting" of accuracy,
and especially in the high-<span style="font-style:italic">x</span> (high-energy) region to wrong
results. This formula corresponds to
</p><table class="display dcenter"><tr style="vertical-align:middle"><td class="dcell">
<a id="eq:budnev_617"></a>
<span style="font-style:italic">f</span>(<span style="font-style:italic">x</span>) = </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center">α</td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center">π</td></tr>
</table></td><td class="dcell"> </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center">1</td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span></td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎡<br>
⎢<br>
⎢<br>
⎣</td><td class="dcell">⎛<br>
⎜<br>
⎜<br>
⎝</td><td class="dcell"><span style="text-decoration:overline">x</span> +
</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center">2</td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎞<br>
⎟<br>
⎟<br>
⎠</td><td class="dcell">log
</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup></td></tr>
</table></td><td class="dcell">
− </td><td class="dcell">⎛<br>
⎜<br>
⎜<br>
⎝</td><td class="dcell">1 − </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center">2</td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎞<br>
⎟<br>
⎟<br>
⎠</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:left">2</td></tr>
<tr><td class="dcell" style="text-align:left"><br>
<br>
<br>
</td></tr>
<tr><td class="dcell" style="text-align:left">&nbsp;</td></tr>
</table></td><td class="dcell">
log</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span><sup>2</sup> + <span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup>/<span style="font-style:italic">E</span><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span><sup>2</sup> +
<span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup>/<span style="font-style:italic">E</span><sup>2</sup></td></tr>
</table></td><td class="dcell">
− </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">m</span><sub><span style="font-style:italic">e</span></sub><sup>2</sup> <span style="font-style:italic">x</span><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup></td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎛<br>
⎜<br>
⎜<br>
⎝</td><td class="dcell">1 −
</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup></td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎞<br>
⎟<br>
⎟<br>
⎠</td><td class="dcell">⎤<br>
⎥<br>
⎥<br>
⎦</td><td class="dcell">     .
    (1)</td></tr>
</table><p>
Here, <span style="font-style:italic">x</span> is the ratio of the photon energy (called frequency ω
in [<a href="#Budnev%3A1974de">23</a>] over the original electron (or positron) beam
energy <span style="font-style:italic">E</span>. The energy of the electron (or positron) after the
splitting is given by <span style="text-decoration:overline">x</span> = 1−<span style="font-style:italic">x</span>.</p><p>The simplified version is the one that corresponds to many
publications about the EPA during SLC and LEP times, and corresponds
to the <span style="font-style:italic">q</span><sup>2</sup> integration of Eq. (6.16e) in [<a href="#Budnev%3A1974de">23</a>], where
<span style="font-style:italic">q</span><sup>2</sup> is the virtuality or momentum transfer of the photon in the EPA:
</p><table class="display dcenter"><tr style="vertical-align:middle"><td class="dcell">
<a id="eq:budnev_616e"></a>
<span style="font-style:italic">f</span>(<span style="font-style:italic">x</span>) = </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center">α</td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center">π</td></tr>
</table></td><td class="dcell"> </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center">1</td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span></td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎡<br>
⎢<br>
⎢<br>
⎣</td><td class="dcell">⎛<br>
⎜<br>
⎜<br>
⎝</td><td class="dcell"><span style="text-decoration:overline">x</span> +
</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center">2</td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎞<br>
⎟<br>
⎟<br>
⎠</td><td class="dcell">log
</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup></td></tr>
</table></td><td class="dcell">
− </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">m</span><sub><span style="font-style:italic">e</span></sub><sup>2</sup> <span style="font-style:italic">x</span><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup></td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎛<br>
⎜<br>
⎜<br>
⎝</td><td class="dcell">1 −
</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup></td></tr>
</table></td><td class="dcell"> </td><td class="dcell">⎞<br>
⎟<br>
⎟<br>
⎠</td><td class="dcell">⎤<br>
⎥<br>
⎥<br>
⎦</td><td class="dcell">     .
    (2)</td></tr>
</table><p>
While Eq. (<a href="#eq%3Abudnev_617">1</a>) is supposed to be the better choice
for simulating hadronic background like low-<span style="font-style:italic">p</span><sub><span style="font-style:italic">T</span></sub> hadrons and should
be applied for the low-<span style="font-style:italic">x</span> region of the EPA,
Eq. (<a href="#eq%3Abudnev_616e">2</a>) seems better suited for high-<span style="font-style:italic">x</span>
simulations like the photoproduction of BSM resonances etc.
Note that the first term in Eqs. (<a href="#eq%3Abudnev_617">1</a>) and
(<a href="#eq%3Abudnev_616e">2</a>) is the standard Altarelli-Parisi QED splitting
function of electron, <span style="font-style:italic">P</span><sub><span style="font-style:italic">e</span>→ <span style="font-style:italic">e</span>γ</sub>(<span style="font-style:italic">x</span>) ∝ 1 + (1−<span style="font-style:italic">x</span>)<sup>2</sup>, while
the last term in both equations is the default power correction.</p><p>The two parameters <span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup> and <span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup> are the
integration boundaries of the photon virtuality integration. Usually,
they are given by the kinematic limits:
</p><table class="display dcenter"><tr style="vertical-align:middle"><td class="dcell">
<span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup> = </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">m</span><sub><span style="font-style:italic">e</span></sub><sup>2</sup> <span style="font-style:italic">x</span><sup>2</sup></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="text-decoration:overline">x</span></td></tr>
</table></td><td class="dcell">         
<span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup> = 4 <span style="font-style:italic">E</span><sup>2</sup> <span style="text-decoration:overline">x</span> = <span style="font-style:italic">s</span> <span style="text-decoration:overline">x</span>      .
    (3)</td></tr>
</table><p>
For low-<span style="font-style:italic">p</span><sub><span style="font-style:italic">T</span></sub> hadron simulations, it is not a good idea to take the
kinematic limit as an upper limit, but one should cut the simulation
off at a hadronic scale like e.g. a multiple of the ρ mass.</p><p>The user can switch between the two different options using the setting
</p><blockquote class="quote">
<span style="font-size:small">
</span><pre class="verbatim"><span style="font-size:small">$epa_mode = "default"
</span></pre><span style="font-size:small">
</span>
</blockquote><p>
or
</p><blockquote class="quote">
<span style="font-size:small">
</span><pre class="verbatim"><span style="font-size:small">$epa_mode = "Budnev_617"
</span></pre><span style="font-size:small">
</span>
</blockquote><p>
for Eq. (<a href="#eq%3Abudnev_617">1</a>), while Eq. (<a href="#eq%3Abudnev_616e">2</a>) can be
chosen with
</p><blockquote class="quote">
<span style="font-size:small">
</span><pre class="verbatim"><span style="font-size:small">$epa_mode = "Budnev_616e"
</span></pre><span style="font-size:small">
</span>
</blockquote><p>
Note that a thorough study for high-energy <span style="font-style:italic">e</span><sup>+</sup><span style="font-style:italic">e</span><sup>−</sup> colliders
regarding the suitability of different EPA options is still lacking.</p><p>For testing purposes also three more variants or simplifications of
Eq. (<a href="#eq%3Abudnev_616e">2</a>) are implemented: the first, steered by
<span style="font-family:monospace">$epa_mode = log_power</span> uses simply <span style="font-style:italic">Q</span><sub>max</sub><sup>2</sup> =
<span style="font-style:italic">s</span>. This is also the case for the two other method. But the switch
<span style="font-family:monospace">$epa_mode = log_simple</span> uses just <span style="font-family:monospace">epa_mass</span> (cf. below)
as <span style="font-style:italic">Q</span><sub>min</sub><sup>2</sup>. The final simplification is to drop the power
correction, which can be chosen with <span style="font-family:monospace">$epa_mode = log</span>. This
corresponds to the simple formula:
</p><table class="display dcenter"><tr style="vertical-align:middle"><td class="dcell">
<span style="font-style:italic">f</span>(<span style="font-style:italic">x</span>) = </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center">α</td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center">2π</td></tr>
</table></td><td class="dcell"> </td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center">1</td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">x</span></td></tr>
</table></td><td class="dcell">   log</td><td class="dcell"><table class="display"><tr><td class="dcell" style="text-align:center"><span style="font-style:italic">s</span></td></tr>
<tr><td class="hbar"></td></tr>
<tr><td class="dcell" style="text-align:center"><span style="font-style:italic">m</span><sup>2</sup></td></tr>
</table></td><td class="dcell">
     . 
    (4)</td></tr>
</table><p>Examples for the application of the EPA in <span style="font-family:monospace">WHIZARD</span> are:
</p><blockquote class="quote">
<span style="font-size:small">
</span><pre class="verbatim"><span style="font-size:small">beams = e1, E1 =&gt; epa
......@@ -7030,7 +7232,7 @@ written.<p>The command is able to print more than one dataset, following the syn
The argument in brackets may also be empty or absent; in this case, all
currently existing datasets are printed.</p><p>The default data format is suitable for compiling analysis data by <span style="font-family:monospace">WHIZARD</span>’s
built-in <span style="font-family:sans-serif">gamelan</span> graphics driver (see below and particularly
Chap. <a href="#chap%3Avisualization">13</a>). Data are written in
Chap. <a href="#chap%3Avisualization">12</a>). Data are written in
blank-separated fixed columns, headlines and comments are initiated by the
<code>#</code> sign, and each data set is terminated by a blank line. However,
external programs often require special formatting.</p><p>The internal graphics driver <span style="font-family:sans-serif">gamelan</span> of <span style="font-family:monospace">WHIZARD</span> is initiated by
......@@ -7038,7 +7240,7 @@ the <span style="font-family:monospace">compile_analysis</span> command. Its syn
contains the <span style="font-family:monospace">write_analysis</span> if that has not been separately
called (which is unnecessary). For more details about the <span style="font-family:sans-serif">gamelan</span>
graphics driver and data visualization within <span style="font-family:monospace">WHIZARD</span>, confer
Chap. <a href="#chap%3Avisualization">13</a>.</p></li><li class="li-enumerate">Custom format. Not yet (re-)implemented in a general form.
Chap. <a href="#chap%3Avisualization">12</a>.</p></li><li class="li-enumerate">Custom format. Not yet (re-)implemented in a general form.
</li></ol>
<!--TOC section id="sec187" 5.10  Custom Input/Output-->
<h2 id="sec187" class="section">5.10  Custom Input/Output</h2><!--SEC END --><p>
......@@ -7423,10 +7625,10 @@ dΦ<sub>rad</sub> Δ<sub>s</sub>(<span style="font-style:italic">k</span><sub>T
⎟<br>
⎟<br>
⎠</td><td class="dcell">   where </td></tr>
</table></td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (1)</td></tr>
</table></td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (5)</td></tr>
<tr><td style="text-align:right;white-space:nowrap" >
<span style="text-decoration:overline">B</span><sub>s</sub></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">B</span> + <span style="color:red"><span style="font-style:italic">V</span></span> + dΦ<sub>rad</sub> 
<span style="color:red"><span style="font-style:italic">R</span></span><sub>s</sub>    and </td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (2)</td></tr>
<span style="color:red"><span style="font-style:italic">R</span></span><sub>s</sub>    and </td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (6)</td></tr>
<tr><td style="text-align:right;white-space:nowrap" >
Δ<sub>s</sub>(<span style="font-style:italic">p</span><sub><span style="font-style:italic">T</span></sub>)</td><td style="text-align:left;white-space:nowrap" ><table class="display"><tr style="vertical-align:middle"><td class="dcell">= exp</td><td class="dcell">⎡<br>
⎢<br>
......@@ -7443,11 +7645,11 @@ dΦ<sub>rad</sub> Δ<sub>s</sub>(<span style="font-style:italic">k</span><sub>T
⎥<br>
⎦</td><td class="dcell"> .
</td></tr>
</table></td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (3)</td></tr>
</table></td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (7)</td></tr>
</table></td></tr>
</table><p>
The subscript s refers to the singular part of the real component, cf.
to the next subsection. Eq. (<a href="#eq%3Apowheg">1</a>) produces either no or one
to the next subsection. Eq. (<a href="#eq%3Apowheg">5</a>) produces either no or one
additional emission. These events can then either be analyzed directly
or passed on to the parton shower<sup><a id="text15" href="#note15">6</a></sup>
for the full simulation. You activate this with
......@@ -7456,7 +7658,7 @@ for the full simulation. You activate this with
?powheg_matching = true
</pre> </td></tr>
</table><p>
The <span style="font-style:italic">p</span><sub><span style="font-style:italic">T</span></sub><sup>min</sup> of Eq. (<a href="#eq%3Apowheg">1</a>) can be set with
The <span style="font-style:italic">p</span><sub><span style="font-style:italic">T</span></sub><sup>min</sup> of Eq. (<a href="#eq%3Apowheg">5</a>) can be set with
<span style="font-family:monospace">powheg_pt_min</span>. It sets the minimal scale for the <span style="font-family:monospace">POWHEG</span>
evolution and should be of order 1 GeV and set accordingly in the
interfaced shower. The maximal scale is currently given by <span style="font-family:monospace">sqrts</span>
......@@ -7480,11 +7682,11 @@ will also only use the singular parts.</p><p>The current implementation uses the
</p><table class="display dcenter"><tr style="vertical-align:middle"><td class="dcell">
     
 
</td><td class="dcell"><table style="border-spacing:6px;border-collapse:separate;" class="cellpading0"><tr><td style="text-align:right;white-space:nowrap" >  <span style="font-style:italic">R</span></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">R</span><sub>fin</sub> + <span style="font-style:italic">R</span><sub>sing</sub>  ,</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (4)</td></tr>
</td><td class="dcell"><table style="border-spacing:6px;border-collapse:separate;" class="cellpading0"><tr><td style="text-align:right;white-space:nowrap" >  <span style="font-style:italic">R</span></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">R</span><sub>fin</sub> + <span style="font-style:italic">R</span><sub>sing</sub>  ,</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (8)</td></tr>
<tr><td style="text-align:right;white-space:nowrap" >
<span style="font-style:italic">R</span><sub>sing</sub></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">R</span> <span style="font-style:italic">F</span>(Φ<sub><span style="font-style:italic">n</span>+1</sub>)  ,</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (5)</td></tr>
<span style="font-style:italic">R</span><sub>sing</sub></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">R</span> <span style="font-style:italic">F</span>(Φ<sub><span style="font-style:italic">n</span>+1</sub>)  ,</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (9)</td></tr>
<tr><td style="text-align:right;white-space:nowrap" >
<span style="font-style:italic">R</span><sub>fin</sub></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">R</span> (1−<span style="font-style:italic">F</span>(Φ<sub><span style="font-style:italic">n</span>+1</sub>))  ,</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (6)</td></tr>
<span style="font-style:italic">R</span><sub>fin</sub></td><td style="text-align:left;white-space:nowrap" >= <span style="font-style:italic">R</span> (1−<span style="font-style:italic">F</span>(Φ<sub><span style="font-style:italic">n</span>+1</sub>))  ,</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (10)</td></tr>
<tr><td style="text-align:right;white-space:nowrap" >
<span style="font-style:italic">F</span>(Φ<sub><span style="font-style:italic">n</span>+1</sub>)</td><td style="text-align:left;white-space:nowrap" ><table class="display"><tr style="vertical-align:middle"><td class="dcell">=
</td><td class="dcell"><table class="display"><tr style="vertical-align:middle"><td class="dcell">⎧<br>
......@@ -7501,7 +7703,7 @@ will also only use the singular parts.</p><p>The current implementation uses the
</table></td></tr>
</table></td><td class="dcell"> .
</td></tr>
</table></td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (7)</td></tr>
</table></td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (11)</td></tr>
</table></td></tr>
</table><p>
Thus, a point is singular (<span style="font-style:italic">F</span>=1), if any of the FKS tuples
......@@ -7522,7 +7724,7 @@ This simplifies in massless case to
</table></td></tr>
</table></td><td class="dcell"> .
</td></tr>
</table></td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (8)</td></tr>
</table></td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >&nbsp;</td><td style="text-align:left;white-space:nowrap" >&nbsp;</td><td style="text-align:right;white-space:nowrap" >    (12)</td></tr>
</table></td></tr>
</table>
<!--BEGIN NOTES chapter-->
......@@ -10146,369 +10348,11 @@ physically irrelevant, but in the output they are representable and
distinguishable from zero.</div></dd><dt class="dt-thefootnotes"><a id="note18" href="#text18">2</a></dt><dd class="dd-thefootnotes"><div class="footnotetext">The FSR is always simulated with
angular ordering enabled.</div></dd></dl>
<!--END NOTES-->
<!--TOC chapter id="sec272" Chapter 12  User Code Plug-Ins-->
<h1 id="sec272" class="chapter">Chapter 12  User Code Plug-Ins</h1><!--SEC END --><p>
<a id="chap:user"></a></p><p><span style="color:red">
Note that the user-code plug-in mechanism has been currently (for
version 2.2.0) disabled, as the huge refactoring of the code between
versions 2.1.X and 2.2.X has completely changed many of the
interfaces. We plan to bring the interface for user code for spectra,
structure functions and event shapes, cuts and observables back online
as soon as possible, at latest for version 2.4.0.
</span></p><p><br>
<br>
<br>
<br>
<br>
<br>
</p>
<!--TOC section id="sec273" 12.1  The plug-in mechanism-->
<h2 id="sec273" class="section">12.1  The plug-in mechanism</h2><!--SEC END --><p>The capabilities of <span style="font-family:monospace">WHIZARD</span> and its <span style="font-family:monospace">SINDARIN</span> command
language are not always sufficient to adapt to all users’ needs. To
make the program more versatile, there are several spots in the
workflow where the user may plug in his/her own code, to enhance or
modify the default behavior.</p><p>User code can be injected, without touching <span style="font-family:monospace">WHIZARD</span>’s source code, in
the following places:
</p><ul class="itemize"><li class="li-itemize">
Cuts, weights, analysis, etc.:
<ul class="itemize"><li class="li-itemize">
Cut functions that operate on a whole subevent.
</li><li class="li-itemize">Observable (e.g., event shapes) calculated from a whole subevent.
</li><li class="li-itemize">Observable calculated for a particle or particle pair.
</li></ul>
</li><li class="li-itemize">Spectra and structure functions.
</li></ul><p>
Additional plug-in locations may be added in the future.</p><p>User code is loaded dynamically by <span style="font-family:monospace">WHIZARD</span>. There are two
possibilities:
</p><ol class="enumerate" type=1><li class="li-enumerate">
The user codes the required procedures in one or more Fortran source
files that are present in the working directory of the <span style="font-family:monospace">WHIZARD</span>