Skip to content

GitLab

Commit e74d0340 authored by Juergen Reuter's avatar Juergen Reuter
Browse files

24e7da9 Merge branch '418-powheg-matching-eeqq-ubf' into 'master'

parent a04f7be2
Expand all Hide whitespace changes
Inline Side-by-side
ChangeLog
View file @ e74d0340
......@@ -4,6 +4,12 @@ Use svn log to see detailed changes.
Version 3.0.0+
2021-07-02
Improvement for POWHEG matching:
- implement massless recoil case
- allow reading in existing POWHEG grids
- support kinematic cuts at generator level
2021-07-01
Distinguish different cases of photons in NLO EW corrections
......
share/doc/manual.tex
View file @ e74d0340
......@@ -940,7 +940,7 @@ separate directory with the appropriate name.
If you want to install the latest development version, you have to
check it out from the \whizard\ SVN repository. Note that since a
couple of years our development is now via a Git revision control
system hosted at the University of Siegen, cf. the next subsection.
system hosted at the University of Siegen, cf. the next subsection.
 
In addition to the prerequisites listed in the previous section, you
need:
......@@ -1022,7 +1022,7 @@ as an official distribution are available from
\url{https://whizard.tp.nt.uni-siegen.de/}
\end{quote}
Building \whizard\ works the way as described in
Sec.~\ref{sec:tarballdistr}.
Sec.~\ref{sec:tarballdistr}.
 
 
%%%%%
......@@ -1327,7 +1327,7 @@ so \whizard\ offers the possibility of backwards compatibility here.
\begin{quote}
configure: --------------------------------------------------------------
configure: --- SHOWERS PYTHIA6 PYTHIA8 MPI ---
configure:
configure:
checking whether we want to enable PYTHIA6... yes
checking for PYTHIA6... (enabled)
checking for PYTHIA6 eh settings... (disabled)
......@@ -1797,7 +1797,7 @@ disable}. However, we do not recommend to do so.
For Windows, from \ttt{Windows 10} onwards, there is the possibility
to install and use an underlying Linux operating system,
e.g. \ttt{Ubuntu}. Installation and usage of \whizard\ works then the
same way as described above.
same way as described above.
 
%%%%%%%%%%%
 
......@@ -9328,7 +9328,7 @@ events in the output file. Thus, the sum of all event weights correctly
reproduces the total cross section.
 
\subsection{\powheg\ matching}
To match the NLO computation with a parton shower, \whizard\ supports
To match the NLO events with a parton shower, \whizard\ supports
the \powheg\ matching. It generates a distribution according to
\begin{align}
\label{eq:powheg}
......@@ -9361,9 +9361,11 @@ interfaced shower. The maximal scale is currently given by \ttt{sqrts}
but should in the future be changeable with \ttt{powheg\_pt\_max}.
 
Note that the \powheg\ event generation needs an additional grid for
efficient event generation that is automatically generated during
integration. Further options that steer the efficiency of this grid are
\ttt{powheg\_grid\_size\_xi} and \ttt{powheg\_grid\_size\_y}.
efficient event generation that is generated during integration if
\ttt{?powheg\_matching = true} is set. Thus, this needs to be set before
the \ttt{integrate} statement.
Further options that steer the efficiency of this grid are
\ttt{powheg\_grid\_size\_xi}, \ttt{powheg\_grid\_size\_y} and \ttt{powheg\_grid\_sampling\_points}.
 
\subsection{Separation of finite and singular contributions}
 
......@@ -9374,7 +9376,7 @@ and singular contributions with the flag
?nlo_use_real_partition = true
\end{code}
The finite contributions, which by definition should not contain soft or
collinear emissions, will then integrate like a ordinary LO integration
collinear emissions, will then integrate like an ordinary LO integration
with one additional particle. Similarly, the event generation will
produce only real events without subtraction terms with Born kinematics
for this additional finite component. The \powheg\ event generation
......@@ -14347,19 +14349,19 @@ must be performed:
 
The \whizard\ subsystem will work with input and output
files in the current working directory, unless asked to do otherwise.
\item
The \ttt{tirpc} library is used by the \ttt{StdHEP} subsystem for \ttt{xdr}
functionality. This
library should be present on the host system.
 
\item
Run the program.
Run the program.
\end{enumerate}
Below is an example program, similar to \whizard's internal unit-test
suite for different external programming languages. The user program
controls the \whizard\ workflow in the same way as a \sindarin\ script
would do. The commands are a mixture of \sindarin\ command
would do. The commands are a mixture of \sindarin\ command
calls and functionality for passing information between the \whizard\
subsystem and the host program.
In particular, the program can process generated events one-by-one.
......@@ -14415,7 +14417,7 @@ for it in range(it_begin, it_end + 1):
 
sample.close()
 
