Fermion flow issue in the loops

Bug #1899608 reported by Benjamin Fuks
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
MadGraph5_aMC@NLO
Undecided
Valentin Hirschi

Bug Description

Hi everyone,

We are performing some new physics calculations in a model where the SM is extended by a squark (and an invisible state that is irrelevant for our purpose). The UFO is provided in attachment. The process of interest consists in dijet production. I use MG5_aMC 2.6.8 and as a process definition, I typed

import model RPVMono_UFO_RG-paper
set complex_mass_scheme True
generate p p > j j NP=2 QCD=0 QED=0 [QCD]
output

In order to solely get QCD corrections to the new physics component, I have modified the madgraph/loop/loop_diagram_generation.py file. I have set the filtering to true (around l388)
    edit_filter_manually = True
and then included
    loop_pdgs = [abs(x) for x in diag.get_loop_lines_pdgs()]
    is_loop_gluon = (21 in loop_pdgs)
    is_loop_sigma = (9000002 in loop_pdgs)
    if is_loop_sigma and not is_loop_gluon:
        valid_diag=False
I get the needed diagrams (so that this works).

However, when launching a cross section calculation, I get a pole non-cancellation issue. The weird thing is that poles are fine for the subprocesses
  d dbar to s sbar or s sbar to d dbar
but not for the subprocess
  dbar d to s sbar
This is weird as the diagrams are exactly the same. I suspect something problematic with the fermion flow here (as we have clashing arrows diagrams). Any thoughts?

Cheers,

Benj

Revision history for this message
Benjamin Fuks (fuks) wrote :
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) wrote :

Did you check the born/loop/real matrix-element with the check function?
If this is an issue related with the fermion flow it will likely spot the issue.

Cheers,

Olivier

Revision history for this message
Benjamin Fuks (fuks) wrote : Re: [Bug 1899608] Re: Fermion flow issue in the loops

Hi Olivier,

I directly ran
   check d d~ > s s~ NP=2 QCD=0 QED=0 [QCD]
   check d~ d > s s~ NP=2 QCD=0 QED=0 [QCD]
and the gauge test failed. It worked for some of the other subprocesses. What I don't understand and how it could work as the exact same vertices/counterterms are used all over the place. Any thoughts?

Cheers,

Benj

> On 13 Oct 2020, at 22:33 , Olivier Mattelaer <email address hidden> wrote:
>
> Did you check the born/loop/real matrix-element with the check function?
> If this is an issue related with the fermion flow it will likely spot the issue.
>
> Cheers,
>
> Olivier
>
> --
> You received this bug notification because you are subscribed to the bug
> report.
> https://bugs.launchpad.net/bugs/1899608
>
> Title:
> Fermion flow issue in the loops
>
> Status in MadGraph5_aMC@NLO:
> New
>
> Bug description:
> Hi everyone,
>
> We are performing some new physics calculations in a model where the
> SM is extended by a squark (and an invisible state that is irrelevant
> for our purpose). The UFO is provided in attachment. The process of
> interest consists in dijet production. I use MG5_aMC 2.6.8 and as a
> process definition, I typed
>
> import model RPVMono_UFO_RG-paper
> set complex_mass_scheme True
> generate p p > j j NP=2 QCD=0 QED=0 [QCD]
> output
>
> In order to solely get QCD corrections to the new physics component, I have modified the madgraph/loop/loop_diagram_generation.py file. I have set the filtering to true (around l388)
> edit_filter_manually = True
> and then included
> loop_pdgs = [abs(x) for x in diag.get_loop_lines_pdgs()]
> is_loop_gluon = (21 in loop_pdgs)
> is_loop_sigma = (9000002 in loop_pdgs)
> if is_loop_sigma and not is_loop_gluon:
> valid_diag=False
> I get the needed diagrams (so that this works).
>
> However, when launching a cross section calculation, I get a pole non-cancellation issue. The weird thing is that poles are fine for the subprocesses
> d dbar to s sbar or s sbar to d dbar
> but not for the subprocess
> dbar d to s sbar
> This is weird as the diagrams are exactly the same. I suspect something problematic with the fermion flow here (as we have clashing arrows diagrams). Any thoughts?
>
> Cheers,
>
> Benj
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1899608/+subscriptions

Revision history for this message
Benjamin Fuks (fuks) wrote :

Hi Olivier,

Any news here? Is there anything I could do on my side?

Cheers,

Benj

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) wrote : Re: [Bug 1899608] Fermion flow issue in the loops

In my case when trying to generate tthose loops they fail to generate with:

Command "check d d~ > s s~ NP=2 QCD=0 QED=0 [QCD]" interrupted with error:
ValueError : -5 is not in list
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

Cheers,

Olivier

