problem with decay
I was generating g g > h with decay h > z z > l l l l and something is wrong with decay .MG5_aMC>import model sm
INFO: Restrict model sm with file models/
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Pass the definition of 'j' and 'p' to 4 flavour scheme.
Kept definitions of multiparticles l+ / l- / vl / vl~ unchanged
Removed obsolete multiparticles all
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
MG5_aMC>generate g g > h [noborn=QCD]
The current model sm does not allow to generate loop corrections of type ['QCD'].
MG5_aMC now loads 'loop_sm'.
import model loop_sm
INFO: Restrict model loop_sm with file models/
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ b t b~ t~ z w+ h w- ta- ta+
1 processes with 4 diagrams generated in 0.020 s
Total: 1 processes with 4 diagrams
MG5_aMC>output ggh
INFO: initialize a new directory: ggh
INFO: remove old information in ggh
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Creating files in directory /home/parallels
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_h
INFO: Generating Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_h
Generated helas calls for 1 subprocesses (4 diagrams) in 0.008 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 7 routines in 0.088 s
save configuration file to /home/parallels
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /home/parallels
Type "launch" to generate events from this process, or see
/home/parallels
Run "open index.html" to see more information about this process.
MG5_aMC>launch
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.4.0.RC2 2022-04-16 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https:/
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /home/parallels
INFO: load configuration from /home/parallels
INFO: load configuration from /home/parallels
Using default text editor "vi". Set another one in ./input/
No valid rivet path found
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 = OFF |
| 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 'shower=Pythia8' 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]
>1
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 = OFF |
| 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 'shower=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.
>3
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 = MadAnalysis5 |
| 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 'shower=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.
>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 = MadAnalysis5 |
| 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 'shower=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.
>
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
| 4. madspin : madspin_card.dat |
| 5. madanalysis5_parton : madanalysis5_
| 6. madanalysis5_hadron : madanalysis5_
| 7. 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, pythia8, 4, enter path, ... ][90s to answer]
>decay h > z z > l- l- l+ l+
INFO: change madspin_card to define the decay of h: h > z z > l- l- l+ l+
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
| 4. madspin : madspin_card.dat |
| 5. madanalysis5_parton : madanalysis5_
| 6. madanalysis5_hadron : madanalysis5_
| 7. 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, pythia8, 4, enter path, ... ]
>
INFO: Update the dependent parameter of the param_card.dat
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
initMadLoop -r -f
ATTENTION: default value of option mesa_glthread overridden by environment.
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_gg_h
INFO: Idle: 2, Running: 0, Completed: 0 [ current time: 14h10 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 0.22s ]
INFO: P0_gg_h/G1 is at 16 +- 0.267 pb. Now submitting iteration #2.
INFO: P0_gg_h/G1 is at 15.85 +- 0.05 pb. Now submitting iteration #3.
INFO: Survey finished for P0_gg_h/G1 at 15.72 +- 0.0483 pb
INFO: Idle: 0, Running: 0, Completed: 9 [ 1.2s ]
INFO: End survey
refine 10000
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
sum of cpu time of last step: 0 second
INFO: Effective Luminosity 731.5056294128872 pb^-1
INFO: need to improve 1 channels
- Current estimate of cross-section: 15.721 +- 0.048287
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 14h11 ]
INFO: P0_gg_h/G1 is at 3884/11470 (34%) event. Resubmit 2 job at iteration 1.
INFO: P0_gg_h/G1 is at 7755/11465 (68%) event. Resubmit 2 job at iteration 2.
INFO: found enough event for P0_gg_h/G1
INFO: Idle: 0, Running: 0, Completed: 9 [ 1.3s ]
INFO: Combining runs
sum of cpu time of last step: 0 second
INFO: finish refine
refine 10000 --treshold=0.9
No need for second refine due to stability of cross-section
INFO: Combining Events
=== Results Summary for run: run_01 tag: tag_1 ===
Cumulative sequential time for this run: 4 seconds
Cross-section : 15.67 +- 0.003413 pb
Nb of events : 10000
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: g g > h
INFO: options:
INFO: detected model: loop_sm. Loading...
No MadAnalysis5 version number could be read from the path supplied '/home/
The specified version of MadAnalysis5 will not be active in your session.
set fastjet to fastjet-config
set ninja to /home/parallels
set collier to /home/parallels
/home/parallels
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/
/home/parallels
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/
Set group_subprocesses to Auto
Note that you need to regenerate all processes
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
decay h > z z > l- l- l+ l+
launch
INFO: Will use seed 724269077
INFO: We need to recalculate the branching fractions for w+,t,z,h,t~,w-
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /home/parallels
INFO: Run "set stdout_level DEBUG" before import for more information.
WARNING: The LO estimate for the width of particle 25
WARNING: differs from the one in the banner by 12 percent
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 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 h : ( width = 0.00564955 GeV )
INFO: BR d1 d2
INFO: 9.541999e-01 b~ b
INFO: 4.580011e-02 ta+ ta-
INFO:
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: generating the production square matrix element
INFO: generate g g > h [sqrvirt=QCD] ;
Command "generate_events run_01" interrupted with error:
InvalidCmd : not supported format standalone_msP
quit
INFO: storing files of previous run
gzipping output file: unweighted_
INFO: Done
INFO:
more information in /home/parallels
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Zhaowei Song
- Solved:
- Last query:
- Last reply: