mgv2.6.1 madspin error
Hi Olivier and Madgraph team
I've been running several test samples with the new released MGv2.6.1
I haven't changed any of the card infos.
Only thing I did was
generate p p > t t~ [QCD]
with MadSpin on (default madspin card decaying t, t~)
Please take a look at this
MG5_aMC>generate p p > t t~ [QCD]
INFO: Generating FKS-subtracted matrix elements for born process: g g > t t~ [ all = QCD ] (1 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: u u~ > t t~ [ all = QCD ] (2 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: c c~ > t t~ [ all = QCD ] (3 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: d d~ > t t~ [ all = QCD ] (4 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: s s~ > t t~ [ all = QCD ] (5 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: u~ u > t t~ [ all = QCD ] (6 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: c~ c > t t~ [ all = QCD ] (7 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: d~ d > t t~ [ all = QCD ] (8 / 9)
INFO: Generating FKS-subtracted matrix elements for born process: s~ s > t t~ [ all = QCD ] (9 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: g g > t t~ QCD<=2 [ all = QCD ] (1 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: u u~ > t t~ QCD<=2 [ all = QCD ] (2 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: c c~ > t t~ QCD<=2 [ all = QCD ] (3 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: d d~ > t t~ QCD<=2 [ all = QCD ] (4 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: s s~ > t t~ QCD<=2 [ all = QCD ] (5 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: u~ u > t t~ QCD<=2 [ all = QCD ] (6 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: c~ c > t t~ QCD<=2 [ all = QCD ] (7 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: d~ d > t t~ QCD<=2 [ all = QCD ] (8 / 9)
INFO: Generating virtual matrix elements using MadLoop:
INFO: Generating virtual matrix element with MadLoop for process: s~ s > t t~ QCD<=2 [ all = QCD ] (9 / 9)
INFO: Generated 9 subprocesses with 136 real emission diagrams, 11 born diagrams and 124 virtual diagrams
MG5_aMC>output TEST
INFO: Writing out the aMC@NLO code, using optimized Loops
INFO: directory /home/shjeon/
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n][60s to answer]
>y
INFO: initialize a new directory: TEST
INFO: remove old information in TEST
INFO: Generating real emission matrix-elements...
INFO: Generating Helas calls for FKS process: g g > t t~ [ all = QCD ] (1 / 9)
INFO: Processing color information for process: g g > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: u u~ > t t~ [ all = QCD ] (2 / 9)
INFO: Processing color information for process: u u~ > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: c c~ > t t~ [ all = QCD ] (3 / 9)
INFO: Generating Helas calls for FKS process: d d~ > t t~ [ all = QCD ] (4 / 9)
INFO: Generating Helas calls for FKS process: s s~ > t t~ [ all = QCD ] (5 / 9)
INFO: Generating Helas calls for FKS process: u~ u > t t~ [ all = QCD ] (6 / 9)
INFO: Processing color information for process: u~ u > t t~ [ all = QCD ]
INFO: Generating Helas calls for FKS process: c~ c > t t~ [ all = QCD ] (7 / 9)
INFO: Generating Helas calls for FKS process: d~ d > t t~ [ all = QCD ] (8 / 9)
INFO: Generating Helas calls for FKS process: s~ s > t t~ [ all = QCD ] (9 / 9)
INFO: ... Done
INFO: Writing files in P0_gg_ttx (1 / 3)
INFO: Creating files in directory V0_gg_ttx
ALOHA: aloha creates FFV1 set of routines with options: L2
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > t t~ QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: g g > t t~ QCD<=2 [ all = QCD ]
INFO: Writing files in P0_uux_ttx (2 / 3)
INFO: Creating files in directory V0_uux_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > t t~ QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: u u~ > t t~ QCD<=2 [ all = QCD ]
INFO: Writing files in P0_uxu_ttx (3 / 3)
INFO: Creating files in directory V0_uxu_ttx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u~ u > t t~ QCD<=2 [ all = QCD ]
INFO: Generating born Feynman diagrams for Process: u~ u > t t~ QCD<=2 [ all = QCD ]
History written to /home/shjeon/
ALOHA: aloha creates FFV1 routines
ALOHA: aloha creates VVVV4 set of routines with options: P0
ALOHA: aloha creates VVV1 routines
ALOHA: aloha creates R2_GG_3 routines
ALOHA: aloha creates R2_GG_2 routines
ALOHA: aloha creates R2_GG_1 routines
ALOHA: aloha creates R2_QQ_2 routines
ALOHA: aloha creates R2_QQ_1 routines
ALOHA: aloha creates R2_GG_1_R2_GG_2 routines
ALOHA: aloha creates R2_GG_1_R2_GG_3 routines
ALOHA: aloha creates R2_QQ_1_R2_QQ_2 routines
save configuration file to /home/shjeon/
INFO: The value for lhapdf in the current configuration does not correspond to a valid executable.
Please set it correctly either in input/mg5_
To avoid regeneration, edit the /home/shjeon/
Note that you can still compile and run aMC@NLO with the built-in PDFs
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Type "launch" to generate events from this process, or see
/home/shjeon/
Run "open index.html" to see more information about this process.
MG5_aMC>launch
INFO: *******
* *
* W E L C O M E to M A D G R A P H 5 *
* a M C @ N L O *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.6.1 2017-12-12 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /home/shjeon/
INFO: load configuration from /home/shjeon/
INFO: load configuration from /home/shjeon/
Using default text editor "vi". Set another one in ./input/
launch auto
The following switches determine which programs are run:
/======
| 1. Type of perturbative computation order = NLO |
| 2. No MC@[N]LO matching / event generation fixed_order = OFF |
| 3. Shower the generated events shower = HERWIG6 |
| 4. Decay onshell particles madspin = OFF |
| 5. Add weights to events for new hypp. reweight = Not Avail. |
| 6. Run MadAnalysis5 on the events generated madanalysis = Not Avail. |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'fixed_order=ON' 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:
/======
| 1. Type of perturbative computation order = NLO |
| 2. No MC@[N]LO matching / event generation fixed_order = OFF |
| 3. Shower the generated events shower = OFF |
| 4. Decay onshell particles madspin = OFF |
| 5. Add weights to events for new hypp. reweight = Not Avail. |
| 6. Run MadAnalysis5 on the events generated madanalysis = Not Avail. |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'fixed_order=ON' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>4
The following switches determine which programs are run:
/======
| 1. Type of perturbative computation order = NLO |
| 2. No MC@[N]LO matching / event generation fixed_order = OFF |
| 3. Shower the generated events shower = OFF |
| 4. Decay onshell particles madspin = ON |
| 5. Add weights to events for new hypp. reweight = Not Avail. |
| 6. Run MadAnalysis5 on the events generated madanalysis = Not Avail. |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'fixed_order=ON' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>0
INFO: will run in mode: noshower
WARNING: You have chosen not to run a parton shower.
NLO events without showering are NOT physical.
Please, shower the LesHouches events before using them for physics analyses.
You have to choose NOW which parton-shower you WILL use and specify it in the run_card.
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_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, 3, madspin, enter path][90s to answer]
>3
open /home/shjeon/
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_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, 3, madspin, enter path]
>1
open /home/shjeon/
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_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, 3, madspin, enter path]
>2
open /home/shjeon/
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_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, 3, madspin, enter path]
>0
INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
INFO: Compiling the code
INFO: For gauge cancellation, the width of 't' has been set to zero.
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
libGL error: No matching fbConfigs or visuals found
libGL error: failed to load driver: swrast
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 8 cores
INFO: Compiling P0_gg_ttx...
INFO: Compiling P0_uux_ttx...
INFO: Compiling P0_uxu_ttx...
INFO: P0_uxu_ttx done.
INFO: P0_uux_ttx done.
INFO: P0_gg_ttx done.
INFO: Checking test output:
INFO: P0_gg_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uux_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxu_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 8 cores
INFO: Cleaning previous results
INFO: Generating events without running the shower.
INFO: Setting up grids
INFO: Idle: 0, Running: 5, Completed: 0 [ current time: 14h11 ]
INFO: Idle: 0, Running: 4, Completed: 1 [ 1.7s ]
INFO: Idle: 0, Running: 2, Completed: 3 [ 3.3s ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 3.4s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 5.7s ]
INFO: Determining the number of unweighted events per channel
Intermediate results:
Random seed: 33
Total cross section: 6.839e+02 +- 6.1e+00 pb
Total abs(cross section): 1.188e+03 +- 1.0e+01 pb
INFO: Computing upper envelope
INFO: Idle: 0, Running: 5, Completed: 0 [ current time: 14h11 ]
INFO: Idle: 0, Running: 4, Completed: 1 [ 5.1s ]
INFO: Idle: 0, Running: 2, Completed: 3 [ 7.6s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 8.2s ]
INFO: Updating the number of unweighted events per channel
Intermediate results:
Random seed: 33
Total cross section: 6.854e+02 +- 4.2e+00 pb
Total abs(cross section): 1.194e+03 +- 5.5e+00 pb
INFO: Generating events
INFO: Idle: 0, Running: 3, Completed: 2 [ current time: 14h11 ]
INFO: Idle: 0, Running: 2, Completed: 3 [ 1.6s ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 10.4s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 11.5s ]
INFO: Doing reweight
INFO: Idle: 0, Running: 2, Completed: 3 [ current time: 14h11 ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 0.41s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 0.45s ]
INFO: Collecting events
INFO:
----
Summary:
Process p p > t t~ [QCD]
Run at p-p collider (6500.0 + 6500.0 GeV)
Number of events generated: 10000
Total cross section: 6.854e+02 +- 4.2e+00 pb
----
Scale variation (computed from LHE events):
----
INFO: The /home/shjeon/
INFO: Events generated
reweight -from_cards
decay_events -from_cards
INFO: Running MadSpin
INFO: This functionality allows for the decay of resonances
INFO: in a .lhe file, keeping track of the spin correlation effets.
INFO: BE AWARE OF THE CURRENT LIMITATIONS:
INFO: (1) Only a succession of 2 body decay are currently allowed
*******
* *
* W E L C O M E to M A D S P I N *
* *
*******
INFO: Extracting the banner ...
INFO: process: p p > t t~
INFO: options:
INFO: detected model: loop_sm. Loading...
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
launch
INFO: Will use seed 637106503
WARNING: set the mass of the b-quark to its value in the param_card.dat: 4.7 GeV
INFO: We need to recalculate the branching fractions for t~,w-,z,w+,t
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /home/shjeon/
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO:
INFO: decay channels for t~ : ( width = 1.4915 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b~ w-
INFO:
INFO:
INFO: decay channels for w- : ( width = 2.04793 GeV )
INFO: BR d1 d2
INFO: 3.333610e-01 d u~
INFO: 3.333610e-01 s c~
INFO: 1.111195e-01 e- ve~
INFO: 1.111195e-01 mu- vm~
INFO: 1.110390e-01 ta- vt~
INFO:
INFO:
INFO: decay channels for z : ( width = 2.441756 GeV )
INFO: BR d1 d2
INFO: 1.523658e-01 d~ d
INFO: 1.523658e-01 s~ s
INFO: 1.507440e-01 b~ b
INFO: 1.188141e-01 u~ u
INFO: 1.188141e-01 c~ c
INFO: 6.793744e-02 ve~ ve
INFO: 6.793744e-02 vm~ vm
INFO: 6.793744e-02 vt~ vt
INFO: 3.438713e-02 e+ e-
INFO: 3.438713e-02 mu+ mu-
INFO: 3.430972e-02 ta+ ta-
INFO:
INFO:
INFO: decay channels for w+ : ( width = 2.04793 GeV )
INFO: BR d1 d2
INFO: 3.333610e-01 d~ u
INFO: 3.333610e-01 s~ c
INFO: 1.111195e-01 e+ ve
INFO: 1.111195e-01 mu+ vm
INFO: 1.110390e-01 ta+ vt
INFO:
INFO:
INFO: decay channels for t : ( width = 1.4915 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b w+
INFO:
INFO: generating the production square matrix element
INFO: generate p p > t t~ --no_warning=
INFO: Done 1.374
INFO: generating the full matrix element squared (with decay)
INFO: generate p p > t t~ --no_warning=
Command "launch auto " interrupted with error:
InvalidCmd : "define" command requires symbols "=" at the second position
quit
INFO:
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply: