error in generating decay channel h1 > a a from (p p > h1 xd xd~ [noborn=QCD])in Pseudoscalar_2HDM

nice answer

Did you check the paper where we introduce the loop-induced mode?
We did use an example very close to yours which has the same limitation and we explain how you can do it with madspin (and the spinmode=none mode).
For the decay of the Higgs, you can either use another model where an effective vertex is used (which is the most efficient way of doing it) or asked for a loop induced decay within your model (but this is quite slow)

Cheers,

Olivier

Did you change the model to use within madspin?
you need to set the model to heft (and be sure to use the same higgs mass obviously).

Cheers,

Olivier

Hi Olivier
i have import model heft but there is no (xd xd~) particle
MG5_aMC>generate p p > h, h > a a
1 processes with 2 diagrams generated in 0.049 s
Total: 1 processes with 2 diagrams
MG5_aMC>import model heft
INFO: Restrict model heft with file models/heft/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 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 p p > h, h > a a
1 processes with 2 diagrams generated in 0.027 s

but the problem is I need (p p > h1 xd xd~) where h1 > a a

Also I saw a paper as u say earlier and I try this

MG5_aMC>import model Pseudoscalar_2HDM UFO
INFO: Restrict model Pseudoscalar_2HDM with file models/Pseudoscalar_2HDM/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
Removed obsolete multiparticles all
Defined multiparticle all = g ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h+ h1 h2 h3 h4 w- ghz~ ghwp~ ghwm~ h- ta- xd ta+ xd~
MG5_aMC>generate p p > h1 xd xd~ , h1 > a a
Command "generate p p > h1 xd xd~ , h1 > a a" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: h1 > a a WEIGHTED=4. Please enter a valid process
MG5_aMC>generate p p > h1 xd xd~ , h1 > a a [noborn=QCD ]
Command "generate p p > h1 xd xd~ , h1 > a a [noborn=QCD ]" interrupted with error:
InvalidCmd : No particle , in model
MG5_aMC>import model Pseudoscalar_2HDM NLO
INFO: Restrict model Pseudoscalar_2HDM with file models/Pseudoscalar_2HDM/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 ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h+ h1 h2 h3 h4 w- ghz~ ghwp~ ghwm~ h- ta- xd ta+ xd~
MG5_aMC>generate p p > h1 xd xd~ , h1 > a a [noborn=QCD ]
Command "generate p p > h1 xd xd~ , h1 > a a [noborn=QCD ]" interrupted with error:
InvalidCmd : No particle , in model
MG5_aMC>

also I change in madspin in
decay h1 > a a
 but h1 > a a doesn't decay

Hi,

In this case you need two different model, one for the loop-induced production and one for the decay
so the production should be
import model Pseudoscalar_2HDM UFO
generate p p > h1 xd xd~ [noborn=QCD ]
output
launch
madspin=ON

and the madspin_card should be
import model heft /path/to/valid_param_card_for_heft.dat
set spinmode none
decay h1 > a a

