From f71a561ea5cba997c8a0850d2d3f7b3b0b3c43bc Mon Sep 17 00:00:00 2001
From: Ian-Woo Kim <ianwookim@gmail.com>
Date: Fri, 10 Jul 2015 14:51:32 +0200
Subject: [PATCH] pseudoscalar of heftNLO now works well!

---
 modelrepo/heftNLO/Fortran/functions.f | 7 +------
 modelrepo/heftNLO/coupling_orders.py  | 4 ++--
 modelrepo/heftNLO/couplings.py        | 2 +-
 modelrepo/heftNLO/vertices.py         | 3 ++-
 template/param_card_heftNLO.dat.st    | 2 +-
 5 files changed, 7 insertions(+), 11 deletions(-)

diff --git a/modelrepo/heftNLO/Fortran/functions.f b/modelrepo/heftNLO/Fortran/functions.f
index efbdc96..eaf3c08 100644
--- a/modelrepo/heftNLO/Fortran/functions.f
+++ b/modelrepo/heftNLO/Fortran/functions.f
@@ -40,7 +40,6 @@ DOUBLE COMPLEX FUNCTION FormFactor2( HALFS )
       DOUBLE COMPLEX beta
       DOUBLE COMPLEX BRA1
       DOUBLE COMPLEX BRA
-      DOUBLE COMPLEX CORRF
       REAL PIE
 
       INCLUDE 'coupl.inc'
@@ -49,13 +48,9 @@ DOUBLE COMPLEX FUNCTION FormFactor2( HALFS )
       PIE=3.14159265358979323846264
       S = 2.0*HALFS
       mt = MDL_MT
-      !wh1 = MDL_WH1
-      wh1 = 11.82
       mh1 = MDL_MP
       tau = 4*mt**2 / S
 
-      CORRF = (S - mh1**2) / ((S-mh1**2)+(0D0,1D0)*mh1*wh1)
- 
       IF (REAL(tau) .LT. 1.0) THEN  
         beta = SQRT(1-tau)
         BRA1 = LOG((1+beta)/(1-beta)) - (0D0,1D0)*PIE
@@ -64,7 +59,7 @@ DOUBLE COMPLEX FUNCTION FormFactor2( HALFS )
         BRA = -2.0*(ASIN(1/SQRT(tau)))**2
       END IF
          
-      FormFactor2 = 4.0*mt / (16.0 * PIE**2 * S)* BRA * CORRF
+      FormFactor2 = 0.5*4.0*mt / (16.0 * PIE**2 * S)* BRA
       !IF (REAL(tau) .LT. 1.0) THEN  
       !  WRITE(0,*) SQRT(S),CORRF !,FormFactor2
       !ENDIF
diff --git a/modelrepo/heftNLO/coupling_orders.py b/modelrepo/heftNLO/coupling_orders.py
index 306740e..84cc43b 100644
--- a/modelrepo/heftNLO/coupling_orders.py
+++ b/modelrepo/heftNLO/coupling_orders.py
@@ -15,10 +15,10 @@
                     hierarchy = 2)
 
 HIG = CouplingOrder(name = 'HIG',
-                    expansion_order = 1,
+                    expansion_order = 99,
                     hierarchy = 2)
 
 HIW = CouplingOrder(name = 'HIW',
-                    expansion_order = 1,
+                    expansion_order = 99,
                     hierarchy = 2)
 
diff --git a/modelrepo/heftNLO/couplings.py b/modelrepo/heftNLO/couplings.py
index 5b4e9fb..debcd00 100644
--- a/modelrepo/heftNLO/couplings.py
+++ b/modelrepo/heftNLO/couplings.py
@@ -78,7 +78,7 @@
 #                 order = {'HIG':1})
 # complex(0,1)*
 GC_16 = Coupling(name = 'GC_16',
-                 value = '-(G**2*yt/cmath.sqrt(2))',
+                 value = '-complex(0,1)*(G**2*yt/cmath.sqrt(2))',
                  order = {'HIG':1})
 
 GC_17 = Coupling(name = 'GC_17',
diff --git a/modelrepo/heftNLO/vertices.py b/modelrepo/heftNLO/vertices.py
index fa91cd6..11ce2a0 100644
--- a/modelrepo/heftNLO/vertices.py
+++ b/modelrepo/heftNLO/vertices.py
@@ -891,9 +891,10 @@
                color = [ '1' ],
                lorentz = [ L.FFV2 ],
                couplings = {(0,0):C.GC_59})
+#####
 
 V_146 = Vertex(name = 'V_146',
                particles = [ P.t__tilde__, P.t, P.h1 ],
                color = [ 'Identity(1,2)' ],
-               lorentz = [ L.FFS2 ],
+               lorentz = [ L.FFS1 ],
                couplings = {(0,0):C.GC_103})
diff --git a/template/param_card_heftNLO.dat.st b/template/param_card_heftNLO.dat.st
index ca1514a..66772fd 100644
--- a/template/param_card_heftNLO.dat.st
+++ b/template/param_card_heftNLO.dat.st
@@ -59,7 +59,7 @@ DECAY   6 1.491500e+00 # WT
 DECAY  23 2.441404e+00 # WZ 
 DECAY  24 2.047600e+00 # WW 
 DECAY  25 $WH$     # WH 
-DECAY 9000006 0.00000e+00  # WH1 $WA$
+DECAY 9000006 $WA$  # WH1 
 ## Dependent parameters, given by model restrictions.
 ## Those values should be edited following the 
 ## analytical expression. MG5 ignores those values