Unable to decay for a specific channel using MadSpin

Hi MG team,
I am trying to generate VBF events at NLO level. I am interested in generating events with Higgs decaying through the 4 lepton channel H->ZZ->4l. Since Madgraph doesnt allow me to decay along with generation at NLO( this wasnt the issue for LO process) I am trying to achieve this decay using MadSpin.

To do that, I specifically mentioned the Higgs Decay H->ZZ in my param_card.dat as shown:
```
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
# BR NDA ID1 ID2
2.619e-02 2 23 23 #H->ZZ
```

The BR was taken from CERNYellowReports. I also specifically mentioned the decay process in madspin_card_default.dat (Only default card was generated during output):
```
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
# set spinmode onshell # Use one of the madspin special mode
 set max_weight_ps_point 400 # number of PS to estimate the maximum for each event

# specify the decay for the final state particles
decay h > z z
#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
# running the actual code
launch
```

With the above two changes, I hoped for MadSpin to identify the process and decay events accordingly. However, MadSpin is completely ignoring the mentioned decay mode:
```
(base) dhruvanshu@dhruvanshu-Inspiron-3501:~/MG5_aMC_v3_1_0/MadSpin$ ./madspin
************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
MadSpin>import /home/dhruvanshu/MG5_aMC_v3_1_0/VBF_NLO_2jets_final/Events/run_09/events.lhe
INFO: process: p p > h j j
INFO: options: $$ w+ w- z
define p = g u c d s u~ c~ d~ s~
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
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
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+
Defined multiparticle p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
Defined multiparticle l+ = e+ mu+
define l- = e- mu-
Defined multiparticle l- = e- mu-
define vl = ve vm vt
Defined multiparticle vl = ve vm vt
define vl~ = ve~ vm~ vt~
Defined multiparticle vl~ = ve~ vm~ vt~
Set group_subprocesses to Auto
Note that you need to regenerate all processes
Complex mass already desactivated.
set loop optimized output to True
set loop color flows to False
INFO: detected model: loop_sm. Loading...
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 l- / j / vl / l+ / p / 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+
Set group_subprocesses to Auto
Note that you need to regenerate all processes
Complex mass already desactivated.
set loop optimized output to True
set loop color flows to False
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
MadSpin>decay h > z z
MadSpin>launch
INFO: Will use seed 383589808
WARNING: set the mass of the b-quark to its value in the param_card.dat: 4.7 GeV
INFO: We need to recalculate the branching fractions for h
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /home/dhruvanshu/MG5_aMC_v3_1_0/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 h : ( width = 0.00564955 GeV )
INFO: BR d1 d2
INFO: 9.541999e-01 b~ b
INFO: 4.580011e-02 ta+ ta-
INFO:
INFO: generating the production square matrix element
INFO: generate p p > h j j $$ w+ w- z --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > h j j pert_QCD $$ w+ w- z --no_warning=duplicate;
INFO: Done 29.19
INFO: generating the full matrix element squared (with decay)
INFO: generate p p > h j j $$ w+ w- z, h > z z QCD=99 --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > h j j pert_QCD $$ w+ w- z, h > z z QCD=99 --no_warning=duplicate;
Command "launch" interrupted with error:
InvalidCmd : No matrix elements generated, check overall coupling orders
Please report this bug to developers
```
So I have the following questions:
1) Is my way correct or I am missing something?
2) Is MadSpin failing to decay because its keeping both Z onshell which would obviously fail this decay?
3) Is there a way to force MadGraph to decay Higgs at NLO without going through the MadSpin path?
4) Can we skip the auto width calculation and rather force MadSpin to take the values given? ( Coz in the process I am trying its showing only bb~ and tata~ modes, completely ignoring the zz mode).

Any help would be appreciated.

Thanks,
Dhruvanshu