基本的にはMG5_aMC_v2_6_0以下のディレクトリで作成したファイルの中に入っているCardsというディレクトリの中に入っているdatファイルをいじれば良い。例えば限定的な領域での散乱断面積を見るためには、mYやmXの範囲(scan:[])をいじる。

以下では、DMsimp_s_spin1 model(FeynRulesの利用を参照)を使用したp p > xd xd~ jというプロセスを計算するとき、mediatorとダークマターの質量を変えて計算した例を紹介する。

$ cd MG5_aMC_v2_6_0
$ ./bin/mg5_aMC
MG5_aMC>import model DMsimp_s_spin1
MG5_aMC>generate p p > xd xd~ j
MG5_aMC>output
MG5_aMC>launch

何か出たら一回[Enter key]を押すと

Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
|  1. param               : param_card.dat                   |
|  2. run                 : run_card.dat                     |
|  3. madanalysis5_parton : madanalysis5_parton_card.dat     |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madanalysis5_parton, enter path][90s to answer]
>

とでるので[1]を押すとparam_card.datが開くので日本語が書いてある部分を編集する。

 GNU nano 2.0.9        File: /data/asai/IntensiveCourses/MG5_aMC_v2_6_0/PROC_DMsimp_s_spin1_1/Cards/param_card.dat                       

######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL   ####
######################################################################
##                                                                  ##
##  Width set on Auto will be computed following the information    ##
##        present in the decay.py files of the model.               ##
##        See  arXiv:1402.1178 for more details.                    ##
##                                                                  ##
######################################################################

###################################
## INFORMATION FOR CKMBLOCK
###################################
Block ckmblock
    1 2.277360e-01 # cabi

###################################
## INFORMATION FOR DMINPUTS
###################################
Block dminputs
    1 0.000000e+00 # gVXc
    2 1.000000e+00 # gVXd
    3 0.000000e+00 # gAXd
    4 2.500000e-01 # gVd11
    5 2.500000e-01 # gVu11
    6 2.500000e-01 # gVd22
    7 2.500000e-01 # gVu22
    8 2.500000e-01 # gVd33
    9 2.500000e-01 # gVu33
   10 0.000000e+00 # gVl11
   11 0.000000e+00 # gVl22
   12 0.000000e+00 # gVl33
   13 0.000000e+00 # gAd11
   14 0.000000e+00 # gAu11
   15 0.000000e+00 # gAd22
   16 0.000000e+00 # gAu22
   17 0.000000e+00 # gAd33
   18 0.000000e+00 # gAu33
   19 0.000000e+00 # gAl11
   20 0.000000e+00 # gAl22
   21 0.000000e+00 # gAl33
   22 0.000000e+00 # gnu11
   23 0.000000e+00 # gnu22
   24 0.000000e+00 # gnu33
   25 0.000000e+00 # gVu31
   26 0.000000e+00 # gAu31
   27 0.000000e+00 # gVd31
   28 0.000000e+00 # gAd31
   29 0.000000e+00 # gVh

###################################
## INFORMATION FOR LOOP
###################################
Block loop
    1 9.118800e+01 # MU_R

###################################
## INFORMATION FOR MASS
###################################
Block mass
    6 1.720000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118760e+01 # MZ
   25 1.250000e+02 # MH
  5000001 scan:[200+200*i for i in range(0,10)] #1.000000e+03 # MY1 #200-2000[GeV]までscan
  5000511 1.000000e+01 # MXr
  5000512 1.000000e+01 # MXc
  5000521 scan:[10+10*i for i in range(0,10)] #1.000000e+01 # MXd #10-100[GeV]までscan
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000 # d : 0.0
  2 0.000000 # u : 0.0
  3 0.000000 # s : 0.0
  4 0.000000 # c : 0.0
  5 0.000000 # b : 0.0
  11 0.000000 # e- : 0.0
  12 0.000000 # ve : 0.0
  13 0.000000 # mu- : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 79.824360 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
  9000002 91.187600 # ghz : MZ
  9000003 79.824360 # ghwp : MW
  9000004 79.824360 # ghwm : MW

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.279000e+02 # aEWM1
    2 1.166370e-05 # Gf
    3 1.184000e-01 # aS

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    6 1.720000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY   6 1.508336e+00 # WT
DECAY  23 2.495200e+00 # WZ
DECAY  24 2.085000e+00 # WW
DECAY  25 4.070000e-03 # WH
DECAY 5000001 auto #1.000000e+01 # WY1 #widthの補正項を自動にするとより厳密に崩壊を再現できるらしい(?)
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY  1 0.000000 # d : 0.0
DECAY  2 0.000000 # u : 0.0
DECAY  3 0.000000 # s : 0.0
DECAY  4 0.000000 # c : 0.0
DECAY  5 0.000000 # b : 0.0
DECAY  11 0.000000 # e- : 0.0
DECAY  12 0.000000 # ve : 0.0
DECAY  13 0.000000 # mu- : 0.0
DECAY  14 0.000000 # vm : 0.0
DECAY  15 0.000000 # ta- : 0.0
DECAY  16 0.000000 # vt : 0.0
DECAY  21 0.000000 # g : 0.0
DECAY  22 0.000000 # a : 0.0
DECAY  5000511 0.000000 # xr : 0.0
DECAY  5000512 0.000000 # xc : 0.0
DECAY  5000521 0.000000 # xd : 0.0
DECAY  9000002 2.495200 # ghz : WZ
DECAY  9000003 2.085000 # ghwp : WW
DECAY  9000004 2.085000 # ghwm : WW
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 5000511  # xr
        1 0  # 3 times electric charge
        2 1  # number of spin states (2S+1)
        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0  # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 5000512  # xc
        1 0  # 3 times electric charge
        2 1  # number of spin states (2S+1)
        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1  # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 5000521  # xd
        1 0  # 3 times electric charge
        2 2  # number of spin states (2S+1)
        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1  # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 5000001  # y1
        1 0  # 3 times electric charge
        2 3  # number of spin states (2S+1)
        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0  # Particle/Antiparticle distinction (0=own anti)

編集し終わったら[Ctrl+x]、[Y]、[Enter key]で終了。また[Enter key]を押すとプログラムが実行される。

いじり方のポイント (最終更新日時 2017-12-01 03:52:50 更新者 KanaeAsai)