c This is the control file for the GEANT simulation.  Parameters defined
c in this file control the kind and extent of simulation that is performed.
c The full list of options is given in section BASE-40 of the GEANT manual.
c
TRIG 1000000

c The following card enables single-track generation (for testing).
c For a single-particle gun, set the momentum (GeV/c), direction
c theta,phi (degrees) and vertex position (cm), and for the particle
c type insert the Geant particle type code plus 100 (eg. 101=gamma,
c 103=electron, 107=pi0, 108=pi+, 109=pi-, 114=proton).  If you use
c the particle code but do not add 100 then theta,phi are ignored
c and the particle direction is generated randomly over 4pi sr.
c For a listing of the Geant particle types, see the following URL.
c http://wwwasdoc.web.cern.ch/wwwasdoc/geant_html3/node72.html
c The meaning of the arguments to KINE are as follows.
c  - particle = GEANT particle type of primary track + 100
c  - momentum = initial track momentum, central value (GeV/c)
c  - theta = initial track polar angle, central value (degrees)
c  - phi = initial track azimuthal angle, central value (degrees)
c  - delta_momentum = spread in initial track momentum, full width (GeV/c)
c  - delta_theta = spread in initial track polar angle, full width (degrees)
c  - delta_phi = spread in initial track azimuthal angle, full width (degrees)
c
c   particle  momentum  theta  phi  delta_momentum delta_theta delta_phi
KINE   50      2e-9      90.   0.      0.             120.        120.

c The SCAP card determines the vertex position for the particle gun.  It
c supports the following three arguments, all of which default to 0.
c
c    vertex_x vertex_y vertex_z
SCAP    1.0       0.      0.

c The following card seeds the random number generator, though it may be
c overridden if seeds are found in the input file (see below). It must be
c unique for each run.  There are two ways to specify the random seed here.
c  1. One argument, must be an integer in the range [1,215]
c  2. Two arguments, must be a pair of positive Integer*4 numbers
c In the first case, one of a limited set of prepared starting seeds is
c chosen from a list.  These seeds have been certified to produce random
c sequences that do not repeat within the first 10^9 or so random numbers.
c For cases where more choices are needed, the two-argument form gives
c access to a total of 2^62 choices, with no guarantees about closed loops.
c
c NOTE: If one uses events read from an HDDM file and that file contains
c random number seeds for the event, those seeds will be used, overwriting
c any value(s) specified here. Most event generators do not include the
c seeds. The seeds are written to the output HDDM file though so if one
c uses the output file for input to another invocation of hdgeant(++)
c then the same seeds will be used. You may check for seeds in the input
c file using hddm-xml file.hddm | grep random .
RNDM 121

c The following line controls the cutoffs for tracking of particles.
c CUTS cutgam cutele cutneu cuthad cutmuo bcute bcutm dcute dcutm ppcutm tofmax
c  - cutgam = Cut for gammas (0.001 GeV)
c  - cutele = Cut for electrons (0.001 GeV)
c  - cutneu = Cut for neutral hadrons (0.01 GeV)
c  - cuthad = Cut for charged hadrons (0.01 GeV)
c  - cutmuo = Cut for muons (0.01 GeV)
c  - bcute  = Cut for electron brems. (CUTGAM)
c  - bcutm  = Cut for muon brems. (CUTGAM)
c  - dcute  = Cut for electron delta-rays. (10 TeV)
c  - dcutm  = Cut for muon delta-rays. (10 TeV)
c  - ppcutm = Cut for e+e- pairs by muons. (0.01 GeV)
c  - tofmax = Time of flight cut (1.E+10 sec)
c  - gcuts  = 5 user words (0.)
CUTS 1e-4 1e-4 1e-3 1e-3 1e-4

c The following line controls a set of generic flags that are used to
c control aspects of the simulation generally related to debugging.
c For normal debugging runs these should be left at zero (or omitted).
c At present the following functionality is defined (assumes debug on).
c  SWIT(2) = 0 turns off trajectory tracing
c          = 2 turns on step-by-step trace during tracking (verbose!)
c          = 3 turns on trajectory plotting after tracking is done
c          = 4 turns on step-by-step plotting during tracking
c  SWIT(3) = 1 stores track trajectories for plotting after tracking is done
c  SWIT(4) = 0 trace trajectories of all particle types
c          = 3 trace only charged particle trajectories
SWIT 0 0 0 0 0 0 0 0 0 0