where you have to replace /path/to/valid_param_card_for_heft.dat as needed (and be sure that your h1 has the same pdg/mass between the two model.

Here is the documentation for the import model command in MG5aMC:
import model ModelNane ParamCardPath [--bypass_check]
Allows to use another model/param_card than the one used for the generation of the events. By default it is checked that the old/new param_card are compatible (and if not the code stops), except if the optional argument --bypass_check is included.

Hi Olivier
i have imported model heft with Pseudoscalar_2HDM UFO using madspin ON
but error same .
NoDiagramException : No amplitudes generated from process Process: h1 > a a. Please enter a valid process

MG5_aMC>import model Pseudoscalar_2HDM UFO
INFO: Restrict model Pseudoscalar_2HDM with file models/Pseudoscalar_2HDM/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 ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h+ h1 h2 h3 h4 w- ghz~ ghwp~ ghwm~ h- ta- xd ta+ xd~
MG5_aMC>generate p p > h1 xd xd~ [noborn=QCD ]
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 QCD+2*QED+DS
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS QCD+2*QED+DS
INFO: Contributing diagrams generated: 0 Born, 60(+4) loops, 0 R2, 0 UV
INFO: Contributing diagrams generated: 0 Born, 2 loops, 2 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 2 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 2 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 2 R2, 0 UV
INFO: Process u~ u > h1 xd xd~ added to mirror process u u~ > h1 xd xd~
INFO: Process c~ c > h1 xd xd~ added to mirror process c c~ > h1 xd xd~
INFO: Process d~ d > h1 xd xd~ added to mirror process d d~ > h1 xd xd~
INFO: Process s~ s > h1 xd xd~ added to mirror process s s~ > h1 xd xd~
5 processes with 68 diagrams generated in 12.817 s
Total: 5 processes with 68 diagrams
MG5_aMC>output final_decay_m
First output using loop matrix-elements has been detected. Now asking for loop reduction:
For loop computations, MadLoop requires dedicated tools to perform the reduction of loop Feynman diagrams using OPP-based and/or TIR approaches.

Which one do you want to install? (this needs to be done only once)
1. cuttools (OPP) [0711.3596] : will be installed (required)
2. iregi (TIR) [1405.0301] : will be installed (required)
3. ninja (OPP) [1403.1229] : /home/faisal/Downloads/MG5_aMC_v3_4_1/ninja/lib
4. collier (TIR) [1604.06792] : /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/collier
5. golem (TIR) [0807.0605] : do not install
You can:
 -> hit 'enter' to proceed
 -> type a number to cycle its options
 -> enter the following command:
    {tool_name} [install|noinstall|{prefixed_installation_path}]

If you are unsure about what this question means, just type enter to proceed. [300s to answer]
>
set ninja /home/faisal/Downloads/MG5_aMC_v3_4_1/ninja/lib
/home/faisal/Downloads/MG5_aMC_v3_4_1/ninja/lib does not seem to correspond to a valid ninja lib . Please enter the full PATH/TO/ninja/lib .
You will NOT be able to run ninja otherwise.

save options ninja
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
save options ninja
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
set collier /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/collier
set collier to /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/collier
save options collier
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
save options collier
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
set golem ''
save options golem
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
save options golem
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
INFO: initialize a new directory: final_decay_m
INFO: remove old information in final_decay_m
WARNING: The ninja reduction library could not be found with PATH:/home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/ninja/bin specified in mg5_configuration.txt. It will not be available.
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Helas calls for process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: u u~ > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Helas calls for process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: c c~ > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Helas calls for process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: d d~ > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Combined process s s~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ] with process d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/SubProcesses/PV0_0_1_gg_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_h1xdxdx
INFO: Generating Feynman diagrams for Process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_h1xdxdx
INFO: Creating files in directory /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/SubProcesses/PV0_1_1_uux_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/SubProcesses/PV0_1_2_ccx_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/SubProcesses/PV0_1_3_ddx_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory P0_qq_h1xdxdx
INFO: Generating Feynman diagrams for Process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Feynman diagrams for Process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Feynman diagrams for Process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group qq_h1xdxdx
Generated helas calls for 4 subprocesses (66 diagrams) in 0.506 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 35 routines in 1.853 s
save configuration file to /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/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/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m done.
Type "launch" to generate events from this process, or see
/home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/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.2 2023-01-20 *
* *
* 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/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/Cards/me5_configuration.txt
INFO: load configuration from /home/faisal/Downloads/MG5_aMC_v3_4_1/input/mg5_configuration.txt
INFO: load configuration from /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/Cards/me5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/=================== Description ===================|============= values ==============|======== other options ========\
| 1. Choose the shower/hadronization program | shower = Not Avail. | Please install module |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = Not Avail. | Please install module |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=======================================================================================================================/
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:
/=================== Description ===================|============= values ==============|======== other options ========\
| 1. Choose the shower/hadronization program | shower = Not Avail. | Please install module |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = Not Avail. | Please install module |
| 4. Decay onshell particles | madspin = ON | onshell|full|OFF |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=======================================================================================================================/
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/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ][90s to answer]
>3
open /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/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 |
| 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/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ]
>0
INFO: Update the dependent parameter of the param_card.dat
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.13
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
initMadLoop -r -f
Gtk-Message: 16:15:16.243: Not loading module "atk-bridge": The functionality is provided by GTK natively. Please try to not load it.
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_h1xdxdx
INFO: P0_qq_h1xdxdx
INFO: Idle: 28, Running: 8, Completed: 0 [ current time: 16h16 ]
INFO: Idle: 27, Running: 8, Completed: 1 [ 8.7s ]
INFO: Idle: 25, Running: 8, Completed: 5 [ 14.1s ]
INFO: Idle: 24, Running: 8, Completed: 8 [ 17.7s ]
INFO: Idle: 23, Running: 8, Completed: 9 [ 26.5s ]
INFO: Idle: 21, Running: 8, Completed: 13 [ 34.3s ]
INFO: Idle: 20, Running: 8, Completed: 14 [ 37.5s ]
INFO: Idle: 19, Running: 8, Completed: 17 [ 44.5s ]
INFO: Idle: 18, Running: 8, Completed: 19 [ 51.8s ]
INFO: Idle: 17, Running: 8, Completed: 21 [ 55.4s ]
INFO: Idle: 15, Running: 8, Completed: 25 [ 58.7s ]
INFO: P0_gg_h1xdxdx/G2.7 is at 0.00142 +- 0.000125 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G1.7 is at 0.0198 +- 0.00136 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G2.8 is at 8.33e-05 +- 3.29e-05 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G1.8 is at 0.000252 +- 0.000182 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G3.5 is at 0.192 +- 0.0161 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G3.7 is at 0.144 +- 0.01 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G3.8 is at 0.028 +- 0.0189 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G4.5 is at 0.000854 +- 3.06e-05 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G4.7 is at 0.0566 +- 0.00387 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G5 is at 245 +- 4.74 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G6.5 is at 0.0524 +- 0.00214 pb. Now submitting iteration #2.
INFO: P0_qq_h1xdxdx/G1.7 is at 0.000422 +- 3.02e-05 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G4.8 is at 0.0107 +- 0.00903 pb. Now submitting iteration #2.
INFO: P0_qq_h1xdxdx/G1.8 is at 2.24e-06 +- 1.1e-06 pb. Now submitting iteration #2.
INFO: Idle: 25, Running: 8, Completed: 47 [ 1m 2s ]
INFO: Idle: 24, Running: 8, Completed: 49 [ 1m 18s ]
INFO: Idle: 23, Running: 8, Completed: 51 [ 1m 25s ]
INFO: P0_qq_h1xdxdx/G2.7 is at 0.00152 +- 0.000124 pb. Now submitting iteration #2.
INFO: Idle: 22, Running: 8, Completed: 57 [ 1m 31s ]
INFO: Idle: 21, Running: 8, Completed: 58 [ 1m 44s ]
INFO: Idle: 20, Running: 8, Completed: 60 [ 1m 49s ]
INFO: Idle: 18, Running: 8, Completed: 64 [ 1m 54s ]
INFO: Idle: 18, Running: 8, Completed: 65 [ 1m 57s ]
INFO: Idle: 17, Running: 8, Completed: 66 [ 2m 9s ]
INFO: Idle: 15, Running: 8, Completed: 69 [ 2m 14s ]
INFO: P0_qq_h1xdxdx/G2.8 is at 0.000137 +- 7.04e-05 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G6.7 is at 0.996 +- 0.083 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G6.8 is at 0.0608 +- 0.0365 pb. Now submitting iteration #2.
INFO: P0_gg_h1xdxdx/G2.7 is at 0.00176 +- 9.89e-05 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G1.7 is at 0.021 +- 0.00037 pb. Now submitting iteration #3.
INFO: Survey finished for P0_gg_h1xdxdx/G2.8 at 0.000209 +- 0.000111 pb
INFO: P0_gg_h1xdxdx/G1.8 is at 0.00549 +- 0.00697 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G3.5 is at 0.224 +- 0.00764 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G3.7 is at 0.1481 +- 0.00123 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G3.8 is at 0.0145 +- 0.00438 pb. Now submitting iteration #3.
INFO: P0_qq_h1xdxdx/G1.7 is at 0.000454 +- 1.04e-05 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G4.7 is at 0.0596 +- 0.000905 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G5 is at 246.1 +- 0.294 pb. Now submitting iteration #3.
INFO: Idle: 24, Running: 8, Completed: 88 [ 2m 17s ]
INFO: Idle: 21, Running: 8, Completed: 93 [ 2m 21s ]
INFO: Survey finished for P0_qq_h1xdxdx/G1.8 at 6.92e-06 +- 2.97e-06 pb
INFO: Idle: 20, Running: 8, Completed: 95 [ 2m 27s ]
INFO: Idle: 19, Running: 8, Completed: 96 [ 2m 40s ]
INFO: Idle: 18, Running: 8, Completed: 98 [ 2m 43s ]
INFO: P0_gg_h1xdxdx/G4.5 is at 0.000831 +- 9.25e-06 pb. Now submitting iteration #3.
INFO: Idle: 17, Running: 8, Completed: 103 [ 3m 11s ]
INFO: Idle: 17, Running: 8, Completed: 104 [ 3m 18s ]
INFO: Idle: 16, Running: 8, Completed: 105 [ 3m 30s ]
INFO: Idle: 14, Running: 8, Completed: 109 [ 3m 34s ]
INFO: Idle: 12, Running: 8, Completed: 113 [ 4m 3s ]
INFO: Survey finished for P0_qq_h1xdxdx/G2.8 at 0.000116 +- 7.1e-06 pb
INFO: P0_qq_h1xdxdx/G2.7 is at 0.00177 +- 7.88e-05 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G6.5 is at 0.0507 +- 0.000584 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G4.8 is at 0.00164 +- 0.00147 pb. Now submitting iteration #3.
INFO: Idle: 14, Running: 8, Completed: 118 [ 4m 7s ]
INFO: Idle: 13, Running: 8, Completed: 119 [ 4m 18s ]
INFO: P0_gg_h1xdxdx/G6.7 is at 1.01 +- 0.00454 pb. Now submitting iteration #3.
INFO: P0_gg_h1xdxdx/G6.8 is at 0.0788 +- 0.00416 pb. Now submitting iteration #3.
INFO: Survey finished for P0_gg_h1xdxdx/G2.7 at 0.00175 +- 2.78e-05 pb
INFO: Survey finished for P0_gg_h1xdxdx/G1.7 at 0.02139 +- 0.000169 pb
INFO: P0_gg_h1xdxdx/G1.8 is at 0.00273 +- 0.00292 pb. Now submitting iteration #4.
INFO: P0_gg_h1xdxdx/G3.5 is at 0.228 +- 0.00594 pb. Now submitting iteration #4.
INFO: Survey finished for P0_qq_h1xdxdx/G1.7 at 0.000444 +- 4.6e-06 pb
INFO: Idle: 13, Running: 8, Completed: 128 [ 4m 23s ]
INFO: Survey finished for P0_gg_h1xdxdx/G3.7 at 0.1445 +- 0.00122 pb
INFO: Idle: 10, Running: 8, Completed: 134 [ 4m 27s ]
INFO: Idle: 9, Running: 8, Completed: 135 [ 5m 3s ]
INFO: Idle: 9, Running: 8, Completed: 136 [ 5m 8s ]
INFO: Idle: 8, Running: 8, Completed: 137 [ 5m 16s ]
INFO: Idle: 7, Running: 8, Completed: 140 [ 5m 20s ]
INFO: Idle: 6, Running: 8, Completed: 141 [ 5m 25s ]
INFO: P0_gg_h1xdxdx/G3.8 is at 0.0214 +- 0.00529 pb. Now submitting iteration #4.
INFO: Survey finished for P0_qq_h1xdxdx/G2.7 at 0.00183 +- 3.07e-05 pb
INFO: Survey finished for P0_gg_h1xdxdx/G4.7 at 0.05832 +- 0.000497 pb
INFO: Survey finished for P0_gg_h1xdxdx/G5 at 247.8 +- 0.608 pb
INFO: Survey finished for P0_gg_h1xdxdx/G4.5 at 0.0008408 +- 5.16e-06 pb
INFO: P0_gg_h1xdxdx/G4.8 is at 0.00333 +- 0.00119 pb. Now submitting iteration #4.
INFO: Survey finished for P0_gg_h1xdxdx/G6.7 at 0.9924 +- 0.00563 pb
INFO: Idle: 2, Running: 8, Completed: 150 [ 5m 51s ]
INFO: Idle: 2, Running: 8, Completed: 151 [ 5m 56s ]
INFO: Idle: 1, Running: 8, Completed: 152 [ 6m 51s ]
INFO: Survey finished for P0_gg_h1xdxdx/G6.5 at 0.05069 +- 0.000196 pb
INFO: P0_gg_h1xdxdx/G6.8 is at 0.068 +- 0.00187 pb. Now submitting iteration #4.
INFO: Idle: 0, Running: 8, Completed: 156 [ 6m 57s ]
INFO: Survey finished for P0_gg_h1xdxdx/G1.8 at 0.00129 +- 0.000444 pb
INFO: Survey finished for P0_gg_h1xdxdx/G3.5 at 0.2277 +- 0.00125 pb
INFO: Idle: 0, Running: 5, Completed: 160 [ 7m 1s ]
INFO: P0_gg_h1xdxdx/G3.8 is at 0.0256 +- 0.00172 pb. Now submitting iteration #5.
INFO: Idle: 0, Running: 5, Completed: 163 [ 7m 26s ]
INFO: P0_gg_h1xdxdx/G4.8 is at 0.00423 +- 0.000551 pb. Now submitting iteration #5.
INFO: Idle: 0, Running: 5, Completed: 166 [ 8m 17s ]
INFO: Survey finished for P0_gg_h1xdxdx/G6.8 at 0.0668 +- 0.000988 pb
INFO: Idle: 0, Running: 3, Completed: 169 [ 9m 11s ]
INFO: Survey finished for P0_gg_h1xdxdx/G3.8 at 0.0271 +- 0.000546 pb
INFO: Idle: 0, Running: 1, Completed: 172 [ 9m 25s ]
INFO: Survey finished for P0_gg_h1xdxdx/G4.8 at 0.00468 +- 0.000155 pb
INFO: Idle: 0, Running: 0, Completed: 174 [ 9m 25s ]
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 46.10325024446438 pb^-1
INFO: need to improve 1 channels
- Current estimate of cross-section: 249.4401140705 +- 0.607850053933017
INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 16h26 ]
INFO: Idle: 0, Running: 3, Completed: 1 [ 7.3s ]
INFO: P0_gg_h1xdxdx/G5 is at 3438/11367 (30%) event. Resubmit 4 job at iteration 1.
INFO: Idle: 0, Running: 3, Completed: 6 [ 14.7s ]
INFO: P0_gg_h1xdxdx/G5 is at 6880/11364 (61%) event. Resubmit 6 job at iteration 2.
INFO: Idle: 0, Running: 5, Completed: 11 [ 24.6s ]
INFO: found enough event for P0_gg_h1xdxdx/G5
INFO: Idle: 0, Running: 0, Completed: 17 [ 26.2s ]
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: 38m22s
     Cross-section : 248.1 +- 0.04692 pb
     Nb of events : 10000

