c LIST - forces FFREAD to echo cards back
c useful to record conditions in log file
LIST

c
c QIE digitization ....
c
MFDI 1

TRIG 10000
FFLX '$FLUXPATH/afs_near_ph2le_v0201.ntup'
GAFN 'gm_near_whole_ph2le.xx_gaf'
HRFN 'gm_near_whole_ph2le.hr'
WGAF 0
ROCK 0

c WGAF  2
c ROCK -1

c PLNL 21 60
c FIDB 50.
c FIDZ 999.


c===============================================================
c general user section
c===============================================================
c
c RUNG - user "run" and "event" numbers
c RNDM - set random number seeds (2 values)
c TRIG - number of events to process
c STAT - if first value = 1 volume statistics are generated (GBSTAT)
c HSTA - standard GEANT histograms (list up to 10)
c        legal values: 'TIME' 'SIZE' 'MULT' 'NTRA' 'STAK'
c PRIN - call gprint() (list up to 10)
c        legal values: 'MATE' 'TMED' 'GEOM' 'VOLU' 'ROTM' 'SETS'
c SAVE - limit information into GAF (not yet implemented)
c        'INIT' 'GEOM' 'EVNT' 'NDIG' 'NHIT' ???
RUNG 102 1
C RNDM 1510229060 988762910
STAT 1
HSTA 'TIME' 'SIZE' 'MULT' 'NTRA' 'STAK'
PRIN 
c SAVE
c
c GAFF - output file format (0=ascii GAF, 1=FZ GAF)
c GAFN - output file name 
c        (if extension is .xx_gaf then "xx" is replaced by "ie" or "fz")
c HRFN - output histogram file name
c WGAF - 0 = write all events to GAF file
c        1 = only those with non-empty hit or digit tables
c        2 = only those with non-empty hit tables
c
GAFF 1
c
c MFLX - flux method 
c        0 = "far" histogram of nu_e,nu_e~,nu_mu,nu_mu~, uniform over face
c        1 = read in Mufson/Miller flux ntuple file created from Jim Hylen's
c            GBEAM output.  This works for both the near and far detectors.
c        2 = read in Peter Litchfield's flux format for cosmic neutrinos
c            (similar to 1, but energy scale is MeV)
c        3 = read in Hatcher flux format; this format avoids trig overhead
c            of method 1 and fixes some coordinate transformation errors
c            in that code.   Flux is derived form reweighting code of J. Hylen
c FFLX - flux method file name (meaningful values depends on method).  For
c        MFLX=1,2,3: input character string denoting name of ntuple file.
c BCEN beam_x0 beam_y0 beam_z0   - set beam spot center in detector co-ord
c                                  defines center of reference plane
c BSLP beam_dxdz beam_dydz       - set beam slope
c ZREF z_tgt2ref                 - desired z from target to reference plane
c FREF ref_z_nominal             - z tgt2ref used in generating flux file
c MFLV - flavor changing method
c        0 = no changes
c        1 = fixed permutations
c PFLV - flavor permutation map (MFLV=1)  6 values
c        integer value = PDG id for nu_e,nu_mu,nu_tau,nu_e~,nu_mu~,nu_tau~
c        nu_e = 12, nu_mu = 14, nu_tau = 16, anti-neutrino is negative
c MKIN - event kinematics method
c        0 = call usr_kine
c        1 = call hepevt_file_kine (reads "hepevt.dat" STDHEP file)
c        2 = call neugen_kine (integrated NEUGEN package)
c        3 = call single_part_kine (uses /GCKINE/ IKINE and PKINE)
c            IKINE = geant particle #
c            PKINE(1..3) = particle momentum
c            set using KINE card, interactively change with KINE command
c        4 = up- and down- going muons
c        5 = allocated for beam (rock) muon method  - NOT YET IMPLEMENTED
c        6 = allocated for radioactive noise method - NOT YET IMPLEMENTED
c FNDE - find event (meaningful only for MKIN=1)
c        0 = start at beginning of file
c       -n = skip abs(n) events
c       +n = look for event # "n"
c MTDK - tau decay method
c        0 = neugen+tauola (no displaced vertex)
c        1 = GEANT GDECAY (not recommended)
c        2 = LUND performs the decay (GLUDKY_MINOS)
c        3 = NEUGEN+tauola; GEANT displace vertex (zombie/reanimate)
c            by tracking the re-animated tau; substitute tauola's decay
c            products (not initially entered) when tau decays
c        4 = tauola called by GUDCAY (not yet implemented)
c ROCK - inclusion of rock in allowed interaction region
c         0 = only the detector
c        -1 = only the upstream rock
c        +1 = upstream rock + detector
c CCNC - CC or NC generation
c        -1 = NC + CC
c         1 = CC only
c         2 = NC only
c DK0F - list of IPDG variables NEUGEN should _not_ decay
c QELR - QE or DIS
c        -1 = all
c         0 = QE-like
c         1 = inelastic
c RSCT - intranuclear rescattering 0=off 1=on (use 0)
c TAUB - choose particular tau branching mode
c         0 = all modes
c         1 = tau -> e   only
c         2 = tau -> mu  only
c         3 = tau -> pi  only
c         4 = tau -> rho only
c         5 = tau -> a1  only
c         6 = tau -> K   only
c         7 = tau -> K*  only
c DISG - Deep Inelastic Scattering generator
c         0 = standard Soudan2 DIS
c         1 = PYTHIA/JETSET for DIS - allows charm production
c HQDK - heavy quark decay method
c         0 = NEUGEN/PYTHIA promptly decays heavy Q states (charm hadrons)
c         1 = GEANT propagates heavy Q hadrons, LUND handles decay kinematics
c ACPT - acceptance cuts on NEUGEN events (see kine/accept_neugen_kine)
c         0 = keep all legal NEUGEN events
c         1 = require event to contain heavy quark (charm or bottom)
c SAVT - list of event # (idevt) in which to save intermediate track info.
c        Saved to StdHep's /HEPEVT/ (limited to 4000 lines).  List can
c        contain 100 events; pairs so negative values denote a range.
c C2ND - threshold (GeV) for saving secondary particle to /HEPEVT/;
c        one threshold per GEANT process.  If the value is positive
c        then any particle with PTOT greater than threshold saves all the
c        secondaries in that interaction; if negative then only those 
c        particles over threshold are recorded.
c        By repeating the card one can overwrite individual values
c        after setting an initial default.  ie.
c            C2ND 60*0.020
c            C2ND 9=0.150 10=0.150 11=0.150
c        sets a limit of 20MeV for all processes except BREM,DRAY,ANNI
c        which use 150MeV.  (default is 0.150).
c        Process are:
c           1='NEXT'  2='MULS'  3='LOSS'  4='FIEL'  5='DCAY'
c           6='PAIR'  7='COMP'  8='PHOT'  9='BREM' 10='DRAY'
c          11='ANNI' 12='HADR' 13='ECOH' 14='EVAP' 15='FISS'
c          16='ABSO' 17='ANNH' 18='CAPT' 19='EINC' 20='INHE'
c          21='MUNU' 22='TOFM' 23='PFIS' 24='SCUT' 25='RAYL'
c          26='PARA' 27='PRED' 28='LOOP' 29='NULL' 30='STOP'
c          31='LABS' 32='LREF' 33='SMAX' 34='SCOR' 35='CKOV'
c          36='REFL' 37='REFR' 38='SYNC' 39='STRA'
c
c ATMD - (MKIN=4 only) direction of muons
c         0 = up
c         non-zero = down
c ATMF - up muon spectrum histogram file
c ADMF - down muon spectrum histogram file
c
MFLX 3

