|
|
|
The user can request any kind of standard FLUKA output for optical photons and also a user specific output, starting for instance from the MGDRAW user routine. Here an example follows, where a few variables are simply recorded in the output "collision tape" (dump file) at each step in the tracking only for particle id = -1 (optical photons). It can be useful, in order to exploit the flags available in this routine, and explained in the different comments, to know that the FLUKA routine which drives the transport of optical photons is KASOPH. The content of the TRACKR common can be used to fully exploit the possibilities offered from the use of MGDRAW routine (see Chap. 13}). Warning: in the present version of FLUKA, there is not yet the possibility of using the User Particle Properties for optical photons (variables SPAUSR, ISPUSR and the STUPFR user routine)* **=== mgdraw ===========================================================** *SUBROUTINE MGDRAW ( ICODE, MREG ) INCLUDE '(DBLPRC)' INCLUDE '(DIMPAR)' INCLUDE '(IOUNIT)'**----------------------------------------------------------------------** ** 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 '(MGDDCM)' INCLUDE '(PAPROP)' INCLUDE '(QUEMGD)' INCLUDE '(SUMCOU)' INCLUDE '(TRACKR)'*DIMENSION DTQUEN ( MXTRCK, MAXQMG )*CHARACTER*20 FILNAM LOGICAL LFCOPE SAVE LFCOPE DATA LFCOPE / .FALSE. /**----------------------------------------------------------------------** ** Icode = 1: call from Kaskad ** Icode = 2: call from Emfsco ** Icode = 3: call from Kasneu ** Icode = 4: call from Kashea ** Icode = 5: call from Kasoph ** **----------------------------------------------------------------------** *IF ( .NOT. LFCOPE ) THEN LFCOPE = .TRUE. IF ( KOMPUT .EQ. 2 ) THEN FILNAM = '/'//CFDRAW(1:8)//' DUMP A' ELSE FILNAM = CFDRAW END IF WRITE(*,*) 'TRAJECTORY OPEN!' WRITE(*,'(A)') 'FILNAM = ',FILNAM OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM = & 'UNFORMATTED' ) END IF C C Write trajectories of optical photons C IF(JTRACK .EQ. -1) THEN WRITE (IODRAW) NTRACK, MTRACK, JTRACK, SNGL (ETRACK), & SNGL (WTRACK) WRITE (IODRAW) ( SNGL (XTRACK (I)), SNGL (YTRACK (I)), & SNGL (ZTRACK (I)), I = 0, NTRACK ), & ( SNGL (DTRACK (I)), I = 1,MTRACK ), & SNGL (CTRACK) WRITE(IODRAW) SNGL(CXTRCK),SNGL(CYTRCK),SNGL(CZTRCK) ENDIF RETURN**======================================================================** ** Boundary-(X)crossing DRAWing: ** ** Icode = 1x: call from Kaskad ** 19: boundary crossing ** Icode = 2x: call from Emfsco ** 29: boundary crossing ** Icode = 3x: call from Kasneu ** 39: boundary crossing ** Icode = 4x: call from Kashea ** 49: boundary crossing ** Icode = 5x: call from Kasoph ** 59: boundary crossing ** **======================================================================** *ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO ) RETURN**======================================================================** ** Event End DRAWing: ** **======================================================================** *ENTRY EEDRAW ( ICODE ) RETURN**======================================================================** ** ENergy deposition DRAWing: ** ** Icode = 1x: call from Kaskad ** 10: elastic interaction recoil ** 11: inelastic interaction recoil ** 12: stopping particle ** 13: pseudo-neutron deposition ** 14: escape ** 15: time kill ** Icode = 2x: call from Emfsco ** 20: local energy deposition (i.e. photoelectric) ** 21: below threshold, iarg=1 ** 22: below threshold, iarg=2 ** 23: escape ** 24: time kill ** Icode = 3x: call from Kasneu ** 30: target recoil ** 31: below threshold ** 32: escape ** 33: time kill ** Icode = 4x: call from Kashea ** 40: escape ** 41: time kill ** 42: delta ray stack overflow ** Icode = 5x: call from Kasoph ** 50: optical photon absorption ** 51: escape ** 52: time kill ** **======================================================================** *ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO ) RETURN**======================================================================** ** SOurce particle DRAWing: ** **======================================================================**ENTRY SODRAW* |* +-------------------------------------------------------------------*RETURN**======================================================================** ** USer dependent DRAWing: ** ** Icode = 10x: call from Kaskad ** 100: elastic interaction secondaries ** 101: inelastic interaction secondaries ** 102: particle decay secondaries ** 103: delta ray generation secondaries ** 104: pair production secondaries ** 105: bremsstrahlung secondaries ** 110: decay products ** Icode = 20x: call from Emfsco ** 208: bremsstrahlung secondaries ** 210: Moller secondaries ** 212: Bhabha secondaries ** 214: in-flight annihilation secondaries ** 215: annihilation at rest secondaries ** 217: pair production secondaries ** 219: Compton scattering secondaries ** 221: photoelectric secondaries ** 225: Rayleigh scattering secondaries ** Icode = 30x: call from Kasneu ** 300: interaction secondaries ** Icode = 40x: call from Kashea ** 400: delta ray generation secondaries ** For all interactions secondaries are put on GENSTK common (kp=1,np) ** but for KASHEA delta ray generation where only the secondary elec- ** tron is present and stacked on FLKSTK common for kp=npflka ** **======================================================================**ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )* No output by default:RETURN*=== End of subroutine Mgdraw ==========================================*END