]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HERWIG/src/hwhssg.f
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
ITS new Geometry files. Not yet ready for uses, committed to allow additional
[u/mrichter/AliRoot.git] / HERWIG / src / hwhssg.f
CommitLineData
3820ca8e 1
2CDECK ID>, HWHSSG.
3
4*CMZ :- -18/05/99 20.33.45 by Kosuke Odagiri
5
6*-- Author : Kosuke Odagiri
7
8C-----------------------------------------------------------------------
9
10 SUBROUTINE HWHSSG
11
12C-----------------------------------------------------------------------
13
14C SUSY 2 PARTON -> 2 GAUGINOS PROCESSES (1 - 3)
15
16C -> GAUGINO + SPARTON PROCESSES (4 - 7)
17
18C-----------------------------------------------------------------------
19
20 INCLUDE 'HERWIG61.INC'
21
22 DOUBLE PRECISION HWR, HWUALF, HWUAEM, EPS, HCS, RCS, DIST,
23
24 & ML(6), ML2(6), MR(6), MR2(6), MCH(2), MCH2(2), MNU(4), MNU2(4),
25
26 & MSQK, MG, MG2, SM, DM, DAB, QPE, SGN, PF, SQPE, EMSC2,
27
28 & FAC0, FACA, FACB, FACC, S, T, T3, U, U4, SN2TH
29
30 DOUBLE PRECISION M1(2,2,6), M2(4,4,6), M3(2,4,6,6),
31
32 & M4(4,6), M5(2,6,6), M6L(4,6), M6R(4,6), M7(2,2,6,6),
33
34 & XA(4), XB(4), XC(4), XD(4), MZ, MW, XW, SQXW, S2W, S22W
35
36 INTEGER I, IQ, IQ1, IQ2, IQ3, IQ4, IG1, IG2, IG3, IG4,
37
38 & ID1, ID2, IGL, SSL, SSR, GLU, SSNU, SSCH, INU, ICH, IWD(6), IPB
39
40 DOUBLE PRECISION DQD(6), DQU(6), HWHSS1
41
42 EXTERNAL HWR, HWUALF, HWUAEM, HWHSS1
43
44 SAVE HCS, M1, M2, M3, M4, M5, M6L, M6R, M7
45
46 PARAMETER (EPS = 1.D-9, IGL = 49, SSL = 400, SSR = 412, GLU = 449)
47
48 PARAMETER (SSNU = 449, SSCH = 453, INU = 49, ICH = 53)
49
50 DOUBLE COMPLEX Z, Z0, C1, C2, C3, GZ, GW, CLL, CLR, CRL, CRR
51
52 PARAMETER (Z = (0.D0,1.D0), Z0 = (0.D0,0.D0))
53
54 EQUIVALENCE (MZ, RMASS(200)), (MW, RMASS(198)), (MG, RMASS(GLU))
55
56 EQUIVALENCE (XA(1), ZMIXSS(1,1)), (XA(2), ZMIXSS(2,1))
57
58 EQUIVALENCE (XA(3), ZMIXSS(3,1)), (XA(4), ZMIXSS(4,1))
59
60 EQUIVALENCE (XB(1), ZMIXSS(1,2)), (XB(2), ZMIXSS(2,2))
61
62 EQUIVALENCE (XB(3), ZMIXSS(3,2)), (XB(4), ZMIXSS(4,2))
63
64 EQUIVALENCE (XC(1), ZMIXSS(1,3)), (XC(2), ZMIXSS(2,3))
65
66 EQUIVALENCE (XC(3), ZMIXSS(3,3)), (XC(4), ZMIXSS(4,3))
67
68 EQUIVALENCE (XD(1), ZMIXSS(1,4)), (XD(2), ZMIXSS(2,4))
69
70 EQUIVALENCE (XD(3), ZMIXSS(3,4)), (XD(4), ZMIXSS(4,4))
71
72 DATA IWD/2,1,4,3,6,5/
73
74 DATA DQD/ONE,ZERO,ONE,ZERO,ONE,ZERO/
75
76 DATA DQU/ZERO,ONE,ZERO,ONE,ZERO,ONE/
77
78C
79
80 CALL HWSGEN(.FALSE.)
81
82 IF (GENEV) THEN
83
84 RCS = HCS*HWR()
85
86 ELSE
87
88 SN2TH = 0.25D0 - 0.25D0*COSTH**2
89
90 S=XX(1)*XX(2)*PHEP(5,3)**2
91
92 EMSC2 = EMSCA**2
93
94 FAC0 = FACTSS*HWUAEM(EMSC2)
95
96c prefactor for pair production, includes 1/Nc colour factor
97
98 FACA = FAC0*HWUAEM(EMSC2) / CAFAC
99
100c prefactor for qq -> gaugino + gluino, includes CF/Nc colour factor
101
102 FACB = FAC0*HWUALF(1,EMSCA) * CFFAC / CAFAC
103
104c prefactor for qg -> gaugino + squark, includes 1/2Nc colour factor
105
106 FACC = FACB / CFFAC / TWO
107
108 MG2 = MG**2
109
110 GZ = S-MZ**2+Z*S/MZ*GAMZ
111
112 GW = S-MW**2+Z*S/MW*GAMW
113
114 DO IQ = 1,6
115
116 IQ1 = SSL + IQ
117
118 IQ2 = SSR + IQ
119
120 ML(IQ) = RMASS(IQ1)
121
122 ML2(IQ) = ML(IQ)**2
123
124 MR(IQ) = RMASS(IQ2)
125
126 MR2(IQ) = MR(IQ)**2
127
128 END DO
129
130 XW = TWO * SWEIN
131
132 SQXW = SQRT(XW)
133
134 S22W = XW * (TWO - XW)
135
136 S2W = SQRT(S22W)
137
138 DO IG1 = 1,4
139
140 MNU(IG1) = RMASS(IG1+SSNU)
141
142 MNU2(IG1) = MNU(IG1)**2
143
144 END DO
145
146 DO IG1 = 1,2
147
148 MCH(IG1) = RMASS(IG1+SSCH)
149
150 MCH2(IG1) = MCH(IG1)**2
151
152 END DO
153
154c _ ~+ ~-
155
156c (1) q q -> X X
157
158c a b
159
160 DO IG1 = 1,2
161
162 DO IG2 = 1,2
163
164 SM = MCH(IG1) + MCH(IG2)
165
166 QPE = S - SM**2
167
168 IF (QPE.GE.ZERO) THEN
169
170 DM = MCH(IG1) - MCH(IG2)
171
172 SQPE = SQRT(QPE*(S-DM**2))
173
174 PF = SQPE/S
175
176 T = (SQPE*COSTH - S + MCH2(IG1) + MCH2(IG2)) / TWO
177
178 U = - T - S + MCH2(IG1) + MCH2(IG2)
179
180 DAB = ABS(FLOAT(IG1+IG2-3))
181
182 C1 = (-WMXVSS(IG1,2)*WMXVSS(IG2,2)+DAB*S22W/XW)/S2W/GZ
183
184 C2 = (-WMXUSS(IG1,2)*WMXUSS(IG2,2)+DAB*S22W/XW)/S2W/GZ
185
186 SGN = WSGNSS(IG1)*WSGNSS(IG2)
187
188 DO IQ = 1,6
189
190 C3 = -DAB*QFCH(IQ)/S
191
192 CLL = C3 - LFCH(IQ)*C1 +
193
194 & DQD(IQ)*WMXUSS(IG1,1)*WMXUSS(IG2,1)/(U-ML2(IWD(IQ)))
195
196 CLR = C3 - LFCH(IQ)*C2 -
197
198 & DQU(IQ)*WMXVSS(IG1,1)*WMXVSS(IG2,1)/(T-ML2(IWD(IQ)))
199
200 CRL = C3 - RFCH(IQ)*C1
201
202 CRR = C3 - RFCH(IQ)*C2
203
204 M1(IG1,IG2,IQ)=FACA*PF*
205
206 & HWHSS1(S,T,U,MCH(IG1),MCH(IG2),SGN,CLL,CLR,CRL,CRR)
207
208 END DO
209
210 ELSE
211
212 DO IQ = 1,6
213
214 M1(IG1,IG2,IQ) = ZERO
215
216 END DO
217
218 END IF
219
220 END DO
221
222 END DO
223
224c _ ~o ~o
225
226c (2) q q -> X X
227
228c i j
229
230 DO IG1 = 1,4
231
232 DO IG2 = 1,4
233
234 SM = MNU(IG1) + MNU(IG2)
235
236 QPE = S - SM**2
237
238 IF (QPE.GE.ZERO) THEN
239
240 DM = MNU(IG1) - MNU(IG2)
241
242 SQPE = SQRT(QPE*(S-DM**2))
243
244 PF = SQPE/S
245
246 T = (SQPE*COSTH - S + MNU2(IG1) + MNU2(IG2)) / TWO
247
248 U = - T - S + MNU2(IG1) + MNU2(IG2)
249
250 C1 = (XD(IG1)*XD(IG2)-XC(IG1)*XC(IG2))/S2W/GZ
251
252 C2 = - C1
253
254 SGN = ZSGNSS(IG1)*ZSGNSS(IG2)
255
256 DO IQ = 1,6
257
258 CLL =LFCH(IQ)*C1+SLFCH(IQ,IG1)*SLFCH(IQ,IG2)/(U-ML2(IQ))
259
260 CLR =LFCH(IQ)*C2-SLFCH(IQ,IG1)*SLFCH(IQ,IG2)/(T-ML2(IQ))
261
262 CRL =RFCH(IQ)*C1-SRFCH(IQ,IG1)*SRFCH(IQ,IG2)/(T-MR2(IQ))
263
264 CRR =RFCH(IQ)*C2+SRFCH(IQ,IG1)*SRFCH(IQ,IG2)/(U-MR2(IQ))
265
266 M2(IG1,IG2,IQ) = FACA*PF*HALF*
267
268 & HWHSS1(S,T,U,MNU(IG1),MNU(IG2),SGN,CLL,CLR,CRL,CRR)
269
270 END DO
271
272 ELSE
273
274 DO IQ = 1,6
275
276 M2(IG1,IG2,IQ) = ZERO
277
278 END DO
279
280 END IF
281
282 END DO
283
284 END DO
285
286c _ ~+ ~o
287
288c (3) U D -> X X
289
290c a i
291
292 DO IG1 = 1,2
293
294 DO IG2 = 1,4
295
296 SM = MCH(IG1) + MNU(IG2)
297
298 QPE = S - SM**2
299
300 IF (QPE.GE.ZERO) THEN
301
302 DM = MCH(IG1) - MNU(IG2)
303
304 SQPE = SQRT(QPE*(S-DM**2))
305
306 PF = SQPE/S
307
308 T = (SQPE*COSTH - S + MCH2(IG1) + MNU2(IG2)) / TWO
309
310 U = - T - S + MCH2(IG1) + MNU2(IG2)
311
312 C1 = XA(IG2)+S2W/XW*XB(IG2)
313
314c note the new s-channel signs below.
315
316 C2 = ( XC(IG2)*WMXUSS(IG1,2)/SQXW+C1*WMXUSS(IG1,1))/GW
317
318 C3 = (-XD(IG2)*WMXVSS(IG1,2)/SQXW+C1*WMXVSS(IG1,1))/GW
319
320 SGN = WSGNSS(IG1)*ZSGNSS(IG2)
321
322 DO IQ1 = 1,3
323
324 IQ3 = IQ1*2
325
326 DO IQ2 = 1,3
327
328 IQ4 = IQ2*2-1
329
330 CLL = C2+WMXVSS(IG1,1)*SLFCH(IQ3,IG2)/(U-ML2(IQ3))
331
332 CLR = C3-WMXUSS(IG1,1)*SLFCH(IQ4,IG2)/(T-ML2(IQ4))
333
334 M3(IG1,IG2,IQ1,IQ2) = FACA*PF*VCKM(IQ1,IQ2)/XW*
335
336 & HWHSS1(S,T,U,MCH(IG1),MNU(IG2),SGN,CLL,CLR,Z0,Z0)
337
338 END DO
339
340 END DO
341
342 ELSE
343
344 DO IQ1 = 1,3
345
346 DO IQ2 = 1,3
347
348 M3(IG1,IG2,IQ1,IQ2) = ZERO
349
350 END DO
351
352 END DO
353
354 END IF
355
356 END DO
357
358 END DO
359
360c _ ~o ~
361
362c (4) q q -> X g
363
364c i
365
366 DO IG1 = 1,4
367
368 SM = MNU(IG1) + MG
369
370 QPE = S - SM**2
371
372 IF (QPE.GE.ZERO) THEN
373
374 DM = MNU(IG1) - MG
375
376 SQPE = SQRT(QPE*(S-DM**2))
377
378 PF = SQPE/S
379
380 T = (SQPE*COSTH - S + MG2 + MNU2(IG1)) / TWO
381
382 U = - T - S + MG2 + MNU2(IG1)
383
384 DO IQ = 1,6
385
386 CLL = SLFCH(IQ,IG1)/(U-ML2(IQ))
387
388 CLR = - SLFCH(IQ,IG1)/(T-ML2(IQ))
389
390 CRL = - SRFCH(IQ,IG1)/(T-MR2(IQ))
391
392 CRR = SRFCH(IQ,IG1)/(U-MR2(IQ))
393
394 M4(IG1,IQ) = FACB*PF*
395
396 & HWHSS1(S,T,U,MNU(IG1),MG,ZSGNSS(IG1),CLL,CLR,CRL,CRR)
397
398 END DO
399
400 ELSE
401
402 DO IQ = 1,6
403
404 M4(IG1,IQ) = ZERO
405
406 END DO
407
408 END IF
409
410 END DO
411
412c _ ~+ ~
413
414c (5) U D -> X g
415
416c a
417
418 DO IG1 = 1,2
419
420 SM = MCH(IG1) + MG
421
422 QPE = S - SM**2
423
424 IF (QPE.GE.ZERO) THEN
425
426 DM = MCH(IG1) - MG
427
428 SQPE = SQRT(QPE*(S-DM**2))
429
430 PF = SQPE/S
431
432 T = (SQPE*COSTH - S + MCH2(IG1) + MG2) / TWO
433
434 U = - T - S + MCH2(IG1) + MG2
435
436 DO IQ1 = 1,3
437
438 IQ3 = IQ1*2
439
440 DO IQ2 = 1,3
441
442 IQ4 = IQ2*2-1
443
444 CLL = WMXVSS(IG1,1)/(U-ML2(IQ3))
445
446 CLR = - WMXUSS(IG1,1)/(T-ML2(IQ4))
447
448 M5(IG1,IQ1,IQ2) = FACB*PF*VCKM(IQ1,IQ2)/XW*
449
450 & HWHSS1(S,T,U,MCH(IG1),MG,WSGNSS(IG1),CLL,CLR,Z0,Z0)
451
452 END DO
453
454 END DO
455
456 ELSE
457
458 DO IQ1 = 1,3
459
460 DO IQ2 = 1,3
461
462 M5(IG1,IQ1,IQ2) = ZERO
463
464 END DO
465
466 END DO
467
468 END IF
469
470 END DO
471
472c ~o ~
473
474c (6) g q -> X q
475
476c i LR
477
478 DO IG1 = 1,4
479
480 DO IQ = 1,6
481
482c left squarks
483
484 SM = MNU(IG1)+ML(IQ)
485
486 QPE = S - SM**2
487
488 IF (QPE.GE.ZERO) THEN
489
490 DM = MNU(IG1)-ML(IQ)
491
492 SQPE = SQRT(QPE*(S-DM**2))
493
494 PF = SQPE/S
495
496 T3 = (SQPE*COSTH - S - SM*DM) / TWO
497
498 U4 = - T3 - S
499
500 M6L(IG1,IQ) = FACC*PF * (
501
502 & QMIXSS(IQ,1,1)*SLFCH(IQ,IG1)+QMIXSS(IQ,2,1)*SRFCH(IQ,IG1)
503
504 & )**2 * T3/S/U4*(-U4+TWO*SM*DM/T3/U4*SQPE*SQPE*SN2TH)
505
506 ELSE
507
508 M6L(IG1,IQ) = ZERO
509
510 END IF
511
512c right squarks
513
514 SM = MNU(IG1)+MR(IQ)
515
516 QPE = S - SM**2
517
518 IF (QPE.GE.ZERO) THEN
519
520 DM = MNU(IG1)-MR(IQ)
521
522 SQPE = SQRT(QPE*(S-DM**2))
523
524 PF = SQPE/S
525
526 T3 = (SQPE*COSTH - S - SM*DM) / TWO
527
528 U4 = - T3 - S
529
530 M6R(IG1,IQ) = FACC*PF * (
531
532 & QMIXSS(IQ,1,2)*SLFCH(IQ,IG1)+QMIXSS(IQ,2,2)*SRFCH(IQ,IG1)
533
534 & )**2 * T3/S/U4*(-U4+TWO*SM*DM/T3/U4*SQPE*SQPE*SN2TH)
535
536 ELSE
537
538 M6R(IG1,IQ) = ZERO
539
540 END IF
541
542 END DO
543
544 END DO
545
546c ~+-~
547
548c (7) g q -> X q'
549
550c a L
551
552 DO IG1 = 1,2
553
554 DO IQ1 = 1,3
555
556 IQ3 = IQ1*2
557
558 DO IQ2 = 1,3
559
560 IQ4 = IQ2*2-1
561
562 DO I = 1,2
563
564c U initiated processes
565
566 IF (I.EQ.1) THEN
567
568 MSQK = ML(IQ4)
569
570 ELSE
571
572 MSQK = MR(IQ4)
573
574 END IF
575
576 SM = MCH(IG1) + MSQK
577
578 QPE = S - SM**2
579
580 IF (((I.EQ.1).OR.(IQ2.EQ.3)).AND.(QPE.GE.ZERO)) THEN
581
582 DM = MCH(IG1) - MSQK
583
584 SQPE = SQRT(QPE*(S-DM**2))
585
586 PF = SQPE/S
587
588 T3 = (SQPE*COSTH - S - SM*DM) / TWO
589
590 U4 = - T3 - S
591
592 M7(I,IG1,IQ3,IQ4)=FACC*PF*WMXUSS(IG1,1)**2*VCKM(IQ1,IQ2)
593
594 & /XW*T3/S/U4*(-U4+TWO*SM*DM/T3/U4*SQPE*SQPE*SN2TH)*
595
596 & QMIXSS(IQ4,1,I)**2
597
598 ELSE
599
600 M7(I,IG1,IQ3,IQ4) = ZERO
601
602 END IF
603
604c D initiated processes
605
606 IF (I.EQ.1) THEN
607
608 MSQK = ML(IQ3)
609
610 ELSE
611
612 MSQK = MR(IQ3)
613
614 END IF
615
616 SM = MCH(IG1) + MSQK
617
618 QPE = S - SM**2
619
620 IF (((I.EQ.1).OR.(IQ1.EQ.3)).AND.(QPE.GE.ZERO)) THEN
621
622 DM = MCH(IG1) - MSQK
623
624 SQPE = SQRT(QPE*(S-DM**2))
625
626 PF = SQPE/S
627
628 T3 = (SQPE*COSTH - S - SM*DM) / TWO
629
630 U4 = - T3 - S
631
632 M7(I,IG1,IQ4,IQ3)=FACC*PF*WMXVSS(IG1,1)**2*VCKM(IQ1,IQ2)
633
634 & /XW*T3/S/U4*(-U4+TWO*SM*DM/T3/U4*SQPE*SQPE*SN2TH)*
635
636 & QMIXSS(IQ3,1,I)**2
637
638 ELSE
639
640 M7(I,IG1,IQ4,IQ3) = ZERO
641
642 END IF
643
644 END DO
645
646 END DO
647
648 END DO
649
650 END DO
651
652 END IF
653
654 HCS = 0.
655
656c _ _ ~+ ~- ~o ~o ~o ~
657
658c q q , q q -> X X , X X , X g
659
660c a b i j i
661
662 DO 1 ID1 = 1,12
663
664 IF (DISF(ID1,1).LT.EPS) GOTO 1
665
666 IF (ID1.GT.6) THEN
667
668 ID2 = ID1 - 6
669
670 IQ = ID2
671
672 IPB = 4132
673
674 ELSE
675
676 ID2 = ID1 + 6
677
678 IQ = ID1
679
680 IPB = 2431
681
682 END IF
683
684 IF (DISF(ID2,2).LT.EPS) GOTO 1
685
686 DIST = DISF(ID1,1)*DISF(ID2,2)
687
688 DO IG1 = 1,2
689
690 IG3 = ICH+IG1
691
692 DO IG2 = 1,2
693
694 IG4 = ICH+IG2+2
695
696 HCS = HCS + DIST*M1(IG1,IG2,IQ)
697
698 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IG4,0,2134,21,*9)
699
700 END DO
701
702 END DO
703
704 DO IG1 = 1,4
705
706 IG3 = INU+IG1
707
708 DO IG2 = 1,4
709
710 IG4 = INU+IG2
711
712 IF (IG2.GE.IG1) HCS = HCS + DIST*M2(IG1,IG2,IQ)
713
714 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IG4,0,2134,22,*9)
715
716 END DO
717
718 HCS = HCS + DIST*M4(IG1,IQ)
719
720 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IGL,0, IPB,24,*9)
721
722 END DO
723
724 1 CONTINUE
725
726c _ _ ~+-~o ~+-~
727
728c q q', q q' -> X X , X g
729
730c a i a
731
732c
733
734c _ _ _ _
735
736c ud(+), ud(-), du(-), du(+)
737
738 DO 2 IQ1 = 1, 3
739
740 DO IQ2 = 1, 3
741
742 IF(VCKM(IQ1,IQ2).GT.EPS) THEN
743
744c _
745
746c ud (+)
747
748 ID1 = IQ1 * 2
749
750 ID2 = IQ2 * 2 + 5
751
752 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
753
754 DIST = DISF(ID1,1)*DISF(ID2,2)
755
756 DO IG1 = 1,2
757
758 IG3 = ICH+IG1
759
760 DO IG2 = 1,4
761
762 IG4 = INU+IG2
763
764 HCS = HCS + DIST*M3(IG1,IG2,IQ1,IQ2)
765
766 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IG4,0,2134,23,*9)
767
768 END DO
769
770 HCS = HCS + DIST*M5(IG1,IQ1,IQ2)
771
772 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IGL,0,2431,25,*9)
773
774 END DO
775
776 END IF
777
778c _
779
780c du (+)
781
782 ID1 = IQ2 * 2 + 5
783
784 ID2 = IQ1 * 2
785
786 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
787
788 DIST = DISF(ID1,1)*DISF(ID2,2)
789
790 DO IG1 = 1,2
791
792 IG3 = ICH+IG1
793
794 DO IG2 = 1,4
795
796 IG4 = INU+IG2
797
798 HCS = HCS + DIST*M3(IG1,IG2,IQ1,IQ2)
799
800 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG4,0,IG3,0,2134,23,*9)
801
802 END DO
803
804 HCS = HCS + DIST*M5(IG1,IQ1,IQ2)
805
806 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IGL,0,IG3,0,3124,25,*9)
807
808 END DO
809
810 END IF
811
812c _
813
814c du (-)
815
816 ID1 = IQ2 * 2 - 1
817
818 ID2 = IQ1 * 2 + 6
819
820 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
821
822 DIST = DISF(ID1,1)*DISF(ID2,2)
823
824 DO IG1 = 1,2
825
826 IG3 = ICH+IG1+2
827
828 DO IG2 = 1,4
829
830 IG4 = INU+IG2
831
832 HCS = HCS + DIST*M3(IG1,IG2,IQ1,IQ2)
833
834 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG4,0,IG3,0,2134,23,*9)
835
836 END DO
837
838 HCS = HCS + DIST*M5(IG1,IQ1,IQ2)
839
840 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IGL,0,IG3,0,2314,25,*9)
841
842 END DO
843
844 END IF
845
846c _
847
848c ud (-)
849
850 ID1 = IQ1 * 2 + 6
851
852 ID2 = IQ2 * 2 - 1
853
854 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
855
856 DIST = DISF(ID1,1)*DISF(ID2,2)
857
858 DO IG1 = 1,2
859
860 IG3 = ICH+IG1+2
861
862 DO IG2 = 1,4
863
864 IG4 = INU+IG2
865
866 HCS = HCS + DIST*M3(IG1,IG2,IQ1,IQ2)
867
868 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IG4,0,2134,23,*9)
869
870 END DO
871
872 HCS = HCS + DIST*M5(IG1,IQ1,IQ2)
873
874 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IGL,0,4132,25,*9)
875
876 END DO
877
878 END IF
879
880 END IF
881
882 END DO
883
884 2 CONTINUE
885
886c _ _ ~o ~ ~+-~
887
888c g q , g q , q g , q g -> X q , X q'
889
890c i LR a L
891
892c neutralino
893
894 DO IQ1 = 1,6
895
896c
897
898c gq
899
900 ID1 = 13
901
902 ID2 = IQ1
903
904 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
905
906 DIST = DISF(ID1,1)*DISF(ID2,2)
907
908 DO IG1 = 1,4
909
910 IG3 = INU+IG1
911
912 HCS = HCS + DIST*M6L(IG1,IQ1)
913
914 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,ID2,0,2431,26,*9)
915
916 HCS = HCS + DIST*M6R(IG1,IQ1)
917
918 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,ID2,2,2431,26,*9)
919
920 END DO
921
922 END IF
923
924c _
925
926c gq
927
928 ID1 = 13
929
930 ID2 = IQ1 + 6
931
932 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
933
934 DIST = DISF(ID1,1)*DISF(ID2,2)
935
936 DO IG1 = 1,4
937
938 IG3 = INU+IG1
939
940 HCS = HCS + DIST*M6L(IG1,IQ1)
941
942 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,ID2,0,4132,26,*9)
943
944 HCS = HCS + DIST*M6R(IG1,IQ1)
945
946 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,ID2,2,4132,26,*9)
947
948 END DO
949
950 END IF
951
952c
953
954c qg
955
956 ID1 = IQ1
957
958 ID2 = 13
959
960 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
961
962 DIST = DISF(ID1,1)*DISF(ID2,2)
963
964 DO IG1 = 1,4
965
966 IG3 = INU+IG1
967
968 HCS = HCS + DIST*M6L(IG1,IQ1)
969
970 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(ID1,0,IG3,0,3124,26,*9)
971
972 HCS = HCS + DIST*M6R(IG1,IQ1)
973
974 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(ID1,2,IG3,0,3124,26,*9)
975
976 END DO
977
978 END IF
979
980c _
981
982c qg
983
984 ID1 = IQ1 + 6
985
986 ID2 = 13
987
988 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
989
990 DIST = DISF(ID1,1)*DISF(ID2,2)
991
992 DO IG1 = 1,4
993
994 IG3 = INU+IG1
995
996 HCS = HCS + DIST*M6L(IG1,IQ1)
997
998 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(ID1,0,IG3,0,2314,26,*9)
999
1000 HCS = HCS + DIST*M6R(IG1,IQ1)
1001
1002 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(ID1,2,IG3,0,2314,26,*9)
1003
1004 END DO
1005
1006 END IF
1007
1008 END DO
1009
1010c chargino
1011
1012 DO IQ1 = 1,3
1013
1014 IQ3 = IQ1*2
1015
1016 DO 3 IQ2 = 1,3
1017
1018 IF (VCKM(IQ1,IQ2).LT.EPS) GOTO 3
1019
1020 IQ4 = IQ2*2-1
1021
1022 DO IG1 = 1,2
1023
1024 IG3 = ICH+IG1
1025
1026 IG4 = ICH+IG1+2
1027
1028c
1029
1030c gq & qg
1031
1032 ID1 = 13
1033
1034 ID2 = IQ3
1035
1036 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ3,IQ4)
1037
1038 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IQ4,0,2431,27,*9)
1039
1040 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ3,IQ4)
1041
1042 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IQ4,2,2431,27,*9)
1043
1044 ID2 = IQ4
1045
1046 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ4,IQ3)
1047
1048 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG4,0,IQ3,0,2431,27,*9)
1049
1050 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ4,IQ3)
1051
1052 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG4,0,IQ3,2,2431,27,*9)
1053
1054 ID1 = IQ3
1055
1056 ID2 = 13
1057
1058 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ3,IQ4)
1059
1060 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ4,0,IG3,0,3124,27,*9)
1061
1062 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ3,IQ4)
1063
1064 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ4,2,IG3,0,3124,27,*9)
1065
1066 ID1 = IQ4
1067
1068 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ4,IQ3)
1069
1070 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ3,0,IG4,0,3124,27,*9)
1071
1072 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ4,IQ3)
1073
1074 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ3,2,IG4,0,3124,27,*9)
1075
1076c _ _
1077
1078c gq & qg
1079
1080 ID1 = 13
1081
1082 ID2 = IQ3 + 6
1083
1084 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ3,IQ4)
1085
1086 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG4,0,IQ4,1,4132,27,*9)
1087
1088 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ3,IQ4)
1089
1090 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG4,0,IQ4,3,4132,27,*9)
1091
1092 ID2 = IQ4 + 6
1093
1094 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ4,IQ3)
1095
1096 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IQ3,1,4132,27,*9)
1097
1098 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ4,IQ3)
1099
1100 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IG3,0,IQ3,3,4132,27,*9)
1101
1102 ID1 = IQ3 + 6
1103
1104 ID2 = 13
1105
1106 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ3,IQ4)
1107
1108 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ4,1,IG4,0,2314,27,*9)
1109
1110 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ3,IQ4)
1111
1112 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ4,3,IG4,0,2314,27,*9)
1113
1114 ID1 = IQ4 + 6
1115
1116 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(1,IG1,IQ4,IQ3)
1117
1118 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ3,1,IG3,0,2314,27,*9)
1119
1120 HCS = HCS + DISF(ID1,1)*DISF(ID2,2)*M7(2,IG1,IQ4,IQ3)
1121
1122 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IQ3,3,IG3,0,2314,27,*9)
1123
1124 END DO
1125
1126 3 CONTINUE
1127
1128 END DO
1129
1130 EVWGT = HCS
1131
1132 RETURN
1133
1134C---GENERATE EVENT
1135
1136 9 IDN(1)=ID1
1137
1138 IDN(2)=ID2
1139
1140 IDCMF=15
1141
1142 CALL HWETWO
1143
1144 IF (AZSPIN) THEN
1145
1146C Calculate coefficients for constructing spin density matrices
1147
1148C Set to zero for now
1149
1150 CALL HWVZRO(7,GCOEF)
1151
1152 END IF
1153
1154 888 END
1155
1156CDECK ID>, HWHSSL.
1157
1158*CMZ :- -18/05/99 20.33.45 by Kosuke Odagiri
1159
1160*-- Author : Kosuke Odagiri
1161
1162C-----------------------------------------------------------------------
1163
1164 SUBROUTINE HWHSSL
1165
1166C-----------------------------------------------------------------------
1167
1168C SUSY 2 PARTON -> 2 SLEPTON PROCESSES
1169
1170C-----------------------------------------------------------------------
1171
1172 INCLUDE 'HERWIG61.INC'
1173
1174 DOUBLE PRECISION HWR, HWUAEM, EPS, HCS, RCS, DIST, S, PF, QPE,
1175
1176 & FACTR, SN2TH, MZ, MW, ME2(2,2,6,2), ME2W(2,3), EMSC2, GW2
1177
1178 INTEGER IQ, IQ1, IQ2, ID1, ID2, IL, IL1, IL2, I, J
1179
1180 EXTERNAL HWR, HWUAEM
1181
1182 SAVE HCS, ME2, ME2W
1183
1184 PARAMETER (EPS = 1.D-9)
1185
1186 DOUBLE COMPLEX Z, GZ, A, BL, BR, CL, CR, D, E
1187
1188 PARAMETER (Z = (0.D0,1.D0))
1189
1190 EQUIVALENCE (MZ, RMASS(200)), (MW, RMASS(198))
1191
1192C
1193
1194 S = XX(1)*XX(2)*PHEP(5,3)**2
1195
1196 EMSC2 = S
1197
1198 EMSCA = SQRT(EMSC2)
1199
1200 CALL HWSGEN(.FALSE.)
1201
1202 IF (GENEV) THEN
1203
1204 RCS = HCS*HWR()
1205
1206 ELSE
1207
1208 SN2TH = 0.25D0 - 0.25D0*COSTH**2
1209
1210 FACTR = FACTSS*HWUAEM(EMSC2)**2/CAFAC*SN2TH
1211
1212 GZ = (S-MZ**2+Z*S*GAMZ/MZ)/S
1213
1214 GW2 = ((ONE-MW**2/S)**2+(GAMW/MW)**2)*(TWO*SWEIN)**2
1215
1216c _ ~ ~*
1217
1218c q q -> l l
1219
1220c
1221
1222 DO IL = 1,6
1223
1224 DO I = 1,2
1225
1226 DO J = 1,2
1227
1228 IF (((I.NE.J).AND.(IL.NE.5)).OR.
1229
1230 & ((I.EQ.2).AND.(((IL/2)*2).EQ.IL))) THEN
1231
1232 QPE = -1.
1233
1234 ELSE
1235
1236 ID1 = 412 + I*12 + IL
1237
1238 ID2 = 412 + J*12 + IL
1239
1240 IL1 = IL + 10
1241
1242 QPE = S-(RMASS(ID1)+RMASS(ID2))**2
1243
1244 END IF
1245
1246 IF (QPE.GT.ZERO) THEN
1247
1248 PF = SQRT(QPE*(S-(RMASS(ID1)-RMASS(ID2))**2))/S
1249
1250 DO IQ = 1,2
1251
1252 A = QFCH(IL1)*QFCH(IQ)
1253
1254 BL = LFCH(IL1)/GZ
1255
1256 BR = RFCH(IL1)/GZ
1257
1258 CL = LMIXSS(IL,1,I)*LMIXSS(IL,1,J)
1259
1260 CR = LMIXSS(IL,2,I)*LMIXSS(IL,2,J)
1261
1262 D = (A+BL*LFCH(IQ))*CL+(A+BR*LFCH(IQ))*CR
1263
1264 E = (A+BL*RFCH(IQ))*CL+(A+BR*RFCH(IQ))*CR
1265
1266 ME2(I,J,IL,IQ)=FACTR*PF**3*DBLE(CONJG(D)*D+CONJG(E)*E)
1267
1268 END DO
1269
1270 ELSE
1271
1272 ME2(I,J,IL,1)=ZERO
1273
1274 ME2(I,J,IL,2)=ZERO
1275
1276 END IF
1277
1278 END DO
1279
1280 END DO
1281
1282 END DO
1283
1284c _ ~ ~*
1285
1286c q q' -> l v
1287
1288c
1289
1290 DO IL = 1,3
1291
1292 DO I = 1,2
1293
1294 IF ((IL.NE.3).AND.(I.EQ.2)) THEN
1295
1296 QPE = -1.
1297
1298 ELSE
1299
1300 ID1 = 411 + IL*2 + I*12
1301
1302 ID2 = 424 + IL*2
1303
1304 QPE = S-(RMASS(ID1)+RMASS(ID2))**2
1305
1306 END IF
1307
1308 IF (QPE.GT.ZERO) THEN
1309
1310 PF = SQRT(QPE*(S-(RMASS(ID1)-RMASS(ID2))**2))/S
1311
1312 ME2W(I,IL)=FACTR*PF**3/GW2
1313
1314 IF (IL.EQ.3) ME2W(I,3)=ME2W(I,3)*LMIXSS(5,1,I)**2
1315
1316 ELSE
1317
1318 ME2W(I,IL)=ZERO
1319
1320 END IF
1321
1322 END DO
1323
1324 END DO
1325
1326
1327
1328 END IF
1329
1330 HCS = 0.
1331
1332C
1333
1334 DO 1 ID1 = 1, 12
1335
1336 IF (DISF(ID1,1).LT.EPS) GOTO 1
1337
1338 IF (ID1.GT.6) THEN
1339
1340 ID2 = ID1 - 6
1341
1342 ELSE
1343
1344 ID2 = ID1 + 6
1345
1346 END IF
1347
1348 IQ = ID1 - ((ID1-1)/2)*2
1349
1350 IF (DISF(ID2,2).LT.EPS) GOTO 1
1351
1352 DIST = DISF(ID1,1)*DISF(ID2,2)
1353
1354 DO IL = 1,6
1355
1356 DO I = 1,2
1357
1358 DO J = 1,2
1359
1360 IL1 = IL+I*12
1361
1362 IL2 = IL+J*12
1363
1364 HCS = HCS + DIST*ME2(I,J,IL,IQ)
1365
1366 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IL1,2,IL2,3,2134,30,*9)
1367
1368 END DO
1369
1370 END DO
1371
1372 END DO
1373
1374 1 CONTINUE
1375
1376c _ _ _ _
1377
1378c ud(+), ud(-), du(-), du(+)
1379
1380 DO 2 IQ1 = 1, 3
1381
1382 DO IQ2 = 1, 3
1383
1384 IF(VCKM(IQ1,IQ2).GT.EPS) THEN
1385
1386c _
1387
1388c ud (+)
1389
1390 ID1 = IQ1 * 2
1391
1392 ID2 = IQ2 * 2 + 5
1393
1394 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
1395
1396 DIST = VCKM(IQ1,IQ2)*DISF(ID1,1)*DISF(ID2,2)
1397
1398 DO IL = 1,3
1399
1400 IL1 = IL*2-1
1401
1402 IL2 = IL1+1
1403
1404 HCS = HCS + DIST*ME2W(1,IL)
1405
1406 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IL1,5,IL2,4,2134,30,*9)
1407
1408 END DO
1409
1410 HCS = HCS + DIST*ME2W(2,3)
1411
1412 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(5,7,6,4,2134,30,*9)
1413
1414 END IF
1415
1416c _
1417
1418c du (+)
1419
1420 ID1 = IQ2 * 2 + 5
1421
1422 ID2 = IQ1 * 2
1423
1424 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
1425
1426 DIST = VCKM(IQ1,IQ2)*DISF(ID1,1)*DISF(ID2,2)
1427
1428 DO IL = 1,3
1429
1430 IL1 = IL*2-1
1431
1432 IL2 = IL1+1
1433
1434 HCS = HCS + DIST*ME2W(1,IL)
1435
1436 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IL1,5,IL2,4,2134,30,*9)
1437
1438 END DO
1439
1440 HCS = HCS + DIST*ME2W(2,3)
1441
1442 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(5,7,6,4,2134,30,*9)
1443
1444 END IF
1445
1446c _
1447
1448c du (-)
1449
1450 ID1 = IQ2 * 2 - 1
1451
1452 ID2 = IQ1 * 2 + 6
1453
1454 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
1455
1456 DIST = VCKM(IQ1,IQ2)*DISF(ID1,1)*DISF(ID2,2)
1457
1458 DO IL = 1,3
1459
1460 IL1 = IL*2-1
1461
1462 IL2 = IL1+1
1463
1464 HCS = HCS + DIST*ME2W(1,IL)
1465
1466 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IL1,4,IL2,5,2134,30,*9)
1467
1468 END DO
1469
1470 HCS = HCS + DIST*ME2W(2,3)
1471
1472 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(5,6,6,5,2134,30,*9)
1473
1474 END IF
1475
1476c _
1477
1478c ud (-)
1479
1480 ID1 = IQ1 * 2 + 6
1481
1482 ID2 = IQ2 * 2 - 1
1483
1484 IF ((DISF(ID1,1).GT.EPS).AND.(DISF(ID2,2).GT.EPS)) THEN
1485
1486 DIST = VCKM(IQ1,IQ2)*DISF(ID1,1)*DISF(ID2,2)
1487
1488 DO IL = 1,3
1489
1490 IL1 = IL*2-1
1491
1492 IL2 = IL1+1
1493
1494 HCS = HCS + DIST*ME2W(1,IL)
1495
1496 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(IL1,4,IL2,5,2134,30,*9)
1497
1498 END DO
1499
1500 HCS = HCS + DIST*ME2W(2,3)
1501
1502 IF (GENEV.AND.HCS.GT.RCS) CALL HWHSSS(5,6,6,5,2134,30,*9)
1503
1504 END IF
1505
1506 END IF
1507
1508 END DO
1509
1510 2 CONTINUE
1511
1512
1513
1514 EVWGT = HCS
1515
1516 RETURN
1517
1518C---GENERATE EVENT
1519
1520 9 IDN(1)=ID1
1521
1522 IDN(2)=ID2
1523
1524 IDCMF=15
1525
1526 CALL HWETWO
1527
1528 IF (AZSPIN) THEN
1529
1530C Calculate coefficients for constructing spin density matrices
1531
1532C Set to zero for now
1533
1534 CALL HWVZRO(7,GCOEF)
1535
1536 END IF
1537
1538 END