problem with heft
I was running g g > h > h h and decay h>zz,z>ll
MG5_aMC>import model heft
INFO: Restrict model heft 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 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 h1 w- ta- ta+
MG5_aMC>generate g g > h > h h , h > z z , z > l+ l-
WARNING: Decay without corresponding particle in core process found.
Decay information for particle(s) z is discarded.
Please check your process definition carefully.
This warning usually means that you forgot parentheses in presence of subdecay.
Example of correct syntax: p p > t t~, ( t > w+ b, w+ > l+ vl)
1 processes with 2 diagrams generated in 0.014 s
Total: 1 processes with 2 diagrams
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
|
Revision history for this message
|
#1 |
Hi,
So as indicated in the warning message, you did forget dome parentheses:
You need to say that the Z are the one coming from the Higgs.
so I guess that what you try to write was
generate g g > h > h h , (h > z z , z > l+ l-)
Now I do see two issues with that syntax:
1) on the " h > z z " this means that you will ask one onshell Higgs decaying into two Z boson.
This is not possible and will create issue, you should rather use
h > z > 2l+ 2l-
2) on the g g > h > h h
Let me apologize, if my following comment here is trivial for you. You might know this already and have plan method to tackle that correctly. In that case just ignore that comment.
You are using the heft model, this method assumes that the top mass is much heavier than the Higgs mass and of any other observable of the system (in particular sqrt(shat) ). While this is kind of working for single Higgs production, the approximation is quite bad for double Higgs production.
So in general for double Higgs production, I would advise to either reweight your sample by the True one loop computation or do directly the loop-induced generation (likely more cpu intensive):
You can use the following script to get the final state that you want.
import model loop_sm
generate g g > h h QED<=99 [noborn=QCD]
output
launch
MadSpin=onshell
decay h > z > 2l+ 2l-
A second comment on your syntax is that you do only include the equivalent of the top triangle and do not include the box (my loop-induced syntax includes both). The two diagrams have strong destructive interference, and therefore it is quite surprising to me that you separate them. Obviously, this is also something that you can correct later and again I apologize if you have a strategy ready to handle such interference in one way or another way (or if your goal is to study such interference).
Cheers,
Olivier
|
Revision history for this message
|
#2 |
I just run as you suggest and there is another error MG5_aMC>generate g g > h h QED<=99 [noborn=QCD]
1 processes with 16 diagrams generated in 0.106 s
Total: 1 processes with 16 diagrams
MG5_aMC>output ggh_loop
INFO: initialize a new directory: ggh_loop
INFO: remove old information in ggh_loop
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Creating files in directory /home/parallels
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_hh
INFO: Generating Feynman diagrams for Process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_hh
Generated helas calls for 1 subprocesses (16 diagrams) in 0.032 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 22 routines in 0.138 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.RC 2022-03-18 *
* *
* 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/
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 = MadAnalysis4 |
| 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 = MadAnalysis4 |
| 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 = MadAnalysis4 |
| 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. plot : plot_card.dat |
| 6. 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 > 2l+ 2l-
INFO: change madspin_card to define the decay of h: h > z > 2l+ 2l-
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. plot : plot_card.dat |
| 6. 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
WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...)
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
initMadLoop -r -f
Error detected in "generate_events run_01"
write debug file /home/parallels
If you need help with this issue please contact us on https:/
str : A compilation Error occurs when trying to compile /home/parallels
The compilation fails with the following output message:
cd MODEL; make
make[1]: Entering directory '/home/
gfortran -w -fPIC -O -ffixed-
..
54 | DOUBLE PRECISION s_qalps(
| 1
Error: Symbol ‘max_particles’ at (1) has no IMPLICIT type
make[1]: *** [<builtin>: couplings.o] Error 1
make[1]: Leaving directory '/home/
make: *** [makefile:53: ../lib/libmodel.a] Error 2
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
quit
INFO:
|
Revision history for this message
|
#3 |
I just run as you suggest and there is another error MG5_aMC>generate g g > h h QED<=99 [noborn=QCD]
1 processes with 16 diagrams generated in 0.106 s
Total: 1 processes with 16 diagrams
MG5_aMC>output ggh_loop
INFO: initialize a new directory: ggh_loop
INFO: remove old information in ggh_loop
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Creating files in directory /home/parallels
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_hh
INFO: Generating Feynman diagrams for Process: g g > h h QED<=99 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_hh
Generated helas calls for 1 subprocesses (16 diagrams) in 0.032 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 22 routines in 0.138 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.RC 2022-03-18 *
* *
* 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/
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 = MadAnalysis4 |
| 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 = MadAnalysis4 |
| 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 = MadAnalysis4 |
| 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. plot : plot_card.dat |
| 6. 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 > 2l+ 2l-
INFO: change madspin_card to define the decay of h: h > z > 2l+ 2l-
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. plot : plot_card.dat |
| 6. 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
WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...)
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
initMadLoop -r -f
Error detected in "generate_events run_01"
write debug file /home/parallels
If you need help with this issue please contact us on https:/
str : A compilation Error occurs when trying to compile /home/parallels
The compilation fails with the following output message:
cd MODEL; make
make[1]: Entering directory '/home/
gfortran -w -fPIC -O -ffixed-
..
54 | DOUBLE PRECISION s_qalps(
| 1
Error: Symbol ‘max_particles’ at (1) has no IMPLICIT type
make[1]: *** [<builtin>: couplings.o] Error 1
make[1]: Leaving directory '/home/
make: *** [makefile:53: ../lib/libmodel.a] Error 2
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
quit
INFO:
|
Revision history for this message
|
#5 |
Thank you Olivier
Actually I'm only a beginner in this field so problem you mentioned I didn't take into account.
What I am trying to do is to generate trilinear Higgs coupling through gluon fusion, I tried to do this through the box.Would you mind telling me how to do it if it is possible.
Besides,would you mind telling me more detail about why heft approximation is quite bad for double Higgs production. I see some people calculating double Higgs production through heft.
Regard
Zhaowei
|
Revision history for this message
|
#6 |
That issue has been fixed in 3.4.0.RC2.
Thanks for reporting a bug on our release candidate.
The syntax that I propose includes both the box and the triangle (you should see both diagrams)
>Besides,would you mind telling me more detail about why heft approximation is quite bad for double Higgs production.
I can not say much more than what I said in my above post, the heft is contracting the loop by using an effective field theory where the High scale is the top mass, so the theory is valid up to the top mass. You can check the literature for more detailed theory explanation and for the phenomenological comparison between the two theories. (see for example: 1401.7340 when di-higgs is done approximatively at NLO. Note that REAL NLO has been computed since then, even if not available in MG5aMC)
|
Revision history for this message
|
#7 |
Dear Olivier
Thank you for your help.
I just upgraded my mg5 to RC2 and run as you suggested but there is still an error.MG5_
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+
MG5_aMC>generate g g > h > h h QED<=99 [noborn=QCD]
INFO: Contributing diagrams generated: 0 Born, 4 loops, 2 R2, 0 UV
1 processes with 4 diagrams generated in 0.102 s
Total: 1 processes with 4 diagrams
MG5_aMC>output ggh_loop
INFO: directory /home/parallels
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: ggh_loop
INFO: remove old information in ggh_loop
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h > h h QED<=99 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h > h h QED<=99 [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h > h h QED<=99 [ noborn = QCD ]
INFO: Creating files in directory /home/parallels
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h > h h QED<=99 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_hh
INFO: Generating Feynman diagrams for Process: g g > h > h h QED<=99 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_hh
Generated helas calls for 1 subprocesses (4 diagrams) in 0.009 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 19 routines in 0.118 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/
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.
>madspin=onshell
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 = onshell |
| 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.
>
INFO: writting at line 22 of the file madspin_card.dat the line: "set spinmode onshell"
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 > 2l+ 2l-
INFO: change madspin_card to define the decay of h: h > z > 2l+ 2l-
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
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_hh
INFO: Idle: 1, Running: 1, Completed: 0 [ current time: 07h32 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 0.22s ]
INFO: P0_gg_hh/G1 is at 0.00409 +- 9.72e-05 pb. Now submitting iteration #2.
INFO: P0_gg_hh/G1 is at 0.004152 +- 2.25e-05 pb. Now submitting iteration #3.
INFO: Survey finished for P0_gg_hh/G1 at 0.004144 +- 7.71e-06 pb
INFO: Idle: 0, Running: 0, Completed: 9 [ 1.3s ]
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 2775163.4933275417 pb^-1
INFO: need to improve 1 channels
- Current estimate of cross-section: 0.0041439 +- 7.7073e-06
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 07h32 ]
INFO: P0_gg_hh/G1 is at 3535/11395 (31%) event. Resubmit 2 job at iteration 1.
INFO: P0_gg_hh/G1 is at 7142/11460 (62%) event. Resubmit 3 job at iteration 2.
INFO: Idle: 1, Running: 0, Completed: 9 [ 2.1s ]
INFO: found enough event for P0_gg_hh/G1
INFO: Idle: 0, Running: 0, Completed: 10 [ 2.1s ]
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: 5 seconds
Cross-section : 0.004122 +- 1.358e-05 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 > h h QED<=99
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
set spinmode onshell
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 > 2l+ 2l-
launch
INFO: Will use seed 650235051
INFO: generate 82000 decay event for particle h
INFO: generating the production square matrix element for onshell
INFO: generating the full matrix element squared (with decay)
INFO: generate matrix element for decay only (1 - > N).
Command "generate_events run_01" interrupted with error:
InvalidCmd : Wrong use of ">" special character.
quit
INFO: storing files of previous run
gzipping output file: unweighted_
INFO: Done
INFO:
Regard
Zhaowei
|
Revision history for this message
|
#8 |
Hi,
Yes indeed madspin=onshell is not supported for loop induced mode (sorry for such mistake)
you need to use the madspin=none mode which is working here.
but looks like you realy need to edit the madspin_card by hand because of a bug in the edition of that file (that I'm going to fix)
Can you help with this problem?
Provide an answer of your own, or ask Zhaowei Song for more information if necessary.
