HelaS: HELAS_diffs.txt

File HELAS_diffs.txt, 6.5 KB (added by anonymous, 6 years ago)
Line 
1---------------------------------------------------------------------------------------
22007-06-02: Differences between the old smadgraph HELAS (from the smadgrapg website)
3and the MG/MEv4 HELAS.
4---------------------------------------------------------------------------------------
5
6These files/routines are in the smadgraph helas, but not in the madgraph.
7fvidmt.F
8fvodmt.F
9hvvxxu.F
10iovdmt.F
11
12
13These files/routines are in the madgraph helas, but not in the smadgraph.
14hstlxx.F  --> 2 Higgs Effective theory (2heft), not tested for the pseudo-scalar higgs's
15httsxx.F  --> 2heft
16httxxx.F  --> Higgs Effective theory (heft)
17hvvhxx.F  --> heft
18hvvshx.F  --> 2heft
19hvvtxx.F  --> old routine, replaced by uvvxxx.F
20hvvvxx.F  --> heft
21jvshxx.F  --> heft
22jvsshx.F  --> 2heft
23jvstxx.F  --> heft
24jvtxxx.F  --> heft
25jvvsxx.F  --> heft
26jvvtlx.F  --> 2heft
27ssstxx.F  --> old routine, replaced by ttsxxx.F
28sstlxx.F  --> 2heft
29ttssxx.F  --> 2heft
30ttsxxx.F  --> heft
31usslxx.F  --> 2heft
32utssxx.F  --> 2heft
33utsxxx.F  --> heft
34uvvvlx.F  --> 2heft
35uvvxxx.F  --> heft
36vvshxx.F  --> heft
37vvsshx.F  --> 2heft
38vvstxx.F  --> old routine, replaced by vvtxx.F
39vvtxxx.F  --> heft
40vvvsxx.F  --> heft
41vvvtlx.F  --> 2heft
42dimensions.inc  --> A file in which the dimensions of the wavefunctions of the above routines are specified.
43
44See also the HELAS_HEFT.ps file with the description/manual of these routines.
45
46
47
48These are the files/routines that are different in the smadgraph and the madgraph helas:
49ggggxx.F
50hioxxx.F
51j3xxxx.F
52jioxxx.F
53jssxxx.F
54jvssxx.F
55jvsxxx.F
56jvvxxx.F
57jw3wxx.F
58jwwwxx.F
59
60**The difference is mainly the implementation of the finite width of the vector bosons.**
61Here are the diff's between the files:
62
63desktop:~$ diff smadgraph/ggggxx.F  MG_ME_V4.1.33/HELAS/ggggxx.F
6429a30,32
65      save dv1,dv2,dv3, dv4
66c      save dv1,dv2,dv3,dv4,dvertx,v12,v13,v14,v23,v24,v34
67
68132a136,141
69
70c      if (abs(dvertx) .gt. 1d40) then
71c         write(*,*) 'Careful',abs(dvertx)
72c         write(*,*) v12,v13,v14
73c         write(*,*) v23,v24,v34
74c      endif
75
76
77desktop:~$ diff smadgraph/hioxxx.F  MG_ME_V4.1.33/HELAS/hioxxx.F
7821c21
79<       double precision rZero
80---
81      double precision rZero, cZero
82
83
84desktop:~$ diff smadgraph/j3xxxx.F  MG_ME_V4.1.33/HELAS/j3xxxx.F
8523a24
86      double complex cm2  ! mass**2- I Gamma mass (Fabio)
87132a134,135
88
89
90146,147d148
91<       csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/zm2
92<       csr = (q(0)*c0r-q(1)*c1r-q(2)*c2r-q(3)*c3r)/zm2
93148a150,156
94c     Fabio's implementation of the fixed width
95      cm2=dcmplx( zm2, -zmw )
96c     csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/zm2
97c     csr = (q(0)*c0r-q(1)*c1r-q(2)*c2r-q(3)*c3r)/zm2
98      csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/cm2
99      csr = (q(0)*c0r-q(1)*c1r-q(2)*c2r-q(3)*c3r)/cm2
100
101
102
103desktop:~$ diff smadgraph/jioxxx.F  MG_ME_V4.1.33/HELAS/jioxxx.F
10419a20,21
105      double complex cm2 ! mass**2- I Gamma mass (Fabio)
106
10792,93c94,96
108< c  For the running width, use below instead of the above d.
109< c         d = rOne/dcmplx( q2-vm2, max(vwidth*q2/vmass,rZero) )
110---
111c     For the running width, use below instead of the above d.
112c     d = rOne/dcmplx( q2-vm2, max(vwidth*q2/vmass,rZero) )
113
114111c114,118
115<          cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/vm2
116---
117
118c     Fabio's implementation of the fixed width
119         cm2=dcmplx( vm2, -vmass*vwidth )
120c     cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/vm2
121         cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/cm2
122
123
124
125
126desktop:~$ diff smadgraph/jssxxx.F  MG_ME_V4.1.33/HELAS/jssxxx.F
12734a35
128      double complex cm2        ! mass**2- I Gamma mass (Fabio)
129111a113
130
131126,129c128,137
132<          jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/vm2)
133<          jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/vm2)
134<          jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/vm2)
135<          jss(4) = adg*( (pp(3)-pa(3)) - q(3)*m2d/vm2)
136---
137c     Fabio's implementation of the fixed width
138         cm2=dcmplx( vm2, -vmass*vwidth )
139c     jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/vm2)
140c     jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/vm2)
141c     jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/vm2)
142c     jss(4) = adg*( (pp(3)-pa(3)) - q(3)*m2d/vm2)
143         jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/cm2)
144         jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/cm2)
145         jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/cm2)
146         jss(4) = adg*( (pp(3)-pa(3)) - q(3)*m2d/cm2)
147
148
149
150
151desktop:~$ diff smadgraph/jvssxx.F  MG_ME_V4.1.33/HELAS/jvssxx.F
15220,21c20,22
153<       double complex vc(6),s1(3),s2(3),gc,jvss(6),dg,vk
154<       double precision q(0:3),vmass,vwidth,q2,vm2
155---
156      double complex vc(6),s1(3),s2(3),gc,jvss(6),dg
157      double complex cm2        ! mass**2- I Gamma mass (Fabio)
158      double precision q(0:3),vmass,vwidth,q2,vk,vm2
159104c105,108
160<          vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/vm2
161---
162c     Fabio's implementation of the fixed width
163         cm2=dcmplx( vm2, -vmass*vwidth )
164c     vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/vm2
165         vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/cm2
166
167
168
169desktop:~$ diff smadgraph/jvsxxx.F  MG_ME_V4.1.33/HELAS/jvsxxx.F
17018a19
171      double complex cm2        ! mass**2- I Gamma mass (Fabio)
17295c96,99
173<          vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/vm2
174---
175c     Fabio's implementation of the fixed width
176         cm2=dcmplx( vm2, -vmass*vwidth )
177c     vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/vm2
178         vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/cm2
179
180
181
182desktop:~$ diff smadgraph/jvvxxx.F  MG_ME_V4.1.33/HELAS/jvvxxx.F
18329a30
184      double complex cm2        ! mass**2- I Gamma mass (Fabio)
185125d125
186<          js = (v12*(-m1+m2) +s11*s12 -s21*s22)/vm2
187126a127,131
188c     Fabio's implementation of the fixed width
189         cm2=dcmplx( vm2, -vmass*vwidth )
190c     js = (v12*(-m1+m2) +s11*s12 -s21*s22)/vm2
191         js = (v12*(-m1+m2) +s11*s12 -s21*s22)/cm2
192
193
194
195
196desktop:~$ diff smadgraph/jw3wxx.F  MG_ME_V4.1.33/HELAS/jw3wxx.F
19753a54
198      double complex cm2        ! mass**2- I Gamma mass (Fabio)
199201,202c202,206
200<          jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mv2
201<
202---
203c     Fabio's implementation of the fixed width
204         cm2=dcmplx( mv2, -dmv*dwv)
205c     jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mv2
206         jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/cm2
207
208
209
210desktop:~$ diff smadgraph/jwwwxx.F  MG_ME_V4.1.33/HELAS/jwwwxx.F
21134a35
212      double complex cm2        ! mass**2- I Gamma mass (Fabio)
213165c166,169
214<       jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mw2
215---
216c     Fabio's implementation of the fixed width
217      cm2=dcmplx( mw2, -dmw*dww )
218c     jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mw2
219      jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/cm2
220