1 #include "isajet/pilot.h"
4 C Calculate angular distributions for W decays from technirho:
5 C d(sigma)/d(qmw**2)d(yw)d(omega)d(omega1)d(omega2)
7 #if defined(CERNLIB_IMPNONE)
11 #include "isajet/itapes.inc"
12 #include "isajet/qcdpar.inc"
13 #include "isajet/jetpar.inc"
14 #include "isajet/pjets.inc"
15 #include "isajet/primar.inc"
16 #include "isajet/q1q2.inc"
17 #include "isajet/jetsig.inc"
18 #include "isajet/wsig.inc"
19 #include "isajet/wwsig.inc"
20 #include "isajet/wcon.inc"
21 #include "isajet/const.inc"
22 #include "isajet/wwpar.inc"
23 #include "isajet/tcpar.inc"
25 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
26 INTEGER I,K,IDADDR(4),IW(2)
27 REAL T12(3,3),T34(3,3),FR(3,3),FI(3,3),CPHI12(3),SPHI12(3),
28 $CPHI34(3),SPHI34(3),PFCM(5,4),PWCM(5,2),CHWW,SHWW,TMP,PTW1,
29 $CPHIW1,SPHIW1,PW1,CTHW1,STHW1,CHW1,SHW1,SHWI,TH12,PHI12,TH34,
30 $PHI34,AMV,GAMV,QMH,A12,B12,A34,B34,TVV12,TVA12,COS12,SIN12,
31 $TVV34,TVA34,COS34,SIN34,TCPHI,TSPHI,TC2PHI,TS2PHI,F0,F1,TOTAL,
36 C Reconstruct W-->FF decay angles
38 C Initialize PFCM and PWCM
48 C Z boost to WW center of mass
49 CHWW=QWJET(4)/QWJET(5)
50 SHWW=QWJET(3)/QWJET(5)
52 TMP=CHWW*PFCM(4,I)-SHWW*PFCM(3,I)
53 PFCM(3,I)=-SHWW*PFCM(4,I)+CHWW*PFCM(3,I)
57 TMP=CHWW*PWCM(4,I)-SHWW*PWCM(3,I)
58 PWCM(3,I)=-SHWW*PWCM(4,I)+CHWW*PWCM(3,I)
62 C Rotate W1 to +z axis
63 PTW1=SQRT(PWCM(1,1)**2+PWCM(2,1)**2)
66 PW1=SQRT(PTW1**2+PWCM(3,1)**2)
71 TMP=CPHIW1*PFCM(1,I)+SPHIW1*PFCM(2,I)
72 PFCM(2,I)=-SPHIW1*PFCM(1,I)+CPHIW1*PFCM(2,I)
77 TMP=CTHW1*PFCM(1,I)-STHW1*PFCM(3,I)
78 PFCM(3,I)=STHW1*PFCM(1,I)+CTHW1*PFCM(3,I)
82 C Boost to W rest frames
83 CHW1=PWCM(4,1)/PWCM(5,1)
91 TMP=CHW1*PFCM(4,I)-SHWI*PFCM(3,I)
92 PFCM(3,I)=-SHWI*PFCM(4,I)+CHW1*PFCM(3,I)
97 TH12=ACOS(PFCM(3,1)/SQRT(PFCM(1,1)**2+PFCM(2,1)**2+PFCM(3,1)**2))
98 PHI12=ATAN2(PFCM(2,1),PFCM(1,1))
99 TH34=ACOS(PFCM(3,3)/SQRT(PFCM(1,3)**2+PFCM(2,3)**2+PFCM(3,3)**2))
100 PHI34=ATAN2(PFCM(2,3),PFCM(1,3))
102 C Compute decay angular distributions.
105 IDADDR(I)=IABS(IDPAIR(I))
106 IF(IDADDR(I).GE.11) IDADDR(I)=IDADDR(I)-4
115 A12=AQ(IDADDR(1),IW(1))
116 B12=BQ(IDADDR(1),IW(1))
117 A34=AQ(IDADDR(3),IW(2))
118 B34=BQ(IDADDR(3),IW(2))
119 C DECAY DISTRIBUTIONS
120 TVV12=8.*PI*ALFA*(A12**2+B12**2)
121 TVA12=16.*PI*ALFA*A12*B12
124 T12(1,1)=TVV12*SIN12**2
125 T12(1,2)=TVV12*SIN12*COS12/SQRT2+TVA12*SIN12/SQRT2
126 T12(1,3)=-TVV12*SIN12*COS12/SQRT2+TVA12*SIN12/SQRT2
128 T12(2,2)=TVV12*(.5+.5*COS12**2)+TVA12*COS12
129 T12(2,3)=TVV12*.5*SIN12**2
132 T12(3,3)=TVV12*(.5+.5*COS12**2)-TVA12*COS12
134 TVV34=8.*PI*ALFA*(A34**2+B34**2)
135 TVA34=16.*PI*ALFA*A34*B34
138 T34(1,1)=TVV34*SIN34**2
139 T34(1,2)=TVV34*SIN34*COS34/SQRT2+TVA34*SIN34/SQRT2
140 T34(1,3)=-TVV34*SIN34*COS34/SQRT2+TVA34*SIN34/SQRT2
142 T34(2,2)=TVV34*(.5+.5*COS34**2)+TVA34*COS34
143 T34(2,3)=TVV34*.5*SIN34**2
146 T34(3,3)=TVV34*(.5+.5*COS34**2)-TVA34*COS34
150 CPHI12(3)=COS(2.*PHI12)
153 SPHI12(3)=SIN(2.*PHI12)
156 CPHI34(3)=COS(2.*PHI34)
159 SPHI34(3)=SIN(2.*PHI34)
161 TCPHI=CPHI12(2)*CPHI34(2)-SPHI12(2)*SPHI34(2)
162 TSPHI=SPHI12(2)*CPHI34(2)+CPHI12(2)*SPHI34(2)
163 TC2PHI=CPHI12(3)*CPHI34(3)-SPHI12(3)*SPHI34(3)
164 TS2PHI=SPHI12(3)*CPHI34(3)+CPHI12(3)*SPHI34(3)
166 C Pure technirho --> WW. Calculate angular distribution for
167 C decay and multiply by cross section.
169 F0=.5*QMH**2/AMV**2-1.
171 TOTAL=(8.*PI/3.)**2*TVV12*TVV34*(F0**2+2.*F1**2)
172 DIFF=F0**2*T12(1,1)*T34(1,1)
173 $+F0*F1*(2.*T12(1,2)*T34(1,2)+2.*T12(1,3)*T34(1,3))*TCPHI
174 $+F1**2*(T12(2,2)*T34(1,2)+T12(3,3)*T34(3,3)
175 $ +2.*T12(2,3)*T34(2,3)*TC2PHI)
176 WWSIG=SIGLLQ*DIFF/TOTAL