MadGraph5_aMC@NLO

Different result between version 3.3.1 and 3.5.1 for some processes

Asked by dong joo kim

Dear Experts

I found an issue when calculating the same NLO process with two different versions of MadGraph5 3.3.1, 3.3.2 and 3.5.1.
As i understand, it should be the same but the result was not same, there was almost 40 times difference at most.

here is a detail,

(define wpm = w+ w-)
i could check only for 5 processes,(p p > z [QCD]/ p p > wpm [QCD] / p p > t t~ b b~ [QCD] / p p > z b b~ j [QCD] / p p > h b b~ j [QCD])
and i changed parameter in 'run_card.dat' like this(for 3.3.1, 3.3.2 and 3.5.1).

ebeam1, ebeam2 = 7000.0
pdlabel = lhapdf
lhaid = 315800
ptj = 30.0
etaj = 4.0
dynamical_scale_choice = 3

for z,

3.3.1 version result > NLO : 5.797 x 10^4 pb
3.5.1 version result > NLO : 2.64 x 10^6 pb

for wpm,

3.3.1 version result > NLO : 1.9 x 10^5 pb
3.5.1 version result > NLO : 7.75 x 10^6 pb

for t t~ b b~,

3.3.2 version result > NLO : 1.885 x 10^1 pb
3.5.1 version result > NLO : 9.436 x 10^1 pb

for z b b~ j,

3.3.1 version result > NLO : 3.582 x 10^2 pb
3.5.1 version result > NLO : error : poles do not cancel

for h b b~ j,

3.3.1 version result > NLO : 1.111 x 10^-1 pb
3.5.1 version result > NLO : error : poles do not cancel

i also checked LO result for 5 same processes,
but for LO results, not like NLO results, there was only few percent difference.

for LO case, i used these condition,

ebeam1, ebeam2 = 7000.0
pdlabel = lhapdf
lhaid = 315000
ptj = 30.0
etaj = 4.0
dynamical_scale_choice = 3

and results are

for z,

3.5.1 version result > LO : 4.213 x 10^4 pb

for wpm,

3.5.1 version result > LO : 1.390 x 10^5 pb

for t t~ b b~,

3.3.2 version result > LO : 7.500 pb
3.5.1 version result > LO : 7.503 pb

for z b b~ j,

3.3.1 version result > LO : 1.891 x 10^2 pb

for h b b~ j,

3.3.1 version result > LO : 8.504 x 10^-2 pb

Although i didn't posted another version results, for the rest processes(z, wpm, z b b~j, h b b~j)
as i checked, there were only few percent difference like t t~ b b~ process as you can see.

that is all for details, LO calculation is not a problem i guess, but NLO.
so can you tell me what should i check or do to fix this issue? or is there some mistake i made?

please Let me know.

*P.S
Above calculation was done in fixed_order condition, so i checked with shower=on with pythia8, in 3.5.1 version
for p p > z process, both NLO, LO, result was same as previous result (LO : 4.216E+04 pb / NLO : 2.617E+06 pb)
which means k-factor(NLO/LO ratio) was "62.07".... too big i guess,,, but for 3.3.1 version , k-factor was 1.371

thank you for your attention
sincerely.

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

Hi,

I have tested the following script:
generate p p > z [QCD]
output
launch
shower=PY8
set ebeam 7k
set pdlabel lhapdf
set lhaid 315800
set ptj = 30.0
set etaj = 4.0
set dynamical_scale_choice = 3

For different version of the code:
First with the Long Term Stable (2.9.16): 5.76(1)e4
For 3.5.2 (not yet released): 5.80(1)e4
for 3.5.0/3.5.1: crash due to a bug in the support of 2>1 matrix-element
For 3.4.2: 5.82(1)e4
For 3.4.1: 5.82(1)e4

So here I do not face the same issue as you and reproduce your number for 3.3 version.
I do face a crash for 3.5.0 and 3.5.1 and not a wrong cross-section. Did you apply a patch to avoid the crash?

