RE: optical photons by scintillation/cherenkov process

From: ercan pilicer <ercanpilicer_at_gmail.com>
Date: Sat, 24 Jan 2009 19:01:00 +0200

hi again

thanks Markus for your comment. but i always get the same energy value
(2.75E-9 eV = 420 nm), not a distribution. that energy value is the
center of wavelength given by ;

* Set Light production/transport properties wavelength limits for PWO
crsytal
OPT-PROP 3.00D-05 4.20D-05 6.50D-05 PWO PWO WV-LIMIT

unfortunately there is still no fluctuation in produced optical photons!

ciao
e.

On Fri, 2009-01-23 at 14:39 +0100, Markus Brugger wrote:
> 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--
>

-- 
==================================================
()".'."()   Ercan Pilicer
(  (T)  )   Uludag University, Turkey
(??)   (??)   High Energy Physics Department
==================================================
Received on Sat Jan 24 2009 - 18:38:56 CET

This archive was generated by hypermail 2.2.0 : Sat Jan 24 2009 - 18:38:59 CET