.. _app-comprehensive-example:

:ref:`4.12節 <comprehensive-example>` で使用するPHITS入力ファイル
--------------------------------------------------------------------------------

::

   [ T i t l e ]
   sample input for dchain: 250MeV protons incident on W/H2O/Fe rod surrounded in Be

   [ P a r a m e t e r s ]
    icntl    =           0     # (D=0) 3:ECH 5:NOR 6:SRC 7,8:GSH 11:DSH 12:DUMP
    maxcas   =        1000     # (D=10) number of particles per one batch
    maxbch   =          10     # (D=10) number of batches
    igamma   =           3     # (D=0) 0:No, 1:Old, 2:EBITEM, 3:EBITEM+Isomer
    file(1)  = c:/phits        # (D=c:/phits) PHITS install folder name
   $  required options for DCHAIN 
    jmout    =           1     # (D=0) Density echo, 0:input, 1:number density
    e-mode   =           0     # (D=0) Event generator mode is not recommended for DCHAIN calculation

   [ S o u r c e ]
      s-type =   1             # mono-energetic axial source
        proj =  proton         # kind of incident particle         
          e0 =   250.00        # energy of beam [MeV/u]
          r0 =   0.5000        # radius [cm]
          x0 =   0.0000        # (D=0.0) center position of x-axis [cm]
          y0 =   0.0000        # (D=0.0) center position of y-axis [cm]
          z0 =  -10.000        # minimum position of z-axis [cm]
          z1 =  -10.000        # maximum position of z-axis [cm]
         dir =   1.0000        # z-direction of beam [cosine]

   [ M a t e r i a l ]
   $ Tungsten, Density (g/cm3) = 19.300000 
   M1      74000  1.000000  $  W
   $ Water, Density (g/cm3) = 0.998207 ~= 1.0
   M2       1000  0.666657  $  H 
            8000  0.333343  $  O
   $ Iron, Density (g/cm3) = 7.874000 
   M3      26000  1.000000  $  Fe           
   $ Beryllium, Density (g/cm3) = 1.848000            
   M4       4000  1.000000  $  Be

   [ S u r f a c e ]
   set:c1[3]  $ thickness of W portion
   set:c2[10] $ thickness of H2O portion
   set:c3[5]  $ thickness of Fe portion
   set:c4[5]  $ thickness of Be back end
   set:c5[4]  $ radius of inner cylinder
   set:c6[8]  $ radius of outer cylinder

   10 pz 0
   11 pz c1
   12 pz c1+c2
   13 pz c1+c2+c3 
   14 pz c1+c2+c3+c4 
   20 cz c5 
   21 cz c6
   99 so  80.0

   [ C e l l ]
   1   1 -19.3    10 -11 -20    $ W 
   2   2 -1.0     11 -12 -20    $ H20
   3   3 -7.874   12 -13 -20    $ Fe 
   4   4 -1.848   (10 -13 20 -21):(13 -14 -21) $ Be 
   99  0         (-99 -10):(-99 14):(-99 21)   $ vacuum 
   999 -1         99            $ particle graveyard

   [ V o l u m e ] 
   $ required section for DCHAIN 
              reg   vol
              1   pi*c1*c5**2
              2   pi*c2*c5**2
              3   pi*c3*c5**2
              4   pi*(((c1+c2+c3+c4)*c6**2)-((c1+c2+c3)*c5**2))

   [ M a t  N a m e  C o l o r ]
      mat     name     size          color   
        0     void     1.00     lightgray   
        1    Tungsten  1.00     darkgreen   
        2    Water     1.00     pastelcyan  
        3    Iron      1.00     orange  
        4    Beryllium 1.00     pastelviolet  

   [ T - D C H A I N ]
   $  required section for DCHAIN 
       title = 250MeV protons incident on Be-encased W/H2O/Fe column
        mesh =  reg               # mesh type is region-wise
         reg =  1 2 3 4
        file = spallation.out     # file name of dchain-phits input file
     timeevo =    3               # time evolution / irradiation schedule
        6.0 m 1.0
        4.0 m 0.5
       50.0 m 0.0
     outtime =    8               # output times
        1.0 m
        3.0 m
        6.0 m
       10.0 m
       15.0 m
       20.0 m
       30.0 m
       40.0 m
   $ beam current (nA)
   set:c11[100.0]
   $ beam power (source/sec) 
   set:c12[ c11 * 1.0e-9 / (1.602177e-19) ]
      amp = c12                   # (D=1.0) Source Intensity(source/sec)

   [ T - T r a c k ]
       title = proton and neutron distributions
     mesh = xyz
   x-type = 2
     xmin = -10
     xmax =  10
       nx = 1
   z-type = 2
     zmin = -10
     zmax =  30
       nz = 200
   y-type = 2
     ymin = -20
     ymax =  20
       ny = 200
   e-type = 2
     emin = 0
     emax =  2000
       ne = 1
     axis = yz
     file = yz-track.out
     part = proton neutron
     unit = 1
   2D-type =  7   # 1:Cont, 2:Clust, 3:Color, 4:xyz, 5:mat, 6:Clust+Cont, 7:Col+Cont
   gshow =    3   # 0: no 1:bnd, 2:bnd+mat, 3:bnd+reg 4:bnd+lat
    epsout=1

   [ E n d ]