RE: optical photons by scintillation/cherenkov process

From: Markus Brugger <Markus.Brugger_at_cern.ch>
Date: Fri, 23 Jan 2009 14:39:00 +0100

Hi Ecran,

I think I understand why you always get 0. for your <wvlngt> parameters.
In your print statement you use the following variables out of the
OPPHCM common, the latter actually being the 'OPtical PHoton CoMmon: '
information (thus constant). Out of the four variables you print, three
actually are in this common:

  Wvmnsn = minimum wavelength for opt. photon sensitivity
  Wvcnsn = central wavelength for opt. photon sensitivity
  Wvmxsn = maximum wavelength for opt. photon sensitivity

What concerns the actual wavelength I a bit puzzled about the variable
name you're using <wvlngt>, which is not included in OOPHCM nor in the
respective stack (OOPHST). Possibly not being initialised nor defined
anywhere else in your mgdraw routine, it's '0' by default.

Concerning the wave length of the optical photon you could theoretically
get it from the OOPHST (OPtical PHoton STack), however as a general rule
it's strongly recommended to use the TRACKR common (as the particle
energy might change during a step, even though not the case for optical
photons).

Thus, please try using TRACKR common, checking for optical photons as
you do and then converting the energy (GeV) into wavelength.

hope this helps, cheers

Markus

