EWK W production @ NLO
Dear all,
I'm trying to produce an EWK W production @ NLO QCD. I have read people trying to do something similar to me but with different problems from the one I'm facing.
I'm splitting the process into two different states of charge. There is no strong reason to do that, but I don't think that can be an issue. So with the w+ I have:
set low_mem_
import model loop_sm-no_b_mass
define l+ = e+ mu+ ta+
define vl = ve vm vt
generate p p > l+ vl j j $$ t t~ h QED=4 QCD=0 [QCD]
output ewkwplvjj_NLO
If I'm not wrong, seems to do all the computations at the matrix element, but then it fails some tests, and finally the error with pilotrun. Now, unfortunately, I would have to re-run the process to get the mentioned log file, because it gets automatically cancelled by the framework (not by Madgraph itself). So if it is needed, I can re-try. In general, if it is possible, I would like to know what these tests are and why they fail. Is it due to some parameter settings?
Is it suggested to try with a stable w and then decay it with madspin? Is there any physical difference from what I'm doing right now?
Cheers
AB
INFO: Generate web pages
INFO: Compiling StdHEP. This has to be done only once.
Start waiting for update. (more info in debug mode)
[1;34mWARNING: StdHep failed to compiled. This forbids to run NLO+PS with PY6 and Herwig6 [0m
INFO: details on the compilation error are available on /home/abulla/
INFO: if you want to retry the compilation automatically, you have to remove that file first
Type "launch" to generate events from this process, or see
/home/abulla/
Run "open index.html" to see more information about this process.
display multiparticles
Multiparticle labels:
p = g u c d s u~ c~ d~ s~ b b~
j = g u c d s b u~ c~ d~ s~ b~
l+ = e+ mu+ ta+
l- = e- mu-
vl = ve vm vt
vl~ = ve~ vm~ vt~
all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
quit
copying generated process to working directory
WARNING: You've chosen not to use the PDF sets recommended for run3 production!
If this isn't intentional, and you prefer to use the recommended sets,
insert the following lines into your process-
copying run_card.dat file
starting NLO mode
No module named 'madgraph'
No hepmc reader since No module named 'madgraph' [1;30m[
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 5.2.9.18 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /home/abulla/
INFO: load configuration from /home/abulla/
INFO: load configuration from /home/abulla/
Using default text editor "vi". Set another one in ./input/
No valid eps viewer found. Please set in ./input/
No valid web browser found. Please set in ./input/
launch -n pilotrun
[1;34mrun_card missed argument pt_min_pdg. Takes default: {}[0m
[1;34mrun_card missed argument pt_max_pdg. Takes default: {}[0m
[1;34mrun_card missed argument mxx_min_pdg. Takes default: {}[0m
stty: standard input: Inappropriate ioctl for device
The following switches determine which programs are run:
/======
1. Type of perturbative computation [1m order[0m = [32mNLO[0m
2. No MC@[N]LO matching / event generation [1mfixed_order[0m = [32mOFF[0m
3. Shower the generated events [1m shower[0m = [31mPYTHIA8[0m
4. Decay onshell particles [1m madspin[0m = [31mOFF[0m
5. Add weights to events for new hypp. [1m reweight[0m = [31mOFF[0m
6. Run MadAnalysis5 on the events generated [1mmadanalysis[0m = [01mNot Avail.[0m
\======
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]
>stty: standard input: Inappropriate ioctl for device
The following switches determine which programs are run:
/======
1. Type of perturbative computation [1m order[0m = [32mNLO[0m
2. No MC@[N]LO matching / event generation [1mfixed_order[0m = [32mOFF[0m
3. Shower the generated events [1m shower[0m = [32mOFF[0m
4. Decay onshell particles [1m madspin[0m = [31mOFF[0m
5. Add weights to events for new hypp. [1m reweight[0m = [31mOFF[0m
6. Run MadAnalysis5 on the events generated [1mmadanalysis[0m = [01mNot Avail.[0m
\======
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.
>stty: standard input: Inappropriate ioctl for device
The following switches determine which programs are run:
/======
1. Type of perturbative computation [1m order[0m = [32mNLO[0m
2. No MC@[N]LO matching / event generation [1mfixed_order[0m = [32mOFF[0m
3. Shower the generated events [1m shower[0m = [32mOFF[0m
4. Decay onshell particles [1m madspin[0m = [31mOFF[0m
5. Add weights to events for new hypp. [1m reweight[0m = [31mOFF[0m
6. Run MadAnalysis5 on the events generated [1mmadanalysis[0m = [01mNot Avail.[0m
\======
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.
>stty: standard input: Inappropriate ioctl for device
INFO: will run in mode: noshower
[1;34mWARNING: 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. [0m
[1;34mrun_card missed argument pt_min_pdg. Takes default: {}[0m
[1;34mrun_card missed argument pt_max_pdg. Takes default: {}[0m
[1;34mrun_card missed argument mxx_min_pdg. Takes default: {}[0m
Do you want to edit a card (press enter to bypass editing)?
/------
| [31m1[0m. param : [32mparam_
| [31m2[0m. run : [32mrun_
\------
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
[[4m0[0m, done, 1, param, 2, run, enter path][90s to answer]
>INFO: Update the dependent parameter of the param_card.dat
[1;34mrun_card missed argument pt_min_pdg. Takes default: {}[0m
[1;34mrun_card missed argument pt_max_pdg. Takes default: {}[0m
[1;34mrun_card missed argument mxx_min_pdg. Takes default: {}[0m
INFO: Starting run
INFO: Compiling the code
[1;34mrun_card missed argument pt_min_pdg. Takes default: {}[0m
[1;34mrun_card missed argument pt_max_pdg. Takes default: {}[0m
[1;34mrun_card missed argument mxx_min_pdg. Takes default: {}[0m
run_card missed argument shower_
run_card missed argument muf_ref_fixed. Takes default: 91.188
run_card missed argument dynamical_
run_card missed argument muf_over_ref. Takes default: 1.0
run_card missed argument rw_rscale. Takes default: [1.0, 2.0, 0.5]
run_card missed argument rw_fscale. Takes default: [1.0, 2.0, 0.5]
run_card missed argument iappl. Takes default: 0
INFO: Using LHAPDF v6.4.0 interface for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 48 cores
INFO: Compiling P0_uxu_tapvtscx...
INFO: Compiling P0_us_epvesd...
INFO: Compiling P0_uc_tapvtcd...
INFO: Compiling P0_udx_tapvtbbx...
INFO: Compiling P0_usx_tapvtucx...
INFO: Compiling P0_udx_tapvtssx...
INFO: Compiling P0_ucx_tapvtcxd...
INFO: Compiling P0_ddx_tapvtscx...
INFO: Compiling P0_ubx_tapvtbxd...
INFO: Compiling P0_dc_tapvtds...
INFO: Compiling P0_ub_tapvtbd...
INFO: Compiling P0_dsx_tapvtdcx...
INFO: Compiling P0_ddx_epvescx...
INFO: Compiling P0_usx_tapvtsxd...
INFO: Compiling P0_uux_tapvtscx...
INFO: Compiling P0_uc_tapvtus...
INFO: Compiling P0_us_tapvtsd...
INFO: Compiling P0_udx_epveccx...
INFO: Compiling P0_uxu_epvescx...
INFO: Compiling P0_dc_epveds...
INFO: Compiling P0_udx_tapvtccx...
INFO: Compiling P0_udx_epvebbx...
INFO: Compiling P0_ub_epvebd...
INFO: Compiling P0_ucx_epvecxd...
INFO: Compiling P0_uc_epvecd...
INFO: Compiling P0_udx_epvessx...
INFO: Compiling P0_dsx_epvedcx...
INFO: Compiling P0_usx_epvesxd...
INFO: Compiling P0_uc_epveus...
INFO: Compiling P0_ubx_epvebxd...
INFO: Compiling P0_uux_epvescx...
INFO: Compiling P0_usx_epveucx...
INFO: Compiling P0_du_tapvtdd...
INFO: Compiling P0_uxu_tapvtuxd...
INFO: Compiling P0_ud_tapvtdd...
INFO: Compiling P0_uux_epveuxd...
INFO: Compiling P0_ddx_tapvtdux...
INFO: Compiling P0_ud_epvedd...
INFO: Compiling P0_du_epvedd...
INFO: Compiling P0_uux_tapvtuxd...
INFO: Compiling P0_udx_epveuux...
INFO: Compiling P0_ddx_epvedux...
INFO: Compiling P0_uxu_epveuxd...
INFO: Compiling P0_udx_tapvtuux...
INFO: Compiling P0_uu_tapvtud...
INFO: Compiling P0_uu_epveud...
INFO: Compiling P0_udx_epveddx...
INFO: Compiling P0_udx_tapvtddx...
INFO: P0_uc_epveus done.
INFO: Compiling P0_uxsx_
INFO: P0_ddx_tapvtscx done.
INFO: Compiling P0_dxu_tapvtccx...
INFO: P0_uux_tapvtscx done.
INFO: Compiling P0_dxu_tapvtbbx...
INFO: P0_us_epvesd done.
INFO: Compiling P0_uxsx_epveuxcx...
INFO: P0_udx_epveccx done.
INFO: Compiling P0_dxcx_
INFO: P0_uc_epvecd done.
INFO: Compiling P0_dxu_epvessx...
INFO: P0_dc_tapvtds done.
INFO: Compiling P0_dxc_tapvtcux...
INFO: P0_udx_epveuux done.
INFO: Compiling P0_dxc_tapvtdxs...
INFO: P0_udx_epvessx done.
INFO: Compiling P0_uxc_tapvtuxs...
INFO: P0_usx_tapvtsxd done.
INFO: Compiling P0_dxu_epvebbx...
INFO: P0_ucx_epvecxd done.
INFO: Compiling P0_dxcx_epvecxux...
INFO: P0_uux_epvescx done.
INFO: Compiling P0_dxu_tapvtssx...
INFO: P0_udx_epvebbx done.
INFO: Compiling P0_dxs_tapvtsux...
INFO: P0_usx_tapvtucx done.
INFO: Compiling P0_dxc_epvedxs...
INFO: P0_ucx_tapvtcxd done.
INFO: Compiling P0_dxd_tapvtscx...
INFO: P0_udx_tapvtddx done.
INFO: Compiling P0_dxs_epvesux...
INFO: P0_usx_epvesxd done.
INFO: Compiling P0_dxu_epveccx...
INFO: P0_du_epvedd done.
INFO: Compiling P0_dxc_epvecux...
INFO: P0_uxu_tapvtscx done.
INFO: Compiling P0_dxd_epvescx...
INFO: P0_udx_tapvtssx done.
INFO: Compiling P0_uxc_epveuxs...
INFO: P0_uux_epveuxd done.
INFO: Compiling P0_dxb_tapvtbux...
INFO: P0_uc_tapvtus done.
INFO: Compiling P0_dxsx_
INFO: P0_dc_epveds done.
INFO: P0_udx_tapvtbbx done.
INFO: Compiling P0_dxsx_epvesxux...
INFO: Compiling P0_dxu_tapvtuux...
INFO: P0_uxu_epveuxd done.
INFO: Compiling P0_bu_tapvtbd...
INFO: P0_ubx_tapvtbxd done.
INFO: Compiling P0_dxsx_
INFO: P0_usx_epveucx done.
INFO: Compiling P0_dxbx_
INFO: P0_ub_epvebd done.
INFO: Compiling P0_dxb_epvebux...
INFO: P0_ddx_epvescx done.
INFO: Compiling P0_dxsx_epvedxcx...
INFO: P0_du_tapvtdd done.
INFO: Compiling P0_uxdx_
INFO: P0_udx_tapvtccx done.
INFO: Compiling P0_bdx_tapvtbux...
INFO: P0_us_tapvtsd done.
INFO: Compiling P0_uxdx_epveuxux...
INFO: P0_uc_tapvtcd done.
INFO: Compiling P0_dxbx_epvebxux...
INFO: P0_ddx_epvedux done.
INFO: Compiling P0_dxux_epveuxux...
INFO: P0_uu_tapvtud done.
INFO: Compiling P0_dxu_epveddx...
INFO: P0_ub_tapvtbd done.
INFO: Compiling P0_bbx_tapvtdux...
INFO: P0_ubx_epvebxd done.
INFO: Compiling P0_bdx_epvebux...
INFO: P0_dsx_tapvtdcx done.
INFO: Compiling P0_dxux_
INFO: P0_udx_epveddx done.
INFO: Compiling P0_dxu_tapvtddx...
INFO: P0_ddx_tapvtdux done.
INFO: Compiling P0_dxdx_epvedxux...
INFO: P0_dsx_epvedcx done.
INFO: Compiling P0_bu_epvebd...
INFO: P0_uxu_epvescx done.
INFO: Compiling P0_dxd_epvedux...
INFO: P0_uu_epveud done.
INFO: Compiling P0_bxu_tapvtbxd...
INFO: P0_ud_tapvtdd done.
INFO: Compiling P0_dxu_epveuux...
INFO: P0_ud_epvedd done.
INFO: Compiling P0_dxd_tapvtdux...
INFO: P0_uxu_tapvtuxd done.
INFO: Compiling P0_dxdx_
INFO: P0_uux_tapvtuxd done.
INFO: Compiling P0_bbx_epvedux...
INFO: P0_udx_tapvtuux done.
INFO: Compiling P0_bxu_epvebxd...
INFO: P0_dxc_tapvtdxs done.
INFO: Compiling P0_bxdx_
INFO: P0_dxu_tapvtccx done.
INFO: Compiling P0_bxdx_epvebxux...
INFO: P0_dxu_tapvtssx done.
INFO: Compiling P0_bxb_tapvtdux...
INFO: P0_dxc_tapvtcux done.
INFO: Compiling P0_bxb_epvedux...
INFO: P0_dxu_epvebbx done.
INFO: P0_dxcx_tapvtcxux done.
INFO: P0_dxs_tapvtsux done.
INFO: P0_uxsx_tapvtuxcx done.
INFO: P0_dxu_tapvtbbx done.
INFO: P0_dxc_epvecux done.
INFO: P0_dxsx_epvedxcx done.
INFO: P0_uxsx_epveuxcx done.
INFO: P0_dxb_epvebux done.
INFO: P0_uxc_tapvtuxs done.
INFO: P0_dxcx_epvecxux done.
INFO: P0_dxs_epvesux done.
INFO: P0_dxc_epvedxs done.
INFO: P0_dxb_tapvtbux done.
INFO: P0_dxu_epvessx done.
INFO: P0_bbx_tapvtdux done.
INFO: P0_dxsx_epvesxux done.
INFO: P0_dxux_epveuxux done.
INFO: P0_dxsx_tapvtsxux done.
INFO: P0_dxd_tapvtscx done.
INFO: P0_dxbx_epvebxux done.
INFO: P0_bu_tapvtbd done.
INFO: P0_bdx_epvebux done.
INFO: P0_dxbx_tapvtbxux done.
INFO: P0_uxc_epveuxs done.
INFO: P0_dxsx_tapvtdxcx done.
INFO: P0_bxu_tapvtbxd done.
INFO: P0_dxd_epvescx done.
INFO: P0_dxux_tapvtuxux done.
INFO: P0_dxdx_tapvtdxux done.
INFO: P0_uxdx_epveuxux done.
INFO: P0_dxu_epveccx done.
INFO: P0_bu_epvebd done.
INFO: P0_bdx_tapvtbux done.
INFO: P0_dxdx_epvedxux done.
INFO: P0_bxu_epvebxd done.
INFO: P0_bbx_epvedux done.
INFO: P0_dxu_tapvtddx done.
INFO: P0_uxdx_tapvtuxux done.
INFO: P0_dxd_epvedux done.
INFO: P0_dxd_tapvtdux done.
INFO: P0_dxu_epveddx done.
INFO: P0_dxu_tapvtuux done.
INFO: P0_dxu_epveuux done.
INFO: P0_bxdx_epvebxux done.
INFO: P0_bxdx_tapvtbxux done.
INFO: P0_bxb_tapvtdux done.
INFO: P0_bxb_epvedux done.
INFO: Checking test output:
INFO: P0_uxu_tapvtscx
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_us_epvesd
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_uc_tapvtcd
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_udx_tapvtbbx
INFO: Result for test_ME:
[1mINFO: Output of the failing test:
==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.4.0 loading /cvmfs/
NNPDF30_
Enter 0 to compute MC/MC(limit)
1 to compute MC/ME(limit)
2 to compute ME/ME(limit)
Enter xi_i, y_ij to be used in coll/soft tests
Enter -2 to generate them randomly
Enter number of tests for soft and collinear limits
A PDF is used, so alpha_s(MZ) is going to be modified
Old value of alpha_s from param_card: 0.11799999999999999
New value of alpha_s from PDF lhapdf : 0.11800222660863767
Give FKS configuration number ("0" loops over all)
======
NEW FKS CONFIGURATION:
FKS configuration number is 1
FKS partons are: i= 7 j= 1
with PDGs: i= 21 j= 2
Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 1 , 26 , 0
with seed 32
Ranmar initialization seeds 10730 9432
tau_min 1 1 : 0.17770E+01 -- 0.80439E+02
tau_min 2 1 : 0.17770E+01 -- 0.80439E+02
tau_min 3 1 : 0.17970E+01 0.17970E+01 0.80439E+02
tau_min 4 1 : 0.17970E+01 0.17970E+01 0.80439E+02
tau_min 5 1 : 0.17770E+01 -- 0.80439E+02
tau_min 6 1 : 0.17770E+01 -- 0.80439E+02
Soft tests done for (Born) config 1
Failures: 24
Soft test 1 PASSED. Fraction of failures: 0.24
Collinear tests done for (Born) config 1
Failures: 0
Collinear test 1 PASSED. Fraction of failures: 0.00
======
NEW FKS CONFIGURATION:
FKS configuration number is 2
FKS partons are: i= 7 j= 2
with PDGs: i= 21 j= -1
Enter graph number (iconfig), '0' loops over all graphs
Soft tests done for (Born) config 1
Failures: 23
Soft test 2 PASSED. Fraction of failures: 0.23
Collinear tests done for (Born) config 1
Failures: 0
Collinear test 2 PASSED. Fraction of failures: 0.00
======
NEW FKS CONFIGURATION:
FKS configuration number is 3
FKS partons are: i= 7 j= 5
with PDGs: i= 21 j= 5
Enter graph number (iconfig), '0' loops over all graphs
Soft tests done for (Born) config 1
Failures: 6
Soft test 3 PASSED. Fraction of failures: 0.06
Collinear tests done for (Born) config 1
Failures: 99
Collinear test 3 FAILED. Fraction of failures: 0.99
======
NEW FKS CONFIGURATION:
FKS configuration number is 4
FKS partons are: i= 7 j= 6
with PDGs: i= 21 j= -5
Enter graph number (iconfig), '0' loops over all graphs
Soft tests done for (Born) config 1
Failures: 5
Soft test 4 PASSED. Fraction of failures: 0.05
Collinear tests done for (Born) config 1
Failures: 99
Collinear test 4 FAILED. Fraction of failures: 0.99
======
NEW FKS CONFIGURATION:
FKS configuration number is 5
FKS partons are: i= 7 j= 1
with PDGs: i= -2 j= 21
Enter graph number (iconfig), '0' loops over all graphs
Soft tests done for (Born) config 1
Failures: 0
Soft test 5 PASSED. Fraction of failures: 0.00
Collinear tests done for (Born) config 1
Failures: 2
Collinear test 5 PASSED. Fraction of failures: 0.02
======
NEW FKS CONFIGURATION:
FKS configuration number is 6
FKS partons are: i= 6 j= 2
with PDGs: i= 1 j= 21
Enter graph number (iconfig), '0' loops over all graphs
Soft tests done for (Born) config 1
Failures: 0
Soft test 6 PASSED. Fraction of failures: 0.00
Thanks for using LHAPDF 6.4.0. Please make sure to cite the paper:
Eur.Phys.J. C75 (2015) 3, 132 (http://
Collinear tests done for (Born) config 1
Failures: 1
Collinear test 6 PASSED. Fraction of failures: 0.01 [0m
[1;31mError detected in "launch -n pilotrun"
write debug file /home/abulla/
If you need help with this issue please contact us on https:/
str : Some tests failed, run cannot continue.
Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.[0m
quit
INFO:
finished pilot run
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
This question was reopened
- by Andrea
Can you help with this problem?
Provide an answer of your own, or ask Andrea for more information if necessary.