Failure in generation DY+2j NLO process in MG5 v3.3.1

Dear Jin Choi,
I see you are using the loop_sm-ckm model. Do you need ckm mixing? was it working in v2?
Also, are you using massive tau's?

Can you try to disable collier as reduction algorithm?
Open Cards/MadLoopParams.dat

and edit
#MLReductionLib
 ! 6|7|1
to
6|1
(remove the exclamation mark at the beginning of the line)

Let me know

Best wishes,

Marco

Dear marco,

Acutually I have used customized loop_sm-ckm_no_b_mass, and I checked that P0_gg_tamtapuux GF8.0 process fails
to work in both loop_sm-ckm and loop_sm-no_b_mass. Could you explain more for v2?

In MGv3.3.1, the model worked fine with process
-----
generate p p > ell+ ell- [QCD]
add process p p > ell+ ell- j [QCD]
add process p p > ell+ ell- j j [QCD]
-----
So I suspect it is not the problem of input model.

Also, I tried to disable the collier and changed to OneLOop for reduction algorithm as you suggested, but it also runs forever with the following message:
-----
########################################################################
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
########################################################################
 VIRT: keeping split order 1

 Splitorders 2
       QCD: 6
       QED: 4
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: -2.3328491646858119E-007 OLP: -2.3328813399993523E-007
  COEFFICIENT SINGLE POLE:
        MadFKS: -5.7059918319874816E-008 OLP: -5.7044497254677827E-008

  SCALE**2: 3546.1736892112958
  Trying another PS point
 REAL 1: keeping split order 1
 REAL 2: keeping split order 1
 REAL 3: keeping split order 1
 REAL 4: keeping split order 1
 REAL 5: keeping split order 1
 REAL 6: keeping split order 1
 REAL 7: keeping split order 1
-----

in the same process.

Thanks for your help!

Best regards,
Jin Choi

Hello,
in the folder where you have that log (should be one among SubProcesses/P*/GF* ), can you please try to do a manual run?
cd in the folder

../madevent_mintMC < input_app.txt
and let me know how it goes?

Thanks!

Best wishes,

Marco

Hello,

The message above is from the manual run as you suggested,
I put the full log in https://choij.web.cern.ch/choij/GenValidation/RUN3/ThirteenTeV/DY2jNLO/log.txt

Thanks for your help!

Best,
Jin

Hi,
I cannot see any error in the log, it looks just truncated. Is there any error printed on screen?
Marco

> On 1 Mar 2022, at 15:01, Jin Choi <email address hidden> wrote:
>
> Question #700727 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/700727
>
> Status: Needs information => Open
>
> Jin Choi gave more information on the question:
> Hello,
>
> The message above is from the manual run as you suggested,
> I put the full log in https://choij.web.cern.ch/choij/GenValidation/RUN3/ThirteenTeV/DY2jNLO/log.txt
>
> Thanks for your help!
>
> Best,
> Jin
>
> --
> You received this question notification because you are assigned to this
> question.

Hello,

There were no error message in the log, and it just stopped at that point.
The message printed on the screen was #2. The process had been running at that point forever.

Best,
Jin

Hi,
I tried to generate e+e- j j

import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model loop_sm-ckm
generate p p > e+ e- j j [QCD]
output lljj-nloqcd

but everything goes through...

Best,

Marco

Hi Marco, Jin

I think the problem occurs when you include massive tau in your lepton container since it always fails with the process 'SubProcesses/P0_gg_tamtapuux/ajob1 8.0 F 0 0', other flavors such as e+e-jj or mu+mu-jj ends normally.

The weirdest thing that I observe here is that this happens only for the exclusive 2j while the inclusive jet process works, i.e.

define ell+ = e+ mu+ ta+
define ell- = e- mu- ta-
generate p p > ell+ ell- j j [QCD] # this doesn't work

while

define ell+ = e+ mu+ ta+
define ell- = e- mu- ta-
generate p p > ell+ ell- [QCD]
add process p p > ell+ ell- j [QCD]
add process p p > ell+ ell- j j [QCD] # this works

although the former should be a complete subset of the latter...

Best,
Sanghyun.

ok.
I will try the tau and let you know

Marco

> On 2 Mar 2022, at 11:01, Sanghyun Ko <email address hidden> wrote:
>
> Question #700727 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/700727
>
> Sanghyun Ko posted a new comment:
> Hi Marco, Jin
>
>
> I think the problem occurs when you include massive tau in your lepton container since it always fails with the process 'SubProcesses/P0_gg_tamtapuux/ajob1 8.0 F 0 0', other flavors such as e+e-jj or mu+mu-jj ends normally.
>
> The weirdest thing that I observe here is that this happens only for the
> exclusive 2j while the inclusive jet process works, i.e.
>
> define ell+ = e+ mu+ ta+
> define ell- = e- mu- ta-
> generate p p > ell+ ell- j j [QCD] # this doesn't work
>
> while
>
> define ell+ = e+ mu+ ta+
> define ell- = e- mu- ta-
> generate p p > ell+ ell- [QCD]
> add process p p > ell+ ell- j [QCD]
> add process p p > ell+ ell- j j [QCD] # this works
>
> although the former should be a complete subset of the latter...
>
>
> Best,
> Sanghyun.
>
> --
> You received this question notification because you are assigned to this
> question.