Options (ENIGMA release 2.8)

Target-related options:

-target <isotope>
-tardens <rho>
-exmax <E>
-exmin <E>
-fermidist <MeV/c>
-fermimom <MeV/c>
-fermimin <MeV/c>
-fermimax <MeV/c>
-fermiwidth <MeV/c>
-thetafermi <d>
-thetafermiwidth <d>
-phifermi <d>
-phifermiwidth <d>
-corpar <cp>

Beam-related options:

-Tbeam <T>
-dTbeam <dT>
-dpoverp <dp>
-time <t>
-lumin <L>

Reaction-related options:

-process <process>
-formalism <theorist>
-omit <term>
-isotropic <on|off>
-xsection <value> [MeV]
-bremsstrahlung <on|off|before|after>
-bremscutoff <kc> [MeV]
-incmframe <on|off>
-thetacm <d> degrees.
-thetacmwidth <d>
-phicm <d>
-phicmwidth <d>

Detector-related options:

-det1 <detname>
-det2 <detname>
-det16 <detname>
-trigmulmin <n>
-trigmulmax <n>

Other options:

-results <filnam>
-eventfile <evfile>
-logfile <logfile>
-tables <path>
-standardinput <on|off>
-incstep <dy>
-maxweight <dy>
-seed <i>
-nsigma <n>
-pmismax <MeV/c>
-pmismin <MeV/c>
-pmisstep <MeV/c>
-Tmismax <MeV>
-Tmismin <MeV>
-Tmisstep <MeV>
-checklimits <on|off>
-condor <on|off>

-target <isotope>

Sets the target nucleus isotope to <isotope>. Conventional notation like 12C, 9Be, 2H etc. is accepted. Mind the Capitals! Default is 12C.

-process <process>

Specifies the reaction type. Available processes are listed in the Description , second column, Mnemonics.

Default is eeppi0.

-Tbeam <T>

Sets the beam kinetic energy to <T> MeV.

