1 | <html> |
---|
2 | <head> |
---|
3 | <title>Error Checks</title> |
---|
4 | <link rel="stylesheet" type="text/css" href="pythia.css"/> |
---|
5 | <link rel="shortcut icon" href="pythia32.gif"/> |
---|
6 | </head> |
---|
7 | <body> |
---|
8 | |
---|
9 | <h2>Error Checks</h2> |
---|
10 | |
---|
11 | There is a few settings related to error checking during program |
---|
12 | execution. Many other checks are performed as well, but do not |
---|
13 | have any specific setting related to themselves. |
---|
14 | |
---|
15 | <p/><code>flag </code><strong> Check:abortIfVeto </strong> |
---|
16 | (<code>default = <strong>off</strong></code>)<br/> |
---|
17 | There are a few ways in which an event can be vetoed, the most |
---|
18 | common being a <a href="UserHooks.html" target="page">User Hooks</a> test. |
---|
19 | Normally this will simply mean that the next parton-level |
---|
20 | configuration is selected inside the <code>Pythia::next()</code> |
---|
21 | routine, without any need for a user intervention. With this |
---|
22 | option switched on, however, <code>Pythia::next()</code> will |
---|
23 | return <code>false</code>. It is then up to the user to decide |
---|
24 | what to do next. |
---|
25 | |
---|
26 | |
---|
27 | <p/><code>flag </code><strong> Check:particleData </strong> |
---|
28 | (<code>default = <strong>off</strong></code>)<br/> |
---|
29 | Check the particle data tables for potential problems during |
---|
30 | initialization. This includes inconsistent use of charge in particle |
---|
31 | names, inconsistent setup of mass, mass range, width and lifetime, |
---|
32 | sum of branching ratios not unity (allowed but discouraged) or charge |
---|
33 | not conserved in a decay channel. Warnings should be viewed as reasons |
---|
34 | to check further, but need not indicate a true problem, and also not all |
---|
35 | problems may be caught. |
---|
36 | The <code>pythia.particleData.checkTable(level)</code> method, |
---|
37 | used for these checks, may also be called directly. |
---|
38 | |
---|
39 | |
---|
40 | <p/><code>mode </code><strong> Check:levelParticleData </strong> |
---|
41 | (<code>default = <strong>1</strong></code>)<br/> |
---|
42 | The level of verbosity and checks of particle data, if switched on. |
---|
43 | <br/><code>option </code><strong> 0</strong> : mimimal amount of checks, e.g. that no channels open. |
---|
44 | |
---|
45 | <br/><code>option </code><strong> 1</strong> : further warning if individual channels closed, |
---|
46 | except for resonances. |
---|
47 | <br/><code>option </code><strong> 2</strong> : also print branching-ratio-averaged threshold mass |
---|
48 | except for resonances. |
---|
49 | <br/><code>option </code><strong> 11</strong> : as 1, but include resonances in detailed checks. |
---|
50 | |
---|
51 | <br/><code>option </code><strong> 12</strong> : as 2, but include resonances in detailed checks. |
---|
52 | |
---|
53 | |
---|
54 | |
---|
55 | <p/><code>flag </code><strong> Check:event </strong> |
---|
56 | (<code>default = <strong>on</strong></code>)<br/> |
---|
57 | When an event has been successfully generated, check that the |
---|
58 | final event record in <code>event</code> does not contain any |
---|
59 | unphysical particles, or nonconserved charge or energy-momentum. |
---|
60 | If this check fails, then <code>pythia.next()</code> obtains the |
---|
61 | value <code>false</code>, i.e. the event is counted as aborted. |
---|
62 | |
---|
63 | |
---|
64 | <p/><code>flag </code><strong> Check:history </strong> |
---|
65 | (<code>default = <strong>on</strong></code>)<br/> |
---|
66 | When <code>Check:event = on</code> and the event is checked as above, |
---|
67 | further checks are made that all mother and daughter pointers are |
---|
68 | consistently set. Specifically that all daughters in the |
---|
69 | <code>daughterlist</code> (or <code>motherList</code>) have the |
---|
70 | particle in their respective <code>motherList</code> |
---|
71 | (or <code>daughterlist</code>). This operation takes a bit more |
---|
72 | time than the other error tests (of the order of 10% of what it takes |
---|
73 | to generate the event in the first place), and so could be switched off |
---|
74 | to save time. |
---|
75 | |
---|
76 | |
---|
77 | |
---|
78 | <p/><code>mode </code><strong> Check:nErrList </strong> |
---|
79 | (<code>default = <strong>0</strong></code>)<br/> |
---|
80 | The number of erroneous events, in the above check, for which |
---|
81 | event listing and other detailed information will be printed. |
---|
82 | After that, only the normal error messages will be issued. |
---|
83 | Error counters are always updated, and accumulated numbers can be |
---|
84 | shown with <code>pythia.statistics()</code> at the end of the run. |
---|
85 | |
---|
86 | |
---|
87 | <p/><code>parm </code><strong> Check:epTolErr </strong> |
---|
88 | (<code>default = <strong>1e-4</strong></code>)<br/> |
---|
89 | Maximum allowed summed deviation of <i>E</i>, <i>p_x</i>, |
---|
90 | <i>p_y</i> and <i>p_z</i> between the incoming beams and the |
---|
91 | final state, as a fraction of the initial energy, above which the |
---|
92 | event is counted as aborted. |
---|
93 | (Unfortunetely roundoff errors do not scale linearly with the energy, |
---|
94 | and also have a very long tail. So while most events at lower energies |
---|
95 | may be correct to better than 1e-10, at LHC it does not have to signal |
---|
96 | any fundamental bug if also the default tolerance above is violated |
---|
97 | occasionally.) |
---|
98 | |
---|
99 | |
---|
100 | <p/><code>parm </code><strong> Check:epTolWarn </strong> |
---|
101 | (<code>default = <strong>1e-6</strong></code>)<br/> |
---|
102 | A check on the same summed deviation as above, but counted as a |
---|
103 | warning rather than an error, and not leading to the event being |
---|
104 | classified as aborted. |
---|
105 | |
---|
106 | |
---|
107 | </body> |
---|
108 | </html> |
---|
109 | |
---|
110 | <!-- Copyright (C) 2012 Torbjorn Sjostrand --> |
---|