For tt~bb~ this is a problematic process at NLO, since it does not make sense in a five flavor computation (since you request to have b)
but this is a process that is not possible to do in four flavor within dedicated diagram substraction (which requires MadSTR PLUGIN which is not yet compatible with version 3.X of the code). So that either does not make sense or fail with all version that you are trying (at least I would not try to produce number for that one).

For the others, are you using a four flavor model?

Cheers,

Olivier

Revision history for this message
dong joo kim (djkdj) said : 		#2

thank you for your reply.

yes i used 4-flavor model for all processes above.

for your information,
although i didn't mentioned, for z process, i checked 3.4.2 version and i got the same result as you and 3.3.1 version
(almost 5.7~5.8 x 10^4 pb).

the main reason i asked is that i got different result only in 3.5.1 version.
(4-body processes like tt~bb~ and z b b~ j can be discussed later, and somehow i understand your point)

i'm now posting p p > z [QCD] log just in case.
can you tell me which is the crash you mentioned?
and also i didn't apply any patch. i just installed pythia8, lhapdf6 and others to run NLO calculation in Madgraph5

i used
generate p p > z [QCD]
output
launch
shower=PY8
set ebeam 7k
set pdlabel lhapdf
set lhaid 315800
and log is like.

************************************************************
* *
* 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.1 2023-07-11 *
* *
* 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
None 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

None 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/kdj/MG5_aMC_v3_5_1/HEPTools/lib
set collier to /home/kdj/MG5_aMC_v3_5_1/HEPTools/lib
set lhapdf to /home/kdj/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config
None 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

Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
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/kdj/MG5_aMC_v3_5_1/NLO_2.txt
The import format was not given, so we guess it as command
generate p p > z [QCD]
No model currently active, so we import the Standard Model
INFO: Restrict model sm with file models/sm/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~
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
The current model sm does not allow to generate loop corrections of type ['QCD'].
MG5_aMC now loads 'loop_sm'.
 import model loop_sm
