How change the mass of charm quark?

Asked by Qiji Xin

Hi,

A quick question. I want to change the mass of charm quark from 1.27 GeV (default value in MadGraph) to 1.5 GeV to match what I am using in Pythia. I am going to edit the restrict_xxx.dat files to realize this.
So is the second line below the only place I need to change? Anywhere else do I need to change?
BLOCK MASS #
      4 1.270000e+00 # mc
      5 4.700000e+00 # mb
      6 1.720000e+02 # mt

Below are the details of my restrict_xxx.dat file:

######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      4 1.270000e+00 # mc
      5 4.700000e+00 # mb
      6 1.720000e+02 # mt
      11 5.110000e-04 # me
      13 1.056600e-01 # mm
      15 1.777000e+00 # mta
      23 9.118800e+01 # mz
      25 1.250000e+02 # mh
      1 1.000000e-105 # d : 0.0
      2 1.000000e-105 # u : 0.0
      3 1.000000e-105 # s : 0.0
      12 1.000000e-105 # ve : 0.0
      14 1.000000e-105 # vm : 0.0
      16 1.000000e-105 # vt : 0.0
      21 1.000000e-105 # g : 0.0
      22 1.000000e-105 # a : 0.0
      24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.325070e+02 # aewm1
      2 1.166390e-05 # gf
      3 1.180000e-01 # as (note that parameter not used if you use a pdf set)
###################################
## INFORMATION FOR WOLFENSTEIN
###################################
BLOCK WOLFENSTEIN #
      1 0.000000e+00 # lamws
      2 0.000000e+00 # aws
      3 0.000000e+00 # rhows
      4 0.000000e+00 # etaws
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      4 0.000000e+00 # ymc
      5 4.200000e+00 # ymb
      6 1.645000e+02 # ymt
      11 0.000000e+00 # yme
      13 0.000000e+00 # ymm
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00 # wt
DECAY 15 2.270000e-12 # wtau
DECAY 23 2.441404e+00 # wz
DECAY 24 2.047600e+00 # ww
DECAY 25 6.382339e-03 # wh
DECAY 1 1.000000e-105 # d : 0.0
DECAY 2 1.000000e-105 # u : 0.0
DECAY 3 1.000000e-105 # s : 0.0
DECAY 4 1.000000e-105 # c : 0.0
DECAY 5 1.000000e-105 # b : 0.0
DECAY 11 1.000000e-105 # e- : 0.0
DECAY 12 1.000000e-105 # ve : 0.0
DECAY 13 1.000000e-105 # mu- : 0.0
DECAY 14 1.000000e-105 # vm : 0.0
DECAY 16 1.000000e-105 # vt : 0.0
DECAY 21 1.000000e-105 # g : 0.0
DECAY 22 1.000000e-105 # a : 0.0

Thank you very much!

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

If you have a param_xxx.dat in your directory, this will become the default value when doing
import model MODELNAME-xxx

So if that file exists, you will need to edit it as well.

Cheers,

Olivier

PS: The best is oviously to check the param_card in the generated directory to check which value is set there (this is the one that is going to be used)

> On 4 May 2022, at 03:11, Qiji Xin <email address hidden> wrote:
>
> New question #701658 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/701658
>
> Hi,
>
> A quick question. I want to change the mass of charm quark from 1.27 GeV (default value in MadGraph) to 1.5 GeV to match what I am using in Pythia. I am going to edit the restrict_xxx.dat files to realize this.
> So is the second line below the only place I need to change? Anywhere else do I need to change?
> BLOCK MASS #
> 4 1.270000e+00 # mc
> 5 4.700000e+00 # mb
> 6 1.720000e+02 # mt
>
> Below are the details of my restrict_xxx.dat file:
>
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 4 1.270000e+00 # mc
> 5 4.700000e+00 # mb
> 6 1.720000e+02 # mt
> 11 5.110000e-04 # me
> 13 1.056600e-01 # mm
> 15 1.777000e+00 # mta
> 23 9.118800e+01 # mz
> 25 1.250000e+02 # mh
> 1 1.000000e-105 # d : 0.0
> 2 1.000000e-105 # u : 0.0
> 3 1.000000e-105 # s : 0.0
> 12 1.000000e-105 # ve : 0.0
> 14 1.000000e-105 # vm : 0.0
> 16 1.000000e-105 # vt : 0.0
> 21 1.000000e-105 # g : 0.0
> 22 1.000000e-105 # a : 0.0
> 24 8.041900e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.325070e+02 # aewm1
> 2 1.166390e-05 # gf
> 3 1.180000e-01 # as (note that parameter not used if you use a pdf set)
> ###################################
> ## INFORMATION FOR WOLFENSTEIN
> ###################################
> BLOCK WOLFENSTEIN #
> 1 0.000000e+00 # lamws
> 2 0.000000e+00 # aws
> 3 0.000000e+00 # rhows
> 4 0.000000e+00 # etaws
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 4 0.000000e+00 # ymc
> 5 4.200000e+00 # ymb
> 6 1.645000e+02 # ymt
> 11 0.000000e+00 # yme
> 13 0.000000e+00 # ymm
> 15 1.777000e+00 # ymtau
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.508336e+00 # wt
> DECAY 15 2.270000e-12 # wtau
> DECAY 23 2.441404e+00 # wz
> DECAY 24 2.047600e+00 # ww
> DECAY 25 6.382339e-03 # wh
> DECAY 1 1.000000e-105 # d : 0.0
> DECAY 2 1.000000e-105 # u : 0.0
> DECAY 3 1.000000e-105 # s : 0.0
> DECAY 4 1.000000e-105 # c : 0.0
> DECAY 5 1.000000e-105 # b : 0.0
> DECAY 11 1.000000e-105 # e- : 0.0
> DECAY 12 1.000000e-105 # ve : 0.0
> DECAY 13 1.000000e-105 # mu- : 0.0
> DECAY 14 1.000000e-105 # vm : 0.0
> DECAY 16 1.000000e-105 # vt : 0.0
> DECAY 21 1.000000e-105 # g : 0.0
> DECAY 22 1.000000e-105 # a : 0.0
>
>
> Thank you very much!
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Qiji Xin (xinqijisuper) said :
#2

Thanks Olivier Mattelaer, that solved my question.