← Back to branch summary

~chirality-flow/chiralityflow/ChiralityFlowMG

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/%..%..%Source%MODEL%intparam_definition.inc

  • Committer: andrew.lifson at lu
  • Date: 2021-09-02 13:57:34 UTC
  • Revision ID: andrew.lifson@thep.lu.se-20210902135734-4eybgli0iljkax9b
added fresh copy of MG5_aMC_v3.2.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/long_ML_SMQCD_optimized/gg_wmtbx/%..%..%Source%MODEL%intparam_definition.inc
 
1
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
2
c      written by the UFO converter
 
3
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
4
 
 
5
C     Parameters that should not be recomputed event by event.
 
6
C     
 
7
      IF(READLHA) THEN
 
8
 
 
9
        G = 2 * DSQRT(AS*PI)  ! for the first init
 
10
 
 
11
        MDL_LHV = 1.000000D+00
 
12
 
 
13
        MDL_CKM33 = 1.000000D+00
 
14
 
 
15
        MDL_CONJG__CKM3X3 = 1.000000D+00
 
16
 
 
17
        MDL_CKM3X3 = 1.000000D+00
 
18
 
 
19
        MDL_CONJG__CKM22 = 1.000000D+00
 
20
 
 
21
        MDL_CKM22 = 1.000000D+00
 
22
 
 
23
        MDL_NCOL = 3.000000D+00
 
24
 
 
25
        MDL_CA = 3.000000D+00
 
26
 
 
27
        MDL_TF = 5.000000D-01
 
28
 
 
29
        MDL_CF = (4.000000D+00/3.000000D+00)
 
30
 
 
31
        MDL_COMPLEXI = DCMPLX(0.000000D+00,1.000000D+00)
 
32
 
 
33
        MDL_MZ__EXP__2 = MDL_MZ**2
 
34
 
 
35
        MDL_MZ__EXP__4 = MDL_MZ**4
 
36
 
 
37
        MDL_SQRT__2 = SQRT(DCMPLX(2.000000D+00))
 
38
 
 
39
        MDL_MH__EXP__2 = MDL_MH**2
 
40
 
 
41
        MDL_NCOL__EXP__2 = MDL_NCOL**2
 
42
 
 
43
        MDL_MB__EXP__2 = MDL_MB**2
 
44
 
 
45
        MDL_MT__EXP__2 = MDL_MT**2
 
46
 
 
47
        MDL_AEW = 1.000000D+00/AEWM1
 
48
 
 
49
        MDL_MW = SQRT(DCMPLX(MDL_MZ__EXP__2/2.000000D+00
 
50
     $   +SQRT(DCMPLX(MDL_MZ__EXP__4/4.000000D+00-(MDL_AEW*PI
 
51
     $   *MDL_MZ__EXP__2)/(MDL_GF*MDL_SQRT__2)))))
 
52
 
 
53
        MDL_SQRT__AEW = SQRT(DCMPLX(MDL_AEW))
 
54
 
 
55
        MDL_EE = 2.000000D+00*MDL_SQRT__AEW*SQRT(DCMPLX(PI))
 
56
 
 
57
        MDL_MW__EXP__2 = MDL_MW**2
 
58
 
 
59
        MDL_SW2 = 1.000000D+00-MDL_MW__EXP__2/MDL_MZ__EXP__2
 
60
 
 
61
        MDL_CW = SQRT(DCMPLX(1.000000D+00-MDL_SW2))
 
62
 
 
63
        MDL_SQRT__SW2 = SQRT(DCMPLX(MDL_SW2))
 
64
 
 
65
        MDL_SW = MDL_SQRT__SW2
 
66
 
 
67
        MDL_G1 = MDL_EE/MDL_CW
 
68
 
 
69
        MDL_GW = MDL_EE/MDL_SW
 
70
 
 
71
        MDL_V = (2.000000D+00*MDL_MW*MDL_SW)/MDL_EE
 
72
 
 
73
        MDL_V__EXP__2 = MDL_V**2
 
74
 
 
75
        MDL_LAM = MDL_MH__EXP__2/(2.000000D+00*MDL_V__EXP__2)
 
76
 
 
77
        MDL_YB = (MDL_YMB*MDL_SQRT__2)/MDL_V
 
78
 
 
79
        MDL_YT = (MDL_YMT*MDL_SQRT__2)/MDL_V
 