fail
Failed to access python version of LHAPDF: If the python interface to LHAPDF is available on your system, try adding its location to the PYTHONPATH environment variable and theLHAPDF library location to LD_LIBRARY_PATH (linux) or DYLD_LIBRARY_PATH (mac os x).The required LD_LIBRARY_PATH is /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/lhapdf6_py3//lib
INFO: can not run systematics since can not link python to lhapdf
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 > h1 xd xd~
INFO: options:
INFO: detected model: Pseudoscalar_2HDM. Loading...
set fastjet to /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/fastjet/bin/fastjet-config
/home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/ninja/bin does not seem to correspond to a valid ninja lib . Please enter the full PATH/TO/ninja/lib .
You will NOT be able to run ninja otherwise.

set collier to /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/collier
set lhapdf to /home/faisal/Downloads/MG5_aMC_v3_4_1/HEPTools/lhapdf6_py3/bin/lhapdf-config
Set group_subprocesses to Auto
Note that you need to regenerate all processes
import model heft PARAM_CARD_PATH [--bypass_check]
set spinmode none
decay h1 > a a
launch
INFO: Will use seed 867728678
INFO: generate 10000 decay event for particle h1
Command "generate_events run_01" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: h1 > a a. Please enter a valid process
quit
INFO: storing files of previous run
gzipping output file: unweighted_events.lhe
INFO: Done
INFO:
more information in /home/faisal/Downloads/MG5_aMC_v3_4_1/final_decay_m/index.html
MG5_aMC>

What is your madspin_card?

You need to adapt the line
import model heft PARAM_CARD_PATH [--bypass_check]
and change
PARAM_CARD_PATH
by a path to a valid (and consistent with your original spectrum --in particular the mass of h1--)
if the automatic check is creating un-reasonable issue then you have to use the "--bypass_check" otherwise you do not have to include it. (in no circumstence you have to keep the "[" and "]" those indicates the the argument is an option.

Cheers,

Olivier