--------------------------------------------------------------------------------------- 2007-06-02: Differences between the old smadgraph HELAS (from the smadgrapg website) and the MG/MEv4 HELAS. --------------------------------------------------------------------------------------- These files/routines are in the smadgraph helas, but not in the madgraph. fvidmt.F fvodmt.F hvvxxu.F iovdmt.F These files/routines are in the madgraph helas, but not in the smadgraph. hstlxx.F --> 2 Higgs Effective theory (2heft), not tested for the pseudo-scalar higgs's httsxx.F --> 2heft httxxx.F --> Higgs Effective theory (heft) hvvhxx.F --> heft hvvshx.F --> 2heft hvvtxx.F --> old routine, replaced by uvvxxx.F hvvvxx.F --> heft jvshxx.F --> heft jvsshx.F --> 2heft jvstxx.F --> heft jvtxxx.F --> heft jvvsxx.F --> heft jvvtlx.F --> 2heft ssstxx.F --> old routine, replaced by ttsxxx.F sstlxx.F --> 2heft ttssxx.F --> 2heft ttsxxx.F --> heft usslxx.F --> 2heft utssxx.F --> 2heft utsxxx.F --> heft uvvvlx.F --> 2heft uvvxxx.F --> heft vvshxx.F --> heft vvsshx.F --> 2heft vvstxx.F --> old routine, replaced by vvtxx.F vvtxxx.F --> heft vvvsxx.F --> heft vvvtlx.F --> 2heft dimensions.inc --> A file in which the dimensions of the wavefunctions of the above routines are specified. See also the HELAS_HEFT.ps file with the description/manual of these routines. These are the files/routines that are different in the smadgraph and the madgraph helas: ggggxx.F hioxxx.F j3xxxx.F jioxxx.F jssxxx.F jvssxx.F jvsxxx.F jvvxxx.F jw3wxx.F jwwwxx.F **The difference is mainly the implementation of the finite width of the vector bosons.** Here are the diff's between the files: desktop:~$ diff smadgraph/ggggxx.F MG_ME_V4.1.33/HELAS/ggggxx.F 29a30,32 save dv1,dv2,dv3, dv4 c save dv1,dv2,dv3,dv4,dvertx,v12,v13,v14,v23,v24,v34 132a136,141 c if (abs(dvertx) .gt. 1d40) then c write(*,*) 'Careful',abs(dvertx) c write(*,*) v12,v13,v14 c write(*,*) v23,v24,v34 c endif desktop:~$ diff smadgraph/hioxxx.F MG_ME_V4.1.33/HELAS/hioxxx.F 21c21 < double precision rZero --- double precision rZero, cZero desktop:~$ diff smadgraph/j3xxxx.F MG_ME_V4.1.33/HELAS/j3xxxx.F 23a24 double complex cm2 ! mass**2- I Gamma mass (Fabio) 132a134,135 146,147d148 < csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/zm2 < csr = (q(0)*c0r-q(1)*c1r-q(2)*c2r-q(3)*c3r)/zm2 148a150,156 c Fabio's implementation of the fixed width cm2=dcmplx( zm2, -zmw ) c csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/zm2 c csr = (q(0)*c0r-q(1)*c1r-q(2)*c2r-q(3)*c3r)/zm2 csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/cm2 csr = (q(0)*c0r-q(1)*c1r-q(2)*c2r-q(3)*c3r)/cm2 desktop:~$ diff smadgraph/jioxxx.F MG_ME_V4.1.33/HELAS/jioxxx.F 19a20,21 double complex cm2 ! mass**2- I Gamma mass (Fabio) 92,93c94,96 < c For the running width, use below instead of the above d. < c d = rOne/dcmplx( q2-vm2, max(vwidth*q2/vmass,rZero) ) --- c For the running width, use below instead of the above d. c d = rOne/dcmplx( q2-vm2, max(vwidth*q2/vmass,rZero) ) 111c114,118 < cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/vm2 --- c Fabio's implementation of the fixed width cm2=dcmplx( vm2, -vmass*vwidth ) c cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/vm2 cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/cm2 desktop:~$ diff smadgraph/jssxxx.F MG_ME_V4.1.33/HELAS/jssxxx.F 34a35 double complex cm2 ! mass**2- I Gamma mass (Fabio) 111a113 126,129c128,137 < jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/vm2) < jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/vm2) < jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/vm2) < jss(4) = adg*( (pp(3)-pa(3)) - q(3)*m2d/vm2) --- c Fabio's implementation of the fixed width cm2=dcmplx( vm2, -vmass*vwidth ) c jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/vm2) c jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/vm2) c jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/vm2) c jss(4) = adg*( (pp(3)-pa(3)) - q(3)*m2d/vm2) jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/cm2) jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/cm2) jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/cm2) jss(4) = adg*( (pp(3)-pa(3)) - q(3)*m2d/cm2) desktop:~$ diff smadgraph/jvssxx.F MG_ME_V4.1.33/HELAS/jvssxx.F 20,21c20,22 < double complex vc(6),s1(3),s2(3),gc,jvss(6),dg,vk < double precision q(0:3),vmass,vwidth,q2,vm2 --- double complex vc(6),s1(3),s2(3),gc,jvss(6),dg double complex cm2 ! mass**2- I Gamma mass (Fabio) double precision q(0:3),vmass,vwidth,q2,vk,vm2 104c105,108 < vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/vm2 --- c Fabio's implementation of the fixed width cm2=dcmplx( vm2, -vmass*vwidth ) c vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/vm2 vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/cm2 desktop:~$ diff smadgraph/jvsxxx.F MG_ME_V4.1.33/HELAS/jvsxxx.F 18a19 double complex cm2 ! mass**2- I Gamma mass (Fabio) 95c96,99 < vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/vm2 --- c Fabio's implementation of the fixed width cm2=dcmplx( vm2, -vmass*vwidth ) c vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/vm2 vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/cm2 desktop:~$ diff smadgraph/jvvxxx.F MG_ME_V4.1.33/HELAS/jvvxxx.F 29a30 double complex cm2 ! mass**2- I Gamma mass (Fabio) 125d125 < js = (v12*(-m1+m2) +s11*s12 -s21*s22)/vm2 126a127,131 c Fabio's implementation of the fixed width cm2=dcmplx( vm2, -vmass*vwidth ) c js = (v12*(-m1+m2) +s11*s12 -s21*s22)/vm2 js = (v12*(-m1+m2) +s11*s12 -s21*s22)/cm2 desktop:~$ diff smadgraph/jw3wxx.F MG_ME_V4.1.33/HELAS/jw3wxx.F 53a54 double complex cm2 ! mass**2- I Gamma mass (Fabio) 201,202c202,206 < jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mv2 < --- c Fabio's implementation of the fixed width cm2=dcmplx( mv2, -dmv*dwv) c jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mv2 jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/cm2 desktop:~$ diff smadgraph/jwwwxx.F MG_ME_V4.1.33/HELAS/jwwwxx.F 34a35 double complex cm2 ! mass**2- I Gamma mass (Fabio) 165c166,169 < jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mw2 --- c Fabio's implementation of the fixed width cm2=dcmplx( mw2, -dmw*dww ) c jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mw2 jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/cm2