Index: trunk/share/tests/ext_tests_nlo/nlo_eettjj.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_eettjj.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_eettjj.sin (revision 8292) @@ -1,28 +1,28 @@ include ("nlo_settings.sin") alpha_power = 2 alphas_power = 2 cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 2 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 2 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]]) / 2 elsif njet == 3 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]]) / 2 endif #?combined_nlo_integration = true process nlo_eettjj_p1 = e1, E1 => t, T, jet, jet { nlo_calculation = full } -integrate (nlo_eettjj_p1) { iterations = 1:100:"gw" } +integrate (nlo_eettjj_p1) { iterations = 1:180:"gw" } Index: trunk/share/tests/ext_tests_nlo/nlo_ee5j.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_ee5j.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_ee5j.sin (revision 8292) @@ -1,30 +1,30 @@ include("nlo_settings.sin") cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 5 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 5 then ( eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]] + eval Pt [extract index 4 [@clustered_jets]] + eval Pt [extract index 5 [@clustered_jets]]) / 2 elsif njet == 6 then ( eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]] + eval Pt [extract index 4 [@clustered_jets]] + eval Pt [extract index 5 [@clustered_jets]] + eval Pt [extract index 6 [@clustered_jets]]) / 2 endif alpha_power = 2 alphas_power = 3 #?combined_nlo_integration = true process nlo_ee5j_p1 = e1, E1 => jet, jet, jet, jet, jet { nlo_calculation = full $restrictions = "!W+:W-" } -integrate (nlo_ee5j_p1) { iterations = 1:100:"gw" } +integrate (nlo_ee5j_p1) { iterations = 1:1040:"gw" } Index: trunk/share/tests/ext_tests_nlo/nlo_eetthjj.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_eetthjj.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_eetthjj.sin (revision 8292) @@ -1,30 +1,30 @@ include ("nlo_settings.sin") alpha_power = 3 alphas_power = 2 cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 2 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 2 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + sqrt (eval Pt [extract index 1 [H]]**2 + eval M [extract index 1 [H]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]]) / 2 elsif njet == 3 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + sqrt (eval Pt [extract index 1 [H]]**2 + eval M [extract index 1 [H]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]]) / 2 endif #?combined_nlo_integration = true process nlo_eetthjj_p1 = e1, E1 => t, T, H, jet, jet {nlo_calculation = full } -integrate (nlo_eetthjj_p1) { iterations = 1:100:"gw" } +integrate (nlo_eetthjj_p1) { iterations = 1:380:"gw" } Index: trunk/share/tests/ext_tests_nlo/nlo_eettjjj.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_eettjjj.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_eettjjj.sin (revision 8292) @@ -1,30 +1,30 @@ include ("nlo_settings.sin") alpha_power = 2 alphas_power = 3 cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 3 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 3 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]]) / 2 elsif njet == 4 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]] + eval Pt [extract index 4 [@clustered_jets]]) / 2 endif #?combined_nlo_integration = true process nlo_eettjjj_p1 = e1, E1 => t, T, jet, jet, jet { nlo_calculation = full } -integrate (nlo_eettjjj_p1) { iterations = 1:100:"gw" } +integrate (nlo_eettjjj_p1) { iterations = 1:780:"gw" } Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettaz.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettaz.ref (revision 8291) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettaz.ref (revision 8292) @@ -1,182 +0,0 @@ -?openmp_logging = false -?vis_history = false -?integration_timer = false -openmp_num_threads = 1 -?debug_decay = false -?debug_process = false -?debug_verbose = false -?write_raw = false -| Including Sindarin from 'nlo_settings.sin' -| Switching to model 'SM', scheme 'GF_MW_MZ' -$blha_ew_scheme = "alpha_qed" -SM.mZ => 9.118800000000E+01 -SM.mW => 8.041900200000E+01 -SM.mH => 1.250000000000E+02 -SM.GF => 1.166390000000E-05 -SM.wZ => 0.000000000000E+00 -SM.wtop => 0.000000000000E+00 -SM.wW => 0.000000000000E+00 -SM.wH => 0.000000000000E+00 -SM.ms => 0.000000000000E+00 -SM.mc => 0.000000000000E+00 -SM.mb => 0.000000000000E+00 -SM.mtop => 1.732000000000E+02 -SM.me => 0.000000000000E+00 -SM.mmu => 0.000000000000E+00 -SM.mtau => 1.777000000000E+00 -SM.alphas => 1.180000000000E-01 -?alphas_is_fixed = false -?alphas_from_mz = true -?alphas_from_lambda_qcd = false -alphas_nf = 5 -alphas_order = 2 -[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) -$method = "openloops" -seed = 8131 -sqrts = 1.000000000000E+03 -jet_algorithm = 2 -jet_r = 5.000000000000E-01 -?use_vamp_equivalences = false -| End of included 'nlo_settings.sin' -alpha_power = 4 -alphas_power = 0 -[user variable] theta_0 = 7.000000000000E-01 -[user variable] tops = PDG(6, -6) -| Process library 'nlo_eettaz_lib': recorded process 'nlo_eettaz_p1' -mult_call_real = 2.000000000000E+00 -| Integrate: current process library needs compilation -| Process library 'nlo_eettaz_lib': compiling ... -| Process library 'nlo_eettaz_lib': writing makefile -| Process library 'nlo_eettaz_lib': removing old files -| Process library 'nlo_eettaz_lib': writing driver -| Process library 'nlo_eettaz_lib': creating source code -| Process library 'nlo_eettaz_lib': compiling sources -| Process library 'nlo_eettaz_lib': linking -| Process library 'nlo_eettaz_lib': loading -| Process library 'nlo_eettaz_lib': ... success. -| Integrate: compilation done -| QCD alpha: using a running strong coupling -| RNG: Initializing TAO random-number generator -| RNG: Setting seed for random-number generator to 8131 -| Initializing integration for process nlo_eettaz_p1: -| Beam structure: [any particles] -| Beam data (collision): -| e- (mass = 0.0000000E+00 GeV) -| e+ (mass = 0.0000000E+00 GeV) -| sqrts = 1.000000000000E+03 GeV -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eettaz_p1.i1.phs' -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eettaz_p1.i3.phs' -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| ------------------------------------------------------------------------ -| Process [scattering]: 'nlo_eettaz_p1' -| Library name = 'nlo_eettaz_lib' -| Process index = 1 -| Process components: -| 1: 'nlo_eettaz_p1_i1': e-, e+ => t, tbar, Z, A [openloops] -| 2: 'nlo_eettaz_p1_i2': e-, e+ => t, tbar, Z, A, gl [openloops], [real] -| 3: 'nlo_eettaz_p1_i3': e-, e+ => t, tbar, Z, A [openloops], [virtual] -| 4: 'nlo_eettaz_p1_i4': e-, e+ => t, tbar, Z, A [inactive], [subtraction] -| ------------------------------------------------------------------------ -| Phase space: 28 channels, 8 dimensions -| Phase space: found 28 channels, collected in 6 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 28 channels, 11 dimensions -| Phase space: found 28 channels, collected in 6 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 28 channels, 8 dimensions -| Phase space: found 28 channels, collected in 6 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Applying user-defined cuts. -| Using user-defined general scale. -| Starting integration for process 'nlo_eettaz_p1' part 'born' -| Integrate: iterations = 4:4000:"gw", 2:5000 -| Integrator: 6 chains, 28 channels, 8 dimensions -| Integrator: 4000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 3976 2.1555586E-01 1.34E-02 6.22 3.92* 3.49 - 2 3962 2.1897091E-01 6.05E-03 2.76 1.74* 9.01 - 3 3948 2.3280371E-01 5.47E-03 2.35 1.48* 11.32 - 4 3932 2.0487806E-01 4.40E-03 2.15 1.35* 12.93 -|-----------------------------------------------------------------------------| - 4 15818 2.1656257E-01 2.91E-03 1.35 1.69 12.93 5.34 4 -|-----------------------------------------------------------------------------| - 5 4988 2.2201713E-01 4.19E-03 1.89 1.33* 12.16 - 6 4988 2.2329813E-01 4.25E-03 1.90 1.34 9.57 -|-----------------------------------------------------------------------------| - 6 9976 2.2264865E-01 2.98E-03 1.34 1.34 9.57 0.05 2 -|=============================================================================| -| Starting integration for process 'nlo_eettaz_p1' part 'real' -| Integrate: iterations = 4:8000:"gw", 2:10000 -| Integrator: 6 chains, 28 channels, 11 dimensions -| Integrator: 8000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 7980 -4.5448605E-02 2.48E-03 5.46 4.87* 1.92 - 2 7962 -4.7085440E-02 1.10E-03 2.34 2.09* 5.95 - 3 7944 -4.8114367E-02 9.88E-04 2.05 1.83* 6.86 - 4 7930 -4.6722717E-02 8.08E-04 1.73 1.54* 8.19 -|-----------------------------------------------------------------------------| - 4 31816 -4.7151433E-02 5.32E-04 1.13 2.01 8.19 0.57 4 -|-----------------------------------------------------------------------------| - 5 9990 -4.6686819E-02 7.66E-04 1.64 1.64 6.59 - 6 9990 -4.5178629E-02 8.23E-04 1.82 1.82 4.77 -|-----------------------------------------------------------------------------| - 6 19980 -4.5986417E-02 5.61E-04 1.22 1.72 4.77 1.80 2 -|=============================================================================| -| Starting integration for process 'nlo_eettaz_p1' part 'virtual' -| Integrate: iterations = 4:4000:"gw", 2:5000 -| Integrator: 6 chains, 28 channels, 8 dimensions -| Integrator: 4000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 3976 6.0684933E-02 4.07E-03 6.70 4.23* 3.26 - 2 3960 6.4939267E-02 1.81E-03 2.79 1.75* 9.19 - 3 3948 6.3222357E-02 1.53E-03 2.42 1.52* 9.94 - 4 3936 6.2701722E-02 1.56E-03 2.49 1.56 10.35 -|-----------------------------------------------------------------------------| - 4 15820 6.3352602E-02 9.12E-04 1.44 1.81 10.35 0.46 4 -|-----------------------------------------------------------------------------| - 5 4986 6.4403587E-02 1.32E-03 2.04 1.44* 10.15 - 6 4986 6.3100881E-02 1.26E-03 2.00 1.42* 8.39 -|-----------------------------------------------------------------------------| - 6 9972 6.3726653E-02 9.12E-04 1.43 1.43 8.39 0.51 2 -|=============================================================================| -| Integrate: sum of all components -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 0 2.4038889E-01 3.17E-03 1.32 0.00* 7.79 -| NLO Correction: [O(alpha_s+1)/O(alpha_s)] -| ( 7.9678 +- 0.49239 ) % -|=============================================================================| -| WHIZARD run finished. -|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eebb.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eebb.ref (revision 8291) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eebb.ref (revision 8292) @@ -1,140 +0,0 @@ -?openmp_logging = false -?vis_history = false -?integration_timer = false -openmp_num_threads = 1 -?debug_decay = false -?debug_process = false -?debug_verbose = false -?write_raw = false -| Including Sindarin from 'nlo_settings.sin' -| Switching to model 'SM', scheme 'GF_MW_MZ' -$blha_ew_scheme = "alpha_qed" -SM.mZ => 9.118800000000E+01 -SM.mW => 8.041900200000E+01 -SM.mH => 1.250000000000E+02 -SM.GF => 1.166390000000E-05 -SM.wZ => 0.000000000000E+00 -SM.wtop => 0.000000000000E+00 -SM.wW => 0.000000000000E+00 -SM.wH => 0.000000000000E+00 -SM.ms => 0.000000000000E+00 -SM.mc => 0.000000000000E+00 -SM.mb => 0.000000000000E+00 -SM.mtop => 1.732000000000E+02 -SM.me => 0.000000000000E+00 -SM.mmu => 0.000000000000E+00 -SM.mtau => 1.777000000000E+00 -SM.alphas => 1.180000000000E-01 -?alphas_is_fixed = false -?alphas_from_mz = true -?alphas_from_lambda_qcd = false -alphas_nf = 5 -alphas_order = 2 -[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) -$method = "openloops" -seed = 8131 -sqrts = 1.000000000000E+03 -jet_algorithm = 2 -jet_r = 5.000000000000E-01 -?use_vamp_equivalences = false -| End of included 'nlo_settings.sin' -alpha_power = 2 -alphas_power = 0 -alphas_nf = 4 -SM.mb => 4.750000000000E+00 -[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 21) -$exclude_gauge_splittings = "t:b" -?combined_nlo_integration = true -| Process library 'nlo_eebb_lib': recorded process 'nlo_eebb_p1' -| Integrate: current process library needs compilation -| Process library 'nlo_eebb_lib': compiling ... -| Process library 'nlo_eebb_lib': writing makefile -| Process library 'nlo_eebb_lib': removing old files -| Process library 'nlo_eebb_lib': writing driver -| Process library 'nlo_eebb_lib': creating source code -| Process library 'nlo_eebb_lib': compiling sources -| Process library 'nlo_eebb_lib': linking -| Process library 'nlo_eebb_lib': loading -| Process library 'nlo_eebb_lib': ... success. -| Integrate: compilation done -| QCD alpha: using a running strong coupling -| RNG: Initializing TAO random-number generator -| RNG: Setting seed for random-number generator to 8131 -| Initializing integration for process nlo_eebb_p1: -| Beam structure: [any particles] -| Beam data (collision): -| e- (mass = 0.0000000E+00 GeV) -| e+ (mass = 0.0000000E+00 GeV) -| sqrts = 1.000000000000E+03 GeV -| Phase space: generating configuration ... -Warning: Intermediate decay of zero-width particle Z may be possible. -Warning: Intermediate decay of zero-width particle H may be possible. -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eebb_p1.i1.phs' -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eebb_p1.i3.phs' -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| ------------------------------------------------------------------------ -| Process [scattering]: 'nlo_eebb_p1' -| Library name = 'nlo_eebb_lib' -| Process index = 1 -| Process components: -| 1: 'nlo_eebb_p1_i1': e-, e+ => b, bbar [openloops] -| 2: 'nlo_eebb_p1_i2': e-, e+ => b, bbar, gl [openloops], [real] -| 3: 'nlo_eebb_p1_i3': e-, e+ => b, bbar [openloops], [virtual] -| 4: 'nlo_eebb_p1_i4': e-, e+ => b, bbar [inactive], [subtraction] -| ------------------------------------------------------------------------ -| Phase space: 1 channels, 5 dimensions -| Phase space: found 1 channel, collected in 1 grove. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 1 channels, 5 dimensions -| Phase space: found 1 channel, collected in 1 grove. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 1 channels, 5 dimensions -| Phase space: found 1 channel, collected in 1 grove. -| Phase space: no equivalences between channels used. -| Phase space: wood -Warning: No cuts have been defined. -| Using user-defined general scale. -| Starting integration for process 'nlo_eebb_p1' -| Integrate: iterations = 5:2000:"gw", 2:1000 -| Integrator: 1 chains, 1 channels, 5 dimensions -| Integrator: 2000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 2000 9.4556570E+01 1.62E+00 1.71 0.77* 28.12 - 2 2000 9.4158530E+01 4.87E-01 0.52 0.23* 48.96 - 3 2000 9.4494604E+01 2.76E-01 0.29 0.13* 56.33 - 4 2000 9.5244719E+01 2.66E-01 0.28 0.13* 53.18 - 5 2000 9.4568866E+01 2.56E-01 0.27 0.12* 54.79 -|-----------------------------------------------------------------------------| - 5 10000 9.4713497E+01 1.46E-01 0.15 0.15 54.79 1.56 5 -|-----------------------------------------------------------------------------| - 6 1000 9.3704669E+01 4.00E-01 0.43 0.13 58.52 - 7 1000 9.4573419E+01 4.22E-01 0.45 0.14 59.15 -|-----------------------------------------------------------------------------| - 7 2000 9.4115881E+01 2.90E-01 0.31 0.14 59.15 2.24 2 -|=============================================================================| -| There were no errors and 3 warning(s). -| WHIZARD run finished. -|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eetta.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eetta.ref (revision 8291) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eetta.ref (revision 8292) @@ -1,193 +0,0 @@ -?openmp_logging = false -?vis_history = false -?integration_timer = false -openmp_num_threads = 1 -?debug_decay = false -?debug_process = false -?debug_verbose = false -?write_raw = false -| Including Sindarin from 'nlo_settings.sin' -| Switching to model 'SM', scheme 'GF_MW_MZ' -$blha_ew_scheme = "alpha_qed" -SM.mZ => 9.118800000000E+01 -SM.mW => 8.041900200000E+01 -SM.mH => 1.250000000000E+02 -SM.GF => 1.166390000000E-05 -SM.wZ => 0.000000000000E+00 -SM.wtop => 0.000000000000E+00 -SM.wW => 0.000000000000E+00 -SM.wH => 0.000000000000E+00 -SM.ms => 0.000000000000E+00 -SM.mc => 0.000000000000E+00 -SM.mb => 0.000000000000E+00 -SM.mtop => 1.732000000000E+02 -SM.me => 0.000000000000E+00 -SM.mmu => 0.000000000000E+00 -SM.mtau => 1.777000000000E+00 -SM.alphas => 1.180000000000E-01 -?alphas_is_fixed = false -?alphas_from_mz = true -?alphas_from_lambda_qcd = false -alphas_nf = 5 -alphas_order = 2 -[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) -$method = "openloops" -seed = 8131 -sqrts = 1.000000000000E+03 -jet_algorithm = 2 -jet_r = 5.000000000000E-01 -?use_vamp_equivalences = false -| End of included 'nlo_settings.sin' -alpha_power = 3 -alphas_power = 0 -[user variable] R_0 = 7.000000000000E-01 -[user variable] tops = PDG(6, -6) -| Process library 'nlo_eetta_lib': recorded process 'nlo_eetta_p1' -| Integrate: current process library needs compilation -| Process library 'nlo_eetta_lib': compiling ... -| Process library 'nlo_eetta_lib': writing makefile -| Process library 'nlo_eetta_lib': removing old files -| Process library 'nlo_eetta_lib': writing driver -| Process library 'nlo_eetta_lib': creating source code -| Process library 'nlo_eetta_lib': compiling sources -| Process library 'nlo_eetta_lib': linking -| Process library 'nlo_eetta_lib': loading -| Process library 'nlo_eetta_lib': ... success. -| Integrate: compilation done -| QCD alpha: using a running strong coupling -| RNG: Initializing TAO random-number generator -| RNG: Setting seed for random-number generator to 8131 -| Initializing integration for process nlo_eetta_p1: -| Beam structure: [any particles] -| Beam data (collision): -| e- (mass = 0.0000000E+00 GeV) -| e+ (mass = 0.0000000E+00 GeV) -| sqrts = 1.000000000000E+03 GeV -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eetta_p1.i1.phs' -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eetta_p1.i3.phs' -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| ------------------------------------------------------------------------ -| Process [scattering]: 'nlo_eetta_p1' -| Library name = 'nlo_eetta_lib' -| Process index = 1 -| Process components: -| 1: 'nlo_eetta_p1_i1': e-, e+ => t, tbar, A [openloops] -| 2: 'nlo_eetta_p1_i2': e-, e+ => t, tbar, A, gl [openloops], [real] -| 3: 'nlo_eetta_p1_i3': e-, e+ => t, tbar, A [openloops], [virtual] -| 4: 'nlo_eetta_p1_i4': e-, e+ => t, tbar, A [inactive], [subtraction] -| ------------------------------------------------------------------------ -| Phase space: 6 channels, 5 dimensions -| Phase space: found 6 channels, collected in 3 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 6 channels, 8 dimensions -| Phase space: found 6 channels, collected in 3 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 6 channels, 5 dimensions -| Phase space: found 6 channels, collected in 3 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Applying user-defined cuts. -| Using user-defined general scale. -| Starting integration for process 'nlo_eetta_p1' part 'born' -| Integrate: iterations = 8:5000:"gw", 2:2000 -| Integrator: 3 chains, 6 channels, 5 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4998 1.2524400E+01 3.72E-01 2.97 2.10* 4.67 - 2 4996 1.2807824E+01 1.73E-01 1.35 0.95* 14.92 - 3 4994 1.2963919E+01 1.33E-01 1.03 0.73* 19.34 - 4 4990 1.2674928E+01 1.27E-01 1.00 0.71* 16.67 - 5 4986 1.2837989E+01 1.20E-01 0.94 0.66* 21.96 - 6 4984 1.2933354E+01 1.18E-01 0.91 0.64* 20.30 - 7 4980 1.2887357E+01 1.17E-01 0.91 0.64 20.70 - 8 4976 1.2790776E+01 1.22E-01 0.96 0.67 15.19 -|-----------------------------------------------------------------------------| - 8 39904 1.2840380E+01 4.77E-02 0.37 0.74 15.19 0.61 8 -|-----------------------------------------------------------------------------| - 9 1996 1.2797864E+01 1.93E-01 1.51 0.67* 18.81 - 10 1996 1.2800378E+01 2.07E-01 1.62 0.72 14.24 -|-----------------------------------------------------------------------------| - 10 3992 1.2799033E+01 1.41E-01 1.10 0.70 14.24 0.00 2 -|=============================================================================| -| Starting integration for process 'nlo_eetta_p1' part 'real' -| Integrate: iterations = 8:5000:"gw", 2:2000 -| Integrator: 3 chains, 6 channels, 8 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4998 -2.9091137E+00 1.08E-01 3.71 2.62* 4.00 - 2 4996 -3.0572272E+00 5.56E-02 1.82 1.29* 7.83 - 3 4994 -2.9804397E+00 4.96E-02 1.67 1.18* 8.52 - 4 4992 -2.9959728E+00 4.78E-02 1.60 1.13* 7.32 - 5 4990 -3.0098998E+00 4.32E-02 1.44 1.01* 11.05 - 6 4986 -2.9203492E+00 4.26E-02 1.46 1.03 10.67 - 7 4982 -3.0464332E+00 3.98E-02 1.31 0.92* 11.54 - 8 4980 -2.9558244E+00 4.12E-02 1.39 0.98 9.36 -|-----------------------------------------------------------------------------| - 8 39918 -2.9903692E+00 1.68E-02 0.56 1.12 9.36 1.09 8 -|-----------------------------------------------------------------------------| - 9 1998 -2.9814727E+00 5.67E-02 1.90 0.85* 17.42 - 10 1998 -3.0120295E+00 5.74E-02 1.91 0.85 14.85 -|-----------------------------------------------------------------------------| - 10 3996 -2.9965661E+00 4.04E-02 1.35 0.85 14.85 0.14 2 -|=============================================================================| -| Starting integration for process 'nlo_eetta_p1' part 'virtual' -| Integrate: iterations = 8:5000:"gw", 2:2000 -| Integrator: 3 chains, 6 channels, 5 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4998 3.8408890E+00 1.13E-01 2.95 2.09* 4.64 - 2 4996 3.9098741E+00 5.08E-02 1.30 0.92* 16.08 - 3 4994 3.9065596E+00 4.17E-02 1.07 0.75* 19.09 - 4 4990 3.9001464E+00 3.79E-02 0.97 0.69* 19.53 - 5 4988 3.8776541E+00 3.82E-02 0.99 0.70 14.41 - 6 4986 3.8941272E+00 3.85E-02 0.99 0.70 14.97 - 7 4982 3.8246674E+00 4.00E-02 1.05 0.74 15.07 - 8 4978 3.9276827E+00 3.63E-02 0.92 0.65* 17.79 -|-----------------------------------------------------------------------------| - 8 39912 3.8906010E+00 1.49E-02 0.38 0.77 17.79 0.63 8 -|-----------------------------------------------------------------------------| - 9 1998 3.8918324E+00 5.66E-02 1.45 0.65* 22.36 - 10 1998 3.9212412E+00 5.78E-02 1.47 0.66 18.98 -|-----------------------------------------------------------------------------| - 10 3996 3.9062326E+00 4.04E-02 1.04 0.65 18.98 0.13 2 -|=============================================================================| -| Integrate: sum of all components -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 0 1.3708699E+01 1.52E-01 1.11 0.00* 12.41 -| NLO Correction: [O(alpha_s+1)/O(alpha_s)] -| ( 7.1073 +- 0.45323 ) % -|=============================================================================| -| WHIZARD run finished. -|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettaa.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettaa.ref (revision 8291) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettaa.ref (revision 8292) @@ -1,180 +0,0 @@ -?openmp_logging = false -?vis_history = false -?integration_timer = false -openmp_num_threads = 1 -?debug_decay = false -?debug_process = false -?debug_verbose = false -?write_raw = false -| Including Sindarin from 'nlo_settings.sin' -| Switching to model 'SM', scheme 'GF_MW_MZ' -$blha_ew_scheme = "alpha_qed" -SM.mZ => 9.118800000000E+01 -SM.mW => 8.041900200000E+01 -SM.mH => 1.250000000000E+02 -SM.GF => 1.166390000000E-05 -SM.wZ => 0.000000000000E+00 -SM.wtop => 0.000000000000E+00 -SM.wW => 0.000000000000E+00 -SM.wH => 0.000000000000E+00 -SM.ms => 0.000000000000E+00 -SM.mc => 0.000000000000E+00 -SM.mb => 0.000000000000E+00 -SM.mtop => 1.732000000000E+02 -SM.me => 0.000000000000E+00 -SM.mmu => 0.000000000000E+00 -SM.mtau => 1.777000000000E+00 -SM.alphas => 1.180000000000E-01 -?alphas_is_fixed = false -?alphas_from_mz = true -?alphas_from_lambda_qcd = false -alphas_nf = 5 -alphas_order = 2 -[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) -$method = "openloops" -seed = 8131 -sqrts = 1.000000000000E+03 -jet_algorithm = 2 -jet_r = 5.000000000000E-01 -?use_vamp_equivalences = false -| End of included 'nlo_settings.sin' -alpha_power = 4 -alphas_power = 0 -| Process library 'nlo_eettaa_lib': recorded process 'nlo_eettaa_p1' -mult_call_real = 2.000000000000E+00 -| Integrate: current process library needs compilation -| Process library 'nlo_eettaa_lib': compiling ... -| Process library 'nlo_eettaa_lib': writing makefile -| Process library 'nlo_eettaa_lib': removing old files -| Process library 'nlo_eettaa_lib': writing driver -| Process library 'nlo_eettaa_lib': creating source code -| Process library 'nlo_eettaa_lib': compiling sources -| Process library 'nlo_eettaa_lib': linking -| Process library 'nlo_eettaa_lib': loading -| Process library 'nlo_eettaa_lib': ... success. -| Integrate: compilation done -| QCD alpha: using a running strong coupling -| RNG: Initializing TAO random-number generator -| RNG: Setting seed for random-number generator to 8131 -| Initializing integration for process nlo_eettaa_p1: -| Beam structure: [any particles] -| Beam data (collision): -| e- (mass = 0.0000000E+00 GeV) -| e+ (mass = 0.0000000E+00 GeV) -| sqrts = 1.000000000000E+03 GeV -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eettaa_p1.i1.phs' -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eettaa_p1.i3.phs' -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| ------------------------------------------------------------------------ -| Process [scattering]: 'nlo_eettaa_p1' -| Library name = 'nlo_eettaa_lib' -| Process index = 1 -| Process components: -| 1: 'nlo_eettaa_p1_i1': e-, e+ => t, tbar, A, A [openloops] -| 2: 'nlo_eettaa_p1_i2': e-, e+ => t, tbar, A, A, gl [openloops], [real] -| 3: 'nlo_eettaa_p1_i3': e-, e+ => t, tbar, A, A [openloops], [virtual] -| 4: 'nlo_eettaa_p1_i4': e-, e+ => t, tbar, A, A [inactive], [subtraction] -| ------------------------------------------------------------------------ -| Phase space: 30 channels, 8 dimensions -| Phase space: found 30 channels, collected in 5 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 30 channels, 11 dimensions -| Phase space: found 30 channels, collected in 5 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 30 channels, 8 dimensions -| Phase space: found 30 channels, collected in 5 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Applying user-defined cuts. -| Using user-defined general scale. -| Starting integration for process 'nlo_eettaa_p1' part 'born' -| Integrate: iterations = 4:5000:"gw", 2:5000 -| Integrator: 5 chains, 30 channels, 8 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4980 3.6991977E-01 2.83E-02 7.64 5.39* 2.18 - 2 4960 3.8294377E-01 1.38E-02 3.61 2.54* 4.85 - 3 4940 3.9224681E-01 1.08E-02 2.76 1.94* 7.14 - 4 4926 3.7558931E-01 9.17E-03 2.44 1.71* 8.26 -|-----------------------------------------------------------------------------| - 4 19806 3.8203117E-01 6.10E-03 1.60 2.25 8.26 0.52 4 -|-----------------------------------------------------------------------------| - 5 4926 3.9202838E-01 9.92E-03 2.53 1.78 6.67 - 6 4926 4.1000778E-01 1.14E-02 2.79 1.96 5.74 -|-----------------------------------------------------------------------------| - 6 9852 3.9976089E-01 7.49E-03 1.87 1.86 5.74 1.41 2 -|=============================================================================| -| Starting integration for process 'nlo_eettaa_p1' part 'real' -| Integrate: iterations = 4:10000:"gw", 2:10000 -| Integrator: 5 chains, 30 channels, 11 dimensions -| Integrator: 10000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 9990 -8.4275061E-02 5.75E-03 6.82 6.82* 1.32 - 2 9972 -9.1266413E-02 2.86E-03 3.13 3.13* 3.21 - 3 9956 -8.8755374E-02 2.29E-03 2.58 2.58* 4.32 - 4 9942 -8.8393569E-02 2.25E-03 2.55 2.54* 4.48 -|-----------------------------------------------------------------------------| - 4 39860 -8.8941246E-02 1.36E-03 1.53 3.06 4.48 0.46 4 -|-----------------------------------------------------------------------------| - 5 9942 -8.6688153E-02 2.30E-03 2.65 2.64 3.10 - 6 9942 -9.0306432E-02 2.01E-03 2.22 2.22* 2.96 -|-----------------------------------------------------------------------------| - 6 19884 -8.8738134E-02 1.51E-03 1.70 2.40 2.96 1.41 2 -|=============================================================================| -| Starting integration for process 'nlo_eettaa_p1' part 'virtual' -| Integrate: iterations = 4:5000:"gw", 2:5000 -| Integrator: 5 chains, 30 channels, 8 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4980 1.3854202E-01 1.28E-02 9.22 6.50* 1.92 - 2 4968 1.1563610E-01 4.27E-03 3.69 2.60* 4.66 - 3 4960 1.1642866E-01 3.21E-03 2.76 1.94* 7.07 - 4 4942 1.1739786E-01 3.01E-03 2.57 1.81* 7.48 -|-----------------------------------------------------------------------------| - 4 19850 1.1717079E-01 1.93E-03 1.65 2.32 7.48 1.00 4 -|-----------------------------------------------------------------------------| - 5 4942 1.1815723E-01 3.21E-03 2.72 1.91 5.88 - 6 4942 1.2353296E-01 3.72E-03 3.01 2.11 4.59 -|-----------------------------------------------------------------------------| - 6 9884 1.2045357E-01 2.43E-03 2.02 2.00 4.59 1.20 2 -|=============================================================================| -| Integrate: sum of all components -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 0 4.3147632E-01 8.02E-03 1.86 0.00* 4.50 -| NLO Correction: [O(alpha_s+1)/O(alpha_s)] -| ( 7.9336 +- 0.73080 ) % -|=============================================================================| -| WHIZARD run finished. -|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettah.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettah.ref (revision 8291) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettah.ref (revision 8292) @@ -1,184 +0,0 @@ -?openmp_logging = false -?vis_history = false -?integration_timer = false -openmp_num_threads = 1 -?debug_decay = false -?debug_process = false -?debug_verbose = false -?write_raw = false -| Including Sindarin from 'nlo_settings.sin' -| Switching to model 'SM', scheme 'GF_MW_MZ' -$blha_ew_scheme = "alpha_qed" -SM.mZ => 9.118800000000E+01 -SM.mW => 8.041900200000E+01 -SM.mH => 1.250000000000E+02 -SM.GF => 1.166390000000E-05 -SM.wZ => 0.000000000000E+00 -SM.wtop => 0.000000000000E+00 -SM.wW => 0.000000000000E+00 -SM.wH => 0.000000000000E+00 -SM.ms => 0.000000000000E+00 -SM.mc => 0.000000000000E+00 -SM.mb => 0.000000000000E+00 -SM.mtop => 1.732000000000E+02 -SM.me => 0.000000000000E+00 -SM.mmu => 0.000000000000E+00 -SM.mtau => 1.777000000000E+00 -SM.alphas => 1.180000000000E-01 -?alphas_is_fixed = false -?alphas_from_mz = true -?alphas_from_lambda_qcd = false -alphas_nf = 5 -alphas_order = 2 -[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) -$method = "openloops" -seed = 8131 -sqrts = 1.000000000000E+03 -jet_algorithm = 2 -jet_r = 5.000000000000E-01 -?use_vamp_equivalences = false -| End of included 'nlo_settings.sin' -alpha_power = 4 -alphas_power = 0 -[user variable] theta_0 = 7.000000000000E-01 -[user variable] tops = PDG(6, -6) -| Process library 'nlo_eettah_lib': recorded process 'nlo_eettah_p1' -| Integrate: current process library needs compilation -| Process library 'nlo_eettah_lib': compiling ... -| Process library 'nlo_eettah_lib': writing makefile -| Process library 'nlo_eettah_lib': removing old files -| Process library 'nlo_eettah_lib': writing driver -| Process library 'nlo_eettah_lib': creating source code -| Process library 'nlo_eettah_lib': compiling sources -| Process library 'nlo_eettah_lib': linking -| Process library 'nlo_eettah_lib': loading -| Process library 'nlo_eettah_lib': ... success. -| Integrate: compilation done -| QCD alpha: using a running strong coupling -| RNG: Initializing TAO random-number generator -| RNG: Setting seed for random-number generator to 8131 -| Initializing integration for process nlo_eettah_p1: -| Beam structure: [any particles] -| Beam data (collision): -| e- (mass = 0.0000000E+00 GeV) -| e+ (mass = 0.0000000E+00 GeV) -| sqrts = 1.000000000000E+03 GeV -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eettah_p1.i1.phs' -| Phase space: generating configuration ... -| Phase space: ... success. -| Phase space: writing configuration file 'nlo_eettah_p1.i3.phs' -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| One-Loop-Provider: Using OpenLoops -| Loading library: [...] -| ------------------------------------------------------------------------ -| Process [scattering]: 'nlo_eettah_p1' -| Library name = 'nlo_eettah_lib' -| Process index = 1 -| Process components: -| 1: 'nlo_eettah_p1_i1': e-, e+ => t, tbar, H, A [openloops] -| 2: 'nlo_eettah_p1_i2': e-, e+ => t, tbar, H, A, gl [openloops], [real] -| 3: 'nlo_eettah_p1_i3': e-, e+ => t, tbar, H, A [openloops], [virtual] -| 4: 'nlo_eettah_p1_i4': e-, e+ => t, tbar, H, A [inactive], [subtraction] -| ------------------------------------------------------------------------ -| Phase space: 12 channels, 8 dimensions -| Phase space: found 12 channels, collected in 2 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 12 channels, 11 dimensions -| Phase space: found 12 channels, collected in 2 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Phase space: 12 channels, 8 dimensions -| Phase space: found 12 channels, collected in 2 groves. -| Phase space: no equivalences between channels used. -| Phase space: wood -| Applying user-defined cuts. -| Using user-defined general scale. -| Starting integration for process 'nlo_eettah_p1' part 'born' -| Integrate: iterations = 4:5000:"gw", 3:3000 -| Integrator: 2 chains, 12 channels, 8 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4992 9.0935267E-02 4.11E-03 4.52 3.19* 2.86 - 2 4986 9.8656204E-02 1.86E-03 1.89 1.33* 10.18 - 3 4980 9.5953441E-02 1.57E-03 1.63 1.15* 13.24 - 4 4974 9.5958323E-02 1.40E-03 1.46 1.03* 15.01 -|-----------------------------------------------------------------------------| - 4 19932 9.6337533E-02 8.90E-04 0.92 1.30 15.01 1.14 4 -|-----------------------------------------------------------------------------| - 5 2994 9.8368415E-02 2.07E-03 2.10 1.15 12.07 - 6 2994 9.7963881E-02 1.99E-03 2.03 1.11* 11.02 - 7 2994 9.8061079E-02 1.91E-03 1.95 1.07* 10.49 -|-----------------------------------------------------------------------------| - 7 8982 9.8123651E-02 1.15E-03 1.17 1.11 10.49 0.01 3 -|=============================================================================| -| Starting integration for process 'nlo_eettah_p1' part 'real' -| Integrate: iterations = 4:5000:"gw", 3:3000 -| Integrator: 2 chains, 12 channels, 11 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4992 -1.8491776E-02 1.02E-03 5.52 3.90* 2.34 - 2 4986 -1.8754982E-02 5.53E-04 2.95 2.08* 5.16 - 3 4980 -1.8382428E-02 3.69E-04 2.01 1.42* 9.84 - 4 4974 -1.9745867E-02 4.06E-04 2.06 1.45 7.51 -|-----------------------------------------------------------------------------| - 4 19932 -1.8926813E-02 2.38E-04 1.26 1.78 7.51 2.18 4 -|-----------------------------------------------------------------------------| - 5 2994 -1.9430246E-02 4.98E-04 2.56 1.40* 10.21 - 6 2994 -1.8719668E-02 5.16E-04 2.76 1.51 7.28 - 7 2994 -1.8622458E-02 4.47E-04 2.40 1.31* 7.34 -|-----------------------------------------------------------------------------| - 7 8982 -1.8905631E-02 2.79E-04 1.48 1.40 7.34 0.82 3 -|=============================================================================| -| Starting integration for process 'nlo_eettah_p1' part 'virtual' -| Integrate: iterations = 4:5000:"gw", 3:3000 -| Integrator: 2 chains, 12 channels, 8 dimensions -| Integrator: 5000 initial calls, 20 bins, stratified = T -| Integrator: VAMP -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 4992 1.5821915E-02 6.42E-04 4.06 2.87* 3.55 - 2 4986 1.6399422E-02 3.34E-04 2.04 1.44* 9.49 - 3 4980 1.6785167E-02 2.79E-04 1.66 1.17* 12.24 - 4 4974 1.6613581E-02 2.49E-04 1.50 1.06* 14.83 -|-----------------------------------------------------------------------------| - 4 19932 1.6572963E-02 1.57E-04 0.95 1.34 14.83 0.75 4 -|-----------------------------------------------------------------------------| - 5 2994 1.6631895E-02 3.13E-04 1.88 1.03* 15.37 - 6 2994 1.6572555E-02 3.24E-04 1.96 1.07 13.92 - 7 2994 1.7279491E-02 3.24E-04 1.87 1.03* 13.72 -|-----------------------------------------------------------------------------| - 7 8982 1.6823922E-02 1.85E-04 1.10 1.04 13.72 1.48 3 -|=============================================================================| -| Integrate: sum of all components -|=============================================================================| -| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | -|=============================================================================| - 1 0 9.6041941E-02 1.19E-03 1.24 0.00* 9.08 -| NLO Correction: [O(alpha_s+1)/O(alpha_s)] -| ( -2.1215 +- 0.34244 ) % -|=============================================================================| -| WHIZARD run finished. -|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4tj.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4tj.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4tj.ref (revision 8292) @@ -0,0 +1,161 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 2 +alphas_power = 3 +| Process library 'nlo_ee4tj_lib': recorded process 'nlo_ee4tj_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_ee4tj_lib': compiling ... +| Process library 'nlo_ee4tj_lib': writing makefile +| Process library 'nlo_ee4tj_lib': removing old files +| Process library 'nlo_ee4tj_lib': writing driver +| Process library 'nlo_ee4tj_lib': creating source code +| Process library 'nlo_ee4tj_lib': compiling sources +| Process library 'nlo_ee4tj_lib': linking +| Process library 'nlo_ee4tj_lib': loading +| Process library 'nlo_ee4tj_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_ee4tj_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +| Phase space: ... failed. Increasing phs_off_shell ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_ee4tj_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... failed. Increasing phs_off_shell ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_ee4tj_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_ee4tj_p1' +| Library name = 'nlo_ee4tj_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_ee4tj_p1_i1': e-, e+ => t, tbar, t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_ee4tj_p1_i2': e-, e+ => t, tbar, t, tbar, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_ee4tj_p1_i3': e-, e+ => t, tbar, t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_ee4tj_p1_i4': e-, e+ => t, tbar, t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 64 channels, 11 dimensions +| Phase space: found 64 channels, collected in 2 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 64 channels, 14 dimensions +| Phase space: found 64 channels, collected in 2 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 64 channels, 11 dimensions +| Phase space: found 64 channels, collected in 2 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_ee4tj_p1' part 'born' +| Integrate: iterations = 1:100:"gw" +| Integrator: 2 chains, 64 channels, 11 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| +Warning: VAMP: too few calls, resetting n_calls to 640 + 1 640 2.6816175E-05 5.43E-06 20.26 5.13* 11.49 +|-----------------------------------------------------------------------------| + 1 640 2.6816175E-05 5.43E-06 20.26 5.13 11.49 +|=============================================================================| +| Starting integration for process 'nlo_ee4tj_p1' part 'real' +| Integrate: iterations = 1:100:"gw" +| Integrator: 2 chains, 64 channels, 14 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| +Warning: VAMP: too few calls, resetting n_calls to 640 + 1 640 2.5351061E-06 2.84E-06 112.13 28.37* 11.72 +|-----------------------------------------------------------------------------| + 1 640 2.5351061E-06 2.84E-06 112.13 28.37 11.72 +|=============================================================================| +| Starting integration for process 'nlo_ee4tj_p1' part 'virtual' +| Integrate: iterations = 1:100:"gw" +| Integrator: 2 chains, 64 channels, 11 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| +Warning: VAMP: too few calls, resetting n_calls to 640 + 1 640 2.8620247E-05 5.66E-06 19.78 5.00* 12.12 +|-----------------------------------------------------------------------------| + 1 640 2.8620247E-05 5.66E-06 19.78 5.00 12.12 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 5.7971527E-05 8.35E-06 14.40 0.00* 11.80 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| (116.1812 +- 33.35126 ) % +|=============================================================================| +| There were no errors and 3 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettzjj.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettzjj.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettzjj.ref (revision 8292) @@ -0,0 +1,158 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 3 +alphas_power = 2 +| Process library 'nlo_eettzjj_lib': recorded process 'nlo_eettzjj_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_eettzjj_lib': compiling ... +| Process library 'nlo_eettzjj_lib': writing makefile +| Process library 'nlo_eettzjj_lib': removing old files +| Process library 'nlo_eettzjj_lib': writing driver +| Process library 'nlo_eettzjj_lib': creating source code +| Process library 'nlo_eettzjj_lib': compiling sources +| Process library 'nlo_eettzjj_lib': linking +| Process library 'nlo_eettzjj_lib': loading +| Process library 'nlo_eettzjj_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_eettzjj_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eettzjj_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eettzjj_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_eettzjj_p1' +| Library name = 'nlo_eettzjj_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_eettzjj_p1_i1': e-, e+ => t, tbar, Z, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_eettzjj_p1_i2': e-, e+ => t, tbar, Z, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_eettzjj_p1_i3': e-, e+ => t, tbar, Z, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_eettzjj_p1_i4': e-, e+ => t, tbar, Z, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 94 channels, 11 dimensions +| Phase space: found 94 channels, collected in 7 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 94 channels, 14 dimensions +| Phase space: found 94 channels, collected in 7 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 94 channels, 11 dimensions +| Phase space: found 94 channels, collected in 7 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_eettzjj_p1' part 'born' +| Integrate: iterations = 1:940:"gw" +| Integrator: 7 chains, 94 channels, 11 dimensions +| Integrator: 940 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 940 5.7893209E-02 1.66E-02 28.61 8.77* 10.18 +|-----------------------------------------------------------------------------| + 1 940 5.7893209E-02 1.66E-02 28.61 8.77 10.18 +|=============================================================================| +| Starting integration for process 'nlo_eettzjj_p1' part 'real' +| Integrate: iterations = 1:940:"gw" +| Integrator: 7 chains, 94 channels, 14 dimensions +| Integrator: 940 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 940 -2.2707819E-03 1.02E-02 447.80 137.29* 10.21 +|-----------------------------------------------------------------------------| + 1 940 -2.2707819E-03 1.02E-02 447.80 137.29 10.21 +|=============================================================================| +| Starting integration for process 'nlo_eettzjj_p1' part 'virtual' +| Integrate: iterations = 1:940:"gw" +| Integrator: 7 chains, 94 channels, 11 dimensions +| Integrator: 940 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 940 3.6978740E-03 2.21E-03 59.74 18.32* 10.46 +|-----------------------------------------------------------------------------| + 1 940 3.6978740E-03 2.21E-03 59.74 18.32 10.46 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 5.9320301E-02 1.96E-02 32.97 0.00* 9.82 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| ( 2.4650 +- 17.98784 ) % +|=============================================================================| +| There were no errors and 2 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettjj.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettjj.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettjj.ref (revision 8292) @@ -0,0 +1,158 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 2 +alphas_power = 2 +| Process library 'nlo_eettjj_lib': recorded process 'nlo_eettjj_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_eettjj_lib': compiling ... +| Process library 'nlo_eettjj_lib': writing makefile +| Process library 'nlo_eettjj_lib': removing old files +| Process library 'nlo_eettjj_lib': writing driver +| Process library 'nlo_eettjj_lib': creating source code +| Process library 'nlo_eettjj_lib': compiling sources +| Process library 'nlo_eettjj_lib': linking +| Process library 'nlo_eettjj_lib': loading +| Process library 'nlo_eettjj_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_eettjj_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eettjj_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eettjj_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_eettjj_p1' +| Library name = 'nlo_eettjj_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_eettjj_p1_i1': e-, e+ => t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_eettjj_p1_i2': e-, e+ => t, tbar, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_eettjj_p1_i3': e-, e+ => t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_eettjj_p1_i4': e-, e+ => t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 18 channels, 8 dimensions +| Phase space: found 18 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 18 channels, 11 dimensions +| Phase space: found 18 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 18 channels, 8 dimensions +| Phase space: found 18 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_eettjj_p1' part 'born' +| Integrate: iterations = 1:180:"gw" +| Integrator: 5 chains, 18 channels, 8 dimensions +| Integrator: 180 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 162 6.8729844E+00 2.28E+00 33.24 4.23* 14.23 +|-----------------------------------------------------------------------------| + 1 162 6.8729844E+00 2.28E+00 33.24 4.23 14.23 +|=============================================================================| +| Starting integration for process 'nlo_eettjj_p1' part 'real' +| Integrate: iterations = 1:180:"gw" +| Integrator: 5 chains, 18 channels, 11 dimensions +| Integrator: 180 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 162 -2.0847511E+01 1.50E+01 71.84 9.14* 12.51 +|-----------------------------------------------------------------------------| + 1 162 -2.0847511E+01 1.50E+01 71.84 9.14 12.51 +|=============================================================================| +| Starting integration for process 'nlo_eettjj_p1' part 'virtual' +| Integrate: iterations = 1:180:"gw" +| Integrator: 5 chains, 18 channels, 8 dimensions +| Integrator: 180 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 162 6.3956075E-01 3.42E-01 53.47 6.81* 12.36 +|-----------------------------------------------------------------------------| + 1 162 6.3956075E-01 3.42E-01 53.47 6.81 12.36 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 -1.3334966E+01 1.52E+01 113.64 0.00* -24.94 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| (******** +- 238.87201 ) % +|=============================================================================| +| There were no errors and 2 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4j.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4j.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4j.ref (revision 8292) @@ -0,0 +1,161 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 2 +alphas_power = 2 +$restrictions = "!W+:W-" +| Process library 'nlo_ee4j_lib': recorded process 'nlo_ee4j_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_ee4j_lib': compiling ... +| Process library 'nlo_ee4j_lib': writing makefile +| Process library 'nlo_ee4j_lib': removing old files +| Process library 'nlo_ee4j_lib': writing driver +| Process library 'nlo_ee4j_lib': creating source code +| Process library 'nlo_ee4j_lib': compiling sources +| Process library 'nlo_ee4j_lib': linking +| Process library 'nlo_ee4j_lib': loading +| Process library 'nlo_ee4j_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_ee4j_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +Warning: Intermediate decay of zero-width particle W- may be possible. +Warning: Intermediate decay of zero-width particle W+ may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_ee4j_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_ee4j_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_ee4j_p1' +| Library name = 'nlo_ee4j_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_ee4j_p1_i1': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_ee4j_p1_i2': e-, e+ => d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_ee4j_p1_i3': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_ee4j_p1_i4': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 48 channels, 8 dimensions +| Phase space: found 48 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 48 channels, 11 dimensions +| Phase space: found 48 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 48 channels, 8 dimensions +| Phase space: found 48 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_ee4j_p1' part 'born' +| Integrate: iterations = 1:480:"gw" +| Integrator: 5 chains, 48 channels, 8 dimensions +| Integrator: 480 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 432 1.7396878E+02 6.71E+01 38.59 8.02* 13.29 +|-----------------------------------------------------------------------------| + 1 432 1.7396878E+02 6.71E+01 38.59 8.02 13.29 +|=============================================================================| +| Starting integration for process 'nlo_ee4j_p1' part 'real' +| Integrate: iterations = 1:480:"gw" +| Integrator: 5 chains, 48 channels, 11 dimensions +| Integrator: 480 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 432 -1.2064627E+01 5.28E+01 437.91 91.02* 12.10 +|-----------------------------------------------------------------------------| + 1 432 -1.2064627E+01 5.28E+01 437.91 91.02 12.10 +|=============================================================================| +| Starting integration for process 'nlo_ee4j_p1' part 'virtual' +| Integrate: iterations = 1:480:"gw" +| Integrator: 5 chains, 48 channels, 8 dimensions +| Integrator: 480 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 432 -9.1758178E+01 5.42E+01 59.06 12.28* 11.47 +|-----------------------------------------------------------------------------| + 1 432 -9.1758178E+01 5.42E+01 59.06 12.28 11.47 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 7.0145979E+01 1.01E+02 144.22 0.00* 5.36 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| (-59.6790 +- 49.22474 ) % +|=============================================================================| +| There were no errors and 4 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eejjj.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eejjj.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eejjj.ref (revision 8292) @@ -0,0 +1,158 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 2 +alphas_power = 1 +| Process library 'nlo_eejjj_lib': recorded process 'nlo_eejjj_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_eejjj_lib': compiling ... +| Process library 'nlo_eejjj_lib': writing makefile +| Process library 'nlo_eejjj_lib': removing old files +| Process library 'nlo_eejjj_lib': writing driver +| Process library 'nlo_eejjj_lib': creating source code +| Process library 'nlo_eejjj_lib': compiling sources +| Process library 'nlo_eejjj_lib': linking +| Process library 'nlo_eejjj_lib': loading +| Process library 'nlo_eejjj_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_eejjj_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eejjj_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eejjj_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_eejjj_p1' +| Library name = 'nlo_eejjj_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_eejjj_p1_i1': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_eejjj_p1_i2': e-, e+ => d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_eejjj_p1_i3': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_eejjj_p1_i4': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 2 channels, 5 dimensions +| Phase space: found 2 channels, collected in 1 grove. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 2 channels, 8 dimensions +| Phase space: found 2 channels, collected in 1 grove. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 2 channels, 5 dimensions +| Phase space: found 2 channels, collected in 1 grove. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_eejjj_p1' part 'born' +| Integrate: iterations = 1:100:"gw" +| Integrator: 1 chains, 2 channels, 5 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 100 3.0454213E+02 1.18E+02 38.83 3.88* 4.87 +|-----------------------------------------------------------------------------| + 1 100 3.0454213E+02 1.18E+02 38.83 3.88 4.87 +|=============================================================================| +| Starting integration for process 'nlo_eejjj_p1' part 'real' +| Integrate: iterations = 1:100:"gw" +| Integrator: 1 chains, 2 channels, 8 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 100 3.1185385E+02 2.40E+02 76.94 7.69* 4.33 +|-----------------------------------------------------------------------------| + 1 100 3.1185385E+02 2.40E+02 76.94 7.69 4.33 +|=============================================================================| +| Starting integration for process 'nlo_eejjj_p1' part 'virtual' +| Integrate: iterations = 1:100:"gw" +| Integrator: 1 chains, 2 channels, 5 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 100 -1.8695383E+01 1.37E+01 73.41 7.34* 2.76 +|-----------------------------------------------------------------------------| + 1 100 -1.8695383E+01 1.37E+01 73.41 7.34 2.76 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 5.9770060E+02 2.68E+02 44.82 0.00* 4.44 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| ( 96.2620 +- 87.32405 ) % +|=============================================================================| +| There were no errors and 2 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee5j.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee5j.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee5j.ref (revision 8292) @@ -0,0 +1,161 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 2 +alphas_power = 3 +$restrictions = "!W+:W-" +| Process library 'nlo_ee5j_lib': recorded process 'nlo_ee5j_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_ee5j_lib': compiling ... +| Process library 'nlo_ee5j_lib': writing makefile +| Process library 'nlo_ee5j_lib': removing old files +| Process library 'nlo_ee5j_lib': writing driver +| Process library 'nlo_ee5j_lib': creating source code +| Process library 'nlo_ee5j_lib': compiling sources +| Process library 'nlo_ee5j_lib': linking +| Process library 'nlo_ee5j_lib': loading +| Process library 'nlo_ee5j_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_ee5j_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +Warning: Intermediate decay of zero-width particle W- may be possible. +Warning: Intermediate decay of zero-width particle W+ may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_ee5j_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_ee5j_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_ee5j_p1' +| Library name = 'nlo_ee5j_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_ee5j_p1_i1': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_ee5j_p1_i2': e-, e+ => d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_ee5j_p1_i3': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_ee5j_p1_i4': e-, e+ => u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 104 channels, 11 dimensions +| Phase space: found 104 channels, collected in 3 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 104 channels, 14 dimensions +| Phase space: found 104 channels, collected in 3 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 104 channels, 11 dimensions +| Phase space: found 104 channels, collected in 3 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_ee5j_p1' part 'born' +| Integrate: iterations = 1:1040:"gw" +| Integrator: 3 chains, 104 channels, 11 dimensions +| Integrator: 1040 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 1040 1.4036172E+01 6.23E+00 44.42 14.32* 10.76 +|-----------------------------------------------------------------------------| + 1 1040 1.4036172E+01 6.23E+00 44.42 14.32 10.76 +|=============================================================================| +| Starting integration for process 'nlo_ee5j_p1' part 'real' +| Integrate: iterations = 1:1040:"gw" +| Integrator: 3 chains, 104 channels, 14 dimensions +| Integrator: 1040 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 1040 3.6383234E+00 2.63E+00 72.27 23.31* 11.65 +|-----------------------------------------------------------------------------| + 1 1040 3.6383234E+00 2.63E+00 72.27 23.31 11.65 +|=============================================================================| +| Starting integration for process 'nlo_ee5j_p1' part 'virtual' +| Integrate: iterations = 1:1040:"gw" +| Integrator: 3 chains, 104 channels, 11 dimensions +| Integrator: 1040 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 1040 -2.6105344E+01 1.45E+01 55.56 17.92* 10.04 +|-----------------------------------------------------------------------------| + 1 1040 -2.6105344E+01 1.45E+01 55.56 17.92 10.04 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 -8.4308487E+00 1.60E+01 189.84 0.00* -5.21 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| (******** +- 126.82245 ) % +|=============================================================================| +| There were no errors and 4 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eetthjj.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eetthjj.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eetthjj.ref (revision 8292) @@ -0,0 +1,161 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 3 +alphas_power = 2 +| Process library 'nlo_eetthjj_lib': recorded process 'nlo_eetthjj_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_eetthjj_lib': compiling ... +| Process library 'nlo_eetthjj_lib': writing makefile +| Process library 'nlo_eetthjj_lib': removing old files +| Process library 'nlo_eetthjj_lib': writing driver +| Process library 'nlo_eetthjj_lib': creating source code +| Process library 'nlo_eetthjj_lib': compiling sources +| Process library 'nlo_eetthjj_lib': linking +| Process library 'nlo_eetthjj_lib': loading +| Process library 'nlo_eetthjj_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_eetthjj_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eetthjj_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eetthjj_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_eetthjj_p1' +| Library name = 'nlo_eetthjj_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_eetthjj_p1_i1': e-, e+ => t, tbar, H, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_eetthjj_p1_i2': e-, e+ => t, tbar, H, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_eetthjj_p1_i3': e-, e+ => t, tbar, H, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_eetthjj_p1_i4': e-, e+ => t, tbar, H, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 38 channels, 11 dimensions +| Phase space: found 38 channels, collected in 2 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 38 channels, 14 dimensions +| Phase space: found 38 channels, collected in 2 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 38 channels, 11 dimensions +| Phase space: found 38 channels, collected in 2 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_eetthjj_p1' part 'born' +| Integrate: iterations = 1:100:"gw" +| Integrator: 2 chains, 38 channels, 11 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| +Warning: VAMP: too few calls, resetting n_calls to 380 + 1 380 3.2695814E-02 1.27E-02 38.75 7.55* 10.18 +|-----------------------------------------------------------------------------| + 1 380 3.2695814E-02 1.27E-02 38.75 7.55 10.18 +|=============================================================================| +| Starting integration for process 'nlo_eetthjj_p1' part 'real' +| Integrate: iterations = 1:100:"gw" +| Integrator: 2 chains, 38 channels, 14 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| +Warning: VAMP: too few calls, resetting n_calls to 380 + 1 380 -3.5561979E-03 3.40E-03 95.52 18.62* 10.49 +|-----------------------------------------------------------------------------| + 1 380 -3.5561979E-03 3.40E-03 95.52 18.62 10.49 +|=============================================================================| +| Starting integration for process 'nlo_eetthjj_p1' part 'virtual' +| Integrate: iterations = 1:100:"gw" +| Integrator: 2 chains, 38 channels, 11 dimensions +| Integrator: 100 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| +Warning: VAMP: too few calls, resetting n_calls to 380 + 1 380 -3.0684326E-03 3.83E-03 124.67 24.30* 10.00 +|-----------------------------------------------------------------------------| + 1 380 -3.0684326E-03 3.83E-03 124.67 24.30 10.00 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 2.6071183E-02 1.37E-02 52.41 0.00* 8.12 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| (-20.2614 +- 17.50661 ) % +|=============================================================================| +| There were no errors and 5 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettjjj.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettjjj.ref (revision 0) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_eettjjj.ref (revision 8292) @@ -0,0 +1,158 @@ +?openmp_logging = false +?vis_history = false +?integration_timer = false +openmp_num_threads = 1 +?debug_decay = false +?debug_process = false +?debug_verbose = false +?write_raw = false +| Switching to model 'SM', scheme 'GF_MW_MZ' +$blha_ew_scheme = "alpha_qed" +SM.mZ => 9.118800000000E+01 +SM.mW => 8.041900200000E+01 +SM.mH => 1.250000000000E+02 +SM.GF => 1.166390000000E-05 +SM.wZ => 0.000000000000E+00 +SM.wtop => 0.000000000000E+00 +SM.wW => 0.000000000000E+00 +SM.wH => 0.000000000000E+00 +SM.ms => 0.000000000000E+00 +SM.mc => 0.000000000000E+00 +SM.mb => 0.000000000000E+00 +SM.mtop => 1.732000000000E+02 +SM.me => 0.000000000000E+00 +SM.mmu => 0.000000000000E+00 +SM.mtau => 1.777000000000E+00 +SM.alphas => 1.180000000000E-01 +?alphas_is_fixed = false +?alphas_from_mz = true +?alphas_from_lambda_qcd = false +alphas_nf = 5 +alphas_order = 2 +[user variable] jet = PDG(2, -2, 1, -1, 3, -3, 4, -4, 5, -5, 21) +$exclude_gauge_splittings = "t" +$method = "openloops" +seed = 8131 +sqrts = 1.000000000000E+03 +jet_algorithm = 2 +jet_r = 5.000000000000E-01 +?use_vamp_equivalences = false +| End of included 'nlo_settings.sin' +alpha_power = 2 +alphas_power = 3 +| Process library 'nlo_eettjjj_lib': recorded process 'nlo_eettjjj_p1' +| Integrate: current process library needs compilation +| Process library 'nlo_eettjjj_lib': compiling ... +| Process library 'nlo_eettjjj_lib': writing makefile +| Process library 'nlo_eettjjj_lib': removing old files +| Process library 'nlo_eettjjj_lib': writing driver +| Process library 'nlo_eettjjj_lib': creating source code +| Process library 'nlo_eettjjj_lib': compiling sources +| Process library 'nlo_eettjjj_lib': linking +| Process library 'nlo_eettjjj_lib': loading +| Process library 'nlo_eettjjj_lib': ... success. +| Integrate: compilation done +| QCD alpha: using a running strong coupling +| RNG: Initializing TAO random-number generator +| RNG: Setting seed for random-number generator to 8131 +| Initializing integration for process nlo_eettjjj_p1: +| Beam structure: [any particles] +| Beam data (collision): +| e- (mass = 0.0000000E+00 GeV) +| e+ (mass = 0.0000000E+00 GeV) +| sqrts = 1.000000000000E+03 GeV +| Phase space: generating configuration ... +Warning: Intermediate decay of zero-width particle Z may be possible. +Warning: Intermediate decay of zero-width particle H may be possible. +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eettjjj_p1.i1.phs' +| Phase space: generating configuration ... +| Phase space: ... success. +| Phase space: writing configuration file 'nlo_eettjjj_p1.i3.phs' +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| One-Loop-Provider: Using OpenLoops +| Loading library: [...] +| ------------------------------------------------------------------------ +| Process [scattering]: 'nlo_eettjjj_p1' +| Library name = 'nlo_eettjjj_lib' +| Process index = 1 +| Process components: +| 1: 'nlo_eettjjj_p1_i1': e-, e+ => t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops] +| 2: 'nlo_eettjjj_p1_i2': e-, e+ => t, tbar, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl, d:dbar:u:ubar:s:sbar:c:cbar:b:bbar:gl [openloops], [real] +| 3: 'nlo_eettjjj_p1_i3': e-, e+ => t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [openloops], [virtual] +| 4: 'nlo_eettjjj_p1_i4': e-, e+ => t, tbar, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl, u:ubar:d:dbar:s:sbar:c:cbar:b:bbar:gl [inactive], [subtraction] +| ------------------------------------------------------------------------ +| Phase space: 78 channels, 11 dimensions +| Phase space: found 78 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 78 channels, 14 dimensions +| Phase space: found 78 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Phase space: 78 channels, 11 dimensions +| Phase space: found 78 channels, collected in 5 groves. +| Phase space: no equivalences between channels used. +| Phase space: wood +| Applying user-defined cuts. +| Using user-defined general scale. +| Starting integration for process 'nlo_eettjjj_p1' part 'born' +| Integrate: iterations = 1:780:"gw" +| Integrator: 5 chains, 78 channels, 11 dimensions +| Integrator: 780 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 780 5.1662429E-01 3.17E-01 61.38 17.14* 10.25 +|-----------------------------------------------------------------------------| + 1 780 5.1662429E-01 3.17E-01 61.38 17.14 10.25 +|=============================================================================| +| Starting integration for process 'nlo_eettjjj_p1' part 'real' +| Integrate: iterations = 1:780:"gw" +| Integrator: 5 chains, 78 channels, 14 dimensions +| Integrator: 780 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 780 1.6439844E-02 3.83E-01 2328.01 650.18* 10.35 +|-----------------------------------------------------------------------------| + 1 780 1.6439844E-02 3.83E-01 2328.01 650.18 10.35 +|=============================================================================| +| Starting integration for process 'nlo_eettjjj_p1' part 'virtual' +| Integrate: iterations = 1:780:"gw" +| Integrator: 5 chains, 78 channels, 11 dimensions +| Integrator: 780 initial calls, 20 bins, stratified = T +| Integrator: VAMP +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 780 9.5871775E-02 4.21E-02 43.94 12.27* 10.04 +|-----------------------------------------------------------------------------| + 1 780 9.5871775E-02 4.21E-02 43.94 12.27 10.04 +|=============================================================================| +| Integrate: sum of all components +|=============================================================================| +| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | +|=============================================================================| + 1 0 6.2893591E-01 4.99E-01 79.31 0.00* 10.22 +| NLO Correction: [O(alpha_s+1)/O(alpha_s)] +| ( 21.7395 +- 75.71365 ) % +|=============================================================================| +| There were no errors and 2 warning(s). +| WHIZARD run finished. +|=============================================================================| Index: trunk/share/tests/ext_tests_nlo/nlo_eettzjj.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_eettzjj.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_eettzjj.sin (revision 8292) @@ -1,30 +1,30 @@ include ("nlo_settings.sin") alpha_power = 3 alphas_power = 2 cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 2 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 2 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + sqrt (eval Pt [extract index 1 [Z]]**2 + eval M [extract index 1 [Z]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]]) / 2 elsif njet == 3 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + sqrt (eval Pt [extract index 1 [Z]]**2 + eval M [extract index 1 [Z]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]]) / 2 endif #?combined_nlo_integration = true process nlo_eettzjj_p1 = e1, E1 => t, T, Z, jet, jet {nlo_calculation = full } -integrate (nlo_eettzjj_p1) { iterations = 1:100:"gw" } +integrate (nlo_eettzjj_p1) { iterations = 1:940:"gw" } Index: trunk/share/tests/ext_tests_nlo/nlo_ee4j.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_ee4j.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_ee4j.sin (revision 8292) @@ -1,28 +1,28 @@ include("nlo_settings.sin") cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 4 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 4 then ( eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]] + eval Pt [extract index 4 [@clustered_jets]]) / 2 elsif njet == 5 then ( eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]] + eval Pt [extract index 3 [@clustered_jets]] + eval Pt [extract index 4 [@clustered_jets]] + eval Pt [extract index 5 [@clustered_jets]]) / 2 endif alpha_power = 2 alphas_power = 2 #?combined_nlo_integration = true process nlo_ee4j_p1 = e1, E1 => jet, jet, jet, jet { nlo_calculation = full $restrictions = "!W+:W-" } -integrate (nlo_ee4j_p1) { iterations = 1:100:"gw" } +integrate (nlo_ee4j_p1) { iterations = 1:480:"gw" } Index: trunk/share/tests/ext_tests_nlo/nlo_ee4tj.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_ee4tj.sin (revision 8291) +++ trunk/share/tests/ext_tests_nlo/nlo_ee4tj.sin (revision 8292) @@ -1,30 +1,30 @@ include ("nlo_settings.sin") alpha_power = 2 alphas_power = 3 cuts = let subevt @clustered_jets = cluster [jet] in let subevt @pt_selected = select if Pt > 30 GeV [@clustered_jets] in let subevt @eta_selected = select if abs(Eta) < 4 [@pt_selected] in count [@eta_selected] >= 1 scale = let subevt @clustered_jets = cluster [jet] in let int njet = count [@clustered_jets] in if njet == 1 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 2 [t]]**2 + eval M [extract index 2 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + sqrt (eval Pt [extract index 2 [T]]**2 + eval M [extract index 2 [T]]**2) + eval Pt [extract index 1 [@clustered_jets]]) / 2 elsif njet == 2 then (sqrt (eval Pt [extract index 1 [t]]**2 + eval M [extract index 1 [t]]**2) + sqrt (eval Pt [extract index 2 [t]]**2 + eval M [extract index 2 [t]]**2) + sqrt (eval Pt [extract index 1 [T]]**2 + eval M [extract index 1 [T]]**2) + sqrt (eval Pt [extract index 2 [T]]**2 + eval M [extract index 2 [T]]**2) + eval Pt [extract index 1 [@clustered_jets]] + eval Pt [extract index 2 [@clustered_jets]]) / 2 endif #?combined_nlo_integration = true process nlo_ee4tj_p1 = e1, E1 => t, T, t, T, jet { nlo_calculation = full } -integrate (nlo_ee4tj_p1) { iterations = 1:100:"gw" } +integrate (nlo_ee4tj_p1) { iterations = 1:640:"gw" }