/
ffs1.in
618 lines (602 loc) · 85.2 KB
/
ffs1.in
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
// real part of fermionic one-loop Z self-energy for k^2=s
double rsz1fs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return((AL*s*real(B0fin(MZs,0,0))*(-98*MWs*MZs + 76*sqr(MWs) + 49*sqr(MZs)))/(36.*MWs*(MWs - MZs)*Pi) - (AL*s*real(B0fin(s,0,0))*(-98*MWs*MZs + 76*sqr(MWs) + 49*sqr(MZs)))/(36.*MWs*(MWs - MZs)*Pi) + (AL*(MZs - s)*A0fin(MBs)*(MBs*(-8*MWs*MZs + 16*sqr(MWs) - 17*sqr(MZs)) + MZs*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))))/(36.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*Pi) + (AL*real(B0fin(s,MBs,MBs))*(s*(4*MWs*MZs - 8*sqr(MWs) - 5*sqr(MZs)) + MBs*(8*MWs*MZs - 16*sqr(MWs) + 17*sqr(MZs))))/(72.*MWs*(MWs - MZs)*Pi) + (AL*(MZs - s)*A0fin(MTs)*(MTs*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs)) + MZs*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))))/(36.*MWs*(MWs - MZs)*MZs*(-4*MTs + MZs)*Pi) - (AL*real(B0fin(s,MTs,MTs))*(MTs*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs)) + s*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))))/(72.*MWs*(MWs - MZs)*Pi) + (AL*(MZs - s)*(-2*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs - 16*sqr(MWs) + 17*sqr(MZs)) + MZs*(MTs*MZs*(800*MWs*MZs - 640*sqr(MWs) - 439*sqr(MZs)) + 24*sqr(MZs)*(-10*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)) - 2*sqr(MTs)*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs))) + MBs*((944*MWs*MZs - 736*sqr(MWs) - 487*sqr(MZs))*sqr(MZs) + 8*sqr(MTs)*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs)) + 8*MTs*MZs*(-392*MWs*MZs + 304*sqr(MWs) + 223*sqr(MZs)))))/(72.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*(-4*MTs + MZs)*Pi) - (AL*real(B0fin(MZs,MBs,MBs))*(2*(3*MZs - s)*sqr(MBs)*(-8*MWs*MZs + 16*sqr(MWs) - 17*sqr(MZs)) + s*(4*MWs*MZs - 8*sqr(MWs) - 5*sqr(MZs))*sqr(MZs) + MBs*MZs*(-8*MWs*MZs*s + 16*s*sqr(MWs) + 10*s*sqr(MZs) + 27*trip(MZs))))/(72.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*Pi) - (AL*real(B0fin(MZs,MTs,MTs))*(s*(40*MWs*MZs - 32*sqr(MWs) - 17*sqr(MZs))*sqr(MZs) + 2*(3*MZs - s)*sqr(MTs)*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs)) + MTs*MZs*(-80*MWs*MZs*s + 64*s*sqr(MWs) + 34*s*sqr(MZs) + 27*trip(MZs))))/(72.*MWs*(MWs - MZs)*MZs*(-4*MTs + MZs)*Pi));
}
// imag. part of fermionic one-loop Z self-energy for k^2=s
double isz1fs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return(-(AL*s*imag(B0fin(s,0,0))*(-98*MWs*MZs + 76*sqr(MWs) + 49*sqr(MZs)))/(36.*MWs*(MWs - MZs)*Pi) - (AL*imag(B0fin(s,MBs,MBs))*(MBs*(-8*MWs*MZs + 16*sqr(MWs) - 17*sqr(MZs)) + s*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))))/(72.*MWs*(MWs - MZs)*Pi) - (AL*imag(B0fin(s,MTs,MTs))*(MTs*(-80*MWs*MZs + 64*sqr(MWs) + 7*sqr(MZs)) + s*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))))/(72.*MWs*(MWs - MZs)*Pi));
}
// real part of bosonic one-loop Z self-energy for k^2=s
double rsz1bs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return((AL*(MZs - s)*A0fin(MWs)*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(24.*MWs*(MWs - MZs)*MZs*Pi) - (AL*(MZs - s)*A0fin(MZs)*(MHs*MZs*(-5*MZs + 9*s) + (MZs - 2*s)*sqr(MHs) + 4*(MZs - 4*s)*sqr(MZs)))/(48.*MWs*(MHs - 4*MZs)*(MWs - MZs)*Pi*s) - (AL*real(B0fin(s,MHs,MZs))*sqr(MZs)*(10*MZs*s - 2*MHs*(MZs + s) + sqr(MHs) + sqr(MZs) + sqr(s)))/(48.*MWs*(MWs - MZs)*Pi*s) + (AL*real(B0fin(s,MWs,MWs))*(4*MWs*MZs*(-5*MZs + s) + 4*(8*MZs + 9*s)*sqr(MWs) - s*sqr(MZs) + 48*trip(MWs)))/(48.*MWs*(MWs - MZs)*Pi) + (AL*real(B0fin(MZs,MWs,MWs))*(4*MZs*(-18*MZs + s)*sqr(MWs) + 2*MWs*(9*MZs - s)*sqr(MZs) + 24*(-3*MZs + s)*trip(MWs) + s*trip(MZs)))/(48.*MWs*(MWs - MZs)*MZs*Pi) + (AL*A0fin(MHs)*(MZs*sqr(MHs)*(16*MZs*s - 5*sqr(MZs) - 11*sqr(s)) + 4*MHs*sqr(MZs)*(-7*MZs*s + sqr(MZs) + 6*sqr(s)) + (-3*MZs*s + sqr(MZs) + 2*sqr(s))*trip(MHs) + 24*(MZs - s)*s*trip(MZs)))/(48.*MHs*MWs*(MHs - 4*MZs)*(MWs - MZs)*Pi*s) + (AL*real(B0fin(MZs,MHs,MZs))*(MZs*(-21*MZs + 13*s)*sqr(MHs) + 4*MHs*(15*MZs - 8*s)*sqr(MZs) + (3*MZs - 2*s)*trip(MHs) + 12*(-7*MZs + 3*s)*trip(MZs)))/(48.*MWs*(MHs - 4*MZs)*(MWs - MZs)*Pi) + (AL*(MZs - s)*(sqr(MHs)*sqr(MZs) + MHs*(28*MZs*sqr(MWs) - 18*MWs*sqr(MZs) + 24*trip(MWs) - 5*trip(MZs)) + 8*MZs*(-14*MZs*sqr(MWs) + 9*MWs*sqr(MZs) - 12*trip(MWs) + 2*trip(MZs))))/(48.*MWs*(MWs - MZs)*MZs*(-MHs + 4*MZs)*Pi));
}
// imag. part of bosonic one-loop Z self-energy for k^2=s
double isz1bs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return(-(AL*imag(B0fin(s,MHs,MZs))*sqr(MZs)*(10*MZs*s - 2*MHs*(MZs + s) + sqr(MHs) + sqr(MZs) + sqr(s)))/(48.*MWs*(MWs - MZs)*Pi*s) + (AL*imag(B0fin(s,MWs,MWs))*(4*MWs*MZs*(-5*MZs + s) + 4*(8*MZs + 9*s)*sqr(MWs) - s*sqr(MZs) + 48*trip(MWs)))/(48.*MWs*(MWs - MZs)*Pi));
}
// real part of fermionic one-loop photon self-energy for k^2=s
double rsg1fs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
double deltaAlpha = input.get(Delal);
return(-(deltaAlpha*s) + (2*AL*(3*MBs*(MZs - s) + 4*MZs*(3*MTs + s)))/(27.*MZs*Pi) - (2*AL*(MZs - s)*A0fin(MBs))/(9.*MZs*Pi) - (4*AL*(2*MTs + s)*A0fin(MTs))/(9.*MTs*Pi) - (19*AL*s*real(B0fin(MZs,0,0)))/(9.*Pi) - (AL*(2*MBs + MZs)*s*real(B0fin(MZs,MBs,MBs)))/(9.*MZs*Pi) + (19*AL*s*real(B0fin(s,0,0)))/(9.*Pi) + (AL*(2*MBs + s)*real(B0fin(s,MBs,MBs)))/(9.*Pi) + (4*AL*(2*MTs + s)*real(B0fin(s,MTs,MTs)))/(9.*Pi));
}
// imag. part of fermionic one-loop photon self-energy for k^2=s
double isg1fs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return((19*AL*s*imag(B0fin(s,0,0)))/(9.*Pi) + (AL*(2*MBs + s)*imag(B0fin(s,MBs,MBs)))/(9.*Pi) + (4*AL*(2*MTs + s)*imag(B0fin(s,MTs,MTs)))/(9.*Pi));
}
// real part of bosonic one-loop photon self-energy for k^2=s
double rsg1bs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return(-(AL*(12*MWs + 7*s))/(12.*Pi) + (AL*(4*MWs + 3*s)*A0fin(MWs))/(4.*MWs*Pi) - (AL*(4*MWs + 3*s)*real(B0fin(s,MWs,MWs)))/(4.*Pi));
}
// imag. part of bosonic one-loop photon self-energy for k^2=s
double isg1bs(double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al);
return(-(AL*(4*MWs + 3*s)*imag(B0fin(s,MWs,MWs)))/(4.*Pi));
}
// real part of fermionic one-loop Z vertex factor for k^2=s
double rz1fs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
deltaAlpha = input.get(Delal),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = (deltaAlpha*(-4*MWs + 3*MZs)*Pi)/(32.*AL) + ((76*MWs - 49*MZs)*real(B0fin(s,0,0)))/144. + ((4*MWs - MZs)*(2*MBs + s)*real(B0fin(s,MBs,MBs)))/(144.*s) + ((8*MWs - 5*MZs)*(2*MTs + s)*real(B0fin(s,MTs,MTs)))/(72.*s) + (3*(2*MWs - 3*MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) - ((2*MWs - 3*MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,0,0))*(1608*MZs*sqr(MWs) - 1446*MWs*sqr(MZs) - 608*trip(MWs) + 473*trip(MZs)))/(1152.*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MWs*sqr(MZs)*(144*MZs*sqr(MWs) - 120*MWs*sqr(MZs) - 64*trip(MWs) + 49*trip(MZs)) + 2*sqr(MBs)*(64*powint(MWs,4) + 153*powint(MZs,4) + 156*sqr(MWs)*sqr(MZs) - 144*MZs*trip(MWs) - 211*MWs*trip(MZs)) + MBs*MZs*(128*powint(MWs,4) - 81*powint(MZs,4) + 240*sqr(MWs)*sqr(MZs) - 288*MZs*trip(MWs) - 44*MWs*trip(MZs))))/(2304.*MWs*MZs*(-4*MBs + MZs)*sqr(MWs - MZs)) - (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-320*MZs*sqr(MWs) + 268*MWs*sqr(MZs) + 128*trip(MWs) - 85*trip(MZs)) + 2*sqr(MTs)*(256*powint(MWs,4) + 63*powint(MZs,4) + 1484*sqr(MWs)*sqr(MZs) - 1024*MZs*trip(MWs) - 797*MWs*trip(MZs)) + MTs*MZs*(-256*powint(MWs,4) + 81*powint(MZs,4) - 536*sqr(MWs)*sqr(MZs) + 640*MZs*trip(MWs) + 116*MWs*trip(MZs))))/(2304.*MWs*MZs*(-4*MTs + MZs)*sqr(MWs - MZs)) - (A0fin(MBs)*(36*(2*MWs - 3*MZs)*s*sqr(MBs)*trip(MZs) + MBs*(288*MZs*(4*MZs - s)*powint(MWs,4) - 128*(4*MZs - s)*powint(MWs,5) + 27*(4*MTs + MZs)*s*powint(MZs,4) + 24*(-32*MZs + 13*s)*sqr(MZs)*trip(MWs) - 36*MWs*(2*MTs - 5*MZs)*s*trip(MZs) + 2*(64*MZs - 175*s)*sqr(MWs)*trip(MZs)) + MZs*(-288*MZs*(MZs - s)*powint(MWs,4) + 128*(MZs - s)*powint(MWs,5) - 27*MTs*s*powint(MZs,4) + 48*(4*MZs - 5*s)*sqr(MZs)*trip(MWs) + 18*MWs*(MTs - 3*MZs)*s*trip(MZs) + 2*(-16*MZs + 67*s)*sqr(MWs)*trip(MZs))))/(2304.*MZs*(-4*MBs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) - (A0fin(MTs)*(32*(4*MWs - 3*MZs)*s*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs - 3*MZs)*s*trip(MTs)*trip(MZs) + MTs*MZs*(-64*MZs*(21*MZs - 32*s)*powint(MWs,4) + 256*(2*MZs - 3*s)*powint(MWs,5) - 27*MBs*s*powint(MZs,4) + 8*(144*MZs - 227*s)*sqr(MZs)*trip(MWs) + 18*MWs*(MBs - 3*MZs)*s*trip(MZs) + 10*(-32*MZs + 59*s)*sqr(MWs)*trip(MZs)) + sqr(MTs)*(256*MZs*(21*MZs + 8*s)*powint(MWs,4) - 512*(4*MZs + s)*powint(MWs,5) + 27*(4*MBs + MZs)*s*powint(MZs,4) - 8*(576*MZs + 371*s)*sqr(MZs)*trip(MWs) - 36*MWs*(2*MBs + 7*MZs)*s*trip(MZs) + 2*(640*MZs + 833*s)*sqr(MWs)*trip(MZs))))/(2304.*MTs*MZs*(-4*MTs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (6*(4*MTs - MZs)*sqr(MBs)*(-80*MZs*(4*MZs - s)*sqr(MWs) + 4*MWs*(16*MZs - 19*s)*sqr(MZs) + 64*(4*MZs - s)*trip(MWs) + 99*s*trip(MZs)) + MZs*(8*s*sqr(MZs)*(-632*MZs*sqr(MWs) + 474*MWs*sqr(MZs) + 320*trip(MWs) - 135*trip(MZs)) - 6*sqr(MTs)*(-128*MZs*(13*MZs + 6*s)*sqr(MWs) + 4*MWs*(160*MZs + 179*s)*sqr(MZs) + 256*(4*MZs + s)*trip(MWs) - 117*s*trip(MZs)) + MTs*MZs*(-64*MZs*(39*MZs - 262*s)*sqr(MWs) + 12*MWs*(80*MZs - 1047*s)*sqr(MZs) + 512*(3*MZs - 17*s)*trip(MWs) + 3807*s*trip(MZs))) + MBs*(24*sqr(MTs)*(-128*MZs*(13*MZs + 6*s)*sqr(MWs) + 4*MWs*(160*MZs + 179*s)*sqr(MZs) + 256*(4*MZs + s)*trip(MWs) - 117*s*trip(MZs)) - 8*MTs*MZs*(8*MZs*(-186*MZs + 1033*s)*sqr(MWs) + 12*MWs*(44*MZs - 523*s)*sqr(MZs) + 32*(30*MZs - 133*s)*trip(MWs) + 1971*s*trip(MZs)) + sqr(MZs)*(16*MZs*(-30*MZs + 1249*s)*sqr(MWs) + 12*MWs*(8*MZs - 1263*s)*sqr(MZs) + 64*(6*MZs - 157*s)*trip(MWs) + 4455*s*trip(MZs))))/(6912.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*(-4*MTs + MZs)*s);
break;
case NUE:
case NUM:
case NUT:
res = (deltaAlpha*MZs*Pi)/(32.*AL) + (3*(2*MWs - MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + (3*(-2*MWs + MZs)*real(B0fin(MZs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) - ((2*MWs - MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) - (18*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs + 16*sqr(MWs) - 11*sqr(MZs)) + MZs*(MTs*MZs*(2368*MWs*MZs - 1792*sqr(MWs) - 1269*sqr(MZs)) - 6*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*sqr(MZs)*(-98*MWs*MZs + 80*sqr(MWs) + 45*sqr(MZs))) + MBs*((3136*MWs*MZs - 2416*sqr(MWs) - 1485*sqr(MZs))*sqr(MZs) + 24*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*MTs*MZs*(-1184*MWs*MZs + 824*sqr(MWs) + 657*sqr(MZs))))/(6912.*MWs*(MWs - MZs)*(-4*MBs + MZs)*(-4*MTs + MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*(40*MWs*MZs - 32*sqr(MWs) - 17*sqr(MZs))*sqr(MZs) + MTs*MZs*(-80*MZs*sqr(MWs) + 88*MWs*sqr(MZs) + 64*trip(MWs) - 27*trip(MZs)) + 2*sqr(MTs)*(-592*MZs*sqr(MWs) + 275*MWs*sqr(MZs) + 320*trip(MWs) - 21*trip(MZs))))/(2304.*MWs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MBs*(16*MWs*MZs + 8*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) + MWs*(-4*MWs + MZs)*trip(MZs) + 2*sqr(MBs)*(-8*MZs*sqr(MWs) - 31*MWs*sqr(MZs) + 16*trip(MWs) + 17*trip(MZs))))/(768.*MWs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(-2*MWs + MZs)*sqr(MZs)*trip(MTs) + MTs*MZs*(-64*powint(MWs,4) - 18*MWs*(MBs - MZs)*sqr(MZs) - 130*sqr(MWs)*sqr(MZs) + 176*MZs*trip(MWs) + 9*MBs*trip(MZs)) + sqr(MTs)*(-640*powint(MWs,4) + 24*MWs*(3*MBs + 4*MZs)*sqr(MZs) - 622*sqr(MWs)*sqr(MZs) + 1184*MZs*trip(MWs) - 9*(4*MBs + MZs)*trip(MZs))))/(2304.*MTs*(-4*MTs + MZs)*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(12*(-2*MWs + MZs)*sqr(MBs)*sqr(MZs) + sqr(MZs)*(6*MWs*MZs*(-MTs + MZs) - 14*MZs*sqr(MWs) + 3*MTs*sqr(MZs) + 8*trip(MWs)) + MBs*(-32*powint(MWs,4) + 8*MWs*(3*MTs - 2*MZs)*sqr(MZs) + 38*sqr(MWs)*sqr(MZs) + 16*MZs*trip(MWs) - 3*(4*MTs + MZs)*trip(MZs))))/(768.*(-4*MBs + MZs)*sqr(MWs)*sqr(MWs - MZs));
break;
case UQU:
case CQU:
res = (deltaAlpha*(8*MWs - 5*MZs)*Pi)/(96.*AL) + ((-76*MWs + 49*MZs)*real(B0fin(s,0,0)))/216. - ((4*MWs - MZs)*(2*MBs + s)*real(B0fin(s,MBs,MBs)))/(216.*s) - ((8*MWs - 5*MZs)*(2*MTs + s)*real(B0fin(s,MTs,MTs)))/(108.*s) + ((-2*MWs + 5*MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + ((2*MWs - 5*MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(768.*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,0,0))*(-3216*MZs*sqr(MWs) + 2838*MWs*sqr(MZs) + 1216*trip(MWs) - 919*trip(MZs)))/(3456.*sqr(MWs - MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-672*MZs*sqr(MWs) + 576*MWs*sqr(MZs) + 256*trip(MWs) - 187*trip(MZs)) + 2*sqr(MTs)*(512*powint(MWs,4) + 105*powint(MZs,4) + 2376*sqr(MWs)*sqr(MZs) - 1728*MZs*trip(MWs) - 1319*MWs*trip(MZs)) + MTs*MZs*(-512*powint(MWs,4) + 135*powint(MZs,4) - 1152*sqr(MWs)*sqr(MZs) + 1344*MZs*trip(MWs) + 320*MWs*trip(MZs))))/(6912.*MWs*MZs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MWs*sqr(MZs)*(-288*MZs*sqr(MWs) + 228*MWs*sqr(MZs) + 128*trip(MWs) - 95*trip(MZs)) + MBs*MZs*(-256*powint(MWs,4) + 135*powint(MZs,4) - 456*sqr(MWs)*sqr(MZs) + 576*MZs*trip(MWs) + 136*MWs*trip(MZs)) + sqr(MBs)*(-256*powint(MWs,4) - 510*powint(MZs,4) - 672*sqr(MWs)*sqr(MZs) + 672*MZs*trip(MWs) + 658*MWs*trip(MZs))))/(6912.*MWs*MZs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(32*(8*MWs - 5*MZs)*s*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs - 5*MZs)*s*trip(MTs)*trip(MZs) + MTs*MZs*(-1344*MZs*(2*MZs - 3*s)*powint(MWs,4) + 512*(2*MZs - 3*s)*powint(MWs,5) - 45*MBs*s*powint(MZs,4) + 1152*(2*MZs - 3*s)*sqr(MZs)*trip(MWs) + 18*MWs*(MBs - 5*MZs)*s*trip(MZs) + 10*(-64*MZs + 105*s)*sqr(MWs)*trip(MZs)) + sqr(MTs)*(384*MZs*(28*MZs + 9*s)*powint(MWs,4) - 1024*(4*MZs + s)*powint(MWs,5) + 45*(4*MBs + MZs)*s*powint(MZs,4) - 144*(64*MZs + 33*s)*sqr(MZs)*trip(MWs) - 24*MWs*(3*MBs + 17*MZs)*s*trip(MZs) + 10*(256*MZs + 271*s)*sqr(MWs)*trip(MZs))))/(6912.*MTs*MZs*(-4*MTs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(36*(2*MWs - 5*MZs)*s*sqr(MBs)*trip(MZs) + MBs*(96*MZs*(24*MZs - 7*s)*powint(MWs,4) - 256*(4*MZs - s)*powint(MWs,5) + 45*(4*MTs + MZs)*s*powint(MZs,4) - 96*(16*MZs - 7*s)*sqr(MZs)*trip(MWs) + 24*MWs*(-3*MTs + 13*MZs)*s*trip(MZs) + 2*(128*MZs - 293*s)*sqr(MWs)*trip(MZs)) + MZs*(-576*MZs*(MZs - s)*powint(MWs,4) + 256*(MZs - s)*powint(MWs,5) - 45*MTs*s*powint(MZs,4) + 24*(16*MZs - 19*s)*sqr(MZs)*trip(MWs) + 18*MWs*(MTs - 5*MZs)*s*trip(MZs) + sqr(MWs)*(-64*powint(MZs,4) + 226*s*trip(MZs)))))/(6912.*MZs*(-4*MBs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (-6*(4*MTs - MZs)*sqr(MBs)*(16*MZs*(-40*MZs + 13*s)*sqr(MWs) + 128*MWs*(MZs - s)*sqr(MZs) + 128*(4*MZs - s)*trip(MWs) + 165*s*trip(MZs)) + MZs*(8*s*sqr(MZs)*(1184*MZs*sqr(MWs) - 850*MWs*sqr(MZs) - 640*trip(MWs) + 225*trip(MZs)) + MTs*MZs*(128*MZs*(39*MZs - 248*s)*sqr(MWs) - 40*MWs*(48*MZs - 569*s)*sqr(MZs) - 1024*(3*MZs - 17*s)*trip(MWs) - 6345*s*trip(MZs)) + 6*sqr(MTs)*(-64*MZs*(52*MZs + 19*s)*sqr(MWs) + 40*MWs*(32*MZs + 29*s)*sqr(MZs) + 512*(4*MZs + s)*trip(MWs) - 195*s*trip(MZs))) + MBs*(sqr(MZs)*(16*MZs*(60*MZs - 2347*s)*sqr(MWs) - 8*MWs*(24*MZs - 3397*s)*sqr(MZs) - 128*(6*MZs - 157*s)*trip(MWs) - 7425*s*trip(MZs)) - 24*sqr(MTs)*(-64*MZs*(52*MZs + 19*s)*sqr(MWs) + 40*MWs*(32*MZs + 29*s)*sqr(MZs) + 512*(4*MZs + s)*trip(MWs) - 195*s*trip(MZs)) + 8*MTs*MZs*(-8*MZs*(372*MZs - 1963*s)*sqr(MWs) + 8*MWs*(132*MZs - 1421*s)*sqr(MZs) + 64*(30*MZs - 133*s)*trip(MWs) + 3285*s*trip(MZs))))/(20736.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*(-4*MTs + MZs)*s);
break;
case DQU:
case SQU:
res = (deltaAlpha*(-4*MWs + MZs)*Pi)/(96.*AL) + ((76*MWs - 49*MZs)*real(B0fin(s,0,0)))/432. + ((4*MWs - MZs)*(2*MBs + s)*real(B0fin(s,MBs,MBs)))/(432.*s) + ((8*MWs - 5*MZs)*(2*MTs + s)*real(B0fin(s,MTs,MTs)))/(216.*s) - ((2*MWs + MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + ((2*MWs + MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(768.*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,0,0))*(1608*MZs*sqr(MWs) - 1338*MWs*sqr(MZs) - 608*trip(MWs) + 419*trip(MZs)))/(3456.*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MWs*sqr(MZs)*(144*MZs*sqr(MWs) - 96*MWs*sqr(MZs) - 64*trip(MWs) + 43*trip(MZs)) + MBs*MZs*(128*powint(MWs,4) - 27*powint(MZs,4) + 192*sqr(MWs)*sqr(MZs) - 288*MZs*trip(MWs) - 140*MWs*trip(MZs)) + 2*sqr(MBs)*(64*powint(MWs,4) + 51*powint(MZs,4) + 204*sqr(MWs)*sqr(MZs) - 240*MZs*trip(MWs) - 25*MWs*trip(MZs))))/(6912.*MWs*MZs*(-4*MBs + MZs)*sqr(MWs - MZs)) - (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-384*MZs*sqr(MWs) + 348*MWs*sqr(MZs) + 128*trip(MWs) - 119*trip(MZs)) + sqr(MTs)*(512*powint(MWs,4) + 42*powint(MZs,4) + 600*sqr(MWs)*sqr(MZs) - 768*MZs*trip(MWs) - 494*MWs*trip(MZs)) + MTs*MZs*(-256*powint(MWs,4) + 27*powint(MZs,4) - 696*sqr(MWs)*sqr(MZs) + 768*MZs*trip(MWs) + 292*MWs*trip(MZs))))/(6912.*MWs*MZs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (A0fin(MBs)*(36*(2*MWs + MZs)*s*sqr(MBs)*trip(MZs) + MBs*(-96*MZs*(12*MZs - 5*s)*powint(MWs,4) + 128*(4*MZs - s)*powint(MWs,5) - 9*(4*MTs + MZs)*s*powint(MZs,4) + 24*(32*MZs - 17*s)*sqr(MZs)*trip(MWs) - 12*MWs*(6*MTs + 7*MZs)*s*trip(MZs) - 2*(64*MZs - 61*s)*sqr(MWs)*trip(MZs)) + MZs*(288*MZs*(MZs - s)*powint(MWs,4) - 128*(MZs - s)*powint(MWs,5) + 9*MTs*s*powint(MZs,4) - 192*(MZs - s)*sqr(MZs)*trip(MWs) + 18*MWs*(MTs + MZs)*s*trip(MZs) + 2*(16*MZs - 25*s)*sqr(MWs)*trip(MZs))))/(6912.*MZs*(-4*MBs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MTs)*(-32*(4*MWs - MZs)*s*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs + MZs)*s*trip(MTs)*trip(MZs) + MTs*MZs*(192*MZs*(7*MZs - 10*s)*powint(MWs,4) - 256*(2*MZs - 3*s)*powint(MWs,5) + 9*MBs*s*powint(MZs,4) - 24*(48*MZs - 61*s)*sqr(MZs)*trip(MWs) + 18*MWs*(MBs + MZs)*s*trip(MZs) + 10*(32*MZs - 33*s)*sqr(MWs)*trip(MZs)) + sqr(MTs)*(-768*MZs*(7*MZs + s)*powint(MWs,4) + 512*(4*MZs + s)*powint(MWs,5) - 9*(4*MBs + MZs)*s*powint(MZs,4) + 24*(192*MZs + 25*s)*sqr(MZs)*trip(MWs) + 12*MWs*(-6*MBs + 5*MZs)*s*trip(MZs) - 2*(640*MZs + 211*s)*sqr(MWs)*trip(MZs))))/(6912.*MTs*MZs*(-4*MTs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (6*(4*MTs - MZs)*sqr(MBs)*(16*MZs*(-20*MZs + 11*s)*sqr(MWs) + 4*MWs*(16*MZs - 7*s)*sqr(MZs) + 64*(4*MZs - s)*trip(MWs) + 33*s*trip(MZs)) + MZs*(8*s*sqr(MZs)*(-472*MZs*sqr(MWs) + 278*MWs*sqr(MZs) + 320*trip(MWs) - 45*trip(MZs)) - 6*sqr(MTs)*(-128*MZs*(13*MZs + s)*sqr(MWs) + 4*MWs*(160*MZs + 43*s)*sqr(MZs) + 256*(4*MZs + s)*trip(MWs) - 39*s*trip(MZs)) + MTs*MZs*(-64*MZs*(39*MZs - 206*s)*sqr(MWs) + 4*MWs*(240*MZs - 1957*s)*sqr(MZs) + 512*(3*MZs - 17*s)*trip(MWs) + 1269*s*trip(MZs))) + MBs*(24*sqr(MTs)*(-128*MZs*(13*MZs + s)*sqr(MWs) + 4*MWs*(160*MZs + 43*s)*sqr(MZs) + 256*(4*MZs + s)*trip(MWs) - 39*s*trip(MZs)) - 8*MTs*MZs*(-8*MZs*(186*MZs - 827*s)*sqr(MWs) + 4*MWs*(132*MZs - 977*s)*sqr(MZs) + 32*(30*MZs - 133*s)*trip(MWs) + 657*s*trip(MZs)) + sqr(MZs)*(-16*MZs*(30*MZs - 947*s)*sqr(MWs) + 4*MWs*(24*MZs - 2221*s)*sqr(MZs) + 64*(6*MZs - 157*s)*trip(MWs) + 1485*s*trip(MZs))))/(20736.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*(-4*MTs + MZs)*s);
break;
case BQU:
res = (deltaAlpha*(-4*MWs + MZs)*Pi)/(96.*AL) + ((76*MWs - 49*MZs)*real(B0fin(s,0,0)))/432. + ((4*MWs - MZs)*(2*MBs + s)*real(B0fin(s,MBs,MBs)))/(432.*s) + ((8*MWs - 5*MZs)*(2*MTs + s)*real(B0fin(s,MTs,MTs)))/(216.*s) - ((2*MWs + MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + ((2*MWs + MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(768.*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,0,0))*(1608*MZs*sqr(MWs) - 1338*MWs*sqr(MZs) - 608*trip(MWs) + 419*trip(MZs)))/(3456.*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MWs*sqr(MZs)*(144*MZs*sqr(MWs) - 96*MWs*sqr(MZs) - 64*trip(MWs) + 43*trip(MZs)) + MBs*MZs*(128*powint(MWs,4) - 27*powint(MZs,4) + 192*sqr(MWs)*sqr(MZs) - 288*MZs*trip(MWs) - 140*MWs*trip(MZs)) + 2*sqr(MBs)*(64*powint(MWs,4) + 51*powint(MZs,4) + 204*sqr(MWs)*sqr(MZs) - 240*MZs*trip(MWs) - 25*MWs*trip(MZs))))/(6912.*MWs*MZs*(-4*MBs + MZs)*sqr(MWs - MZs)) - (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-384*MZs*sqr(MWs) + 348*MWs*sqr(MZs) + 128*trip(MWs) - 119*trip(MZs)) + sqr(MTs)*(512*powint(MWs,4) + 42*powint(MZs,4) + 600*sqr(MWs)*sqr(MZs) - 768*MZs*trip(MWs) - 494*MWs*trip(MZs)) + MTs*MZs*(-256*powint(MWs,4) + 27*powint(MZs,4) - 696*sqr(MWs)*sqr(MZs) + 768*MZs*trip(MWs) + 292*MWs*trip(MZs))))/(6912.*MWs*MZs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (A0fin(MBs)*(36*(2*MWs + MZs)*s*sqr(MBs)*trip(MZs) + MBs*(-96*MZs*(12*MZs - 5*s)*powint(MWs,4) + 128*(4*MZs - s)*powint(MWs,5) - 9*(4*MTs + MZs)*s*powint(MZs,4) + 24*(32*MZs - 17*s)*sqr(MZs)*trip(MWs) - 12*MWs*(6*MTs + 7*MZs)*s*trip(MZs) - 2*(64*MZs - 61*s)*sqr(MWs)*trip(MZs)) + MZs*(288*MZs*(MZs - s)*powint(MWs,4) - 128*(MZs - s)*powint(MWs,5) + 9*MTs*s*powint(MZs,4) - 192*(MZs - s)*sqr(MZs)*trip(MWs) + 18*MWs*(MTs + MZs)*s*trip(MZs) + 2*(16*MZs - 25*s)*sqr(MWs)*trip(MZs))))/(6912.*MZs*(-4*MBs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MTs)*(-32*(4*MWs - MZs)*s*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs + MZs)*s*trip(MTs)*trip(MZs) + MTs*MZs*(192*MZs*(7*MZs - 10*s)*powint(MWs,4) - 256*(2*MZs - 3*s)*powint(MWs,5) + 9*MBs*s*powint(MZs,4) - 24*(48*MZs - 61*s)*sqr(MZs)*trip(MWs) + 18*MWs*(MBs + MZs)*s*trip(MZs) + 10*(32*MZs - 33*s)*sqr(MWs)*trip(MZs)) + sqr(MTs)*(-768*MZs*(7*MZs + s)*powint(MWs,4) + 512*(4*MZs + s)*powint(MWs,5) - 9*(4*MBs + MZs)*s*powint(MZs,4) + 24*(192*MZs + 25*s)*sqr(MZs)*trip(MWs) + 12*MWs*(-6*MBs + 5*MZs)*s*trip(MZs) - 2*(640*MZs + 211*s)*sqr(MWs)*trip(MZs))))/(6912.*MTs*MZs*(-4*MTs + MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (6*(4*MTs - MZs)*sqr(MBs)*(16*MZs*(-20*MZs + 11*s)*sqr(MWs) + 4*MWs*(16*MZs - 7*s)*sqr(MZs) + 64*(4*MZs - s)*trip(MWs) + 33*s*trip(MZs)) + MZs*(8*s*sqr(MZs)*(-472*MZs*sqr(MWs) + 278*MWs*sqr(MZs) + 320*trip(MWs) - 45*trip(MZs)) - 6*sqr(MTs)*(-128*MZs*(13*MZs + s)*sqr(MWs) + 4*MWs*(160*MZs + 43*s)*sqr(MZs) + 256*(4*MZs + s)*trip(MWs) - 39*s*trip(MZs)) + MTs*MZs*(-64*MZs*(39*MZs - 206*s)*sqr(MWs) + 4*MWs*(240*MZs - 1957*s)*sqr(MZs) + 512*(3*MZs - 17*s)*trip(MWs) + 1269*s*trip(MZs))) + MBs*(24*sqr(MTs)*(-128*MZs*(13*MZs + s)*sqr(MWs) + 4*MWs*(160*MZs + 43*s)*sqr(MZs) + 256*(4*MZs + s)*trip(MWs) - 39*s*trip(MZs)) - 8*MTs*MZs*(-8*MZs*(186*MZs - 827*s)*sqr(MWs) + 4*MWs*(132*MZs - 977*s)*sqr(MZs) + 32*(30*MZs - 133*s)*trip(MWs) + 657*s*trip(MZs)) + sqr(MZs)*(-16*MZs*(30*MZs - 947*s)*sqr(MWs) + 4*MWs*(24*MZs - 2221*s)*sqr(MZs) + 64*(6*MZs - 157*s)*trip(MWs) + 1485*s*trip(MZs))))/(20736.*MWs*(MWs - MZs)*MZs*(-4*MBs + MZs)*(-4*MTs + MZs)*s);
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = -(deltaAlpha*MZs*Pi)/(32.*AL) + (3*(-2*MWs + MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + (3*(2*MWs - MZs)*real(B0fin(MZs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + ((2*MWs - MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) - (18*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs + 16*sqr(MWs) - 11*sqr(MZs)) + MZs*(MTs*MZs*(2368*MWs*MZs - 1792*sqr(MWs) - 1269*sqr(MZs)) - 6*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*sqr(MZs)*(-98*MWs*MZs + 80*sqr(MWs) + 45*sqr(MZs))) + MBs*((3136*MWs*MZs - 2416*sqr(MWs) - 1485*sqr(MZs))*sqr(MZs) + 24*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*MTs*MZs*(-1184*MWs*MZs + 824*sqr(MWs) + 657*sqr(MZs))))/(6912.*MWs*(-4*MBs + MZs)*(-4*MTs + MZs)*(-MWs + MZs)) - (real(B0fin(MZs,MBs,MBs))*(MBs*(16*MWs*MZs + 8*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) + MWs*(-4*MWs + MZs)*trip(MZs) + 2*sqr(MBs)*(-8*MZs*sqr(MWs) - 31*MWs*sqr(MZs) + 16*trip(MWs) + 17*trip(MZs))))/(768.*MWs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs)) + MTs*MZs*(80*MZs*sqr(MWs) - 88*MWs*sqr(MZs) - 64*trip(MWs) + 27*trip(MZs)) + sqr(MTs)*(1184*MZs*sqr(MWs) - 550*MWs*sqr(MZs) - 640*trip(MWs) + 42*trip(MZs))))/(2304.*MWs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(-32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs - MZs)*sqr(MZs)*trip(MTs) + MTs*MZs*(64*powint(MWs,4) + 18*MWs*(MBs - MZs)*sqr(MZs) + 130*sqr(MWs)*sqr(MZs) - 176*MZs*trip(MWs) - 9*MBs*trip(MZs)) + sqr(MTs)*(640*powint(MWs,4) - 24*MWs*(3*MBs + 4*MZs)*sqr(MZs) + 622*sqr(MWs)*sqr(MZs) - 1184*MZs*trip(MWs) + 9*(4*MBs + MZs)*trip(MZs))))/(2304.*MTs*(-4*MTs + MZs)*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(12*(-2*MWs + MZs)*sqr(MBs)*sqr(MZs) + sqr(MZs)*(6*MWs*MZs*(-MTs + MZs) - 14*MZs*sqr(MWs) + 3*MTs*sqr(MZs) + 8*trip(MWs)) + MBs*(-32*powint(MWs,4) + 8*MWs*(3*MTs - 2*MZs)*sqr(MZs) + 38*sqr(MWs)*sqr(MZs) + 16*MZs*trip(MWs) - 3*(4*MTs + MZs)*trip(MZs))))/(768.*(4*MBs - MZs)*sqr(MWs)*sqr(MWs - MZs));
break;
case NUE:
case NUM:
case NUT:
res = (deltaAlpha*MZs*Pi)/(32.*AL) + (3*(2*MWs - MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + (3*(-2*MWs + MZs)*real(B0fin(MZs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) - ((2*MWs - MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) - (18*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs + 16*sqr(MWs) - 11*sqr(MZs)) + MZs*(MTs*MZs*(2368*MWs*MZs - 1792*sqr(MWs) - 1269*sqr(MZs)) - 6*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*sqr(MZs)*(-98*MWs*MZs + 80*sqr(MWs) + 45*sqr(MZs))) + MBs*((3136*MWs*MZs - 2416*sqr(MWs) - 1485*sqr(MZs))*sqr(MZs) + 24*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*MTs*MZs*(-1184*MWs*MZs + 824*sqr(MWs) + 657*sqr(MZs))))/(6912.*MWs*(MWs - MZs)*(-4*MBs + MZs)*(-4*MTs + MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*(40*MWs*MZs - 32*sqr(MWs) - 17*sqr(MZs))*sqr(MZs) + MTs*MZs*(-80*MZs*sqr(MWs) + 88*MWs*sqr(MZs) + 64*trip(MWs) - 27*trip(MZs)) + 2*sqr(MTs)*(-592*MZs*sqr(MWs) + 275*MWs*sqr(MZs) + 320*trip(MWs) - 21*trip(MZs))))/(2304.*MWs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MBs*(16*MWs*MZs + 8*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) + MWs*(-4*MWs + MZs)*trip(MZs) + 2*sqr(MBs)*(-8*MZs*sqr(MWs) - 31*MWs*sqr(MZs) + 16*trip(MWs) + 17*trip(MZs))))/(768.*MWs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(-2*MWs + MZs)*sqr(MZs)*trip(MTs) + MTs*MZs*(-64*powint(MWs,4) - 18*MWs*(MBs - MZs)*sqr(MZs) - 130*sqr(MWs)*sqr(MZs) + 176*MZs*trip(MWs) + 9*MBs*trip(MZs)) + sqr(MTs)*(-640*powint(MWs,4) + 24*MWs*(3*MBs + 4*MZs)*sqr(MZs) - 622*sqr(MWs)*sqr(MZs) + 1184*MZs*trip(MWs) - 9*(4*MBs + MZs)*trip(MZs))))/(2304.*MTs*(-4*MTs + MZs)*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(12*(-2*MWs + MZs)*sqr(MBs)*sqr(MZs) + sqr(MZs)*(6*MWs*MZs*(-MTs + MZs) - 14*MZs*sqr(MWs) + 3*MTs*sqr(MZs) + 8*trip(MWs)) + MBs*(-32*powint(MWs,4) + 8*MWs*(3*MTs - 2*MZs)*sqr(MZs) + 38*sqr(MWs)*sqr(MZs) + 16*MZs*trip(MWs) - 3*(4*MTs + MZs)*trip(MZs))))/(768.*(-4*MBs + MZs)*sqr(MWs)*sqr(MWs - MZs));
break;
case UQU:
case CQU:
res = (deltaAlpha*MZs*Pi)/(32.*AL) + (3*(2*MWs - MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + (3*(-2*MWs + MZs)*real(B0fin(MZs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) - ((2*MWs - MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) - (18*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs + 16*sqr(MWs) - 11*sqr(MZs)) + MZs*(MTs*MZs*(2368*MWs*MZs - 1792*sqr(MWs) - 1269*sqr(MZs)) - 6*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*sqr(MZs)*(-98*MWs*MZs + 80*sqr(MWs) + 45*sqr(MZs))) + MBs*((3136*MWs*MZs - 2416*sqr(MWs) - 1485*sqr(MZs))*sqr(MZs) + 24*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*MTs*MZs*(-1184*MWs*MZs + 824*sqr(MWs) + 657*sqr(MZs))))/(6912.*MWs*(MWs - MZs)*(-4*MBs + MZs)*(-4*MTs + MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*(40*MWs*MZs - 32*sqr(MWs) - 17*sqr(MZs))*sqr(MZs) + MTs*MZs*(-80*MZs*sqr(MWs) + 88*MWs*sqr(MZs) + 64*trip(MWs) - 27*trip(MZs)) + 2*sqr(MTs)*(-592*MZs*sqr(MWs) + 275*MWs*sqr(MZs) + 320*trip(MWs) - 21*trip(MZs))))/(2304.*MWs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MBs,MBs))*(MBs*(16*MWs*MZs + 8*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) + MWs*(-4*MWs + MZs)*trip(MZs) + 2*sqr(MBs)*(-8*MZs*sqr(MWs) - 31*MWs*sqr(MZs) + 16*trip(MWs) + 17*trip(MZs))))/(768.*MWs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(-2*MWs + MZs)*sqr(MZs)*trip(MTs) + MTs*MZs*(-64*powint(MWs,4) - 18*MWs*(MBs - MZs)*sqr(MZs) - 130*sqr(MWs)*sqr(MZs) + 176*MZs*trip(MWs) + 9*MBs*trip(MZs)) + sqr(MTs)*(-640*powint(MWs,4) + 24*MWs*(3*MBs + 4*MZs)*sqr(MZs) - 622*sqr(MWs)*sqr(MZs) + 1184*MZs*trip(MWs) - 9*(4*MBs + MZs)*trip(MZs))))/(2304.*MTs*(-4*MTs + MZs)*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(12*(-2*MWs + MZs)*sqr(MBs)*sqr(MZs) + sqr(MZs)*(6*MWs*MZs*(-MTs + MZs) - 14*MZs*sqr(MWs) + 3*MTs*sqr(MZs) + 8*trip(MWs)) + MBs*(-32*powint(MWs,4) + 8*MWs*(3*MTs - 2*MZs)*sqr(MZs) + 38*sqr(MWs)*sqr(MZs) + 16*MZs*trip(MWs) - 3*(4*MTs + MZs)*trip(MZs))))/(768.*(-4*MBs + MZs)*sqr(MWs)*sqr(MWs - MZs));
break;
case DQU:
case SQU:
res = -(deltaAlpha*MZs*Pi)/(32.*AL) + (3*(-2*MWs + MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + (3*(2*MWs - MZs)*real(B0fin(MZs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + ((2*MWs - MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) - (18*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs + 16*sqr(MWs) - 11*sqr(MZs)) + MZs*(MTs*MZs*(2368*MWs*MZs - 1792*sqr(MWs) - 1269*sqr(MZs)) - 6*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*sqr(MZs)*(-98*MWs*MZs + 80*sqr(MWs) + 45*sqr(MZs))) + MBs*((3136*MWs*MZs - 2416*sqr(MWs) - 1485*sqr(MZs))*sqr(MZs) + 24*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*MTs*MZs*(-1184*MWs*MZs + 824*sqr(MWs) + 657*sqr(MZs))))/(6912.*MWs*(-4*MBs + MZs)*(-4*MTs + MZs)*(-MWs + MZs)) - (real(B0fin(MZs,MBs,MBs))*(MBs*(16*MWs*MZs + 8*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) + MWs*(-4*MWs + MZs)*trip(MZs) + 2*sqr(MBs)*(-8*MZs*sqr(MWs) - 31*MWs*sqr(MZs) + 16*trip(MWs) + 17*trip(MZs))))/(768.*MWs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs)) + MTs*MZs*(80*MZs*sqr(MWs) - 88*MWs*sqr(MZs) - 64*trip(MWs) + 27*trip(MZs)) + sqr(MTs)*(1184*MZs*sqr(MWs) - 550*MWs*sqr(MZs) - 640*trip(MWs) + 42*trip(MZs))))/(2304.*MWs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(-32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs - MZs)*sqr(MZs)*trip(MTs) + MTs*MZs*(64*powint(MWs,4) + 18*MWs*(MBs - MZs)*sqr(MZs) + 130*sqr(MWs)*sqr(MZs) - 176*MZs*trip(MWs) - 9*MBs*trip(MZs)) + sqr(MTs)*(640*powint(MWs,4) - 24*MWs*(3*MBs + 4*MZs)*sqr(MZs) + 622*sqr(MWs)*sqr(MZs) - 1184*MZs*trip(MWs) + 9*(4*MBs + MZs)*trip(MZs))))/(2304.*MTs*(-4*MTs + MZs)*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(12*(-2*MWs + MZs)*sqr(MBs)*sqr(MZs) + sqr(MZs)*(6*MWs*MZs*(-MTs + MZs) - 14*MZs*sqr(MWs) + 3*MTs*sqr(MZs) + 8*trip(MWs)) + MBs*(-32*powint(MWs,4) + 8*MWs*(3*MTs - 2*MZs)*sqr(MZs) + 38*sqr(MWs)*sqr(MZs) + 16*MZs*trip(MWs) - 3*(4*MTs + MZs)*trip(MZs))))/(768.*(4*MBs - MZs)*sqr(MWs)*sqr(MWs - MZs));
break;
case BQU:
res = -(deltaAlpha*MZs*Pi)/(32.*AL) + (3*(-2*MWs + MZs)*real(B0fin(MWs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + (3*(2*MWs - MZs)*real(B0fin(MZs,0,0))*sqr(MZs))/(128.*sqr(MWs - MZs)) + ((2*MWs - MZs)*real(B0fin(MWs,MTs,MBs))*(MTs*MWs + MBs*(-2*MTs + MWs) + sqr(MBs) + sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/(256.*sqr(MWs)*sqr(MWs - MZs)) - (18*(4*MTs - MZs)*sqr(MBs)*(8*MWs*MZs + 16*sqr(MWs) - 11*sqr(MZs)) + MZs*(MTs*MZs*(2368*MWs*MZs - 1792*sqr(MWs) - 1269*sqr(MZs)) - 6*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*sqr(MZs)*(-98*MWs*MZs + 80*sqr(MWs) + 45*sqr(MZs))) + MBs*((3136*MWs*MZs - 2416*sqr(MWs) - 1485*sqr(MZs))*sqr(MZs) + 24*sqr(MTs)*(-272*MWs*MZs + 320*sqr(MWs) + 39*sqr(MZs)) + 8*MTs*MZs*(-1184*MWs*MZs + 824*sqr(MWs) + 657*sqr(MZs))))/(6912.*MWs*(-4*MBs + MZs)*(-4*MTs + MZs)*(-MWs + MZs)) - (real(B0fin(MZs,MBs,MBs))*(MBs*(16*MWs*MZs + 8*sqr(MWs) - 9*sqr(MZs))*sqr(MZs) + MWs*(-4*MWs + MZs)*trip(MZs) + 2*sqr(MBs)*(-8*MZs*sqr(MWs) - 31*MWs*sqr(MZs) + 16*trip(MWs) + 17*trip(MZs))))/(768.*MWs*(-4*MBs + MZs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MTs,MTs))*(MWs*sqr(MZs)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs)) + MTs*MZs*(80*MZs*sqr(MWs) - 88*MWs*sqr(MZs) - 64*trip(MWs) + 27*trip(MZs)) + sqr(MTs)*(1184*MZs*sqr(MWs) - 550*MWs*sqr(MZs) - 640*trip(MWs) + 42*trip(MZs))))/(2304.*MWs*(-4*MTs + MZs)*sqr(MWs - MZs)) + (A0fin(MTs)*(-32*sqr(MWs)*sqr(MWs - MZs)*sqr(MZs) + 36*(2*MWs - MZs)*sqr(MZs)*trip(MTs) + MTs*MZs*(64*powint(MWs,4) + 18*MWs*(MBs - MZs)*sqr(MZs) + 130*sqr(MWs)*sqr(MZs) - 176*MZs*trip(MWs) - 9*MBs*trip(MZs)) + sqr(MTs)*(640*powint(MWs,4) - 24*MWs*(3*MBs + 4*MZs)*sqr(MZs) + 622*sqr(MWs)*sqr(MZs) - 1184*MZs*trip(MWs) + 9*(4*MBs + MZs)*trip(MZs))))/(2304.*MTs*(-4*MTs + MZs)*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MBs)*(12*(-2*MWs + MZs)*sqr(MBs)*sqr(MZs) + sqr(MZs)*(6*MWs*MZs*(-MTs + MZs) - 14*MZs*sqr(MWs) + 3*MTs*sqr(MZs) + 8*trip(MWs)) + MBs*(-32*powint(MWs,4) + 8*MWs*(3*MTs - 2*MZs)*sqr(MZs) + 38*sqr(MWs)*sqr(MZs) + 16*MZs*trip(MWs) - 3*(4*MTs + MZs)*trip(MZs))))/(768.*(4*MBs - MZs)*sqr(MWs)*sqr(MWs - MZs));
break;
}
break;
}
return(res*(8*pow(AL,1.5))/(sqrt(MWs)*sqrt(-MWs + MZs)*sqrt(Pi)));
}
// imag. part of fermionic one-loop Z vertex factor for k^2=s
double iz1fs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = ((76*MWs - 49*MZs)*imag(B0fin(s,0,0)))/144. + ((4*MWs - MZs)*(2*MBs + s)*imag(B0fin(s,MBs,MBs)))/(144.*s) + ((8*MWs - 5*MZs)*(2*MTs + s)*imag(B0fin(s,MTs,MTs)))/(72.*s);
break;
case NUE:
case NUM:
case NUT:
res = 0;
break;
case UQU:
case CQU:
res = ((-76*MWs + 49*MZs)*imag(B0fin(s,0,0)))/216. - ((4*MWs - MZs)*(2*MBs + s)*imag(B0fin(s,MBs,MBs)))/(216.*s) - ((8*MWs - 5*MZs)*(2*MTs + s)*imag(B0fin(s,MTs,MTs)))/(108.*s);
break;
case DQU:
case SQU:
res = ((76*MWs - 49*MZs)*imag(B0fin(s,0,0)))/432. + ((4*MWs - MZs)*(2*MBs + s)*imag(B0fin(s,MBs,MBs)))/(432.*s) + ((8*MWs - 5*MZs)*(2*MTs + s)*imag(B0fin(s,MTs,MTs)))/(216.*s);
break;
case BQU:
res = ((76*MWs - 49*MZs)*imag(B0fin(s,0,0)))/432. + ((4*MWs - MZs)*(2*MBs + s)*imag(B0fin(s,MBs,MBs)))/(432.*s) + ((8*MWs - 5*MZs)*(2*MTs + s)*imag(B0fin(s,MTs,MTs)))/(216.*s);
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = 0;
break;
case NUE:
case NUM:
case NUT:
res = 0;
break;
case UQU:
case CQU:
res = 0;
break;
case DQU:
case SQU:
res = 0;
break;
case BQU:
res = 0;
break;
}
break;
}
return(res*(8*pow(AL,1.5))/(sqrt(MWs)*sqrt(-MWs + MZs)*sqrt(Pi)));
}
// real part of bosonic one-loop Z vertex factor for k^2=s
double rz1bs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) + ((2*MWs - 3*MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) - ((2*MWs - 3*MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + (real(C0(s,0,0,0,0,MWs))*sqr(MZs)*sqr(MWs + s))/(64.*(MWs - MZs)*s) + (real(B0fin(s,MWs,MWs))*(MWs*MZs*(16*MZs + 31*s) + (26*MZs - 36*s)*sqr(MWs) + 2*s*sqr(MZs) - 48*trip(MWs)))/(192.*(MWs - MZs)*s) + (real(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-36*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 9*trip(MZs)))/(128.*MWs*(-MWs + MZs)*s) + (real(B0fin(MZs,MWs,MWs))*(-96*powint(MWs,4) + 59*powint(MZs,4) - 44*sqr(MWs)*sqr(MZs) + 248*MZs*trip(MWs) - 266*MWs*trip(MZs)))/(1536.*MZs*sqr(MWs - MZs)) + (A0fin(MHs)*(MZs*sqr(MHs)*(41*MWs*MZs - 44*sqr(MWs) - 12*sqr(MZs)) + 24*MHs*MWs*(4*MWs - 3*MZs)*sqr(MZs) + (-8*MWs*MZs + 8*sqr(MWs) + 3*sqr(MZs))*trip(MHs) + 24*MWs*(-4*MWs + 3*MZs)*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(-MWs + MZs)*sqr(MWs)) + (real(B0fin(s,0,0))*(-4*MZs*(19*MZs + 27*s)*sqr(MWs) + 12*MWs*(5*MZs + 7*s)*sqr(MZs) + 16*(2*MZs + 3*s)*trip(MWs) - 9*(2*MZs + 3*s)*trip(MZs)))/(256.*MWs*(-MWs + MZs)*s) + (A0fin(MWs)*(-8*MZs*(74*MZs - 5*s)*powint(MWs,4) + 96*(4*MZs + s)*powint(MWs,5) - 3*(MHs + MZs)*s*powint(MZs,4) + 4*(14*MZs - 99*s)*sqr(MZs)*trip(MWs) + 2*MWs*(MHs - 53*MZs)*s*trip(MZs) + 2*(76*MZs + 193*s)*sqr(MWs)*trip(MZs)))/(1536.*MZs*s*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MHs,MZs))*(MZs*sqr(MHs)*(107*MWs*MZs - 52*sqr(MWs) - 63*sqr(MZs)) + 4*MHs*sqr(MZs)*(-70*MWs*MZs + 32*sqr(MWs) + 45*sqr(MZs)) + (-16*MWs*MZs + 8*sqr(MWs) + 9*sqr(MZs))*trip(MHs) - 12*(-29*MWs*MZs + 12*sqr(MWs) + 21*sqr(MZs))*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + ((4*MWs - 3*MZs)*s*sqr(MHs)*sqr(MZs) + 2*MHs*(48*(4*MZs + s)*powint(MWs,4) - 6*s*powint(MZs,4) - 8*(8*MZs + 41*s)*sqr(MWs)*sqr(MZs) + 4*MZs*(-32*MZs + 13*s)*trip(MWs) + 293*MWs*s*trip(MZs)) - 4*MZs*(96*(4*MZs + s)*powint(MWs,4) - 15*s*powint(MZs,4) - 16*(8*MZs + 41*s)*sqr(MWs)*sqr(MZs) - 8*MZs*(32*MZs - 13*s)*trip(MWs) + 590*MWs*s*trip(MZs)))/(1536.*MWs*(MWs - MZs)*MZs*(-MHs + 4*MZs)*s) + (A0fin(MZs)*(MWs*MZs*s*sqr(MHs)*(16*MWs*MZs - 8*sqr(MWs) - 9*sqr(MZs)) - MHs*(-312*MZs*(2*MZs + 3*s)*powint(MWs,4) + 96*(2*MZs + 3*s)*powint(MWs,5) + MWs*(108*MZs + 139*s)*powint(MZs,4) + 3*s*powint(MZs,5) + 24*(33*MZs + 50*s)*sqr(MZs)*trip(MWs) - 13*(36*MZs + 55*s)*sqr(MWs)*trip(MZs)) + 4*MZs*(-312*MZs*(2*MZs + 3*s)*powint(MWs,4) + 96*(2*MZs + 3*s)*powint(MWs,5) + 2*MWs*(54*MZs + 83*s)*powint(MZs,4) + 3*s*powint(MZs,5) + 4*(198*MZs + 305*s)*sqr(MZs)*trip(MWs) - 2*(234*MZs + 379*s)*sqr(MWs)*trip(MZs))))/(1536.*MZs*(-MHs + 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case NUE:
case NUM:
case NUT:
res = (MWs*MZs*(2*MWs + s)*real(B0fin(s,MWs,MWs)))/(64.*(MWs - MZs)*s) + (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + ((2*MWs - MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) - ((2*MWs - MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + ((2*MWs - MZs)*MZs*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(64.*(-MWs + MZs)*s) + (MZs*real(B0fin(s,0,0))*(6*MWs*MZs*s + 4*(MZs - 3*s)*sqr(MWs) + (2*MZs + 3*s)*sqr(MZs) - 8*trip(MWs)))/(256.*MWs*(MWs - MZs)*s) + (real(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*trip(MZs))/(128.*MWs*(MWs - MZs)*s) + (real(B0fin(MZs,MWs,MWs))*(68*MZs*sqr(MWs) - 106*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*sqr(MWs - MZs)) + (MZs*A0fin(MHs)*((19*MWs - 4*MZs)*MZs*sqr(MHs) - 24*MHs*MWs*sqr(MZs) + (-4*MWs + MZs)*trip(MHs) + 24*MWs*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(-MWs + MZs)*sqr(MWs)) - (MZs*real(B0fin(MZs,MHs,MZs))*(MZs*(-29*MWs + 21*MZs)*sqr(MHs) + 4*MHs*(22*MWs - 15*MZs)*sqr(MZs) + (4*MWs - 3*MZs)*trip(MHs) + 12*(-11*MWs + 7*MZs)*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) - (A0fin(MWs)*(120*s*powint(MWs,4) - 2*MWs*(MHs - 5*MZs)*s*sqr(MZs) - 2*(12*MZs + 19*s)*sqr(MWs)*sqr(MZs) + 4*MZs*(6*MZs - 19*s)*trip(MWs) + (MHs + MZs)*s*trip(MZs)))/(1536.*s*sqr(MWs)*sqr(MWs - MZs)) + (sqr(MHs)*sqr(MZs) + 2*MHs*(84*MZs*sqr(MWs) - 84*MWs*sqr(MZs) + 60*trip(MWs) - trip(MZs)) + 4*MZs*(-168*MZs*sqr(MWs) + 168*MWs*sqr(MZs) - 120*trip(MWs) + trip(MZs)))/(1536.*MWs*(-MWs + MZs)*(-MHs + 4*MZs)) + (MZs*A0fin(MZs)*(MWs*(-4*MWs + 3*MZs)*s*sqr(MHs) - 4*MZs*(-2*MZs*(6*MZs + 25*s)*sqr(MWs) + 6*MWs*(2*MZs + 3*s)*sqr(MZs) + 48*s*trip(MWs) + s*trip(MZs)) + MHs*(-(MZs*(12*MZs + 37*s)*sqr(MWs)) + 3*MWs*(4*MZs + 3*s)*sqr(MZs) + 48*s*trip(MWs) + s*trip(MZs))))/(1536.*(MHs - 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case UQU:
case CQU:
res = (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) - ((2*MWs - 5*MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + ((2*MWs - 5*MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*(2*MWs + MZs)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(-MWs + MZs)*s) + (real(B0fin(s,MWs,MWs))*(-(MWs*MZs*(32*MZs + 59*s)) + (-46*MZs + 72*s)*sqr(MWs) - 4*s*sqr(MZs) + 96*trip(MWs)))/(576.*(MWs - MZs)*s) + (real(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-240*MZs*sqr(MWs) + 204*MWs*sqr(MZs) + 128*trip(MWs) - 65*trip(MZs)))/(3456.*MWs*(MWs - MZs)*s) + (real(B0fin(MZs,MWs,MWs))*(192*powint(MWs,4) - 101*powint(MZs,4) + 156*sqr(MWs)*sqr(MZs) - 376*MZs*trip(MWs) + 426*MWs*trip(MZs)))/(4608.*MZs*sqr(MWs - MZs)) + (A0fin(MHs)*(MZs*sqr(MHs)*(63*MWs*MZs - 88*sqr(MWs) - 20*sqr(MZs)) + 24*MHs*MWs*(8*MWs - 5*MZs)*sqr(MZs) + (-12*MWs*MZs + 16*sqr(MWs) + 5*sqr(MZs))*trip(MHs) + 24*MWs*(-8*MWs + 5*MZs)*trip(MZs)))/(4608.*MHs*(MHs - 4*MZs)*(MWs - MZs)*sqr(MWs)) + (real(B0fin(s,0,0))*(-12*MZs*(43*MZs + 69*s)*sqr(MWs) + 6*MWs*(68*MZs + 93*s)*sqr(MZs) + 8*(23*MZs + 48*s)*trip(MWs) - 65*(2*MZs + 3*s)*trip(MZs)))/(6912.*MWs*(MWs - MZs)*s) - (A0fin(MWs)*(-8*MZs*(148*MZs - 25*s)*powint(MWs,4) + 192*(4*MZs + s)*powint(MWs,5) - 5*(MHs + MZs)*s*powint(MZs,4) + 4*(34*MZs - 217*s)*sqr(MZs)*trip(MWs) + 2*MWs*(MHs - 101*MZs)*s*trip(MZs) + 2*(140*MZs + 367*s)*sqr(MWs)*trip(MZs)))/(4608.*MZs*s*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MHs,MZs))*(-4*MHs*sqr(MZs)*(-118*MWs*MZs + 64*sqr(MWs) + 75*sqr(MZs)) + MZs*sqr(MHs)*(-185*MWs*MZs + 104*sqr(MWs) + 105*sqr(MZs)) + (28*MWs*MZs - 16*sqr(MWs) - 15*sqr(MZs))*trip(MHs) + 12*(-47*MWs*MZs + 24*sqr(MWs) + 35*sqr(MZs))*trip(MZs)))/(4608.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + (3*(8*MWs - 5*MZs)*s*sqr(MHs)*sqr(MZs) + 2*MHs*(288*(4*MZs + s)*powint(MWs,4) + 5*s*powint(MZs,4) - 24*(16*MZs + 63*s)*sqr(MWs)*sqr(MZs) + (412*MZs*s - 768*sqr(MZs))*trip(MWs) + 1338*MWs*s*trip(MZs)) - 4*MZs*(576*(4*MZs + s)*powint(MWs,4) - 5*s*powint(MZs,4) - 48*(16*MZs + 63*s)*sqr(MWs)*sqr(MZs) + 8*MZs*(-192*MZs + 103*s)*trip(MWs) + 2700*MWs*s*trip(MZs)))/(13824.*MWs*MZs*(-MWs + MZs)*(-MHs + 4*MZs)*s) + (A0fin(MZs)*(3*MWs*MZs*s*sqr(MHs)*(-28*MWs*MZs + 16*sqr(MWs) + 15*sqr(MZs)) - 4*MZs*(-736*MZs*(2*MZs + 3*s)*powint(MWs,4) + 256*(2*MZs + 3*s)*powint(MWs,5) + 2*MWs*(130*MZs + 207*s)*powint(MZs,4) + 15*s*powint(MZs,5) + 24*(74*MZs + 113*s)*sqr(MZs)*trip(MWs) - 2*(538*MZs + 927*s)*sqr(MWs)*trip(MZs)) + MHs*(-736*MZs*(2*MZs + 3*s)*powint(MWs,4) + 256*(2*MZs + 3*s)*powint(MWs,5) + MWs*(260*MZs + 279*s)*powint(MZs,4) + 15*s*powint(MZs,5) + 48*(37*MZs + 54*s)*sqr(MZs)*trip(MWs) - (1076*MZs + 1635*s)*sqr(MWs)*trip(MZs))))/(13824.*MZs*(-MHs + 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case DQU:
case SQU:
res = (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) - ((2*MWs + MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + ((4*MWs - MZs)*MZs*(2*MWs + MZs)*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + ((4*MWs - MZs)*real(B0fin(s,0,0))*(4*MWs*MZs*(MZs + 15*s) + 4*(7*MZs - 3*s)*sqr(MWs) - 7*(2*MZs + 3*s)*sqr(MZs)))/(6912.*MWs*(MWs - MZs)*s) + ((4*MWs - MZs)*MZs*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(MWs - MZs)*s) + (real(B0fin(s,MWs,MWs))*(MWs*MZs*(16*MZs + 25*s) + 2*(7*MZs - 18*s)*sqr(MWs) + 2*s*sqr(MZs) - 48*trip(MWs)))/(576.*(MWs - MZs)*s) + (real(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-12*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 7*trip(MZs)))/(3456.*MWs*(-MWs + MZs)*s) - ((4*MWs - MZs)*real(B0fin(MZs,MWs,MWs))*(4*MZs*sqr(MWs) + 46*MWs*sqr(MZs) + 24*trip(MWs) + 25*trip(MZs)))/(4608.*MZs*sqr(MWs - MZs)) + (A0fin(MHs)*(MZs*sqr(MHs)*(3*MWs*MZs - 44*sqr(MWs) - 4*sqr(MZs)) + 24*MHs*MWs*(4*MWs - MZs)*sqr(MZs) + (8*sqr(MWs) + sqr(MZs))*trip(MHs) + 24*MWs*(-4*MWs + MZs)*trip(MZs)))/(4608.*MHs*(MHs - 4*MZs)*(-MWs + MZs)*sqr(MWs)) + (A0fin(MWs)*(-8*MZs*(74*MZs - 35*s)*powint(MWs,4) + 96*(4*MZs + s)*powint(MWs,5) - (MHs + MZs)*s*powint(MZs,4) + 4*(26*MZs - 137*s)*sqr(MZs)*trip(MWs) - 2*MWs*(MHs + 43*MZs)*s*trip(MZs) + 2*(52*MZs + 155*s)*sqr(MWs)*trip(MZs)))/(4608.*MZs*s*sqr(MWs)*sqr(MWs - MZs)) + (real(B0fin(MZs,MHs,MZs))*(MZs*sqr(MHs)*(49*MWs*MZs - 52*sqr(MWs) - 21*sqr(MZs)) + 4*MHs*sqr(MZs)*(-26*MWs*MZs + 32*sqr(MWs) + 15*sqr(MZs)) + (-8*MWs*MZs + 8*sqr(MWs) + 3*sqr(MZs))*trip(MHs) - 12*(-7*MWs*MZs + 12*sqr(MWs) + 7*sqr(MZs))*trip(MZs)))/(4608.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + (3*(4*MWs - MZs)*s*sqr(MHs)*sqr(MZs) + 2*MHs*(144*(4*MZs + s)*powint(MWs,4) + 4*s*powint(MZs,4) - 12*(16*MZs + 27*s)*sqr(MWs)*sqr(MZs) + 4*MZs*(-96*MZs + 113*s)*trip(MWs) + 255*MWs*s*trip(MZs)) - 4*MZs*(288*(4*MZs + s)*powint(MWs,4) + 5*s*powint(MZs,4) - 24*(16*MZs + 27*s)*sqr(MWs)*sqr(MZs) + 8*MZs*(-96*MZs + 113*s)*trip(MWs) + 522*MWs*s*trip(MZs)))/(13824.*MWs*(MWs - MZs)*MZs*(-MHs + 4*MZs)*s) - (A0fin(MZs)*(3*MWs*MZs*s*sqr(MHs)*(-8*MWs*MZs + 8*sqr(MWs) + 3*sqr(MZs)) + MHs*(4*MWs - MZs)*(8*(2*MZs + 3*s)*powint(MWs,4) - 3*s*powint(MZs,4) + 9*(4*MZs - s)*sqr(MWs)*sqr(MZs) - 12*MZs*(2*MZs + 3*s)*trip(MWs) - MWs*(28*MZs + 39*s)*trip(MZs)) - 4*MZs*(-56*MZs*(2*MZs + 3*s)*powint(MWs,4) + 32*(2*MZs + 3*s)*powint(MWs,5) + 2*MWs*(14*MZs + 27*s)*powint(MZs,4) + 3*s*powint(MZs,5) + 12*(14*MZs + 5*s)*sqr(MZs)*trip(MWs) - 2*(74*MZs + 99*s)*sqr(MWs)*trip(MZs))))/(13824.*MZs*(-MHs + 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case BQU:
res = (MZs*real(B0fin(s,MTs,MTs))*(-(MWs*(4*MWs - MZs)*(2*MWs + 3*s)) + 4*(MWs - MZs)*sqr(MTs) + 2*MTs*(-(MZs*s) + MWs*(MZs + s) + 2*sqr(MWs))))/(384.*MWs*(-MWs + MZs)*s) - ((2*MWs + MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + ((4*MWs - MZs)*MZs*(2*MWs + MZs)*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*(2*MWs + MZs)*A0fin(MTs)*(-6*MTs*MWs*(MWs + s) + s*sqr(MTs) + 4*(MWs + 2*s)*sqr(MWs) + 2*trip(MTs)))/(768.*MWs*(-MWs + MZs)*s*sqr(MTs - MWs)) + (MZs*real(C0(s,0,0,MTs,MTs,MWs))*((-6*MWs*MZs - 4*MWs*s + MZs*s)*sqr(MTs) + 4*MTs*MWs*(2*MWs*s + MZs*s + 3*sqr(MWs)) - 2*MWs*(4*MWs - MZs)*sqr(MWs + s) - 4*(MWs - MZs)*trip(MTs)))/(384.*MWs*(-MWs + MZs)*s) + (real(B0fin(s,MWs,MWs))*(-192*powint(MWs,4) + 2*MWs*MZs*(4*MZs*s + 9*MTs*(MZs + s) + 18*sqr(MTs)) + 4*MZs*(9*MTs + 16*MZs + 25*s)*sqr(MWs) - 9*MTs*(2*MTs + s)*sqr(MZs) + 8*(7*MZs - 18*s)*trip(MWs)))/(2304.*MWs*(MWs - MZs)*s) + (MZs*real(C0(0,s,0,MTs,MWs,MWs))*((MZs*s - 2*MWs*(MZs + s))*sqr(MTs) + MTs*MWs*(-4*MZs*s + MWs*(MZs + 4*s) + 6*sqr(MWs)) + (-2*MWs + MZs)*trip(MTs) - 4*(MWs + 2*s)*trip(MWs)))/(128.*MWs*(-MWs + MZs)*s) + ((2*MZs + 3*s)*real(B0fin(s,0,0))*(-12*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 7*trip(MZs)))/(6912.*MWs*(-MWs + MZs)*s) + (real(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-12*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 7*trip(MZs)))/(3456.*MWs*(-MWs + MZs)*s) - ((4*MWs - MZs)*real(B0fin(MZs,MWs,MWs))*(4*MZs*sqr(MWs) + 46*MWs*sqr(MZs) + 24*trip(MWs) + 25*trip(MZs)))/(4608.*MZs*sqr(MWs - MZs)) + (A0fin(MWs)*(136*powint(MWs,4) + (96*powint(MWs,5))/MZs - powint(MZs,4) + 252*MWs*MZs*sqr(MWs - MZs) + (6*MWs*MZs*(6*MTs*MWs + 4*sqr(MTs) - 4*sqr(MWs))*sqr(MWs - MZs))/sqr(MTs - MWs) + (192*(MZs - s)*sqr(MWs)*sqr(MWs - MZs))/s - 2*MHs*MWs*sqr(MZs) + (9*MWs*(MWs - MZs)*(6*MTs*MWs + 4*sqr(MTs) - 4*sqr(MWs))*sqr(MZs))/sqr(MTs - MWs) + 34*sqr(MWs)*sqr(MZs) - 248*MZs*trip(MWs) + (4*MWs*(MWs - MZs)*(3*MZs*(2*MWs + MZs)*sqr(MTs) + MTs*MWs*(-106*MWs*MZs + 96*sqr(MWs) + 19*sqr(MZs)) - 2*MWs*(MWs*MZs*(11*MZs + 12*s) - 50*MZs*sqr(MWs) + 6*s*sqr(MZs) + 48*trip(MWs))))/((MTs - MWs)*s) + 336*MWs*trip(MWs - MZs) - MHs*trip(MZs) + 34*MWs*trip(MZs)))/(4608.*sqr(MWs)*sqr(MWs - MZs)) + (A0fin(MHs)*(MZs*sqr(MHs)*(3*MWs*MZs - 44*sqr(MWs) - 4*sqr(MZs)) + 24*MHs*MWs*(4*MWs - MZs)*sqr(MZs) + (8*sqr(MWs) + sqr(MZs))*trip(MHs) + 24*MWs*(-4*MWs + MZs)*trip(MZs)))/(4608.*MHs*(MHs - 4*MZs)*(-MWs + MZs)*sqr(MWs)) + (real(B0fin(MZs,MHs,MZs))*(MZs*sqr(MHs)*(49*MWs*MZs - 52*sqr(MWs) - 21*sqr(MZs)) + 4*MHs*sqr(MZs)*(-26*MWs*MZs + 32*sqr(MWs) + 15*sqr(MZs)) + (-8*MWs*MZs + 8*sqr(MWs) + 3*sqr(MZs))*trip(MHs) - 12*(-7*MWs*MZs + 12*sqr(MWs) + 7*sqr(MZs))*trip(MZs)))/(4608.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) - (A0fin(MZs)*(3*MWs*MZs*s*sqr(MHs)*(-8*MWs*MZs + 8*sqr(MWs) + 3*sqr(MZs)) + MHs*(4*MWs - MZs)*(8*(2*MZs + 3*s)*powint(MWs,4) - 3*s*powint(MZs,4) + 9*(4*MZs - s)*sqr(MWs)*sqr(MZs) - 12*MZs*(2*MZs + 3*s)*trip(MWs) - MWs*(28*MZs + 39*s)*trip(MZs)) - 4*MZs*(-56*MZs*(2*MZs + 3*s)*powint(MWs,4) + 32*(2*MZs + 3*s)*powint(MWs,5) + 2*MWs*(14*MZs + 27*s)*powint(MZs,4) + 3*s*powint(MZs,5) + 12*(14*MZs + 5*s)*sqr(MZs)*trip(MWs) - 2*(74*MZs + 99*s)*sqr(MWs)*trip(MZs))))/(13824.*MZs*(-MHs + 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + (3*(MTs - MWs)*(4*MWs - MZs)*s*sqr(MHs)*sqr(MZs) + MHs*(-27*(6*MWs - 5*MZs)*s*sqr(MTs)*sqr(MZs) - 2*MWs*(144*(4*MZs + s)*powint(MWs,4) + 4*s*powint(MZs,4) - 12*(16*MZs + 27*s)*sqr(MWs)*sqr(MZs) + 4*MZs*(-96*MZs + 113*s)*trip(MWs) + 255*MWs*s*trip(MZs)) + MTs*(288*(4*MZs + s)*powint(MWs,4) + 8*s*powint(MZs,4) - 6*(64*MZs + 63*s)*sqr(MWs)*sqr(MZs) + 8*MZs*(-96*MZs + 113*s)*trip(MWs) + 429*MWs*s*trip(MZs))) + 4*MZs*(27*(6*MWs - 5*MZs)*s*sqr(MTs)*sqr(MZs) - MTs*(288*(4*MZs + s)*powint(MWs,4) + 5*s*powint(MZs,4) - 6*(64*MZs + 63*s)*sqr(MWs)*sqr(MZs) + 8*MZs*(-96*MZs + 113*s)*trip(MWs) + 441*MWs*s*trip(MZs)) + MWs*(288*(4*MZs + s)*powint(MWs,4) + 5*s*powint(MZs,4) - 24*(16*MZs + 27*s)*sqr(MWs)*sqr(MZs) + 8*MZs*(-96*MZs + 113*s)*trip(MWs) + 522*MWs*s*trip(MZs))))/(13824.*MWs*(-MTs + MWs)*MZs*(-MWs + MZs)*(-MHs + 4*MZs)*s);
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = (MWs*MZs*(2*MWs + s)*real(B0fin(s,MWs,MWs)))/(64.*(-MWs + MZs)*s) + (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) - ((2*MWs - MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + ((2*MWs - MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*real(B0fin(s,0,0))*(-12*MWs*MZs*(3*MZs + 5*s) + 4*(5*MZs + 9*s)*sqr(MWs) + 7*(2*MZs + 3*s)*sqr(MZs)))/(256.*MWs*(-MWs + MZs)*s) + (real(C0(s,0,0,0,0,MWs))*sqr(MZs)*sqr(MWs + s))/(64.*(MWs - MZs)*s) + (MZs*real(C0(s,0,0,0,0,MZs))*(-18*MWs*MZs + 12*sqr(MWs) + 7*sqr(MZs))*sqr(MZs + s))/(128.*MWs*(-MWs + MZs)*s) - (real(B0fin(MZs,MWs,MWs))*(68*MZs*sqr(MWs) - 106*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*sqr(MWs - MZs)) + (MZs*A0fin(MHs)*((19*MWs - 4*MZs)*MZs*sqr(MHs) - 24*MHs*MWs*sqr(MZs) + (-4*MWs + MZs)*trip(MHs) + 24*MWs*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(MWs - MZs)*sqr(MWs)) + (MZs*real(B0fin(MZs,MHs,MZs))*(MZs*(-29*MWs + 21*MZs)*sqr(MHs) + 4*MHs*(22*MWs - 15*MZs)*sqr(MZs) + (4*MWs - 3*MZs)*trip(MHs) + 12*(-11*MWs + 7*MZs)*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + (A0fin(MWs)*(24*(2*MZs + 5*s)*powint(MWs,4) + 2*MWs*(-MHs + 41*MZs)*s*sqr(MZs) + 2*(12*MZs - 91*s)*sqr(MWs)*sqr(MZs) - 4*MZs*(18*MZs + s)*trip(MWs) + (MHs + MZs)*s*trip(MZs)))/(1536.*s*sqr(MWs)*sqr(MWs - MZs)) + (sqr(MHs)*sqr(MZs) + 2*MHs*(108*MZs*sqr(MWs) - 117*MWs*sqr(MZs) + 60*trip(MWs) + 8*trip(MZs)) - 4*MZs*(216*MZs*sqr(MWs) - 234*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*MWs*(MWs - MZs)*(-MHs + 4*MZs)) + (A0fin(MZs)*(MWs*(4*MWs - 3*MZs)*MZs*s*sqr(MHs) + 4*MZs*(-72*(2*MZs + 3*s)*powint(MWs,4) + s*powint(MZs,4) - 2*(150*MZs + 241*s)*sqr(MWs)*sqr(MZs) + 12*MZs*(30*MZs + 49*s)*trip(MWs) + 42*MWs*(2*MZs + 3*s)*trip(MZs)) + MHs*(72*(2*MZs + 3*s)*powint(MWs,4) - s*powint(MZs,4) + (300*MZs + 469*s)*sqr(MWs)*sqr(MZs) - 12*MZs*(30*MZs + 49*s)*trip(MWs) - 3*MWs*(28*MZs + 39*s)*trip(MZs))))/(1536.*(MHs - 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case NUE:
case NUM:
case NUT:
res = (MWs*MZs*(2*MWs + s)*real(B0fin(s,MWs,MWs)))/(64.*(MWs - MZs)*s) + (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + ((2*MWs - MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) - ((2*MWs - MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + ((2*MWs - MZs)*MZs*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(64.*(-MWs + MZs)*s) + (MZs*real(B0fin(s,0,0))*(6*MWs*MZs*s + 4*(MZs - 3*s)*sqr(MWs) + (2*MZs + 3*s)*sqr(MZs) - 8*trip(MWs)))/(256.*MWs*(MWs - MZs)*s) + (real(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*trip(MZs))/(128.*MWs*(MWs - MZs)*s) + (real(B0fin(MZs,MWs,MWs))*(68*MZs*sqr(MWs) - 106*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*sqr(MWs - MZs)) + (MZs*A0fin(MHs)*((19*MWs - 4*MZs)*MZs*sqr(MHs) - 24*MHs*MWs*sqr(MZs) + (-4*MWs + MZs)*trip(MHs) + 24*MWs*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(-MWs + MZs)*sqr(MWs)) - (MZs*real(B0fin(MZs,MHs,MZs))*(MZs*(-29*MWs + 21*MZs)*sqr(MHs) + 4*MHs*(22*MWs - 15*MZs)*sqr(MZs) + (4*MWs - 3*MZs)*trip(MHs) + 12*(-11*MWs + 7*MZs)*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) - (A0fin(MWs)*(120*s*powint(MWs,4) - 2*MWs*(MHs - 5*MZs)*s*sqr(MZs) - 2*(12*MZs + 19*s)*sqr(MWs)*sqr(MZs) + 4*MZs*(6*MZs - 19*s)*trip(MWs) + (MHs + MZs)*s*trip(MZs)))/(1536.*s*sqr(MWs)*sqr(MWs - MZs)) + (sqr(MHs)*sqr(MZs) + 2*MHs*(84*MZs*sqr(MWs) - 84*MWs*sqr(MZs) + 60*trip(MWs) - trip(MZs)) + 4*MZs*(-168*MZs*sqr(MWs) + 168*MWs*sqr(MZs) - 120*trip(MWs) + trip(MZs)))/(1536.*MWs*(-MWs + MZs)*(-MHs + 4*MZs)) + (MZs*A0fin(MZs)*(MWs*(-4*MWs + 3*MZs)*s*sqr(MHs) - 4*MZs*(-2*MZs*(6*MZs + 25*s)*sqr(MWs) + 6*MWs*(2*MZs + 3*s)*sqr(MZs) + 48*s*trip(MWs) + s*trip(MZs)) + MHs*(-(MZs*(12*MZs + 37*s)*sqr(MWs)) + 3*MWs*(4*MZs + 3*s)*sqr(MZs) + 48*s*trip(MWs) + s*trip(MZs))))/(1536.*(MHs - 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case UQU:
case CQU:
res = (MWs*MZs*(2*MWs + s)*real(B0fin(s,MWs,MWs)))/(64.*(MWs - MZs)*s) + (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + ((2*MWs - MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) - ((2*MWs - MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*(2*MWs + MZs)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(-MWs + MZs)*s) + (MZs*real(C0(s,0,0,0,0,MZs))*(-20*MWs*MZs + 16*sqr(MWs) + 7*sqr(MZs))*sqr(MZs + s))/(384.*MWs*(MWs - MZs)*s) + (MZs*real(B0fin(s,0,0))*(-2*MWs*MZs*(20*MZs + 33*s) + 4*(7*MZs + 9*s)*sqr(MWs) + 7*(2*MZs + 3*s)*sqr(MZs) - 8*trip(MWs)))/(768.*MWs*(MWs - MZs)*s) + (real(B0fin(MZs,MWs,MWs))*(68*MZs*sqr(MWs) - 106*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*sqr(MWs - MZs)) + (MZs*A0fin(MHs)*((19*MWs - 4*MZs)*MZs*sqr(MHs) - 24*MHs*MWs*sqr(MZs) + (-4*MWs + MZs)*trip(MHs) + 24*MWs*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(-MWs + MZs)*sqr(MWs)) - (MZs*real(B0fin(MZs,MHs,MZs))*(MZs*(-29*MWs + 21*MZs)*sqr(MHs) + 4*MHs*(22*MWs - 15*MZs)*sqr(MZs) + (4*MWs - 3*MZs)*trip(MHs) + 12*(-11*MWs + 7*MZs)*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + (A0fin(MWs)*(-8*(4*MZs + 15*s)*powint(MWs,4) + 2*MWs*(MHs - 29*MZs)*s*sqr(MZs) + 2*(-4*MZs + 67*s)*sqr(MWs)*sqr(MZs) + 4*MZs*(10*MZs + 7*s)*trip(MWs) - (MHs + MZs)*s*trip(MZs)))/(1536.*s*sqr(MWs)*sqr(MWs - MZs)) + (sqr(MHs)*sqr(MZs) + 2*MHs*(96*MZs*sqr(MWs) - 98*MWs*sqr(MZs) + 60*trip(MWs) + trip(MZs)) - 4*MZs*(192*MZs*sqr(MWs) - 196*MWs*sqr(MZs) + 120*trip(MWs) + 3*trip(MZs)))/(1536.*MWs*(-MWs + MZs)*(-MHs + 4*MZs)) + (A0fin(MZs)*(MWs*MZs*(-4*MWs + 3*MZs)*s*sqr(MHs) - 4*MZs*(-32*(2*MZs + 3*s)*powint(MWs,4) + s*powint(MZs,4) - 2*(54*MZs + 97*s)*sqr(MWs)*sqr(MZs) + 24*MZs*(6*MZs + 11*s)*trip(MWs) + 14*MWs*(2*MZs + 3*s)*trip(MZs)) + MHs*(-32*(2*MZs + 3*s)*powint(MWs,4) + s*powint(MZs,4) - (108*MZs + 181*s)*sqr(MWs)*sqr(MZs) + 24*MZs*(6*MZs + 11*s)*trip(MWs) + MWs*(28*MZs + 33*s)*trip(MZs))))/(1536.*(MHs - 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case DQU:
case SQU:
res = (MWs*MZs*(2*MWs + s)*real(B0fin(s,MWs,MWs)))/(64.*(-MWs + MZs)*s) + (MZs*(MWs + 2*s)*real(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) - ((2*MWs - MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + ((2*MWs - MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + ((4*MWs - MZs)*MZs*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(MWs - MZs)*s) + (MZs*real(C0(s,0,0,0,0,MZs))*(-2*MWs*MZs + 4*sqr(MWs) + sqr(MZs))*sqr(MZs + s))/(384.*MWs*(-MWs + MZs)*s) + (MZs*real(B0fin(s,0,0))*(-12*(MZs - s)*sqr(MWs) + 4*MWs*sqr(MZs) - (2*MZs + 3*s)*sqr(MZs) + 16*trip(MWs)))/(768.*MWs*(MWs - MZs)*s) - (real(B0fin(MZs,MWs,MWs))*(68*MZs*sqr(MWs) - 106*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*sqr(MWs - MZs)) + (MZs*A0fin(MHs)*((19*MWs - 4*MZs)*MZs*sqr(MHs) - 24*MHs*MWs*sqr(MZs) + (-4*MWs + MZs)*trip(MHs) + 24*MWs*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(MWs - MZs)*sqr(MWs)) + (MZs*real(B0fin(MZs,MHs,MZs))*(MZs*(-29*MWs + 21*MZs)*sqr(MHs) + 4*MHs*(22*MWs - 15*MZs)*sqr(MZs) + (4*MWs - 3*MZs)*trip(MHs) + 12*(-11*MWs + 7*MZs)*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + (A0fin(MWs)*(8*(2*MZs + 15*s)*powint(MWs,4) + 2*MWs*(-MHs + 17*MZs)*s*sqr(MZs) - 2*(4*MZs + 43*s)*sqr(MWs)*sqr(MZs) - 4*MZs*(2*MZs + 13*s)*trip(MWs) + (MHs + MZs)*s*trip(MZs)))/(1536.*s*sqr(MWs)*sqr(MWs - MZs)) + (sqr(MHs)*sqr(MZs) + 2*MHs*(88*MZs*sqr(MWs) - 87*MWs*sqr(MZs) + 60*trip(MWs) - 2*trip(MZs)) + 4*MZs*(-176*MZs*sqr(MWs) + 174*MWs*sqr(MZs) - 120*trip(MWs) + 3*trip(MZs)))/(1536.*MWs*(MWs - MZs)*(-MHs + 4*MZs)) + (A0fin(MZs)*(MWs*(4*MWs - 3*MZs)*MZs*s*sqr(MHs) + MHs*(8*(2*MZs + 3*s)*powint(MWs,4) - s*powint(MZs,4) + (12*MZs + 37*s)*sqr(MWs)*sqr(MZs) - 12*MZs*(2*MZs + 7*s)*trip(MWs) + MWs*(-4*MZs + 3*s)*trip(MZs)) + 4*MZs*(-8*(2*MZs + 3*s)*powint(MWs,4) + s*powint(MZs,4) - 2*(6*MZs + 25*s)*sqr(MWs)*sqr(MZs) + 12*MZs*(2*MZs + 7*s)*trip(MWs) + 2*MWs*(2*MZs + 3*s)*trip(MZs))))/(1536.*(MHs - 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs));
break;
case BQU:
res = (MZs*real(B0fin(s,MWs,MWs))*((-4*MWs + 2*MZs)*sqr(MTs) + MTs*(MZs*s - 2*MWs*(MZs + s) - 4*sqr(MWs)) + 4*(2*MWs + s)*sqr(MWs)))/(256.*MWs*(-MWs + MZs)*s) + (MZs*real(B0fin(s,MTs,MTs))*(-(MWs*(4*MWs - MZs)*(2*MWs + 3*s)) + 4*(MWs - MZs)*sqr(MTs) + 2*MTs*(-(MZs*s) + MWs*(MZs + s) + 2*sqr(MWs))))/(384.*MWs*(-MWs + MZs)*s) - ((2*MWs - MZs)*real(B0fin(MWs,MHs,MWs))*(-4*MHs*MWs + sqr(MHs) + 12*sqr(MWs))*sqr(MZs))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*(2*MZs + 3*s)*real(B0fin(s,0,0))*(-2*MWs*MZs + 4*sqr(MWs) + sqr(MZs)))/(768.*MWs*(-MWs + MZs)*s) + ((2*MWs - MZs)*(4*MWs - MZs)*MZs*real(B0fin(MWs,MZs,MWs))*(20*MWs*MZs + 12*sqr(MWs) + sqr(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*real(C0(s,0,0,0,0,MZs))*(-2*MWs*MZs + 4*sqr(MWs) + sqr(MZs))*sqr(MZs + s))/(384.*MWs*(-MWs + MZs)*s) + (MZs*(2*MWs + MZs)*A0fin(MTs)*(-6*MTs*MWs*(MWs + s) + s*sqr(MTs) + 4*(MWs + 2*s)*sqr(MWs) + 2*trip(MTs)))/(768.*MWs*(-MWs + MZs)*s*sqr(MTs - MWs)) + (MZs*real(C0(s,0,0,MTs,MTs,MWs))*((-6*MWs*MZs - 4*MWs*s + MZs*s)*sqr(MTs) + 4*MTs*MWs*(2*MWs*s + MZs*s + 3*sqr(MWs)) - 2*MWs*(4*MWs - MZs)*sqr(MWs + s) - 4*(MWs - MZs)*trip(MTs)))/(384.*MWs*(-MWs + MZs)*s) + (MZs*real(C0(0,s,0,MTs,MWs,MWs))*((MZs*s - 2*MWs*(MZs + s))*sqr(MTs) + MTs*MWs*(-4*MZs*s + MWs*(MZs + 4*s) + 6*sqr(MWs)) + (-2*MWs + MZs)*trip(MTs) - 4*(MWs + 2*s)*trip(MWs)))/(128.*MWs*(-MWs + MZs)*s) - (real(B0fin(MZs,MWs,MWs))*(68*MZs*sqr(MWs) - 106*MWs*sqr(MZs) + 120*trip(MWs) + 17*trip(MZs)))/(1536.*sqr(MWs - MZs)) + (A0fin(MWs)*(120*powint(MWs,4) + powint(MZs,4) + (4*MWs*(MWs - MZs)*MZs*(2*MWs + MZs)*(MTs*MWs - 2*MWs*(MWs + 2*s) + sqr(MTs)))/((MTs - MWs)*s) + 84*MWs*MZs*sqr(MWs - MZs) + (4*MWs*MZs*(3*MTs*MWs + 2*sqr(MTs) - 2*sqr(MWs))*sqr(MWs - MZs))/sqr(MTs - MWs) - 2*MHs*MWs*sqr(MZs) + (6*MWs*(MWs - MZs)*(3*MTs*MWs + 2*sqr(MTs) - 2*sqr(MWs))*sqr(MZs))/sqr(MTs - MWs) + 94*sqr(MWs)*sqr(MZs) - 160*MZs*trip(MWs) + MHs*trip(MZs) - 38*MWs*trip(MZs)))/(1536.*sqr(MWs)*sqr(MWs - MZs)) + (MZs*A0fin(MHs)*((19*MWs - 4*MZs)*MZs*sqr(MHs) - 24*MHs*MWs*sqr(MZs) + (-4*MWs + MZs)*trip(MHs) + 24*MWs*trip(MZs)))/(1536.*MHs*(MHs - 4*MZs)*(MWs - MZs)*sqr(MWs)) + (MZs*real(B0fin(MZs,MHs,MZs))*(MZs*(-29*MWs + 21*MZs)*sqr(MHs) + 4*MHs*(22*MWs - 15*MZs)*sqr(MZs) + (4*MWs - 3*MZs)*trip(MHs) + 12*(-11*MWs + 7*MZs)*trip(MZs)))/(1536.*MWs*(MHs - 4*MZs)*sqr(MWs - MZs)) + (A0fin(MZs)*(MWs*(4*MWs - 3*MZs)*MZs*s*sqr(MHs) + MHs*(8*(2*MZs + 3*s)*powint(MWs,4) - s*powint(MZs,4) + (12*MZs + 37*s)*sqr(MWs)*sqr(MZs) - 12*MZs*(2*MZs + 7*s)*trip(MWs) + MWs*(-4*MZs + 3*s)*trip(MZs)) + 4*MZs*(-8*(2*MZs + 3*s)*powint(MWs,4) + s*powint(MZs,4) - 2*(6*MZs + 25*s)*sqr(MWs)*sqr(MZs) + 12*MZs*(2*MZs + 7*s)*trip(MWs) + 2*MWs*(2*MZs + 3*s)*trip(MZs))))/(1536.*(MHs - 4*MZs)*s*sqr(MWs)*sqr(MWs - MZs)) + ((MTs - MWs)*sqr(MHs)*sqr(MZs) + MHs*(-3*(6*MWs - 5*MZs)*MZs*sqr(MTs) + MTs*(206*MZs*sqr(MWs) - 183*MWs*sqr(MZs) + 120*trip(MWs) - 4*trip(MZs)) + 2*MWs*(-88*MZs*sqr(MWs) + 87*MWs*sqr(MZs) - 60*trip(MWs) + 2*trip(MZs))) + 4*MZs*(3*(6*MWs - 5*MZs)*MZs*sqr(MTs) + MWs*(176*MZs*sqr(MWs) - 174*MWs*sqr(MZs) + 120*trip(MWs) - 3*trip(MZs)) + MTs*(-206*MZs*sqr(MWs) + 183*MWs*sqr(MZs) - 120*trip(MWs) + 3*trip(MZs))))/(1536.*(MTs - MWs)*MWs*(MHs - 4*MZs)*(-MWs + MZs));
break;
}
break;
}
return(res*(8*pow(AL,1.5))/(sqrt(MWs)*sqrt(-MWs + MZs)*sqrt(Pi)));
}
// imag. part of bosonic one-loop Z vertex factor for k^2=s
double iz1bs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) + (imag(C0(s,0,0,0,0,MWs))*sqr(MZs)*sqr(MWs + s))/(64.*(MWs - MZs)*s) + (imag(B0fin(s,MWs,MWs))*(MWs*MZs*(16*MZs + 31*s) + (26*MZs - 36*s)*sqr(MWs) + 2*s*sqr(MZs) - 48*trip(MWs)))/(192.*(MWs - MZs)*s) + (imag(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-36*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 9*trip(MZs)))/(128.*MWs*(-MWs + MZs)*s) + (imag(B0fin(s,0,0))*(-4*MZs*(19*MZs + 27*s)*sqr(MWs) + 12*MWs*(5*MZs + 7*s)*sqr(MZs) + 16*(2*MZs + 3*s)*trip(MWs) - 9*(2*MZs + 3*s)*trip(MZs)))/(256.*MWs*(-MWs + MZs)*s);
break;
case NUE:
case NUM:
case NUT:
res = (MWs*MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs)))/(64.*(MWs - MZs)*s) + (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + ((2*MWs - MZs)*MZs*imag(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(64.*(-MWs + MZs)*s) + (MZs*imag(B0fin(s,0,0))*(6*MWs*MZs*s + 4*(MZs - 3*s)*sqr(MWs) + (2*MZs + 3*s)*sqr(MZs) - 8*trip(MWs)))/(256.*MWs*(MWs - MZs)*s) + (imag(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*trip(MZs))/(128.*MWs*(MWs - MZs)*s);
break;
case UQU:
case CQU:
res = (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + (MZs*(2*MWs + MZs)*imag(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(-MWs + MZs)*s) + (imag(B0fin(s,MWs,MWs))*(-(MWs*MZs*(32*MZs + 59*s)) + (-46*MZs + 72*s)*sqr(MWs) - 4*s*sqr(MZs) + 96*trip(MWs)))/(576.*(MWs - MZs)*s) + (imag(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-240*MZs*sqr(MWs) + 204*MWs*sqr(MZs) + 128*trip(MWs) - 65*trip(MZs)))/(3456.*MWs*(MWs - MZs)*s) + (imag(B0fin(s,0,0))*(-12*MZs*(43*MZs + 69*s)*sqr(MWs) + 6*MWs*(68*MZs + 93*s)*sqr(MZs) + 8*(23*MZs + 48*s)*trip(MWs) - 65*(2*MZs + 3*s)*trip(MZs)))/(6912.*MWs*(MWs - MZs)*s);
break;
case DQU:
case SQU:
res = (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) + ((4*MWs - MZs)*imag(B0fin(s,0,0))*(4*MWs*MZs*(MZs + 15*s) + 4*(7*MZs - 3*s)*sqr(MWs) - 7*(2*MZs + 3*s)*sqr(MZs)))/(6912.*MWs*(MWs - MZs)*s) + ((4*MWs - MZs)*MZs*imag(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(MWs - MZs)*s) + (imag(B0fin(s,MWs,MWs))*(MWs*MZs*(16*MZs + 25*s) + 2*(7*MZs - 18*s)*sqr(MWs) + 2*s*sqr(MZs) - 48*trip(MWs)))/(576.*(MWs - MZs)*s) + (imag(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-12*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 7*trip(MZs)))/(3456.*MWs*(-MWs + MZs)*s);
break;
case BQU:
res = (MZs*imag(B0fin(s,MTs,MTs))*(MWs*(4*MWs - MZs)*(2*MWs + 3*s) - 4*(MWs - MZs)*sqr(MTs) - 2*MTs*(-(MZs*s) + MWs*(MZs + s) + 2*sqr(MWs))))/(384.*MWs*(MWs - MZs)*s) + (MZs*imag(C0(s,0,0,MTs,MTs,MWs))*((-6*MWs*MZs - 4*MWs*s + MZs*s)*sqr(MTs) + 4*MTs*MWs*(2*MWs*s + MZs*s + 3*sqr(MWs)) - 2*MWs*(4*MWs - MZs)*sqr(MWs + s) - 4*(MWs - MZs)*trip(MTs)))/(384.*MWs*(-MWs + MZs)*s) + (imag(B0fin(s,MWs,MWs))*(-192*powint(MWs,4) + 2*MWs*MZs*(4*MZs*s + 9*MTs*(MZs + s) + 18*sqr(MTs)) + 4*MZs*(9*MTs + 16*MZs + 25*s)*sqr(MWs) - 9*MTs*(2*MTs + s)*sqr(MZs) + 8*(7*MZs - 18*s)*trip(MWs)))/(2304.*MWs*(MWs - MZs)*s) + (MZs*imag(C0(0,s,0,MTs,MWs,MWs))*((MZs*s - 2*MWs*(MZs + s))*sqr(MTs) + MTs*MWs*(-4*MZs*s + MWs*(MZs + 4*s) + 6*sqr(MWs)) + (-2*MWs + MZs)*trip(MTs) - 4*(MWs + 2*s)*trip(MWs)))/(128.*MWs*(-MWs + MZs)*s) + ((2*MZs + 3*s)*imag(B0fin(s,0,0))*(-12*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 7*trip(MZs)))/(6912.*MWs*(-MWs + MZs)*s) + (imag(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*(-12*MZs*sqr(MWs) + 30*MWs*sqr(MZs) + 16*trip(MWs) - 7*trip(MZs)))/(3456.*MWs*(-MWs + MZs)*s);
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = (MWs*MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs)))/(64.*(-MWs + MZs)*s) + (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) + (MZs*imag(B0fin(s,0,0))*(-12*MWs*MZs*(3*MZs + 5*s) + 4*(5*MZs + 9*s)*sqr(MWs) + 7*(2*MZs + 3*s)*sqr(MZs)))/(256.*MWs*(-MWs + MZs)*s) + (imag(C0(s,0,0,0,0,MWs))*sqr(MZs)*sqr(MWs + s))/(64.*(MWs - MZs)*s) + (MZs*imag(C0(s,0,0,0,0,MZs))*(-18*MWs*MZs + 12*sqr(MWs) + 7*sqr(MZs))*sqr(MZs + s))/(128.*MWs*(-MWs + MZs)*s);
break;
case NUE:
case NUM:
case NUT:
res = (MWs*MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs)))/(64.*(MWs - MZs)*s) + (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + ((2*MWs - MZs)*MZs*imag(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(64.*(-MWs + MZs)*s) + (MZs*imag(B0fin(s,0,0))*(6*MWs*MZs*s + 4*(MZs - 3*s)*sqr(MWs) + (2*MZs + 3*s)*sqr(MZs) - 8*trip(MWs)))/(256.*MWs*(MWs - MZs)*s) + (imag(C0(s,0,0,0,0,MZs))*sqr(MZs + s)*trip(MZs))/(128.*MWs*(MWs - MZs)*s);
break;
case UQU:
case CQU:
res = (MWs*MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs)))/(64.*(MWs - MZs)*s) + (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(-MWs + MZs)*s) + (MZs*(2*MWs + MZs)*imag(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(-MWs + MZs)*s) + (MZs*imag(C0(s,0,0,0,0,MZs))*(-20*MWs*MZs + 16*sqr(MWs) + 7*sqr(MZs))*sqr(MZs + s))/(384.*MWs*(MWs - MZs)*s) + (MZs*imag(B0fin(s,0,0))*(-2*MWs*MZs*(20*MZs + 33*s) + 4*(7*MZs + 9*s)*sqr(MWs) + 7*(2*MZs + 3*s)*sqr(MZs) - 8*trip(MWs)))/(768.*MWs*(MWs - MZs)*s);
break;
case DQU:
case SQU:
res = (MWs*MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs)))/(64.*(-MWs + MZs)*s) + (MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*sqr(MWs))/(32.*(MWs - MZs)*s) + ((4*MWs - MZs)*MZs*imag(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(192.*(MWs - MZs)*s) + (MZs*imag(C0(s,0,0,0,0,MZs))*(-2*MWs*MZs + 4*sqr(MWs) + sqr(MZs))*sqr(MZs + s))/(384.*MWs*(-MWs + MZs)*s) + (MZs*imag(B0fin(s,0,0))*(-12*(MZs - s)*sqr(MWs) + 4*MWs*sqr(MZs) - (2*MZs + 3*s)*sqr(MZs) + 16*trip(MWs)))/(768.*MWs*(MWs - MZs)*s);
break;
case BQU:
res = (MZs*imag(B0fin(s,MWs,MWs))*((-4*MWs + 2*MZs)*sqr(MTs) + MTs*(MZs*s - 2*MWs*(MZs + s) - 4*sqr(MWs)) + 4*(2*MWs + s)*sqr(MWs)))/(256.*MWs*(-MWs + MZs)*s) + (MZs*imag(B0fin(s,MTs,MTs))*(MWs*(4*MWs - MZs)*(2*MWs + 3*s) - 4*(MWs - MZs)*sqr(MTs) - 2*MTs*(-(MZs*s) + MWs*(MZs + s) + 2*sqr(MWs))))/(384.*MWs*(MWs - MZs)*s) + (MZs*(2*MZs + 3*s)*imag(B0fin(s,0,0))*(-2*MWs*MZs + 4*sqr(MWs) + sqr(MZs)))/(768.*MWs*(-MWs + MZs)*s) + (MZs*imag(C0(s,0,0,0,0,MZs))*(-2*MWs*MZs + 4*sqr(MWs) + sqr(MZs))*sqr(MZs + s))/(384.*MWs*(-MWs + MZs)*s) + (MZs*imag(C0(s,0,0,MTs,MTs,MWs))*((-6*MWs*MZs - 4*MWs*s + MZs*s)*sqr(MTs) + 4*MTs*MWs*(2*MWs*s + MZs*s + 3*sqr(MWs)) - 2*MWs*(4*MWs - MZs)*sqr(MWs + s) - 4*(MWs - MZs)*trip(MTs)))/(384.*MWs*(-MWs + MZs)*s) + (MZs*imag(C0(0,s,0,MTs,MWs,MWs))*((MZs*s - 2*MWs*(MZs + s))*sqr(MTs) + MTs*MWs*(-4*MZs*s + MWs*(MZs + 4*s) + 6*sqr(MWs)) + (-2*MWs + MZs)*trip(MTs) - 4*(MWs + 2*s)*trip(MWs)))/(128.*MWs*(-MWs + MZs)*s);
break;
}
break;
}
return(res*(8*pow(AL,1.5))/(sqrt(MWs)*sqrt(-MWs + MZs)*sqrt(Pi)));
}
// real part of fermionic one-loop photon vertex factor for k^2=s
double rg1fs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = 0;
break;
case NUE:
case NUM:
case NUT:
res = 0;
break;
case UQU:
case CQU:
res = 0;
break;
case DQU:
case SQU:
res = 0;
break;
case BQU:
res = 0;
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = 0;
break;
case NUE:
case NUM:
case NUT:
res = 0;
break;
case UQU:
case CQU:
res = 0;
break;
case DQU:
case SQU:
res = 0;
break;
case BQU:
res = 0;
break;
}
break;
}
return(res*1);
}
// imag. part of fermionic one-loop photon vertex factor for k^2=s
double ig1fs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = 0;
break;
case NUE:
case NUM:
case NUT:
res = 0;
break;
case UQU:
case CQU:
res = 0;
break;
case DQU:
case SQU:
res = 0;
break;
case BQU:
res = 0;
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = 0;
break;
case NUE:
case NUM:
case NUT:
res = 0;
break;
case UQU:
case CQU:
res = 0;
break;
case DQU:
case SQU:
res = 0;
break;
case BQU:
res = 0;
break;
}
break;
}
return(res*1);
}
// real part of bosonic one-loop photon vertex factor for k^2=s
double rg1bs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = -(MZs*(2*MWs + s)*A0fin(MWs)*pow(AL,1.5))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(-MWs + MZs)*s*sqrt(Pi)) + (pow(AL,1.5)*(-26*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(32.*MWs*(MWs - MZs)*sqrt(Pi)) - ((2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5)*(-12*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(16.*MWs*(MWs - MZs)*MZs*s*sqrt(Pi)) + ((2*MZs + 3*s)*pow(AL,1.5)*real(B0fin(s,0,0))*(-12*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*(-12*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))*sqr(MZs + s))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case NUE:
case NUM:
case NUT:
res = (MZs*pow(AL,1.5))/((2*MWs - 2*MZs)*sqrt(Pi)) + (MZs*A0fin(MWs)*pow(AL,1.5))/(4.*MWs*(-MWs + MZs)*sqrt(Pi)) + (MZs*(2*MWs + 3*s)*pow(AL,1.5)*real(B0fin(s,0,0)))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(4.*(MWs - MZs)*s*sqrt(Pi));
break;
case UQU:
case CQU:
res = (MZs*A0fin(MWs)*pow(AL,1.5))/((6*MWs*s - 6*MZs*s)*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(MWs - MZs)*s*sqrt(Pi)) - (pow(AL,1.5)*(-238*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs)))/(432.*MWs*(MWs - MZs)*sqrt(Pi)) + ((2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs)))/(216.*MWs*(MWs - MZs)*MZs*s*sqrt(Pi)) + (pow(AL,1.5)*real(B0fin(s,0,0))*(-46*MWs*MZs - 96*MWs*s + 51*MZs*s + 34*sqr(MZs)))/(216.*MWs*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(12.*(MWs - MZs)*s*sqrt(Pi)) - (pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))*sqr(MZs + s))/(108.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case DQU:
case SQU:
res = -(MZs*(2*MWs - 3*s)*A0fin(MWs)*pow(AL,1.5))/(24.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(-MWs + MZs)*s*sqrt(Pi)) + (pow(AL,1.5)*(-418*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(864.*MWs*(MWs - MZs)*sqrt(Pi)) - ((2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5)*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(432.*MWs*(MWs - MZs)*MZs*s*sqrt(Pi)) - (pow(AL,1.5)*real(B0fin(s,0,0))*(8*MWs*MZs*(MZs + 15*s) + 8*(7*MZs - 3*s)*sqr(MWs) - 5*(2*MZs + 3*s)*sqr(MZs)))/(432.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(6.*(MWs - MZs)*s*sqrt(Pi)) + (pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))*sqr(MZs + s))/(216.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case BQU:
res = -(MZs*pow(AL,1.5)*real(B0fin(s,MWs,MWs))*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + s) + 2*sqr(MTs)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*real(B0fin(s,MTs,MTs))*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + 3*s) + 2*sqr(MTs)))/(24.*MWs*(MWs - MZs)*s*sqrt(Pi)) - ((2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5)*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(432.*MWs*(MWs - MZs)*MZs*s*sqrt(Pi)) + ((2*MZs + 3*s)*pow(AL,1.5)*real(B0fin(s,0,0))*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(432.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (pow(AL,1.5)*(81*MZs*sqr(MTs) + MTs*(283*MWs*MZs - 8*sqr(MWs) - 5*sqr(MZs)) + MWs*(-418*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))))/(864.*MWs*(-MTs + MWs)*(MWs - MZs)*sqrt(Pi)) + (pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))*sqr(MZs + s))/(216.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(C0(0,s,0,MTs,MWs,MWs))*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*(MWs + 2*s)*sqr(MWs) + trip(MTs)))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(C0(s,0,0,MTs,MTs,MWs))*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*MWs*sqr(MWs + s) + trip(MTs)))/(12.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*A0fin(MWs)*pow(AL,1.5)*(MTs*MWs*(-6*MWs + 19*s) - 10*s*sqr(MTs) + 2*(2*MWs - 3*s)*sqr(MWs) + 2*trip(MTs)))/(48.*MWs*(MWs - MZs)*s*sqr(MTs - MWs)*sqrt(Pi)) + (MZs*A0fin(MTs)*pow(AL,1.5)*(-6*MTs*MWs*(MWs + s) + s*sqr(MTs) + 4*(MWs + 2*s)*sqr(MWs) + 2*trip(MTs)))/(48.*MWs*(MWs - MZs)*s*sqr(MTs - MWs)*sqrt(Pi));
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = -(MZs*(10*MWs + 3*MZs)*pow(AL,1.5))/(32.*MWs*(MWs - MZs)*sqrt(Pi)) - (MZs*(2*MWs + s)*A0fin(MWs)*pow(AL,1.5))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi)) - ((4*MWs - 3*MZs)*(2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + ((4*MWs - 3*MZs)*MZs*(2*MZs + 3*s)*pow(AL,1.5)*real(B0fin(s,0,0)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(-MWs + MZs)*s*sqrt(Pi)) + ((4*MWs - 3*MZs)*MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*sqr(MZs + s))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case NUE:
case NUM:
case NUT:
res = (MZs*pow(AL,1.5))/((2*MWs - 2*MZs)*sqrt(Pi)) + (MZs*A0fin(MWs)*pow(AL,1.5))/(4.*MWs*(-MWs + MZs)*sqrt(Pi)) + (MZs*(2*MWs + 3*s)*pow(AL,1.5)*real(B0fin(s,0,0)))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(4.*(MWs - MZs)*s*sqrt(Pi));
break;
case UQU:
case CQU:
res = (MZs*(58*MWs + 5*MZs)*pow(AL,1.5))/(144.*MWs*(MWs - MZs)*sqrt(Pi)) + (MZs*A0fin(MWs)*pow(AL,1.5))/((6*MWs*s - 6*MZs*s)*sqrt(Pi)) + ((8*MWs - 5*MZs)*(2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5))/(72.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(6*MWs - 5*(2*MZs + 3*s))*pow(AL,1.5)*real(B0fin(s,0,0)))/(72.*MWs*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(12.*(MWs - MZs)*s*sqrt(Pi)) - ((8*MWs - 5*MZs)*MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*sqr(MZs + s))/(36.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case DQU:
case SQU:
res = -(MZs*(134*MWs + MZs)*pow(AL,1.5))/(288.*MWs*(MWs - MZs)*sqrt(Pi)) - (MZs*(2*MWs - 3*s)*A0fin(MWs)*pow(AL,1.5))/(24.*MWs*(MWs - MZs)*s*sqrt(Pi)) - ((4*MWs - MZs)*(2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5))/(144.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*pow(AL,1.5)*real(B0fin(s,MWs,MWs)))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*pow(AL,1.5)*real(C0(0,s,0,0,MWs,MWs)))/(4.*(-MWs + MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(B0fin(s,0,0))*(-8*MWs*(MZs - 3*s) + MZs*(2*MZs + 3*s) + 24*sqr(MWs)))/(144.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MWs))*sqr(MWs + s))/(6.*(MWs - MZs)*s*sqrt(Pi)) + ((4*MWs - MZs)*MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*sqr(MZs + s))/(72.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case BQU:
res = -((4*MWs - MZs)*(2*MZs + 3*s)*A0fin(MZs)*pow(AL,1.5))/(144.*MWs*(MWs - MZs)*s*sqrt(Pi)) + ((4*MWs - MZs)*MZs*(2*MZs + 3*s)*pow(AL,1.5)*real(B0fin(s,0,0)))/(144.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(B0fin(s,MWs,MWs))*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + s) + 2*sqr(MTs)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*real(B0fin(s,MTs,MTs))*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + 3*s) + 2*sqr(MTs)))/(24.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*pow(AL,1.5)*(MTs*(89*MWs + MZs) - MWs*(134*MWs + MZs) + 27*sqr(MTs)))/(288.*MWs*(-MTs + MWs)*(MWs - MZs)*sqrt(Pi)) + ((4*MWs - MZs)*MZs*pow(AL,1.5)*real(C0(s,0,0,0,0,MZs))*sqr(MZs + s))/(72.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(C0(0,s,0,MTs,MWs,MWs))*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*(MWs + 2*s)*sqr(MWs) + trip(MTs)))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*pow(AL,1.5)*real(C0(s,0,0,MTs,MTs,MWs))*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*MWs*sqr(MWs + s) + trip(MTs)))/(12.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*A0fin(MWs)*pow(AL,1.5)*(MTs*MWs*(-6*MWs + 19*s) - 10*s*sqr(MTs) + 2*(2*MWs - 3*s)*sqr(MWs) + 2*trip(MTs)))/(48.*MWs*(MWs - MZs)*s*sqr(MTs - MWs)*sqrt(Pi)) + (MZs*A0fin(MTs)*pow(AL,1.5)*(-6*MTs*MWs*(MWs + s) + s*sqr(MTs) + 4*(MWs + 2*s)*sqr(MWs) + 2*trip(MTs)))/(48.*MWs*(MWs - MZs)*s*sqr(MTs - MWs)*sqrt(Pi));
break;
}
break;
}
return(res*1);
}
// imag. part of bosonic one-loop photon vertex factor for k^2=s
double ig1bs(int type, int form, double s, const inval& input)
{
double MZs = sqr(input.get(MZ)),
MWs = sqr(input.get(MW)),
MHs = sqr(input.get(MH)),
MTs = sqr(input.get(MT)),
MBs = sqr(input.get(MB)),
AL = input.get(al),
res;
switch(form) {
case VEC:
switch(type) {
case ELE:
case MUO:
case TAU:
res = (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(-MWs + MZs)*s*sqrt(Pi)) + ((2*MZs + 3*s)*imag(B0fin(s,0,0))*pow(AL,1.5)*(-12*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*(-12*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))*sqr(MZs + s))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case NUE:
case NUM:
case NUT:
res = (MZs*(2*MWs + 3*s)*imag(B0fin(s,0,0))*pow(AL,1.5))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*imag(C0(s,0,0,0,0,MWs))*pow(AL,1.5)*sqr(MWs + s))/(4.*(MWs - MZs)*s*sqrt(Pi));
break;
case UQU:
case CQU:
res = (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (imag(B0fin(s,0,0))*pow(AL,1.5)*(-46*MWs*MZs - 96*MWs*s + 51*MZs*s + 34*sqr(MZs)))/(216.*MWs*s*sqrt(Pi)) + (MZs*imag(C0(s,0,0,0,0,MWs))*pow(AL,1.5)*sqr(MWs + s))/(12.*(MWs - MZs)*s*sqrt(Pi)) - (imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*(-40*MWs*MZs + 32*sqr(MWs) + 17*sqr(MZs))*sqr(MZs + s))/(108.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case DQU:
case SQU:
res = (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(-MWs + MZs)*s*sqrt(Pi)) - (imag(B0fin(s,0,0))*pow(AL,1.5)*(8*MWs*MZs*(MZs + 15*s) + 8*(7*MZs - 3*s)*sqr(MWs) - 5*(2*MZs + 3*s)*sqr(MZs)))/(432.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(C0(s,0,0,0,0,MWs))*pow(AL,1.5)*sqr(MWs + s))/(6.*(MWs - MZs)*s*sqrt(Pi)) + (imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))*sqr(MZs + s))/(216.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case BQU:
res = -(MZs*imag(B0fin(s,MWs,MWs))*pow(AL,1.5)*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + s) + 2*sqr(MTs)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*imag(B0fin(s,MTs,MTs))*pow(AL,1.5)*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + 3*s) + 2*sqr(MTs)))/(24.*MWs*(MWs - MZs)*s*sqrt(Pi)) + ((2*MZs + 3*s)*imag(B0fin(s,0,0))*pow(AL,1.5)*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs)))/(432.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*(-4*MWs*MZs + 8*sqr(MWs) + 5*sqr(MZs))*sqr(MZs + s))/(216.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(C0(0,s,0,MTs,MWs,MWs))*pow(AL,1.5)*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*(MWs + 2*s)*sqr(MWs) + trip(MTs)))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(C0(s,0,0,MTs,MTs,MWs))*pow(AL,1.5)*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*MWs*sqr(MWs + s) + trip(MTs)))/(12.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
}
break;
case AXV:
switch(type) {
case ELE:
case MUO:
case TAU:
res = ((4*MWs - 3*MZs)*MZs*(2*MZs + 3*s)*imag(B0fin(s,0,0))*pow(AL,1.5))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(-MWs + MZs)*s*sqrt(Pi)) + ((4*MWs - 3*MZs)*MZs*imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*sqr(MZs + s))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case NUE:
case NUM:
case NUT:
res = (MZs*(2*MWs + 3*s)*imag(B0fin(s,0,0))*pow(AL,1.5))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*imag(C0(s,0,0,0,0,MWs))*pow(AL,1.5)*sqr(MWs + s))/(4.*(MWs - MZs)*s*sqrt(Pi));
break;
case UQU:
case CQU:
res = (MZs*(6*MWs - 5*(2*MZs + 3*s))*imag(B0fin(s,0,0))*pow(AL,1.5))/(72.*MWs*s*sqrt(Pi)) + (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(-MWs + MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(MWs - MZs)*s*sqrt(Pi)) + (MZs*imag(C0(s,0,0,0,0,MWs))*pow(AL,1.5)*sqr(MWs + s))/(12.*(MWs - MZs)*s*sqrt(Pi)) - ((8*MWs - 5*MZs)*MZs*imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*sqr(MZs + s))/(36.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case DQU:
case SQU:
res = (MZs*(2*MWs + s)*imag(B0fin(s,MWs,MWs))*pow(AL,1.5))/(8.*(MWs - MZs)*s*sqrt(Pi)) + (MWs*MZs*(MWs + 2*s)*imag(C0(0,s,0,0,MWs,MWs))*pow(AL,1.5))/(4.*(-MWs + MZs)*s*sqrt(Pi)) - (MZs*imag(B0fin(s,0,0))*pow(AL,1.5)*(-8*MWs*(MZs - 3*s) + MZs*(2*MZs + 3*s) + 24*sqr(MWs)))/(144.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(C0(s,0,0,0,0,MWs))*pow(AL,1.5)*sqr(MWs + s))/(6.*(MWs - MZs)*s*sqrt(Pi)) + ((4*MWs - MZs)*MZs*imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*sqr(MZs + s))/(72.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
case BQU:
res = ((4*MWs - MZs)*MZs*(2*MZs + 3*s)*imag(B0fin(s,0,0))*pow(AL,1.5))/(144.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(B0fin(s,MWs,MWs))*pow(AL,1.5)*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + s) + 2*sqr(MTs)))/(16.*MWs*(MWs - MZs)*s*sqrt(Pi)) + (MZs*imag(B0fin(s,MTs,MTs))*pow(AL,1.5)*(MTs*(2*MWs + s) - 2*MWs*(2*MWs + 3*s) + 2*sqr(MTs)))/(24.*MWs*(MWs - MZs)*s*sqrt(Pi)) + ((4*MWs - MZs)*MZs*imag(C0(s,0,0,0,0,MZs))*pow(AL,1.5)*sqr(MZs + s))/(72.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(C0(0,s,0,MTs,MWs,MWs))*pow(AL,1.5)*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*(MWs + 2*s)*sqr(MWs) + trip(MTs)))/(8.*MWs*(MWs - MZs)*s*sqrt(Pi)) - (MZs*imag(C0(s,0,0,MTs,MTs,MWs))*pow(AL,1.5)*(-(MTs*MWs*(3*MWs + 2*s)) + s*sqr(MTs) + 2*MWs*sqr(MWs + s) + trip(MTs)))/(12.*MWs*(MWs - MZs)*s*sqrt(Pi));
break;
}
break;
}
return(res*1);
}