next up previous contents
Next: A.5 EXTREAD steering example Up: A. Examples of steering-files Previous: A.3 PYTHIA6 steering example   Contents


A.4 SUSYGEN steering example

  LIST
C
C In this example, all keys know by default by SGV are shown.
C The commented value is the default.

C ========== General steerings : ============


C  Max # of events
C
C         MAXEV 1000

C  Max # of events to list
C
C         MAXPR 1


C  Unit to list events on
C
C         LUNPR 6

C          Starting seed for detector simulation random number generator.
C          On VMS, setting ISEED to 0 yields new seeds
C          each run (the seed will be calculated from the
C          wall-clock
C
C         DSEED 123345678

C          Starting seed for Lund evet generator random number generator.
C          Same comment. Note that other generators might or might not
C          care about this value...
C
C         GSEED 19780503

C  Save random-number generator seeds on file ( to be able to
C  restart after a crash. Only the most current set will be saved)
C
C         SEED_BACKUP TRUE

C  Frequency of seed-backups. Don't put it much less than this default:
C  I/O take time !!!
C
C         BACKUP_FREQ  100

C  Should the analysing code do initialisation ? (normally not, but if
C  You saved simulated events on a file and are reading them back, it
C  should.
C
C         NEED_ANA_INI  FALSE

C  Skip the event and continue with next on errors (TRUE) or stop
C  execution ?
C
C         SKIP TRUE

C  Should histograms read in be updated (1) or reset (0) ?
C
C         UPDF 1

C  Use column-wise ntuples ?
C
C         CWN FALSE

C  Use any HBOOK external file at all ?
C
C         USE_HBOOK TRUE

C  I/O generated and simulated events :
C           RZ   : output to RZ file
C           FZUT : output to FZ file
C           FZIN : input from FZ file
C           NO   : No I/O (normally what You want)
C  (if 'FZIN' is selected NEED_ANA_INI is set to TRUE)
C
C         IOMODE 'NO  '

C========= Event generator steering ======


C  CMS energy (GeV)
C
C          CMS_ENE 92.0


C  JETSET steering arrays :
C
C    MSTU
C    MSTJ
C    PARU
C    PARJ
C    PMAS
C    MDME

C  PYTHIA steering arrays :
C
C    MSTP
C    PARP
C    MSTI
C    PARI

C  more PYTHIA steering
C    MSEL
C    MSUB
C    CKIN
C    KFIN

C     (Note that KFIN is declared as KFIN(2,-40:40), PMAS as PMAS(4,500),
C     and MDME as MDME(2,2000), so You must do
C     some algebra to find the corresponding number when in a 1-dim
C     array)

C Average interaction point
C
C MEAN_VERTEX 0.0 0.0 0.0

C Beam-spot size
C
C VERTEX_SPREAD 0.015 0.0010 1.000

C Simulate interaction point ? (by default it will be a 0,0,0)
C
C PRIMARY_VERTEX_SIM FALSE


C      The following is example of how to set up a SUSYGEN
C      run. Refere to the SUSYGEN documentation for all the
C      possible settings and their defaults!


C---------- MSSM INPUTS
  MODES 5
  M 0.
  MU -100.
  M0 1000.
  TANB 2.
  A 0.
  MA 50.


C ---- ARE RELEVANT FOR MODES 2

MSQUARK 1000.
MLSTOP 1000.
MRSTOP 1000.
MLSEL 1000.
MRSEL 1000.
 MLSEL 40.0
 MRSEL 40.0
MSNU 1000.
MIX 0
C-- MAS1 0. 0. 60.5
PHIMX 0. 0. 0.

C ---- WHAT PROCESSES TO GENERATE

ZINO FALSE
WINO FALSE
SELECTRON FALSE
SMUON FALSE
STAU TRUE
SNU FALSE
SQUARK FALSE
SSTOP FALSE
SBOTTOM FALSE
HIGGS FALSE
C HIGSEL 0 0 0 0 0

C    Note that SGV needs to have the PROCESL array filled in order to
C    ramdomly choos between processes. If only the card STAU TRUE is
C    used, PROCSEL will not be filled. So, we must also select the
C    processes (stau_R and Stau_L) explicitly by :

PROCSEL 2 22 23

C---- R Parity switches
RPARITY TRUE
INDIC 3 1 2 3
LAMDA 0.1

C DECSEL 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0
C---- INCLUDE OR NOT RADIATIVE CORRECTIONS
ISR 1

DEBUG 0

C  Do parameter scan or not
SCAN FALSE
VSCAN 2 -200. 200. 2 0. 200.

C------- RUN CONDITIONS

C NB !!!! The following not used within SGV! Use the corresponding
C SGV steering card instead (CMS_ENE)
C --- ECM 183
C---    GENERATE EVENTS OR NOT (ignored by SGV: You always want to
C       generate evets in SGV.
GENER 1

C--- write out events to for012.dat or not
LUWRIT FALSE
LEPI FALSE


C========= Detector simulation steering ======

C  Generate hits or track parameters at perigee
C
C         VDHITS FALSE

C  If hits to be generated : In how many layers (not counting the beam-pipe)
C
C         VDLAYS 1

C  Make full set of track parameters with Error matrix, or just
C  3-momentum and production vertex ? If You write the event
C  to a file, putting MTKR to FALSE might be interesting: The size
C  of the file is much reduced. The parameters can the be reconstructed
C  by the analysing code. If You're not writing events out, putting
C  MTKR = false is just a waste of CPU-cycles (the same thing will be
C  done both in the simulation and the analysis)
C
C         MTKR TRUE

C  Number of detectors (MAX=3)
C
C         NDETS  1


C  Generate showers in calorimeters also for charged
C  particles (else only for neutrals).
C
C         CHSHOW TRUE

C          Generate brems and photon conversions in the detector
C          material
C        GENBC  FALSE


C          Minimum electron momentum to generate brems.
C        PMINBR  0.3

C          Minimum photon momentum to generate pair-production.
C        PMINPA  0.0

C          Minimum fraction of pt after to pt before the brems for the
C          original electron to be kept for the tracking
C        PTLOSLIM  0.9

C          Send particles with these codes to analysis
C          By default, all particles with non-zero lifetime
C          are sent and need not be specified here.
C          Use the LUND partcle codes.

C          Here, as an example, we want to save staus and Chi^0_1

       SAVE_PARTICLES  55 -55 59 -59 65 -65 68 -68  71 -71

C          Print the geometry of the dectector after loading.
C         PRDET FALSE

C          Print the geometry of the dectector after loading, in a
C          way usefull for displaying the detector with SHOWDET.KUMAC
C          (if true, will also set PRDET to TRUE)
C         PLDET FALSE

C========= Analysis steering ======

C  Make Event ntuple ? (no-op with the default ZAUSER)
C
C         MEVNT FALSE

C  Make Jet ntuple ? (no-op with the default ZAUSER)
C
C         MJETNT FALSE

C  How to analyse the VD hits ? (no-op with the default ZAUSER)
C
C         VDMET  1

C  Minimum number of jets to accept the event.
C
C         MINJET  0

C  Make primary vertex by adding tracks until the Xi^2 of the vertex
C  gets bad (DOWNUP TRUE) or by by removing tracks until it gets good
C  (DOWNUP = FALSE)
C
C         DOWNUP FALSE


C  Limiting probability for two showers in the calorimeter to be
C  separated.
C  If the separartion of the two showers is such that the
C  probability >= SHOW_SEP_LIM that two independent meassuerments
C  of the SAME shower would give that same observed separartion (or
C  less), the showers are merged. Hence, big SHOW_SEP_LIM -> many
C  showers merged, and v.v. In particular: SHOW_SEP_LIM=0. means
C  no showers are merged, and the cluster-information will essentially
C  be a copy of the calorimeter information in ZATRS.
C  NB. This parameter only takes effect if ZAUSHO is called as in
C  the example ZAUSER!
C
C         SHOW_SEP_LIM 0.9995

C  Lowest distance between shower start-point and track extrapolation
C  to the calorimeter at which the shower will not be attached
C  to the track (in cm).
C  NB. This parameter only takes effect if ZAUSHO is called as in
C  the example ZAUSER!
C
C         MIN_SEP_CLU_TRK 5.0



Mikael Berggren 2003-03-27