error detected in generating events
| Affects | Status | Importance | Assigned to | Milestone | |
|---|---|---|---|---|---|
| MadGraph5_aMC@NLO |
Invalid
|
Undecided
|
Unassigned | ||
Bug Description
Hi,
On trying to generate process p p > h > mu+ mu- [noborn = NP], I get the diagrams correct but after launch I get some error.
The debug file is attached.
MG5_aMC>generate p p > h > mu+ mu-[noborn = NP]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS NP+2*QED+QCD
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 6 loops, 0 R2, 0 UV
INFO: Contributing diagrams generated: 0 Born, 3 loops, 0 R2, 0 UV
INFO: Contributing diagrams generated: 0 Born, 3 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 6 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 3 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (NP)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 3 loops, 0 R2, 0 UV
INFO: Process u~ u > mu+ mu- added to mirror process u u~ > mu+ mu-
INFO: Process u~ c > mu+ mu- added to mirror process c u~ > mu+ mu-
INFO: Process c~ u > mu+ mu- added to mirror process u c~ > mu+ mu-
INFO: Process c~ c > mu+ mu- added to mirror process c c~ > mu+ mu-
INFO: Process d~ d > mu+ mu- added to mirror process d d~ > mu+ mu-
INFO: Process s~ s > mu+ mu- added to mirror process s s~ > mu+ mu-
6 processes with 24 diagrams generated in 47.253 s
Total: 6 processes with 24 diagrams
MG5_aMC>output new
INFO: initialize a new directory: new
INFO: remove old information in new
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: u u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Processing color information for loop process: u u~ > h > mu+ mu- [ noborn = NP ]
INFO: Creating color matrix loop process: u u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Generating Helas calls for process: c c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Processing color information for loop process: c c~ > h > mu+ mu- [ noborn = NP ]
INFO: Creating color matrix loop process: c c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Generating Helas calls for process: u c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Processing color information for loop process: u c~ > h > mu+ mu- [ noborn = NP ]
INFO: Creating color matrix loop process: u c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Generating Helas calls for process: c u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Processing color information for loop process: c u~ > h > mu+ mu- [ noborn = NP ]
INFO: Creating color matrix loop process: c u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Generating Helas calls for process: d d~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Processing color information for loop process: d d~ > h > mu+ mu- [ noborn = NP ]
INFO: Creating color matrix loop process: d d~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Generating Helas calls for process: s s~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Processing color information for loop process: s s~ > h > mu+ mu- [ noborn = NP ]
INFO: Creating color matrix loop process: s s~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory /home/ipsita/
ALOHA: aloha creates FFV8 set of routines with options: L2,P0
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory P0_uux_mupmum
INFO: Generating Feynman diagrams for Process: u u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Finding symmetric diagrams for subprocess group uux_mupmum
INFO: Creating files in directory /home/ipsita/
ALOHA: aloha creates FFV8 set of routines with options: L2,P0
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: c c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory P0_ccx_mupmum
INFO: Generating Feynman diagrams for Process: c c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Finding symmetric diagrams for subprocess group ccx_mupmum
INFO: Creating files in directory /home/ipsita/
ALOHA: aloha creates FFV8 set of routines with options: L2,P0
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory P0_ucx_mupmum
INFO: Generating Feynman diagrams for Process: u c~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Finding symmetric diagrams for subprocess group ucx_mupmum
INFO: Creating files in directory /home/ipsita/
ALOHA: aloha creates FFV8 set of routines with options: L2,P0
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: c u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory P0_cux_mupmum
INFO: Generating Feynman diagrams for Process: c u~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Finding symmetric diagrams for subprocess group cux_mupmum
INFO: Creating files in directory /home/ipsita/
ALOHA: aloha creates FFV8 set of routines with options: L2,P0
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d d~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory P0_ddx_mupmum
INFO: Generating Feynman diagrams for Process: d d~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Finding symmetric diagrams for subprocess group ddx_mupmum
INFO: Creating files in directory /home/ipsita/
ALOHA: aloha creates FFV8 set of routines with options: L2,P0
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: s s~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Creating files in directory P0_ssx_mupmum
INFO: Generating Feynman diagrams for Process: s s~ > h > mu+ mu- WEIGHTED<=20 [ noborn = NP ]
INFO: Finding symmetric diagrams for subprocess group ssx_mupmum
Generated helas calls for 6 subprocesses (24 diagrams) in 0.541 s
ALOHA: aloha creates FFV8 set of routines with options: L3
The option gauge is modified [unitary] but will not be written in the configuration files.
If you want to make this value the default for future session, you can run 'save options --all'
save configuration file to /home/ipsita/
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /home/ipsita/
Type "launch" to generate events from this process, or see
/home/ipsita/
Run "open index.html" to see more information about this process.
MG5_aMC>launch new
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.6.4 2018-11-09 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https:/
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /home/ipsita/
INFO: load configuration from /home/ipsita/
INFO: load configuration from /home/ipsita/
Using default text editor "vi". Set another one in ./input/
generate_events run_01
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program detector = Not Avail. |
| 3. Choose an analysis package (plot/convert) analysis = Not Avail. |
| 4. Decay onshell particles madspin = OFF |
| 5. Add weights to events for new hypp. reweight = Not Avail. |
\======
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'madspin=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]
>4
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program detector = Not Avail. |
| 3. Choose an analysis package (plot/convert) analysis = Not Avail. |
| 4. Decay onshell particles madspin = ON |
| 5. Add weights to events for new hypp. reweight = Not Avail. |
\======
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'madspin=OFF' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program detector = Not Avail. |
| 3. Choose an analysis package (plot/convert) analysis = Not Avail. |
| 4. Decay onshell particles madspin = ON |
| 5. Add weights to events for new hypp. reweight = Not Avail. |
\======
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'madspin=OFF' 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
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 |
| 4. MadLoopParams : MadLoopParams.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, 4, enter path, ... ][90s to answer]
>0
INFO: Update the dependent parameter of the param_card.dat
Generating 1000 events with run name run_01
survey run_01
INFO: compile directory
initMadLoop -r -f
Initializing MadLoop loop-induced matrix elements (this can take some time)...
MadLoop initialization finished.
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: P0_uux_mupmum
INFO: P0_ccx_mupmum
INFO: P0_ucx_mupmum
INFO: P0_cux_mupmum
INFO: P0_ddx_mupmum
INFO: P0_ssx_mupmum
INFO: Idle: 3, Running: 4, Completed: 10 [ 0.13s ]
WARNING: Channel 1 has encountered a fraction of 0.00338
of numerically unstable loop matrix element computations
(which could not be rescued using quadruple precision).
The results might not be trusted.
WARNING: Channel 1 has encountered a fraction of 0.00338
of numerically unstable loop matrix element computations
(which could not be rescued using quadruple precision).
The results might not be trusted.
INFO: P0_ccx_mupmum/G1 is at 4e-18 +- 8.42e-20 pb. Now submitting iteration #2.
INFO: P0_uux_mupmum/G1 is at 2.96e-23 +- 6.76e-25 pb. Now submitting iteration #2.
INFO: Idle: 4, Running: 4, Completed: 11 [ 0.18s ]
WARNING: Channel 1 has encountered a fraction of 0.00338
of numerically unstable loop matrix element computations
(which could not be rescued using quadruple precision).
The results might not be trusted.
WARNING: Channel 1 has encountered a fraction of 0.00338
of numerically unstable loop matrix element computations
(which could not be rescued using quadruple precision).
The results might not be trusted.
INFO: P0_ucx_mupmum/G1 is at 7.54e-20 +- 1.81e-21 pb. Now submitting iteration #2.
INFO: P0_cux_mupmum/G1 is at 3.79e-20 +- 8.01e-22 pb. Now submitting iteration #2.
INFO: P0_ddx_mupmum/G1 is at 9.01e-23 +- 1.93e-24 pb. Now submitting iteration #2.
INFO: Idle: 6, Running: 4, Completed: 18 [ 4.8s ]
INFO: Idle: 4, Running: 4, Completed: 23 [ 10.5s ]
INFO: P0_ssx_mupmum/G1 is at 1.07e-20 +- 2.23e-22 pb. Now submitting iteration #2.
INFO: P0_ccx_mupmum/G1 is at 4.052e-18 +- 2.72e-20 pb. Now submitting iteration #3.
INFO: P0_uux_mupmum/G1 is at 2.9569e-23 +- 7.77e-27 pb. Now submitting iteration #3.
INFO: P0_ucx_mupmum/G1 is at 7.463e-20 +- 3.83e-22 pb. Now submitting iteration #3.
INFO: Idle: 6, Running: 4, Completed: 30 [ 13.7s ]
INFO: P0_cux_mupmum/G1 is at 3.837e-20 +- 2.15e-22 pb. Now submitting iteration #3.
INFO: Idle: 5, Running: 4, Completed: 35 [ 22.2s ]
INFO: P0_ddx_mupmum/G1 is at 9.079e-23 +- 3.89e-25 pb. Now submitting iteration #3.
INFO: Idle: 4, Running: 4, Completed: 40 [ 33.8s ]
INFO: P0_ssx_mupmum/G1 is at 1.073e-20 +- 3.43e-23 pb. Now submitting iteration #3.
INFO: Survey finished for P0_ccx_mupmum/G1 at 4.06e-18 +- 1.25e-20 pb
INFO: Survey finished for P0_uux_mupmum/G1 at 2.9584e-23 +- 8.6e-27 pb
INFO: Survey finished for P0_ucx_mupmum/G1 at 7.564e-20 +- 5.58e-22 pb
INFO: Idle: 0, Running: 4, Completed: 48 [ 39.9s ]
INFO: Survey finished for P0_cux_mupmum/G1 at 3.806e-20 +- 1.94e-22 pb
INFO: Survey finished for P0_ssx_mupmum/G1 at 1.072e-20 +- 1.57e-23 pb
INFO: Survey finished for P0_ddx_mupmum/G1 at 9.048e-23 +- 2.41e-25 pb
INFO: Idle: 0, Running: 0, Completed: 54 [ 40s ]
INFO: End survey
refine 1000
Creating Jobs
INFO: Refine results to 1000
INFO: Generating 1000.0 unweigthed events.
INFO: Effective Luminosity 2.74833995797e+20 pb^-1
INFO: need to improve 0 channels
Current estimate of cross-section: 4.184344068e-18 +- 1.24869775523e-20
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 12h30 ]
INFO: Combining runs
INFO: finish refine
refine 1000
Creating Jobs
INFO: Refine results to 1000
INFO: Generating 1000.0 unweigthed events.
INFO: Effective Luminosity 2.74833995797e+20 pb^-1
INFO: need to improve 0 channels
Current estimate of cross-section: 4.184344068e-18 +- 1.24869775523e-20
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 12h30 ]
INFO: Combining runs
INFO: finish refine
INFO: Combining Events
Some PS with numerical instability have been set to a zero matrix-element (0.0446%)
=== Results Summary for run: run_01 tag: tag_1 ===
Cumulative sequential time for this run: 3m01s
Cross-section : 4.184e-18 +- 1.249e-20 pb
Nb of events : 1000
INFO: No version of lhapdf. Can not run systematics computation
store_events
INFO: Storing parton level results
INFO: End Parton
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 > h > mu+ mu-
INFO: options:
INFO: detected model: /home/ipsita/
Set group_subprocesses to Auto
Note that you need to regenerate all processes
save options auto_update
save configuration file to /home/ipsita/
save options cluster_type
save configuration file to /home/ipsita/
save options OLP
save configuration file to /home/ipsita/
save options applgrid
save configuration file to /home/ipsita/
lhapdf-config 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/
save options lhapdf
save configuration file to /home/ipsita/
set ninja to /home/ipsita/
save options ninja
save configuration file to /home/ipsita/
save options amcfast
save configuration file to /home/ipsita/
save options output_dependencies
save configuration file to /home/ipsita/
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 27434403
INFO: We need to recalculate the branching fractions for t~,w-,z,w+,t
INFO: using the FeynRules formula present in the model (arXiv:1402.1178)
INFO: load particles
INFO: load vertices
INFO:
INFO: decay channels for t~ : ( width = 1.466877 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b~ w-
INFO:
INFO:
INFO: decay channels for w- : ( width = 2.004505 GeV )
INFO: BR d1 d2
INFO: 3.163626e-01 d u~
INFO: 3.162417e-01 s c~
INFO: 1.110247e-01 e- ve~
INFO: 1.110244e-01 mu- vm~
INFO: 1.109421e-01 ta- vt~
INFO: 1.692634e-02 s u~
INFO: 1.691995e-02 d c~
INFO: 5.543055e-04 b c~
INFO: 3.880363e-06 b u~
INFO:
INFO:
INFO: decay channels for z : ( width = 2.411596 GeV )
INFO: BR d1 d2
INFO: 1.520727e-01 d~ d
INFO: 1.520719e-01 s~ s
INFO: 1.504304e-01 b~ b
INFO: 1.178265e-01 u~ u
INFO: 1.177064e-01 c~ c
INFO: 6.878490e-02 vt~ vt
INFO: 6.878490e-02 vm~ vm
INFO: 6.878490e-02 ve~ ve
INFO: 3.453867e-02 e+ e-
INFO: 3.453839e-02 mu+ mu-
INFO: 3.446033e-02 ta+ ta-
INFO:
INFO:
INFO: decay channels for w+ : ( width = 2.004505 GeV )
INFO: BR d1 d2
INFO: 3.163626e-01 d~ u
INFO: 3.162417e-01 s~ c
INFO: 1.110247e-01 e+ ve
INFO: 1.110244e-01 mu+ vm
INFO: 1.109421e-01 ta+ vt
INFO: 1.692634e-02 s~ u
INFO: 1.691995e-02 d~ c
INFO: 5.543055e-04 b~ c
INFO: 3.880363e-06 b~ u
INFO:
INFO:
INFO: decay channels for t : ( width = 1.466877 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b w+
INFO:
INFO: generating the production square matrix element
INFO: generate p p > h > mu+ mu- [noborn = NP] ;
Error detected in "generate_events run_01"
write debug file /home/ipsita/
If you need help with this issue please contact us on https:/
MadGraph5Error : The analytic information wavefunction_rank has not been computed yet for this wavefunction and an alohaModel was not specified, so that the information cannot be retrieved.
quit
INFO: storing files of previous run
gzipping output file: unweighted_
INFO: Done
INFO:
The debug file is attached.
Thanks in advance
| IPSITA RAY (ray176121005) wrote | #1 |
| Olivier Mattelaer (olivier-mattelaer) wrote | #2 |
| IPSITA RAY (ray176121005) wrote | #3 |
| Olivier Mattelaer (olivier-mattelaer) wrote Re: [Bug 1816293] Re: error detected in generating events | #4 |
| Changed in mg5amcnlo: | |
| status: | New → Invalid |
| IPSITA RAY (ray176121005) wrote | #5 |
| Olivier Mattelaer (olivier-mattelaer) wrote | #6 |
Hi,
Using MadSpin here does not make any sense.
What particle are you trying to decay?
In top of that MadSpin is not able to be used for loop-induced process.
(you have to pass to spinmode=none for madspin to work)
Cheers,
Olivier