Model parameters for sepics.
Any comment may be placed before ---- line.
you must give a value for each row at a column corresponding to
  "specified"   if you give / only there, the standard is taken.
  No items  should be deleted or added.  If the column for the
  standard has only ",", no  standard value is defined.

The name field must start from the 2nd column.
Blank lines, if any,  are ignored. If the first 2 columns are blank,
the line is regarded as a comment line and ignored.



name    standard   specified    meaning    options & note
----------------------------------------------------------------------
 BaseTime     1.             / average cpu time(sec) for 1 GeV.

 CosNormal  (.6,  1.)      / The range of the cosine of the angle 
	        between the "normal vector" and the incident particle
       	     	direction.  Used only if InputA = 'is' or 'cos *'
	        The normal vector is directed to +z axis in many cases
	        but actually it depends on InputP as follows.
                InputP = 'fix':
                     +z axis   
	        InputP(2:4) = 'xy'      'xy2'  
                     +z axis            -z axis
                     In this case, Zinp is used as the incident Z-pos.
                InputP(2:4) = 'yz'      'yz2'
                     +x axis            -x axis
		     In this case, Xinp is used as the incident X-pos.
                InputP(2:3) = 'xz'      'xz2'
	             +y axis            -y axis
	             In this case, Yinp is used as the incident X-pos.
	        InputP  = 'usph'
                     +z axis
                     In this case, incident point is distributed uniformly
                     on the world sphere surface.  If the incident  particle 
		     is found not directed to the inner part of the sphere,
                     the particle direction is reversed.


                      
 DCInpX     0.d1           /  These three are x, y, z components of the 
 DCInpY     0.d1           /  direction cosine of  the incident ptcl.
 DCInpZ     1.d0,          /  Used only if InputA is 'fix'.
                              The square sum must be equal to 1. within 
                      the accuracy of 4~5 digits. Since it is bothersome
                      to give higher accuracy, sepics renormalize them.

                      The value is for the world  coordinate.

 DeadLine  '00.00.00'        /     Deadline time  (not used yet)
 E0Max     300000.         / If the sampled incident energy > E0Max, it will
                           be discarded.
 Hwhm       ,            /   Half-width at half-max of a Gaussian spread beam.
                             Used when InputP='gxy' etc. (cm) (see also ProfR)
 InputA    'is'   'cos 0'      / 'is' --> Incident particle angle is isotropic
                                   (See CosNormal).
                          'fix'--> Incident particle angle is a fixed one
                                   given by the direction cosine,
                                   (DCInpX, DCInpY, DCInpZ).
                          'cos pw'--> cos^pw dcos angular distribution
                           if pw = 0 equivalent to 'is'.  
 InputP    'fix'   'uxy'   /  Specify the incident partilce position in the
                           world coordinate.  
        'fix'--> (Xinp, Yinp, Zinp) (cm) is used.
        'uxy'-->  uniform on an x-y plane. The z-value is dependent on  Zinp.
                  z=Zinp if it is 0 (strictly, min z of the detector) or
                  z=max z of the detector if Zinp> max z.  'uxy'
                  can be used  only if
	        1) the world is given by box
                    or
                2) the first component is a box ( Zinp= min z), or the last
                   component is a box (Zinp = max z).
        'uxy2'  The same as 'uxy' but the direction of the incident is
                 affected. See CosNormal.
        'uxz'-->  uniform on an x-z plane. The y-value is dependent on Yinp as
                  z on Zinp for 'uxy'.  To use this, the world must be given
                  by box.
        'uxz2'   The same as 'uxz' but see CosNormal.
        'uyz'-->  uniform on a y-z plane. The x-value is dependent on  Xinp
                  as  z on Zinp for 'uxy'. To use this,  the world must be
                  given by box.
        'uyz2'  The same as 'uyz' but see CosNormal.
        'uxyz'--> uniform in (x,y,z) space. The world must be given by box.

        'usph'--> uniform on the world sphere.  Therefore, the world must
                  be given by a sphere.  The region of the input position 
                  can be controled by Xinp, Yinp and Zinp.
                  Xinp and Yinp are regarded as a polar angle (theta deg.)
                  and azimuthal angle (phi deg.) around which the incident
                  particles are distributed within an half opening angle of
                  Zinp (deg.).  See also CosNormal

              ****************
              Gassian beam cases follow: The Gasussian beam spread is
              controled by Hwhm (Half-width at half-max) and ProfR.
              The latter is the radius of the beam beyond which no particle
              is distributed.
        'gxy'     Gaussian beam on an x-y plane.  Otheres are the same as
                  'uxy'. However, the world must be given by box always.
        'gxy2'    See CosNormal
        'gxz'     The same as 'gxy' but the plane is on 'x-z'
        'gxy2'    See CosNormal
        'gyz'     The same as 'gxy' but the plane is on 'y-z'
        'gyz2'    See CosNormal.
        'gsph'    Gaussian beam on the world sphere. Xinp, Yinp and Zinp
                  serve as in the case of 'usph'.  ProfR is not referred
	          since Zinp can do the same. Gaussian means here that
                  the beam has a Gaussian density if it is  projected to
                  the plane which is tangent to the sphere at the  Gaussian 
                  center. See also CosNormal.
      
            If these input positions do not satisfy you, you have to make
            the incident particle in your userhook rouitne.

 Ir1(1)    123   789          / The random number generator is the same as
 Ir1(2)    50000011        / the one used in Cosmos. It can be initialized
                             by two 32-bit integers.  
 JobTime   1190           /  cpu time alloted (not used yet)
 Nevent    1    20          /   # of events to be generated.
 PrimaryFile  'primary'   / primary spectrum data file.
                            For details, see $EPiCSTOP/Data/Primary/sample.d.
                            The format is completely the same as Cosmos.
 ProfR      ,            /   Limit of gaussian beam radius (cm).
                             Used when InputP is 'gxy' etc.
 Xinp     1    14,    /   Incident pos of x. Used when InputP='fix'

 Yinp     1    14,    /   see Xinp. (see InputP in document)

 Zinp      0.          /   see Xinp
                           The values are in cm in the world coordinate.
                           These are used quite a different purpose if
                           InputP='usph'  or 'gsph'. See InputP.

 LogIr     t           / Take log of intial seed of random num. etc 
                         for each event
