C++ Interface to Tauola
main_2q_2tau_NLO.F
1 program main_upup_tautau_nlo
2 implicit none!
3!#include "s2n_declare.h"
4
5 integer*4 l1,l2,l3,l4,iz
6 integer*4 iqed,iew,iborn,gfscheme,ifgg
7 real*8 sum,sig(2),har,hai,nc,fc,tlmu2,cosf
8 real*8 s,t,u,betaf,sigma,mta,conhc,pi
9 complex*16 UpUp_TauTau_HA_NLO,ha
10 external UpUp_TauTau_HA_NLO
11
12**********************************************************
13* iqed = 0 without QED corrections
14* 1 with QED corrections
15* 2 ISR
16* 3 IFI
17* 4 FSR
18* 5 [IFI+FSR]
19*
20* iew = 0 without Weak corrections
21* 1 with Weak corrections
22*
23* iborn = 0 Born + 1-loop
24* 1 Born
25*
26* gfscheme = 0 alpha(0) calculation scheme
27* 1 GF
28* 2 GF prime
29*
30* ifgg = -1 ffgg = 0
31* 0 zoro
32* 1 [zoro+coeff*ff_gg]
33* 2 [zoro/(1d0 - coeff*ff_gg)]
34*
35********************************************************
36
37*///////////////////////
38*-- Flag settings
39 iqed = 0
40 iew = 1
41 iborn = 0
42 gfscheme = 0
43 ifgg = 1
44*///////////////////////
45 nc = 1d0
46 fc = 3d0
47*///////////////////////
48* 'tlmu2' is IR singularity regulator
49*(photon mass sqared). Affects NLO
50* amlitude in case QED is on(iqed=1)
51 tlmu2 = 1d-5
52 call flagset(iqed,iew,iborn,gfscheme,ifgg,nc,fc,tlmu2)
53*///////////////////////
54
55 call paraget(mta,conhc,pi)
56
57*///////////////////////
58 call printconsts(0)
59
60*///////////////////////
61* Mandelstam variables
62 s = 4900.d0
63*- banchmark values for sigma
64 cosf = 0.896067987d0 ! sigma = 0.314137724
65 cosf = 0.515891036d0 ! sigma = 0.563049574
66 cosf = 0.103059961d0 ! sigma = 1.37836685
67 cosf = 0.000969684d0 ! sigma = 1.66750761
68 cosf =-0.255370913d0 ! sigma = 2.54648731
69*-
70 betaf = sqrt(1d0-4d0*mta**2/s)
71 t = mta**2 - s/2*(1d0-betaf*cosf)
72 u = mta**2 - s/2*(1d0+betaf*cosf)
73*//////////////////////
74
75 do iz = 0,1
76 sum = 0d0
77 do l1 = 1,2
78 do l2 = 1,2
79 do l3 = 1,2
80 do l4 = 1,2
81 call upup(l1,l2,l3,l4,s,t,u,iz,har,hai);
82 sum = sum +har*har+hai*hai! ha*dconjg(ha)
83 enddo
84 enddo
85 enddo
86 enddo
87 sig(iz+1) = sum
88 enddo
89 sigma = conhc* ! to pbarn
90 & nc/fc*1d0/2/s *
91 & 1d0/4 * ! spin sum
92 & (sig(2) - sig(1)) * ! |Amp|^2 - linearized
93 & betaf/16/pi ! phase_space/dcos{theta}
94 print*,''
95 print*, "d_sigma/d_cos{theta} = ",sigma
96 print*,''
97
98 stop
99 end