Dear FLUKA users and experts,
I want to simulate the wave length shifted photon production in selected
region but I'm not good in FORTRAN. Can I calculate a number of absorbed
photons in region also? I don't really understand how WVLNSH subroutine is
working. I ask you help in this thing.
Best regards,
Viktor
*$ CREATE WVLNSH.FOR
*COPY WVLNSH
*
*=== Wvlnsh ===========================================================*
*
SUBROUTINE WVLNSH ( EKPHOT, WVLNGT, MREG , MMAT , NWVSHP,
& WVSHPH, DWVSHP, MXWSPH )
INCLUDE '(DBLPRC)'
INCLUDE '(DIMPAR)'
INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1997-2010 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* WaVe LeNgth SHifted photon production: *
* *
* Created on 21 september 1997 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 16-Dec-10 by Alfredo Ferrari *
* *
* Note the wavelength is always assumed to be the vacuum one *
* *
* Input variables: *
* *
* Ekphot = absorbed photon energy (GeV) *
* Wvlngt = absorbed photon wave-length (cm) *
* Mreg = current region *
* Mmat = current material *
* Mxwsph = maximum number of secondary photons *
* *
* *
* Output variables: *
* *
* Nwvshp = number of produced secondary photons *
* Wvshph(i) > 0 -> energy (GeV) of the i_th produced photon *
* < 0 -> wavelength (cm) of the i_th produced photon *
* Dwvshp(i) = production delay (s) of the i_th produced photon *
* *
*----------------------------------------------------------------------*
*
INCLUDE '(FLKMAT)'
DIMENSION WVSHPH (MXWSPH), DWVSHP (MXWSPH)
* IF ( MMAT .NE. 20 ) THEN
CALL FLABRT ( 'WVLNSH', 'MMAT IS NOT SHIFTER!' )
END IF
NWVSHP = 1
Wvshph = -4.2D-05
RETURN
*=== End of subroutine Wvlnsh =========================================*
END
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sun Dec 18 2016 - 20:02:15 CET