[fluka-discuss]: Optical Photon Example - Cherenkov

From: Justin Raudabaugh <justin.raudabaugh_at_duke.edu>
Date: Wed, 8 Feb 2017 17:55:47 +0000

Dear FLUKA experts,

I am trying to understand generation of optical photons, and attempting to recreate the code from the Cherenkov example in section 12.2.3 of the FLUKA manual. I have typed up the code, but I am not seeing any photons generated. The code runs (suspiciously quickly for a large number of histories) , but doesn’t appear to do anything. Any thoughts on what the problem could be? Code attached


TITLE
scint
* Set the defaults for precision simulations
DEFAULTS PRECISIO
* Define the beam characteristics
BEAM -10.00 MUON+
* Define the beam position
BEAMPOS 0.0 0.0 199.0 NEGATIVE
DELTARAY -1.0 ARGON ARGON NOPRINT
PAIRBREM -3. ARGON ARGON
MUPHOTON -1.0 ARGON ARGON
PHOTONUC -1. HYDROGEN NICKEL
IONTRANS 4-HELIUM
DISCARD NEUTRIM ANEUTRIM NEUTRIT ANEUTRIT NEUTRIE ANEUTRIE
GEOBEGIN COMBNAME
    0 0
RPP blkhole -999.0 999.0 -999.0 999.0 -999.0 999.0
RPP argon -200.0 200.0 -200.0 200.0 -200.0 200.0
END
* Black hole
BLKBODY 5 +blkhole -argon
* Target
liquid 5 +argon
END
GEOEND
EMF EMF-OFF
MATERIAL 18.0 39.948 0.00166 20. ARGON
LOW-MAT ARGON 18. -2. 87. ARGON
OPT-PROP 1.000E-05 3.5E-05 6.000E-05 HYDROGEN NICKEL WV-LIMIT
OPT-PROP 1.0 HYDROGEN NICKEL METAL
OPT-PROP ARGON ARGON RESET
OPT-PROD ARGON ARGON SCIN-OFF
OPT-PROD 2.0E-09 3.0E-09 ARGON ARGON CERENKOV
MAT-PROP 1. ARGON ARGON USERDIRE
OPT-PROD 1.100E-05 5.500E-05 ARGON ARGON CEREN-WV
OPT-PROP 1.000E-05 3.500E-05 6.000E-05 ARGON ARGON WV-LIMIT
OPT-PROP -100.0 0.001 0.01111 ARGON ARGON
OPT-PROP 0.1 0.1 SENSITIV
SCORE ENERGY OPTIPHOT ALL-PART BEAMPART
* Set the random number seed
RANDOMIZ 1.0
USRBDX 1. OPTIPHOT -55. liquid BLKBODY 1.0opt.phot
USRBDX 12.0E-09 0.0 120.0 &
USRBIN 10. ENERGY -56. 200. 200. 200.optphotbin
USRBIN -200. -200. -200. 200. 200. 200. &
*USERDUMP 100. 49. 2. MGDRAW
* Set the number of primary histories to be simulated in the run
START 10000.0
STOP

Thanks for your time,
Justin Raudabaugh��������������������������������������������������������(�Ƨ��ځ맲��r���*'~�&������x%��Z�m����0��n���ߖ���i�'i���w�
Received on Wed Feb 08 2017 - 20:24:16 CET