]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION SMSZG(Q,QB,K,E,EB,AQ,BQ,AE,BE) | |
3 | IMPLICIT NONE | |
4 | C | |
5 | C This does squared matrix element for q+qb -> Z+gamma | |
6 | C where Z-> e+eb | |
7 | C I have factored out 128*e^6*Q_q^2*|D_Z(z^2)| from | |
8 | C the expression. Also 1/12 from spin/color ave. is out. | |
9 | C | |
10 | REAL Q(5),QB(5),K(5),E(5),EB(5),AQ,BQ,AE,BE,SMSZG | |
11 | REAL M1S,M2S,M12 | |
12 | REAL EDQ,EBDK,EBDQ,EDK,QBDK,EDQB,EBDQB,QDK,QDQB | |
13 | EDQ=E(4)*Q(4)-E(1)*Q(1)-E(2)*Q(2)-E(3)*Q(3) | |
14 | EBDK=EB(4)*K(4)-EB(1)*K(1)-EB(2)*K(2)-EB(3)*K(3) | |
15 | EBDQ=EB(4)*Q(4)-EB(1)*Q(1)-EB(2)*Q(2)-EB(3)*Q(3) | |
16 | EDK=E(4)*K(4)-E(1)*K(1)-E(2)*K(2)-E(3)*K(3) | |
17 | QBDK=QB(4)*K(4)-QB(1)*K(1)-QB(2)*K(2)-QB(3)*K(3) | |
18 | EDQB=E(4)*QB(4)-E(1)*QB(1)-E(2)*QB(2)-E(3)*QB(3) | |
19 | EBDQB=EB(4)*QB(4)-EB(1)*QB(1)-EB(2)*QB(2)-EB(3)*QB(3) | |
20 | QDK=Q(4)*K(4)-Q(1)*K(1)-Q(2)*K(2)-Q(3)*K(3) | |
21 | QDQB=Q(4)*QB(4)-Q(1)*QB(1)-Q(2)*QB(2)-Q(3)*QB(3) | |
22 | M1S=(((AQ**2+BQ**2)*(AE**2+BE**2)-4*AQ*BQ*AE*BE)*EDQ*EBDK+ | |
23 | $((AQ**2+BQ**2)*(AE**2+BE**2)+4*AQ*BQ*AE*BE)*EBDQ*EDK)/ | |
24 | $4./QBDK | |
25 | M2S=(((AQ**2+BQ**2)*(AE**2+BE**2)+4*AQ*BQ*AE*BE)*EDQB*EBDK+ | |
26 | $((AQ**2+BQ**2)*(AE**2+BE**2)-4*AQ*BQ*AE*BE)*EBDQB*EDK)/ | |
27 | $4./QDK | |
28 | M12=(2*(AQ**2+BQ**2)*(AE**2+BE**2)*(EDQ*EBDQ*QBDK+EDQB*EBDQB*QDK) | |
29 | $+((AQ**2+BQ**2)*(AE**2+BE**2)-4*AQ*BQ*AE*BE)*(2*EDQ*EBDQB*QDQB+ | |
30 | $EDQ*EBDK*QDQB-EDK*EBDQB*QDQB+EDQ*EBDQB*QDK-EDQ*EBDQB*QBDK)+ | |
31 | $((AQ**2+BQ**2)*(AE**2+BE**2)+4*AQ*BQ*AE*BE)*(2*EDQB*EBDQ*QDQB- | |
32 | $EDQB*EBDK*QDQB+EDK*EBDQ*QDQB+EDQB*EBDQ*QDK-EDQB*EBDQ*QBDK))/ | |
33 | $4./QBDK/QDK | |
34 | SMSZG=M1S+M2S+M12 | |
35 | RETURN | |
36 | END |