[u/mrichter/AliRoot.git] / ISAJET / code / sigwh.F

#include "isajet/pilot.h"
      SUBROUTINE SIGWH
C
C          Calculate d(sigma)/d(pt**2)d(y1)d(y2) for WH and ZH
C          associated production.
C
C          SIGMA    = cross section summed over types allowed by
C                     JETTYPE cards.
C          SIGS(I)  = partial cross section for I1 + I2 --> I3 + I4
C          INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 +I1
C
C          Extra factor of 1/2 needed for nonidentical final jets.
C          Y=-log(tan(theta/2)) gives jacobean P1*P2/E1*E2
C
C
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
#include "isajet/itapes.inc"
#include "isajet/const.inc"
#include "isajet/jetpar.inc"
#include "isajet/jetsig.inc"
#include "isajet/primar.inc"
#include "isajet/q1q2.inc"
#include "isajet/qcdpar.inc"
#include "isajet/wcon.inc"
#include "isajet/hcon.inc"
#include "isajet/xmssm.inc"
C
      REAL X(2)
      EQUIVALENCE (X(1),X1)
      EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
      REAL SIG,S,T,U,FAC,AMW,AMZ,AMW2,AMZ2,E1,E2
      REAL QFCN,STRUC,SIGHW
      REAL PROPZ,PROPW,GV(2),GA(2),AMH,AMH2,GAMW,GAMZ
      INTEGER IS2UD(25),IQ,IH,I,IQ1,IQ2,IFLQ
      SAVE IS2UD
C
C          IS2UD: Susy jettype -> u/d code
      DATA IS2UD/0,1,1,2,2,2,2,1,1,2,2,1,1,1,1,2,2,2,2,1,1,2,2,1,1/