> On 16 Oct 2020, at 08:09, Benjamin Fuks <email address hidden> wrote:
>
> Hi Olivier,
>
> Any news here? Is there anything I could do on my side?
>
> Cheers,
>
> Benj
>
> --
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1899608
>
> Title:
> Fermion flow issue in the loops
>
> Status in MadGraph5_aMC@NLO:
> New
>
> Bug description:
> Hi everyone,
>
> We are performing some new physics calculations in a model where the
> SM is extended by a squark (and an invisible state that is irrelevant
> for our purpose). The UFO is provided in attachment. The process of
> interest consists in dijet production. I use MG5_aMC 2.6.8 and as a
> process definition, I typed
>
> import model RPVMono_UFO_RG-paper
> set complex_mass_scheme True
> generate p p > j j NP=2 QCD=0 QED=0 [QCD]
> output
>
> In order to solely get QCD corrections to the new physics component, I have modified the madgraph/loop/loop_diagram_generation.py file. I have set the filtering to true (around l388)
> edit_filter_manually = True
> and then included
> loop_pdgs = [abs(x) for x in diag.get_loop_lines_pdgs()]
> is_loop_gluon = (21 in loop_pdgs)
> is_loop_sigma = (9000002 in loop_pdgs)
> if is_loop_sigma and not is_loop_gluon:
> valid_diag=False
> I get the needed diagrams (so that this works).
>
> However, when launching a cross section calculation, I get a pole non-cancellation issue. The weird thing is that poles are fine for the subprocesses
> d dbar to s sbar or s sbar to d dbar
> but not for the subprocess
> dbar d to s sbar
> This is weird as the diagrams are exactly the same. I suspect something problematic with the fermion flow here (as we have clashing arrows diagrams). Any thoughts?
>
> Cheers,
>
> Benj
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1899608/+subscriptions

Revision history for this message
Benjamin Fuks (fuks) wrote :
Download full text (4.3 KiB)

That's weird. Are you using the same version of mg5 as me (version 2.6.8)?

Cheers,

B>

> On 16 Oct 2020, at 09:33 , Olivier Mattelaer <email address hidden> wrote:
>
> In my case when trying to generate tthose loops they fail to generate
> with:
>
> Command "check d d~ > s s~ NP=2 QCD=0 QED=0 [QCD]" interrupted with error:
> ValueError : -5 is not in list
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> Cheers,
>
> Olivier
>
>> On 16 Oct 2020, at 08:09, Benjamin Fuks <email address hidden> wrote:
>>
>> Hi Olivier,
>>
>> Any news here? Is there anything I could do on my side?
>>
>> Cheers,
>>
>> Benj
>>
>> --
>> You received this bug notification because you are subscribed to
>> MadGraph5_aMC@NLO.
>> https://bugs.launchpad.net/bugs/1899608
>>
>> Title:
>> Fermion flow issue in the loops
>>
>> Status in MadGraph5_aMC@NLO:
>> New
>>
>> Bug description:
>> Hi everyone,
>>
>> We are performing some new physics calculations in a model where the
>> SM is extended by a squark (and an invisible state that is irrelevant
>> for our purpose). The UFO is provided in attachment. The process of
>> interest consists in dijet production. I use MG5_aMC 2.6.8 and as a
>> process definition, I typed
>>
>> import model RPVMono_UFO_RG-paper
>> set complex_mass_scheme True
>> generate p p > j j NP=2 QCD=0 QED=0 [QCD]
>> output
>>
>> In order to solely get QCD corrections to the new physics component, I have modified the madgraph/loop/loop_diagram_generation.py file. I have set the filtering to true (around l388)
>> edit_filter_manually = True
>> and then included
>> loop_pdgs = [abs(x) for x in diag.get_loop_lines_pdgs()]
>> is_loop_gluon = (21 in loop_pdgs)
>> is_loop_sigma = (9000002 in loop_pdgs)
>> if is_loop_sigma and not is_loop_gluon:
>> valid_diag=False
>> I get the needed diagrams (so that this works).
>>
>> However, when launching a cross section calculation, I get a pole non-cancellation issue. The weird thing is that poles are fine for the subprocesses
>> d dbar to s sbar or s sbar to d dbar
>> but not for the subprocess
>> dbar d to s sbar
>> This is weird as the diagrams are exactly the same. I suspect something problematic with the fermion flow here (as we have clashing arrows diagrams). Any thoughts?
>>
>> Cheers,
>>
>> Benj
>>
>> To manage notifications about this bug go to:
>> https://bugs.launchpad.net/mg5amcnlo/+bug/1899608/+subscriptions
>
> --
> You received this bug notification because you are subscribed to the bug
> report.
> https://bugs.launchpad.net/bugs/1899608
>
> Title:
> Fermion flow issue in the loops
>
> Status in MadGraph5_aMC@NLO:
> New
>
> Bug description:
> Hi everyone,
>
> We are performing some new physics calculations in a model where the
> SM is extended by a squark (and an invisible state that is irrelevant
> for our purpose). The UFO is provided in attachment. The process of
> interest consists in dijet production. I use MG5_aMC 2.6.8 and as a
> process definition, I typed
>
> import model RPVMono_UFO_RG-paper
> set complex_...

