Hallo Viktor,
as indicated in the routine header you pasted, you can set the number of
*produced* photons as a function of the absorbing/emitting material and of
the energy (wave-length) of the absorbed photon.
If you set NWVSHP = 1, you need to fill WVSHPH(1)= and DWVSHP(1)=;
if the secondary photons are two (NWVSHP = 2), you need to define
WVSHPH(1), DWVSHP(1), WVSHPH(2), and DWVSHP(2).
Then, if you wish, you can also keep track of the absorption event by
writing something.
Kind regards
Francesco
**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854
On Sun, 18 Dec 2016, ๗ฆหิฯา ๒ฯฤฆฮ wrote:
> 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 Mon Dec 19 2016 - 12:15:58 CET