<html> <head> <title>Error Checks</title> <link rel="stylesheet" type="text/css" href="pythia.css"/> <link rel="shortcut icon" href="pythia32.gif"/> </head> <body> <h2>Error Checks</h2> There is a few settings related to error checking during program execution. Many other checks are performed as well, but do not have any specific setting related to themselves. <p/><code>flag </code><strong> Check:abortIfVeto </strong> (<code>default = <strong>off</strong></code>)<br/> There are a few ways in which an event can be vetoed, the most common being a <a href="UserHooks.html" target="page">User Hooks</a> test. Normally this will simply mean that the next parton-level configuration is selected inside the <code>Pythia::next()</code> routine, without any need for a user intervention. With this option switched on, however, <code>Pythia::next()</code> will return <code>false</code>. It is then up to the user to decide what to do next. <p/><code>flag </code><strong> Check:particleData </strong> (<code>default = <strong>off</strong></code>)<br/> Check the particle data tables for potential problems during initialization. This includes inconsistent use of charge in particle names, inconsistent setup of mass, mass range, width and lifetime, sum of branching ratios not unity (allowed but discouraged) or charge not conserved in a decay channel. Warnings should be viewed as reasons to check further, but need not indicate a true problem, and also not all problems may be caught. The <code>pythia.particleData.checkTable(level)</code> method, used for these checks, may also be called directly. <p/><code>mode </code><strong> Check:levelParticleData </strong> (<code>default = <strong>1</strong></code>)<br/> The level of verbosity and checks of particle data, if switched on. <br/><code>option </code><strong> 0</strong> : minimal amount of checks, e.g. that no channels open. <br/><code>option </code><strong> 1</strong> : further warning if individual channels closed, except for resonances. <br/><code>option </code><strong> 2</strong> : also print branching-ratio-averaged threshold mass except for resonances. <br/><code>option </code><strong> 11</strong> : as 1, but include resonances in detailed checks. <br/><code>option </code><strong> 12</strong> : as 2, but include resonances in detailed checks. <p/><code>flag </code><strong> Check:event </strong> (<code>default = <strong>on</strong></code>)<br/> When an event has been successfully generated, check that the final event record in <code>event</code> does not contain any unphysical particles, or nonconserved charge or energy-momentum. If this check fails, then <code>pythia.next()</code> obtains the value <code>false</code>, i.e. the event is counted as aborted. <p/><code>flag </code><strong> Check:history </strong> (<code>default = <strong>on</strong></code>)<br/> When <code>Check:event = on</code> and the event is checked as above, further checks are made that all mother and daughter pointers are consistently set. Specifically that all daughters in the <code>daughterlist</code> (or <code>motherList</code>) have the particle in their respective <code>motherList</code> (or <code>daughterlist</code>). This operation takes a bit more time than the other error tests (of the order of 10% of what it takes to generate the event in the first place), and so could be switched off to save time. <p/><code>mode </code><strong> Check:nErrList </strong> (<code>default = <strong>0</strong></code>)<br/> The number of erroneous events, in the above check, for which event listing and other detailed information will be printed. After that, only the normal error messages will be issued. Error counters are always updated, and accumulated numbers can be shown with <code>pythia.statistics()</code> at the end of the run. <p/><code>parm </code><strong> Check:epTolErr </strong> (<code>default = <strong>1e-4</strong></code>)<br/> Maximum allowed summed deviation of <i>E</i>, <i>p_x</i>, <i>p_y</i> and <i>p_z</i> between the incoming beams and the final state, as a fraction of the initial energy, above which the event is counted as aborted. (Unfortunately roundoff errors do not scale linearly with the energy, and also have a very long tail. So while most events at lower energies may be correct to better than 1e-10, at LHC it does not have to signal any fundamental bug if also the default tolerance above is violated occasionally.) <p/><code>parm </code><strong> Check:epTolWarn </strong> (<code>default = <strong>1e-6</strong></code>)<br/> A check on the same summed deviation as above, but counted as a warning rather than an error, and not leading to the event being classified as aborted. </body> </html> <!-- Copyright (C) 2013 Torbjorn Sjostrand -->