MFLV 1
PFLV 12 14 16 -12 -14 -16
MKIN 2
c FNDE 10
MTDK 3
RSCT 0
c ROCK 0
c BCEN +131.06 52.56
c BCEN +131.06 326.05

c
c restrict the generated vertex position
c     PLNL  - plane limits first & last
c     FIDB  - max radius around the beam axis
c     FIDZ  - max radius around the detector z axis (ie. coil)
c

c
c SXYZ - save step information in JXYZ structure during GUSTEP call
c        0 = nothing stored (use this during batch runs)
c       -x = store only charged particles
c        1 = store only while entering active detector volume
c        2 = store when entering any volume
SXYZ 0
c
c===============================================================
c detector geometry section
c===============================================================
c
c SUPR - overall supermodule organization
c        -1 = near LST detector   +1 = far LST detector
c        -2 = near FLS detector   +2 = far FLS detector
c        -3 = near RPC detector   +3 = far RPC detector
c NORF - override Near/Far as set by SUPR, leave geometry unaffected
c        (near=negative values, far=positive value)
SUPR -2
C NORF 2
c
c VHAL - detector hall volume size (half lengths in cm)
c VHAL VHAL 700. 600. 5250.

c VHAL  700. 700. 3784

c HXY0 0. 0.

c up to 15 independent supermodules can be defined
c NPLS - total number of active planes in each supermodule (end list w/ 0)
c SSPA - additional space up(down)-stream of the first(last) plane 
c        of a supermodule (allows for magnet coils); 2 values/supermodule.
c        The default grants 2m extra space upstream of first module to
c        allow for "walk space".
c SX0,SY0 - x,y offsets of supermodule relative to centerline of HALL volume
c NPLS    240    240    0    0
c SSPA 250. 50. 2*50.  2*50. 

