Problem with NLO correction

[Hridoy Debnath]

Hello,

 I have created a NLO UFO file using feynrules for a new model. It works perfectly okay for LO contributions. To check the NLO contributions , I have tried to generate p p > t t~ [QCD] process.

I got this error when I was trying generate events:

Error detected in "launch auto "

Also it shows as warning "WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD." and "WARNING: Could not compile StdHEP because its source directory could not be found in the SOURCE folder.
 Check the MG5_aMC option 'output_dependencies'.
 This will prevent the use of HERWIG6/Pythia6 shower"

I would really appreciate it if one could help me regarding what might be the possible source of errors/ solutions ?

Here I have attached the terminal output

MG5_aMC>generate p p > t t~ [QCD]
WARNING: Process: g/u/c/d/s/u~/c~/d~/s~ g/u/c/d/s/u~/c~/d~/s~ > t t~ WEIGHTED=2 [ all = QCD ] can have real emission processes which are not finite.
To avoid this, please use multiparticles when generating the process and be sure to include all the following particles in the multiparticle definition:
 g
INFO: Generating FKS-subtracted matrix elements for born process: g g > t t~ [ all = QCD ] (1 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: u u~ > t t~ [ all = QCD ] (2 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: c c~ > t t~ [ all = QCD ] (3 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: d d~ > t t~ [ all = QCD ] (4 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: s s~ > t t~ [ all = QCD ] (5 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: u~ u > t t~ [ all = QCD ] (6 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: c~ c > t t~ [ all = QCD ] (7 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: d~ d > t t~ [ all = QCD ] (8 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: s~ s > t t~ [ all = QCD ] (9 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: g g > t t~ NP<=2 QCD<=2 [ all = QCD ] (1 / 9)
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: u u~ > t t~ NP<=2 QCD<=2 [ all = QCD ] (2 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: c c~ > t t~ NP<=2 QCD<=2 [ all = QCD ] (3 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: d d~ > t t~ NP<=2 QCD<=2 [ all = QCD ] (4 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: s s~ > t t~ NP<=2 QCD<=2 [ all = QCD ] (5 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: u~ u > t t~ NP<=2 QCD<=2 [ all = QCD ] (6 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: c~ c > t t~ NP<=2 QCD<=2 [ all = QCD ] (7 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: d~ d > t t~ NP<=2 QCD<=2 [ all = QCD ] (8 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: s~ s > t t~ NP<=2 QCD<=2 [ all = QCD ] (9 / 9)
WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Generated 9 subprocesses with 232 real emission diagrams, 19 born diagrams and 185 virtual diagrams
MG5_aMC>output try4
INFO: Writing out the aMC@NLO code, using optimized Loops
INFO: initialize a new directory: try4
INFO: remove old information in try4
INFO: Generating real emission matrix-elements...
INFO: Generating Helas calls for FKS process: g g > t t~ [ all = QCD ] (1 / 9)
INFO: Processing color information for process: g g > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: u u~ > t t~ [ all = QCD ] (2 / 9)
INFO: Processing color information for process: u u~ > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: c c~ > t t~ [ all = QCD ] (3 / 9)
INFO: Reusing existing color information for process: c c~ > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: d d~ > t t~ [ all = QCD ] (4 / 9)
INFO: Reusing existing color information for process: d d~ > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: s s~ > t t~ [ all = QCD ] (5 / 9)
INFO: Reusing existing color information for process: s s~ > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: u~ u > t t~ [ all = QCD ] (6 / 9)
INFO: Processing color information for process: u~ u > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: c~ c > t t~ [ all = QCD ] (7 / 9)
INFO: Reusing existing color information for process: c~ c > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: d~ d > t t~ [ all = QCD ] (8 / 9)
INFO: Reusing existing color information for process: d~ d > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: s~ s > t t~ [ all = QCD ] (9 / 9)
INFO: Reusing existing color information for process: s~ s > t t~ [ all = QCD ]
INFO: ... Done
INFO: Writing files in P0_gg_ttx (1 / 9)
INFO: Creating files in directory V0_gg_ttx
ALOHA: aloha creates FFV1 set of routines with options: L1
aloha creates FFV1 set of routines with options: L1
aloha creates FFV1 set of routines with options: L1
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: g g > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_uux_ttx (2 / 9)
INFO: Creating files in directory V0_uux_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: u u~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_ccx_ttx (3 / 9)
INFO: Creating files in directory V0_ccx_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: c c~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: c c~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_ddx_ttx (4 / 9)
INFO: Creating files in directory V0_ddx_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d d~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: d d~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_ssx_ttx (5 / 9)
INFO: Creating files in directory V0_ssx_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: s s~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: s s~ > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_uxu_ttx (6 / 9)
INFO: Creating files in directory V0_uxu_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u~ u > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: u~ u > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_cxc_ttx (7 / 9)
INFO: Creating files in directory V0_cxc_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: c~ c > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: c~ c > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_dxd_ttx (8 / 9)
INFO: Creating files in directory V0_dxd_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d~ d > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: d~ d > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Writing files in P0_sxs_ttx (9 / 9)
INFO: Creating files in directory V0_sxs_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: s~ s > t t~ NP<=2 QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: s~ s > t t~ NP<=2 QCD<=2 [ all = QCD ]
History written to /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/Cards/proc_card_mg5.dat
ALOHA: aloha creates UUV1 set of routines with options: L2
aloha creates UUV1 set of routines with options: L2
aloha creates UUV1 set of routines with options: L2
ALOHA: aloha creates VVV2 routines
aloha creates VVV2 routines
aloha creates VVV2 routines
ALOHA: aloha creates FFV1 routines
aloha creates FFV1 routines
aloha creates FFV1 routines
ALOHA: aloha creates FFV6 routines
aloha creates FFV6 routines
aloha creates FFV6 routines
ALOHA: aloha creates VVV1 routines
aloha creates VVV1 routines
aloha creates VVV1 routines
ALOHA: aloha creates VV3 routines
aloha creates VV3 routines
aloha creates VV3 routines
ALOHA: aloha creates VV2 routines
aloha creates VV2 routines
aloha creates VV2 routines
ALOHA: aloha creates VV5 routines
aloha creates VV5 routines
aloha creates VV5 routines
ALOHA: aloha creates VV4 routines
aloha creates VV4 routines
aloha creates VV4 routines
ALOHA: aloha creates FF1 routines
aloha creates FF1 routines
aloha creates FF1 routines
ALOHA: aloha creates FF2 routines
aloha creates FF2 routines
aloha creates FF2 routines
ALOHA: aloha creates VV1 routines
aloha creates VV1 routines
aloha creates VV1 routines
ALOHA: aloha creates FFV5 routines
aloha creates FFV5 routines
aloha creates FFV5 routines
ALOHA: aloha creates VVV3 routines
aloha creates VVV3 routines
aloha creates VVV3 routines
ALOHA: aloha creates FFV1_6 routines
aloha creates FFV1_6 routines
aloha creates FFV1_6 routines
ALOHA: aloha creates FF1_2 routines
aloha creates FF1_2 routines
aloha creates FF1_2 routines
ALOHA: aloha creates FFV1_5 routines
aloha creates FFV1_5 routines
aloha creates FFV1_5 routines
The option max_npoint_for_channel is modified [0] but will not be written in the configuration files.
If you want to make this value the default for future session, you can run 'save options --all'
The option default_unset_couplings is modified [99] but will not be written in the configuration files.
If you want to make this value the default for future session, you can run 'save options --all'
The option max_t_for_channel is modified [99] but will not be written in the configuration files.
If you want to make this value the default for future session, you can run 'save options --all'
The option auto_update is modified [7] but will not be written in the configuration files.
If you want to make this value the default for future session, you can run 'save options --all'
save configuration file to /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/Cards/amcatnlo_configuration.txt
INFO: The value for lhapdf in the current configuration does not correspond to a valid executable.
Please set it correctly either in input/mg5_configuration or with "set lhapdf /path/to/lhapdf-config" and regenrate the process.
To avoid regeneration, edit the /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/Cards/amcatnlo_configuration.txt file.
Note that you can still compile and run aMC@NLO with the built-in PDFs

INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
INFO: Generating MadAnalysis5 default cards tailored to this process
Type "launch" to generate events from this process, or see
/home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch try4
INFO: ************************************************************
* *
* W E L C O M E to M A D G R A P H 5 *
* a M C @ N L O *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.9.17 2023-10-19 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/Cards/amcatnlo_configuration.txt
INFO: load configuration from /home/hridoy/Desktop/MG5_aMC_v2_9_16/input/mg5_configuration.txt
INFO: load configuration from /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/Cards/amcatnlo_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
launch auto
The following switches determine which programs are run:
/================== Description ===================|================ values ================|====== other options =======\
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = HERWIG6 | OFF|PYTHIA8 |
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\========================================================================================================================/
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>0
INFO: will run in mode: noshower
WARNING: You have chosen not to run a parton shower.
    NLO events without showering are NOT physical.
    Please, shower the LesHouches events before using them for physics analyses.
    You have to choose NOW which parton-shower you WILL use and specify it in the run_card.
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.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, run, enter path][90s to answer]
>0
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/...)
INFO: Starting run
INFO: Compiling the code
INFO: For gauge cancellation, the width of 't' has been set to zero.
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
WARNING: Could not compile StdHEP because its source directory could not be found in the SOURCE folder.
 Check the MG5_aMC option 'output_dependencies'.
 This will prevent the use of HERWIG6/Pythia6 shower.
INFO: Compiling directories...
INFO: Compiling on 8 cores
INFO: Compiling P0_gg_ttx...
INFO: Compiling P0_uux_ttx...
INFO: Compiling P0_ccx_ttx...
INFO: Compiling P0_ddx_ttx...
INFO: Compiling P0_ssx_ttx...
INFO: Compiling P0_uxu_ttx...
INFO: Compiling P0_cxc_ttx...
INFO: Compiling P0_dxd_ttx...
INFO: P0_uux_ttx done.
INFO: Compiling P0_sxs_ttx...
INFO: P0_uxu_ttx done.
INFO: P0_dxd_ttx done.
INFO: P0_ssx_ttx done.
INFO: P0_cxc_ttx done.
INFO: P0_ddx_ttx done.
INFO: P0_ccx_ttx done.
INFO: P0_gg_ttx done.
INFO: P0_sxs_ttx done.
INFO: Checking test output:
INFO: P0_gg_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uux_ttx
INFO: Result for test_ME:
INFO: Output of the failing test:
 Enter 0 to compute MC/MC(limit)
       1 to compute MC/ME(limit)
       2 to compute ME/ME(limit)
 Enter xi_i, y_ij to be used in coll/soft tests
  Enter -2 to generate them randomly
 Enter number of tests for soft and collinear limits
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11839999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
 Give FKS configuration number ("0" loops over all)

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 5 j= 1
 with PDGs: i= 21 j= 2

 Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 1 , 2 , 0
  with seed 32
 Ranmar initialization seeds 10730 9408
tau_min 1 1 : 0.34400E+03 -- 0.34400E+03
tau_min 2 1 : 0.34400E+03 -- 0.34400E+03
tau_min 3 1 : 0.34400E+03 0.34400E+03 0.34400E+03
tau_min 4 1 : 0.34400E+03 0.34400E+03 0.34400E+03
tau_min 5 1 : 0.34400E+03 -- 0.34400E+03
tau_min 6 1 : 0.34400E+03 -- 0.34400E+03

 Soft tests done for (Born) config 1
 Failures: 34
     Soft test 1 FAILED. Fraction of failures: 0.34

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 2
 FKS partons are: i= 5 j= 2
 with PDGs: i= 21 j= -2

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 26
     Soft test 2 PASSED. Fraction of failures: 0.26

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 3
 FKS partons are: i= 5 j= 3
 with PDGs: i= 21 j= 6

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 3 PASSED. Fraction of failures: 0.00

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 4
 FKS partons are: i= 5 j= 4
 with PDGs: i= 21 j= -6

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 0
     Soft test 4 PASSED. Fraction of failures: 0.00

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 5
 FKS partons are: i= 5 j= 1
 with PDGs: i= -2 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 2
     Soft test 5 PASSED. Fraction of failures: 0.02

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 6
 FKS partons are: i= 5 j= 2
 with PDGs: i= 2 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 1
     Soft test 6 PASSED. Fraction of failures: 0.01

 No collinear test for massive j_fks
Error detected in "launch auto "
write debug file /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str : Some tests failed, run cannot continue.
 Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.

[Sihyun Jeon (shjeon)]

PYTHIA8=parton_shower

change this in the run_card and see if that works.
maybe MC@NLO for herwig6/pp is not supported anymore because it's too old.

then what you have to assure yourself is that you want to shower LHE file with pythia8, not herwig

[Hridoy Debnath (hridoy-1)]

Hi Sihyun

Thanks a lot. I have tried to run it with pythia8, but still having the same issue.

No collinear test for massive j_fks
Error detected in "launch auto "
write debug file /home/hridoy/Desktop/MG5_aMC_v2_9_16/try4/run_02_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str : Some tests failed, run cannot continue.
 Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.

[Olivier Mattelaer (olivier-mattelaer)]

This will avoid the issue that you do not have stdhep installed (or even installable on your system).

But you have here at least two others issue:

1)WARNING: All Born diagrams do not factorize the same sum of power(s) of the the perturbed order(s) QCD.
This is potentially dangerous as the real-emission diagrams from aMC@NLO will not be consistent with these virtual contributions.

This warning is quite serious, and is typically followed later in the code by a crash (in version 2.9.x) when computing the actual power of alpha_s for the born. (version 3 is more flexible here).
If this warning is what you expect then 2.9.x is not for you and you need to use 3.x for such kind of NLO computation (but you need to use the framework of double expansion in that case and be clear what type of correction you want to include --3.1 will likely complain that your syntax is ambiguous on that prospect).

Note that you actually do not see such crash yet in 2.9, since you face another crash first (see next point).

2) str : Some tests failed, run cannot continue.

So here, you have an issue with the behaviour of your amplitude in the soft limit:
technical detail in the configuration where the soft limit is not correct:
FKS configuration number is 1
 FKS partons are: i= 5 j= 1
 with PDGs: i= 21 j= 2
for config =1

In all NLO compuation, we first check for a couple of phase-space pont the soft and collinear limit of the amplitude in order to check that substraction of the divergence will be done properly for the process.
For your BSM model, looks like one soft limit of your amplitude is not the compatible with the FKS substraction scheme and therefore the result will not be NLO accurate.

Here I do not know if this is related to a bad handling of the double expansion (like pointed above) or a more serious issue.

Cheers,

Olivier

[Hridoy Debnath (hridoy-1)]

Hi Olivier

Thanks a lot for your response. I really appreciate it.

Do you have any suggestions or comments about how to solve these issues, please?

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

The first suggestion is to try 3.5.X release.
This should fix the first issue. For the second issue, it might also be fixed with that version, if not I do not have any real comment/suggestion for you, expect studying the soft limit of your amplitude and understand what the issue is with such limit.

Cheers,

Olivier