Index: trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4j.ref =================================================================== --- trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4j.ref (revision 8856) +++ trunk/share/tests/ext_tests_nlo/ref-output/nlo_ee4j.ref (revision 8857) @@ -1,177 +1,178 @@ ?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 $integration_method = "vamp2" $rng_method = "rng_stream" | End of included 'nlo_settings.sin' +$openloops_allowed_libs = "pplljj" 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 RNG Stream 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: Using 320 equivalences between channels. | Phase space: wood | Phase space: 48 channels, 11 dimensions | Phase space: found 48 channels, collected in 5 groves. | Phase space: Using 320 equivalences between channels. | Phase space: wood | Phase space: 48 channels, 8 dimensions | Phase space: found 48 channels, collected in 5 groves. | Phase space: Using 320 equivalences between channels. | 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: Using VAMP2 channel equivalences | Integrator: Write grid header and grids to 'nlo_ee4j_p1.m1.vg2' | Integrator: Grid checkpoint after each iteration | Integrator: 480 initial calls, 20 max. bins, stratified = T | Integrator: VAMP2 |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| | VAMP2: Initialize new grids and write to file 'nlo_ee4j_p1.m1.vg2'. | VAMP2: set chain: use chained weights. 1 494 8.8893730E+01 4.52E+01 50.81 11.29* 12.59 |-----------------------------------------------------------------------------| 1 494 8.8893730E+01 4.52E+01 50.81 11.29 12.59 |=============================================================================| | Starting integration for process 'nlo_ee4j_p1' part 'real' | Integrate: iterations = 1:480:"gw" | Integrator: 5 chains, 48 channels, 11 dimensions | Integrator: Using VAMP2 channel equivalences | Integrator: Write grid header and grids to 'nlo_ee4j_p1.m2.vg2' | Integrator: Grid checkpoint after each iteration | Integrator: 480 initial calls, 20 max. bins, stratified = T | Integrator: VAMP2 |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| | VAMP2: Initialize new grids and write to file 'nlo_ee4j_p1.m2.vg2'. | VAMP2: set chain: use chained weights. 1 494 4.2862344E+01 1.13E+01 26.34 5.85* 14.36 |-----------------------------------------------------------------------------| 1 494 4.2862344E+01 1.13E+01 26.34 5.85 14.36 |=============================================================================| | Starting integration for process 'nlo_ee4j_p1' part 'virtual' | Integrate: iterations = 1:480:"gw" | Integrator: 5 chains, 48 channels, 8 dimensions | Integrator: Using VAMP2 channel equivalences | Integrator: Write grid header and grids to 'nlo_ee4j_p1.m3.vg2' | Integrator: Grid checkpoint after each iteration | Integrator: 480 initial calls, 20 max. bins, stratified = T | Integrator: VAMP2 |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| | VAMP2: Initialize new grids and write to file 'nlo_ee4j_p1.m3.vg2'. | VAMP2: set chain: use chained weights. 1 488 -2.3587867E+01 9.36E+00 39.69 8.77* 10.96 |-----------------------------------------------------------------------------| 1 488 -2.3587867E+01 9.36E+00 39.69 8.77 10.96 |=============================================================================| | Integrate: sum of all components |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 0 1.0816821E+02 4.75E+01 43.91 0.00* 10.77 | NLO Correction: [O(alpha_s+1)/O(alpha_s)] | ( 21.68 +- 19.84 ) % |=============================================================================| | There were no errors and 4 warning(s). | WHIZARD run finished. |=============================================================================| Index: trunk/share/tests/ext_tests_nlo/nlo_ee4j.sin =================================================================== --- trunk/share/tests/ext_tests_nlo/nlo_ee4j.sin (revision 8856) +++ trunk/share/tests/ext_tests_nlo/nlo_ee4j.sin (revision 8857) @@ -1,16 +1,18 @@ include("nlo_settings.sin") +$openloops_allowed_libs = "pplljj" + 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 sum Pt/2 [@clustered_jets] 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:480:"gw" }