Default is 730 MeV (The O'Connell energy).

-dTbeam <dT>

Sets the beam kinetic energy interval to: Tbeam +/- 0.5*<dT> This gives a uniform distribution in beam energies, and is used e.g. for real photon beams.

Default <dT>=0.

-dpoverp <dp>

Sets a Gaussian distribution with variance <dp>*|p| on every component p_i of the beam momentum. This is used for electron beams and can be combined with dTbeam, to get a flat distribution with Gaussian edges.

Default <dp>=0.

-time <t>

Sets the beamtime for your simulated run to <t> seconds. (This is NOT the CPU time, but actual simulated beamtime!).

Default is <t>=1

-logfile <logfile>

Send the logging messages of Enigma to a file named <logfile>, instead of piping them to

Default: standard output.

-lumin <L>

Sets the luminosity to <L> beam+targetnucleus encounters per cm^-2 per second. This option, together with -time, influences the CPU-time needed for the ENIGMA program. Of course the detector acceptances influence the CPU time as well.

Default <L>=1e33 cm^-2.

-tardens <rho>

Specifies the target density as <rho> targetnuclei cm^-2. This parameter does not influence countrates or CPU time. It is only used to calculate the beam intensity from the specified luminosity.

Default td=1e22.

-eventfile <evfile>

Send the Enigma list-mode events to the file <evfile>, instead of sending them to:

Default: standard output

-exmin <E>

Sets the mimimal excitation energy (positive), left in the recoil nucleus ,to <E> MeV. An excitation energy for the recoil nucleus is uniformly drawn between exmin and exmax, and added to the recoiling restmass.

Default <E>=0.

-exmax <E>

see exmin, but now maximal.

Default <E>=0

-maxweight <m>

Freezes the value of maxweight to <m>. No attempt is made to adapt the value of maxweight dynamically. This option is ment for test purposes. Using it may cause a large rejection ratio, or on the other side, a large number of underweighted events.

-seed <i>

Sets the random seed to <i>. If this is set to a positive integer, it will reproduce exactly the statistics of another run with the same seed. If set to a negative integer, a different seed will be derived from the CPU's date and time, for every run. In that case, different runs can be added to improve statistics.

Default <i>=-1.

-tables <path>

Specifies the search path that ENIGMA will use to find the tables of isotope masses, cross sections etc.

Default: <path>=/local/eminutil1/Enigma/tables

-bremsstrahlung <on|off|before|after>

Generate the internal bremsstrahlung tails for incoming and outgoing electrons, according to the factorization theorem see e.g. Borie&Drechsel Nucl Phys A167(1971)369 or earlier and better and not referenced in Borie&Drechsel: Ngoc&Jorba PR136(1964)B1036.
The <on>parameter generates events either in the radiative tail or in the main peak, depending on the setting of the -bremscutoff option. The <before> parameter generates only events where the incoming electron has emitted a real photon with k>k_cutoff, while the <after> parameter does the same for the outgoing electron. In this sense, the parameters should be understood as <only-before> and <only-after>. Correct normalisation is kept, so the <before> and <after> runs should together have a statistically equal number of `radiation tail events' as the <on> run. (neglecting the few events with <before AND after>, they are only present in the <on> run).
Warning: '-bremsstrahlung' may sample some anomalously high cross-sections for large electron energy losses. This may give rise to exceedingly long CPU-times. So restrict your electron energy acceptance as much as possible and keep an eye on the "maxweight" convergence.

Default: off

-bremscutoff <kc> [MeV]

Set the cutoff photon energy for emission of real photons by internal bremsstrahlung. When the foton energy would be below cutoff, and bremsstrahlung is switched on, a non-bremsstrahlung event is generated. In combination with the -bremsstrahlun before or after option, this can lead to big savings in CPU time.

Default: kc=1 MeV

-nsigma <n>

When normalizing, we need to know the maximum possible weight of an event over the experimental acceptance. We define it as the mean weight plus a real number <n> times its variance. Events that exceed this are counted and the result xx is given in the ENIGMA output: 'underweighted hits: xx'. This number xx should not exceed a small fraction of the total number of hits, because these events will not be correctly distributed. If more than a few percent of the hits are underweighted, try to increase <n> from its Default value of 3 to 4 or 5. For crude runs, one may try values of <n> equal to 2.5 or 2, saving some CPU time.

Default: <n>=3

-thetacm <d> degrees.

Always take <d> as theta in the center of mass of the interacting subsystem, instead of (Default) drawing a uniform -1<cos(thetacm)<+1. Also specify thetacmwidth. This and the next four options can be used to restrict the CM-solid angle where events are generated. CM angles are defined in the conventional way: Theta is the polar angle with respect to the 3-momentum transfer direction q. Phicm is the azimuthal angle of the rotation of the decay-plane with respect to the collision plane. Phicm defined to be zero when the both planes are the same, but the pion is forward of q. In case of (ee',pn) decay, phicm=0 when the proton is in plane and more forward than q. This option may speed up production runs for small solid-angle set-ups.

Default: draw a uniform -1<cos(thetacm)<+1.

-thetacmwidth <d>

With thetacm set to <th>, suppress events that are outside the interval: <th>+<d>/2<thetacm<<th>-<d>/2 If <d> is put to zero, correct normalisation is lost.

Default is <d>=180.

-phicm <d>

Always take <d> as phi in the CM of the interacting subsystem instead of the Default.

Default: draw a uniform 0 < phicm < 360.

-phicmwidth <d>

With phicm set to <phi>, suppress events that are outside the interval: <phi>-<d>/2 <phicm< <phi>+<d>/2 If put to zero, correct normalisation is lost.

Default is <d>=360 (full phi coverage)

-incmframe <on|off>

Output the event 4-vectors and the histograms in the CM frame of the interacting cluster. The events are still detected in the Lab frame, so a detection-trigger can still be applied. This option can be used to determine CM-acceptance before restricting the thetacm and phicm ranges.

Default is off.

-incstep <dy>

Set the increment step for histogramming to the (possibly real or negative) value <dy>

Default: <dy>=+1.

-isotropic <on|off>

make the angular decay distribution in CM isotropic by choosing the cross section always equal to the thetacm=90 degrees value.

Default: off.

-xsection <value> [MeV]

Take the cross-section energy-independent and give it the value that it would have at energy<value> MeV. (for the "elementary" process, that is to say before applying fermi motion, and for electron beam it is the virtual gamma energy...)

Default: take energy dependency from formalism.

-thetafermi <d>

always take <d> as theta of the direction of fermimotion of the interacting cluster, instead of drawing a uniform -1< cos(thetafermi) <+1

-thetafermiwidth <d>

draw a uniform cos(thetafermi-0.5*thetafermiwidth)< cos(thetafermi) < thetafermi+0.5*thetafermiwidth)

-phifermi <d>

always take <d> as phi of the fermimotion of the interacting cluster, instead of drawing a uniform 0 < phifermi <2pi.

-phifermiwidth <d>

draw a uniform phifermi-0.5*phifermiwidth < phifermi <phifermi+0.5*phifermiwidth.

-fermidist <MeV/c>

specify the type of fermimomentum distribution for the interacting nucleon (or nucleon cluster). Present choices:
constant
Constant as a function of the magnitude of |p|, over the interval fermimom-0.5*fermiwidth < |p| < fermimom+0.5*fermiwidth.
Gaussian
has a default width of 110 MeV/c (sigma), but this may be changed with the -fermiwidth option.
Urbana2
the momentum distribution in 2H, as given by Sciavilla et al [NuclPhysA449(1986)219, table6.
Urbana3
the proton distribution in 3He [ibidem, table 8].
pwave
is a distribution, due to Leon Kester, for the ~4 p-wave protons in 12C.

Defaults:

2H: Urbana2
3H, 3He: Urbana3
4He: Gaussian with sigma=85MeV/c
other: Gaussian with sigma=110MeV/c

-fermimin <MeV/c>

only accept fermimomenta, drawn according to a given distribution, if they are at least equal to <MeV/c>.

-fermimax <MeV/c>

only accept fermimomenta, drawn according to a given distribution, if they are at most equal to <MeV/c>.

-fermimom <MeV/c>

always take the magnitude of the fermimomentum equal to <MeV/c>, instead of drawing it according to a Gaussian or Urbana type of distribution.

-fermiwidth <MeV/c>

change the variance of the Gaussian magnitude distribution of the fermi momentum.

Default value: 110 to <MeV/c> per nucleon of the interacting cluster. In case of 2H and 3He where Urbana distributions are Default, setting this option will force Gaussian distributions.

-trigmulmin <n>

set the least number of detectors, required for a trigger to <n>.

Default is the total number of triggers having their trigger-mode bit equal to 1 in their detector description file.

-trigmulmax <n>

set the highest number of detectors, allowed for a trigger to <n>.

Default is the total number of triggers having their trigger-mode bit equal to 1 in their detector description file.

-standardinput <on|off>

Instead of generating simulated events, read events from standard input in the ENIGMA-output format. In this way one can compare ENIGMA histograms with experimental histograms, without resorting to PAW or other packages. Some of the normal ENIGMA options will not work (like the % of beamtime done etc).

Default = off.

-det1 <detname>

Installs a detector named <detname> in the simulated experiment. A Detector File named <detname> should exist in the current directory, specifying the detector properties. When <detname>=none is specified, an already defined detector is removed and no new one installed.

Default: none.

-det2 <detname>

see -det1. It is advised to not use the same detectorname twice. I don't know what happens.

So better use a different one here. ... And so on. The present maximum is:

-det16 <detname>

-formalism <theorist>

Select a theoretician to calculate your cross sections:

  1. Present choices for pion production are:
    laget
    Blomqvist-Laget [NuclPhysA280(1977)405] tabulated via the FreeSig code.
    dressler
    E.T.Dressler [CanJPhys66(1988)279] This option uses cross sections tabulated by F.Blaazer(VuA)
    dressler.eps
    The same formalism, but with inline code, and no tables. Code due to E.P.Sichterman (RUU)

    Default: <theorist>=dressler.eps.

  2. For eepn and gampn we offer also several formalisms:
    lorentzian
    An isotropic cross section with a lorentzian energy dependence, adapted to the missing strength in the 'dip-region' in the inclusive missing energy spectra.
    chant
    D(gamma,pn) photodissociation cross section as parametrized by Chant & Roos in their ThreeDee code.
    rossi
    The 2H(gamma,p)n cross section parametrization by Rossi et al [PhysRevC 40(1989)2412].
    We multiply with the number of Quasi-deuterons as given by the Levinger constant [Tavares e.a, JPhysG:NucPartPhys18(1992)521].

    Default theorist for eepn and gampn: <theorist>=rossi.

  3. For pion absorption, we have only one formalism:
    ashery
    the differential pi+,pp cross section is parametrized as a Legendre expansion with energy-dependent coefficients as in Ashery[Ann.Rev.Nucl.Part.Sci26(1986)207].
    3- and 4 nucleon pion absorption processes are assumed to have the same energy dependence, but are taken to be isotropic in their (3N or 4N) CM system.

-omit <term>

Omit a response function term (T or L or P or I). See Dressler formalism. e.g. "-omit P -omit I" will omit the Polarisation and Interference terms which are not needed for inclusive experiments. For real photons only the T term is nonzero, and so Enigma only uses the T term.

Default is to include all four terms for virtual photons, and only the T term for real photons.

-condor <on|off>

When running enigma.condor, this option should be included. Deleting and re-creating the results-subdirectory is suppressed. This is necessary because Condor tolerates no shell commands from within the application. N.B. Be aware that the results/ directory has to pre-exist, if you are generating histograms.

Default = off.

-corpar <cp>

This parameter determines the amount of correlation between nucleons in a cluster that interacts with the beam particle. It may vary between 0 and A, where A is the number of nucleons in the cluster. cp=0 means the cluster has no momentum, the sum of nucleon momenta is zero (completely anti-correlated) cp=A means the sum of nucleon momenta is A*(single nucleon) (completely correlated motion). cp=1 means independent motion: sum=sqrt(A)*individual momentum distribution width.

Default: cp=1

-checklimits <on|off>

Give warnings if histogram or table limits are exceeded. This can occur with read, lookup, increment or write operations on tables, histograms, scatterplots etc. In all cases (on or off) the actual histogram channel is truncated to the lowest (for overflow) or the highest (for overflow) valid value.

Default: <on>.

-results <filnam>

Instead of creating histograms in the directory with the name ./results, do it in ./<filnam>

-Tmismin <MeV>

Set lower bound of missing energy spectra equal to <MeV>.

Default: -Tbeam.

-Tmismax <MeV>

Set upper bound of missing energy spectra equal to <MeV>.

Default: +Tbeam.

-Tmisstep <MeV>

Set lower bound of missing energy spectra equal to <MeV>.

Default: 2*Tbeam/36.

-pmismin <MeV/c>

Set lower bound of missing energy spectra equal to <MeV/c>.

Default: -2*Tbeam.

-pmismax <MeV/c>

Set upper bound of missing energy spectra equal to <MeV/c>.

Default: +2*Tbeam.

-pmisstep <MeV/c>

Set lower bound of missing energy spectra equal to <MeV/c>.

Default: 4*Tbeam/36.


Back to NIKHEF home page.
JV / 6 Jan 1994