HC_NLO_X0_UFO in MadGraph 3.X
Dear experts,
I would like to ask a follow-up question on the topic I posted here:
https:/
I am still aiming to simulate t-channel tH production at NLO, with a separation of quadratic and interference terms in order to model the process continuously in the parameter space of CP-even and CP-odd top Yukawa coupling modifiers.
In the last thread, I mentioned the problem of the pole cancellation error for the interference terms at NLO when using version 2.9.18 and it was suggested to use version 3.X and try with that.
I tried your suggestion. However, using version 3.5.3, I do not manage to get any of the Higgs production processes described in the HC_NLO_X0_UFO model’s readme simulated.
Trying the tHq process with the syntax from the model readme (without specifying coupling orders), it does not find the right diagrams (It sets the born squared orders automatically to “QED^2<=6 QCD^2<=2 QNP=0”, while QNP is needed for the X0-WW coupling).
So, I assume I need to run it via:
`import model HC_NLO_X0_UFO
generate p p > t x0 b~ j $$ w+ w- [QCD] QED^2<=6 QCD^2<=2 QNP^2<=2`
However, in both cases it seems to freeze in the compiling step:
INFO: Compiling directories...
INFO: Compiling on 64 cores
INFO: Compiling P0_gu_x0tbxd...
INFO: Compiling P0_gdx_x0tbxux...
INFO: Compiling P0_ug_x0tbxd...
INFO: Compiling P0_dxg_x0tbxux…
Also, trying to run the simpler process ggH does not work in my setup with the HC model and version 3.5.3:
`import model HC_NLO_X0_UFO-heft
generate p p > x0 / t [QCD] QCD^2<=2 QNP^2<=2`
(.....)
INFO: Compiling directories...
INFO: Compiling on 96 cores
INFO: Compiling P0_gg_x0_no_t...
INFO: P0_gg_x0_no_t done.
INFO: Checking test output:
INFO: P0_gg_x0_no_t
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
Error detected in "launch auto "
write debug file /.automount/
If you need help with this issue, please, contact us on https:/
str : Poles do not cancel, run cannot continue
quit
So I would like to know your impression on these issues:
Is it possibly due to the model generally being incompatible with version 3.X?
Is there anything wrong with the syntax I used? Possibly anything odd in my environment?
I attach the full output I receive for the tHq process here.
I would greatly appreciate to hear your opinion on this.
Thank you very much for your help.
Cheers
Florian
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- marco zaro Edit question
- Solved by:
- Florian Mausolf
- Solved:
- Last query:
- Last reply:
|
Revision history for this message
|
#1 |
Full output for tHq is this:
./mg5_aMC
launch*
* *
* 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.3 2023-12-23 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* and *
* http://
* *
* 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/
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/
MG5_aMC> set fastjet /PATH/TO/
eMELA-config does not seem to correspond to a valid eMELA-config executable.
Please set the 'fastjet'variable to the full (absolute) /PATH/TO/
MG5_aMC> set eMELA /PATH/TO/
set ninja to /.automount/
set collier to /.automount/
set lhapdf to lhapdf-config
None does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/
Note that you can still compile and run aMC@NLO with the built-in PDFs
MG5_aMC> set lhapdf /PATH/TO/
Using default text editor "vi". Set another one in ./input/
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
Loading default model: sm
INFO: Restrict model sm with file ../models/
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+
MG5_aMC>launch tHq_4f_NLO_MH125_SM
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.3 2023-12-23 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /.automount/
INFO: load configuration from /.automount/
INFO: load configuration from /.automount/
Using default text editor "vi". Set another one in ./input/
launch auto
The following switches determine which programs are run:
/================== Description =======
| 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.[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.
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
\------
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path][90s to answer]
>
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
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
/usr/sbin/
/usr/sbin/
XPCOMGlueLoad error for file /usr/lib64/
/cvmfs/
Couldn't load XPCOM.
INFO: For gauge cancellation, the width of 't' has been set to zero.
INFO: For gauge cancellation, the width of 'X0' has been set to zero.
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 96 cores
INFO: Compiling P0_gu_x0tbxd...
INFO: Compiling P0_gdx_x0tbxux...
INFO: Compiling P0_ug_x0tbxd...
INFO: Compiling P0_dxg_x0tbxux...
When interrupting hanging tHq process:
WARNING: fct <function compile_dir at 0x7f3a947819d0> does not return 0. Stopping the code in a clean way. The error was:
A compilation Error occurs when trying to compile /.automount/
The compilation fails with the following output message:
gfortran -O -fno-automatic -ffixed-
check_
311 | 101 Z(K)=Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 101 at (1)
check_
313 | 102 Z(K)=(Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 102 at (1)
check_
325 | 105 XMT=XMT+ABS(XM(I))
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 105 at (1)
check_
365 | 202 Q(1,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 202 at (1)
check_
369 | 203 R(I)=0.
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 203 at (1)
check_
371 | DO 204 K=1,4
| 1
Warning: Fortran 2018 deleted feature: Shared DO termination label 204 at (1)
check_
372 | 204 R(K)=R(K)+Q(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 204 at (1)
check_
375 | 205 B(K)=-R(K)/RMAS
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 205 at (1)
check_
384 | 206 P(K,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 206 at (1)
check_
385 | 207 P(0,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 207 at (1)
check_
407 | 301 P2(I)=P(0,I)**2
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 301 at (1)
check_
417 | 303 G0=G0+P2(I)/E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 303 at (1)
check_
428 | 306 P(K,I)=X*P(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 306 at (1)
check_
429 | 307 P(0,I)=E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 307 at (1)
check_
436 | 308 WT3=WT3+
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 308 at (1)
run.inc:79:21:
79 | common/to_rwgt/ do_rwgt_scale, rw_Fscale_down, rw_Fscale_up, rw_Rscale_down, rw_Rscale_up,
| 1
Warning: Padding of 4 bytes required before 'rw_fscale_down' in COMMON 'to_rwgt' at (1); reorder elements or use '-fno-align-
run.inc:10:22:
10 | common/
| 1
Warning: Padding of 4 bytes required before 'mue_ref_fixed' in COMMON 'to_scale' at (1); reorder elements or use '-fno-align-
gfortran -O -fno-automatic -ffixed-
for dir in `ls -d V*`; do cd $dir; make; cd ../; done
make[1]: Entering directory '/.automount/
rm -f ../libMadLoop.a
gfortran -fno-automatic -ffixed-
make[1]: *** [makefile:34: polynomial.o] Interrupt
make: *** [makefile:130: libMadLoop.a] Interrupt
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
WARNING: fct <function compile_dir at 0x7f3a947819d0> does not return 0. Stopping the code in a clean way. The error was:
A compilation Error occurs when trying to compile /.automount/
The compilation fails with the following output message:
gfortran -O -fno-automatic -ffixed-
check_
311 | 101 Z(K)=Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 101 at (1)
check_
313 | 102 Z(K)=(Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 102 at (1)
check_
325 | 105 XMT=XMT+ABS(XM(I))
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 105 at (1)
check_
365 | 202 Q(1,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 202 at (1)
check_
369 | 203 R(I)=0.
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 203 at (1)
check_
371 | DO 204 K=1,4
| 1
Warning: Fortran 2018 deleted feature: Shared DO termination label 204 at (1)
check_
372 | 204 R(K)=R(K)+Q(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 204 at (1)
check_
375 | 205 B(K)=-R(K)/RMAS
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 205 at (1)
check_
384 | 206 P(K,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 206 at (1)
check_
385 | 207 P(0,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 207 at (1)
check_
407 | 301 P2(I)=P(0,I)**2
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 301 at (1)
check_
417 | 303 G0=G0+P2(I)/E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 303 at (1)
check_
428 | 306 P(K,I)=X*P(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 306 at (1)
check_
429 | 307 P(0,I)=E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 307 at (1)
check_
436 | 308 WT3=WT3+
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 308 at (1)
run.inc:79:21:
79 | common/to_rwgt/ do_rwgt_scale, rw_Fscale_down, rw_Fscale_up, rw_Rscale_down, rw_Rscale_up,
| 1
Warning: Padding of 4 bytes required before 'rw_fscale_down' in COMMON 'to_rwgt' at (1); reorder elements or use '-fno-align-
run.inc:10:22:
10 | common/
| 1
Warning: Padding of 4 bytes required before 'mue_ref_fixed' in COMMON 'to_scale' at (1); reorder elements or use '-fno-align-
gfortran -O -fno-automatic -ffixed-
for dir in `ls -d V*`; do cd $dir; make; cd ../; done
make[1]: Entering directory '/.automount/
rm -f ../libMadLoop.a
gfortran -fno-automatic -ffixed-
make[1]: *** [makefile:34: polynomial.o] Interrupt
make: *** [makefile:130: libMadLoop.a] Interrupt
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
WARNING: fct <function compile_dir at 0x7f3a947819d0> does not return 0. Stopping the code in a clean way. The error was:
A compilation Error occurs when trying to compile /.automount/
The compilation fails with the following output message:
gfortran -O -fno-automatic -ffixed-
check_
311 | 101 Z(K)=Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 101 at (1)
check_
313 | 102 Z(K)=(Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 102 at (1)
check_
325 | 105 XMT=XMT+ABS(XM(I))
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 105 at (1)
check_
365 | 202 Q(1,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 202 at (1)
check_
369 | 203 R(I)=0.
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 203 at (1)
check_
371 | DO 204 K=1,4
| 1
Warning: Fortran 2018 deleted feature: Shared DO termination label 204 at (1)
check_
372 | 204 R(K)=R(K)+Q(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 204 at (1)
check_
375 | 205 B(K)=-R(K)/RMAS
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 205 at (1)
check_
384 | 206 P(K,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 206 at (1)
check_
385 | 207 P(0,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 207 at (1)
check_
407 | 301 P2(I)=P(0,I)**2
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 301 at (1)
check_
417 | 303 G0=G0+P2(I)/E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 303 at (1)
check_
428 | 306 P(K,I)=X*P(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 306 at (1)
check_
429 | 307 P(0,I)=E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 307 at (1)
check_
436 | 308 WT3=WT3+
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 308 at (1)
run.inc:79:21:
79 | common/to_rwgt/ do_rwgt_scale, rw_Fscale_down, rw_Fscale_up, rw_Rscale_down, rw_Rscale_up,
| 1
Warning: Padding of 4 bytes required before 'rw_fscale_down' in COMMON 'to_rwgt' at (1); reorder elements or use '-fno-align-
run.inc:10:22:
10 | common/
| 1
Warning: Padding of 4 bytes required before 'mue_ref_fixed' in COMMON 'to_scale' at (1); reorder elements or use '-fno-align-
gfortran -O -fno-automatic -ffixed-
for dir in `ls -d V*`; do cd $dir; make; cd ../; done
make[1]: Entering directory '/.automount/
rm -f ../libMadLoop.a
gfortran -fno-automatic -ffixed-
make[1]: *** [makefile:34: polynomial.o] Interrupt
make: *** [makefile:130: libMadLoop.a] Interrupt
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
WARNING: fct <function compile_dir at 0x7f3a947819d0> does not return 0. Stopping the code in a clean way. The error was:
A compilation Error occurs when trying to compile /.automount/
The compilation fails with the following output message:
gfortran -O -fno-automatic -ffixed-
check_
311 | 101 Z(K)=Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 101 at (1)
check_
313 | 102 Z(K)=(Z(
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 102 at (1)
check_
325 | 105 XMT=XMT+ABS(XM(I))
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 105 at (1)
check_
365 | 202 Q(1,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 202 at (1)
check_
369 | 203 R(I)=0.
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 203 at (1)
check_
371 | DO 204 K=1,4
| 1
Warning: Fortran 2018 deleted feature: Shared DO termination label 204 at (1)
check_
372 | 204 R(K)=R(K)+Q(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 204 at (1)
check_
375 | 205 B(K)=-R(K)/RMAS
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 205 at (1)
check_
384 | 206 P(K,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 206 at (1)
check_
385 | 207 P(0,I)=
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 207 at (1)
check_
407 | 301 P2(I)=P(0,I)**2
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 301 at (1)
check_
417 | 303 G0=G0+P2(I)/E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 303 at (1)
check_
428 | 306 P(K,I)=X*P(K,I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 306 at (1)
check_
429 | 307 P(0,I)=E(I)
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 307 at (1)
check_
436 | 308 WT3=WT3+
| 1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 308 at (1)
run.inc:79:21:
79 | common/to_rwgt/ do_rwgt_scale, rw_Fscale_down, rw_Fscale_up, rw_Rscale_down, rw_Rscale_up,
| 1
Warning: Padding of 4 bytes required before 'rw_fscale_down' in COMMON 'to_rwgt' at (1); reorder elements or use '-fno-align-
run.inc:10:22:
10 | common/
| 1
Warning: Padding of 4 bytes required before 'mue_ref_fixed' in COMMON 'to_scale' at (1); reorder elements or use '-fno-align-
gfortran -O -fno-automatic -ffixed-
for dir in `ls -d V*`; do cd $dir; make; cd ../; done
make[1]: Entering directory '/.automount/
rm -f ../libMadLoop.a
gfortran -fno-automatic -ffixed-
make[1]: *** [makefile:34: polynomial.o] Interrupt
make: *** [makefile:130: libMadLoop.a] Interrupt
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
stopping all operation
in order to quit MadGraph5_aMC@NLO please enter exit
q
Cheers
Florian
|
Revision history for this message
|
#2 |
Hi,
the compilation of polynomial.f is cumbersome on some systems, in particular if the flag -O3 is given.
For this, I suggest you to edit SubProcesses/
POLYNOMIAL_
instead of O3.
Normally, it should work OK
Let us know
Marco
|
Revision history for this message
|
#3 |
Hi Marco,
thank you for your reply!
I tried it with your suggestion using the same syntax for tHq as detailed above.
I get such errors now:
`Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP No qp loop library can deal with this integral` with a crash afterwards.
Posting the full output and content of one of the log files below.
Do you have an idea how to solve this?
Thanks!
Cheers,
Florian
./mg5_aMC
*******
* *
* 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.3 2023-12-23 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* and *
* http://
* *
* 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/
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/
MG5_aMC> set fastjet /PATH/TO/
eMELA-config does not seem to correspond to a valid eMELA-config executable.
Please set the 'fastjet'variable to the full (absolute) /PATH/TO/
MG5_aMC> set eMELA /PATH/TO/
set ninja to /.automount/
set collier to /.automount/
set lhapdf to lhapdf-config
None does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/
Note that you can still compile and run aMC@NLO with the built-in PDFs
MG5_aMC> set lhapdf /PATH/TO/
Using default text editor "vi". Set another one in ./input/
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
Loading default model: sm
INFO: Restrict model sm with file ../models/
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+
MG5_aMC>launch tHq_4f_NLO_MH125_SM
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.3 2023-12-23 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /.automount/
INFO: load configuration from /.automount/
INFO: load configuration from /.automount/
Using default text editor "vi". Set another one in ./input/
launch auto
The following switches determine which programs are run:
/================== Description =======
| 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.[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.
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
\------
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path][90s to answer]
>
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
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
INFO: For gauge cancellation, the width of 't' has been set to zero.
INFO: For gauge cancellation, the width of 'X0' has been set to zero.
INFO: Compiling source...
/usr/sbin/
/usr/sbin/
XPCOMGlueLoad error for file /usr/lib64/
/cvmfs/
Couldn't load XPCOM.
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 96 cores
INFO: Compiling P0_gu_x0tbxd...
INFO: Compiling P0_gdx_x0tbxux...
INFO: Compiling P0_ug_x0tbxd...
INFO: Compiling P0_dxg_x0tbxux...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP No qp loop library can deal with this integral
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP No qp loop library can deal with this integral
INFO: P0_dxg_x0tbxux done.
INFO: P0_gdx_x0tbxux done.
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP No qp loop library can deal with this integral
INFO: P0_gu_x0tbxd done.
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP No qp loop library can deal with this integral
INFO: P0_ug_x0tbxd done.
INFO: Checking test output:
INFO: P0_gu_x0tbxd
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_gdx_x0tbxux
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_ug_x0tbxd
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_dxg_x0tbxux
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: Starting run
INFO: Using 96 cores
INFO: Cleaning previous results
INFO: Generating events without running the shower.
INFO: Setting up grids
INFO: Idle: 0, Running: 0, Completed: 32 [ current time: 12h10 ]
Error detected in "launch auto "
write debug file /.automount/
If you need help with this issue, please, contact us on https:/
str : An error occurred during the collection of results.
Please check the .log files inside the directories which failed:
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
/.automount/
quit
the first of the mentioned log files contains this:
304549
======
INFO: MadFKS read these parameters from FKS_params.dat
======
> IRPoleCheckThre
> PrecisionVirtua
> SelectedContrib
> VetoedContribut
> QCD_squared_
> QED_squared_
> SelectedCouplin
> NHelForMCoverHels = 4
> VirtualFraction = 1.0000000000000000
> MinVirtualFraction = 5.0000000000000
> SeparateFlavour
> UsePolyVirtual = F
======
SPLIT TYPE USED: T F T
A PDF is used, so alpha_s(MZ) is going to be modified
Old value of alpha_s from param_card: 0.11900000000000002
*****
NNPDFDriver version 1.0.3
Grid: NNPDF23nlo_
*****
New value of alpha_s from PDF nn23nlo: 0.11899999999999999
******
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
******
External Params
-----
MU_R = 91.188000000000002
aEWM1 = 132.50700000000001
mdl_Gf = 1.1663900000000
aS = 0.11900000000000002
mdl_ymb = 4.7000000000000002
mdl_ymt = 173.00000000000000
mdl_ymtau = 1.7769999999999999
mdl_MT = 173.00000000000000
mdl_MB = 4.7000000000000002
mdl_MZ = 91.188000000000002
mdl_MTA = 1.7769999999999999
mdl_WT = 0.0000000000000000
mdl_WZ = 2.4414039999999999
mdl_WW = 2.0476000000000001
mdl_MX0 = 125.00000000000000
mdl_WX0 = 0.0000000000000000
mdl_Lambda = 1000.0000000000000
mdl_cosa = 1.0000000000000000
mdl_kSM = 1.0000000000000000
mdl_kHtt = 1.0000000000000000
mdl_kAtt = 1.0000000000000000
mdl_kHbb = 1.0000000000000000
mdl_kAbb = 1.0000000000000000
mdl_kHll = 1.0000000000000000
mdl_kAll = 1.0000000000000000
mdl_kHaa = 1.0000000000000000
mdl_kAaa = 1.0000000000000000
mdl_kHza = 1.0000000000000000
mdl_kAza = 1.0000000000000000
mdl_kHzz = 0.0000000000000000
mdl_kAzz = 0.0000000000000000
mdl_kHww = 0.0000000000000000
mdl_kAww = 0.0000000000000000
mdl_kHda = 0.0000000000000000
mdl_kHdz = 0.0000000000000000
mdl_kHdwR = 0.0000000000000000
mdl_kHdwI = 0.0000000000000000
Internal Params
-----
mdl_CKM22 = 1.0000000000000000
mdl_lhv = 1.0000000000000000
mdl_conjg__CKM22 = 1.0000000000000000
mdl_CKM3x3 = 1.0000000000000000
mdl_conjg__CKM3x3 = 1.0000000000000000
mdl_CKM33 = 1.0000000000000000
mdl_Ncol = 3.0000000000000000
mdl_CA = 3.0000000000000000
mdl_TF = 0.50000000000000000
mdl_CF = 1.3333333333333333
mdl_complexi = (0.000000000000
mdl_MZ__exp__2 = 8315.2513440000002
mdl_sqrt__2 = 1.4142135623730951
mdl_MX0__exp__2 = 15625.000000000000
mdl_cosa__exp__2 = 1.0000000000000000
mdl_sina = 0.0000000000000000
mdl_kHdw = (0.000000000000
mdl_nb_
mdl_Ncol__exp__2 = 9.0000000000000000
mdl_MB__exp__2 = 22.090000000000003
mdl_MT__exp__2 = 29929.000000000000
mdl_conjg__kHdw = (0.000000000000
mdl_aEW = 7.5467711139788
mdl_MW = 80.419002445756163
mdl_sqrt__aEW = 8.6872153846781
mdl_ee = 0.30795376724436879
mdl_MW__exp__2 = 6467.2159543705357
mdl_sw2 = 0.22224648578577766
mdl_cw = 0.88190334743339216
mdl_sqrt__sw2 = 0.47143025548407230
mdl_sw = 0.47143025548407230
mdl_g1 = 0.34919219678733299
mdl_gw = 0.65323293034757990
mdl_vev = 246.21845810181637
mdl_vev__exp__2 = 60623.529110035903
mdl_lam = 0.12886910601690263
mdl_yb = 2.6995554250465
mdl_yt = 0.99366614581500623
mdl_ytau = 1.0206617000654
mdl_muH = 88.388347648318430
mdl_AxialZUp = -0.185177018617
mdl_AxialZDown = 0.18517701861793787
mdl_VectorZUp = 7.5430507588273
mdl_VectorZDown = -0.130303763103
mdl_VectorAUp = 0.20530251149624587
mdl_VectorADown = -0.102651255748
mdl_VectorWmDxU = 0.23095271737156670
mdl_AxialWmDxU = -0.230952717371
mdl_VectorWpUxD = 0.23095271737156670
mdl_AxialWpUxD = -0.230952717371
mdl_I1x33 = (2.699555425046
mdl_I2x33 = (0.993666145815
mdl_I3x33 = (0.993666145815
mdl_I4x33 = (2.699555425046
mdl_Vector_tbGp = (-0.96667059156
mdl_Axial_tbGp = (-1.02066170006
mdl_Vector_tbGm = (0.966670591564
mdl_Axial_tbGm = (-1.02066170006
mdl_ee__exp__2 = 9.4835522759998
mdl_gAaa = 1.3008566310666
mdl_cw__exp__2 = 0.77775351421422245
mdl_gAza = 4.7794971072590
mdl_gHaa = 2.5475109025056
mdl_gHza = 3.9182129211851
mdl_gw__exp__2 = 0.42671326129048615
mdl_sw__exp__2 = 0.22224648578577769
Internal Params evaluated point by point
-----
mdl_sqrt__aS = 0.34496376621320685
mdl_G__exp__2 = 1.4953981031087418
mdl_G__exp__4 = 2.2362154867812229
mdl_R2MixedFac
mdl_G_UVg_1EPS_ = -5.208345512682
mdl_G_UVb_1EPS_ = 3.1565730379892
mdl_bMass_UV_1EPS_ = (0.000000000000
mdl_tMass_UV_1EPS_ = (0.000000000000
mdl_UV_yuk_b_1EPS_ = -3.787887645587
mdl_UV_yuk_t_1EPS_ = -3.787887645587
mdl_G__exp__3 = 1.8286695702272273
mdl_MU_R__exp__2 = 8315.2513440000002
mdl_gAgg = 7.6921277039691
mdl_gHgg = -5.128085135979
mdl_G_UVb_FIN_ = 1.8720755903176
mdl_G_UVt_FIN_ = -4.042738525768
mdl_bMass_UV_FIN_ = (0.000000000000
mdl_tMass_UV_FIN_ = (0.000000000000
mdl_UV_yuk_b_FIN_ = -0.275154239445
mdl_UV_yuk_t_FIN_ = -1.992306298605
Couplings of HC_NLO_X0_UFO
-----
GC_4 = -0.12229E+01 0.00000E+00
GC_5 = 0.00000E+00 0.12229E+01
R2_GQQ = -0.00000E+00 -0.30880E-01
R2_QQq = 0.00000E+00 0.12626E-01
R2_QQb = 0.00000E+00 0.11869E+00
R2_QQt = 0.00000E+00 0.43687E+01
UV_GQQg_1eps = 0.00000E+00 -0.63691E-01
UV_GQQb_1eps = 0.00000E+00 0.38601E-02
UV_bMass_1eps = 0.00000E+00 0.17803E+00
UV_tMass_1eps = 0.00000E+00 0.65530E+01
UVWfct_b_0_1eps -0.18939E-01 0.00000E+00
UVWfct_G_1_1eps -0.31566E-02 0.00000E+00
R2_bxtW = -0.00000E+00 -0.11664E-01
R2_bbX0_h = 0.00000E+00 0.96408E-03
R2_ttX0_h = 0.00000E+00 0.35486E-01
R2_bbX0_a = -0.00000E+00 0.00000E+00
R2_ttX0_a = -0.00000E+00 0.00000E+00
R2_ggX0_hb = 0.00000E+00 -0.16992E-02
R2_ggX0_ht = 0.00000E+00 -0.23022E+01
UV_X0bb_h_1eps 0.00000E+00 0.72306E-03
UV_X0tt_h_1eps 0.00000E+00 0.26615E-01
UV_X0bb_a_1eps -0.00000E+00 0.00000E+00
UV_X0tt_a_1eps -0.00000E+00 0.00000E+00
UV_GQQb = 0.00000E+00 0.22893E-01
UV_GQQt = 0.00000E+00 -0.49437E-02
UV_bMass = 0.00000E+00 0.12932E+01
UV_tMass = 0.00000E+00 0.34467E+00
UVWfct_b_0 = -0.13758E+00 -0.00000E+00
UVWfct_t_0 = -0.99615E-03 -0.00000E+00
UVWfct_G_1 = -0.18721E-01 0.00000E+00
UVWfct_G_2 = 0.40427E-02 0.00000E+00
UV_X0bb_h = 0.00000E+00 0.52523E-02
UV_X0tt_h = 0.00000E+00 0.13999E-02
UV_X0bb_a = -0.00000E+00 0.00000E+00
UV_X0tt_a = -0.00000E+00 0.00000E+00
GC_1 = -0.00000E+00 -0.10265E+00
GC_2 = 0.00000E+00 0.20530E+00
GC_7 = 0.00000E+00 0.57609E+00
GC_8 = -0.00000E+00 -0.42671E+00
GC_9 = 0.00000E+00 0.33188E+00
GC_11 = 0.00000E+00 0.46191E+00
GC_21 = -0.00000E+00 -0.28804E+00
GC_23 = -0.00000E+00 -0.27437E-01
GC_25 = 0.00000E+00 0.30795E+00
GC_26 = -0.00000E+00 -0.35482E+00
GC_27 = 0.00000E+00 0.94836E-01
GC_3001h = -0.00000E+00 -0.25475E-04
GC_3001a = -0.00000E+00 -0.00000E+00
GC_3005h1 = 0.00000E+00 0.52532E+02
GC_3006a = -0.00000E+00 -0.00000E+00
GC_3006h2 = -0.00000E+00 -0.39182E-04
GC_3007h1 = 0.00000E+00 0.67544E+02
GC_3009a = 0.00000E+00 0.00000E+00
GC_3009h = -0.00000E+00 -0.19089E-01
GC_3010a = 0.00000E+00 0.00000E+00
GC_3010h = -0.00000E+00 -0.70263E+00
GC_3005a = -0.00000E+00 -0.00000E+00
GC_3005h2 = -0.00000E+00 -0.00000E+00
GC_3005h3 = -0.00000E+00 -0.00000E+00
GC_3005h4 = -0.00000E+00 -0.00000E+00
GC_3005h4 = -0.00000E+00 -0.00000E+00
GC_3006h1 = -0.00000E+00 -0.00000E+00
GC_3007a = -0.00000E+00 -0.00000E+00
GC_3007h2 = -0.00000E+00 -0.00000E+00
GC_3007h3 = -0.00000E+00 -0.00000E+00
Collider parameters:
------
Running at P P machine @ 13000.000000000000 GeV
PDF set = nn23nlo
alpha_s(Mz)= 0.1190 running at 2 loops.
alpha_s(Mz)= 0.1190 running at 2 loops.
Renormalization scale set on event-by-event basis
Factorization scale set on event-by-event basis
INFO: orders_tag_plot is computed as: + QCD * 1 + QNP * 100 + QED * 10000
orders_tag_plot= 60002 for QCD,QNP,QED, = 2 , 0 , 6 ,
AMP_SPLIT: 1 correspond to S.O. 2 0 6
orders_tag_plot= 50102 for QCD,QNP,QED, = 2 , 1 , 5 ,
AMP_SPLIT: 2 correspond to S.O. 2 1 5
orders_tag_plot= 40202 for QCD,QNP,QED, = 2 , 2 , 4 ,
AMP_SPLIT: 3 correspond to S.O. 2 2 4
orders_tag_plot= 60004 for QCD,QNP,QED, = 4 , 0 , 6 ,
AMP_SPLIT: 4 correspond to S.O. 4 0 6
orders_tag_plot= 50104 for QCD,QNP,QED, = 4 , 1 , 5 ,
AMP_SPLIT: 5 correspond to S.O. 4 1 5
orders_tag_plot= 40204 for QCD,QNP,QED, = 4 , 2 , 4 ,
AMP_SPLIT: 6 correspond to S.O. 4 2 4
orders_tag_plot= 60202 for QCD,QNP,QED, = 2 , 2 , 6 ,
AMP_SPLIT: 7 correspond to S.O. 2 2 6
orders_tag_plot= 70102 for QCD,QNP,QED, = 2 , 1 , 7 ,
AMP_SPLIT: 8 correspond to S.O. 2 1 7
orders_tag_plot= 80200 for QCD,QNP,QED, = 0 , 2 , 8 ,
AMP_SPLIT: 9 correspond to S.O. 0 2 8
orders_tag_plot= 50302 for QCD,QNP,QED, = 2 , 3 , 5 ,
AMP_SPLIT: 10 correspond to S.O. 2 3 5
getting user params
Enter number of events and iterations:
Number of events and iterations -1 12
Enter desired fractional accuracy:
Desired fractional accuracy: 2.9999999999999
Enter alpha, beta for G_soft
Enter alpha<0 to set G_soft=1 (no ME soft)
for G_soft: alpha= 1.0000000000000000 , beta= -0.100000000000
Enter alpha, beta for G_azi
Enter alpha>0 to set G_azi=0 (no azi corr)
for G_azi: alpha= -1.0000000000000000 , beta= -0.100000000000
Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 1 0
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
MINT running mode: 0
Set the three folding parameters for MINT
xi_i, y_ij, phi_i
1 1 1
'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
Enter 'born0' or 'virt0' to perform
a pure n-body integration (no S functions)
doing the all of this channel
Normal integration (Sfunction != 1)
about to integrate 13 -1 12 1
imode is 0
channel 1 : 1 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
------- iteration 1
Update # PS points (even_rn): 1040 --> 1040
Using random seed offsets: 1 , 1 , 0
with seed 33
Ranmar initialization seeds 11949 9408
Total number of FKS directories is 17
Initial state 2 g->gg already found in driver_mintMC
21 21 21 21 21 22 22 22 22 0 0 0 0 0 0 0 0
21 2 6 -5 1 2 6 -5 1 0 0 0 0 0 0 0 0
1 2 4 5 6 2 4 5 6 0 0 0 0 0 0 0 0
Time in seconds: 0
|
Revision history for this message
|
#4 |
Hi,
can you try to re-generate the process, issuing
set nlo_mixed_expansion False
bofore the generate command ?
Thanks
Marco
|
Revision history for this message
|
#5 |
Hi,
thanks for your quick reply!
Setting nlo_mixed_expansion False makes it work for the SM configuration of tHq.
Trying to constrain coupling orders such that it will only simulate the Htt-HWW interference, as formulated in the originial request in https:/
The syntax I used trying to constrain it to the Htt-HWW interference with the original HC model is this:
```
import model HC_NLO_X0_UFO
set nlo_mixed_expansion False
generate p p > t x0 b~ j $$ w+ w- QNP^2<=1 QED^2<=5 QCD^2<=2 [QCD]
add process p p > t~ x0 b j $$ w+ w- QNP^2<=1 QED^2<=5 QCD^2<=2 [QCD]
```
It fails with such messages:
WARNING: program /.automount/
Do you know how this can be handled?
I am attaching the output and the log file content below.
Thanks and best wishes,
Florian
Full output:
(base) ./mg5_aMC
*******
* *
* 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.3 2023-12-23 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* and *
* http://
* *
* 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/
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/
MG5_aMC> set fastjet /PATH/TO/
eMELA-config does not seem to correspond to a valid eMELA-config executable.
Please set the 'fastjet'variable to the full (absolute) /PATH/TO/
MG5_aMC> set eMELA /PATH/TO/
set ninja to /.automount/
set collier to /.automount/
set lhapdf to lhapdf-config
None does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/
Note that you can still compile and run aMC@NLO with the built-in PDFs
MG5_aMC> set lhapdf /PATH/TO/
Using default text editor "vi". Set another one in ./input/
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
Loading default model: sm
INFO: Restrict model sm with file ../models/
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+
MG5_aMC>launch tHq_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.3 2023-12-23 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /.automount/
INFO: load configuration from /.automount/
INFO: load configuration from /.automount/
Using default text editor "vi". Set another one in ./input/
launch auto
The following switches determine which programs are run:
/================== Description =======
| 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|
| 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 = Not Avail | Please install module |
\======
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]
>3
The following switches determine which programs are run:
/================== Description =======
| 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 = OFF | PYTHIA6Q|
| 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 = Not Avail | Please install module |
\======
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: noshower
WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical. Please, shower the LesHouches events before using them for physics analyses. You have to choose NOW which parton-shower you WILL use and specify it in the run_card.
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
\------
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path][90s to answer]
>2
open run_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
\------
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path]
>
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
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
/usr/sbin/
/usr/sbin/
INFO: For gauge cancellation, the width of 't' has been set to zero.
XPCOMGlueLoad error for file /usr/lib64/
/cvmfs/
Couldn't load XPCOM.
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 96 cores
INFO: Compiling P0_gu_x0tbxd...
INFO: Compiling P0_gdx_x0tbxux...
INFO: Compiling P0_ug_x0tbxd...
INFO: Compiling P0_dxg_x0tbxux...
INFO: Compiling P0_gd_x0btxu...
INFO: Compiling P0_gux_x0btxdx...
INFO: Compiling P0_dg_x0btxu...
INFO: Compiling P0_uxg_x0btxdx...
INFO: P0_dg_x0btxu done.
INFO: P0_uxg_x0btxdx done.
INFO: P0_gd_x0btxu done.
INFO: P0_gu_x0tbxd done.
INFO: P0_dxg_x0tbxux done.
INFO: P0_gdx_x0tbxux done.
INFO: P0_gux_x0btxdx done.
INFO: P0_ug_x0tbxd done.
INFO: Checking test output:
INFO: P0_gu_x0tbxd
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_gdx_x0tbxux
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_ug_x0tbxd
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_dxg_x0tbxux
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_gd_x0btxu
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_gux_x0btxdx
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_dg_x0btxu
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_uxg_x0btxdx
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 96 cores
INFO: Cleaning previous results
INFO: Generating events without running the shower.
INFO: Setting up grids
WARNING: program /.automount/
WARNING: program /.automount/
WARNING: program /.automount/
Exception in thread Thread-53:
Traceback (most recent call last):
File "/cvmfs/
self.run()
File "/cvmfs/
self.
File "/.automount/
WARNING: program /.automount/
self.
ValueError: list.remove(x): x not in list
Terminated
WARNING: program /.automount/
WARNING: program /.automount/
WARNING: program /.automount/
WARNING: program /.automount/
Exception in thread Thread-60:
Traceback (most recent call last):
File "/cvmfs/
self.run()
File "/cvmfs/
self.
File "/.automount/
self.
ValueError: list.remove(x): x not in list
WARNING: program /.automount/
Terminated
Terminated
WARNING: program /.automount/
INFO: Idle: 0, Running: 1, Completed: 63 [ current time: 17h10 ]
INFO: Idle: 0, Running: 0, Completed: 64 [ 0.38s ]
Command "launch auto " interrupted with error:
Exception : program /.automount/
Please report this bug on https:/
More information is found in '/.automount/
Please attach this file to your report.
INFO:
quit
The content of the log file is this:
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.3 2023-12-23 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* and *
#* http://
#* *
#******
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#******
launch auto
Traceback (most recent call last):
File "/.automount/
return self.onecmd_
File "/.automount/
return func(arg, **opt)
File "/.automount/
evt_file = self.run(mode, options)
File "/.automount/
self.
File "/.automount/
self.
File "/.automount/
self.
File "/.automount/
raise Exception(
Exception: program /.automount/
Value of current Options:
exrootana
cluster_
mg5amc_
output_
acknowledged_
ignore_
low_mem_
complex_
include_
loop_
max_
default_
nlo_
automatic_
notificat
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.3 2023-12-23 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set low_mem_
set complex_mass_scheme False
set include_
set gauge unitary
set loop_optimized_
set loop_color_flows False
set max_npoint_
set default_
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
set auto_convert_model T
import model HC_NLO_X0_UFO
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 nlo_mixed_expansion False
generate p p > t x0 b~ j $$ w+ w- QNP^2<=1 QED^2<=5 QCD^2<=2 [QCD]
add process p p > t~ x0 b j $$ w+ w- QNP^2<=1 QED^2<=5 QCD^2<=2 [QCD]
output tHq_4f_
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## INFORMATION FOR FRBLOCK
#######
BLOCK FRBLOCK #
1 1.000000e+03 # lambda
2 1.000000e+00 # cosa
3 1.000000e+00 # ksm
4 1.000000e+00 # khtt
5 1.000000e+00 # katt
6 1.000000e+00 # khbb
7 1.000000e+00 # kabb
8 1.000000e+00 # khll
9 1.000000e+00 # kall
10 1.000000e+00 # khaa
11 1.000000e+00 # kaaa
12 1.000000e+00 # khza
13 1.000000e+00 # kaza
16 0.000000e+00 # khzz
17 0.000000e+00 # kazz
18 0.000000e+00 # khww
19 0.000000e+00 # kaww
20 0.000000e+00 # khda
21 0.000000e+00 # khdz
22 0.000000e+00 # khdwr (real part of khdw)
23 0.000000e+00 # khdwi (imaginary part of khdw)
#######
## INFORMATION FOR LOOP
#######
BLOCK LOOP #
1 9.118800e+01 # mu_r
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
1 1.325070e+02 # aewm1
2 1.166390e-05 # gf
3 1.190000e-01 # as
#######
## INFORMATION FOR MASS
#######
BLOCK MASS #
5 4.700000e+00 # mb
6 1.730000e+02 # mt
15 1.777000e+00 # mta
23 9.118800e+01 # mz
25 1.250000e+02 # mx0
1 0.000000e+00 # d : 0.0
2 0.000000e+00 # u : 0.0
3 0.000000e+00 # s : 0.0
4 0.000000e+00 # c : 0.0
11 0.000000e+00 # e- : 0.0
12 0.000000e+00 # ve : 0.0
13 0.000000e+00 # mu- : 0.0
14 0.000000e+00 # vm : 0.0
16 0.000000e+00 # vt : 0.0
21 0.000000e+00 # g : 0.0
22 0.000000e+00 # a : 0.0
24 8.041900e+01 # w+ : cmath.sqrt(
82 0.000000e+00 # gh : 0.0
#######
## INFORMATION FOR YUKAWA
#######
BLOCK YUKAWA #
5 4.700000e+00 # ymb
6 1.730000e+02 # ymt
15 1.777000e+00 # ymtau
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.491500e+00 # wt
DECAY 23 2.441404e+00 # wz
DECAY 24 2.047600e+00 # ww
DECAY 25 4.070000e-03 # wx0
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 82 0.000000e+00 # gh : 0.0
#######
## INFORMATION FOR QNUMBERS 82
#######
BLOCK QNUMBERS 82 # gh
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# Some of the values of variables can be list. These can either be *
# comma or space separated. *
# *
# To display additional parameter, you can use the command: *
# update to_full *
#******
#
#******
# Running parameters
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#******
1000 = nevents ! Number of unweighted events requested
-1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
-1 = nevt_job! Max number of events per job in event generation.
! (-1= no split).
#******
# Output format
#******
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
average = event_norm ! average/sum/bias. Normalization of the weight in the LHEF
#******
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#******
0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
# These numbers are ignored except if req_acc_FO is equal to -1
5000 = npoints_FO_grid ! number of points to setup grids
4 = niters_FO_grid ! number of iter. to setup grids
10000 = npoints_FO ! number of points to compute Xsec
6 = niters_FO ! number of iter. to compute Xsec
#******
# Random number seed *
#******
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# 0 = no PDF *
# 1/-1 = proton/antiproton *
# 3/-3 = electron/positron with ISR/Beamstrahlung; *
# 4/-4 = muon/antimuon with ISR/Beamstrahlung; *
#******
1 = lpp1 ! beam 1 type (0 = no PDF)
1 = lpp2 ! beam 2 type (0 = no PDF)
6500.0 = ebeam1 ! beam 1 energy in GeV
6500.0 = ebeam2 ! beam 2 energy in GeV
#******
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#******
nn23nlo = pdlabel ! PDF set
244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
! numbers for central PDF sets are allowed. Can be a list;
! PDF sets beyond the first are included via reweighting.
0 = pdfscheme ! the scheme of the input PDFs. 0->MSbar; 1->DIS
! 2->eta (leptonic); 3->beta (leptonic)
! 4->mixed (leptonic); 6->delta (leptonic)
#******
# The following block is specific to lepton collisions (lpp=+-3) *
#******
True = photons_from_lepton ! whether to include or not photons from
#******
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#******
HERWIG6 = parton_shower
1.0 = shower_scale_factor ! multiply default shower starting
False = MCatNLO_DELTA ! use MC@NLO-Delta matching, arXiv:2002.12716
#******
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses divided by two of all final state particles *
# and partons. This can be changed in SubProcesses/
# dynamical_
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.118 = muR_ref_fixed ! fixed ren reference scale
91.118 = muF_ref_fixed ! fixed fact reference scale
-1 = dynamical_
! dynamical choices. Can be a list; scale choices beyond the
! first are included via reweighting
1.0 = muR_over_ref ! ratio of current muR over reference muR
1.0 = muF_over_ref ! ratio of current muF over reference muF
#******
# Reweight variables for scale dependence and PDF uncertainty *
#******
1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
True = reweight_scale ! Reweight to get scale variation using the
! rw_rscale and rw_fscale factors. Should be a list of
! booleans of equal length to dynamical_
! specify for which choice to include scale dependence.
False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
! list booleans of equal length to lhaid to specify for
! which PDF set to include the uncertainties.
#******
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#******
False = store_rwgt_info ! Store info for reweighting in LHE file
#******
# Customization of the code. List of files containing user hook function
#******
= custom_fcts ! List of files containing user hook function
#******
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#******
0 = ickkw
#******
#
#******
# BW cutoff (M+/-bwcutoff*
# written in the LHE event file *
#******
15.0 = bwcutoff
#******
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - If gamma_is_j, photons are also clustered with jets. *
# Otherwise, they will be treated as tagged particles and photon *
# isolation will be applied. Note that photons in the real emission *
# will always be clustered with QCD partons. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#******
1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
0.7 = jetradius ! The radius parameter for the jet algorithm
10.0 = ptj ! Min jet transverse momentum
-1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
False = gamma_is_j! Wether to cluster photons as jets or not
#******
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#******
0.0 = ptl ! Min lepton transverse momentum
-1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
0.0 = drll ! Min distance between opposite sign lepton pairs
0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#******
# Fermion-photon recombination parameters *
# If Rphreco=0, no recombination is performed *
#******
0.1 = Rphreco ! Minimum fermion-photon distance for recombination
-1.0 = etaphreco ! Maximum abs(pseudo-rap) for photons to be recombined (a value .lt.0 means no cut)
False = lepphreco ! Recombine photons and leptons together
False = quarkphreco ! Recombine photons and quarks together
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# Not applied if gamma_is_j *
# When ptgmin=0, all the other parameters are ignored *
# More specific cuts can be specified in SubProcesses/cuts.f *
#******
20.0 = ptgmin ! Min photon transverse momentum
-1.0 = etagamma ! Max photon abs(pseudo-rap)
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# Cuts associated to MASSIVE particles identified by their PDG codes. *
# All cuts are applied to both particles and anti-particles, so use *
# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
# {6:100, 25:200} for multiple particles *
#******
{} = pt_min_pdg ! Min pT for a massive particle
{} = pt_max_pdg ! Max pT for a massive particle
{} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
#******
# Use PineAPPL to generate PDF-independent fast-interpolation grid *
# (https:/
#******
False = pineappl ! PineAPPL switch
#******
# Folding parameters for S-events to reduce the number of negatively *
# weighted events. Allowed values are 1, 2, 4 or 8 for each of the *
# three variables. Typically, folding in xi_i or y_ij results in the *
# largest reduction of negatively weighted events. (arXiv:2002.12716) *
#******
1, 1, 1 = folding ! correspond to folding in xi_i, y_ij, and phi_i
#******
|
Revision history for this message
|
#6 |
Would you have a directory
/.automount/
and if yes, would you have a log file in that directory?
If yes, could you copy/paste his content here?
Thanks,
Olivier
|
Revision history for this message
|
#7 |
Hi,
yes, the log file exists, I paste the content below.
Thanks!
Cheers,
Florian
1266736
======
INFO: MadFKS read these parameters from FKS_params.dat
======
> IRPoleCheckThre
> PrecisionVirtua
> SelectedContrib
> VetoedContribut
> QCD_squared_
> QED_squared_
> SelectedCouplin
> NHelForMCoverHels = 4
> VirtualFraction = 1.0000000000000000
> MinVirtualFraction = 5.0000000000000
> SeparateFlavour
> UsePolyVirtual = F
======
SPLIT TYPE USED: T F F
A PDF is used, so alpha_s(MZ) is going to be modified
Old value of alpha_s from param_card: 0.11900000000000002
*****
NNPDFDriver version 1.0.3
Grid: NNPDF23nlo_
*****
New value of alpha_s from PDF nn23nlo: 0.11899999999999999
******
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
******
External Params
-----
MU_R = 91.188000000000002
aEWM1 = 132.50700000000001
mdl_Gf = 1.1663900000000
aS = 0.11900000000000002
mdl_ymb = 4.7000000000000002
mdl_ymt = 173.00000000000000
mdl_ymtau = 1.7769999999999999
mdl_MT = 173.00000000000000
mdl_MB = 4.7000000000000002
mdl_MZ = 91.188000000000002
mdl_MTA = 1.7769999999999999
mdl_WT = 0.0000000000000000
mdl_WZ = 2.4414039999999999
mdl_WW = 2.0476000000000001
mdl_MX0 = 125.00000000000000
mdl_WX0 = 4.0699999999999
mdl_Lambda = 1000.0000000000000
mdl_cosa = 1.0000000000000000
mdl_kSM = 1.0000000000000000
mdl_kHtt = 1.0000000000000000
mdl_kAtt = 1.0000000000000000
mdl_kHbb = 1.0000000000000000
mdl_kAbb = 1.0000000000000000
mdl_kHll = 1.0000000000000000
mdl_kAll = 1.0000000000000000
mdl_kHaa = 1.0000000000000000
mdl_kAaa = 1.0000000000000000
mdl_kHza = 1.0000000000000000
mdl_kAza = 1.0000000000000000
mdl_kHzz = 0.0000000000000000
mdl_kAzz = 0.0000000000000000
mdl_kHww = 0.0000000000000000
mdl_kAww = 0.0000000000000000
mdl_kHda = 0.0000000000000000
mdl_kHdz = 0.0000000000000000
mdl_kHdwR = 0.0000000000000000
mdl_kHdwI = 0.0000000000000000
Internal Params
-----
mdl_lhv = 1.0000000000000000
mdl_CKM22 = 1.0000000000000000
mdl_conjg__CKM22 = 1.0000000000000000
mdl_conjg__CKM3x3 = 1.0000000000000000
mdl_CKM33 = 1.0000000000000000
mdl_CKM3x3 = 1.0000000000000000
mdl_Ncol = 3.0000000000000000
mdl_CA = 3.0000000000000000
mdl_TF = 0.50000000000000000
mdl_CF = 1.3333333333333333
mdl_complexi = (0.000000000000
mdl_MZ__exp__2 = 8315.2513440000002
mdl_MZ__exp__4 = 69143404.913893804
mdl_sqrt__2 = 1.4142135623730951
mdl_MX0__exp__2 = 15625.000000000000
mdl_cosa__exp__2 = 1.0000000000000000
mdl_sina = 0.0000000000000000
mdl_kHdw = (0.000000000000
mdl_nb_
mdl_Ncol__exp__2 = 9.0000000000000000
mdl_MB__exp__2 = 22.090000000000003
mdl_MT__exp__2 = 29929.000000000000
mdl_conjg__kHdw = (0.000000000000
mdl_aEW = 7.5467711139788
mdl_MW = 80.419002445756163
mdl_sqrt__aEW = 8.6872153846781
mdl_ee = 0.30795376724436879
mdl_MW__exp__2 = 6467.2159543705357
mdl_sw2 = 0.22224648578577766
mdl_cw = 0.88190334743339216
mdl_sqrt__sw2 = 0.47143025548407230
mdl_sw = 0.47143025548407230
mdl_g1 = 0.34919219678733299
mdl_gw = 0.65323293034757990
mdl_vev = 246.21845810181637
mdl_vev__exp__2 = 60623.529110035903
mdl_lam = 0.12886910601690263
mdl_yb = 2.6995554250465
mdl_yt = 0.99366614581500623
mdl_ytau = 1.0206617000654
mdl_muH = 88.388347648318430
mdl_AxialZUp = -0.185177018617
mdl_AxialZDown = 0.18517701861793787
mdl_VectorZUp = 7.5430507588273
mdl_VectorZDown = -0.130303763103
mdl_VectorAUp = 0.20530251149624587
mdl_VectorADown = -0.102651255748
mdl_VectorWmDxU = 0.23095271737156670
mdl_AxialWmDxU = -0.230952717371
mdl_VectorWpUxD = 0.23095271737156670
mdl_AxialWpUxD = -0.230952717371
mdl_I1x33 = (2.699555425046
mdl_I2x33 = (0.993666145815
mdl_I3x33 = (0.993666145815
mdl_I4x33 = (2.699555425046
mdl_Vector_tbGp = (-0.96667059156
mdl_Axial_tbGp = (-1.02066170006
mdl_Vector_tbGm = (0.966670591564
mdl_Axial_tbGm = (-1.02066170006
mdl_ee__exp__2 = 9.4835522759998
mdl_gAaa = 1.3008566310666
mdl_cw__exp__2 = 0.77775351421422245
mdl_gAza = 4.7794971072590
mdl_gHaa = 2.5475109025056
mdl_gHza = 3.9182129211851
mdl_gw__exp__2 = 0.42671326129048615
mdl_sw__exp__2 = 0.22224648578577769
Internal Params evaluated point by point
-----
mdl_sqrt__aS = 0.34496376621320685
mdl_G__exp__2 = 1.4953981031087418
mdl_G__exp__4 = 2.2362154867812229
mdl_R2MixedFac
mdl_G_UVg_1EPS_ = -5.208345512682
mdl_G_UVb_1EPS_ = 3.1565730379892
mdl_bMass_UV_1EPS_ = (0.000000000000
mdl_tMass_UV_1EPS_ = (0.000000000000
mdl_UV_yuk_b_1EPS_ = -3.787887645587
mdl_UV_yuk_t_1EPS_ = -3.787887645587
mdl_G__exp__3 = 1.8286695702272273
mdl_MU_R__exp__2 = 8315.2513440000002
mdl_gAgg = 7.6921277039691
mdl_gHgg = -5.128085135979
mdl_G_UVb_FIN_ = 1.8720755903176
mdl_G_UVt_FIN_ = -4.042738525768
mdl_bMass_UV_FIN_ = (0.000000000000
mdl_tMass_UV_FIN_ = (0.000000000000
mdl_UV_yuk_b_FIN_ = -0.275154239445
mdl_UV_yuk_t_FIN_ = -1.992306298605
Couplings of HC_NLO_X0_UFO
-----
GC_4 = -0.12229E+01 0.00000E+00
GC_5 = 0.00000E+00 0.12229E+01
R2_GQQ = -0.00000E+00 -0.30880E-01
R2_QQq = 0.00000E+00 0.12626E-01
R2_QQb = 0.00000E+00 0.11869E+00
R2_QQt = 0.00000E+00 0.43687E+01
UV_GQQg_1eps = 0.00000E+00 -0.63691E-01
UV_GQQb_1eps = 0.00000E+00 0.38601E-02
UV_bMass_1eps = 0.00000E+00 0.17803E+00
UV_tMass_1eps = 0.00000E+00 0.65530E+01
UVWfct_b_0_1eps -0.18939E-01 0.00000E+00
UVWfct_G_1_1eps -0.31566E-02 0.00000E+00
R2_bxtW = -0.00000E+00 -0.11664E-01
R2_bbX0_h = 0.00000E+00 0.96408E-03
R2_ttX0_h = 0.00000E+00 0.35486E-01
R2_bbX0_a = -0.00000E+00 0.00000E+00
R2_ttX0_a = -0.00000E+00 0.00000E+00
R2_ggX0_hb = 0.00000E+00 -0.16992E-02
R2_ggX0_ht = 0.00000E+00 -0.23022E+01
UV_X0bb_h_1eps 0.00000E+00 0.72306E-03
UV_X0tt_h_1eps 0.00000E+00 0.26615E-01
UV_X0bb_a_1eps -0.00000E+00 0.00000E+00
UV_X0tt_a_1eps -0.00000E+00 0.00000E+00
UV_GQQb = 0.00000E+00 0.22893E-01
UV_GQQt = 0.00000E+00 -0.49437E-02
UV_bMass = 0.00000E+00 0.12932E+01
UV_tMass = 0.00000E+00 0.34467E+00
UVWfct_b_0 = -0.13758E+00 -0.00000E+00
UVWfct_t_0 = -0.99615E-03 -0.00000E+00
UVWfct_G_1 = -0.18721E-01 0.00000E+00
UVWfct_G_2 = 0.40427E-02 0.00000E+00
UV_X0bb_h = 0.00000E+00 0.52523E-02
UV_X0tt_h = 0.00000E+00 0.13999E-02
UV_X0bb_a = -0.00000E+00 0.00000E+00
UV_X0tt_a = -0.00000E+00 0.00000E+00
GC_11 = 0.00000E+00 0.46191E+00
GC_3005h1 = 0.00000E+00 0.52532E+02
GC_3009a = 0.00000E+00 0.00000E+00
GC_3009h = -0.00000E+00 -0.19089E-01
GC_3010a = 0.00000E+00 0.00000E+00
GC_3010h = -0.00000E+00 -0.70263E+00
GC_3005a = -0.00000E+00 -0.00000E+00
GC_3005h2 = -0.00000E+00 -0.00000E+00
GC_3005h3 = -0.00000E+00 -0.00000E+00
GC_3005h4 = -0.00000E+00 -0.00000E+00
Collider parameters:
------
Running at P P machine @ 13000.000000000000 GeV
PDF set = nn23nlo
alpha_s(Mz)= 0.1190 running at 2 loops.
alpha_s(Mz)= 0.1190 running at 2 loops.
Renormalization scale set on event-by-event basis
Factorization scale set on event-by-event basis
INFO: orders_tag_plot is computed as: + QCD * 1 + QNP * 100 + QED * 10000
orders_tag_plot= 60002 for QCD,QNP,QED, = 2 , 0 , 6 ,
AMP_SPLIT: 1 correspond to S.O. 2 0 6
orders_tag_plot= 50102 for QCD,QNP,QED, = 2 , 1 , 5 ,
AMP_SPLIT: 2 correspond to S.O. 2 1 5
orders_tag_plot= 40202 for QCD,QNP,QED, = 2 , 2 , 4 ,
AMP_SPLIT: 3 correspond to S.O. 2 2 4
orders_tag_plot= 60004 for QCD,QNP,QED, = 4 , 0 , 6 ,
AMP_SPLIT: 4 correspond to S.O. 4 0 6
orders_tag_plot= 50104 for QCD,QNP,QED, = 4 , 1 , 5 ,
AMP_SPLIT: 5 correspond to S.O. 4 1 5
orders_tag_plot= 40204 for QCD,QNP,QED, = 4 , 2 , 4 ,
AMP_SPLIT: 6 correspond to S.O. 4 2 4
getting user params
Enter number of events and iterations:
Number of events and iterations -1 12
Enter desired fractional accuracy:
Desired fractional accuracy: 2.9999999999999
Enter alpha, beta for G_soft
Enter alpha<0 to set G_soft=1 (no ME soft)
for G_soft: alpha= 1.0000000000000000 , beta= -0.100000000000
Enter alpha, beta for G_azi
Enter alpha>0 to set G_azi=0 (no azi corr)
for G_azi: alpha= -1.0000000000000000 , beta= -0.100000000000
Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Do MC over helicities for the virtuals
Enter Configuration Number:
Running Configuration Number: 3 0
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
MINT running mode: 0
Set the three folding parameters for MINT
xi_i, y_ij, phi_i
1 1 1
'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
Enter 'born0' or 'virt0' to perform
a pure n-body integration (no S functions)
doing the all of this channel
Normal integration (Sfunction != 1)
about to integrate 13 -1 12 3
imode is 0
channel 1 : 3 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
#------
# FastJet release 3.1.3 [fjcore]
# M. Cacciari, G.P. Salam and G. Soyez
# A software package for jet finding and analysis at colliders
# http://
#
# Please cite EPJC72(2012)1896 [arXiv:1111.6097] if you use this package
# for scientific work and optionally PLB641(2006)57 [hep-ph/0512210].
#
# FastJet is provided without warranty under the terms of the GNU GPLv2.
# It uses T. Chan's closest pair algorithm, S. Fortune's Voronoi code
# and 3rd party plugin jet algorithms. See COPYING file for details.
#------
------- iteration 1
Update # PS points (even_rn): 1040 --> 1040
Using random seed offsets: 3 , 1 , 0
with seed 33
Ranmar initialization seeds 11951 9408
Total number of FKS directories is 10
FKS process map (sum= 3 ) :
1 --> 5 : 1 6 7 8 9
2 --> 2 : 2 10
3 --> 1 : 3
4 --> 1 : 4
5 --> 1 : 5
======
process combination map (specified per FKS dir):
1 map 1 2
1 inv. map 1 2
2 map 1 2
2 inv. map 1 2
3 map 1 2
3 inv. map 1 2
4 map 1 2
4 inv. map 1 2
5 map 1 2
5 inv. map 1 2
6 map 1 1 1 1 2 2 2 2
6 inv. map 4 8
7 map 1 2
7 inv. map 1 2
8 map 1 2
8 inv. map 1 2
9 map 1 1 2 2
9 inv. map 2 4
10 map 1 2
10 inv. map 1 2
======
tau_min 1 1 : 0.30270E+03 -- 0.31270E+03
tau_min 2 1 : 0.30270E+03 -- 0.31270E+03
tau_min 3 1 : 0.31270E+03 0.31270E+03 0.31270E+03
tau_min 4 1 : 0.31270E+03 0.31270E+03 0.31270E+03
tau_min 5 1 : 0.31270E+03 0.31270E+03 0.31270E+03
tau_min 6 1 : 0.30270E+03 -- 0.31270E+03
tau_min 7 1 : 0.30270E+03 -- 0.31270E+03
tau_min 8 1 : 0.30270E+03 -- 0.31270E+03
tau_min 9 1 : 0.30270E+03 -- 0.31270E+03
tau_min 10 1 : 0.30270E+03 -- 0.31270E+03
Scale values (may change event by event):
muR, muR_reference: 0.183667D+03 0.183667D+03 1.00
muF1, muF1_reference: 0.183667D+03 0.183667D+03 1.00
muF2, muF2_reference: 0.183667D+03 0.183667D+03 1.00
QES, QES_reference: 0.183667D+03 0.183667D+03 1.00
muR_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
muF1_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
muF2_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
QES_reference [functional form]:
H_T/2 := sum_i mT(i)/2, i=final state
alpha_s= 0.10756719615258985
BORN: not keeping split order 1
BORN: keeping split order 2
BORN: not keeping split order 3
counterterm S.O 1 QCD
BORN: not keeping split order 1
BORN: keeping split order 2
BORN: not keeping split order 3
counterterm S.O 2 QNP
BORN: not keeping split order 1
BORN: not keeping split order 2
BORN: not keeping split order 3
counterterm S.O 3 QED
BORN: not keeping split order 1
BORN: not keeping split order 2
BORN: not keeping split order 3
Charge-linked born are not used
Color-linked born are used
alpha_s value used for the virtuals is (for the first PS point): 0.10756719615258985
=====
{ }
{ [32m [0m }
{ [32m ,, [0m }
{ [32m`7MMM. ,MMF' `7MM `7MMF' [0m }
{ [32m MMMb dPMM MM MM [0m }
{ [32m M YM ,M MM ,6"Yb. ,M""bMM MM ,pW"Wq. ,pW"Wq.`7MMpdMAo. [0m }
{ [32m M Mb M' MM 8) MM ,AP MM MM 6W' `Wb 6W' `Wb MM `Wb [0m }
{ [32m M YM.P' MM ,pm9MM 8MI MM MM , 8M M8 8M M8 MM M8 [0m }
{ [32m M `YM' MM 8M MM `Mb MM MM ,M YA. ,A9 YA. ,A9 MM ,AP [0m }
{ [32m.JML. `' .JMML.`
{ [32m MM [0m }
{ [32m .JMML. [0m }
{ [32m[0mv3.5.3 (2023-12-23), Ref: arXiv:1103.0621v2, arXiv:1405.
{ [32m [0m }
{ }
=====
======
INFO: MadLoop read these parameters from ../MadLoop5_
======
> MLReductionLib = 6|7|1
> CTModeRun = -1
> MLStabThres = 1.0000000000000
> NRotations_DP = 0
> NRotations_QP = 0
> CTStabThres = 1.0000000000000
> CTLoopLibrary = 2
> CTModeInit = 1
> CheckCycle = 3
> MaxAttempts = 10
> UseLoopFilter = F
> HelicityFilterLevel = 2
> ImprovePSPoint = 2
> DoubleCheckHeli
> LoopInitStartOver = F
> HelInitStartOver = F
> ZeroThres = 1.0000000000000
> OSThres = 1.0000000000000
> WriteOutFilters = T
> UseQPIntegrandF
> UseQPIntegrandF
> IREGIMODE = 2
> IREGIRECY = T
> COLLIERMode = 1
> COLLIERRequired
> COLLIERCanOutput = F
> COLLIERComputeU
> COLLIERComputeI
> COLLIERGlobalCache = -1
> COLLIERUseCache
> COLLIERUseInter
======
-------
| You are using CutTools - Version 1.9.3 |
| Authors: G. Ossola, C. Papadopoulos, R. Pittau |
| Published in JHEP 0803:042,2008 |
| http://
| |
| Compiler with 34 significant digits detetected |
------
#######
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #
+----
| |
| Ninja - version 1.1.0 |
| |
| Author: Tiziano Peraro |
| |
| Based on: |
| |
| P. Mastrolia, E. Mirabella and T. Peraro, |
| "Integrand reduction of one-loop scattering amplitudes |
| through Laurent series expansion," |
| JHEP 1206 (2012) 095 [arXiv:1203.0291 [hep-ph]]. |
| |
| T. Peraro, |
| "Ninja: Automated Integrand Reduction via Laurent |
| Expansion for One-Loop Amplitudes," |
| Comput.Phys.Commun. 185 (2014) [arXiv:1403.1229 [hep-ph]] |
| |
+----
STOP 1
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
#######
#######
# #
# You are using OneLOop-3.6 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 18-02-2015 #
# #
# Please cite #
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
#######
VIRT: not keeping split order 1
VIRT: keeping split order 2
VIRT: not keeping split order 3
Splitorders 5
QCD: 4
QNP: 1
QED: 5
---- POLES CANCELLED ----
COEFFICIENT DOUBLE POLE:
MadFKS: 2.5778610373175
COEFFICIENT SINGLE POLE:
MadFKS: 4.9403200259598
Parameters related to the different treatment of alpha in ME and PDFs
alphascheme 0
nlep_run -1
nupq_run -1
ndnq_run -1
w_run 1
Negative sumMCsec -1.746870865354
Time in seconds: 0
|
Revision history for this message
|
#8 |
This is the issue
w_run 1
Negative sumMCsec -1.746870865354
Time in seconds: 0
I add Rikkert to the thread, as he is more familiar than me with montecarlocounter.f
Thanks a lot
Marco
|
Revision history for this message
|
#9 |
Hello,
This error message means that the MC subtraction term is negative. This should not happen. If the Born matrix element is positive, so are the MC subtraction terms. Are you trying to compute an interference here that could make the Born negative for some phase-space points?
Best,
Rikkert
|
Revision history for this message
|
#10 |
Hi Rikkert,
thanks for your reply!
This interference is destructive / negative on the total cross section indeed.
So, is this approach generally not recommended in such cases?
I was pointed to this older thread https:/
My alternative approach would be to pick several points in the phase space (kHtt & kAtt in the model) and interpolate subsequently. This would also be ok unless issues appear, but possibly less handy for the analysis.
Please let me know your suggestion. Thanks!
Cheers
Florian
|
Revision history for this message
|
#11 |
Hello Florian,
Event generation for pure interference contributions is not supported at NLO.
Best,
Rikkert
|
Revision history for this message
|
#12 |
Hello Rikkert,
alright, I will not go on with this approach then.
Thanks for the information!
Best,
Florian