80
 
 
81
        MDL_YTAU = (MDL_YMTAU*MDL_SQRT__2)/MDL_V
 
82
 
 
83
        MDL_MUH = SQRT(DCMPLX(MDL_LAM*MDL_V__EXP__2))
 
84
 
 
85
        MDL_AXIALZUP = (3.000000D+00/2.000000D+00)*(-(MDL_EE*MDL_SW)
 
86
     $   /(6.000000D+00*MDL_CW))-(1.000000D+00/2.000000D+00)*((MDL_CW
 
87
     $   *MDL_EE)/(2.000000D+00*MDL_SW))
 
88
 
 
89
        MDL_AXIALZDOWN = (-1.000000D+00/2.000000D+00)*(-(MDL_CW*MDL_EE)
 
90
     $   /(2.000000D+00*MDL_SW))+(-3.000000D+00/2.000000D+00)*(
 
91
     $   -(MDL_EE*MDL_SW)/(6.000000D+00*MDL_CW))
 
92
 
 
93
        MDL_VECTORZUP = (1.000000D+00/2.000000D+00)*((MDL_CW*MDL_EE)
 
94
     $   /(2.000000D+00*MDL_SW))+(5.000000D+00/2.000000D+00)*(-(MDL_EE
 
95
     $   *MDL_SW)/(6.000000D+00*MDL_CW))
 
96
 
 
97
        MDL_VECTORZDOWN = (1.000000D+00/2.000000D+00)*(-(MDL_CW*MDL_EE)
 
98
     $   /(2.000000D+00*MDL_SW))+(-1.000000D+00/2.000000D+00)*(
 
99
     $   -(MDL_EE*MDL_SW)/(6.000000D+00*MDL_CW))
 
100
 
 
101
        MDL_VECTORAUP = (2.000000D+00*MDL_EE)/3.000000D+00
 
102
 
 
103
        MDL_VECTORADOWN = -(MDL_EE)/3.000000D+00
 
104
 
 
105
        MDL_VECTORWMDXU = (1.000000D+00/2.000000D+00)*((MDL_EE)
 
106
     $   /(MDL_SW*MDL_SQRT__2))
 
107
 
 
108
        MDL_AXIALWMDXU = (-1.000000D+00/2.000000D+00)*((MDL_EE)
 
109
     $   /(MDL_SW*MDL_SQRT__2))
 
110
 
 
111
        MDL_VECTORWPUXD = (1.000000D+00/2.000000D+00)*((MDL_EE)
 
112
     $   /(MDL_SW*MDL_SQRT__2))
 
113
 
 
114
        MDL_AXIALWPUXD = -(1.000000D+00/2.000000D+00)*((MDL_EE)
 
115
     $   /(MDL_SW*MDL_SQRT__2))
 
116
 
 
117
        MDL_I1X33 = MDL_YB*MDL_CONJG__CKM3X3
 
118
 
 
119
        MDL_I2X33 = MDL_YT*MDL_CONJG__CKM3X3
 
120
 
 
121
        MDL_I3X33 = MDL_CKM3X3*MDL_YT
 
122
 
 
123
        MDL_I4X33 = MDL_CKM3X3*MDL_YB
 
124
 
 
125
        MDL_VECTOR_TBGP = MDL_I1X33-MDL_I2X33
 
126
 
 
127
        MDL_AXIAL_TBGP = -MDL_I2X33-MDL_I1X33
 
128
 
 
129
        MDL_VECTOR_TBGM = MDL_I3X33-MDL_I4X33
 
130
 
 
131
        MDL_AXIAL_TBGM = -MDL_I4X33-MDL_I3X33
 
132
 
 
133
        MDL_GW__EXP__2 = MDL_GW**2
 
134
 
 
135
        MDL_CW__EXP__2 = MDL_CW**2
 
136
 
 
137
        MDL_EE__EXP__2 = MDL_EE**2
 
138
 
 
139
        MDL_SW__EXP__2 = MDL_SW**2
 
140
 
 
141
        MDL_YB__EXP__2 = MDL_YB**2
 
142
 
 
143
        MDL_YT__EXP__2 = MDL_YT**2
 
144
 
 
145
      ENDIF
 
146
C     
 
147
C     Parameters that should be recomputed at an event by even basis.
 
148
C     
 
149
      AS = G**2/4/PI
 
