Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pysfdc.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pysfdc.f
1 
2 C*********************************************************************
3 
4 C...PYSFDC
5 C...Calculates decays of sfermions.
6 
7  SUBROUTINE pysfdc(KFIN,XLAM,IDLAM,IKNT)
8 
9 C...Double precision and integer declarations.
10  IMPLICIT DOUBLE PRECISION(a-h, o-z)
11  IMPLICIT INTEGER(i-n)
12  INTEGER pyk,pychge,pycomp
13 C...Parameter statement to help give large particle numbers.
14  parameter(ksusy1=1000000,ksusy2=2000000,ktechn=3000000,
15  &kexcit=4000000,kdimen=5000000)
16 C...Commonblocks.
17  common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
18  common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
19  common/pymssm/imss(0:99),rmss(0:99)
20  common/pyssmt/zmix(4,4),umix(2,2),vmix(2,2),smz(4),smw(2),
21  &sfmix(16,4),zmixi(4,4),umixi(2,2),vmixi(2,2)
22  SAVE /pydat1/,/pydat2/,/pymssm/,/pyssmt/
23 
24 C...Local variables.
25  COMPLEX*16 zmixc(4,4),vmixc(2,2),umixc(2,2)
26  COMPLEX*16 cal,car,cbl,cbr,calp,carp,cblp,cbrp,ca,cb
27  INTEGER kfin,kcin
28  DOUBLE PRECISION xmi,xmj,xmf,xmsf1,xmsf2,xmw,xmw2,xmz,axmj
29  DOUBLE PRECISION xmi2,xmi3,xma2,xmb2,xmfp
30  DOUBLE PRECISION pylamf,xl
31  DOUBLE PRECISION tanw,xw,aem,c1,as
32  DOUBLE PRECISION al,ar,bl,br
33  DOUBLE PRECISION ch1,ch2,ch3,ch4
34  DOUBLE PRECISION xmbot,xmtop
35  DOUBLE PRECISION xlam(0:400)
36  INTEGER idlam(400,3)
37  INTEGER lknt,ix,ilr,idu,j,i,iknt,ifl,ii
38  DOUBLE PRECISION sr2
39  DOUBLE PRECISION cbeta,sbeta
40  DOUBLE PRECISION cw
41  DOUBLE PRECISION beta,alfa,xmu,at,ab,atrit,atrib,atril
42  DOUBLE PRECISION cosa,sina,tanb
43  DOUBLE PRECISION pyalem,pi,pyalps,ei
44  DOUBLE PRECISION ghrr,ghll,ghlr,xmb,blr
45  INTEGER ig,kf1,kf2
46  INTEGER igg(4),kfnchi(4),kfcchi(2)
47  DATA igg/23,25,35,36/
48  DATA pi/3.141592654d0/
49  DATA sr2/1.4142136d0/
50  DATA kfnchi/1000022,1000023,1000025,1000035/
51  DATA kfcchi/1000024,1000037/
52 
53 C...COUNT THE NUMBER OF DECAY MODES
54  lknt=0
55 
56 C...NO NU_R DECAYS
57  IF(kfin.EQ.ksusy2+12.OR.kfin.EQ.ksusy2+14.OR.
58  &kfin.EQ.ksusy2+16) RETURN
59 
60  xmw=pmas(24,1)
61  xmw2=xmw**2
62  xmz=pmas(23,1)
63  xw=paru(102)
64  tanw = sqrt(xw/(1d0-xw))
65  cw=sqrt(1d0-xw)
66 
67  DO 110 i=1,4
68  DO 100 j=1,4
69  zmixc(j,i)=dcmplx(zmix(j,i),zmixi(j,i))
70  100 CONTINUE
71  110 CONTINUE
72  DO 130 i=1,2
73  DO 120 j=1,2
74  vmixc(j,i)=dcmplx(vmix(j,i),vmixi(j,i))
75  umixc(j,i)=dcmplx(umix(j,i),umixi(j,i))
76  120 CONTINUE
77  130 CONTINUE
78 
79 C...KCIN
80  kcin=pycomp(kfin)
81 C...ILR is 1 for left and 2 for right.
82  ilr=kfin/ksusy1
83 C...IFL is matching non-SUSY flavour.
84  ifl=mod(kfin,ksusy1)
85 C...IDU is weak isospin, 1 for down and 2 for up.
86  idu=2-mod(ifl,2)
87 
88  xmi=pmas(kcin,1)
89  xmi2=xmi**2
90  aem=pyalem(xmi2)
91  as =pyalps(xmi2)
92  c1=aem/xw
93  xmi3=xmi**3
94  ei=kchg(ifl,1)/3d0
95 
96  xmbot=pymrun(5,xmi2)
97  xmtop=pymrun(6,xmi2)
98 
99  tanb=rmss(5)
100  beta=atan(tanb)
101  alfa=rmss(18)
102  cbeta=cos(beta)
103  sbeta=tanb*cbeta
104  sina=sin(alfa)
105  cosa=cos(alfa)
106  xmu=-rmss(4)
107  atrit=rmss(16)
108  atrib=rmss(15)
109  atril=rmss(17)
110 
111 C...2-BODY DECAYS OF SFERMION -> GRAVITINO + FERMION
112 
113  IF(imss(11).EQ.1) THEN
114  xmp=rmss(29)
115  idg=39+ksusy1
116  xmgr=pmas(pycomp(idg),1)
117  xfac=(xmi2/(xmp*xmgr))**2*xmi/48d0/pi
118  IF(ifl.EQ.5) THEN
119  xmf=xmbot
120  ELSEIF(ifl.EQ.6) THEN
121  xmf=xmtop
122  ELSE
123  xmf=pmas(ifl,1)
124  ENDIF
125  IF(xmi.GT.xmgr+xmf) THEN
126  lknt=lknt+1
127  idlam(lknt,1)=idg
128  idlam(lknt,2)=ifl
129  idlam(lknt,3)=0
130  xlam(lknt)=xfac*(1d0-xmf**2/xmi2)**4
131  ENDIF
132  ENDIF
133 
134 C...2-BODY DECAYS OF SFERMION -> FERMION + GAUGE/GAUGINO
135 
136 C...CHARGED DECAYS:
137  DO 140 ix=1,2
138 C...DI -> U CHI1-,CHI2-
139  IF(idu.EQ.1) THEN
140  xmfp=pmas(ifl+1,1)
141  xmf =pmas(ifl,1)
142 C...UI -> D CHI1+,CHI2+
143  ELSE
144  xmfp=pmas(ifl-1,1)
145  xmf =pmas(ifl,1)
146  ENDIF
147  xmj=smw(ix)
148  axmj=abs(xmj)
149  IF(xmi.GE.axmj+xmfp) THEN
150  xma2=xmj**2
151  xmb2=xmfp**2
152  IF(idu.EQ.2) THEN
153  IF(ifl.EQ.6) THEN
154  xmfp=xmbot
155  xmf =xmtop
156  ELSEIF(ifl.LT.6) THEN
157  xmf=0d0
158  xmfp=0d0
159  ENDIF
160  cbl=vmixc(ix,1)
161  cal=-xmfp*umixc(ix,2)/sr2/xmw/cbeta
162  cbr=-xmf*vmixc(ix,2)/sr2/xmw/sbeta
163  car=0d0
164  ELSE
165  IF(ifl.EQ.5) THEN
166  xmf =xmbot
167  xmfp=xmtop
168  ELSEIF(ifl.LT.5) THEN
169  xmf=0d0
170  xmfp=0d0
171  ENDIF
172  cbl=umixc(ix,1)
173  cal=-xmfp*vmixc(ix,2)/sr2/xmw/sbeta
174  cbr=-xmf*umixc(ix,2)/sr2/xmw/cbeta
175  car=0d0
176  ENDIF
177 
178  calp=sfmix(ifl,1)*cal + sfmix(ifl,2)*car
179  cblp=sfmix(ifl,1)*cbl + sfmix(ifl,2)*cbr
180  carp=sfmix(ifl,4)*car + sfmix(ifl,3)*cal
181  cbrp=sfmix(ifl,4)*cbr + sfmix(ifl,3)*cbl
182  cal=calp
183  cbl=cblp
184  car=carp
185  cbr=cbrp
186 
187 C...F1 -> F` CHI
188  IF(ilr.EQ.1) THEN
189  ca=cal
190  cb=cbl
191 C...F2 -> F` CHI
192  ELSE
193  ca=car
194  cb=cbr
195  ENDIF
196  lknt=lknt+1
197  xl=pylamf(xmi2,xma2,xmb2)
198 C...SPIN AVERAGE = 1/1 NOT 1/2....NO COLOR ENHANCEMENT
199  xlam(lknt)=2d0*c1/8d0/xmi3*sqrt(xl)*((xmi2-xmb2-xma2)*
200  & (abs(ca)**2+abs(cb)**2)-4d0*dble(ca*dconjg(cb))*xmj*xmfp)
201  idlam(lknt,3)=0
202  IF(idu.EQ.1) THEN
203  idlam(lknt,1)=-kfcchi(ix)
204  idlam(lknt,2)=ifl+1
205  ELSE
206  idlam(lknt,1)=kfcchi(ix)
207  idlam(lknt,2)=ifl-1
208  ENDIF
209  ENDIF
210  140 CONTINUE
211 
212 C...NEUTRAL DECAYS
213  DO 150 ix=1,4
214 C...DI -> D CHI10
215  xmf=pmas(ifl,1)
216  xmj=smz(ix)
217  axmj=abs(xmj)
218  IF(xmi.GE.axmj+xmf) THEN
219  xma2=xmj**2
220  xmb2=xmf**2
221  IF(idu.EQ.1) THEN
222  IF(ifl.EQ.5) THEN
223  xmf=xmbot
224  ELSEIF(ifl.LT.5) THEN
225  xmf=0d0
226  ENDIF
227  cbl=-zmixc(ix,2)+tanw*zmixc(ix,1)*(2d0*ei+1)
228  cal=xmf*zmixc(ix,3)/xmw/cbeta
229  car=-2d0*ei*tanw*zmixc(ix,1)
230  cbr=cal
231  ELSE
232  IF(ifl.EQ.6) THEN
233  xmf=xmtop
234  ELSEIF(ifl.LT.5) THEN
235  xmf=0d0
236  ENDIF
237  cbl=zmixc(ix,2)+tanw*zmixc(ix,1)*(2d0*ei-1)
238  cal=xmf*zmixc(ix,4)/xmw/sbeta
239  car=-2d0*ei*tanw*zmixc(ix,1)
240  cbr=cal
241  ENDIF
242 
243  calp=sfmix(ifl,1)*cal + sfmix(ifl,2)*car
244  cblp=sfmix(ifl,1)*cbl + sfmix(ifl,2)*cbr
245  carp=sfmix(ifl,4)*car + sfmix(ifl,3)*cal
246  cbrp=sfmix(ifl,4)*cbr + sfmix(ifl,3)*cbl
247  cal=calp
248  cbl=cblp
249  car=carp
250  cbr=cbrp
251 
252 C...F1 -> F CHI
253  IF(ilr.EQ.1) THEN
254  ca=cal
255  cb=cbl
256 C...F2 -> F CHI
257  ELSE
258  ca=car
259  cb=cbr
260  ENDIF
261  lknt=lknt+1
262  xl=pylamf(xmi2,xma2,xmb2)
263 C...SPIN AVERAGE = 1/1 NOT 1/2....NO COLOR ENHANCEMENT
264  xlam(lknt)=c1/8d0/xmi3*sqrt(xl)*((xmi2-xmb2-xma2)*
265  & (abs(ca)**2+abs(cb)**2)-4d0*dble(ca*dconjg(cb))*xmj*xmf)
266  idlam(lknt,1)=kfnchi(ix)
267  idlam(lknt,2)=ifl
268  idlam(lknt,3)=0
269  ENDIF
270  150 CONTINUE
271 
272 C...2-BODY DECAYS TO SM GAUGE AND HIGGS BOSONS
273 C...IG=23,25,35,36
274  DO 160 ii=1,4
275  ig=igg(ii)
276  IF(ilr.EQ.1) goto 160
277  xmb=pmas(ig,1)
278  xmsf1=pmas(pycomp(kfin-ksusy1),1)
279  IF(xmi.LT.xmsf1+xmb) goto 160
280  IF(ig.EQ.23) THEN
281  bl=-sign(.5d0,ei)/cw+ei*xw/cw
282  br=ei*xw/cw
283  blr=0d0
284  ELSEIF(ig.EQ.25) THEN
285  IF(ifl.EQ.5) THEN
286  xmf=xmbot
287  ELSEIF(ifl.EQ.6) THEN
288  xmf=xmtop
289  ELSEIF(ifl.LT.5) THEN
290  xmf=0d0
291  ELSE
292  xmf=pmas(ifl,1)
293  ENDIF
294  IF(idu.EQ.2) THEN
295  ghll=xmz/cw*(0.5d0-ei*xw)*(-sin(alfa+beta))+
296  & xmf**2/xmw*cosa/sbeta
297  ghrr=xmz/cw*(ei*xw)*(-sin(alfa+beta))+
298  & xmf**2/xmw*cosa/sbeta
299  ELSE
300  ghll=xmz/cw*(0.5d0-ei*xw)*(-sin(alfa+beta))+
301  & xmf**2/xmw*(-sina)/cbeta
302  ghrr=xmz/cw*(ei*xw)*(-sin(alfa+beta))+
303  & xmf**2/xmw*(-sina)/cbeta
304  ENDIF
305  IF(ifl.EQ.5) THEN
306  at=atrib
307  ELSEIF(ifl.EQ.6) THEN
308  at=atrit
309  ELSEIF(ifl.EQ.15) THEN
310  at=atril
311  ELSE
312  at=0d0
313  ENDIF
314 C.........need to complexify
315  IF(idu.EQ.2) THEN
316  ghlr=xmf/2d0/xmw/sbeta*(-xmu*sina+
317  & at*cosa)
318  ELSE
319  ghlr=xmf/2d0/xmw/cbeta*(xmu*cosa-
320  & at*sina)
321  ENDIF
322  bl=ghll
323  br=ghrr
324  blr=-ghlr
325  ELSEIF(ig.EQ.35) THEN
326  IF(ifl.EQ.5) THEN
327  xmf=xmbot
328  ELSEIF(ifl.EQ.6) THEN
329  xmf=xmtop
330  ELSEIF(ifl.LT.5) THEN
331  xmf=0d0
332  ELSE
333  xmf=pmas(ifl,1)
334  ENDIF
335  IF(idu.EQ.2) THEN
336  ghll=xmz/cw*(0.5d0-ei*xw)*cos(alfa+beta)+
337  & xmf**2/xmw*sina/sbeta
338  ghrr=xmz/cw*(ei*xw)*cos(alfa+beta)+
339  & xmf**2/xmw*sina/sbeta
340  ELSE
341  ghll=xmz/cw*(0.5d0-ei*xw)*cos(alfa+beta)+
342  & xmf**2/xmw*cosa/cbeta
343  ghrr=xmz/cw*(ei*xw)*cos(alfa+beta)+
344  & xmf**2/xmw*cosa/cbeta
345  ENDIF
346  IF(ifl.EQ.5) THEN
347  at=atrib
348  ELSEIF(ifl.EQ.6) THEN
349  at=atrit
350  ELSEIF(ifl.EQ.15) THEN
351  at=atril
352  ELSE
353  at=0d0
354  ENDIF
355 C.........Need to complexify
356  IF(idu.EQ.2) THEN
357  ghlr=xmf/2d0/xmw/sbeta*(xmu*cosa+
358  & at*sina)
359  ELSE
360  ghlr=xmf/2d0/xmw/cbeta*(xmu*sina+
361  & at*cosa)
362  ENDIF
363  bl=ghll
364  br=ghrr
365  blr=ghlr
366  ELSEIF(ig.EQ.36) THEN
367  ghll=0d0
368  ghrr=0d0
369  IF(ifl.EQ.5) THEN
370  xmf=xmbot
371  ELSEIF(ifl.EQ.6) THEN
372  xmf=xmtop
373  ELSEIF(ifl.LT.5) THEN
374  xmf=0d0
375  ELSE
376  xmf=pmas(ifl,1)
377  ENDIF
378  IF(ifl.EQ.5) THEN
379  at=atrib
380  ELSEIF(ifl.EQ.6) THEN
381  at=atrit
382  ELSEIF(ifl.EQ.15) THEN
383  at=atril
384  ELSE
385  at=0d0
386  ENDIF
387 C.........Need to complexify
388  IF(idu.EQ.2) THEN
389  ghlr=xmf/2d0/xmw*(-xmu+at/tanb)
390  ELSE
391  ghlr=xmf/2d0/xmw/(-xmu+at*tanb)
392  ENDIF
393  bl=ghll
394  br=ghrr
395  blr=ghlr
396  ENDIF
397  al=sfmix(ifl,1)*sfmix(ifl,3)*bl+
398  & sfmix(ifl,2)*sfmix(ifl,4)*br+
399  & (sfmix(ifl,1)*sfmix(ifl,4)+sfmix(ifl,3)*sfmix(ifl,2))*blr
400  xl=pylamf(xmi2,xmsf1**2,xmb**2)
401  lknt=lknt+1
402  IF(ig.EQ.23) THEN
403  xlam(lknt)=c1/4d0/xmi3*xl**1.5d0/xmb**2*al**2
404  ELSE
405  xlam(lknt)=c1/4d0/xmi3*sqrt(xl)*al**2
406  ENDIF
407  idlam(lknt,3)=0
408  idlam(lknt,1)=kfin-ksusy1
409  idlam(lknt,2)=ig
410  160 CONTINUE
411 
412 C...SF -> SF' + W
413  xmb=pmas(24,1)
414  IF(mod(ifl,2).EQ.0) THEN
415  kf1=ksusy1+ifl-1
416  ELSE
417  kf1=ksusy1+ifl+1
418  ENDIF
419  kf2=kf1+ksusy1
420  xmsf1=pmas(pycomp(kf1),1)
421  xmsf2=pmas(pycomp(kf2),1)
422  IF(xmi.GT.xmb+xmsf1) THEN
423  IF(mod(ifl,2).EQ.0) THEN
424  IF(ilr.EQ.1) THEN
425  al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl-1,1)
426  ELSE
427  al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl-1,1)
428  ENDIF
429  ELSE
430  IF(ilr.EQ.1) THEN
431  al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl+1,1)
432  ELSE
433  al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl+1,1)
434  ENDIF
435  ENDIF
436  xl=pylamf(xmi2,xmsf1**2,xmb**2)
437  lknt=lknt+1
438  xlam(lknt)=c1/4d0/xmi3*xl**1.5d0/xmb**2*al**2
439  idlam(lknt,3)=0
440  idlam(lknt,1)=kf1
441  idlam(lknt,2)=sign(24,kchg(ifl,1))
442  ENDIF
443  IF(xmi.GT.xmb+xmsf2) THEN
444  IF(mod(ifl,2).EQ.0) THEN
445  IF(ilr.EQ.1) THEN
446  al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl-1,3)
447  ELSE
448  al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl-1,3)
449  ENDIF
450  ELSE
451  IF(ilr.EQ.1) THEN
452  al=1d0/sr2*sfmix(ifl,1)*sfmix(ifl+1,3)
453  ELSE
454  al=1d0/sr2*sfmix(ifl,3)*sfmix(ifl+1,3)
455  ENDIF
456  ENDIF
457  xl=pylamf(xmi2,xmsf2**2,xmb**2)
458  lknt=lknt+1
459  xlam(lknt)=c1/4d0/xmi3*xl**1.5d0/xmb**2*al**2
460  idlam(lknt,3)=0
461  idlam(lknt,1)=kf2
462  idlam(lknt,2)=sign(24,kchg(ifl,1))
463  ENDIF
464 
465 C...SF -> SF' + HC
466  xmb=pmas(37,1)
467  IF(mod(ifl,2).EQ.0) THEN
468  kf1=ksusy1+ifl-1
469  ELSE
470  kf1=ksusy1+ifl+1
471  ENDIF
472  kf2=kf1+ksusy1
473  xmsf1=pmas(pycomp(kf1),1)
474  xmsf2=pmas(pycomp(kf2),1)
475  IF(xmi.GT.xmb+xmsf1) THEN
476  xmf=0d0
477  xmfp=0d0
478  at=0d0
479  ab=0d0
480  IF(mod(ifl,2).EQ.0) THEN
481 C...T1-> B1 HC
482  IF(ilr.EQ.1) THEN
483  ch1=-sfmix(ifl,1)*sfmix(ifl-1,1)
484  ch2= sfmix(ifl,2)*sfmix(ifl-1,2)
485  ch3=-sfmix(ifl,1)*sfmix(ifl-1,2)
486  ch4=-sfmix(ifl,2)*sfmix(ifl-1,1)
487 C...T2-> B1 HC
488  ELSE
489  ch1= sfmix(ifl,3)*sfmix(ifl-1,1)
490  ch2=-sfmix(ifl,4)*sfmix(ifl-1,2)
491  ch3= sfmix(ifl,3)*sfmix(ifl-1,2)
492  ch4= sfmix(ifl,4)*sfmix(ifl-1,1)
493  ENDIF
494  IF(ifl.EQ.6) THEN
495  xmf=xmtop
496  xmfp=xmbot
497  at=atrit
498  ab=atrib
499  ENDIF
500  ELSE
501 C...B1 -> T1 HC
502  IF(ilr.EQ.1) THEN
503  ch1=-sfmix(ifl+1,1)*sfmix(ifl,1)
504  ch2= sfmix(ifl+1,2)*sfmix(ifl,2)
505  ch3=-sfmix(ifl+1,1)*sfmix(ifl,2)
506  ch4=-sfmix(ifl+1,2)*sfmix(ifl,1)
507 C...B2-> T1 HC
508  ELSE
509  ch1= sfmix(ifl,3)*sfmix(ifl+1,1)
510  ch2=-sfmix(ifl,4)*sfmix(ifl+1,2)
511  ch3= sfmix(ifl,4)*sfmix(ifl+1,1)
512  ch4= sfmix(ifl,3)*sfmix(ifl+1,2)
513  ENDIF
514  IF(ifl.EQ.5) THEN
515  xmf=xmtop
516  xmfp=xmbot
517  at=atrit
518  ab=atrib
519  ENDIF
520  ENDIF
521  xl=pylamf(xmi2,xmsf1**2,xmb**2)
522  lknt=lknt+1
523 C.......Need to complexify
524  al=ch1*(xmw2*2d0*cbeta*sbeta-xmfp**2*tanb-xmf**2/tanb)+
525  & ch2*2d0*xmf*xmfp/(2d0*cbeta*sbeta)+
526  & ch3*xmfp*(-xmu+ab*tanb)+ch4*xmf*(-xmu+at/tanb)
527  xlam(lknt)=c1/8d0/xmi3*sqrt(xl)/xmw2*al**2
528  idlam(lknt,3)=0
529  idlam(lknt,1)=kf1
530  idlam(lknt,2)=sign(37,kchg(ifl,1))
531  ENDIF
532  IF(xmi.GT.xmb+xmsf2) THEN
533  xmf=0d0
534  xmfp=0d0
535  at=0d0
536  ab=0d0
537  IF(mod(ifl,2).EQ.0) THEN
538 C...T1-> B2 HC
539  IF(ilr.EQ.1) THEN
540  ch1= sfmix(ifl-1,3)*sfmix(ifl,1)
541  ch2=-sfmix(ifl-1,4)*sfmix(ifl,2)
542  ch3= sfmix(ifl-1,4)*sfmix(ifl,1)
543  ch4= sfmix(ifl-1,3)*sfmix(ifl,2)
544 C...T2-> B2 HC
545  ELSE
546  ch1= -sfmix(ifl,3)*sfmix(ifl-1,3)
547  ch2= sfmix(ifl,4)*sfmix(ifl-1,4)
548  ch3= -sfmix(ifl,3)*sfmix(ifl-1,4)
549  ch4= -sfmix(ifl,4)*sfmix(ifl-1,3)
550  ENDIF
551  IF(ifl.EQ.6) THEN
552  xmf=xmtop
553  xmfp=xmbot
554  at=atrit
555  ab=atrib
556  ENDIF
557  ELSE
558 C...B1 -> T2 HC
559  IF(ilr.EQ.1) THEN
560  ch1= sfmix(ifl+1,3)*sfmix(ifl,1)
561  ch2=-sfmix(ifl+1,4)*sfmix(ifl,2)
562  ch3= sfmix(ifl+1,3)*sfmix(ifl,2)
563  ch4= sfmix(ifl+1,4)*sfmix(ifl,1)
564 C...B2-> T2 HC
565  ELSE
566  ch1= -sfmix(ifl+1,3)*sfmix(ifl,3)
567  ch2= sfmix(ifl+1,4)*sfmix(ifl,4)
568  ch3= -sfmix(ifl+1,3)*sfmix(ifl,4)
569  ch4= -sfmix(ifl+1,4)*sfmix(ifl,3)
570  ENDIF
571  IF(ifl.EQ.5) THEN
572  xmf=xmtop
573  xmfp=xmbot
574  at=atrit
575  ab=atrib
576  ENDIF
577  ENDIF
578  xl=pylamf(xmi2,xmsf1**2,xmb**2)
579  lknt=lknt+1
580 C.......Need to complexify
581  al=ch1*(xmw2*2d0*cbeta*sbeta-xmfp**2*tanb-xmf**2/tanb)+
582  & ch2*2d0*xmf*xmfp/(2d0*cbeta*sbeta)+
583  & ch3*xmfp*(-xmu+ab*tanb)+ch4*xmf*(-xmu+at/tanb)
584  xlam(lknt)=c1/8d0/xmi3*sqrt(xl)/xmw2*al**2
585  idlam(lknt,3)=0
586  idlam(lknt,1)=kf2
587  idlam(lknt,2)=sign(37,kchg(ifl,1))
588  ENDIF
589 
590 C...2-BODY DECAYS OF SQUARK -> QUARK GLUINO
591 
592  IF(ifl.LE.6) THEN
593  xmfp=0d0
594  xmf=0d0
595  IF(ifl.EQ.6) xmf=pmas(6,1)
596  IF(ifl.EQ.5) xmf=pmas(5,1)
597  xmj=pmas(pycomp(ksusy1+21),1)
598  axmj=abs(xmj)
599  IF(xmi.GE.axmj+xmf) THEN
600  al=-sfmix(ifl,3)
601  bl=sfmix(ifl,1)
602  ar=-sfmix(ifl,4)
603  br=sfmix(ifl,2)
604 C...F1 -> F CHI
605  IF(ilr.EQ.1) THEN
606  xca=al
607  xcb=bl
608 C...F2 -> F CHI
609  ELSE
610  xca=ar
611  xcb=br
612  ENDIF
613  lknt=lknt+1
614  xma2=xmj**2
615  xmb2=xmf**2
616  xl=pylamf(xmi2,xma2,xmb2)
617  xlam(lknt)=4d0/3d0*as/2d0/xmi3*sqrt(xl)*((xmi2-xmb2-xma2)*
618  & (xca**2+xcb**2)+4d0*xca*xcb*xmj*xmf)
619  idlam(lknt,1)=ksusy1+21
620  idlam(lknt,2)=ifl
621  idlam(lknt,3)=0
622  ENDIF
623  ENDIF
624 
625 C...IF NOTHING ELSE FOR T1, THEN T1* -> C+CHI0
626  IF(kfin.EQ.ksusy1+6.AND.pmas(kcin,1).GT.
627  &pmas(pycomp(ksusy1+22),1)+pmas(4,1)) THEN
628 C...THIS IS A BACK-OF-THE-ENVELOPE ESTIMATE
629 C...M = 1/(16PI**2)G**3 = G*2/(4PI) G/(4PI) = C1 * G/(4PI)
630 C...M*M = C1**2 * G**2/(16PI**2)
631 C...G = 1/(8PI)P/MI**2 * M*M = C1**3/(32PI**2)*LAM/(2*MI**3)
632  lknt=lknt+1
633  xl=pylamf(xmi2,0d0,pmas(pycomp(ksusy1+22),1)**2)
634  xlam(lknt)=c1**3/64d0/pi**2/xmi3*sqrt(xl)
635  IF(xlam(lknt).EQ.0) xlam(lknt)=1d-3
636  idlam(lknt,1)=ksusy1+22
637  idlam(lknt,2)=4
638  idlam(lknt,3)=0
639  ENDIF
640 
641 C...R-violating sfermion decays (SKANDS).
642  CALL pyrvsf(kfin,xlam,idlam,lknt)
643 
644  iknt=lknt
645  xlam(0)=0d0
646  DO 170 i=1,iknt
647  IF(xlam(i).LT.0d0) xlam(i)=0d0
648  xlam(0)=xlam(0)+xlam(i)
649  170 CONTINUE
650  IF(xlam(0).EQ.0d0) xlam(0)=1d-3
651 
652  RETURN
653  END