> -----Original Message-----
> From: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-
> discuss_at_mi.infn.it] On Behalf Of ercan pilicer
> Sent: Donnerstag, 22. Jšnner 2009 10:04
> To: fluka-discuss_at_fluka.org
> Subject: optical photons by scintillation/cherenkov process
>
>
> --=3D-fyAaG3HVQvBjw1SihXoi
> Content-Type: text/plain; charset=3DUTF-8
> Content-Transfer-Encoding: 8bit
>
> dear users
>
> i want to get a distribution for optical photons generated by
> scintillation or cherenkov process inside and at the exit of a pwo
> crystal. before starting using usrbdx and usryield cards i wanted to
> check wavelengths of generated optical photons. i always get 0. values
> for wvlngt parameter like;
>
> =2E..
> 0. 2.5E-06 5.89E-05 0.0006
> 0. 2.5E-06 5.89E-05 0.0006
> 0. 2.5E-06 5.89E-05 0.0006
> 0. 2.5E-06 5.89E-05 0.0006
> 0. 2.5E-06 5.89E-05 0.0006
> =2E..
>
> part of modified mgdraw.f routine is ;
>
> ...
> INCLUDE '(OPPHCM)'
> INCLUDE '(OPPHST)'
> ...
> if(jtrack.eq.-1)then
> write(81,*) wvlngt,wvmnsn,wvcnsn,wvmxsn
> c if(Lcrnkv(27).eqv..true.) mykcrn=3Dmykcrn+1
> c if(Lscntl(27).eqv..true.) mykscn=3Dmykscn+1
> c write(82,*) mykcrn,mykscn
> endif
> ...
> ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO )
> if(jtrack.eq.-1)then
> write(88,*) wvlngt,wvmnsn,wvcnsn,wvmxsn
> c if(Lcrnkv(27).eqv..true.) mykbdxcrn=3Dmykbdxcrn+1
> c if(Lscntl(27).eqv..true.) mykbdxscn=3Dmykbdxscn+1
> c write(82,*) mykcrn,mykscn
> endif
> RETURN
> ...
>
> i attach my inp file as well.
>
> any comments would be deeply appreciated
>
> ciao
> e.
>
>
> --=20
> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3
> D=3D=
> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3
> D=3D
> ()".'."() Ercan Pilicer
> ( (T) ) Uludag University, Turkey
> (=C3=B6) (=C3=B6) High Energy Physics Department
> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3
> D=3D=
> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3
> D=3D
>
>
> --=3D-fyAaG3HVQvBjw1SihXoi
> Content-Disposition: attachment; filename=3Dpwo.inp
> Content-Type: text/x-fluka; name=3Dpwo.inp; charset=3DUTF-8
> Content-Transfer-Encoding: 7bit
>
> TITLE
> Test of scintillation light production in PbWO4 crystal
> DEFAULTS
> PRECI=
> SIO
> BEAM -0.01
> ELECT=
> RON
> BEAMPOS 0.0 0.0 -11.5
> DISCARD NEUTRIM ANEUTRIM NEUTRIT ANEUTRIT NEUTRIE ANEUTRIE
> GEOBEGIN
> COMBN=
> AME
> 0 0 Single PWO crystal (3x3x23cm) with Silicon (5mm)
> RPP B1 -9999999. +9999999. -9999999. +9999999. -9999999.
> +9999999.
> RPP B2 -20. 20. -20. 20. -20. 20.
> RPP B3 -1.6 1.6 -1.6 1.6 -11.5 11.5
> RPP B4 -1.5 1.5 -1.5 1.5 -11.5 11.5
> RPP B5 -1.5 1.5 -1.5 1.5 11.5 12.
> END
> blkHole 5 +B1 -B2
> expHole 5 +B2 -B3 -B5
> Wrap 5 +B3 -B4
> Crystal 5 +B4
> SiDet 5 +B5
> END
> GEOEND
> MATERIAL 0.0012048 26.0
> AIR
> MATERIAL 8.28 27.0
> PWO
> MATERIAL 0.331 28.0
> TYVEK
> COMPOUND -0.000124 CARBON -0.755267 NITROGEN -0.231781
> OXYGENAIR
> COMPOUND -0.012827 ARGON 0.0 0.0
> AIR
> COMPOUND 1.0 LEAD 1.0 TUNGSTEN 4.0
> OXYGENPWO
> COMPOUND 2.0 HYDROGEN 1.0 OXYGEN
> TYVEK
> ASSIGNMA BLCKHOLE blkHole 0.0
> ASSIGNMA AIR expHole
> ASSIGNMA TYVEK Wrap
> ASSIGNMA PWO Crystal
> ASSIGNMA SILICON SiDet
> * switches on Scintillation light production in PWO crystal
> OPT-PROD PWO
> SCINT=
> ILL
> #if 0
> OPT-PROD PWO
> CERE-=
> OFF
> #endif
> * Set Light production/transport properties: from 100 to 600 nm in PWO
> crs=
> ytal
> OPT-PROP 3.00D-05 4.20D-05 6.50D-05 HYDROGEN @LASTMAT
> WV-LI=
> MIT
> * The following card restores the wave-length limits for PWO crsytal
> OPT-PROP 3.00D-05 4.20D-05 6.50D-05 PWO
> WV-LI=
> MIT
> * Set wrapping material to "metal" with 0 reflectivity:
> OPT-PROP 1.0 TYVEK
> METAL
> #if 0
> OPT-PROP PWO
> RESET
> #endif
> * defines Scint. light production for PWO crystal parameters:
> * a) wavelength (cm) of first scintillation line.
> * b) fraction of deposited energy going into scint. light
> * (in PWO ~ 150 photons/MeV)
> OPT-PROD 4.50D-05 2.755D-2 PWO
> SCINT=
> -WV
> * The following card, for PWO crystal
> * a) sets to 2.29 for refraction index
> * b) sets to 1000 cm the mean free path for absorption.
> * c) sets to 90 cm the mean free path for Rayleigh scattering
> OPT-PROP 2.29 0.001 0.01111 PWO
> * The following card defines the "Sensitivity" in order to introduce
> the
> * maximum Quantum Efficiency at generation level. Here 1/10 of photons
> are
> * actually generated.
> * Fluctuations are properly sampled
> OPT-PROP 0.1
> SENSI=
> TIV
> SCORE ENERGY OPTIPHOT
> #if 0
> USRBDX 101. OPTIPHOT -55. Crystal SiDet
> OptPh=
> ot1
> USRBDX 10. 1D-12 200. &
> #endif
> #if 0
> USRBDX 99. OPTIPHOT -55. Crystal SiDet
> OptPh=
> ot2
> USRBDX 10. 1D-12 200. &
> #endif
> #if 0
> USRBDX 1. OPTIPHOT -56. Crystal SiDet
> OptPh=
> ot1
> USRBDX 200. &
> #endif
> #if 0
> USRBDX -1. OPTIPHOT -56. Crystal SiDet
> OptPh=
> ot2
> USRBDX 200. &
> #endif
> #if 0
> USRBIN 10. OPTIPHOT -57. 1.5 1.5
> 11.5OptPh=
> o
> USRBIN -1.5 -1.5 -11.5 100. 100. 100.&
> #endif
> #if 0
> USRBIN 10. ENERGY -57. 1.5 1.5
> 11.5Edep
> USRBIN -1.5 -1.5 -11.5 100. 100. 100.&
> #endif
> #if 0
> USRBIN 10. ENERGY -57. 20. 20.
> 20.EdepA=
> ll
> USRBIN -20. -20. -20. 100. 100. 100.&
> #endif
> #if 0
> USRYIELD 2401. OPTIPHOT -58. -1. -2.
> Scint=
> 1
> USRYIELD 4.13D-9 1.55D-9 100. 6.2831853 0.0 2703.&
> #endif
> #if 0
> USRYIELD 2501. OPTIPHOT -58. -1. -2.
> Scint=
> 2
> USRYIELD 4.13D-9 1.55D-9 100. 6.2831853 0.0 2703.&
> #endif
> #if 0
> USRTRACK 1. OPTIPHOT -59. Crystal
> 100.Scint=
> 1
> USRTRACK 4.13D-9 1.55D-9 &
> #endif
> #if 0
> USRTRACK -1. OPTIPHOT -59. Crystal
> 100.Scint=
> 2
> USRTRACK 4.13D-9 1.55D-9 &
> #endif
> USERDUMP 100. 21. 2.
> MGDRA=
> W
> RANDOMIZ 1.0
> START 10.
> STOP
>
> --=3D-fyAaG3HVQvBjw1SihXoi--
>
> --3D--
Received on Fri Jan 23 2009 - 15:11:33 CET

This archive was generated by hypermail 2.2.0 : Fri Jan 23 2009 - 15:11:37 CET