* * $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.29 by Alan Weinstein *-- Author : SUBROUTINE DEXEL(MODE,ISGN,POL,PNU,PWB,Q1,Q2,PH) C ---------------------------------------------------------------------- C THIS SIMULATES TAU DECAY IN TAU REST FRAME C INTO ELECTRON AND TWO NEUTRINOS C C called by : DEXAY,DEXAY1 C ---------------------------------------------------------------------- REAL POL(4),HV(4),PWB(4),PNU(4),Q1(4),Q2(4),PH(4) DATA IWARM/0/ C IF(MODE.EQ.-1) THEN C =================== IWARM=1 CALL DADMEL( -1,ISGN,HV,PNU,PWB,Q1,Q2,PH) CC CALL HBOOK1(813,'WEIGHT DISTRIBUTION DEXEL $',100,0,2) C ELSEIF(MODE.EQ. 0) THEN C ======================= 300 CONTINUE IF(IWARM.EQ.0) GOTO 902 CALL DADMEL( 0,ISGN,HV,PNU,PWB,Q1,Q2,PH) WT=(1+POL(1)*HV(1)+POL(2)*HV(2)+POL(3)*HV(3))/2. CC CALL HFILL(813,WT) CALL RANMAR(RN,1) IF(RN.GT.WT) GOTO 300 C ELSEIF(MODE.EQ. 1) THEN C ======================= CALL DADMEL( 1,ISGN,HV,PNU,PWB,Q1,Q2,PH) CC CALL HPRINT(813) ENDIF C ===== RETURN 902 PRINT 9020 9020 FORMAT(' ----- DEXEL: LACK OF INITIALISATION') STOP END