problem with decay

Asked by Zhaowei Song

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/sm/restrict_default.dat .
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/loop_sm/restrict_default.dat .
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/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh/SubProcesses/PV0_0_1_gg_h
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/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh done.
Type "launch" to generate events from this process, or see
/home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh/README
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://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh/Cards/me5_configuration.txt
INFO: load configuration from /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/input/mg5_configuration.txt
INFO: load configuration from /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
No valid rivet path found
generate_events run_01
The following switches determine which programs are run:
/==============================================================================\
| 1. Choose the shower/hadronization program shower = OFF |
| 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/hadronization program shower = Pythia8 |
| 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/hadronization program shower = Pythia8 |
| 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/hadronization program shower = Pythia8 |
| 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_parton_card.dat |
| 6. madanalysis5_hadron : madanalysis5_hadron_card.dat |
| 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/MadWidth/...).
     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_parton_card.dat |
| 6. madanalysis5_hadron : madanalysis5_hadron_card.dat |
| 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/MadWidth/...).
     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/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/HEPTools/madanalysis5/madanalysis5'.
The specified version of MadAnalysis5 will not be active in your session.
set fastjet to fastjet-config
set ninja to /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/HEPTools/lib
set collier to /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/HEPTools/lib
/home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/HEPTools/lhapdf6_py3/bin/lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

/home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/HEPTools/lhapdf6_py3/bin/lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

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/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/models/loop_sm with file models/loop_sm/restrict_default.dat .
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_events.lhe
INFO: Done
INFO:
more information in /home/parallels/Downloads/MG5_aMC_v3.4.0.RC2/MG5_aMC_v3_4_0_RC2/ggh/index.html

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Zhaowei Song
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

You are here facing a couple of the limitation of the main MadSpin algorithm.
1) that algorithm can not handle to decay for loop-induced matrix. (since MG5aMC is not able to generate such matrix-element)
2) that algorithm can only handle sequence of two body-decay (and your syntax is a four body decay)

You need to use another algorithm:
in this case have the line in madspin
set spinmode none

This mode is perfect for Higgs production but very bad for any other particle.
It does not have any limitation but you do not have spin-correlation between production and decay, nor off-shell effect (none of those are relevant for Higgs).

Please look at 1507.00020 for all details concerning loop-induced processes.

Revision history for this message
Zhaowei Song (zhaowei12345678) said :
#2

Thank you Olivier,that solved my problem