| 1 | ---------------------------------------------------------------------------------------
|
|---|
| 2 | 2007-06-02: Differences between the old smadgraph HELAS (from the smadgrapg website)
|
|---|
| 3 | and the MG/MEv4 HELAS.
|
|---|
| 4 | ---------------------------------------------------------------------------------------
|
|---|
| 5 |
|
|---|
| 6 | These files/routines are in the smadgraph helas, but not in the madgraph.
|
|---|
| 7 | fvidmt.F
|
|---|
| 8 | fvodmt.F
|
|---|
| 9 | hvvxxu.F
|
|---|
| 10 | iovdmt.F
|
|---|
| 11 |
|
|---|
| 12 |
|
|---|
| 13 | These files/routines are in the madgraph helas, but not in the smadgraph.
|
|---|
| 14 | hstlxx.F --> 2 Higgs Effective theory (2heft), not tested for the pseudo-scalar higgs's
|
|---|
| 15 | httsxx.F --> 2heft
|
|---|
| 16 | httxxx.F --> Higgs Effective theory (heft)
|
|---|
| 17 | hvvhxx.F --> heft
|
|---|
| 18 | hvvshx.F --> 2heft
|
|---|
| 19 | hvvtxx.F --> old routine, replaced by uvvxxx.F
|
|---|
| 20 | hvvvxx.F --> heft
|
|---|
| 21 | jvshxx.F --> heft
|
|---|
| 22 | jvsshx.F --> 2heft
|
|---|
| 23 | jvstxx.F --> heft
|
|---|
| 24 | jvtxxx.F --> heft
|
|---|
| 25 | jvvsxx.F --> heft
|
|---|
| 26 | jvvtlx.F --> 2heft
|
|---|
| 27 | ssstxx.F --> old routine, replaced by ttsxxx.F
|
|---|
| 28 | sstlxx.F --> 2heft
|
|---|
| 29 | ttssxx.F --> 2heft
|
|---|
| 30 | ttsxxx.F --> heft
|
|---|
| 31 | usslxx.F --> 2heft
|
|---|
| 32 | utssxx.F --> 2heft
|
|---|
| 33 | utsxxx.F --> heft
|
|---|
| 34 | uvvvlx.F --> 2heft
|
|---|
| 35 | uvvxxx.F --> heft
|
|---|
| 36 | vvshxx.F --> heft
|
|---|
| 37 | vvsshx.F --> 2heft
|
|---|
| 38 | vvstxx.F --> old routine, replaced by vvtxx.F
|
|---|
| 39 | vvtxxx.F --> heft
|
|---|
| 40 | vvvsxx.F --> heft
|
|---|
| 41 | vvvtlx.F --> 2heft
|
|---|
| 42 | dimensions.inc --> A file in which the dimensions of the wavefunctions of the above routines are specified.
|
|---|
| 43 |
|
|---|
| 44 | See also the HELAS_HEFT.ps file with the description/manual of these routines.
|
|---|
| 45 |
|
|---|
| 46 |
|
|---|
| 47 |
|
|---|
| 48 | These are the files/routines that are different in the smadgraph and the madgraph helas:
|
|---|
| 49 | ggggxx.F
|
|---|
| 50 | hioxxx.F
|
|---|
| 51 | j3xxxx.F
|
|---|
| 52 | jioxxx.F
|
|---|
| 53 | jssxxx.F
|
|---|
| 54 | jvssxx.F
|
|---|
| 55 | jvsxxx.F
|
|---|
| 56 | jvvxxx.F
|
|---|
| 57 | jw3wxx.F
|
|---|
| 58 | jwwwxx.F
|
|---|
| 59 |
|
|---|
| 60 | **The difference is mainly the implementation of the finite width of the vector bosons.**
|
|---|
| 61 | Here are the diff's between the files:
|
|---|
| 62 |
|
|---|
| 63 | desktop:~$ diff smadgraph/ggggxx.F MG_ME_V4.1.33/HELAS/ggggxx.F
|
|---|
| 64 | 29a30,32
|
|---|
| 65 | save dv1,dv2,dv3, dv4
|
|---|
| 66 | c save dv1,dv2,dv3,dv4,dvertx,v12,v13,v14,v23,v24,v34
|
|---|
| 67 |
|
|---|
| 68 | 132a136,141
|
|---|
| 69 |
|
|---|
| 70 | c if (abs(dvertx) .gt. 1d40) then
|
|---|
| 71 | c write(*,*) 'Careful',abs(dvertx)
|
|---|
| 72 | c write(*,*) v12,v13,v14
|
|---|
| 73 | c write(*,*) v23,v24,v34
|
|---|
| 74 | c endif
|
|---|
| 75 |
|
|---|
| 76 |
|
|---|
| 77 | desktop:~$ diff smadgraph/hioxxx.F MG_ME_V4.1.33/HELAS/hioxxx.F
|
|---|
| 78 | 21c21
|
|---|
| 79 | < double precision rZero
|
|---|
| 80 | ---
|
|---|
| 81 | double precision rZero, cZero
|
|---|
| 82 |
|
|---|
| 83 |
|
|---|
| 84 | desktop:~$ diff smadgraph/j3xxxx.F MG_ME_V4.1.33/HELAS/j3xxxx.F
|
|---|
| 85 | 23a24
|
|---|
| 86 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 87 | 132a134,135
|
|---|
| 88 |
|
|---|
| 89 |
|
|---|
| 90 | 146,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
|
|---|
| 93 | 148a150,156
|
|---|
| 94 | c Fabio's implementation of the fixed width
|
|---|
| 95 | cm2=dcmplx( zm2, -zmw )
|
|---|
| 96 | c csl = (q(0)*c0l-q(1)*c1l-q(2)*c2l-q(3)*c3l)/zm2
|
|---|
| 97 | c 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 |
|
|---|
| 103 | desktop:~$ diff smadgraph/jioxxx.F MG_ME_V4.1.33/HELAS/jioxxx.F
|
|---|
| 104 | 19a20,21
|
|---|
| 105 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 106 |
|
|---|
| 107 | 92,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 | ---
|
|---|
| 111 | c For the running width, use below instead of the above d.
|
|---|
| 112 | c d = rOne/dcmplx( q2-vm2, max(vwidth*q2/vmass,rZero) )
|
|---|
| 113 |
|
|---|
| 114 | 111c114,118
|
|---|
| 115 | < cs = (q(0)*c0-q(1)*c1-q(2)*c2-q(3)*c3)/vm2
|
|---|
| 116 | ---
|
|---|
| 117 |
|
|---|
| 118 | c Fabio's implementation of the fixed width
|
|---|
| 119 | cm2=dcmplx( vm2, -vmass*vwidth )
|
|---|
| 120 | c 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 |
|
|---|
| 126 | desktop:~$ diff smadgraph/jssxxx.F MG_ME_V4.1.33/HELAS/jssxxx.F
|
|---|
| 127 | 34a35
|
|---|
| 128 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 129 | 111a113
|
|---|
| 130 |
|
|---|
| 131 | 126,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 | ---
|
|---|
| 137 | c Fabio's implementation of the fixed width
|
|---|
| 138 | cm2=dcmplx( vm2, -vmass*vwidth )
|
|---|
| 139 | c jss(1) = adg*( (pp(0)-pa(0)) - q(0)*m2d/vm2)
|
|---|
| 140 | c jss(2) = adg*( (pp(1)-pa(1)) - q(1)*m2d/vm2)
|
|---|
| 141 | c jss(3) = adg*( (pp(2)-pa(2)) - q(2)*m2d/vm2)
|
|---|
| 142 | c 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 |
|
|---|
| 151 | desktop:~$ diff smadgraph/jvssxx.F MG_ME_V4.1.33/HELAS/jvssxx.F
|
|---|
| 152 | 20,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
|
|---|
| 159 | 104c105,108
|
|---|
| 160 | < vk = (q(0)*vc(1)-q(1)*vc(2)-q(2)*vc(3)-q(3)*vc(4))/vm2
|
|---|
| 161 | ---
|
|---|
| 162 | c Fabio's implementation of the fixed width
|
|---|
| 163 | cm2=dcmplx( vm2, -vmass*vwidth )
|
|---|
| 164 | c 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 |
|
|---|
| 169 | desktop:~$ diff smadgraph/jvsxxx.F MG_ME_V4.1.33/HELAS/jvsxxx.F
|
|---|
| 170 | 18a19
|
|---|
| 171 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 172 | 95c96,99
|
|---|
| 173 | < vk = (-q(0)*vc(1)+q(1)*vc(2)+q(2)*vc(3)+q(3)*vc(4))/vm2
|
|---|
| 174 | ---
|
|---|
| 175 | c Fabio's implementation of the fixed width
|
|---|
| 176 | cm2=dcmplx( vm2, -vmass*vwidth )
|
|---|
| 177 | c 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 |
|
|---|
| 182 | desktop:~$ diff smadgraph/jvvxxx.F MG_ME_V4.1.33/HELAS/jvvxxx.F
|
|---|
| 183 | 29a30
|
|---|
| 184 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 185 | 125d125
|
|---|
| 186 | < js = (v12*(-m1+m2) +s11*s12 -s21*s22)/vm2
|
|---|
| 187 | 126a127,131
|
|---|
| 188 | c Fabio's implementation of the fixed width
|
|---|
| 189 | cm2=dcmplx( vm2, -vmass*vwidth )
|
|---|
| 190 | c js = (v12*(-m1+m2) +s11*s12 -s21*s22)/vm2
|
|---|
| 191 | js = (v12*(-m1+m2) +s11*s12 -s21*s22)/cm2
|
|---|
| 192 |
|
|---|
| 193 |
|
|---|
| 194 |
|
|---|
| 195 |
|
|---|
| 196 | desktop:~$ diff smadgraph/jw3wxx.F MG_ME_V4.1.33/HELAS/jw3wxx.F
|
|---|
| 197 | 53a54
|
|---|
| 198 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 199 | 201,202c202,206
|
|---|
| 200 | < jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mv2
|
|---|
| 201 | <
|
|---|
| 202 | ---
|
|---|
| 203 | c Fabio's implementation of the fixed width
|
|---|
| 204 | cm2=dcmplx( mv2, -dmv*dwv)
|
|---|
| 205 | c 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 |
|
|---|
| 210 | desktop:~$ diff smadgraph/jwwwxx.F MG_ME_V4.1.33/HELAS/jwwwxx.F
|
|---|
| 211 | 34a35
|
|---|
| 212 | double complex cm2 ! mass**2- I Gamma mass (Fabio)
|
|---|
| 213 | 165c166,169
|
|---|
| 214 | < jq = (jj(0)*q(0)-jj(1)*q(1)-jj(2)*q(2)-jj(3)*q(3))/mw2
|
|---|
| 215 | ---
|
|---|
| 216 | c Fabio's implementation of the fixed width
|
|---|
| 217 | cm2=dcmplx( mw2, -dmw*dww )
|
|---|
| 218 | c 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 |
|
|---|
