diff --git a/HHVBF_UFO/restrict_vbf_c2v0_kl0.dat b/HHVBF_UFO/restrict_vbf_c2v0_kl0.dat
new file mode 100644
index 0000000000000000000000000000000000000000..881dc061fa2698c0755fa05ef9cccd7c4da3f75d
--- /dev/null
+++ b/HHVBF_UFO/restrict_vbf_c2v0_kl0.dat
@@ -0,0 +1,143 @@
+######################################################################
+## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO  #####################
+######################################################################
+
+###################################
+## INFORMATION FOR SMINPUTS
+###################################
+Block SMINPUTS 
+    1 1.322330e+02 # aEWM1 
+    2 1.166370e-05 # Gf 
+    3 1.184000e-01 # aS 
+
+###################################
+## INFORMATION FOR MASS
+###################################
+Block MASS 
+    1 5.040000e-03 # MD 
+    2 2.550000e-03 # MU 
+    3 1.010000e-01 # MS 
+    4 1.270000e+00 # MC 
+    5 4.700000e+00 # MB 
+    6 1.720000e+02 # MT 
+   11 5.110000e-04 # Me 
+   13 1.056600e-01 # MMU 
+   15 1.777000e+00 # MTA 
+   23 9.118760e+01 # MZ 
+   25 1.250000e+02 # MH 
+##  Not dependent paramater.
+## Those values should be edited following analytical the 
+## analytical expression. Some generator could simply ignore 
+## those values and use the analytical expression
+  22 0.000000 # a : 0.0 
+  24 80.384851 # W+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  21 0.000000 # g : 0.0 
+  9000001 0.000000 # ghA : 0.0 
+  9000003 80.384851 # ghWp : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  9000004 80.384851 # ghWm : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  9000005 0.000000 # ghG : 0.0 
+  12 0.000000 # ve : 0.0 
+  14 0.000000 # vm : 0.0 
+  16 0.000000 # vt : 0.0 
+  251 80.384851 # G+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  9000002 91.187600 # ghZ : MZ 
+  250 91.187600 # G0 : MZ 
+
+###################################
+## INFORMATION FOR DECAY
+###################################
+DECAY   6 1.508336e+00 
+DECAY  23 2.495200e+00 
+DECAY  24 2.085000e+00 
+DECAY  25 4.070000e-03 
+##  Not dependent paramater.
+## Those values should be edited following analytical the 
+## analytical expression. Some generator could simply ignore 
+## those values and use the analytical expression
+DECAY  22 0.000000 # a : 0.0 
+DECAY  21 0.000000 # g : 0.0 
+DECAY  9000001 0.000000 # ghA : 0.0 
+DECAY  9000005 0.000000 # ghG : 0.0 
+DECAY  12 0.000000 # ve : 0.0 
+DECAY  14 0.000000 # vm : 0.0 
+DECAY  16 0.000000 # vt : 0.0 
+DECAY  11 0.000000 # e- : 0.0 
+DECAY  13 0.000000 # mu- : 0.0 
+DECAY  15 0.000000 # ta- : 0.0 
+DECAY  2 0.000000 # u : 0.0 
+DECAY  4 0.000000 # c : 0.0 
+DECAY  1 0.000000 # d : 0.0 
+DECAY  3 0.000000 # s : 0.0 
+DECAY  5 0.000000 # b : 0.0 
+DECAY  9000002 2.495200 # ghZ : WZ 
+DECAY  9000003 2.085000 # ghWp : WW 
+DECAY  9000004 2.085000 # ghWm : WW 
+DECAY  250 2.495200 # G0 : WZ 
+DECAY  251 2.085000 # G+ : WW 
+
+###################################
+## INFORMATION FOR CKMBLOCK
+###################################
+Block CKMBLOCK 
+    1 2.277360e-01 # cabi 
+
+###################################
+## INFORMATION FOR NEW
+###################################
+Block NEW 
+    1 1.000000 # CV 
+    2 0.000000 # C2V 
+    3 0.000000 # C3 
+
+###################################
+## INFORMATION FOR YUKAWA
+###################################
+Block YUKAWA 
+    1 5.040000e-03 # ymdo 
+    2 2.550000e-03 # ymup 
+    3 1.010000e-01 # yms 
+    4 1.270000e+00 # ymc 
+    5 4.700000e+00 # ymb 
+    6 1.720000e+02 # ymt 
+   11 5.110000e-04 # yme 
+   13 1.056600e-01 # ymm 
+   15 1.777000e+00 # ymtau 
+#===========================================================
+# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
+#===========================================================
+
+Block QNUMBERS 9000001  # ghA 
+        1 0  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000002  # ghZ 
+        1 0  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000003  # ghWp 
+        1 3  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000004  # ghWm 
+        1 -3  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000005  # ghG 
+        1 0  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 8  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 250  # G0 
+        1 0  # 3 times electric charge
+        2 1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 0  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 251  # G+ 
+        1 3  # 3 times electric charge
+        2 1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
diff --git a/HHVBF_UFO/restrict_vbf_kl0.dat b/HHVBF_UFO/restrict_vbf_kl0.dat
new file mode 100644
index 0000000000000000000000000000000000000000..6916db72ab8a227d48397d30c76659c4fb4a86f9
--- /dev/null
+++ b/HHVBF_UFO/restrict_vbf_kl0.dat
@@ -0,0 +1,143 @@
+######################################################################
+## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO  #####################
+######################################################################
+
+###################################
+## INFORMATION FOR SMINPUTS
+###################################
+Block SMINPUTS 
+    1 1.322330e+02 # aEWM1 
+    2 1.166370e-05 # Gf 
+    3 1.184000e-01 # aS 
+
+###################################
+## INFORMATION FOR MASS
+###################################
+Block MASS 
+    1 5.040000e-03 # MD 
+    2 2.550000e-03 # MU 
+    3 1.010000e-01 # MS 
+    4 1.270000e+00 # MC 
+    5 4.700000e+00 # MB 
+    6 1.720000e+02 # MT 
+   11 5.110000e-04 # Me 
+   13 1.056600e-01 # MMU 
+   15 1.777000e+00 # MTA 
+   23 9.118760e+01 # MZ 
+   25 1.250000e+02 # MH 
+##  Not dependent paramater.
+## Those values should be edited following analytical the 
+## analytical expression. Some generator could simply ignore 
+## those values and use the analytical expression
+  22 0.000000 # a : 0.0 
+  24 80.384851 # W+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  21 0.000000 # g : 0.0 
+  9000001 0.000000 # ghA : 0.0 
+  9000003 80.384851 # ghWp : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  9000004 80.384851 # ghWm : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  9000005 0.000000 # ghG : 0.0 
+  12 0.000000 # ve : 0.0 
+  14 0.000000 # vm : 0.0 
+  16 0.000000 # vt : 0.0 
+  251 80.384851 # G+ : cmath.sqrt(MZ**2/2. + cmath.sqrt(MZ**4/4. - (aEW*cmath.pi*MZ**2)/(Gf*cmath.sqrt(2)))) 
+  9000002 91.187600 # ghZ : MZ 
+  250 91.187600 # G0 : MZ 
+
+###################################
+## INFORMATION FOR DECAY
+###################################
+DECAY   6 1.508336e+00 
+DECAY  23 2.495200e+00 
+DECAY  24 2.085000e+00 
+DECAY  25 4.070000e-03 
+##  Not dependent paramater.
+## Those values should be edited following analytical the 
+## analytical expression. Some generator could simply ignore 
+## those values and use the analytical expression
+DECAY  22 0.000000 # a : 0.0 
+DECAY  21 0.000000 # g : 0.0 
+DECAY  9000001 0.000000 # ghA : 0.0 
+DECAY  9000005 0.000000 # ghG : 0.0 
+DECAY  12 0.000000 # ve : 0.0 
+DECAY  14 0.000000 # vm : 0.0 
+DECAY  16 0.000000 # vt : 0.0 
+DECAY  11 0.000000 # e- : 0.0 
+DECAY  13 0.000000 # mu- : 0.0 
+DECAY  15 0.000000 # ta- : 0.0 
+DECAY  2 0.000000 # u : 0.0 
+DECAY  4 0.000000 # c : 0.0 
+DECAY  1 0.000000 # d : 0.0 
+DECAY  3 0.000000 # s : 0.0 
+DECAY  5 0.000000 # b : 0.0 
+DECAY  9000002 2.495200 # ghZ : WZ 
+DECAY  9000003 2.085000 # ghWp : WW 
+DECAY  9000004 2.085000 # ghWm : WW 
+DECAY  250 2.495200 # G0 : WZ 
+DECAY  251 2.085000 # G+ : WW 
+
+###################################
+## INFORMATION FOR CKMBLOCK
+###################################
+Block CKMBLOCK 
+    1 2.277360e-01 # cabi 
+
+###################################
+## INFORMATION FOR NEW
+###################################
+Block NEW 
+    1 1.000000 # CV 
+    2 1.000000 # C2V 
+    3 0.000000 # C3 
+
+###################################
+## INFORMATION FOR YUKAWA
+###################################
+Block YUKAWA 
+    1 5.040000e-03 # ymdo 
+    2 2.550000e-03 # ymup 
+    3 1.010000e-01 # yms 
+    4 1.270000e+00 # ymc 
+    5 4.700000e+00 # ymb 
+    6 1.720000e+02 # ymt 
+   11 5.110000e-04 # yme 
+   13 1.056600e-01 # ymm 
+   15 1.777000e+00 # ymtau 
+#===========================================================
+# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
+#===========================================================
+
+Block QNUMBERS 9000001  # ghA 
+        1 0  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000002  # ghZ 
+        1 0  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000003  # ghWp 
+        1 3  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000004  # ghWm 
+        1 -3  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 9000005  # ghG 
+        1 0  # 3 times electric charge
+        2 -1  # number of spin states (2S+1)
+        3 8  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 250  # G0 
+        1 0  # 3 times electric charge
+        2 1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 0  # Particle/Antiparticle distinction (0=own anti)
+Block QNUMBERS 251  # G+ 
+        1 3  # 3 times electric charge
+        2 1  # number of spin states (2S+1)
+        3 1  # colour rep (1: singlet, 3: triplet, 8: octet)
+        4 1  # Particle/Antiparticle distinction (0=own anti)