基本的にはMG5_aMC_v2_6_0以下のディレクトリで作成したファイルの中に入っているCardsというディレクトリの中に入っているdatファイルをいじれば良い。例えば限定的な領域での散乱断面積を見るためには、mYやmXの範囲(scan:[])をいじる。 以下では、DMsimp_s_spin1 model（[[http://hpx.phys.ocha.ac.jp/wiki/InternalPages/集中講義・物理学特論2017#FeynRules.2BMG5SKXUo-|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]を押すとプログラムが実行される。