MadGraph5_aMC@NLO

Integration channel issue while generating e- e+ > 4j [QCD]

Asked by Swapan Kumar Majhi

Dear experts,
We are trying to generate four jet events at QCD NLO in electron-postiron collider. We are using loop_qcd_qed_sm_Gmu model. We have found there is no pole cancellation issues but lhe file is not getting as an output. Runs only avialable for fixed-order. Integraation channle issues are there. We paste the output below,

************************************************************
* *
* 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 3.5.2 2023-11-08 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://madgraph.phys.ucl.ac.be/ *
* and *
* http://amcatnlo.web.cern.ch/amcatnlo/ *
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
************************************************************
load MG5 configuration from ../input/mg5_configuration.txt
fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
 Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
 MG5_aMC> set fastjet /PATH/TO/fastjet-config

eMELA-config does not seem to correspond to a valid eMELA-config executable.
 Please set the 'fastjet'variable to the full (absolute) /PATH/TO/eMELA-config (including eMELA-config).
 MG5_aMC> set eMELA /PATH/TO/eMELA-config

set ninja to /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/HEPTools/lib
set collier to /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/HEPTools/lib
lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

set lhapdf to /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config
Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
import /home/q/Documents/Computations/HEP/Projects/ee4j/errors/4j
The import format was not given, so we guess it as command
import model loop_qcd_qed_sm_Gmu
Switching to Feynman gauge because it is the only one supported by the model loop_qcd_qed_sm_Gmu.
INFO: Restrict model loop_qcd_qed_sm_Gmu with file ../models/loop_qcd_qed_sm_Gmu/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
Pass the definition of 'j' and 'p' to 5 flavour scheme, including the photon.
Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt e- mu- ta- ve~ vm~ vt~ e+ mu+ ta+ t t~ z w+ ghz ghwp ghwm h g0 g+ w- ghz~ ghwp~ ghwm~ g-
set complex_mass_scheme True
Activate complex mass scheme.
INFO: Restrict model loop_qcd_qed_sm_Gmu with file ../models/loop_qcd_qed_sm_Gmu/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g ghg ghg~ u c d s b u~ c~ d~ s~ b~ a gha gha~ ve vm vt e- mu- ta- ve~ vm~ vt~ e+ mu+ ta+ t t~ z w+ ghz ghwp ghwm h g0 g+ w- ghz~ ghwp~ ghwm~ g-
set gauge Feynman
define j = g u u~ c c~ d d~ s s~ b b~
Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
generate e- e+ > 4j aEW=4 aS=2 [QCD]
change syntax aS=2 to QCD^2=4 to correspond to UFO model convention
Interpreting 'QCD^2=4' as 'QCD^2<=4'
change syntax aEW=4 to QED^2=8 to correspond to UFO model convention
Interpreting 'QED^2=8' as 'QED^2<=8'
The process definition has leptons in the initial state
The flag 'include_lepton_initiated_processes' will be set to True
WARNING: Use of multiparticles is non-trivial for NLO process generation and depends on the orders included, the process considered, as well as the PDF set chosen. See appendix D of arXiv:1804.10017 [hep-ph] for some guidance. 
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > g g u u~ QCD^2<=6 QED^2<=8 (1 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > g g c c~ QCD^2<=6 QED^2<=8 (2 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > g g d d~ QCD^2<=6 QED^2<=8 (3 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > g g s s~ QCD^2<=6 QED^2<=8 (4 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > g g b b~ QCD^2<=6 QED^2<=8 (5 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > u u u~ u~ QCD^2<=6 QED^2<=8 (6 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > u c u~ c~ QCD^2<=6 QED^2<=8 (7 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > u d u~ d~ QCD^2<=6 QED^2<=8 (8 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > u s u~ s~ QCD^2<=6 QED^2<=8 (9 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > u s c~ d~ QCD^2<=6 QED^2<=8 (10 / 22)
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > u b u~ b~ QCD^2<=6 QED^2<=8 (11 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > c c c~ c~ QCD^2<=6 QED^2<=8 (12 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > c d u~ s~ QCD^2<=6 QED^2<=8 (13 / 22)
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > c d c~ d~ QCD^2<=6 QED^2<=8 (14 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > c s c~ s~ QCD^2<=6 QED^2<=8 (15 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > c b c~ b~ QCD^2<=6 QED^2<=8 (16 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > d d d~ d~ QCD^2<=6 QED^2<=8 (17 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > d s d~ s~ QCD^2<=6 QED^2<=8 (18 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > d b d~ b~ QCD^2<=6 QED^2<=8 (19 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > s s s~ s~ QCD^2<=6 QED^2<=8 (20 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > s b s~ b~ QCD^2<=6 QED^2<=8 (21 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating FKS-subtracted matrix elements for born process: e- e+ > b b b~ b~ QCD^2<=6 QED^2<=8 (22 / 22)
INFO: link_born_real: skipping because not all splittings are QCD
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > u u~ g g QCD^2<=6 QED^2<=8 (1 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > c c~ g g QCD^2<=6 QED^2<=8 (2 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > d d~ g g QCD^2<=6 QED^2<=8 (3 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > s s~ g g QCD^2<=6 QED^2<=8 (4 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > b b~ g g QCD^2<=6 QED^2<=8 (5 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > u u u~ u~ QCD^2<=6 QED^2<=8 (6 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > c u u~ c~ QCD^2<=6 QED^2<=8 (7 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > u d d~ u~ QCD^2<=6 QED^2<=8 (8 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > s u u~ s~ QCD^2<=6 QED^2<=8 (9 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > s u d~ c~ QCD^2<=6 QED^2<=8 (10 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > b u u~ b~ QCD^2<=6 QED^2<=8 (11 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > c c c~ c~ QCD^2<=6 QED^2<=8 (12 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > c d u~ s~ QCD^2<=6 QED^2<=8 (13 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > c d d~ c~ QCD^2<=6 QED^2<=8 (14 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > c s s~ c~ QCD^2<=6 QED^2<=8 (15 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > b c c~ b~ QCD^2<=6 QED^2<=8 (16 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > d d d~ d~ QCD^2<=6 QED^2<=8 (17 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > s d d~ s~ QCD^2<=6 QED^2<=8 (18 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > b d d~ b~ QCD^2<=6 QED^2<=8 (19 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > s s s~ s~ QCD^2<=6 QED^2<=8 (20 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > b s s~ b~ QCD^2<=6 QED^2<=8 (21 / 22)
INFO: Generating virtual matrix element with MadLoop for process: e- e+ > b b b~ b~ QCD^2<=6 QED^2<=8 (22 / 22)
INFO: Generated 22 subprocesses with 24300 real emission diagrams, 764 born diagrams and 79627 virtual diagrams
------------------------------------------------------------------------
This computation can involve not only purely SM-QCD corrections at NLO.
Please also cite ref. 'arXiv:1804.10017' when using results from this code.
------------------------------------------------------------------------

display processes
Born processes:
Process: e- e+ > u u~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u u~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c u u~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u d d~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b u u~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c c~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c s s~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b c c~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d d~ d~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s d d~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b d d~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s s~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b s s~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b b~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Real processes:
Process: e- e+ > u u~ g g a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > u u~ g g e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > u u~ g g e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u~ g g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u~ d d~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u~ s s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u~ b b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u~ u u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u~ c c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ g g a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > c c~ g g e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > c c~ g g e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ g g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ d d~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ s s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ b b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ u u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ c c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ g g a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > d d~ g g e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > d d~ g g e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ g g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ d d~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ s s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ b b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ u u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ c c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ g g a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > s s~ g g e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > s s~ g g e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ g g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ d d~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ s s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ b b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ u u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ c c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ g g a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > b b~ g g e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > b b~ g g e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ g g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ d d~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ s s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ b b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ u u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ c c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u u~ u~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > u u u~ u~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > u u u~ u~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u u~ u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c u u~ c~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > c u u~ c~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > c u u~ c~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c u u~ c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u d d~ u~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > u d d~ u~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > u d d~ u~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u d d~ u~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u u~ s~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > s u u~ s~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > s u u~ s~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u u~ s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u d~ c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b u u~ b~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > b u u~ b~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > b u u~ b~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b u u~ b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c c~ c~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > c c c~ c~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > c c c~ c~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c c~ c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d u~ s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d d~ c~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > c d d~ c~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > c d d~ c~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d d~ c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c s s~ c~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > c s s~ c~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > c s s~ c~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c s s~ c~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b c c~ b~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > b c c~ b~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > b c c~ b~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b c c~ b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d d~ d~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > d d d~ d~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > d d d~ d~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d d~ d~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s d d~ s~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > s d d~ s~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > s d d~ s~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s d d~ s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b d d~ b~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > b d d~ b~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > b d d~ b~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b d d~ b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s s~ s~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > s s s~ s~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > s s s~ s~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s s~ s~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b s s~ b~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > b s s~ b~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > b s s~ b~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b s s~ b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b b~ b~ a [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: a e+ > b b b~ b~ e+ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- a > b b b~ b~ e- [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b b~ b~ g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Loop processes:
Process: e- e+ > u u~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u u u~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c u u~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > u d d~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s u d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b u u~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c c c~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c d d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > c s s~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b c c~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > d d d~ d~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s d d~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b d d~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > s s s~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b s s~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
Process: e- e+ > b b b~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
output 4j_NLO
INFO: Writing out the aMC@NLO code, using optimized Loops
INFO: directory /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO already exists.
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n]
found line : launch
This answer is not valid for current question. Keep it for next question and use here default: y
INFO: initialize a new directory: 4j_NLO
INFO: remove old information in 4j_NLO
INFO: Generating real emission matrix-elements...
INFO: Generating Helas calls for FKS process: e- e+ > u u~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8 (1 / 22)
INFO: Processing color information for process: e- e+ > u u~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > c c~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8 (2 / 22)
INFO: Reusing existing color information for process: e- e+ > c c~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > d d~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8 (3 / 22)
INFO: Reusing existing color information for process: e- e+ > d d~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > s s~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8 (4 / 22)
INFO: Reusing existing color information for process: e- e+ > s s~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > b b~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8 (5 / 22)
INFO: Reusing existing color information for process: e- e+ > b b~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > u u u~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (6 / 22)
INFO: Processing color information for process: e- e+ > u u u~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > c u u~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (7 / 22)
INFO: Processing color information for process: e- e+ > c u u~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > u d d~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (8 / 22)
INFO: Processing color information for process: e- e+ > u d d~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > s u u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (9 / 22)
INFO: Reusing existing color information for process: e- e+ > s u u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > s u d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (10 / 22)
INFO: Processing color information for process: e- e+ > s u d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > b u u~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (11 / 22)
INFO: Reusing existing color information for process: e- e+ > b u u~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > c c c~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (12 / 22)
INFO: Generating Helas calls for FKS process: e- e+ > c d u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (13 / 22)
INFO: Reusing existing color information for process: e- e+ > c d u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > c d d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (14 / 22)
INFO: Reusing existing color information for process: e- e+ > c d d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > c s s~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (15 / 22)
INFO: Generating Helas calls for FKS process: e- e+ > b c c~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (16 / 22)
INFO: Generating Helas calls for FKS process: e- e+ > d d d~ d~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (17 / 22)
INFO: Reusing existing color information for process: e- e+ > d d d~ d~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > s d d~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (18 / 22)
INFO: Reusing existing color information for process: e- e+ > s d d~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > b d d~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (19 / 22)
INFO: Reusing existing color information for process: e- e+ > b d d~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating Helas calls for FKS process: e- e+ > s s s~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (20 / 22)
INFO: Generating Helas calls for FKS process: e- e+ > b s s~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (21 / 22)
INFO: Generating Helas calls for FKS process: e- e+ > b b b~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8 (22 / 22)
INFO: Reusing existing color information for process: e- e+ > b b b~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: ... Done
INFO: Writing files in P0_emep_uuxgg (1 / 17)
INFO: Creating files in directory V0_emep_uuxgg
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > u u~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > u u~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_ccxgg (2 / 17)
INFO: Creating files in directory V0_emep_ccxgg
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > c c~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > c c~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_ddxgg (3 / 17)
INFO: Creating files in directory V0_emep_ddxgg
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > d d~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > d d~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_ssxgg (4 / 17)
INFO: Creating files in directory V0_emep_ssxgg
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > s s~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > s s~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_bbxgg (5 / 17)
INFO: Creating files in directory V0_emep_bbxgg
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > b b~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > b b~ g g [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_uuuxux (6 / 17)
INFO: Creating files in directory V0_emep_uuuxux
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > u u u~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > u u u~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_cuuxcx (7 / 17)
INFO: Creating files in directory V0_emep_cuuxcx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > c u u~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > c u u~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_uddxux (8 / 17)
INFO: Creating files in directory V0_emep_uddxux
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > u d d~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > u d d~ u~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_suuxsx (9 / 17)
INFO: Creating files in directory V0_emep_suuxsx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > s u u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > s u u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_sudxcx (10 / 17)
INFO: Creating files in directory V0_emep_sudxcx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > s u d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > s u d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_buuxbx (11 / 17)
INFO: Creating files in directory V0_emep_buuxbx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > b u u~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > b u u~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_cduxsx (12 / 17)
INFO: Creating files in directory V0_emep_cduxsx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > c d u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > c d u~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_cddxcx (13 / 17)
INFO: Creating files in directory V0_emep_cddxcx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > c d d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > c d d~ c~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_dddxdx (14 / 17)
INFO: Creating files in directory V0_emep_dddxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > d d d~ d~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > d d d~ d~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_sddxsx (15 / 17)
INFO: Creating files in directory V0_emep_sddxsx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > s d d~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > s d d~ s~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_bddxbx (16 / 17)
INFO: Creating files in directory V0_emep_bddxbx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > b d d~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > b d d~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Writing files in P0_emep_bbbxbx (17 / 17)
INFO: Creating files in directory V0_emep_bbbxbx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: e- e+ > b b b~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
INFO: Generating born Feynman diagrams for Process: e- e+ > b b b~ b~ [ all = QCD QED ] QCD^2<=6 QED^2<=8
History written to /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO/Cards/proc_card_mg5.dat
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 138 routines in 2.162 s
The option complex_mass_scheme is modified [False] 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 include_lepton_initiated_processes is modified [False] 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 gauge is modified [unitary] 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/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO/Cards/amcatnlo_configuration.txt
INFO: Use Fortran compiler /usr/bin/gfortran
INFO: Use c++ compiler /usr/bin/g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Type "launch" to generate events from this process, or see
/home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO/README
Run "open index.html" to see more information about this process.
launch
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 3.5.2 2023-11-08 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO/Cards/amcatnlo_configuration.txt
INFO: load configuration from /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/input/mg5_configuration.txt
INFO: load configuration from /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO/Cards/amcatnlo_configuration.txt
launch auto
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
WARNING: NLO+PS mode is not allowed for processes including electroweak corrections 
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 | No NLO+PS available for EW correction |
\================================================================================================================================/
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.
INFO: will run in mode: NLO
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
set nevents 1000
INFO: modify parameter nevents of the run_card.dat to 1000 
INFO: Update the dependent parameter of the param_card.dat
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.119 
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
INFO: Starting run
INFO: Compiling the code
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 16 cores
INFO: Compiling P0_emep_uuxgg...
INFO: Compiling P0_emep_ccxgg...
INFO: Compiling P0_emep_ddxgg...
INFO: Compiling P0_emep_ssxgg...
INFO: Compiling P0_emep_bbxgg...
INFO: Compiling P0_emep_uuuxux...
INFO: Compiling P0_emep_cuuxcx...
INFO: Compiling P0_emep_uddxux...
INFO: Compiling P0_emep_suuxsx...
INFO: Compiling P0_emep_sudxcx...
INFO: Compiling P0_emep_buuxbx...
INFO: Compiling P0_emep_cduxsx...
INFO: Compiling P0_emep_cddxcx...
INFO: Compiling P0_emep_dddxdx...
INFO: Compiling P0_emep_sddxsx...
INFO: Compiling P0_emep_bddxbx...
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
INFO: P0_emep_cduxsx done.
INFO: Compiling P0_emep_bbbxbx...
INFO: P0_emep_sudxcx done.
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
INFO: P0_emep_cddxcx done.
STOP 1
INFO: P0_emep_suuxsx done.
INFO: P0_emep_sddxsx done.
INFO: P0_emep_bddxbx done.
INFO: P0_emep_cuuxcx done.
INFO: P0_emep_buuxbx done.
STOP 1
INFO: P0_emep_uddxux done.
STOP 1
STOP 1
INFO: P0_emep_dddxdx done.
INFO: P0_emep_uuuxux done.
STOP 1
INFO: P0_emep_bbbxbx done.
STOP 1
INFO: P0_emep_ccxgg done.
STOP 1
INFO: P0_emep_uuxgg done.
STOP 1
STOP 1
STOP 1
INFO: P0_emep_ssxgg done.
INFO: P0_emep_ddxgg done.
INFO: P0_emep_bbxgg done.
INFO: Checking test output:
INFO: P0_emep_uuxgg
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_emep_ccxgg
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_emep_ddxgg
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_emep_ssxgg
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_emep_bbxgg
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_emep_uuuxux
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_emep_cuuxcx
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_emep_uddxux
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_emep_suuxsx
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_emep_sudxcx
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_emep_buuxbx
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_emep_cduxsx
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_emep_cddxcx
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_emep_dddxdx
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_emep_sddxsx
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_emep_bddxbx
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_emep_bbbxbx
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: Starting run
INFO: Using 16 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
WARNING: No integration channels found for contribution P0_emep_uuxgg 
WARNING: No integration channels found for contribution P0_emep_ccxgg 
WARNING: No integration channels found for contribution P0_emep_ddxgg 
WARNING: No integration channels found for contribution P0_emep_ssxgg 
WARNING: No integration channels found for contribution P0_emep_bbxgg 
WARNING: No integration channels found for contribution P0_emep_uuuxux 
WARNING: No integration channels found for contribution P0_emep_cuuxcx 
WARNING: No integration channels found for contribution P0_emep_uddxux 
WARNING: No integration channels found for contribution P0_emep_suuxsx 
WARNING: No integration channels found for contribution P0_emep_buuxbx 
WARNING: No integration channels found for contribution P0_emep_cddxcx 
WARNING: No integration channels found for contribution P0_emep_dddxdx 
WARNING: No integration channels found for contribution P0_emep_sddxsx 
WARNING: No integration channels found for contribution P0_emep_bddxbx 
WARNING: No integration channels found for contribution P0_emep_bbbxbx 
INFO: Setting up grids
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 15h17 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 22.8s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 23.6s ]
sum of cpu time of last step: 0 second
INFO:
      Results after grid setup:
      Total cross section: 4.250e-01 +- 1.7e-02 pb

INFO: Refining results, step 1
INFO: Idle: 8, Running: 16, Completed: 0 [ current time: 15h18 ]
INFO: Idle: 7, Running: 16, Completed: 1 [ 1m 48s ]
INFO: Idle: 6, Running: 16, Completed: 2 [ 1m 49s ]
INFO: Idle: 5, Running: 16, Completed: 3 [ 1m 51s ]
INFO: Idle: 4, Running: 16, Completed: 4 [ 1m 51s ]
INFO: Idle: 3, Running: 16, Completed: 5 [ 1m 51s ]
INFO: Idle: 2, Running: 16, Completed: 6 [ 1m 51s ]
INFO: Idle: 1, Running: 16, Completed: 7 [ 1m 54s ]
INFO: Idle: 0, Running: 16, Completed: 8 [ 1m 57s ]
INFO: Idle: 0, Running: 15, Completed: 9 [ 2m 5s ]
INFO: Idle: 0, Running: 14, Completed: 10 [ 2m 5s ]
INFO: Idle: 0, Running: 13, Completed: 11 [ 2m 5s ]
INFO: Idle: 0, Running: 12, Completed: 12 [ 2m 5s ]
INFO: Idle: 0, Running: 11, Completed: 13 [ 2m 5s ]
INFO: Idle: 0, Running: 10, Completed: 14 [ 2m 5s ]
INFO: Idle: 0, Running: 9, Completed: 15 [ 2m 6s ]
INFO: Idle: 0, Running: 8, Completed: 16 [ 2m 7s ]
INFO: Idle: 0, Running: 7, Completed: 17 [ 2m 58s ]
INFO: Idle: 0, Running: 6, Completed: 18 [ 2m 59s ]
INFO: Idle: 0, Running: 5, Completed: 19 [ 2m 59s ]
INFO: Idle: 0, Running: 4, Completed: 20 [ 3m 0s ]
INFO: Idle: 0, Running: 3, Completed: 21 [ 3m 0s ]
INFO: Idle: 0, Running: 2, Completed: 22 [ 3m 0s ]
INFO: Idle: 0, Running: 1, Completed: 23 [ 3m 0s ]
INFO: Idle: 0, Running: 0, Completed: 24 [ 3m 2s ]
sum of cpu time of last step: 0 second
INFO:
      Current results:
      Total cross section: 4.374e-01 +- 6.1e-03 pb

INFO: Refining results, step 2
INFO: Idle: 0, Running: 16, Completed: 0 [ current time: 15h21 ]
INFO: Idle: 0, Running: 0, Completed: 16 [ 5m 36s ]
sum of cpu time of last step: 0 second
INFO:
   --------------------------------------------------------------
      Final results and run summary:
      Process e- e+ > 4j aEW=4 aS=2 [QCD]
      Run at l-l collider (500.0 + 500.0 GeV)
      Total cross section: 4.335e-01 +- 2.4e-03 pb
   --------------------------------------------------------------
      Scale variation (computed from histogram information):
          Dynamical_scale_choice -1 (envelope of 9 values):
              4.335e-01 pb +0.4% -0.3%
   --------------------------------------------------------------

INFO: The results of this run and the HwU and GnuPlot files with the plots have been saved in /home/q/Documents/Computations/HEP/MG5_aMC_v3_5_1/bin/4j_NLO/Events/run_01
INFO: Run complete
INFO:
quit
INFO:
quit

Why we are getting "STOP 1" message and "No integration channel issues"? Kindly suggest on way of getting the lhe output!

Any suggestions and comments are highly appreciated.

Thank you

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi,
you have this warning
WARNING: NLO+PS mode is not allowed for processes including electroweak corrections 
and the reason is likely that you generate the process with
generate e- e+ > 4j aEW=4 aS=2 [QCD]
 This includes contributions at the born of the order (squared-ME level) aEW=2 aS=2 ; aEW=3 aS=1 ; aEW=4 aS=0.

The proper syntax, to be able to get events (involving only QCD corrections) is to generate only the term with most powers of alphaS
generate e- e+ > 4j aEW=2 aS=2 [QCD]

LEt me know if this solves your issue

Cheers,

Marco

Revision history for this message
user1905 (userinhep) said : 		#2

Dear Marco,

Thanks for the quick response.

We did checked with the syntax you suggested but we do have a different goal. the idea was to consider the NLO QCD to WW production and for that the minimum syntax needed was,

generate e- e+ > 4j aEW<=4 aS=2 [QCD],

and this syntax is having all the error. Can you suggest if there is any safe way to handle the above process.

regards

amir

Revision history for this message
user1905 (userinhep) said : 		#3

Just to add an additional info to this,

generate e- e+ > 4j aEW<=4 [QCD] works fine if we define j to be quarks of different flavor, but if j are taken to be of same flavor there is an integration channel issues.

Also we have noted that event though we have explicitly defined j devoid of photon, we obtain photon at the final states.

Regards
amir

Revision history for this message
marco zaro (marco-zaro) said : 		#4

Hi,
in the case of 4j you have big troubles in the purely-EW contribution. In particular, in principle, since you have photons which can be jets (they are needed for IR safety), they can split (in principle) also into leptons.

For contributions to a given process beyond the LO with the maximum power of alpha_s, event generation is not possible for QCD corrections.

Best wishes,

Marco

> On 23 May 2024, at 12:20, user1905 <email address hidden> wrote:
>
> Question #816013 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/816013
>
> user1905 proposed the following answer:
> Dear Marco,
>
> Thanks for the quick response.
>
> We did checked with the syntax you suggested but we do have a different
> goal. the idea was to consider the NLO QCD to WW production and for that
> the minimum syntax needed was,
>
> generate e- e+ > 4j aEW<=4 aS=2 [QCD],
>
> and this syntax is having all the error. Can you suggest if there is any
> safe way to handle the above process.
>
> regards
>
> amir
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
Swapan Kumar Majhi (majhi-majhi) said : 		#5

Dear Marco,

              Thanks for your prompt reply. I need some clarification regarding event generation of as follows:

import model loop_qcd_qed_sm
generate e- e+ > 4j aEW<=4 [QCD]

It generates the events with some warning as we pointed out earlier.
My question is whether this results are reliable or not?

Here is my proc_card:

#****************
set group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
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~
set stdout_level DEBUG
import model loop_qcd_qed_sm
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 22 # pass to 5 flavors
define j = p
set loop_optimized_output True
set include_lepton_initiated_processes True
set complex_mass_scheme True
set max_npoint_for_channel 99
set max_t_for_channel 99
define j = g u u~ d d~ s s~ c c~
generate e+ e- > 4j aEW<=4 [QCD]

output ee4jNLOQCDsm

MG5_aMC>launch

The following switches determine which programs are run:
/====================================================================\
| 1. Type of perturbative computation order = NLO |
| 2. No MC@[N]LO matching / event generation fixed_order = OFF |
 \====================================================================/
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]
>
 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.
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat | #### Here I choose PYTHIA8 parton-shower ####
 \------------------------------------------------------------/

##### After that I got these messages ########

INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 4 cores
INFO: Compiling P0_epem_uuuxux...
INFO: Compiling P0_epem_uddxux...
INFO: Compiling P0_epem_suuxsx...
INFO: Compiling P0_epem_sudxcx...
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
STOP 1
INFO: P0_epem_sudxcx done.
INFO: Compiling P0_epem_cuuxcx...
STOP 1
STOP 1
STOP 1
INFO: P0_epem_suuxsx done.
INFO: Compiling P0_epem_dddxdx...
STOP 1
STOP 1
INFO: P0_epem_uddxux done.
INFO: Compiling P0_epem_sddxsx...
STOP 1
STOP 1
STOP 1
INFO: P0_epem_uuuxux done.
INFO: Compiling P0_epem_cduxsx...
STOP 1
STOP 1
INFO: P0_epem_cuuxcx done.
INFO: Compiling P0_epem_cddxcx...
STOP 1
STOP 1
INFO: P0_epem_cduxsx done.
STOP 1
INFO: P0_epem_sddxsx done.
STOP 1
STOP 1
INFO: P0_epem_dddxdx done.
INFO: P0_epem_cddxcx done.
INFO: Checking test output:
INFO: P0_epem_uuuxux
INFO: Result for test_ME:
INFO: Passed.
DEBUG:

INFO: Result for test_MC:
INFO: Passed.
DEBUG:
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 4 cores
INFO: Cleaning previous results
INFO: Generating events without running the shower.
WARNING: No integration channels found for contribution P0_epem_uuuxux
WARNING: No integration channels found for contribution P0_epem_uddxux
WARNING: No integration channels found for contribution P0_epem_suuxsx
WARNING: No integration channels found for contribution P0_epem_cuuxcx
WARNING: No integration channels found for contribution P0_epem_dddxdx
WARNING: No integration channels found for contribution P0_epem_sddxsx
WARNING: No integration channels found for contribution P0_epem_cddxcx
INFO: Setting up grids
INFO: Idle: 18, Running: 4, Completed: 0 [ current time: 18h07 ]
INFO: Idle: 17, Running: 4, Completed: 1 [ 36.6s ]
INFO: Idle: 16, Running: 4, Completed: 2 [ 38.4s ]
INFO: Idle: 15, Running: 4, Completed: 3 [ 1m 10s ]
INFO: Idle: 14, Running: 4, Completed: 4 [ 1m 19s ]
INFO: Idle: 13, Running: 4, Completed: 5 [ 1m 22s ]
INFO: Idle: 12, Running: 4, Completed: 6 [ 1m 41s ]
INFO: Idle: 11, Running: 4, Completed: 7 [ 1m 44s ]
INFO: Idle: 10, Running: 4, Completed: 8 [ 1m 48s ]
INFO: Idle: 9, Running: 4, Completed: 9 [ 2m 10s ]
INFO: Idle: 8, Running: 4, Completed: 10 [ 2m 21s ]
INFO: Idle: 7, Running: 4, Completed: 11 [ 2m 27s ]
INFO: Idle: 6, Running: 4, Completed: 12 [ 2m 44s ]
INFO: Idle: 5, Running: 4, Completed: 13 [ 2m 48s ]
INFO: Idle: 4, Running: 4, Completed: 14 [ 3m 3s ]
INFO: Idle: 3, Running: 4, Completed: 15 [ 3m 6s ]
INFO: Idle: 2, Running: 4, Completed: 16 [ 3m 11s ]
INFO: Idle: 1, Running: 4, Completed: 17 [ 3m 13s ]
INFO: Idle: 0, Running: 4, Completed: 18 [ 3m 36s ]
INFO: Idle: 0, Running: 3, Completed: 19 [ 3m 42s ]
INFO: Idle: 0, Running: 2, Completed: 20 [ 3m 46s ]
INFO: Idle: 0, Running: 1, Completed: 21 [ 4m 4s ]
INFO: Idle: 0, Running: 0, Completed: 22 [ 4m 5s ]
sum of cpu time of last step: 0 second
INFO: Determining the number of unweighted events per channel

      Intermediate results:
      Random seed: 33
      Total cross section: 4.135e-01 +- 2.4e-02 pb
      Total abs(cross section): 2.102e+00 +- 2.5e-02 pb

INFO: Computing upper envelope
INFO: Idle: 18, Running: 4, Completed: 0 [ current time: 18h11 ]
INFO: Idle: 17, Running: 4, Completed: 1 [ 14.8s ]
INFO: Idle: 16, Running: 4, Completed: 2 [ 15.1s ]
INFO: Idle: 15, Running: 4, Completed: 3 [ 15.6s ]
INFO: Idle: 14, Running: 4, Completed: 4 [ 16.2s ]
INFO: Idle: 13, Running: 4, Completed: 5 [ 31.4s ]
INFO: Idle: 12, Running: 4, Completed: 6 [ 32.3s ]
INFO: Idle: 11, Running: 4, Completed: 7 [ 33.8s ]
INFO: Idle: 10, Running: 4, Completed: 8 [ 37.1s ]
INFO: Idle: 9, Running: 4, Completed: 9 [ 47.3s ]
INFO: Idle: 8, Running: 4, Completed: 10 [ 48.6s ]
INFO: Idle: 7, Running: 4, Completed: 11 [ 53.5s ]
INFO: Idle: 6, Running: 4, Completed: 12 [ 54.9s ]
INFO: Idle: 5, Running: 4, Completed: 13 [ 1m 3s ]
INFO: Idle: 4, Running: 4, Completed: 14 [ 1m 4s ]
INFO: Idle: 3, Running: 4, Completed: 15 [ 1m 9s ]
INFO: Idle: 2, Running: 4, Completed: 16 [ 1m 16s ]
INFO: Idle: 1, Running: 4, Completed: 17 [ 1m 21s ]
INFO: Idle: 0, Running: 4, Completed: 18 [ 1m 22s ]
INFO: Idle: 0, Running: 3, Completed: 19 [ 1m 25s ]
INFO: Idle: 0, Running: 2, Completed: 20 [ 1m 30s ]
INFO: Idle: 0, Running: 1, Completed: 21 [ 1m 34s ]
INFO: Idle: 0, Running: 0, Completed: 22 [ 1m 36s ]
sum of cpu time of last step: 0 second
INFO: Updating the number of unweighted events per channel

      Intermediate results:
      Random seed: 33
      Total cross section: 4.745e-01 +- 1.8e-02 pb
      Total abs(cross section): 2.179e+00 +- 1.8e-02 pb

INFO: Generating events
INFO: Idle: 0, Running: 0, Completed: 6 [ current time: 18h13 ]
sum of cpu time of last step: 0 second
INFO: Doing reweight
INFO: Idle: 0, Running: 0, Completed: 6 [ current time: 18h13 ]
INFO: Collecting events
DEBUG:
  Number of loop ME evaluations (by MadLoop): 34646
    Stability unknown: 0
    Stable PS point: 34646
    Unstable PS point (and rescued): 0
    Unstable PS point (and not rescued): 0
    Only double precision used: 34646
    Quadruple precision used: 0
    Initialization phase-space points: 0
    Reduction methods used:
      > Ninja (double precision) 34646

  Statistics on virtual integration optimization :
  Could not find statistics on the integration optimization.

 Inclusive timing profile :
    Overall slowest channel 0:01:41 (P0_epem_sudxcx/GF1.0 [step 0])
    Average channel running time 0:00:25
    Aggregated total running time 0:20:51
  Timing profile for <Other_tasks> :
    Overall fraction of time 25.284 %
    Largest fraction of time 49.243 % (P0_epem_sudxcx/GF5.0 [step 2])
    Smallest fraction of time 9.935 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Virtuals> :
    Overall fraction of time 20.821 %
    Largest fraction of time 81.665 % (P0_epem_cduxsx/GF6.0 [step 2])
    Smallest fraction of time 0.000 % (P0_epem_cduxsx/GF11.0 [step 2])
  Timing profile for <Nbody_prefactor> :
    Overall fraction of time 17.404 %
    Largest fraction of time 22.930 % (P0_epem_sudxcx/GF11.0 [step 1])
    Smallest fraction of time 1.892 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Reals_evaluation> :
    Overall fraction of time 8.847 %
    Largest fraction of time 10.719 % (P0_epem_sudxcx/GF9.0 [step 1])
    Smallest fraction of time 0.883 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Applying_cuts> :
    Overall fraction of time 8.789 %
    Largest fraction of time 10.719 % (P0_epem_sudxcx/GF7.0 [step 1])
    Smallest fraction of time 0.807 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <PS_Generation> :
    Overall fraction of time 4.912 %
    Largest fraction of time 5.915 % (P0_epem_sudxcx/GF7.0 [step 1])
    Smallest fraction of time 0.539 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Counter_terms> :
    Overall fraction of time 3.609 %
    Largest fraction of time 5.184 % (P0_epem_sudxcx/GF11.0 [step 1])
    Smallest fraction of time 0.348 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Integrated_CT> :
    Overall fraction of time 3.338 %
    Largest fraction of time 4.451 % (P0_epem_sudxcx/GF11.0 [step 1])
    Smallest fraction of time 0.446 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <AlphaS_dependencies> :
    Overall fraction of time 3.267 %
    Largest fraction of time 4.405 % (P0_epem_sudxcx/GF11.0 [step 1])
    Smallest fraction of time 0.364 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <MCsubtraction> :
    Overall fraction of time 1.306 %
    Largest fraction of time 1.792 % (P0_epem_sudxcx/GF11.0 [step 1])
    Smallest fraction of time 0.124 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Born> :
    Overall fraction of time 0.713 %
    Largest fraction of time 0.956 % (P0_epem_sudxcx/GF11.0 [step 1])
    Smallest fraction of time 0.064 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Adding_alphas_pdf> :
    Overall fraction of time 0.668 %
    Largest fraction of time 0.841 % (P0_epem_cduxsx/GF11.0 [step 1])
    Smallest fraction of time 0.087 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <N1body_prefactor> :
    Overall fraction of time 0.526 %
    Largest fraction of time 0.651 % (P0_epem_cduxsx/GF10.0 [step 1])
    Smallest fraction of time 0.043 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Sum_ident_contr> :
    Overall fraction of time 0.502 %
    Largest fraction of time 0.616 % (P0_epem_cduxsx/GF9.0 [step 1])
    Smallest fraction of time 0.061 % (P0_epem_cduxsx/GF6.0 [step 2])
  Timing profile for <Write_events> :
    Overall fraction of time 0.014 %
    Largest fraction of time 13.725 % (P0_epem_sudxcx/GF5.0 [step 2])
    Smallest fraction of time 0.000 % (P0_epem_cduxsx/GF1.0 [step 0])
  Timing profile for <Pick_unwgt> :
    Overall fraction of time 0.000 %
    Largest fraction of time 0.017 % (P0_epem_cduxsx/GF11.0 [step 2])
    Smallest fraction of time 0.000 % (P0_epem_cduxsx/GF1.0 [step 0])
  Zero time record for FxFx_cluster.
  Zero time record for Reweight_scale.
  Zero time record for Reweight_pdf.
  Zero time record for Filling_plots.
INFO:
   --------------------------------------------------------------
      Summary:
      Process e+ e- > 4j aEW<=4 [QCD]
      Run at l-l collider (500.0 + 500.0 GeV)
      Number of events generated: 100
      Total cross section: 4.737e-01 +- 1.8e-02 pb
   --------------------------------------------------------------

INFO: The /home/swapan/hep_packages/MG5_aMC_v3.5.4/MG5_aMC_v3_5_4/ee4jNLOQCDsm/Events/run_01/events.lhe.gz file has been generated.

INFO: Events generated
reweight -from_cards
decay_events -from_cards
WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical. Please, shower the Les Houches events before using them for physics analyses.
INFO:
quit
INFO:

Revision history for this message
marco zaro (marco-zaro) said : 		#6

Hi,
to answer your question, they are not reliable. What you are generating is not only QCD corrections to WW decay, there is a whole extra lot in it (including EW corrections on different LO contributions).
this is explained in appendix D of https://inspirehep.net/literature/1670211 <https://inspirehep.net/literature/1670211>
Best wishes,

Marco

> On 23 May 2024, at 16:35, Swapan Kumar Majhi <email address hidden> wrote:
>
> Question #816013 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/816013
>
> Swapan Kumar Majhi posted a new comment:
> Dear Marco,
>
> Thanks for your prompt reply. I need some clarification
> regarding event generation of as follows:
>
> import model loop_qcd_qed_sm
> generate e- e+ > 4j aEW<=4 [QCD]
>
> It generates the events with some warning as we pointed out earlier.
> My question is whether this results are reliable or not?
>
> Here is my proc_card:
>
> #****************
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set low_mem_multicore_nlo_generation False
> set complex_mass_scheme False
> set include_lepton_initiated_processes False
> set gauge unitary
> set loop_optimized_output True
> set loop_color_flows False
> set max_npoint_for_channel 0
> set default_unset_couplings 99
> set max_t_for_channel 99
> set zerowidth_tchannel True
> set nlo_mixed_expansion True
> 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~
> set stdout_level DEBUG
> import model loop_qcd_qed_sm
> define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 22 # pass to 5 flavors
> define j = p
> set loop_optimized_output True
> set include_lepton_initiated_processes True
> set complex_mass_scheme True
> set max_npoint_for_channel 99
> set max_t_for_channel 99
> define j = g u u~ d d~ s s~ c c~
> generate e+ e- > 4j aEW<=4 [QCD]
>
> output ee4jNLOQCDsm
>
> MG5_aMC>launch
>
> The following switches determine which programs are run:
> /====================================================================\
> | 1. Type of perturbative computation order = NLO |
> | 2. No MC@[N]LO matching / event generation fixed_order = OFF |
> \====================================================================/
> 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]
>>
> 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.
> Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
> Lepton-lepton collisions: ignoring PDF related parameters in the run_card.dat (pdlabel, lhaid, reweight_pdf, ...)
> Do you want to edit a card (press enter to bypass editing)?
> /------------------------------------------------------------\
> | 1. param : param_card.dat |
> | 2. run : run_card.dat | #### Here I choose PYTHIA8 parton-shower ####
> \------------------------------------------------------------/
>
> ##### After that I got these messages ########
>
> INFO: ...done, continuing with P* directories
> INFO: Compiling directories...
> INFO: Compiling on 4 cores
> INFO: Compiling P0_epem_uuuxux...
> INFO: Compiling P0_epem_uddxux...
> INFO: Compiling P0_epem_suuxsx...
> INFO: Compiling P0_epem_sudxcx...
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> INFO: P0_epem_sudxcx done.
> INFO: Compiling P0_epem_cuuxcx...
> STOP 1
> STOP 1
> STOP 1
> INFO: P0_epem_suuxsx done.
> INFO: Compiling P0_epem_dddxdx...
> STOP 1
> STOP 1
> INFO: P0_epem_uddxux done.
> INFO: Compiling P0_epem_sddxsx...
> STOP 1
> STOP 1
> STOP 1
> INFO: P0_epem_uuuxux done.
> INFO: Compiling P0_epem_cduxsx...
> STOP 1
> STOP 1
> INFO: P0_epem_cuuxcx done.
> INFO: Compiling P0_epem_cddxcx...
> STOP 1
> STOP 1
> INFO: P0_epem_cduxsx done.
> STOP 1
> INFO: P0_epem_sddxsx done.
> STOP 1
> STOP 1
> INFO: P0_epem_dddxdx done.
> INFO: P0_epem_cddxcx done.
> INFO: Checking test output:
> INFO: P0_epem_uuuxux
> INFO: Result for test_ME:
> INFO: Passed.
> DEBUG:
>
> INFO: Result for test_MC:
> INFO: Passed.
> DEBUG:
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> INFO: Starting run
> INFO: Using 4 cores
> INFO: Cleaning previous results
> INFO: Generating events without running the shower.
> WARNING: No integration channels found for contribution P0_epem_uuuxux
> WARNING: No integration channels found for contribution P0_epem_uddxux
> WARNING: No integration channels found for contribution P0_epem_suuxsx
> WARNING: No integration channels found for contribution P0_epem_cuuxcx
> WARNING: No integration channels found for contribution P0_epem_dddxdx
> WARNING: No integration channels found for contribution P0_epem_sddxsx
> WARNING: No integration channels found for contribution P0_epem_cddxcx
> INFO: Setting up grids
> INFO: Idle: 18, Running: 4, Completed: 0 [ current time: 18h07 ]
> INFO: Idle: 17, Running: 4, Completed: 1 [ 36.6s ]
> INFO: Idle: 16, Running: 4, Completed: 2 [ 38.4s ]
> INFO: Idle: 15, Running: 4, Completed: 3 [ 1m 10s ]
> INFO: Idle: 14, Running: 4, Completed: 4 [ 1m 19s ]
> INFO: Idle: 13, Running: 4, Completed: 5 [ 1m 22s ]
> INFO: Idle: 12, Running: 4, Completed: 6 [ 1m 41s ]
> INFO: Idle: 11, Running: 4, Completed: 7 [ 1m 44s ]
> INFO: Idle: 10, Running: 4, Completed: 8 [ 1m 48s ]
> INFO: Idle: 9, Running: 4, Completed: 9 [ 2m 10s ]
> INFO: Idle: 8, Running: 4, Completed: 10 [ 2m 21s ]
> INFO: Idle: 7, Running: 4, Completed: 11 [ 2m 27s ]
> INFO: Idle: 6, Running: 4, Completed: 12 [ 2m 44s ]
> INFO: Idle: 5, Running: 4, Completed: 13 [ 2m 48s ]
> INFO: Idle: 4, Running: 4, Completed: 14 [ 3m 3s ]
> INFO: Idle: 3, Running: 4, Completed: 15 [ 3m 6s ]
> INFO: Idle: 2, Running: 4, Completed: 16 [ 3m 11s ]
> INFO: Idle: 1, Running: 4, Completed: 17 [ 3m 13s ]
> INFO: Idle: 0, Running: 4, Completed: 18 [ 3m 36s ]
> INFO: Idle: 0, Running: 3, Completed: 19 [ 3m 42s ]
> INFO: Idle: 0, Running: 2, Completed: 20 [ 3m 46s ]
> INFO: Idle: 0, Running: 1, Completed: 21 [ 4m 4s ]
> INFO: Idle: 0, Running: 0, Completed: 22 [ 4m 5s ]
> sum of cpu time of last step: 0 second
> INFO: Determining the number of unweighted events per channel
>
> Intermediate results:
> Random seed: 33
> Total cross section: 4.135e-01 +- 2.4e-02 pb
> Total abs(cross section): 2.102e+00 +- 2.5e-02 pb
>
> INFO: Computing upper envelope
> INFO: Idle: 18, Running: 4, Completed: 0 [ current time: 18h11 ]
> INFO: Idle: 17, Running: 4, Completed: 1 [ 14.8s ]
> INFO: Idle: 16, Running: 4, Completed: 2 [ 15.1s ]
> INFO: Idle: 15, Running: 4, Completed: 3 [ 15.6s ]
> INFO: Idle: 14, Running: 4, Completed: 4 [ 16.2s ]
> INFO: Idle: 13, Running: 4, Completed: 5 [ 31.4s ]
> INFO: Idle: 12, Running: 4, Completed: 6 [ 32.3s ]
> INFO: Idle: 11, Running: 4, Completed: 7 [ 33.8s ]
> INFO: Idle: 10, Running: 4, Completed: 8 [ 37.1s ]
> INFO: Idle: 9, Running: 4, Completed: 9 [ 47.3s ]
> INFO: Idle: 8, Running: 4, Completed: 10 [ 48.6s ]
> INFO: Idle: 7, Running: 4, Completed: 11 [ 53.5s ]
> INFO: Idle: 6, Running: 4, Completed: 12 [ 54.9s ]
> INFO: Idle: 5, Running: 4, Completed: 13 [ 1m 3s ]
> INFO: Idle: 4, Running: 4, Completed: 14 [ 1m 4s ]
> INFO: Idle: 3, Running: 4, Completed: 15 [ 1m 9s ]
> INFO: Idle: 2, Running: 4, Completed: 16 [ 1m 16s ]
> INFO: Idle: 1, Running: 4, Completed: 17 [ 1m 21s ]
> INFO: Idle: 0, Running: 4, Completed: 18 [ 1m 22s ]
> INFO: Idle: 0, Running: 3, Completed: 19 [ 1m 25s ]
> INFO: Idle: 0, Running: 2, Completed: 20 [ 1m 30s ]
> INFO: Idle: 0, Running: 1, Completed: 21 [ 1m 34s ]
> INFO: Idle: 0, Running: 0, Completed: 22 [ 1m 36s ]
> sum of cpu time of last step: 0 second
> INFO: Updating the number of unweighted events per channel
>
> Intermediate results:
> Random seed: 33
> Total cross section: 4.745e-01 +- 1.8e-02 pb
> Total abs(cross section): 2.179e+00 +- 1.8e-02 pb
>
>
> INFO: Generating events
> INFO: Idle: 0, Running: 0, Completed: 6 [ current time: 18h13 ]
> sum of cpu time of last step: 0 second
> INFO: Doing reweight
> INFO: Idle: 0, Running: 0, Completed: 6 [ current time: 18h13 ]
> INFO: Collecting events
> DEBUG:
> Number of loop ME evaluations (by MadLoop): 34646
> Stability unknown: 0
> Stable PS point: 34646
> Unstable PS point (and rescued): 0
> Unstable PS point (and not rescued): 0
> Only double precision used: 34646
> Quadruple precision used: 0
> Initialization phase-space points: 0
> Reduction methods used:
>> Ninja (double precision) 34646
>
> Statistics on virtual integration optimization :
> Could not find statistics on the integration optimization.
>
> Inclusive timing profile :
> Overall slowest channel 0:01:41 (P0_epem_sudxcx/GF1.0 [step 0])
> Average channel running time 0:00:25
> Aggregated total running time 0:20:51
> Timing profile for <Other_tasks> :
> Overall fraction of time 25.284 %
> Largest fraction of time 49.243 % (P0_epem_sudxcx/GF5.0 [step 2])
> Smallest fraction of time 9.935 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Virtuals> :
> Overall fraction of time 20.821 %
> Largest fraction of time 81.665 % (P0_epem_cduxsx/GF6.0 [step 2])
> Smallest fraction of time 0.000 % (P0_epem_cduxsx/GF11.0 [step 2])
> Timing profile for <Nbody_prefactor> :
> Overall fraction of time 17.404 %
> Largest fraction of time 22.930 % (P0_epem_sudxcx/GF11.0 [step 1])
> Smallest fraction of time 1.892 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Reals_evaluation> :
> Overall fraction of time 8.847 %
> Largest fraction of time 10.719 % (P0_epem_sudxcx/GF9.0 [step 1])
> Smallest fraction of time 0.883 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Applying_cuts> :
> Overall fraction of time 8.789 %
> Largest fraction of time 10.719 % (P0_epem_sudxcx/GF7.0 [step 1])
> Smallest fraction of time 0.807 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <PS_Generation> :
> Overall fraction of time 4.912 %
> Largest fraction of time 5.915 % (P0_epem_sudxcx/GF7.0 [step 1])
> Smallest fraction of time 0.539 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Counter_terms> :
> Overall fraction of time 3.609 %
> Largest fraction of time 5.184 % (P0_epem_sudxcx/GF11.0 [step 1])
> Smallest fraction of time 0.348 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Integrated_CT> :
> Overall fraction of time 3.338 %
> Largest fraction of time 4.451 % (P0_epem_sudxcx/GF11.0 [step 1])
> Smallest fraction of time 0.446 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <AlphaS_dependencies> :
> Overall fraction of time 3.267 %
> Largest fraction of time 4.405 % (P0_epem_sudxcx/GF11.0 [step 1])
> Smallest fraction of time 0.364 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <MCsubtraction> :
> Overall fraction of time 1.306 %
> Largest fraction of time 1.792 % (P0_epem_sudxcx/GF11.0 [step 1])
> Smallest fraction of time 0.124 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Born> :
> Overall fraction of time 0.713 %
> Largest fraction of time 0.956 % (P0_epem_sudxcx/GF11.0 [step 1])
> Smallest fraction of time 0.064 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Adding_alphas_pdf> :
> Overall fraction of time 0.668 %
> Largest fraction of time 0.841 % (P0_epem_cduxsx/GF11.0 [step 1])
> Smallest fraction of time 0.087 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <N1body_prefactor> :
> Overall fraction of time 0.526 %
> Largest fraction of time 0.651 % (P0_epem_cduxsx/GF10.0 [step 1])
> Smallest fraction of time 0.043 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Sum_ident_contr> :
> Overall fraction of time 0.502 %
> Largest fraction of time 0.616 % (P0_epem_cduxsx/GF9.0 [step 1])
> Smallest fraction of time 0.061 % (P0_epem_cduxsx/GF6.0 [step 2])
> Timing profile for <Write_events> :
> Overall fraction of time 0.014 %
> Largest fraction of time 13.725 % (P0_epem_sudxcx/GF5.0 [step 2])
> Smallest fraction of time 0.000 % (P0_epem_cduxsx/GF1.0 [step 0])
> Timing profile for <Pick_unwgt> :
> Overall fraction of time 0.000 %
> Largest fraction of time 0.017 % (P0_epem_cduxsx/GF11.0 [step 2])
> Smallest fraction of time 0.000 % (P0_epem_cduxsx/GF1.0 [step 0])
> Zero time record for FxFx_cluster.
> Zero time record for Reweight_scale.
> Zero time record for Reweight_pdf.
> Zero time record for Filling_plots.
> INFO:
> --------------------------------------------------------------
> Summary:
> Process e+ e- > 4j aEW<=4 [QCD]
> Run at l-l collider (500.0 + 500.0 GeV)
> Number of events generated: 100
> Total cross section: 4.737e-01 +- 1.8e-02 pb
> --------------------------------------------------------------
>
> INFO: The /home/swapan/hep_packages/MG5_aMC_v3.5.4/MG5_aMC_v3_5_4/ee4jNLOQCDsm/Events/run_01/events.lhe.gz file has been generated.
>
> INFO: Events generated
> reweight -from_cards
> decay_events -from_cards
> WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical. Please, shower the Les Houches events before using them for physics analyses.
> INFO:
> quit
> INFO:
>
> --
> You received this question notification because you are subscribed to
> the question.

Can you help with this problem?

Provide an answer of your own, or ask Swapan Kumar Majhi for more information if necessary.

To post a message you must log in.