c The following card enables the GelHad package (from BaBar)
c   on/off  ecut   scale   mode   thresh
GELH  1     0.2     1.0     4     0.160

c The following card selects the hadronic physics package
c   HADR 0	no hadronic interactions
c   HADR 1	GHEISHA only (default)
c   HADR 2	GHEISHA only, with no generation of secondaries
c   HADR 3      FLUKA (with GHEISHA for neutrons below 20MeV)
c   HADR 4	FLUKA (with MICAP for neutrons below 20MeV)
HADR 1

c The following cards are needed if optical photons are being
c being generated and tracked in the simulation.  The CKOV directive
c enables Cerenkov generation in materials for which the refractive
c index table has been specified.  The LABS card enables absorption
c of optical photons.  The ABAN directive controls a special feature
c of Geant which allows it to "abandon" tracking of charged particles
c once their remaining range drops below the distance to the next
c discrete interaction or geometric boundary.  Particles abandoned
c during tracking are stopped immediately and dump all remaining energy
c where they lie.  The remaining energy is dumped in the correct volume
c so this is OK in most cases, but it can cut into the yield of
c Cerenkov photons (eg. in a lead glass calorimeter) at the end of
c a particle track.  If this might be important, set ABAN to 0.
CKOV 1
LABS 1

c The following card prevents GEANT tracking code from abandoning the
c tracking of particles near the end of their range, once it determines
c that their fate is just to stop (i.e. electrons and protons).  This
c behaviour is normal in most cases, but in the case of Cerenkov light
c generation it leads to an underestimate for the yields.
c   ABAN 1	abandon stopping tracks (default)
c   ABAN 0	do not abandon stopping tracks
ABAN 0

c The following card sets up the simulation to perform debugging on
c a subset of the simulated events.
c DEBUG first last step
c  - first (int) = event number of first event to debug
c  - last (int) = event number of last event to debug
c  - step (int) = only debug one event every step events
DEBU 1 10 1000000

c The following tracking parameters are defined for each tracking medium
c   TMAXFD (REAL) maximum angular deviation due to the magnetic field
c                 permitted in one step (degrees)
c   DEEMAX (REAL) maximum fractional energy loss in one step (0< DEEMAX <=0.1)
c   STEMAX (REAL) maximum step permitted (cm)
c   STMIN  (REAL) minimum value for the maximum step imposed by energy loss,
c                 multiple scattering, Cerenkov or magnetic field effects (cm)
c Normally they are assigned appropriate values calculated automatically by
c Geant when the geometry is defined, overwriting the values declared by
c the user code in the GSTMED() call.  Users who know what they are doing can
c force Geant to instead use the values passed in the arguments to GSTMED()
c by removing the comment in front of the following card.  Any parameters with
c zero values are still assigned automatic values even when AUTO is turned off.
cAUTO 0

c The following cards allow one to switch on/off some physics processes in GEANT:
c MULS 0 no multiple scattering
c      1 Moliere or Coulomb scattering (default)  
c 
c BREM 0 no bremsstrahlung 
c      1 bremsstrahlung (default)
c
c COMP 0 no Compton
c      1 Compton scattering (default)
c 
c PAIR 0 no pair production
c      1 pair production (default)
c 
c LOSS 0 (controls energy losses) no energy loss
c      1  delta-rays are produced above the threshold. Reduced fluctuations from 
c         delta-rays below the threshold are added to the energy losses. The threshold
c         energies for delta-ray production can be set using the CUTS card (see above). 
c         The fields 'dcute' and 'dcutm' in the CUTS card correspond to energy thresholds 
c         for electron and muon delta-rays, respectively. The default energy threshold 
c         value is 100 keV (default  see uginit.F  12/16/2011  DL). 
c      2  no delta-rays are produced. Complete fluctuations are calculated .
c
c DCAY 0 no decay in flight
c      1 decay in flight with generation of secondaries (default)
c      2 decay in flight without generation of secondaries
c
c DRAY 0 no delta ray production
c      1 delta ray production with generation of secondaries (default)
c      2 delta ray production without generation of secondaries

END