*$ CREATE MGDRAW.FOR *COPY MGDRAW * * *=== mgdraw ===========================================================* * * SUBROUTINE MGDRAW ( ICODE, MREG ) INCLUDE '(DBLPRC)' INCLUDE '(DIMPAR)' INCLUDE '(IOUNIT)' * *----------------------------------------------------------------------* * * * Copyright (C) 1990-2006 by Alfredo Ferrari * * All Rights Reserved. * * * * * * MaGnetic field trajectory DRAWing: actually this entry manages * * all trajectory dumping for * * drawing * * * * Created on 01 march 1990 by Alfredo Ferrari * * INFN - Milan * * Last change 05-may-06 by Alfredo Ferrari * * INFN - Milan * * * *----------------------------------------------------------------------* * INCLUDE '(CASLIM)' INCLUDE '(COMPUT)' INCLUDE '(SOURCM)' INCLUDE '(FHEAVY)' INCLUDE '(FLKSTK)' INCLUDE '(GENSTK)' INCLUDE '(EMFSTK)' INCLUDE '(MGDDCM)' INCLUDE '(PAPROP)' INCLUDE '(QUEMGD)' INCLUDE '(SUMCOU)' INCLUDE '(TRACKR)' * DIMENSION DTQUEN ( MXTRCK, MAXQMG ) * CHARACTER*20 FILNAM CHARACTER*8 MRGNAM, NRGNAM LOGICAL LFCOPE SAVE LFCOPE DATA LFCOPE / .FALSE. / * LOGICAL UDFIRST DATA UDFIRST / .TRUE. / SAVE UDFIRST * RETURN * *======================================================================* * * ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO ) IF ( .NOT. LFCOPE ) THEN LFCOPE = .TRUE. OPEN ( UNIT = 80, FILE = "collfile1", STATUS = 'UNKNOWN') END IF CALL GEOR2N ( MREG, MRGNAM, IERR1 ) CALL GEOR2N ( NEWREG, NRGNAM, IERR2 ) IF(IERR1 .NE. 0 .OR. IERR2 .NE. 0) STOP "Error in name conversion" IF(MRGNAM .EQ. "TARG1" .AND. nrgnam .EQ. "TARG2" .AND. & JTRACK .EQ. -6) THEN WRITE(80,'(1P,10G25.15)') JTRACK, ZSCO,XSCO,YCSO, & ETRACK-AM(JTRACK), WTRACK END IF RETURN * *======================================================================* * * ENTRY EEDRAW ( ICODE ) RETURN * * *======================================================================* * * ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO ) RETURN * *======================================================================* * * ENTRY SODRAW RETURN * *======================================================================* * * ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO ) IF ( UDFIRST ) THEN UDFIRST = .FALSE. OPEN ( UNIT = 70, FILE = "filejuju_usd", STATUS = 'UNKNOWN') ENDIF do 20 ip=1,NPHEAV WRITE(80,'(1P,12G25.15)') 'cHeav', NCASE, JTRACK, & ip,ICODE,KHEAVY(ip),ICHEAV(KHEAVY(ip)),IBHEAV(KHEAVY(ip)), & TKHEAV(ip),CXR(ip),CYR(ip),CZR(ip) 20 continue RETURN *=== End of subrutine Mgdraw ==========================================* END