C          Functions
      QFCN(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
C
      IF (GOMSSM) THEN
        CALL SIGWHS
        RETURN
      END IF
C          Initialize
      DO 10 I=1,MXSIGS
10    SIGS(I)=0.
      SIGMA=0.
      NSIGS=0
C
      AMW=WMASS(2)
      AMW2=AMW**2
      AMZ=WMASS(4)
      AMZ2=AMZ**2
      AMH=HMASS
      AMH2=AMH**2
      GAMW=WGAM(2)
      GAMZ=WGAM(4)
      GV(1)=.25-2*SIN2W/3.
      GV(2)=-.25+SIN2W/3.
      GA(1)=-.25
      GA(2)=.25
C
C          WH production via W-*
C
      IF (GOQ(28,1).AND.GOQ(30,2)) THEN
          CALL TWOKIN(0.,0.,AMW,AMH)
          IF(X1.GE.1..OR.X2.GE.1.) GO TO 100
          E1=SQRT(P(1)**2+AMW**2)
          E2=SQRT(P(2)**2+AMH**2)
          FAC=1./(12.*PI*S**2)
          FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
          PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
          SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
     $          PROPW*TBRWW(3,1)
          SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
          CALL SIGFIL(SIG,3,4,28,30)
          SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
          CALL SIGFIL(SIG,4,3,28,30)
          SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
          CALL SIGFIL(SIG,9,6,28,30)
          SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
          CALL SIGFIL(SIG,6,9,28,30)
100       CONTINUE
      END IF
C
      IF (GOQ(30,1).AND.GOQ(28,2)) THEN
          CALL TWOKIN(0.,0.,AMH,AMW)
          IF(X1.GE.1..OR.X2.GE.1.) GO TO 110
          E1=SQRT(P(1)**2+AMH**2)
          E2=SQRT(P(2)**2+AMW**2)
          FAC=1./(12.*PI*S**2)
          FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
          PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
          SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
     $          PROPW*TBRWW(3,2)
          SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
          CALL SIGFIL(SIG,3,4,30,28)
          SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
          CALL SIGFIL(SIG,4,3,30,28)
          SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
          CALL SIGFIL(SIG,9,6,30,28)
          SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
          CALL SIGFIL(SIG,6,9,30,28)
110       CONTINUE
      END IF
C
C
C          WH production via W+*
C
      IF (GOQ(27,1).AND.GOQ(30,2)) THEN
          CALL TWOKIN(0.,0.,AMW,AMH)
          IF(X1.GE.1..OR.X2.GE.1.) GO TO 120
          E1=SQRT(P(1)**2+AMW**2)
          E2=SQRT(P(2)**2+AMH**2)
          FAC=1./(12.*PI*S**2)
          FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
          PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
          SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
     $          PROPW*TBRWW(2,1)
          SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
          CALL SIGFIL(SIG,2,5,27,30)
          SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
          CALL SIGFIL(SIG,5,2,27,30)
          SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
          CALL SIGFIL(SIG,8,7,27,30)
          SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
          CALL SIGFIL(SIG,7,8,27,30)
120       CONTINUE
      END IF
C
      IF (GOQ(30,1).AND.GOQ(27,2)) THEN
          CALL TWOKIN(0.,0.,AMH,AMW)
          IF(X1.GE.1..OR.X2.GE.1.) GO TO 130
          E1=SQRT(P(1)**2+AMH**2)
          E2=SQRT(P(2)**2+AMW**2)
          FAC=1./(12.*PI*S**2)
          FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
          PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
          SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
     $          PROPW*TBRWW(2,2)
          SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
          CALL SIGFIL(SIG,2,5,30,27)
          SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
          CALL SIGFIL(SIG,5,2,30,27)
          SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
          CALL SIGFIL(SIG,8,7,30,27)
          SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
          CALL SIGFIL(SIG,7,8,30,27)
130       CONTINUE
      END IF
C
C          ZH production via Z*
C          
      IF (GOQ(29,1).AND.GOQ(30,2)) THEN
          CALL TWOKIN(0.,0.,AMZ,AMH)
          IF(X1.GE.1..OR.X2.GE.1.) GO TO 200
          E1=SQRT(P(1)**2+AMZ2)
          E2=SQRT(P(2)**2+AMH2)
          FAC=1./(3.*PI*S**2)
          FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
          PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
          DO 210 IQ1=2,11
            IFLQ=IS2UD(IQ1)
            IQ2=MATCH(IQ1,4)
            IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 210
            SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
     $     (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,1)
            SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
            CALL SIGFIL(SIG,IQ1,IQ2,29,30)
210         CONTINUE
200       CONTINUE
      END IF
C          HZ production via Z*
C          
      IF (GOQ(30,1).AND.GOQ(29,2)) THEN
          CALL TWOKIN(0.,0.,AMH,AMZ)
          IF(X1.GE.1..OR.X2.GE.1.) GO TO 220
          E1=SQRT(P(1)**2+AMH2)
          E2=SQRT(P(2)**2+AMZ2)
          FAC=1./(3.*PI*S**2)
          FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
          PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
          DO 230 IQ1=2,11
            IFLQ=IS2UD(IQ1)
            IQ2=MATCH(IQ1,4)
            IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 230
            SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
     $     (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,2)
            SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
            CALL SIGFIL(SIG,IQ1,IQ2,30,29)
230         CONTINUE
220       CONTINUE
      END IF
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	SUBROUTINE SIGWH
	3	C
	4	C Calculate d(sigma)/d(pt**2)d(y1)d(y2) for WH and ZH
	5	C associated production.
	6	C
	7	C SIGMA = cross section summed over types allowed by
	8	C JETTYPE cards.
	9	C SIGS(I) = partial cross section for I1 + I2 --> I3 + I4
	10	C INOUT(I) = IOPAK*3I4 + IOPAK*2I3 + IOPAK*I2 +I1
	11	C
	12	C Extra factor of 1/2 needed for nonidentical final jets.
	13	C Y=-log(tan(theta/2)) gives jacobean P1P2/E1E2
	14	C
	15	C
	16	#if defined(CERNLIB_IMPNONE)
	17	IMPLICIT NONE
	18	#endif
	19	#include "isajet/itapes.inc"
	20	#include "isajet/const.inc"
	21	#include "isajet/jetpar.inc"
	22	#include "isajet/jetsig.inc"
	23	#include "isajet/primar.inc"
	24	#include "isajet/q1q2.inc"
	25	#include "isajet/qcdpar.inc"
	26	#include "isajet/wcon.inc"
	27	#include "isajet/hcon.inc"
	28	#include "isajet/xmssm.inc"
	29	C
	30	REAL X(2)
	31	EQUIVALENCE (X(1),X1)
	32	EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
	33	REAL SIG,S,T,U,FAC,AMW,AMZ,AMW2,AMZ2,E1,E2
	34	REAL QFCN,STRUC,SIGHW
	35	REAL PROPZ,PROPW,GV(2),GA(2),AMH,AMH2,GAMW,GAMZ
	36	INTEGER IS2UD(25),IQ,IH,I,IQ1,IQ2,IFLQ
	37	SAVE IS2UD
	38	C
	39	C IS2UD: Susy jettype -> u/d code
	40	DATA IS2UD/0,1,1,2,2,2,2,1,1,2,2,1,1,1,1,2,2,2,2,1,1,2,2,1,1/
	41
	42	C Functions
	43	QFCN(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
	44	C
	45	IF (GOMSSM) THEN
	46	CALL SIGWHS
	47	RETURN
	48	END IF
	49	C Initialize
	50	DO 10 I=1,MXSIGS
	51	10 SIGS(I)=0.
	52	SIGMA=0.
	53	NSIGS=0
	54	C
	55	AMW=WMASS(2)
	56	AMW2=AMW**2
	57	AMZ=WMASS(4)
	58	AMZ2=AMZ**2
	59	AMH=HMASS
	60	AMH2=AMH**2
	61	GAMW=WGAM(2)
	62	GAMZ=WGAM(4)
	63	GV(1)=.25-2*SIN2W/3.
	64	GV(2)=-.25+SIN2W/3.
65	GA(1)=-.25
66	GA(2)=.25
67	C
68	C WH production via W-*
69	C
70	IF (GOQ(28,1).AND.GOQ(30,2)) THEN
71	CALL TWOKIN(0.,0.,AMW,AMH)
72	IF(X1.GE.1..OR.X2.GE.1.) GO TO 100
73	E1=SQRT(P(1)2+AMW2)
74	E2=SQRT(P(2)2+AMH2)
75	FAC=1./(12.PIS**2)
76	FAC=FACS/SCM(P(1)P(2)/(E1E2))*UNITS
77	PROPW=(S-AMW2)2+AMW*2GAMW**2
78	SIGHW=GF*2AMW*8(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
79	$ PROPW*TBRWW(3,1)
80	SIG=.5SIGHWFACQFCN(3,1)QFCN(4,2)
81	CALL SIGFIL(SIG,3,4,28,30)
82	SIG=.5SIGHWFACQFCN(4,1)QFCN(3,2)
83	CALL SIGFIL(SIG,4,3,28,30)
84	SIG=.5SIGHWFACQFCN(9,1)QFCN(6,2)
85	CALL SIGFIL(SIG,9,6,28,30)
86	SIG=.5SIGHWFACQFCN(6,1)QFCN(9,2)
87	CALL SIGFIL(SIG,6,9,28,30)
88	100 CONTINUE
89	END IF
90	C
91	IF (GOQ(30,1).AND.GOQ(28,2)) THEN
92	CALL TWOKIN(0.,0.,AMH,AMW)
93	IF(X1.GE.1..OR.X2.GE.1.) GO TO 110
94	E1=SQRT(P(1)2+AMH2)
95	E2=SQRT(P(2)2+AMW2)
96	FAC=1./(12.PIS**2)
97	FAC=FACS/SCM(P(1)P(2)/(E1E2))*UNITS
98	PROPW=(S-AMW2)2+AMW*2GAMW**2
99	SIGHW=GF*2AMW*8(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
100	$ PROPW*TBRWW(3,2)
101	SIG=.5SIGHWFACQFCN(3,1)QFCN(4,2)
102	CALL SIGFIL(SIG,3,4,30,28)
103	SIG=.5SIGHWFACQFCN(4,1)QFCN(3,2)
104	CALL SIGFIL(SIG,4,3,30,28)
105	SIG=.5SIGHWFACQFCN(9,1)QFCN(6,2)
106	CALL SIGFIL(SIG,9,6,30,28)
107	SIG=.5SIGHWFACQFCN(6,1)QFCN(9,2)
108	CALL SIGFIL(SIG,6,9,30,28)
109	110 CONTINUE
110	END IF
111	C
112	C
113	C WH production via W+*
114	C
115	IF (GOQ(27,1).AND.GOQ(30,2)) THEN
116	CALL TWOKIN(0.,0.,AMW,AMH)
117	IF(X1.GE.1..OR.X2.GE.1.) GO TO 120
118	E1=SQRT(P(1)2+AMW2)
119	E2=SQRT(P(2)2+AMH2)
120	FAC=1./(12.PIS**2)
121	FAC=FACS/SCM(P(1)P(2)/(E1E2))*UNITS
122	PROPW=(S-AMW2)2+AMW*2GAMW**2
123	SIGHW=GF*2AMW*8(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
124	$ PROPW*TBRWW(2,1)
125	SIG=.5SIGHWFACQFCN(2,1)QFCN(5,2)
126	CALL SIGFIL(SIG,2,5,27,30)
127	SIG=.5SIGHWFACQFCN(5,1)QFCN(2,2)
128	CALL SIGFIL(SIG,5,2,27,30)
129	SIG=.5SIGHWFACQFCN(8,1)QFCN(7,2)
130	CALL SIGFIL(SIG,8,7,27,30)
131	SIG=.5SIGHWFACQFCN(7,1)QFCN(8,2)
132	CALL SIGFIL(SIG,7,8,27,30)
133	120 CONTINUE
134	END IF
135	C
136	IF (GOQ(30,1).AND.GOQ(27,2)) THEN
137	CALL TWOKIN(0.,0.,AMH,AMW)
138	IF(X1.GE.1..OR.X2.GE.1.) GO TO 130
139	E1=SQRT(P(1)2+AMH2)
140	E2=SQRT(P(2)2+AMW2)
141	FAC=1./(12.PIS**2)
142	FAC=FACS/SCM(P(1)P(2)/(E1E2))*UNITS
143	PROPW=(S-AMW2)2+AMW*2GAMW**2
144	SIGHW=GF*2AMW*8(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
145	$ PROPW*TBRWW(2,2)
146	SIG=.5SIGHWFACQFCN(2,1)QFCN(5,2)
147	CALL SIGFIL(SIG,2,5,30,27)
148	SIG=.5SIGHWFACQFCN(5,1)QFCN(2,2)
149	CALL SIGFIL(SIG,5,2,30,27)
150	SIG=.5SIGHWFACQFCN(8,1)QFCN(7,2)
151	CALL SIGFIL(SIG,8,7,30,27)
152	SIG=.5SIGHWFACQFCN(7,1)QFCN(8,2)
153	CALL SIGFIL(SIG,7,8,30,27)
154	130 CONTINUE
155	END IF
156	C
157	C ZH production via Z*
158	C
159	IF (GOQ(29,1).AND.GOQ(30,2)) THEN
160	CALL TWOKIN(0.,0.,AMZ,AMH)
161	IF(X1.GE.1..OR.X2.GE.1.) GO TO 200
162	E1=SQRT(P(1)**2+AMZ2)
163	E2=SQRT(P(2)**2+AMH2)
164	FAC=1./(3.PIS**2)
165	FAC=FACS/SCM(P(1)P(2)/(E1E2))*UNITS
166	PROPZ=(S-AMZ2)2+AMZ*2GAMZ**2
167	DO 210 IQ1=2,11
168	IFLQ=IS2UD(IQ1)
169	IQ2=MATCH(IQ1,4)
170	IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 210
171	SIG=GF*2AMZ*8(GV(IFLQ)2+GA(IFLQ)2)*
172	$ (S/AMZ2+(1.-T/AMZ2)(1.-U/AMZ2))/PROPZTBRWW(4,1)
173	SIG=.5SIGFACQFCN(IQ1,1)QFCN(IQ2,2)
174	CALL SIGFIL(SIG,IQ1,IQ2,29,30)
175	210 CONTINUE
176	200 CONTINUE
177	END IF
178	C HZ production via Z*
179	C
180	IF (GOQ(30,1).AND.GOQ(29,2)) THEN
181	CALL TWOKIN(0.,0.,AMH,AMZ)
182	IF(X1.GE.1..OR.X2.GE.1.) GO TO 220
183	E1=SQRT(P(1)**2+AMH2)
184	E2=SQRT(P(2)**2+AMZ2)
185	FAC=1./(3.PIS**2)
186	FAC=FACS/SCM(P(1)P(2)/(E1E2))*UNITS
187	PROPZ=(S-AMZ2)2+AMZ*2GAMZ**2
188	DO 230 IQ1=2,11
189	IFLQ=IS2UD(IQ1)
190	IQ2=MATCH(IQ1,4)
191	IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 230
192	SIG=GF*2AMZ*8(GV(IFLQ)2+GA(IFLQ)2)*
193	$ (S/AMZ2+(1.-T/AMZ2)(1.-U/AMZ2))/PROPZTBRWW(4,2)
194	SIG=.5SIGFACQFCN(IQ1,1)QFCN(IQ2,2)
195	CALL SIGFIL(SIG,IQ1,IQ2,30,29)
196	230 CONTINUE
197	220 CONTINUE
198	END IF
199	RETURN
200	END