|=============================================================================| | | | WW WW WW WW WW WWWWWW WW WWWWW WWWW | | WW WW WW WW WW WW WW WWWW WW WW WW WW | | WW WW WW WW WWWWWWW WW WW WW WW WWWWW WW WW | | WWWW WWWW WW WW WW WW WWWWWWWW WW WW WW WW | | WW WW WW WW WW WWWWWW WW WW WW WW WWWW | | | | | | W | | sW | | WW | | sWW | | WWW | | wWWW | | wWWWW | | WW WW | | WW WW | | wWW WW | | wWW WW | | WW WW | | WW WW | | WW WW | | WW WW | | WW WW | | WW WW | | wwwwww WW WW | | WWWWWww WW WW | | WWWWWwwwww WW WW | | wWWWwwwwwWW WW | | wWWWWWWWWWWwWWW WW | | wWWWWW wW WWWWWWW | | WWWW wW WW wWWWWWWWwww | | WWWW wWWWWWWWwwww | | WWWW WWWW WWw | | WWWWww WWWW | | WWWwwww WWWW | | wWWWWwww wWWWWW | | WwwwwwwwwWWW | | | | | | | | by: Wolfgang Kilian | | Thorsten Ohl | | Juergen Reuter | | | | if you use WHIZARD please cite: | | W. Kilian, T. Ohl, J. Reuter, arXiv: 0708.4233 [hep-ph] | | M. Moretti, T. Ohl, J. Reuter, arXiv: hep-ph/0102195 | | | |=============================================================================| | WHIZARD 2.0.4 |=============================================================================| | Initializing process library 'processes' | Loading process library 'processes' | Reading model file '/opt/whizard/share/whizard/models/SM.mdl' | Using model: SM | Reading commands from file 'LEP_higgs.sin' [user variable] n = PDG(12, 14, 16) [user variable] N = PDG(-12, -14, -16) [user variable] q = PDG(2, 1, 3, 4) [user variable] Q = PDG(-2, -1, -3, -4) | Process 'zh': keeping configuration | [L] zh = e-, e+ => Z, H | Process 'nnbb': keeping configuration | [L] nnbb = e-, e+ => nue:numu:nutau, nuebar:numubar:nutaubar, b, bbar | Process 'qqbb': keeping configuration | [L] qqbb = e-, e+ => u:d:s:c, ubar:dbar:sbar:cbar, b, bbar | Process 'bbbb': keeping configuration | [L] bbbb = e-, e+ => b, bbar, b, bbar | Process 'eebb': keeping configuration | [L] eebb = e-, e+ => e-, e+, b, bbar | Process 'qqtt': keeping configuration | [L] qqtt = e-, e+ => u:d:s:c, ubar:dbar:sbar:cbar, tau-, tau+ | Process 'bbtt': keeping configuration | [L] bbtt = e-, e+ => b, bbar, tau-, tau+ | Generating code for process library 'processes' | Skipping process 'zh' (source code exists) | Skipping process 'nnbb' (source code exists) | Skipping process 'qqbb' (source code exists) | Skipping process 'bbbb' (source code exists) | Skipping process 'eebb' (source code exists) | Skipping process 'qqtt' (source code exists) | Skipping process 'bbtt' (source code exists) | Writing interface code for process library 'processes' | Unloading process library 'processes' | Compiling process library 'processes' | Skipping process 'zh' (object code exists) | Skipping process 'nnbb' (object code exists) | Skipping process 'qqbb' (object code exists) | Skipping process 'bbbb' (object code exists) | Skipping process 'eebb' (object code exists) | Skipping process 'qqtt' (object code exists) | Skipping process 'bbtt' (object code exists) | Skipping library 'processes' (no processes have been recompiled) | Loading process library 'processes' | Process 'zh': keeping configuration | Process 'nnbb': keeping configuration | Process 'qqbb': keeping configuration | Process 'bbbb': keeping configuration | Process 'eebb': keeping configuration | Process 'qqtt': keeping configuration | Process 'bbtt': keeping configuration sqrts = 209.00000000000000 SM.mH = 115.00000000000000 SM.wH = 3.22800000000000040E-003 SM.mb = 2.8999999999999999 SM.me = 0.0000000000000000 SM.ms = 0.0000000000000000 SM.mc = 0.0000000000000000 | Initializating integration for process zh: | Read phase-space configuration from file 'zh.phs'... | ... found 1 phase space channels, collected in 1 groves. | Phase space: found 1 equivalences between channels. | iterations = 5:5000 | Applying user-defined cuts. | Reading integration grids and results from file 'zh.vg': | 5000 calls, 1 channels, 2 dimensions, 0 bins, stratified = T | Integrating process 'zh': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 5000 1.6114132E+02 5.36E-04 0.00 0.00* 99.43 2 5000 1.6114176E+02 3.44E-04 0.00 0.00* 45.49 3 5000 1.6114156E+02 2.70E-04 0.00 0.00* 45.49 4 5000 1.6114164E+02 2.30E-04 0.00 0.00* 45.49 5 5000 1.6114172E+02 2.03E-04 0.00 0.00* 45.49 |-----------------------------------------------------------------------------| 5 25000 1.6114165E+02 1.21E-04 0.00 0.00 45.49 0.18 5 |-----------------------------------------------------------------------------| |=============================================================================| 5 25000 1.6114165E+02 1.21E-04 0.00 0.00 45.49 0.18 5 |=============================================================================| | Initializating integration for process nnbb: | Read phase-space configuration from file 'nnbb.phs'... | ... found 51 phase space channels, collected in 12 groves. | Phase space: found 71 equivalences between channels. | iterations = 12:20000 | Applying user-defined cuts. | Reading integration grids and results from file 'nnbb.vg': | 20000 calls, 51 channels, 8 dimensions, 0 bins, stratified = T | Integrating process 'nnbb': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 20000 1.0435458E+02 2.33E+00 2.23 3.16* 4.31 2 20000 1.0621482E+02 1.06E+00 0.99 1.41* 4.33 3 20000 1.0625233E+02 6.43E-01 0.61 0.86* 4.29 4 20000 1.0626800E+02 4.52E-01 0.43 0.60* 4.28 5 20000 1.0638305E+02 3.52E-01 0.33 0.47* 4.29 6 20000 1.0631941E+02 2.93E-01 0.28 0.39* 4.26 7 20000 1.0648900E+02 2.53E-01 0.24 0.34* 4.29 8 20000 1.0685861E+02 2.26E-01 0.21 0.30* 4.33 9 20000 1.0665972E+02 2.06E-01 0.19 0.27* 4.22 10 20000 1.0667204E+02 1.90E-01 0.18 0.25* 4.27 11 20000 1.0664100E+02 1.77E-01 0.17 0.24* 4.24 12 20000 1.0674402E+02 1.67E-01 0.16 0.22* 4.29 |-----------------------------------------------------------------------------| 12 240000 1.0662719E+02 7.41E-02 0.07 0.34 4.29 0.51 12 |-----------------------------------------------------------------------------| |=============================================================================| 12 240000 1.0662719E+02 7.41E-02 0.07 0.34 4.29 0.51 12 |=============================================================================| | Initializating integration for process qqbb: | Read phase-space configuration from file 'qqbb.phs'... | ... found 100 phase space channels, collected in 12 groves. | Phase space: found 164 equivalences between channels. | iterations = 12:20000 | Applying user-defined cuts. | Reading integration grids and results from file 'qqbb.vg': | 20000 calls, 100 channels, 8 dimensions, 0 bins, stratified = T | Integrating process 'qqbb': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 20000 4.7173507E+02 1.19E+01 2.52 3.57* 4.09 2 20000 4.8616899E+02 5.56E+00 1.14 1.62* 4.00 3 20000 4.9380742E+02 3.50E+00 0.71 1.00* 4.01 4 20000 4.9105816E+02 2.46E+00 0.50 0.71* 3.91 5 20000 4.9305346E+02 1.94E+00 0.39 0.56* 3.96 6 20000 4.9414758E+02 1.61E+00 0.33 0.46* 3.95 7 20000 4.9389827E+02 1.39E+00 0.28 0.40* 3.92 8 20000 4.9358988E+02 1.24E+00 0.25 0.35* 3.89 9 20000 4.9411605E+02 1.11E+00 0.23 0.32* 3.91 10 20000 4.9335069E+02 1.01E+00 0.21 0.29* 3.86 11 20000 4.9374536E+02 9.32E-01 0.19 0.27* 3.90 12 20000 4.9346985E+02 8.66E-01 0.18 0.25* 3.87 |-----------------------------------------------------------------------------| 12 240000 4.9353448E+02 3.96E-01 0.08 0.39 3.87 0.62 12 |-----------------------------------------------------------------------------| |=============================================================================| 12 240000 4.9353448E+02 3.96E-01 0.08 0.39 3.87 0.62 12 |=============================================================================| | Initializating integration for process bbbb: | Read phase-space configuration from file 'bbbb.phs'... | ... found 100 phase space channels, collected in 12 groves. | Phase space: found 640 equivalences between channels. | iterations = 12:20000 | Applying user-defined cuts. | Reading integration grids and results from file 'bbbb.vg': | 20000 calls, 100 channels, 8 dimensions, 0 bins, stratified = T | Integrating process 'bbbb': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 20000 6.4007188E+01 1.38E+00 2.15 3.04* 4.61 2 20000 6.5070484E+01 6.13E-01 0.94 1.33* 4.51 3 20000 6.5731606E+01 3.89E-01 0.59 0.84* 4.56 4 20000 6.5967644E+01 2.84E-01 0.43 0.61* 4.56 5 20000 6.5810692E+01 2.25E-01 0.34 0.48* 4.52 6 20000 6.5780562E+01 1.89E-01 0.29 0.41* 4.53 7 20000 6.5901951E+01 1.65E-01 0.25 0.35* 4.57 8 20000 6.5908997E+01 1.48E-01 0.22 0.32* 4.54 9 20000 6.5862913E+01 1.35E-01 0.20 0.29* 4.52 10 20000 6.5786926E+01 1.25E-01 0.19 0.27* 4.50 11 20000 6.5798729E+01 1.16E-01 0.18 0.25* 4.54 12 20000 6.5816762E+01 1.10E-01 0.17 0.24* 4.54 |-----------------------------------------------------------------------------| 12 240000 6.5825683E+01 4.83E-02 0.07 0.36 4.54 0.40 12 |-----------------------------------------------------------------------------| |=============================================================================| 12 240000 6.5825683E+01 4.83E-02 0.07 0.36 4.54 0.40 12 |=============================================================================| | Initializating integration for process eebb: | Read phase-space configuration from file 'eebb.phs'... | ... found 65 phase space channels, collected in 14 groves. | Phase space: found 97 equivalences between channels. | iterations = 12:20000 | Applying user-defined cuts. | Reading integration grids and results from file 'eebb.vg': | 20000 calls, 65 channels, 8 dimensions, 0 bins, stratified = T | Integrating process 'eebb': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 20000 1.0055936E+03 6.18E+02 61.50 86.97* 0.41 2 20000 9.8602291E+02 2.76E+02 28.02 39.63* 0.24 3 20000 1.0493581E+03 2.71E+02 25.84 36.54* 0.32 4 20000 1.3279782E+03 2.22E+02 16.69 23.60* 0.21 5 20000 1.3522964E+03 2.20E+02 16.31 23.06* 0.23 6 20000 1.6082186E+03 1.85E+02 11.53 16.31* 0.13 7 20000 1.8012326E+03 1.51E+02 8.40 11.89* 0.13 8 20000 2.0743178E+03 1.27E+02 6.14 8.69* 0.16 9 20000 2.3474142E+03 1.03E+02 4.40 6.23* 0.17 10 20000 2.3541329E+03 1.03E+02 4.38 6.20* 0.23 11 20000 2.4191849E+03 1.02E+02 4.21 5.95* 0.20 12 20000 2.4302107E+03 1.02E+02 4.18 5.91* 0.20 |-----------------------------------------------------------------------------| 12 240000 2.1308397E+03 4.13E+01 1.94 9.50 0.20 9.15 12 |-----------------------------------------------------------------------------| |=============================================================================| 12 240000 2.1308397E+03 4.13E+01 1.94 9.50 0.20 9.15 12 |=============================================================================| | Initializating integration for process qqtt: | Read phase-space configuration from file 'qqtt.phs'... | ... found 50 phase space channels, collected in 12 groves. | Phase space: found 82 equivalences between channels. | iterations = 12:20000 | Applying user-defined cuts. | Reading integration grids and results from file 'qqtt.vg': | 20000 calls, 50 channels, 8 dimensions, 0 bins, stratified = T | Integrating process 'qqtt': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 20000 1.0412719E+02 1.80E+00 1.73 2.44* 5.88 2 20000 1.0543703E+02 7.75E-01 0.74 1.04* 5.86 3 20000 1.0465655E+02 4.86E-01 0.46 0.66* 5.74 4 20000 1.0433101E+02 3.54E-01 0.34 0.48* 5.71 5 20000 1.0441592E+02 2.94E-01 0.28 0.40* 5.54 6 20000 1.0413029E+02 2.49E-01 0.24 0.34* 5.43 7 20000 1.0417818E+02 2.17E-01 0.21 0.29* 5.47 8 20000 1.0401866E+02 1.94E-01 0.19 0.26* 5.41 9 20000 1.0401723E+02 1.75E-01 0.17 0.24* 5.44 10 20000 1.0404310E+02 1.68E-01 0.16 0.23* 5.19 11 20000 1.0392141E+02 1.57E-01 0.15 0.21* 5.11 12 20000 1.0382861E+02 1.46E-01 0.14 0.20* 5.11 |-----------------------------------------------------------------------------| 12 240000 1.0404050E+02 6.39E-02 0.06 0.30 5.11 0.95 12 |-----------------------------------------------------------------------------| |=============================================================================| 12 240000 1.0404050E+02 6.39E-02 0.06 0.30 5.11 0.95 12 |=============================================================================| | Initializating integration for process bbtt: | Read phase-space configuration from file 'bbtt.phs'... | ... found 50 phase space channels, collected in 12 groves. | Phase space: found 82 equivalences between channels. | iterations = 12:20000 | Applying user-defined cuts. | Reading integration grids and results from file 'bbtt.vg': | 20000 calls, 50 channels, 8 dimensions, 0 bins, stratified = T | Integrating process 'bbtt': |=============================================================================| | It Calls Integral[fb] Error[fb] Err[%] Acc Eff[%] Chi2 N[It] | |=============================================================================| 1 20000 3.0839637E+01 5.40E-01 1.75 2.48* 6.06 2 20000 3.0185227E+01 2.16E-01 0.72 1.01* 5.78 3 20000 3.0220917E+01 1.33E-01 0.44 0.62* 5.80 4 20000 3.0199659E+01 9.41E-02 0.31 0.44* 5.79 5 20000 3.0221076E+01 7.35E-02 0.24 0.34* 5.80 6 20000 3.0241962E+01 6.13E-02 0.20 0.29* 5.80 7 20000 3.0268504E+01 5.36E-02 0.18 0.25* 5.79 8 20000 3.0244126E+01 4.81E-02 0.16 0.22* 5.75 9 20000 3.0260689E+01 4.41E-02 0.15 0.21* 5.73 10 20000 3.0276479E+01 4.21E-02 0.14 0.20* 5.47 11 20000 3.0280962E+01 3.93E-02 0.13 0.18* 5.44 12 20000 3.0262671E+01 3.69E-02 0.12 0.17* 5.40 |-----------------------------------------------------------------------------| 12 240000 3.0260165E+01 1.60E-02 0.05 0.26 5.40 0.25 12 |-----------------------------------------------------------------------------| |=============================================================================| 12 240000 3.0260165E+01 1.60E-02 0.05 0.26 5.40 0.25 12 |=============================================================================| $description = "A WHIZARD 2.0 Example. Using weighted events to speed things up." $y_label = "$N_{\textrm{events}}$" $title = "Invisible mass distribution in $e^+e^- \to \nu\bar\nu b \bar b$" $x_label = "$M_{\nu\nu}$/GeV" $title = "$bb$ invariant mass distribution in $e^+e^- \to \nu\bar\nu b \bar b$" $x_label = "$M_{b\bar b}$/GeV" luminosity = 10.000000000000000 | Initializing simulation for processes nnbb: | Simulation mode = unweighted, event_normalization = '1', unpolarized_events | Requested number of events = 1066 | This corresponds to luminosity [fb-1] = 10.00 | Reading events from file 'nnbb.evx' ... | Applying user-defined analysis setup. | ... 1066 events read, 1066 total. | Simulation finished. $title = "Dijet invariant mass distribution in $e^+e^- \to q \bar q b \bar b$" $x_label = "$M_{q\bar q}$/GeV" | Initializing simulation for processes qqbb: | Simulation mode = unweighted, event_normalization = '1', unpolarized_events | Requested number of events = 4935 | This corresponds to luminosity [fb-1] = 10.00 | Reading events from file 'qqbb.evx' ... | Applying user-defined analysis setup. | ... 4935 events read, 4935 total. | Simulation finished. $out_file = "lep_higgs.dat" | Writing analysis data to file 'lep_higgs.dat' | Compiling analysis results in file 'lep_higgs.tex' | Compiling analysis results display in 'lep_higgs.tex' | WHIZARD run finished. |=============================================================================|