del(wz)
del(wz)
\end{code}
 
\subsubsection{Python module import}
......@@ -14464,11 +14466,11 @@ and string values map to \python\ strings.
To set a \sindarin\ variable of appropriate type:
\begin{quote}
\tt wz.set\_int ( \textit{name}, \textit{value} );
\\
\\
\tt wz.set\_double ( \textit{name}, \textit{value} );
\\
\\
\tt wz.set\_bool ( \textit{name}, \textit{value} );
\\
\\
\tt wz.set\_string ( \textit{name}, \textit{value} );
\end{quote}
\textit{name} is a \python\ string value. It must match the corresponding
......@@ -19020,9 +19022,9 @@ doi:10.1016/j.cpc.2013.05.021
C.~Bierlich, A.~Buckley, J.~Butterworth, C.~H.~Christensen, L.~Corpe,
D.~Grellscheid, J.~F.~Grosse-Oetringhaus, C.~Gutschow,
P.~Karczmarczyk, J.~Klein, L.~L\"onnblad, C.~S.~Pollard, P.~Richardson,
H.~Schulz and F.~Siegert,
H.~Schulz and F.~Siegert,
{\em Robust Independent Validation of Experiment and Theory: Rivet
version 3},
version 3},
SciPost Phys. \textbf{8}, 026 (2020)
doi:10.21468/SciPostPhys.8.2.026
[arXiv:1912.05451 [hep-ph]].
......
share/tests/ext_tests_nlo_add/nlo_qq_powheg.sin
View file @ e74d0340
......@@ -39,7 +39,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 2
powheg_grid_size_y = 2
powheg_grid_sampling_points = 10000
powheg_pt_min = 1.0
?powheg_use_singular_jacobian = false
......
share/tests/ext_tests_nlo_add/nlo_tt_powheg.sin
View file @ e74d0340
......@@ -41,7 +41,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000000
powheg_pt_min = 1
?powheg_use_singular_jacobian = false
......
share/tests/ext_tests_nlo_add/nlo_tt_powheg_sudakov.sin
View file @ e74d0340
......@@ -41,7 +41,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000000
powheg_pt_min = 1
?powheg_use_singular_jacobian = false
......
share/tests/ext_tests_nlo_add/nlo_uu_powheg.sin
View file @ e74d0340
......@@ -41,7 +41,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000000
powheg_pt_min = 1
?powheg_use_singular_jacobian = false
......
share/tests/functional_tests/powheg_1.sin
View file @ e74d0340
# SINDARIN input for WHIZARD self-test
# Testing POWHEG matching for NLO matrix elements with shower
# Testing POWHEG matching for the process ee -> tt
# The second part of this test tests powheg matching together with a
# separation of the real contribution into finite and singular part.
model = SM
?logging = true
......@@ -25,7 +27,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000
powheg_pt_min = 1
?powheg_use_singular_jacobian = false
......
share/tests/functional_tests/ref-output-double/powheg_1.ref
View file @ e74d0340
......@@ -16,7 +16,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000
powheg_pt_min = 1.00000E+00
?powheg_use_singular_jacobian = false
| Process library 'powheg_1_lib': recorded process 'powheg_1_p1'
......@@ -103,13 +102,14 @@ n_events = 1
| Events: writing to ASCII file 'powheg_1_p1.debug'
| Events: generating 1 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p1_1000.pg'
| POWHEG: using grids from file 'powheg_1_p1.pg'
| Grid: Mean value of the grid: 1.702238298E-02
| Grid: Max value of the grid: 3.968172459E-02
| Grid: Mean/Max value of the grid: 4.289728624E-01
| ... event sample complete.
| Events: actual unweighting efficiency = 20.00 %
| Events: closing ASCII file 'powheg_1_p1.debug'
Warning: There have been 1 ( 50.00%) POWHEG grid excesses.
n_events = 2
| Starting simulation for process 'powheg_1_p1'
| Simulate: using integration grids from file 'powheg_1_p1.m1.vg'
......@@ -124,13 +124,14 @@ n_events = 2
| Events: writing to ASCII file 'powheg_1_p1.debug'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p1_1000.pg'
| POWHEG: using grids from file 'powheg_1_p1.pg'
| Grid: Mean value of the grid: 1.702238298E-02
| Grid: Max value of the grid: 3.968172459E-02
| Grid: Mean/Max value of the grid: 4.289728624E-01
| ... event sample complete.
| Events: actual unweighting efficiency = 100.00 %
| Events: closing ASCII file 'powheg_1_p1.debug'
Warning: There have been 1 ( 11.11%) POWHEG grid excesses.
?nlo_use_real_partition = true
real_partition_scale = 5.00000E+00
| Process library 'powheg_1_lib': unloading
......@@ -241,7 +242,7 @@ tolerance = 1.08344E+02
| Events: writing to ASCII file 'powheg_1_p2.debug'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p2_1000.pg'
| POWHEG: using grids from file 'powheg_1_p2.pg'
| Grid: Mean value of the grid: 1.509359941E-03
| Grid: Max value of the grid: 1.846389141E-02
| Grid: Mean/Max value of the grid: 8.174657807E-02
......@@ -250,7 +251,7 @@ tolerance = 1.08344E+02
| Events: closing ASCII file 'powheg_1_p2.debug'
| Summary of value checks:
| Failures: 0 / Total: 1
| There were no errors and 2 warning(s).
| There were no errors and 4 warning(s).
| WHIZARD run finished.
|=============================================================================|
Contents of powheg_1_p1.debug:
......@@ -321,9 +322,9 @@ Contents of powheg_1_p1.debug:
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [v] t 1 0 1-2 5-7 250.000 96.085 -123.677 -85.712 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -96.085 123.677 85.712 30625.000
5 [o] t 2 0 3-4 [none] 160.829 104.543 -121.712 11.079 30625.000
6 [o] tbar 0 1 3-4 [none] 244.953 -92.243 118.732 82.285 30625.000
7 [o] gl 1 2 3-4 [none] 94.218 -12.299 2.981 -93.364 0.000
5 [o] t 2 0 3-4 [none] 214.865 125.308 -146.439 -94.976 30625.000
6 [o] tbar 0 1 3-4 [none] 265.700 -107.599 138.497 95.983 30625.000
7 [o] gl 1 2 3-4 [none] 19.435 -17.709 7.942 -1.007 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -367,9 +368,9 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 1.6082874E+02; 1.0454291E+02,-1.2171243E+02, 1.1079156E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.4495304E+02;-9.2243471E+01, 1.1873192E+02, 8.2285257E+01| 3.0625000E+04| 4)
5 prt(o:21| 9.4218222E+01;-1.2299435E+01, 2.9805095E+00,-9.3364414E+01| 0.0000000E+00| 5)
3 prt(o:6| 2.1486485E+02; 1.2530850E+02,-1.4643856E+02,-9.4975950E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.6570049E+02;-1.0759908E+02, 1.3849702E+02, 9.5983140E+01| 3.0625000E+04| 4)
5 prt(o:21| 1.9434661E+01;-1.7709417E+01, 7.9415454E+00,-1.0071898E+00| 0.0000000E+00| 5)
========================================================================
========================================================================
Event #2
......@@ -436,8 +437,11 @@ event_excess* => 0.00000E+00
Nr Status Flavor Col ACol Parents Children P(0) P(1) P(2) P(3) P^2
1 [i] e+ 0 0 [none] 3-4 250.000 0.000 0.000 250.000 0.000
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [o] t 1 0 1-2 [none] 250.000 168.867 30.937 -49.010 30625.000
4 [o] tbar 0 1 1-2 [none] 250.000 -168.867 -30.937 49.010 30625.000
3 [v] t 1 0 1-2 5-7 250.000 168.867 30.937 -49.010 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -168.867 -30.937 49.010 30625.000
5 [o] t 1 0 3-4 [none] 251.490 170.834 31.297 -49.581 30625.000
6 [o] tbar 0 2 3-4 [none] 198.733 -166.680 -59.028 90.710 30625.000
7 [o] gl 2 1 3-4 [none] 49.778 -4.154 27.731 -41.128 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -467,8 +471,8 @@ Sum of outgoing momenta: p(0:3) = 500.000 0.000 0.000
sqrts* = 5.00000E+02
sqrts_hat* => 5.00000E+02
n_in* => 2
n_out* => 2
n_tot* => 4
n_out* => 3
n_tot* => 5
$process_id* => "powheg_1_p1"
process_num_id* => [unknown integer]
sqme* => 2.76793E-02
......@@ -481,8 +485,9 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 2.5000000E+02; 1.6886659E+02, 3.0936816E+01,-4.9010068E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.5000000E+02;-1.6886659E+02,-3.0936816E+01, 4.9010068E+01| 3.0625000E+04| 4)
3 prt(o:6| 2.5148982E+02; 1.7083419E+02, 3.1297285E+01,-4.9581122E+01| 3.0625000E+04| 3)
4 prt(o:-6| 1.9873263E+02;-1.6667970E+02,-5.9028135E+01, 9.0709516E+01| 3.0625000E+04| 4)
5 prt(o:21| 4.9777552E+01;-4.1544875E+00, 2.7730851E+01,-4.1128394E+01| 0.0000000E+00| 5)
========================================================================
Contents of powheg_1_p2.debug:
========================================================================
......@@ -670,8 +675,11 @@ event_excess* => 0.00000E+00
Nr Status Flavor Col ACol Parents Children P(0) P(1) P(2) P(3) P^2
1 [i] e+ 0 0 [none] 3-4 250.000 0.000 0.000 250.000 0.000
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [o] t 1 0 1-2 [none] 250.000 114.315 2.823 -137.110 30625.000
4 [o] tbar 0 1 1-2 [none] 250.000 -114.315 -2.823 137.110 30625.000
3 [v] t 1 0 1-2 5-7 250.000 114.315 2.823 -137.110 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -114.315 -2.823 137.110 30625.000
5 [o] t 2 0 3-4 [none] 218.688 140.896 4.102 -167.200 30625.000
6 [o] tbar 0 1 3-4 [none] 279.991 -139.945 -3.456 167.850 30625.000
7 [o] gl 1 2 3-4 [none] 1.320 -0.951 -0.646 -0.650 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -701,8 +709,8 @@ Sum of outgoing momenta: p(0:3) = 500.000 0.000 0.000
sqrts* = 5.00000E+02
sqrts_hat* => 5.00000E+02
n_in* => 2
n_out* => 2
n_tot* => 4
n_out* => 3
n_tot* => 5
$process_id* => "powheg_1_p2"
process_num_id* => [unknown integer]
sqme* => 3.51646E-02
......@@ -715,14 +723,15 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 2.5000000E+02; 1.1431539E+02, 2.8231580E+00,-1.3710953E+02| 3.0625000E+04| 3)
4 prt(o:-6| 2.5000000E+02;-1.1431539E+02,-2.8231580E+00, 1.3710953E+02| 3.0625000E+04| 4)
3 prt(o:6| 2.1868819E+02; 1.4089623E+02, 4.1019634E+00,-1.6720033E+02| 3.0625000E+04| 3)
4 prt(o:-6| 2.7999141E+02;-1.3994522E+02,-3.4561180E+00, 1.6784987E+02| 3.0625000E+04| 4)
5 prt(o:21| 1.3203955E+00;-9.5101202E-01,-6.4584531E-01,-6.4954149E-01| 0.0000000E+00| 5)
========================================================================
Contents of powheg_1_p1_1000.pg:
Contents of powheg_1_p1.pg:
3
5 5 2
2.538608113E-02 3.429289805E-02 2.985681667E-02 2.158620136E-02 1.879672179E-02 3.711364208E-02 3.824737015E-02 3.256180673E-02 2.916146820E-02 1.778642354E-02 2.207064596E-02 2.618968673E-02 2.480053435E-02 1.952907052E-02 1.107019080E-02 1.008267321E-02 1.126668430E-02 8.570795734E-03 6.226312737E-03 1.322538637E-03 1.046999756E-03 1.760766185E-03 1.612180666E-03 1.471504660E-03 5.574793422E-04 2.526311554E-02 3.491502186E-02 2.956010450E-02 2.935704230E-02 8.353646513E-03 3.832502345E-02 3.968172459E-02 3.016863471E-02 2.385262853E-02 1.816909461E-02 2.179073657E-02 2.445579634E-02 2.209286125E-02 1.269046939E-02 1.663231625E-02 1.033479411E-02 1.102169006E-02 9.497679173E-03 3.893308452E-03 2.479969318E-03 1.121037866E-03 1.538909193E-03 1.772675648E-03 1.201811393E-03 5.815643085E-04
Contents of powheg_1_p2_1000.pg:
2.538608113E-02 3.429289805E-02 2.985681667E-02 2.158620136E-02 1.879672179E-02 3.711364208E-02 3.824737015E-02 3.256180673E-02 2.916146820E-02 1.778642354E-02 2.207064596E-02 2.618968673E-02 2.480053435E-02 1.952907052E-02 1.107019080E-02 1.008267321E-02 1.126668430E-02 8.570795734E-03 6.226312737E-03 1.322538637E-03 1.046999756E-03 1.760766185E-03 1.612180666E-03 1.471504660E-03 5.574793422E-04 2.526311554E-02 3.491502186E-02 2.956010450E-02 2.935704230E-02 8.353646513E-03 3.832502345E-02 3.968172459E-02 3.016863471E-02 2.385262853E-02 1.816909461E-02 2.179073657E-02 2.445579634E-02 2.209286125E-02 1.269046939E-02 1.663231625E-02 1.033479411E-02 1.102169006E-02 9.497679173E-03 3.893308452E-03 2.479969318E-03 1.121037866E-03 1.538909193E-03 1.772675648E-03 1.201811393E-03 5.815643085E-04
Contents of powheg_1_p2.pg:
3
5 5 2
1.845110429E-02 0.000000000E+00 3.846058726E-04 0.000000000E+00 0.000000000E+00 1.703644571E-02 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.630346817E-03 1.959254406E-04 3.587897845E-05 0.000000000E+00 2.654261351E-06 1.846389141E-02 0.000000000E+00 3.802521236E-04 0.000000000E+00 0.000000000E+00 1.704619181E-02 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.588756600E-03 2.100838567E-04 3.905383913E-05 0.000000000E+00 2.806018250E-06
1.845110429E-02 0.000000000E+00 3.846058726E-04 0.000000000E+00 0.000000000E+00 1.703644571E-02 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.630346817E-03 1.959254406E-04 3.587897845E-05 0.000000000E+00 2.654261351E-06 1.846389141E-02 0.000000000E+00 3.802521236E-04 0.000000000E+00 0.000000000E+00 1.704619181E-02 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.588756600E-03 2.100838567E-04 3.905383913E-05 0.000000000E+00 2.806018250E-06
share/tests/functional_tests/ref-output-ext/powheg_1.ref
View file @ e74d0340
......@@ -16,7 +16,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000
powheg_pt_min = 1.00000E+00
?powheg_use_singular_jacobian = false
| Process library 'powheg_1_lib': recorded process 'powheg_1_p1'
......@@ -103,13 +102,14 @@ n_events = 1
| Events: writing to ASCII file 'powheg_1_p1.debug'
| Events: generating 1 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p1_1000.pg'
| POWHEG: using grids from file 'powheg_1_p1.pg'
| Grid: Mean value of the grid: 1.702238736E-02
| Grid: Max value of the grid: 3.968173850E-02
| Grid: Mean/Max value of the grid: 4.289728224E-01
| ... event sample complete.
| Events: actual unweighting efficiency = 20.00 %
| Events: closing ASCII file 'powheg_1_p1.debug'
Warning: There have been 1 ( 50.00%) POWHEG grid excesses.
n_events = 2
| Starting simulation for process 'powheg_1_p1'
| Simulate: using integration grids from file 'powheg_1_p1.m1.vg'
......@@ -124,13 +124,14 @@ n_events = 2
| Events: writing to ASCII file 'powheg_1_p1.debug'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p1_1000.pg'
| POWHEG: using grids from file 'powheg_1_p1.pg'
| Grid: Mean value of the grid: 1.702238736E-02
| Grid: Max value of the grid: 3.968173850E-02
| Grid: Mean/Max value of the grid: 4.289728224E-01
| ... event sample complete.
| Events: actual unweighting efficiency = 100.00 %
| Events: closing ASCII file 'powheg_1_p1.debug'
Warning: There have been 1 ( 11.11%) POWHEG grid excesses.
?nlo_use_real_partition = true
real_partition_scale = 5.00000E+00
| Process library 'powheg_1_lib': unloading
......@@ -241,7 +242,7 @@ tolerance = 1.08344E+02
| Events: writing to ASCII file 'powheg_1_p2.debug'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p2_1000.pg'
| POWHEG: using grids from file 'powheg_1_p2.pg'
| Grid: Mean value of the grid: 1.509352292E-03
| Grid: Max value of the grid: 1.846388060E-02
| Grid: Mean/Max value of the grid: 8.174621169E-02
......@@ -250,7 +251,7 @@ tolerance = 1.08344E+02
| Events: closing ASCII file 'powheg_1_p2.debug'
| Summary of value checks:
| Failures: 0 / Total: 1
| There were no errors and 2 warning(s).
| There were no errors and 4 warning(s).
| WHIZARD run finished.
|=============================================================================|
Contents of powheg_1_p1.debug:
......@@ -321,9 +322,9 @@ Contents of powheg_1_p1.debug:
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [v] t 1 0 1-2 5-7 250.000 96.085 -123.677 -85.712 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -96.085 123.677 85.712 30625.000
5 [o] t 2 0 3-4 [none] 160.829 104.543 -121.712 11.079 30625.000
6 [o] tbar 0 1 3-4 [none] 244.953 -92.243 118.732 82.285 30625.000
7 [o] gl 1 2 3-4 [none] 94.218 -12.299 2.981 -93.364 0.000
5 [o] t 2 0 3-4 [none] 214.865 125.308 -146.439 -94.976 30625.000
6 [o] tbar 0 1 3-4 [none] 265.700 -107.599 138.497 95.983 30625.000
7 [o] gl 1 2 3-4 [none] 19.435 -17.709 7.942 -1.007 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -367,9 +368,9 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 1.6082874E+02; 1.0454291E+02,-1.2171243E+02, 1.1079156E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.4495304E+02;-9.2243471E+01, 1.1873192E+02, 8.2285257E+01| 3.0625000E+04| 4)
5 prt(o:21| 9.4218222E+01;-1.2299435E+01, 2.9805095E+00,-9.3364414E+01| 0.0000000E+00| 5)
3 prt(o:6| 2.1486485E+02; 1.2530850E+02,-1.4643856E+02,-9.4975950E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.6570049E+02;-1.0759908E+02, 1.3849702E+02, 9.5983140E+01| 3.0625000E+04| 4)
5 prt(o:21| 1.9434661E+01;-1.7709417E+01, 7.9415454E+00,-1.0071898E+00| 0.0000000E+00| 5)
========================================================================
========================================================================
Event #2
......@@ -436,8 +437,11 @@ event_excess* => 0.00000E+00
Nr Status Flavor Col ACol Parents Children P(0) P(1) P(2) P(3) P^2
1 [i] e+ 0 0 [none] 3-4 250.000 0.000 0.000 250.000 0.000
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [o] t 1 0 1-2 [none] 250.000 168.867 30.937 -49.010 30625.000
4 [o] tbar 0 1 1-2 [none] 250.000 -168.867 -30.937 49.010 30625.000
3 [v] t 1 0 1-2 5-7 250.000 168.867 30.937 -49.010 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -168.867 -30.937 49.010 30625.000
5 [o] t 1 0 3-4 [none] 251.490 170.834 31.297 -49.581 30625.000
6 [o] tbar 0 2 3-4 [none] 198.733 -166.680 -59.028 90.710 30625.000
7 [o] gl 2 1 3-4 [none] 49.778 -4.154 27.731 -41.128 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -467,8 +471,8 @@ Sum of outgoing momenta: p(0:3) = 500.000 0.000 0.000
sqrts* = 5.00000E+02
sqrts_hat* => 5.00000E+02
n_in* => 2
n_out* => 2
n_tot* => 4
n_out* => 3
n_tot* => 5
$process_id* => "powheg_1_p1"
process_num_id* => [unknown integer]
sqme* => 2.76793E-02
......@@ -481,8 +485,9 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 2.5000000E+02; 1.6886659E+02, 3.0936816E+01,-4.9010068E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.5000000E+02;-1.6886659E+02,-3.0936816E+01, 4.9010068E+01| 3.0625000E+04| 4)
3 prt(o:6| 2.5148982E+02; 1.7083419E+02, 3.1297285E+01,-4.9581122E+01| 3.0625000E+04| 3)
4 prt(o:-6| 1.9873263E+02;-1.6667970E+02,-5.9028135E+01, 9.0709516E+01| 3.0625000E+04| 4)
5 prt(o:21| 4.9777552E+01;-4.1544875E+00, 2.7730851E+01,-4.1128394E+01| 0.0000000E+00| 5)
========================================================================
Contents of powheg_1_p2.debug:
========================================================================
......@@ -670,8 +675,11 @@ event_excess* => 0.00000E+00
Nr Status Flavor Col ACol Parents Children P(0) P(1) P(2) P(3) P^2
1 [i] e+ 0 0 [none] 3-4 250.000 0.000 0.000 250.000 0.000
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [o] t 1 0 1-2 [none] 250.000 114.315 2.823 -137.110 30625.000
4 [o] tbar 0 1 1-2 [none] 250.000 -114.315 -2.823 137.110 30625.000
3 [v] t 1 0 1-2 5-7 250.000 114.315 2.823 -137.110 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -114.315 -2.823 137.110 30625.000
5 [o] t 2 0 3-4 [none] 218.688 140.896 4.102 -167.200 30625.000
6 [o] tbar 0 1 3-4 [none] 279.991 -139.945 -3.456 167.850 30625.000
7 [o] gl 1 2 3-4 [none] 1.320 -0.951 -0.646 -0.650 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -701,8 +709,8 @@ Sum of outgoing momenta: p(0:3) = 500.000 0.000 0.000
sqrts* = 5.00000E+02
sqrts_hat* => 5.00000E+02
n_in* => 2
n_out* => 2
n_tot* => 4
n_out* => 3
n_tot* => 5
$process_id* => "powheg_1_p2"
process_num_id* => [unknown integer]
sqme* => 3.51646E-02
......@@ -715,14 +723,15 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 2.5000000E+02; 1.1431539E+02, 2.8231580E+00,-1.3710953E+02| 3.0625000E+04| 3)
4 prt(o:-6| 2.5000000E+02;-1.1431539E+02,-2.8231580E+00, 1.3710953E+02| 3.0625000E+04| 4)
3 prt(o:6| 2.1868819E+02; 1.4089623E+02, 4.1019634E+00,-1.6720033E+02| 3.0625000E+04| 3)
4 prt(o:-6| 2.7999141E+02;-1.3994522E+02,-3.4561180E+00, 1.6784987E+02| 3.0625000E+04| 4)
5 prt(o:21| 1.3203955E+00;-9.5101202E-01,-6.4584531E-01,-6.4954149E-01| 0.0000000E+00| 5)
========================================================================
Contents of powheg_1_p1_1000.pg:
Contents of powheg_1_p1.pg:
3
5 5 2
2.538608541E-02 3.429290764E-02 2.985682734E-02 2.158620934E-02 1.879672724E-02 3.711365002E-02 3.824738385E-02 3.256181610E-02 2.916147563E-02 1.778642711E-02 2.207065111E-02 2.618969324E-02 2.480054100E-02 1.952907677E-02 1.107019458E-02 1.008267038E-02 1.126668736E-02 8.570798820E-03 6.226314717E-03 1.322539021E-03 1.047000011E-03 1.760766746E-03 1.612181296E-03 1.471505200E-03 5.574795380E-04 2.526311995E-02 3.491503111E-02 2.956011518E-02 2.935705104E-02 8.353648802E-03 3.832503227E-02 3.968173850E-02 3.016864348E-02 2.385263511E-02 1.816910045E-02 2.179073172E-02 2.445580372E-02 2.209286763E-02 1.269047307E-02 1.663232008E-02 1.033479040E-02 1.102169336E-02 9.497682272E-03 3.893309967E-03 2.479970072E-03 1.121038054E-03 1.538909782E-03 1.772676298E-03 1.201811894E-03 5.815645111E-04
Contents of powheg_1_p2_1000.pg:
Contents of powheg_1_p2.pg:
3
5 5 2
1.845109355E-02 0.000000000E+00 3.846060261E-04 0.000000000E+00 0.000000000E+00 1.703626528E-02 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.630346745E-03 1.959255008E-04 3.587899256E-05 0.000000000E+00 2.654262470E-06 1.846388060E-02 0.000000000E+00 3.802522765E-04 0.000000000E+00 0.000000000E+00 1.704601098E-02 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 0.000000000E+00 1.588756609E-03 2.100839106E-04 3.905385338E-05 0.000000000E+00 2.806019413E-06
share/tests/functional_tests/ref-output-quad/powheg_1.ref
View file @ e74d0340
......@@ -16,7 +16,6 @@ alphas_power = 0
?powheg_matching = true
powheg_grid_size_xi = 5
powheg_grid_size_y = 5
powheg_grid_sampling_points = 1000
powheg_pt_min = 1.00000E+00
?powheg_use_singular_jacobian = false
| Process library 'powheg_1_lib': recorded process 'powheg_1_p1'
......@@ -103,13 +102,14 @@ n_events = 1
| Events: writing to ASCII file 'powheg_1_p1.debug'
| Events: generating 1 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p1_1000.pg'
| POWHEG: using grids from file 'powheg_1_p1.pg'
| Grid: Mean value of the grid: 1.702238737E-02
| Grid: Max value of the grid: 3.968173851E-02
| Grid: Mean/Max value of the grid: 4.289728224E-01
| ... event sample complete.
| Events: actual unweighting efficiency = 20.00 %
| Events: closing ASCII file 'powheg_1_p1.debug'
Warning: There have been 1 ( 50.00%) POWHEG grid excesses.
n_events = 2
| Starting simulation for process 'powheg_1_p1'
| Simulate: using integration grids from file 'powheg_1_p1.m1.vg'
......@@ -124,13 +124,14 @@ n_events = 2
| Events: writing to ASCII file 'powheg_1_p1.debug'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p1_1000.pg'
| POWHEG: using grids from file 'powheg_1_p1.pg'
| Grid: Mean value of the grid: 1.702238737E-02
| Grid: Max value of the grid: 3.968173851E-02
| Grid: Mean/Max value of the grid: 4.289728224E-01
| ... event sample complete.
| Events: actual unweighting efficiency = 100.00 %
| Events: closing ASCII file 'powheg_1_p1.debug'
Warning: There have been 1 ( 11.11%) POWHEG grid excesses.
?nlo_use_real_partition = true
real_partition_scale = 5.00000E+00
| Process library 'powheg_1_lib': unloading
......@@ -241,7 +242,7 @@ tolerance = 1.08344E+02
| Events: writing to ASCII file 'powheg_1_p2.debug'
| Events: generating 2 unweighted, unpolarized events ...
| Events: event normalization mode '1'
| POWHEG: using grids from file 'powheg_1_p2_1000.pg'
| POWHEG: using grids from file 'powheg_1_p2.pg'
| Grid: Mean value of the grid: 1.509352289E-03
| Grid: Max value of the grid: 1.846388060E-02
| Grid: Mean/Max value of the grid: 8.174621150E-02
......@@ -250,7 +251,7 @@ tolerance = 1.08344E+02
| Events: closing ASCII file 'powheg_1_p2.debug'
| Summary of value checks:
| Failures: 0 / Total: 1
| There were no errors and 2 warning(s).
| There were no errors and 4 warning(s).
| WHIZARD run finished.
|=============================================================================|
Contents of powheg_1_p1.debug:
......@@ -321,9 +322,9 @@ Contents of powheg_1_p1.debug:
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [v] t 1 0 1-2 5-7 250.000 96.085 -123.677 -85.712 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -96.085 123.677 85.712 30625.000
5 [o] t 2 0 3-4 [none] 160.829 104.543 -121.712 11.079 30625.000
6 [o] tbar 0 1 3-4 [none] 244.953 -92.243 118.732 82.285 30625.000
7 [o] gl 1 2 3-4 [none] 94.218 -12.299 2.981 -93.364 0.000
5 [o] t 2 0 3-4 [none] 214.865 125.308 -146.439 -94.976 30625.000
6 [o] tbar 0 1 3-4 [none] 265.700 -107.599 138.497 95.983 30625.000
7 [o] gl 1 2 3-4 [none] 19.435 -17.709 7.942 -1.007 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -367,9 +368,9 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 1.6082874E+02; 1.0454291E+02,-1.2171243E+02, 1.1079156E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.4495304E+02;-9.2243471E+01, 1.1873192E+02, 8.2285257E+01| 3.0625000E+04| 4)
5 prt(o:21| 9.4218222E+01;-1.2299435E+01, 2.9805095E+00,-9.3364414E+01| 0.0000000E+00| 5)
3 prt(o:6| 2.1486485E+02; 1.2530850E+02,-1.4643856E+02,-9.4975950E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.6570049E+02;-1.0759908E+02, 1.3849702E+02, 9.5983140E+01| 3.0625000E+04| 4)
5 prt(o:21| 1.9434661E+01;-1.7709417E+01, 7.9415454E+00,-1.0071898E+00| 0.0000000E+00| 5)
========================================================================
========================================================================
Event #2
......@@ -436,8 +437,11 @@ event_excess* => 0.00000E+00
Nr Status Flavor Col ACol Parents Children P(0) P(1) P(2) P(3) P^2
1 [i] e+ 0 0 [none] 3-4 250.000 0.000 0.000 250.000 0.000
2 [i] e- 0 0 [none] 3-4 250.000 0.000 0.000 -250.000 0.000
3 [o] t 1 0 1-2 [none] 250.000 168.867 30.937 -49.010 30625.000
4 [o] tbar 0 1 1-2 [none] 250.000 -168.867 -30.937 49.010 30625.000
3 [v] t 1 0 1-2 5-7 250.000 168.867 30.937 -49.010 30625.000
4 [v] tbar 0 1 1-2 5-7 250.000 -168.867 -30.937 49.010 30625.000
5 [o] t 1 0 3-4 [none] 251.490 170.834 31.297 -49.581 30625.000
6 [o] tbar 0 2 3-4 [none] 198.733 -166.680 -59.028 90.710 30625.000
7 [o] gl 2 1 3-4 [none] 49.778 -4.154 27.731 -41.128 0.000
------------------------------------------------------------------------
Sum of incoming momenta: p(0:3) = 500.000 0.000 0.000 0.000
Sum of beam remnant momenta: p(0:3) = 0.000 0.000 0.000 0.000
......@@ -467,8 +471,8 @@ Sum of outgoing momenta: p(0:3) = 500.000 0.000 0.000
sqrts* = 5.00000E+02
sqrts_hat* => 5.00000E+02
n_in* => 2
n_out* => 2
n_tot* => 4
n_out* => 3
n_tot* => 5
$process_id* => "powheg_1_p1"
process_num_id* => [unknown integer]
sqme* => 2.76793E-02
......@@ -481,8 +485,9 @@ event_excess* => 0.00000E+00
subevent:
1 prt(i:-11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00,-2.5000000E+02| 2.6111793E-07| 1)
2 prt(i:11|-2.5000000E+02; 0.0000000E+00, 0.0000000E+00, 2.5000000E+02| 2.6111793E-07| 2)
3 prt(o:6| 2.5000000E+02; 1.6886659E+02, 3.0936816E+01,-4.9010068E+01| 3.0625000E+04| 3)
4 prt(o:-6| 2.5000000E+02;-1.6886659E+02,-3.0936816E+01, 4.9010068E+01| 3.0625000E+04| 4)
3 prt(o:6| 2.5148982E+02; 1.7083419E+02, 3.1297285E+01,-4.9581122E+01| 3.0625000E+04| 3)
4 prt(o:-6| 1.9873263E+02;-1.6667970E+02,-5.9028135E+01, 9.0709516E+01| 3.0625000E+04| 4)
5 prt(o:21| 4.9777552E+01;-4.1544875E+00, 2.7730851E+01,-4.1128394E+01| 0.0000000E+00| 5)
========================================================================