mgv2.6.1 madspin error

Asked by Sihyun Jeon on 2017-12-20

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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST already exists.
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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/proc_card_mg5.dat
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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt
INFO: The value for lhapdf in the current configuration does not correspond to a valid executable.
Please set it correctly either in input/mg5_configuration or with "set lhapdf /path/to/lhapdf-config" and regenrate the process.
To avoid regeneration, edit the /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt file.
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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/README
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://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt
INFO: load configuration from /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/input/mg5_configuration.txt
INFO: load configuration from /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
launch auto
The following switches determine which programs are run:
/=================================================================================\
| 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/MadWidth/...).
     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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/madspin_card.dat
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/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madspin, enter path]
>1
open /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/param_card.dat
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/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madspin, enter path]
>2
open /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/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 |
| 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/MadWidth/...).
     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):
          Dynamical_scale_choice -1 (envelope of 9 values):
              6.695e+02 pb +10.2% -11.9%
   --------------------------------------------------------------

INFO: The /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Events/run_01/events.lhe.gz file has been generated.

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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/models/loop_sm with file models/loop_sm/restrict_default.dat .
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=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > t t~ pert_QCD --no_warning=duplicate;
INFO: Done 1.374
INFO: generating the full matrix element squared (with decay)
INFO: generate p p > t t~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21, (t~ > b~ w- , w- > all all QCD=99), (t > b w+ , w+ > all all QCD=99) --no_warning=duplicate;add process p p > t t~ pert_QCD --no_warning=duplicate, (t~ > b~ w- , w- > all all QCD=99), (t > b w+ , w+ > all all QCD=99) --no_warning=duplicate;
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
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2017-12-22
Last query:
2017-12-22
Last reply:
2017-12-21

Hi,

Yeah I discover that couple of days ago.
This is already fixed in 2.6.2 that will be release early January.
I have also added some unitest and acceptance test to avoid such problem to occur again in the future.

Sorry for that,

Olivier

PS: You can use bzr to download the development version via

bzr branch lp:~mg5core1/mg5amcnlo/2.6.2<https://code.launchpad.net/~mg5core1/mg5amcnlo/2.6.2>

On 20 Dec 2017, at 06:17, Sihyun Jeon <<email address hidden><mailto:<email address hidden>>> wrote:

New question #661988 on MadGraph5_aMC@NLO:
https://answers.launchpad.net/mg5amcnlo/+question/661988

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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST already exists.
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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/proc_card_mg5.dat
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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt
INFO: The value for lhapdf in the current configuration does not correspond to a valid executable.
Please set it correctly either in input/mg5_configuration or with "set lhapdf /path/to/lhapdf-config" and regenrate the process.
To avoid regeneration, edit the /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt file.
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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/README
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://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt
INFO: load configuration from /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/input/mg5_configuration.txt
INFO: load configuration from /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/amcatnlo_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
launch auto
The following switches determine which programs are run:
/=================================================================================\
| 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/MadWidth/...).
    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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/madspin_card.dat
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/MadWidth/...).
    Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, enter path]
1
open /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/param_card.dat
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/MadWidth/...).
    Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, enter path]
2
open /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Cards/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 |
| 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/MadWidth/...).
    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):
         Dynamical_scale_choice -1 (envelope of 9 values):
             6.695e+02 pb +10.2% -11.9%
  --------------------------------------------------------------

INFO: The /home/shjeon/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/TEST/Events/run_01/events.lhe.gz file has been generated.

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/Desktop/Physics/MadGraph/MG5_aMC_v2_6_1/models/loop_sm with file models/loop_sm/restrict_default.dat .
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=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > t t~ pert_QCD --no_warning=duplicate;
INFO: Done 1.374
INFO: generating the full matrix element squared (with decay)
INFO: generate p p > t t~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21, (t~ > b~ w- , w- > all all QCD=99), (t > b w+ , w+ > all all QCD=99) --no_warning=duplicate;add process p p > t t~ pert_QCD --no_warning=duplicate, (t~ > b~ w- , w- > all all QCD=99), (t > b w+ , w+ > all all QCD=99) --no_warning=duplicate;
Command "launch auto " interrupted with error:
InvalidCmd : "define" command requires symbols "=" at the second position
quit
INFO:

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

Sihyun Jeon (shjeon) said : #2

Thanks Olivier Mattelaer, that solved my question.