電子を水に打ち込む

• 幅100cm、長さ400cmの直方体の水(cell=1)に、左から 1 GeV の電子を1発打ち込んでみた。

---

電子 (2d_xz_e-.eps)

陽電子 (2d_xz_e+.eps)

光子 (2d_xz_photon.eps)

インプットファイル

$======================================
[ Parameters ]
  maxcas = 1
  maxbch = 1

  negs = 2

$======================================
[ source ]
   s-type = 1             $ mono-energetic axial source
     proj = electron
       z0 = -100
       z1 = -100
       e0 = 1000

$======================================
[ surface ]
 1  rpp    -50 50  -50 50  0 400
 99  so    2000

$======================================
[ material ]
 mat[1] H 2 O 1

$======================================
[ cell ]
 1    1  -1.00  -1
 98   0         -99 #1
 99  -1          99

$======================================
[ t-track ]
     file = 2d_xz_e-.out
     part = electron

     mesh =  xyz

   x-type = 2
     xmin = -60
     xmax = 60
       nx = 99

   y-type = 2
     ymin = -60
     ymax = 60
       ny = 1

   z-type = 2
     zmin = -100
     zmax = 420
       nz = 200

   e-type = 2
     emin = 0
     emax = 10000000
       ne = 1

   t-type = 2
     tmin = 0
     tmax = 15
       nt = 15

     unit = 1        $ 1: [1/cm^2/source]
     axis = xz

   epsout = 1
    gshow = 3

$======================================
[ t-track ]
     file = 2d_xz_e+.out
     part = positron

     mesh =  xyz

   x-type = 2
       nx = 99
     xmin = -60
     xmax = 60

   y-type = 2
       ny = 1
     ymin = -60
     ymax = 60

   z-type = 2
       nz = 200
     zmin = -100
     zmax = 420

   e-type = 2
       ne = 1
     emin = 0
     emax = 10000000

$ 時間情報は省略可
   t-type = 2
       nt = 15
     tmin = 0
     tmax = 15

     unit = 1        $ 1: [1/cm^2/source]
     axis = xz

   epsout = 1
    gshow = 3

$======================================
[ t-track ]
     file = 2d_xz_photon.out
     part = photon

     mesh =  xyz

   x-type = 2
       nx = 99
     xmin = -60
     xmax = 60

   y-type = 2
       ny = 1
     ymin = -60
     ymax = 60

   z-type = 2
       nz = 200
     zmin = -100
     zmax = 420

   e-type = 2
       ne = 1
     emin = 0
     emax = 10000000

$ 時間情報は省略可
   t-type = 2
       nt = 15
     tmin = 0
     tmax = 15

     unit = 1        $ 1: [1/cm^2/source]
     axis = xz

   epsout = 1
    gshow = 1