New boson

• ミューオンがターゲットと相互作用することで暗黒光子が生成され、その暗黒光子が電子陽電子ペアに崩壊する。

• μ N → μ N A’, A’ → e+ e-

---

ミューオン (2d_xz_muon.eps)

---

新粒子 (2d_xz_900000.eps)

---

電子・陽電子 (2d_xz_e.eps)

---

main.inp

[ Parameters ]
  maxcas = 1
  maxbch = 1
  negs = 2
  file(1) = /Users/sakaki/programs/phits/332

  $ Turn off default muon reactions
  imucap = 0
  imuint = 0
  imubrm = 0
  imuppd = 0

[ User Defined Interaction ]
  number = 1
$ Name                 Bias     Parameters
$                               kf      boson_type  mass  g_e   g_mu  Kmax nZ  Z_array [option]
  new_boson_radiation  1.0e+0   900000  1           50    1e-6  1e-6  1000 2   8 26    $ -1

[ User Defined Particle ]
  number = 3
$ Name              kfcode   Parameters
$                            boson_type  mass  g_e   g_mu
  new_boson         900000   1           50    1e-6  1e-6
  MuonMinus_NoDecay     13
  MuonPlus_NoDecay     -13

$ (boson_type) 1:Scalar, 2:Pseudoscalar, 3:Vector, 4:Axial-vector

$==============================================================
[ source ]
 s-type = 1             # mono-energetic axial source
   proj = muon-
     z0 = -2 $ cm
     z1 = -2 $ cm
     e0 = 1000 $ MeV

$==============================================================
[ surface ]
 1  rpp  -1 1  -1 1  0 1
 99 so   1000

$==============================================================
[ material ]
 mat[1] Fe 1  O 1

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

$==============================================================
[ forced collisions ]
 part = muon- muon+
 reg   fcl
   1   -1.0

$==============================================================
[ t-track ]
   file = 2d_xz.out
   part = muon- 900000 (electron positron)
   mesh =  xyz

 x-type =  2
     nx =  199
   xmin = -2
   xmax =  2

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

 z-type =  2
     nz =  200
   zmin = -2
   zmax =  10

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

   unit =  1            # 1: [1/cm^2/source]
   axis =  xz
 epsout =  1
  gshow =  3