c SSPA    4907.3 0. 0.

c SX0      0.     0.    0.
c SY0      0.     0.    0.
c
c CSET - coil winding placements for supermodule
c        0 = no coils
c        1 = center bus; return bus in (-x,-y) corner
c        4 = center bus; return bus bars in 4 corners
c        8 = center bus; return bus bars every 45degrees
c       10 = dipole configuration; 4 bus bars at 80% of half-length;
c            horizontal connects at front/back faces
c CCUR - B field map multiplier for supermodule (-1=focus mu-,+1=focus mu+)
c RMAG - coil bus bar cylinder radius 
c CMAG - coil bus bar clearance
c BFLD - maximum expected field in KGauss
c CSET  15*1
CCUR  15*-1.
RMAG 15.
CMAG 22.5
BFLD 20.0
c
c planes in the supermodule are grouped as (passive,active) planes
c each supermodule is made up of identical integral "modules"
c each module can be made of up to 16 pairs
c the number of pairs set in SiPA should integrally divide NPLS value
c SiPA - (Super)module "i" (hex value 1 to F) "pair" list (up to 16 pairs)
c        terminate the list with the string '....' (can be 17th value)
c        a "pair" takes the form 'PjAk' where:
c           P - passive plane type
c               'A'=aluminum       'B'=magnetized iron
c               'C'=concrete       'G'=glass
c               'K'=skip (air)     'P'=Pb (lead)
c               'S'=steel (no B)   'Z'=user specified (cf ABSx cards)
c           j - "instance" character (1..9,A..Z)
c               distinquish between different configurations
c           A - active plane type
c                'L' = LST = limited streamer tube (also proportional tubes)
c                'F' = FLS = fiber in liq scint (also plastic scint)
c                'R' = RPC = resistive plate chambers
c                'T' = TST = test planes (placeholder for user geometry)
c           k - "instance" character (1..9,A..Z)
c        any combination of passives and actives is permissible
c        different active can be included in the same module
c SiRM - (Super)module "i" rotation matrix string (one char/pair)
c        only the active planes are rotated
c              cell direction   readout end   increasing cell #
c        'X' = horizontal           -x         (-x,-y) -> (-x,+y)
c        'Y' = vertical             +y         (-x,+y) -> (+x,+y)
c        'U' = 45deg             (-x,+y)       (-x, 0) -> ( 0,+y)
c        'V' = 45deg             (+x,+y)       ( 0,+y) -> (+x, 0)
c
c example:
C S1PA 'B4F4' 'B4F4' 'B4F4' 'B4F4' 'B4F4' 'B4F4' 'B4F4' 'B4F4' '....'
C S1RM 'XYUVXYUV'
C S2PA 'B1L2' 'B3F4' 'B5R6' 'A1T1' '....'
C S2RM 'XYXY'
c
c Volumes are labelled 'PjTT' where:
c    'Pj' is the plane type and instance (either passive or active)
c    'TT' is one of:
c       'PL'  plane      clipping volume
c       'BX'  box        artifact for extrusion division (no user adjustment)
c       'XT'  extrusion  extrusion volume divided out of BX volume
c       'CB'  comb       forms cell walls within extrusion
c       'CV'  cover      forms cell top
c       'CL'  cell       cell volume displaces comb, touches cover
c       'GV'  gap CV     dead region between cover and extrusion
c       'GB'  gap CB     dead region between comb and extrusion
c       passive planes only have 'PL' specifications
c Each volume has a number of characteristics
c    name             type       volume  description/legal values
c    'SHAP' shape     hollerith  '??PL'  'BOX' 'TUBE' 'PGON' (PGON=octogon)
c    'MEDI' medium    hollerith   any    first 4 char of tracking medium
c    'WIDT' width     real        any    transverse size of volume
c    'THIC' thickness real        any    depth of volume in z
c    'AIR1' airgap(1) real       '??PL'  air space on front side of plane
c    'AIR2' airgap(2) real       '??PL'  extra space on back side of plane
c    'BMAP' bmap      integer    'B?PL'  magnetic field map (0=none)
c    'CHIL' children  integer     ----   # of sub-volumes
c                                '??XT'  1=comb alone,2=CB,CV,3+=CB,CV,GB,GV
c    'DIGT' digittype integer    '??PL'  digitization scheme
c                                        1=LST,2=FLS,3=RPC
c Characteristics are set as triplets (volume name,attribute,value)
c Real valued attributes are set using the RVOL,RTAG,RVAL cards
c Integer/Hollerith attributes use IVOL,ITAG,IVAL cards.
c Up to 1500 triplets of each type can be set.
c An instance value of '*' sets the default.
c Note: re-issuing the same card name restarts the sequence and potentially 
c overwrites values.  Users must explicitly give an offset.  Alternatively,
c values can be continued by placing them on lines immediately following
c the named card with no intervening keys.
c example:
C RVOL 'F4PL' 'F4XT' 'F4XT' 'F4CB' 'F4CB' 'F4CL' 'F4CL'
C RTAG 'THIC' 'THIC' 'WIDT' 'THIC' 'WIDT' 'THIC' 'WIDT'
C RVAL   2.0    2.0   66.5    1.7   66.4    1.7   1.971875 
C IVOL 'F4CB' 'F4XT'
C ITAG 'CHIL' 'CHIL'
C IVAL    32      1
C RVOL 10='R2CL' 'R2CB' 'R2GV' 'R2GV' 'R2GB' 'R2GB'
C RTAG 10='THIC' 'THIC' 'THIC' 'WIDT' 'THIC' 'WIDT'
C RVAL 10=0.2      0.3    0.2    7.8    0.5    8.0
C IVOL 10='R2GV' 'R2GB' 'R2XT' 'R2GV' 'R2GB' 
C ITAG 10='MEDI' 'MEDI' 'CHIL' 'SHAP' 'SHAP' 
C IVAL 10='RPC*' 'RPC*'     4  'BOX ' 'BOX ' 
c RVOL 'F4PL'
c RTAG 'AIR2'
c RVAL 0.4
c
c Of special interest to the "average" user might be how one modifies
c the B field map used by that plane.  All "B"-type planes of the
c same "instance" use the same field map; if you need to mix different
c maps in the same setup you'll need to construct a duplicate physical
c configuration.  The `?' gets replaced by the instance value, and
c the "bmap#" must be a positive integer.  Map 0 is the z-independent,
c azimuthally symmetric, phi-only field.  Additional maps should be
c added interfaced through the GUFLD routine (not replace it).
C IVOL 'B?PL'
C ITAG 'BMAP'
C IVAL  bmap#
c
c===============================================================
c detector materials/tracking section
c===============================================================
c AUTO - override automagic computation of tracking parameters 
c         0 = use manual tracking limits (see LBL_MTRK)
c         1 = use normal computation
c SLST,SFLS,SRPC,SEMU,SFE - tracking cuts for active detector volumes (and Fe)
c  these do not take effect unless AUTO = 0
c  negative values for 1,2,3,5 are still computed using auto scheme
c  unless explicitly set 6,7,8,9 are taken from general GEANT physics cards
c    SXXX( 1) = tmaxfd - max ang deviation due to B field (in degrees)
c    SXXX( 2) = stemax - max step permitted (cm)
c    SXXX( 3) = deemax - max fractional energy loss in 1 step [0..1.]
c    SXXX( 4) = epsil  - boundary crossing precision (cm)
c    SXXX( 5) = stmin  - minimum step size from (e-loss,scattering,Bfield)
c    SXXX( 6) = idray  - delta ray mechanism (see PHYS001)
c    SXXX( 7) = iloss  - continuous energy loss mechanism
c    SXXX( 8) = dcute  - cutoff between continuous & delta-ray by e
c    SXXX( 9) = dcutm  -   "      "       "            "       by mu
c CLST,CFLS,CRPC,CEMU,CFE - cutoff values for active volumes (and iron)
c    cutgam cutele cutneu cuthad cutmuo bcute bcutm dcute dcutm ppcutm
c FLST,FFLS,FRPC,FEMU,FFE - physics flags for active volumes (and iron)
c    pair comp phot pfis dray anni brem hadr munu dcay loss muls stra 
C SXXX tmaxfd stemax deemax epsil stmin idray iloss dcute dcutm
C CXXX cutgam cutele cutneu cuthad cutmuo bcute bcutm dcute dcutm ppcutm
C PXXX pair comp phot pfis dray anni brem hadr munu dcay loss muls stra
c LSTD,RPCD - gas density of LST or RPC mixture
c LSTW,RPCW - gas mixture fractions (by volume)
c  XXXW(1) = isobutane
c  XXXW(2) = CO2
c  XXXW(3) = argon
c  XXXW(4) = SF6
c  XXXW(5) = Freon13b1
c MXST - maximum steps GEANT will track a particle
C AUTO 0
C SLST -5.0 -0.01 -0.01 0.005 -0.005 DRAY LOSS CUTS(8) CUTS(9)
C SFLS -5.0 -0.01 -0.01 0.005 -0.005 DRAY LOSS CUTS(8) CUTS(9)
C SRPC -5.0 -0.01 -0.01 0.005 -0.005 DRAY LOSS CUTS(8) CUTS(9)
C LSTD 1.205e-3
C LSTW .090 .700 .210 .000 .000 
C RPCD 1.205e-3
C RPCD .300 .000 .700 .000 .000
SFE 4=0.001
MXST 200000
c
c Allow the user to create their on passive absorber material
c ABSN - mixture name (*must* be 20 char or less, FFREAD apparently doesn't
c        catch overlong strings; on some machines this causes core dumps)
c ABSA - list of element's A values (up to 10, end list w/ 0)
c ABSZ - list of element's Z values (up to 10, end list w/ 0)
c ABSW - list of element's fraction by weight (up to 10)
c ABSD - mixture density
c     12345678901234567890
ABSN 'LUCITE CH2CCH3CO2CH3'
ABSA 1.   12.   16.    0.
ABSZ 1.    6.    8.    0.
ABSW  .08   .60   .32  0.
ABSD 1.16
c
c LIQD - liquid scintillator density (Bicron sez: .86, MACRO sez: .8222)
LIQD 0.86
c===============================================================
c active detector section
c===============================================================
c
c FATL - fiber attenuation lengths (short+long) in cm
c FATN - fiber attenuation normalizations (short+long)
c STUB - fiber pig-tail beyond extrusion end
c PEMU - photo electrons in normalization calculation
c ELMU - dedx energy loss for muon in normalization calculation
c FATL 135. 1100.
c FATN 1.0  0.66
c STUB 100.
c PEMU 15.6
c ELMU 0.0028
c
c===============================================================
c standard GEANT physics data card section
c===============================================================
c
c  Turn on various PHYSICS processes.
c
c -- avoid double counting of delta-rays and energy loss
c -- restricted fluctuations DRAY=1,LOSS=1
DRAY 1
LOSS 1
c -- no explicit delta rays DRAY=0,LOSS=2
c DRAY 0
c LOSS 2
c -- no fluctuations DRAY=1,LOSS=4
c DRAY 1
c LOSS 4
ANNI 1
BREM 1
COMP 1
DCAY 1
HADR 1
MULS 1
MUNU 1
PAIR 1
PFIS 1
PHOT 1
c
c 100KeV cuts 
c
CUTS 10*.0001
c
c===============================================================
c signal end off FFREAD input
c===============================================================
END