Process with loop_qcd_qed_sm model

Asked by SUBHASISH BEHERA on 2019-02-18

Dear MG5 Team,

I was trying to generate process with "Z > h a" or "a > h a" by using "loop_qcd_qed_sm model".
I am having the following problem while "launch".

launch PROC_SA_loop_qcd_qed_sm_0

Do you want to edit a card (press enter to bypass editing)?

/------------------------------------------------------------\

| 1. param : param_card.dat |

| 2. MadLoopParams : MadLoopParams.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, MadLoopParams, enter path]

>

INFO: Update the dependent parameter of the param_card.dat

WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...)

Initializing process emep > ah.

root: Error while executing ./check in /home/subhasish/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_2/PROC_SA_loop_qcd_qed_sm_0/SubProcesses/P0_emep_ah

root: Failed at running the process in /home/subhasish/RESEARCH/EVENT_GENERATORS/MADGRAPH5/MG5_aMC_v2_6_2/PROC_SA_loop_qcd_qed_sm_0/SubProcesses/P0_emep_ah.

Command "launch PROC_SA_loop_qcd_qed_sm_0" interrupted with error:

NameError : global name 'MadGraph5Error' is not defined

Please report this bug on https://bugs.launchpad.net/mg5amcnlo

More information is found in 'MG5_debug'.

Please attach this file to your report.

MG5_aMC>quit

** I didn't change any of the parameters in the param or MadLoopParam files.

Thanks,
Subhasish

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Valentin Hirschi Edit question
Last query:
2019-02-21
Last reply:
2019-02-20

Hi,

None of those process ( z > h a and a > a h ) does make sense since they are not enough energy to produce the final state.
So I doubt that this is what you are trying to run.

Independently, it seems that you have a bug in ... the handling of the bug.
So we will need to debug this in two steps:
1) fix the issue with the "MadGraph5Error" not defined
2) then rerun your code and find the correct bug report

In order to move on 1).
Could you create a bug report (this is a question) and attach the MG5_debug file.
(attachement are not allowed for question)

If you could also provide full details on how to reproduce the bug, this could also help us to investigate this directly and then simplify the amount of request that we will need you to do.

Thanks,

Olivier

Hi,

Ok so do I understand correctly that the problem is actually solved in the latest version?

> ** Now I don't know how to get the cross-section out of this for onshell
> Higgs production.

This cross-section should be zero due to Furry's theorem.
This is why it is not included in this paper:
https://arxiv.org/pdf/1507.00020.pdf

Cheers,

Olivier

> On 19 Feb 2019, at 23:07, SUBHASISH BEHERA <email address hidden> wrote:
>
> Question #678635 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/678635
>
> Status: Answered => Open
>
> SUBHASISH BEHERA is still having a problem:
> Dear Olivier,
>
> Thanks for the info.
>
> I was trying to check the cross-section for the process : e- e+ > h a
> [virt=QED] /W+ W- g+ g- ghwp ghwm e-
> i.e. the process will only have t or tbar in the vertex correction.
>
> I want to check for ILC@250 GeV.
>
>
> After posting yesterday I tried with "MG5_aMC_v2_6_5" which work with the
> above process and give me the following Subprocess/result.dat as :
>
> PS 5.000000000000000E+002 0.000000000000000E+000
> 0.000000000000000E+000 5.000000000000000E+002
>
> PS 5.000000000000000E+002 0.000000000000000E+000
> 0.000000000000000E+000
> -5.000000000000000E+002
>
> PS 5.078124999999999E+002 1.091910924993979E+002
> 4.378803084023696E+002
> -1.964349154007088E+002
>
> PS 4.921875000000001E+002 -1.091910924993979E+002
> -4.378803084023696E+002 1.964349154007088E+002
>
> EXP 0
>
> BORN 0.000000000000000E+000
>
> FIN 4.134354580331755E-007
>
> 1EPS 6.390188840651834E-009
>
> 2EPS 5.257571624808558E-010
>
> ASO2PI 1.878028328484365E-002
>
> Export_Format LoopInduced
>
> RETCODE 130
>
> ACC -.1000E+01
>
> Born_kept F
>
> Loop_kept T
>
> Split_Orders_Names QED
>
> Loop_SO_Results 8
>
> SO_Loop ACC -1.0000000000000000
>
> SO_Loop FIN 4.1343545803317549E-007
>
> SO_Loop 1EPS 6.3901888406518337E-009
>
> SO_Loop 2EPS 5.2575716248085581E-010
>
>
> ** Now I don't know how to get the cross-section out of this for onshell
> Higgs production.
>
>
> Thanks,
>
> Subhasish
>
> On Wed, Feb 20, 2019 at 3:27 AM Olivier Mattelaer <
> <email address hidden>> wrote:
>
>> Your question #678635 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/678635
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Hi,
>>
>> None of those process ( z > h a and a > a h ) does make sense since they
>> are not enough energy to produce the final state.
>> So I doubt that this is what you are trying to run.
>>
>> Independently, it seems that you have a bug in ... the handling of the
>> bug.
>> So we will need to debug this in two steps:
>> 1) fix the issue with the "MadGraph5Error" not defined
>> 2) then rerun your code and find the correct bug report
>>
>> In order to move on 1).
>> Could you create a bug report (this is a question) and attach the
>> MG5_debug file.
>> (attachement are not allowed for question)
>>
>> If you could also provide full details on how to reproduce the bug, this
>> could also help us to investigate this directly and then simplify the
>> amount of request that we will need you to do.
>>
>>
>> Thanks,
>>
>> Olivier
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/mg5amcnlo/+question/678635/+confirm?answer_id=0
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/mg5amcnlo/+question/678635
>>
>> You received this question notification because you asked the question.
>>
>
>
> --
>
> Thanks,
> Subhasish
>
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Dear Olivier,

in our case its e- e+ > h a, so the possible vertices are :

             a(from e- e+ fusion) h a, Z h a : something related to h > a a etc.

So as far as I understood Furry's theorem, tells that "the vacuum expectation value of the time-ordered product of any odd number of electromagnetic fields and/or currents vanish"
or

"one-loop contribution to the photon one-point and three-point functions vanish".

But here it is not the case.

Can you comment more ?

Thanks,
Subhasish

Can you help with this problem?

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

To post a message you must log in.