* * $Id$ * * $Log$ * Revision 1.1 2000/06/19 20:00:31 eugenio * Initial revision * * Revision 1.1.1.1 1994/11/22 16:57:03 zfiles * first version of korb in CVS * * #include "sys/CLEO_machine.h" #include "pilot.h" *CMZ : 2.00/00 21/01/93 15.42.31 by Alan Weinstein *-- Author : SUBROUTINE DEXRO(MODE,ISGN,POL,PNU,PRO,PIC,PIZ) C ---------------------------------------------------------------------- C THIS SIMULATES TAU DECAY IN TAU REST FRAME C INTO NU RHO, NEXT RHO DECAYS INTO PION PAIR. C OUTPUT FOUR MOMENTA: PNU TAUNEUTRINO, C PRO RHO C PIC PION CHARGED C PIZ PION ZERO C ---------------------------------------------------------------------- COMMON / INOUT / INUT,IOUT REAL POL(4),HV(4),PRO(4),PNU(4),PIC(4),PIZ(4) DATA IWARM/0/ C IF(MODE.EQ.-1) THEN C =================== IWARM=1 CALL DADMRO( -1,ISGN,HV,PNU,PRO,PIC,PIZ) CC CALL HBOOK1(816,'WEIGHT DISTRIBUTION DEXRO $',100,0,2) CC CALL HBOOK1(916,'ABS2 OF HV IN ROUTINE DEXRO $',100,0,2) C ELSEIF(MODE.EQ. 0) THEN C ======================= 300 CONTINUE IF(IWARM.EQ.0) GOTO 902 CALL DADMRO( 0,ISGN,HV,PNU,PRO,PIC,PIZ) WT=(1+POL(1)*HV(1)+POL(2)*HV(2)+POL(3)*HV(3))/2. CC CALL HFILL(816,WT) CC XHELP=HV(1)**2+HV(2)**2+HV(3)**2 CC CALL HFILL(916,XHELP) CALL RANMAR(RN,1) IF(RN.GT.WT) GOTO 300 C ELSEIF(MODE.EQ. 1) THEN C ======================= CALL DADMRO( 1,ISGN,HV,PNU,PRO,PIC,PIZ) CC CALL HPRINT(816) CC CALL HPRINT(916) ENDIF C ===== RETURN 902 WRITE(IOUT, 9020) 9020 FORMAT(' ----- DEXRO: LACK OF INITIALISATION') STOP END