INFO: Restrict model loop_sm with file models/loop_sm/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 gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ b t b~ t~ z w+ h w- ta- ta+
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED
Setting the born squared orders automatically in the process definition to QED^2<=2 QCD=0.
If this is not what you need, please regenerate with the correct orders.
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: u u~ > z QCD^2=2 QED^2=2 (1 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: c c~ > z QCD^2=2 QED^2=2 (2 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: d d~ > z QCD^2=2 QED^2=2 (3 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: s s~ > z QCD^2=2 QED^2=2 (4 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: u~ u > z QCD^2=2 QED^2=2 (5 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: c~ c > z QCD^2=2 QED^2=2 (6 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: d~ d > z QCD^2=2 QED^2=2 (7 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: s~ s > z QCD^2=2 QED^2=2 (8 / 8)
INFO: Generating virtual matrix element with MadLoop for process: u u~ > z QCD^2=2 QED^2=2 (1 / 8)
INFO: Generating virtual matrix element with MadLoop for process: c c~ > z QCD^2=2 QED^2=2 (2 / 8)
INFO: Generating virtual matrix element with MadLoop for process: d d~ > z QCD^2=2 QED^2=2 (3 / 8)
INFO: Generating virtual matrix element with MadLoop for process: s s~ > z QCD^2=2 QED^2=2 (4 / 8)
INFO: Generating virtual matrix element with MadLoop for process: u~ u > z QCD^2=2 QED^2=2 (5 / 8)
INFO: Generating virtual matrix element with MadLoop for process: c~ c > z QCD^2=2 QED^2=2 (6 / 8)
INFO: Generating virtual matrix element with MadLoop for process: d~ d > z QCD^2=2 QED^2=2 (7 / 8)
INFO: Generating virtual matrix element with MadLoop for process: s~ s > z QCD^2=2 QED^2=2 (8 / 8)
INFO: Generated 8 subprocesses with 48 real emission diagrams, 8 born diagrams and 8 virtual diagrams
output NLO_z_4f
INFO: Writing out the aMC@NLO code, using optimized Loops
INFO: directory /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f already exists.
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n]
found line : launch NLO_z_4f
This answer is not valid for current question. Keep it for next question and use here default: y
INFO: initialize a new directory: NLO_z_4f
INFO: remove old information in NLO_z_4f
INFO: Generating real emission matrix-elements...
INFO: Generating Helas calls for FKS process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (1 / 8)
INFO: Processing color information for process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: c c~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (2 / 8)
INFO: Generating Helas calls for FKS process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (3 / 8)
INFO: Reusing existing color information for process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: s s~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (4 / 8)
INFO: Generating Helas calls for FKS process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2 (5 / 8)
INFO: Processing color information for process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: c~ c > z [ all = QCD QED ] QCD^2=2 QED^2=2 (6 / 8)
INFO: Generating Helas calls for FKS process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2 (7 / 8)
INFO: Reusing existing color information for process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: s~ s > z [ all = QCD QED ] QCD^2=2 QED^2=2 (8 / 8)
INFO: ... Done
INFO: Writing files in P0_uux_z (1 / 4)
INFO: Creating files in directory V0_uux_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Writing files in P0_ddx_z (2 / 4)
INFO: Creating files in directory V0_ddx_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Writing files in P0_uxu_z (3 / 4)
INFO: Creating files in directory V0_uxu_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Writing files in P0_dxd_z (4 / 4)
INFO: Creating files in directory V0_dxd_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 16 routines in 0.441 s
save configuration file to /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/Cards/amcatnlo_configuration.txt
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/kdj/MG5_aMC_v3_5_1/NLO_z_4f/README
Run "open index.html" to see more information about this process.
launch NLO_z_4f
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.1 2023-07-11 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/Cards/amcatnlo_configuration.txt
INFO: load configuration from /home/kdj/MG5_aMC_v3_5_1/input/mg5_configuration.txt
INFO: load configuration from /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/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|PYTHIA6Q|PYTHIA6PT|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.
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|PYTHIA6Q|PYTHIA6PT|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.
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|PYTHIA6Q|PYTHIA6PT|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.
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|PYTHIA6Q|PYTHIA6PT|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.
INFO: will run in mode: aMC@NLO
INFO: modify parameter parton_shower of the run_card.dat to PYTHIA8
set ebeam1 7000
INFO: modify parameter ebeam1 of the run_card.dat to 7000.0
set ebeam2 7000
INFO: modify parameter ebeam2 of the run_card.dat to 7000.0
set pdlabel lhapdf
INFO: modify parameter pdlabel of the run_card.dat to lhapdf
set lhaid 315800
INFO: modify parameter lhaid of the run_card.dat to [315800]
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. shower : shower_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, 3, shower, enter path]
The answer to the previous question is not set in your input file
Use 0 value
INFO: modify parameter extralibs of the shower_card.dat to pythia8 z dl stdc++
INFO: modify parameter extrapaths of the shower_card.dat to ../lib /home/kdj/MG5_aMC_v3_5_1/HEPTools/pythia8//lib l,-rpath,/home/kdj/MG5_aMC_v3_5_1/HEPTools/pythia8//lib /home/kdj/MG5_aMC_v3_5_1/HEPTools/zlib/lib l,-rpath,/home/kdj/MG5_aMC_v3_5_1/HEPTools/zlib/lib
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.11800216951342939
INFO: Starting run
INFO: Compiling the code
Error: no DISPLAY environment variable specified
INFO: Using LHAPDF v6.3.0 interface for PDFs
INFO: For gauge cancellation, the width of 'Z' has been set to zero.
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 32 cores
INFO: Compiling P0_uux_z...
INFO: Compiling P0_ddx_z...
INFO: Compiling P0_uxu_z...
INFO: Compiling P0_dxd_z...
INFO: P0_uxu_z done.
INFO: P0_dxd_z done.
INFO: P0_uux_z done.
INFO: P0_ddx_z done.
INFO: Checking test output:
INFO: P0_uux_z
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_ddx_z
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_uxu_z
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_dxd_z
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: Starting run
INFO: Using 32 cores
INFO: Cleaning previous results
INFO: Doing NLO matched to parton shower
INFO: Setting up grids
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 12h41 ]
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: 2.623e+06 +- 1.4e+04 pb
      Total abs(cross section): 2.965e+06 +- 1.3e+04 pb

INFO: Computing upper envelope
INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 12h41 ]
INFO: Idle: 0, Running: 3, Completed: 1 [ 0.3s ]
INFO: Idle: 0, Running: 2, Completed: 2 [ 0.31s ]
INFO: Idle: 0, Running: 1, Completed: 3 [ 0.31s ]
INFO: Idle: 0, Running: 0, Completed: 4 [ 0.34s ]
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: 2.646e+06 +- 7.8e+03 pb
      Total abs(cross section): 2.992e+06 +- 7.4e+03 pb

INFO: Generating events
INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 12h41 ]
INFO: Idle: 0, Running: 3, Completed: 1 [ 0.45s ]
INFO: Idle: 0, Running: 2, Completed: 2 [ 0.49s ]
INFO: Idle: 0, Running: 1, Completed: 3 [ 0.56s ]
INFO: Idle: 0, Running: 0, Completed: 4 [ 0.58s ]
sum of cpu time of last step: 0 second
INFO: Doing reweight
INFO: Idle: 0, Running: 2, Completed: 2 [ current time: 12h41 ]
INFO: Idle: 0, Running: 1, Completed: 3 [ 0.068s ]
INFO: Idle: 0, Running: 0, Completed: 4 [ 0.079s ]
INFO: Collecting events
INFO:
   --------------------------------------------------------------
      Summary:
      Process p p > z [QCD]
      Run at p-p collider (7000.0 + 7000.0 GeV)
      Number of events generated: 10000
      Total cross section: 2.646e+06 +- 7.8e+03 pb
   --------------------------------------------------------------
      Scale variation (computed from LHE events):
          Dynamical_scale_choice -1 (envelope of 9 values):
              2.626e+06 pb +4.5% -8.2%
   --------------------------------------------------------------

INFO: The /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/Events/run_01/events.lhe.gz file has been generated.

INFO: Events generated
reweight -from_cards
decay_events -from_cards
INFO: Preparing MCatNLO run
INFO: Using LHAPDF v6.3.0 interface for PDFs
INFO: Compiling MCatNLO for PYTHIA8...
INFO: ... done
INFO: Showering events...
INFO: (Running in /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/MCatNLO/RUN_PYTHIA8_1)
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 12h41 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 57.8s ]
INFO: The file /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/Events/run_01/events_PYTHIA8_0.hepmc.gz has been generated.
It contains showered and hadronized events in the HEPMC format obtained by showering the parton-level event file /home/kdj/MG5_aMC_v3_5_1/NLO_z_4f/Events/run_01/events.lhe.gz with PYTHIA8
INFO: Run complete
INFO:
quit
INFO:
GDBus.Error:org.freedesktop.DBus.Error.ServiceUnknown: The name org.freedesktop.Notifications was not provided by any .service files
quit

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#3

Hi,

I do not have the message of the crash right now, and my computer is busy to run all the validation to relase 2.9.17 (and then will move to run all the check for 3.5.2).

Looks like that your compiler is less secure than mine and does not correctly detect the issue in the code (the one introduced in 3.5.0) and that such bug is likely the reason of your wrong cross-section in this case. So this is even more serious than estimated so far since it is not always a crash but only a wrong reported cross-section.

I hope to be able to release 2.9.17 today and 3.5.2 early next week.

Thanks,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask dong joo kim for more information if necessary.

To post a message you must log in.