pseudoscalar of heftNLO now works well!

wavewave · Jul 10, 2015 · f71a561 · f71a561
1 parent 3ea0e98
commit f71a561
Show file tree

Hide file tree

Showing 5 changed files with 7 additions and 11 deletions.
diff --git 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
@@ -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
@@ -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
@@ -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
@@ -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