Read more...

Changed in mg5amcnlo:
assignee: nobody → marco zaro (marco-zaro)
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) wrote :

The issue does not seem to come from the loop. But Benj is going to re-generate the model and check if the problem persist

Changed in mg5amcnlo:
status: New → Incomplete
Revision history for this message
marco zaro (marco-zaro) wrote : Re: [Bug 1899608] Re: Fermion flow issue in the loops

Hi,
in case, please attach the check_poles.log file from the folder that fails.

Cheers,

Marco

> On 16 Oct 2020, at 11:58, Olivier Mattelaer <email address hidden> wrote:
>
> The issue does not seem to come from the loop. But Benj is going to re-
> generate the model and check if the problem persist
>
> ** Changed in: mg5amcnlo
> Status: New => Incomplete
>
> --
> You received this bug notification because you are a bug assignee.
> https://bugs.launchpad.net/bugs/1899608
>
> Title:
> Fermion flow issue in the loops
>
> Status in MadGraph5_aMC@NLO:
> Incomplete
>
> Bug description:
> Hi everyone,
>
> We are performing some new physics calculations in a model where the
> SM is extended by a squark (and an invisible state that is irrelevant
> for our purpose). The UFO is provided in attachment. The process of
> interest consists in dijet production. I use MG5_aMC 2.6.8 and as a
> process definition, I typed
>
> import model RPVMono_UFO_RG-paper
> set complex_mass_scheme True
> generate p p > j j NP=2 QCD=0 QED=0 [QCD]
> output
>
> In order to solely get QCD corrections to the new physics component, I have modified the madgraph/loop/loop_diagram_generation.py file. I have set the filtering to true (around l388)
> edit_filter_manually = True
> and then included
> loop_pdgs = [abs(x) for x in diag.get_loop_lines_pdgs()]
> is_loop_gluon = (21 in loop_pdgs)
> is_loop_sigma = (9000002 in loop_pdgs)
> if is_loop_sigma and not is_loop_gluon:
> valid_diag=False
> I get the needed diagrams (so that this works).
>
> However, when launching a cross section calculation, I get a pole non-cancellation issue. The weird thing is that poles are fine for the subprocesses
> d dbar to s sbar or s sbar to d dbar
> but not for the subprocess
> dbar d to s sbar
> This is weird as the diagrams are exactly the same. I suspect something problematic with the fermion flow here (as we have clashing arrows diagrams). Any thoughts?
>
> Cheers,
>
> Benj
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1899608/+subscriptions

Revision history for this message
Benjamin Fuks (fuks) wrote :

Hi Marco,

I have tried with a fresh UFO and nothing changes. The UFO seems to be therefore OK. The checks are the following:
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ddx_ssx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dsx_dsx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_sd_sd
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_sdx_sdx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ssx_ddx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxd_ssx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
Error detected in "calculate_xsect "
write debug file /Users/fuks/Work/tools/madgraph/2.6.8/PROCNLO_RPVMono_UFO_RG-paper_7/run_03_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Poles do not cancel, run cannot continue

I will attach the last log below, since all the others are fine.

Revision history for this message
Benjamin Fuks (fuks) wrote :

Hi Valentin,

After offline discussions on skype with Marco, it turns out that the issue is related to the treatment of the fermion flow in the loops. Therefore, I assign this issue to you :-)

If you want to treat it faster, feel free to contact me on skype.

Cheers,

Benj

PS: I have tested both with a fresh and with the old UFO (from our common 2015 paper), and both with 2.6.8 (python2) and 2.8.1 (python3). The problem is the same.

Revision history for this message
marco zaro (marco-zaro) wrote :

Hi Valentin,
Benjamin checked that when doing a standalone output of the virtuals, they return same born but different poles/finite part for
d d~ > s s~
and
d~ d > s s~.

I fear I can't help more than this....

Thanks!

Cheers,

Marco

Changed in mg5amcnlo:
assignee: marco zaro (marco-zaro) → Valentin Hirschi (valentin-hirschi)
Changed in mg5amcnlo:
status: Incomplete → Confirmed
Revision history for this message
Benjamin Fuks (fuks) wrote :

Hi everyone,

Any news?

Cheers,

Benj

Revision history for this message
Benjamin Fuks (fuks) wrote :

Hi everyone,

Has this been finally addressed?

Cheers,

Benj

To post a comment you must log in.
This report contains Public information  Edit
Everyone can see this information.

Other bug subscribers