Issue with compute_widths not updating parameter values in param_card.dat

Hi,

Concerning:
./bin/mg5_aMC
MG5_aMC> import model ScotoSinglet_Complex_UFO/
MG5_aMC> compute_widths h

In this case they are no param_card.dat associated to the model, so it makes sense that it is using the default.
if a file restrict_default.dat exists in the UFO model, this is the file that we will use for the default BUT if a file param_default.dat also exists.
If restrict_default.dat does not exists then we will use the hardcoded value from the model.

So here you do not face an unexpected behavior.

So as you do at the end, you need to specify the param_card that you want to use with
compute_widths h --path=/home/sebastian/Documents/scoto_hTOmumu_and_g2/MadWidth_Computations/Non-Processed_Param_Cards/param_card_1.dat

The error:
InvalidParamCard : frblock (8,) is already define to 0.1172 impossible to assign 495688.0509465932.

is typically happening when a syntax error is present in the param_card and that lead to a double definition of a parameter.

Now I see that the name of that parameter is "aS" which needs to be handle in a special way within the code (due to it's running).
In principle that parameter should be in the SMinputs block (entry #3).
Not sure if this is the reason of your issue, but in any case, I would advocate to follow the standard convention for the meaning of the parameter
(and therefore remove entry 1-8 from the frblock of your model).

If you send me the model by email (<email address hidden> <mailto:<email address hidden>>) and the param_cad (and a link to this page) then i will be able to take a deeper look.

Olivier

> On 7 Nov 2024, at 00:40, Sebastián Ariel Norero Cárdenas <email address hidden> wrote:
>
> New question #819285 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/819285
>
> I am currently using MadGraph 5 (v3.4.0) with a custom model defined via FeynRules, named ScotoSinglet_Complex_UFO. I've encountered a problem where updates to the param_card.dat are not being recognized by the compute_widths command.
>
> Next are two steps to reproduce my issue:
>
> 1.- Here is my original param_card.dat file:
>
> ######################################################################
> ## 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 FRBLOCK
> ###################################
> Block frblock
> 1 7.757952e-03 # aEWM1
> 2 1.166370e-05 # Gf
> 3 9.118760e+01 # MMZ
> 4 7.994700e+01 # MMW
> 5 1.200000e+00 # MMC
> 6 4.230000e+00 # MMB
> 7 1.750000e+02 # MMT
> 8 1.172000e-01 # aS
> 9 1.000000e+02 # QS
> 10 2.532000e+00 # lam2
> 11 2.054300e+00 # lam3
> 12 2.050000e-10 # lam4
> 13 7.950000e-08 # lam5
> 14 1.180000e-02 # lamphi
> 15 -1.216000e-07 # lametaphi
> 16 7.950000e-08 # lamHphi
> 17 1.713900e+00 # muparam
> 18 8.235480e+00 # mphiSQ
> 19 -2.184172e+00 # metaSQ
> 20 9.982780e-01 # sint
> 21 5.866680e-02 # cost
>
> ###################################
> ## INFORMATION FOR FRBLOCK10
> ###################################
> Block frblock10
> 1 1 7.892700e-12 # varthetaR1x1
> 1 2 5.974833e-12 # varthetaR1x2
> 1 3 1.416562e-12 # varthetaR1x3
> 2 1 4.727221e-12 # varthetaR2x1
> 2 2 5.339249e-12 # varthetaR2x2
> 2 3 7.010406e-12 # varthetaR2x3
> 3 1 3.919012e-12 # varthetaR3x1
> 3 2 5.982790e-12 # varthetaR3x2
> 3 3 6.989103e-12 # varthetaR3x3
>
> ###################################
> ## INFORMATION FOR FRBLOCK6
> ###################################
> Block frblock6
> 1 1 3.590000e-03 # yv1x1
> 1 2 2.092000e-02 # yv1x2
> 1 3 2.156000e-02 # yv1x3
> 2 1 1.140000e-02 # yv2x1
> 2 2 1.350000e-02 # yv2x2
> 2 3 -1.002000e-02 # yv2x3
> 3 1 -1.656400e-02 # yv3x1
> 3 2 -6.970000e-03 # yv3x2
> 3 3 1.029000e-02 # yv3x3
>
> ###################################
> ## INFORMATION FOR FRBLOCK7
> ###################################
> Block frblock7
> 1 1 -1.946180e-05 # kappa1x1
> 1 2 1.196800e-09 # kappa1x2
> 1 3 4.023900e-11 # kappa1x3
> 2 1 5.464252e-01 # kappa2x1
> 2 2 -3.893000e-10 # kappa2x2
> 2 3 -8.847189e-06 # kappa2x3
> 3 1 -2.915510e-04 # kappa3x1
> 3 2 3.464147e+00 # kappa3x2
> 3 3 6.418249e-05 # kappa3x3
>
> ###################################
> ## INFORMATION FOR FRBLOCK8
> ###################################
> Block frblock8
> 1 1 7.892700e-01 # PMNS1x1
> 1 2 5.974833e-01 # PMNS1x2
> 1 3 1.416562e-01 # PMNS1x3
> 2 1 4.727221e-01 # PMNS2x1
> 2 2 5.339249e-01 # PMNS2x2
> 2 3 7.010406e-01 # PMNS2x3
> 3 1 3.919012e-01 # PMNS3x1
> 3 2 5.982790e-01 # PMNS3x2
> 3 3 6.989103e-01 # PMNS3x3
>
> ###################################
> ## INFORMATION FOR FRBLOCK9
> ###################################
> Block frblock9
> 1 1 7.892700e-12 # varthetaL1x1
> 1 2 5.974833e-12 # varthetaL1x2
> 1 3 1.416562e-12 # varthetaL1x3
> 2 1 4.727221e-12 # varthetaL2x1
> 2 2 5.339249e-12 # varthetaL2x2
> 2 3 7.010406e-12 # varthetaL2x3
> 3 1 3.919012e-12 # varthetaL3x1
> 3 2 5.982790e-12 # varthetaL3x2
> 3 3 6.989103e-12 # varthetaL3x3
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 1 0.000000e+00 # MD
> 2 0.000000e+00 # MU
> 3 1.000000e-01 # MS
> 11 5.110000e-04 # Me
> 12 2.856963e-11 # Mve
> 13 1.057000e-01 # MM
> 14 2.978972e-11 # Mvm
> 15 1.777000e+00 # MTA
> 16 5.833830e-11 # Mvt
> 25 1.252000e+02 # mh
> 312000 1.077820e+03 # MN1
> 314000 4.635148e+02 # MN2
> 316000 6.511449e+02 # MN3
> ## 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.
> 4 1.200000e+00 # c : MMC
> 5 4.230000e+00 # b : MMB
> 6 1.750000e+02 # t : MMT
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 23 9.118760e+01 # Z : MMZ
> 24 7.994700e+01 # W+ : MMW
> 3500 2.495363e+02 # etaR : cmath.sqrt(metaSQ + 0.5*(lam3 + lam4 + lam5)*v__exp__2)
> 3600 2.495362e+02 # etaI : cmath.sqrt(metaSQ + 0.5*(lam3 + lam4 - lam5)*v__exp__2)
> 4500 2.609172e+00 # S1p : cmath.sqrt(M11SQ + M22SQ - cmath.sqrt(4*M12SQ__exp__2 + (M11SQ - M22SQ)**2))/sqrt__2
> 4600 2.495391e+02 # S2p : cmath.sqrt(M11SQ + M22SQ + cmath.sqrt(4*M12SQ__exp__2 + (M11SQ - M22SQ)**2))/sqrt__2
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
>
> Block QNUMBERS 312000 # N1
> 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 314000 # N2
> 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 316000 # N3
> 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 3500 # etaR
> 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 3600 # etaI
> 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 4500 # S1p
> 1 3 # 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 4600 # S2p
> 1 3 # 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)
> #
> #*************************
> # Decay widths *
> #*************************
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 2.000000e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.495200e+00
> #
> # PDG Width
> DECAY 24 2.085000e+00
> #
> # PDG Width
> DECAY 25 0.000000e+00
> #
> # PDG Width
> DECAY 3500 0.000000e+00
> #
> # PDG Width
> DECAY 3600 0.000000e+00
> #
> # PDG Width
> DECAY 4500 1.000000e+00
> #
> # PDG Width
> DECAY 4600 1.000000e+00
> #
> # PDG Width
> DECAY 312000 0.000000e+00
> #
> # PDG Width
> DECAY 314000 0.000000e+00
> #
> # PDG Width
> DECAY 316000 0.000000e+00
>
>
> So, now I will manually update the parameter lam2 in param_card.dat from its default value of 2.532000e+00 to, say, 1.7 (and save the file, of course)
>
> 2.- Now I run MadWidth to compute the branching ratios of the 125-Higgs (h) as follows:
>
> ./bin/mg5_aMC
> MG5_aMC> import model ScotoSinglet_Complex_UFO/
> MG5_aMC> compute_widths h
>
> However, when opening the param_card.dat file, I can notice that the value of lam2 has been reverted back to the default 2.532000e+00 and it seems like the calculations were done using this old value. What could be causing this kind of issue? By the way, I created the "ScotoSinglet_Complex_UFO" model using FeynRules and I'm unsure if there might be a setting within that process that's leading to this unexpected behavior.
>
> Additionally, when I attempted to specify a path to an updated param card using:
>
> MG5_aMC> compute_widths h --path=/home/sebastian/Documents/scoto_hTOmumu_and_g2/MadWidth_Computations/Non-Processed_Param_Cards/param_card_1.dat
>
> I ran into an error saying:
>
> InvalidParamCard : frblock (8,) is already define to 0.1172 impossible to assign 495688.0509465932.
>
> It feels like there’s a setting in the UFO model that I'm missing.
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.