* * $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 DEKAY(KTO,HX) C *********************** C THIS DEKAY IS IN SPIRIT OF THE 'DECAY' WHICH C WAS INCLUDED IN KORAL-B PROGRAM, COMP. PHYS. COMMUN. C VOL. 36 (1985) 191, SEE COMMENTS ON GENERAL PHILOSOPHY THERE. C KTO=0 INITIALISATION (OBLIGATORY) C KTO=1,11 DENOTES TAU+ AND KTO=2,12 TAU- C DEKAY(1,H) AND DEKAY(2,H) IS CALLED INTERNALLY BY MC GENERATOR. C H DENOTES THE POLARIMETRIC VECTOR, USED BY THE HOST PROGRAM FOR C CALCULATION OF THE SPIN WEIGHT. C USER MAY OPTIONALLY CALL DEKAY(11,H) DEKAY(12,H) IN ORDER C TO TRANSFORM DECAY PRODUCTS TO CMS AND WRITE LUND RECORD IN /LUJETS/. C KTO=100, PRINT FINAL REPORT (OPTIONAL). C DECAY MODES: C JAK=1 ELECTRON DECAY C JAK=2 MU DECAY C JAK=3 PI DECAY C JAK=4 RHO DECAY C JAK=5 A1 DECAY C JAK=6 K DECAY C JAK=7 K* DECAY C JAK=8 NPI DECAY C JAK=0 INCLUSIVE: JAK=1,2,3,4,5,6,7,8 REAL H(4) REAL*8 HX(4) COMMON / JAKI / JAK1,JAK2,JAKP,JAKM,KTOM COMMON / IDFC / IDF COMMON / TAUBMC / GAMPMC(30),GAMPER(30),NEVDEC(30) REAL*4 GAMPMC ,GAMPER COMMON / INOUT / INUT,IOUT REAL PDUM1(4),PDUM2(4),PDUM3(4),PDUM4(4),PDUM5(4),HDUM(4) REAL PDUMX(4,6) DATA IWARM/0/ KTOM=KTO IF(KTO.EQ.-1) THEN C ================== C INITIALISATION OR REINITIALISATION KTOM=1 IF (IWARM.EQ.1) X=5/(IWARM-1) IWARM=1 WRITE(IOUT,7001) JAK1,JAK2 NEVTOT=0 NEV1=0 NEV2=0 IF(JAK1.NE.-1.OR.JAK2.NE.-1) THEN CALL DADMEL(-1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5) CALL DADMMU(-1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5) CALL DADMPI(-1,IDUM,PDUM,PDUM1,PDUM2) CALL DADMRO(-1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4) CALL DADMAA(-1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5,JDUM) CALL DADMKK(-1,IDUM,PDUM,PDUM1,PDUM2) CALL DADMKS(-1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,JDUM) CALL DAD4PI(-1,IDUM,HDUM,PDUM1,PDUM2,PDUMX,JDUM) CALL DADNPI(-1,IDUM, PDUM1,PDUM2,PDUMX,JDUM) CALL DADMPK(-1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5,JDUM) ENDIF ELSEIF(KTO.EQ.1) THEN C ===================== C DECAY OF TAU+ IN THE TAU REST FRAME NEVTOT=NEVTOT+1 IF(IWARM.EQ.0) GOTO 902 ISGN= IDF/IABS(IDF) C AJWMOD to change BRs depending on sign: CALL TAURDF(KTO) C end AJWMOD CALL DEKAY1(0,H,ISGN) ELSEIF(KTO.EQ.2) THEN C ================================= C DECAY OF TAU- IN THE TAU REST FRAME NEVTOT=NEVTOT+1 IF(IWARM.EQ.0) GOTO 902 ISGN=-IDF/IABS(IDF) C AJWMOD to change BRs depending on sign: CALL TAURDF(KTO) C end AJWMOD CALL DEKAY2(0,H,ISGN) ELSEIF(KTO.EQ.11) THEN C ====================== C REST OF DECAY PROCEDURE FOR ACCEPTED TAU+ DECAY NEV1=NEV1+1 ISGN= IDF/IABS(IDF) CALL DEKAY1(1,H,ISGN) ELSEIF(KTO.EQ.12) THEN C ====================== C REST OF DECAY PROCEDURE FOR ACCEPTED TAU- DECAY NEV2=NEV2+1 ISGN=-IDF/IABS(IDF) CALL DEKAY2(1,H,ISGN) ELSEIF(KTO.EQ.100) THEN C ======================= IF(JAK1.NE.-1.OR.JAK2.NE.-1) THEN CALL DADMEL( 1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5) CALL DADMMU( 1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5) CALL DADMPI( 1,IDUM,PDUM,PDUM1,PDUM2) CALL DADMRO( 1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4) CALL DADMAA( 1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5,JDUM) CALL DADMKK( 1,IDUM,PDUM,PDUM1,PDUM2) CALL DADMKS( 1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,JDUM) CALL DAD4PI( 1,IDUM,HDUM,PDUM1,PDUM2,PDUMX,JDUM) CALL DADNPI( 1,IDUM, PDUM1,PDUM2,PDUMX,JDUM) CALL DADMPK( 1,IDUM,HDUM,PDUM1,PDUM2,PDUM3,PDUM4,PDUM5,JDUM) WRITE(IOUT,7010) NEV1,NEV2,NEVTOT WRITE(IOUT,7011) (NEVDEC(I),GAMPMC(I),GAMPER(I),I= 1,11) WRITE(IOUT,7012) (NEVDEC(I),GAMPMC(I),GAMPER(I),I=12,19) ENDIF ELSE C ==== GOTO 910 ENDIF C ===== DO 78 K=1,4 78 HX(K)=H(K) RETURN 7001 FORMAT(///1X,15(5H*****) $ /,' *', 25X,'*****TAUOLA LIBRARY: VERSION 2.4 ******',9X,1H*, $ /,' *', 25X,'***********NOVEMBER 1992***************',9X,1H*, $ /,' *', 25X,'**AUTHORS: S.JADACH, Z.WAS*************',9X,1H*, $ /,' *', 25X,'**R. DECKER, M. JEZABEK, J.H.KUEHN*****',9X,1H*, $ /,' *', 25X,'**AVAILABLE FROM: WASM AT CERNVM ******',9X,1H*, $ /,' *', 25X,'***** PUBLISHED IN COMP. PHYS. COMM.***',9X,1H*, $ /,' *', 25X,'*******CERN-TH-5856 SEPTEMBER 1990*****',9X,1H*, $ /,' *', 25X,'*******CERN-TH-6195 SEPTEMBER 1991*****',9X,1H*, $ /,' *', 25X,'*******CERN TH NOVEMBER 1992*****',9X,1H*, $ /,' *', 25X,'****DEKAY ROUTINE: INITIALIZATION******',9X,1H*, $ /,' *',I20 ,5X,'JAK1 = DECAY MODE TAU+ ',9X,1H*, $ /,' *',I20 ,5X,'JAK2 = DECAY MODE TAU- ',9X,1H*, $ /,1X,15(5H*****)/) 7010 FORMAT(///1X,15(5H*****) $ /,' *', 25X,'*****TAUOLA LIBRARY: VERSION 2.4 ******',9X,1H*, $ /,' *', 25X,'***********NOVEMBER 1992***************',9X,1H*, $ /,' *', 25X,'**AUTHORS: S.JADACH, Z.WAS*************',9X,1H*, $ /,' *', 25X,'**R. DECKER, M. JEZABEK, J.H.KUEHN*****',9X,1H*, $ /,' *', 25X,'**AVAILABLE FROM: WASM AT CERNVM ******',9X,1H*, $ /,' *', 25X,'***** PUBLISHED IN COMP. PHYS. COMM.***',9X,1H*, $ /,' *', 25X,'*******CERN-TH-5856 SEPTEMBER 1990*****',9X,1H*, $ /,' *', 25X,'*******CERN-TH-6195 SEPTEMBER 1991*****',9X,1H*, $ /,' *', 25X,'*******CERN TH NOVEMBER 1992*****',9X,1H*, $ /,' *', 25X,'*****DEKAY ROUTINE: FINAL REPORT*******',9X,1H*, $ /,' *',I20 ,5X,'NEV1 = NO. OF TAU+ DECS. ACCEPTED ',9X,1H*, $ /,' *',I20 ,5X,'NEV2 = NO. OF TAU- DECS. ACCEPTED ',9X,1H*, $ /,' *',I20 ,5X,'NEVTOT = SUM ',9X,1H*, $ /,' *',' NOEVTS ', $ ' PART.WIDTH ERROR ROUTINE DECAY MODE ',9X,1H*) 7011 FORMAT(1X,'*' $ ,I10,2F12.7 ,' DADMEL ELECTRON ',9X,1H* $ /,' *',I10,2F12.7 ,' DADMMU MUON ',9X,1H* $ /,' *',I10,2F12.7 ,' DADMPI PION ',9X,1H* $ /,' *',I10,2F12.7, ' DADMRO RHO (->2PI) ',9X,1H* $ /,' *',I10,2F12.7, ' DADMAA A1 (->3PI) ',9X,1H* $ /,' *',I10,2F12.7, ' DADMKK KAON ',9X,1H* $ /,' *',I10,2F12.7, ' DADMKS K* ',9X,1H* $ /,' *',I10,2F12.7, ' DAD4PI PI- 2PI+ PI0 ',9X,1H* $ /,' *',I10,2F12.7, ' DAD4PI PI- 3PI0 ',9X,1H* $ /,' *',I10,2F12.7, ' DADNPI 3PI- 2PI+ ',9X,1H* $ /,' *',I10,2F12.7, ' DADNPI 3PI- 2PI+ PI0 ',9X,1H*) 7012 FORMAT(1X,'*' $ ,I10,2F12.7, ' DADMPK K+ K- PI+ ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK K0B K0 PI+ ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK K+ K0B PI0 ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK K+ PI0 PI0 ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK K+ PI- PI+ ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK K0 PI0 PI+ ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK ET PI- PI0 ',9X,1H* $ /,' *',I10,2F12.7, ' DADMPK PI-PI0 GAM ',9X,1H* $ /,' *',20X,'THE ERROR IS RELATIVE AND PART.WIDTH ',10X,1H* $ /,' *',20X,'IN UNITS GFERMI**2*MASS**5/192/PI**3 ',10X,1H* $ /,1X,15(5H*****)/) 902 PRINT 9020 9020 FORMAT(' ----- DEKAY: LACK OF INITIALISATION') STOP 910 PRINT 9100 9100 FORMAT(' ----- DEKAY: WRONG VALUE OF KTO ') STOP END