* * $Id$ * * $Log$ * Revision 1.1 2000/06/19 20:00:31 eugenio * Initial revision * * Revision 1.1.1.1 1994/11/22 16:57:02 zfiles * first version of korb in CVS * * #include "sys/CLEO_machine.h" #include "pilot.h" *CMZ : 2.00/00 21/01/93 15.42.28 by Alan Weinstein *-- Author : SUBROUTINE JAKER(JAK) C ********************* C C ********************************************************************** C * C *********TAUOLA LIBRARY: VERSION 2.4 ******** * C **************NOVEMBER 1992****************** * C ** AUTHORS: S.JADACH, Z.WAS ***** * C ** R. DECKER, M. JEZABEK, J.H.KUEHN, ***** * C ********AVAILABLE FROM: WASM AT CERNVM ****** * C *******PUBLISHED IN COMP. PHYS. COMM.******** * C *** PREPRINT CERN-TH-5856 SEPTEMBER 1990 **** * C *** PREPRINT CERN-TH-6195 OCTOBER 1991 **** * C *** PREPRINT CERN-TH- NOVEMBER 1992 **** * C ********************************************************************** C ---------------------------------------------------------------------- c SUBROUTINE JAKER, C CHOOSES DECAY MODE ACCORDING TO LIST OF BRANCHING RATIOS C JAK=1 ELECTRON MODE C JAK=2 MUON MODE C JAK=3 PION MODE C JAK=4 RHO MODE C JAK=5 A1 MODE C JAK=6 K MODE C JAK=7 K* MODE C JAK=8 nPI MODE C C called by : DEXAY C ---------------------------------------------------------------------- COMMON / TAUBRA / GAMPRT(30),JLIST(30),NCHAN REAL CUMUL(20) C IF(NCHAN.LE.0.OR.NCHAN.GT.30) GOTO 902 CALL RANMAR(RRR,1) SUM=0 DO 20 I=1,NCHAN SUM=SUM+GAMPRT(I) 20 CUMUL(I)=SUM DO 25 I=NCHAN,1,-1 IF(RRR.LT.CUMUL(I)/CUMUL(NCHAN)) JI=I 25 CONTINUE JAK=JLIST(JI) RETURN 902 PRINT 9020 9020 FORMAT(' ----- JAKER: WRONG NCHAN') STOP END