150
 
 
151
      MDL_SQRT__AS = SQRT(DCMPLX(AS))
 
152
 
 
153
      MDL_G__EXP__4 = G**4
 
154
 
 
155
      MDL_G__EXP__2 = G**2
 
156
 
 
157
      MDL_G__EXP__3 = G**3
 
158
 
 
159
      MDL_MU_R__EXP__2 = MU_R**2
 
160
 
 
161
C     
 
162
C     Parameters that should be updated for the loops.
 
163
C     
 
164
      IF (UPDATELOOP) THEN
 
165
        MDL_R2MIXEDFACTOR_FIN_ = -(MDL_G__EXP__2*(1.000000D+00+MDL_LHV)
 
166
     $   *(MDL_NCOL__EXP__2-1.000000D+00))/(2.000000D+00*MDL_NCOL
 
167
     $   *1.600000D+01*PI**2)
 
168
 
 
169
        MDL_G_UVG_1EPS_ = -((MDL_G__EXP__2)/(2.000000D+00*4.800000D+01
 
170
     $   *PI**2))*1.100000D+01*MDL_CA
 
171
 
 
172
        MDL_G_UVB_1EPS_ = ((MDL_G__EXP__2)/(2.000000D+00*4.800000D+01
 
173
     $   *PI**2))*4.000000D+00*MDL_TF
 
174
 
 
175
        MDL_BMASS_UV_1EPS_ = COND(DCMPLX(MDL_MB),DCMPLX(0.000000D+00)
 
176
     $   ,DCMPLX(MDL_COMPLEXI*((MDL_G__EXP__2)/(1.600000D+01*PI**2))
 
177
     $   *(3.000000D+00*MDL_CF)*MDL_MB))
 
178
 
 
179
        MDL_TMASS_UV_1EPS_ = COND(DCMPLX(MDL_MT),DCMPLX(0.000000D+00)
 
180
     $   ,DCMPLX(MDL_COMPLEXI*((MDL_G__EXP__2)/(1.600000D+01*PI**2))
 
181
     $   *3.000000D+00*MDL_CF*MDL_MT))
 
182
 
 
183
        MDL_G_UVB_FIN_ = COND(DCMPLX(MDL_MB),DCMPLX(0.000000D+00)
 
184
     $   ,DCMPLX(-((MDL_G__EXP__2)/(2.000000D+00*4.800000D+01*PI**2))
 
185
     $   *4.000000D+00*MDL_TF*REGLOG(DCMPLX(MDL_MB__EXP__2
 
186
     $   /MDL_MU_R__EXP__2))))
 
187
 
 
188
        MDL_G_UVT_FIN_ = COND(DCMPLX(MDL_MT),DCMPLX(0.000000D+00)
 
189
     $   ,DCMPLX(-((MDL_G__EXP__2)/(2.000000D+00*4.800000D+01*PI**2))
 
190
     $   *4.000000D+00*MDL_TF*REGLOG(DCMPLX(MDL_MT__EXP__2
 
191
     $   /MDL_MU_R__EXP__2))))
 
192
 
 
193
        MDL_BMASS_UV_FIN_ = COND(DCMPLX(MDL_MB),DCMPLX(0.000000D+00)
 
194
     $   ,DCMPLX(MDL_COMPLEXI*((MDL_G__EXP__2)/(1.600000D+01*PI**2))
 
195
     $   *MDL_CF*(4.000000D+00-3.000000D+00
 
196
     $   *REGLOG(DCMPLX(MDL_MB__EXP__2/MDL_MU_R__EXP__2)))*MDL_MB))
 
197
 
 
198
        MDL_TMASS_UV_FIN_ = COND(DCMPLX(MDL_MT),DCMPLX(0.000000D+00)
 
199
     $   ,DCMPLX(MDL_COMPLEXI*((MDL_G__EXP__2)/(1.600000D+01*PI**2))
 
200
     $   *MDL_CF*(4.000000D+00-3.000000D+00
 
201
     $   *REGLOG(DCMPLX(MDL_MT__EXP__2/MDL_MU_R__EXP__2)))*MDL_MT))
 
202
 
 
203
      ENDIF
 
204
C     
 
205
C     Definition of the EW coupling used in the write out of aqed
 
206
C     
 
207
      GAL(1) = 3.5449077018110318D0 / DSQRT(ABS(AEWM1))
 
208
      GAL(2) = 1D0
 
209