MadDM

newton1 crash

Asked by Lianyou SHAN

Dear Experts,

Please pardon a naive question.

I saw job failure like below when launching my calculation :

  File "/home/shanly/newPrime/MG5aMCv320/madgraph/various/lhe_parser.py", line 866, in get_alphas
    self.alpsrunner = Alphas_Runner(asmz, nloop, zmass, cmass, bmass)
  File "/home/shanly/newPrime/MG5aMCv320/models/model_reader.py", line 321, in __init__
    self.amb = self.newton1(t, asmz, 5)
  File "/home/shanly/newPrime/MG5aMCv320/models/model_reader.py", line 375, in newton1
    a_out = alphas/(1+b0*alphas*t+c1*alphas*math.log(1+alphas*b0*t))

The failure insist even I tried

  False = running_as

in my Cards/maddm_card.dat

In principle the beta functions should not be disturbed since my model change nothing with color.

Thanks a lot in advance !

   Lianyou

Question information

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

What is the Z, B, C mass for your model?

Revision history for this message
Lianyou SHAN (lyshan) said : 		#2

Dear Olivier,

Thanks for attention. Below is my setup :

BLOCK MASS #
............
      3 9.500000e-02 # ms
      4 1.275000e+00 # mc
      5 4.180000e+00 # mb
      6 1.730000e+02 # mt
............
      23 9.118760e+01 # z : mz

in Cards/param_card.dat .

Best !

   Lianyou

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#3

so those sounds reasonable.
What is the value of a_S in the param_card?

Revision history for this message
Lianyou SHAN (lyshan) said : 		#4

Dear Olivier,

In the card I have "as" as below:

BLOCK SMINPUTS #
      1 7.812500e-03 # aew
      2 1.166379e-05 # gf
      3 9.118760e+01 # mz (note that parameter not used if you use a pdf set)
      4 1.184000e-01 # as

Best !

    Lianyou

Revision history for this message
Lianyou SHAN (lyshan) said : 		#5

Dear Olivier,

I guess the

BlockOrder = 4

 was wronly used in my FR

I am changing back to 3 and retry. Thanks !

Lianyou

Revision history for this message
Lianyou SHAN (lyshan) said : 		#6

Dear Olivier,

Thanking to your suggestions to look into, the "newton1" issue has been fixed.

Best ! Lianyou