changing coupling

Asked by Sarah

Good morning

I am working in 2HDMtII_NLO model and the results I got are slightly different than the results in the paper.
I am thinking to change the couplings.
can I change couplings in Madgraph of this model?

best regards

Question information

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

You can change the value of the param_card, which should impact the value of the coupling(s).

Cheers,

Olivier

Revision history for this message
Sarah (sarahwarad) said :
#2

I did, I used the same benchmark points in the paper, but I did not get the same results.

so cannot I change the couplings?

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

You can hack the MODEL directory if you want to, but then you need to be careful (especially at NLO) that you edit the coupling value in a consistent way. (especially at NLO, you need to be carefull about all counterterm, finite piece of the loop, ...)

Cheers,

Olivier

Revision history for this message
Sarah (sarahwarad) said :
#4

Thank you Olivier

Revision history for this message
Sarah (sarahwarad) said :
#5

Good morning
I have a question related to this point.

I am working on 2HDMtII_NLO model, studying the decay h3 > h1 z
( h3 = A Higgs boson, h1 = the SM Higgs)

in the coupling.py file of the model there are many couplings their values are functions. For example
GC_140 = Coupling(name = 'GC_140',
                  value = '-(cw*g1*TH2x1*TH3x3)/2. + (gw*sw*TH2x1*TH3x3)/2.',
                  order = {'QED':1})

Can I change the value to a number? I mean like this
Coupling(name = 'GC_140',
                  value = '1.36074468e-01',
                  order = {'QED':1})
or Madgraph will calculate them while it is running the process?

I will calculate the coupling values in 2HDMC code.

best regards

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

Hi,

> Can I change the value to a number? I mean like this
> Coupling(name = 'GC_140',
> value = '1.36074468e-01',
> order = {'QED':1})

This is a valid UFO syntax so you can do that.

Now It is likely not a smart thing to do
1) It is not clear to me that the value that you have to put here match the one of the 2HDMC code.
it is highly possible that the internal of the UFO model does not follow the convention of 2HDMC.
espeically since the splitting of an interaction into three piece (coupling, lorentz,color) is not unique
and you can easily move factor from one to the other.

2) you are using a NLO model, then a lot of those formula are correlated to each other. Just changing one coupling might break gauge invariance, nlo accuracy, lorentz invariance, .... (the same issue can also occur at LO but their the relation are simpler and the risk is less critical in general)

Cheers,

Olivier

PS: Obviously this will force you to regenerate the code since you change the UFO model.

> On 11 Nov 2021, at 10:30, Sarah <email address hidden> wrote:
>
> Question #699394 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/699394
>
> Sarah posted a new comment:
> Good morning
> I have a question related to this point.
>
> I am working on 2HDMtII_NLO model, studying the decay h3 > h1 z
> ( h3 = A Higgs boson, h1 = the SM Higgs)
>
> in the coupling.py file of the model there are many couplings their values are functions. For example
> GC_140 = Coupling(name = 'GC_140',
> value = '-(cw*g1*TH2x1*TH3x3)/2. + (gw*sw*TH2x1*TH3x3)/2.',
> order = {'QED':1})
>
> Can I change the value to a number? I mean like this
> Coupling(name = 'GC_140',
> value = '1.36074468e-01',
> order = {'QED':1})
> or Madgraph will calculate them while it is running the process?
>
> I will calculate the coupling values in 2HDMC code.
>
>
> best regards
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Sarah (sarahwarad) said :
#7

thank you
sorry I made a mistake, the value 1.36074468e-01 is the coupling of h3 h1 z in 2HDMC for CG_138 not for GC_140'

Can you help with this problem?

Provide an answer of your own, or ask Sarah for more information if necessary.

To post a message you must log in.