Page MenuHomeHEPForge
Files F7879273

No OneTemporary
Actions

View Options

Index: trunk/share/tests/ref-output-double/helicity.ref
===================================================================
--- trunk/share/tests/ref-output-double/helicity.ref (revision 5384)
+++ trunk/share/tests/ref-output-double/helicity.ref (revision 5385)
@@ -1,277 +1,272 @@
?vis_history = false
?integration_timer = false
+?pacify = true
| Process library 'helicity_lib': recorded process 'helicity_ee_mm'
| Process library 'helicity_lib': recorded process 'helicity_aa_ww'
| Process library 'helicity_lib': compiling ...
| Process library 'helicity_lib': writing makefile
| Process library 'helicity_lib': removing old files
| Process library 'helicity_lib': writing driver
| Process library 'helicity_lib': creating source code
| Process library 'helicity_lib': compiling sources
| Process library 'helicity_lib': linking
| Process library 'helicity_lib': loading
| Process library 'helicity_lib': ... success.
-SM.me = 0.000000000000E+00
-SM.mmu = 0.000000000000E+00
-sqrts = 5.000000000000E+02
+SM.me = 0.0000000000E+00
+SM.mmu = 0.0000000000E+00
+sqrts = 5.0000000000E+02
openmp_num_threads = 1
?helicity_selection_active = false
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_ee_mm:
| Beam structure: e-, e+
| Beam data (collision):
| e- (mass = 0.0000000E+00 GeV)
| e+ (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+02 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_ee_mm_i1.phs'
| Phase space: 2 channels, 2 dimensions
| Phase space: found 2 channels, collected in 2 groves.
| Phase space: Using 2 equivalences between channels.
| Phase space: wood
Warning: No cuts have been defined.
| Starting integration for process 'helicity_ee_mm'
| Integrator: 2 chains, 2 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
1 9216 4.2823234E+02 6.92E-02 0.02 0.02* 40.29
|-----------------------------------------------------------------------------|
1 9216 4.2823234E+02 6.92E-02 0.02 0.02 40.29
|=============================================================================|
-[user variable] without_helicity_selection = 4.282323386776E+02
+[user variable] without_helicity_selection = 4.2823233868E+02
?helicity_selection_active = true
helicity_selection_cutoff = 100
-helicity_selection_threshold = 1.000000000000E+10
+helicity_selection_threshold = 1.0000000000E+10
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_ee_mm:
| Beam structure: e-, e+
| Beam data (collision):
| e- (mass = 0.0000000E+00 GeV)
| e+ (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+02 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_ee_mm_i1.phs'
| Phase space: 2 channels, 2 dimensions
| Phase space: found 2 channels, collected in 2 groves.
| Phase space: Using 2 equivalences between channels.
| Phase space: wood
Warning: No cuts have been defined.
| Starting integration for process 'helicity_ee_mm'
| Integrator: 2 chains, 2 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
1 9216 4.2823234E+02 6.92E-02 0.02 0.02* 40.29
|-----------------------------------------------------------------------------|
1 9216 4.2823234E+02 6.92E-02 0.02 0.02 40.29
|=============================================================================|
-[user variable] with____helicity_selection = 4.282323386776E+02
-helicity_selection_threshold = 1.000000000000E+17
+[user variable] with____helicity_selection = 4.2823233868E+02
+helicity_selection_threshold = 1.0000000000E+17
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_ee_mm:
| Beam structure: e-, e+
| Beam data (collision):
| e- (mass = 0.0000000E+00 GeV)
| e+ (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+02 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_ee_mm_i1.phs'
| Phase space: 2 channels, 2 dimensions
| Phase space: found 2 channels, collected in 2 groves.
| Phase space: Using 2 equivalences between channels.
| Phase space: wood
Warning: No cuts have been defined.
| Starting integration for process 'helicity_ee_mm'
| Integrator: 2 chains, 2 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
Warning: ! vamp: Function identically zero in channel 2
1 9216 4.5261726E+00 2.38E-02 0.52 0.50* 0.87
|-----------------------------------------------------------------------------|
1 9216 4.5261726E+00 2.38E-02 0.52 0.50 0.87
|=============================================================================|
-[user variable] crazy___helicity_selection = 4.526172644731E+00
+[user variable] crazy___helicity_selection = 4.5261726447E+00
Comparing results:
-tolerance = 1.000000000000E-12
-?pacify = true
+tolerance = 1.0000000000E-12
[user variable] without_helicity_selection = 4.2823233868E+02
-?pacify = true
[user variable] with____helicity_selection = 4.2823233868E+02
| expect: success
-[user variable] crazy___helicity_selection = 4.526172644731E+00
+[user variable] crazy___helicity_selection = 4.5261726447E+00
| expect: success
********************************************************
* A A => W+ W- @ 5 TeV, w/o cuts, but w/ higher cutoff *
********************************************************
-sqrts = 5.000000000000E+03
+sqrts = 5.0000000000E+03
?helicity_selection_active = false
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_aa_ww:
| Beam structure: A, A
| Beam data (collision):
| A (mass = 0.0000000E+00 GeV)
| A (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_aa_ww_i1.phs'
| Phase space: 4 channels, 2 dimensions
| Phase space: found 4 channels, collected in 1 grove.
| Phase space: Using 8 equivalences between channels.
| Phase space: wood
Warning: No cuts have been defined.
| Starting integration for process 'helicity_aa_ww'
| Integrator: 1 chains, 4 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
1 8192 8.5843811E+04 1.04E+02 0.12 0.11* 24.14
|-----------------------------------------------------------------------------|
1 8192 8.5843811E+04 1.04E+02 0.12 0.11 24.14
|=============================================================================|
-[user variable] without_helicity_selection = 8.584381148085E+04
+[user variable] without_helicity_selection = 8.5843811481E+04
?helicity_selection_active = true
helicity_selection_cutoff = 1000
-helicity_selection_threshold = 1.000000000000E+10
+helicity_selection_threshold = 1.0000000000E+10
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_aa_ww:
| Beam structure: A, A
| Beam data (collision):
| A (mass = 0.0000000E+00 GeV)
| A (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_aa_ww_i1.phs'
| Phase space: 4 channels, 2 dimensions
| Phase space: found 4 channels, collected in 1 grove.
| Phase space: Using 8 equivalences between channels.
| Phase space: wood
Warning: No cuts have been defined.
| Starting integration for process 'helicity_aa_ww'
| Integrator: 1 chains, 4 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
1 8192 8.5843811E+04 1.04E+02 0.12 0.11* 24.14
|-----------------------------------------------------------------------------|
1 8192 8.5843811E+04 1.04E+02 0.12 0.11 24.14
|=============================================================================|
-[user variable] with____helicity_selection = 8.584381148085E+04
+[user variable] with____helicity_selection = 8.5843811481E+04
Comparing results:
-tolerance = 1.000000000000E-12
-?pacify = true
+tolerance = 1.0000000000E-12
[user variable] without_helicity_selection = 8.5843811481E+04
-?pacify = true
[user variable] with____helicity_selection = 8.5843811481E+04
| expect: success
**************************************************
* A A => W+ W- @ 5 TeV, w/ cuts and lower cutoff *
**************************************************
-sqrts = 5.000000000000E+03
+sqrts = 5.0000000000E+03
?helicity_selection_active = false
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_aa_ww:
| Beam structure: A, A
| Beam data (collision):
| A (mass = 0.0000000E+00 GeV)
| A (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_aa_ww_i1.phs'
| Phase space: 4 channels, 2 dimensions
| Phase space: found 4 channels, collected in 1 grove.
| Phase space: Using 8 equivalences between channels.
| Phase space: wood
| Applying user-defined cuts.
| Starting integration for process 'helicity_aa_ww'
| Integrator: 1 chains, 4 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
1 8192 3.4087979E+04 2.77E+02 0.81 0.74* 15.22
|-----------------------------------------------------------------------------|
1 8192 3.4087979E+04 2.77E+02 0.81 0.74 15.22
|=============================================================================|
-[user variable] without_helicity_selection = 3.408797900700E+04
+[user variable] without_helicity_selection = 3.4087979007E+04
?helicity_selection_active = true
helicity_selection_cutoff = 100
-helicity_selection_threshold = 1.000000000000E+10
+helicity_selection_threshold = 1.0000000000E+10
seed = 0
| RNG: Initializing TAO random-number generator
| RNG: Setting seed for random-number generator to 0
| Initializing integration for process helicity_aa_ww:
| Beam structure: A, A
| Beam data (collision):
| A (mass = 0.0000000E+00 GeV)
| A (mass = 0.0000000E+00 GeV)
| sqrts = 5.000000000000E+03 GeV
| Phase space: generating configuration ...
| Phase space: ... success.
| Phase space: writing configuration file 'helicity_aa_ww_i1.phs'
| Phase space: 4 channels, 2 dimensions
| Phase space: found 4 channels, collected in 1 grove.
| Phase space: Using 8 equivalences between channels.
| Phase space: wood
| Applying user-defined cuts.
| Starting integration for process 'helicity_aa_ww'
| Integrator: 1 chains, 4 channels, 2 dimensions
| Integrator: Using VAMP channel equivalences
| Integrator: 10000 initial calls, 20 bins, stratified = T
| Integrator: VAMP
|=============================================================================|
| It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] |
|=============================================================================|
1 8192 3.4087979E+04 2.77E+02 0.81 0.74* 15.22
|-----------------------------------------------------------------------------|
1 8192 3.4087979E+04 2.77E+02 0.81 0.74 15.22
|=============================================================================|
-[user variable] with____helicity_selection = 3.408797900700E+04
+[user variable] with____helicity_selection = 3.4087979007E+04
Comparing results:
-tolerance = 1.000000000000E-10
-?pacify = true
+tolerance = 1.0000000000E-10
[user variable] without_helicity_selection = 3.4087979007E+04
-?pacify = true
[user variable] with____helicity_selection = 3.4087979007E+04
| expect: success
| Summary of value checks:
| Failures: 0 / Total: 4
| There were no errors and 6 warning(s).
| WHIZARD run finished.
|=============================================================================|
Index: trunk/share/tests/helicity.sin
===================================================================
--- trunk/share/tests/helicity.sin (revision 5384)
+++ trunk/share/tests/helicity.sin (revision 5385)
@@ -1,110 +1,111 @@
# $Id$
# Validating the helicity selection rule heuristics
?logging = true
?vis_history = false
?integration_timer = false
+?pacify = true
process helicity_ee_mm = e1, E1 => e2, E2
process helicity_aa_ww = A, A => Wp, Wm
compile
me = 0
mmu = 0
sqrts = 500 GeV
beams = e1, E1
iterations = 1:10000
!!! Tests should be run single-threaded
openmp_num_threads = 1
?helicity_selection_active = false
seed = 0
integrate (helicity_ee_mm)
real without_helicity_selection = integral(helicity_ee_mm)
?helicity_selection_active = true
helicity_selection_cutoff = 100
helicity_selection_threshold = 1E10
seed = 0
integrate (helicity_ee_mm)
real with____helicity_selection = integral(helicity_ee_mm)
helicity_selection_threshold = 1E17
seed = 0
integrate (helicity_ee_mm)
real crazy___helicity_selection = integral(helicity_ee_mm)
printf ("Comparing results:")
tolerance = 1E-12
-show(without_helicity_selection) { ?pacify = true }
-show(with____helicity_selection) { ?pacify = true }
+show(without_helicity_selection)
+show(with____helicity_selection)
expect (with____helicity_selection == without_helicity_selection)
show(crazy___helicity_selection)
if (real_precision < 18) then
expect (crazy___helicity_selection <> without_helicity_selection)
endif
! ----------------------------------------------------------------------
printf ("********************************************************")
printf ("* A A => W+ W- @ 5 TeV, w/o cuts, but w/ higher cutoff *")
printf ("********************************************************")
sqrts = 5 TeV
beams = A, A
iterations = 1:10000
?helicity_selection_active = false
seed = 0
integrate (helicity_aa_ww)
without_helicity_selection = integral(helicity_aa_ww)
?helicity_selection_active = true
helicity_selection_cutoff = 1000
helicity_selection_threshold = 1E10
seed = 0
integrate (helicity_aa_ww)
with____helicity_selection = integral(helicity_aa_ww)
printf ("Comparing results:")
tolerance = 1E-12
-show(without_helicity_selection) { ?pacify = true }
-show(with____helicity_selection) { ?pacify = true }
+show(without_helicity_selection)
+show(with____helicity_selection)
expect (with____helicity_selection == without_helicity_selection)
! ----------------------------------------------------------------------
printf ("**************************************************")
printf ("* A A => W+ W- @ 5 TeV, w/ cuts and lower cutoff *")
printf ("**************************************************")
sqrts = 5 TeV
beams = A, A
cuts = all Pt > 100 GeV [Wp:Wm]
iterations = 1:10000
?helicity_selection_active = false
seed = 0
integrate (helicity_aa_ww)
without_helicity_selection = integral(helicity_aa_ww)
?helicity_selection_active = true
helicity_selection_cutoff = 100
helicity_selection_threshold = 1E10
seed = 0
integrate (helicity_aa_ww)
with____helicity_selection = integral(helicity_aa_ww)
printf ("Comparing results:")
tolerance = 1E-10
-show(without_helicity_selection) { ?pacify = true }
-show(with____helicity_selection) { ?pacify = true }
+show(without_helicity_selection)
+show(with____helicity_selection)
expect (with____helicity_selection == without_helicity_selection)

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 7:48 PM (1 d, 7 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3805916
Default Alt Text
(17 KB)

Event Timeline