RE: optical photons by scintillation/cherenkov process

From: Markus Brugger <Markus.Brugger_at_cern.ch>
Date: Sun, 25 Jan 2009 20:35:23 +0100

Hi Ercan,

did you use TRACKR to retrieve the energy (wavelength) of the photon? If
yes, could you please attach (please not inline) your input-file as well
as the respective user routine.

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: 24 January 2009 18:01
> To: Markus Brugger
> Cc: fluka-discuss_at_fluka.org
> Subject: RE: optical photons by scintillation/cherenkov process
>=20
> hi again
>=20
> thanks Markus for your comment. but i always get the same energy
> value
> (2.75E-9 eV =3D 420 nm), not a distribution. that energy value is the
> center of wavelength given by ;
>=20
> * 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
>=20
> unfortunately there is still no fluctuation in produced optical
> photons!
>=20
> ciao
> e.
>=20
>=20
>=20
>=20
> 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 =3D minimum wavelength for opt. photon sensitivity
> > Wvcnsn =3D central wavelength for opt. photon sensitivity
> > Wvmxsn =3D 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
> > >
> > >
> > > --=3D3D-fyAaG3HVQvBjw1SihXoi
> > > Content-Type: text/plain; charset=3D3DUTF-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;
> > >
> > > =3D2E..
> > > 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
> > > =3D2E..
> > >
> > > 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=3D3Dmykcrn+1
> > > c if(Lscntl(27).eqv..true.) mykscn=3D3Dmykscn+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=3D3Dmykbdxcrn+1
> > > c if(Lscntl(27).eqv..true.) mykbdxscn=3D3Dmykbdxscn+1
> > > c write(82,*) mykcrn,mykscn
> > > endif
> > > RETURN
> > > ...
> > >
> > > i attach my inp file as well.
> > >
> > > any comments would be deeply appreciated
> > >
> > > ciao
> > > e.
> > >
> > >
> > > --=3D20
> > >
> =
=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D=
3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D
> =3D3
> > > D=3D3D=3D
> > >
> =
=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D=
3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D
> =3D3
> > > D=3D3D
> > > ()".'."() Ercan Pilicer
> > > ( (T) ) Uludag University, Turkey
> > > (=3DC3=3DB6) (=3DC3=3DB6) High Energy Physics Department
> > >
> =
=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D=
3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D
> =3D3
> > > D=3D3D=3D
> > >
> =
=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D=
3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D=3D3D
> =3D3
> > > D=3D3D
> > >
> > >
> > > --=3D3D-fyAaG3HVQvBjw1SihXoi
> > > Content-Disposition: attachment; filename=3D3Dpwo.inp
> > > Content-Type: text/x-fluka; name=3D3Dpwo.inp; charset=3D3DUTF-8
> > > Content-Transfer-Encoding: 7bit
> > >
> > > TITLE
> > > Test of scintillation light production in PbWO4 crystal
> > > DEFAULTS
> > > PRECI=3D
> > > SIO
> > > BEAM -0.01
> > > ELECT=3D
> > > RON
> > > BEAMPOS 0.0 0.0 -11.5
> > > DISCARD NEUTRIM ANEUTRIM NEUTRIT ANEUTRIT NEUTRIE
> ANEUTRIE
> > > GEOBEGIN
> > > COMBN=3D
> > > 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=3D
> > > ILL
> > > #if 0
> > > OPT-PROD PWO
> > > CERE-=3D
> > > OFF
> > > #endif
> > > * Set Light production/transport properties: from 100 to 600 nm
> in PWO
> > > crs=3D
> > > ytal
> > > OPT-PROP 3.00D-05 4.20D-05 6.50D-05 HYDROGEN @LASTMAT
> > > WV-LI=3D
> > > 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=3D
> > > 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=3D
> > > -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=3D
> > > TIV
> > > SCORE ENERGY OPTIPHOT
> > > #if 0
> > > USRBDX 101. OPTIPHOT -55. Crystal SiDet
> > > OptPh=3D
> > > ot1
> > > USRBDX 10. 1D-12 200.
> &
> > > #endif
> > > #if 0
> > > USRBDX 99. OPTIPHOT -55. Crystal SiDet
> > > OptPh=3D
> > > ot2
> > > USRBDX 10. 1D-12 200.
> &
> > > #endif
> > > #if 0
> > > USRBDX 1. OPTIPHOT -56. Crystal SiDet
> > > OptPh=3D
> > > ot1
> > > USRBDX 200.
> &
> > > #endif
> > > #if 0
> > > USRBDX -1. OPTIPHOT -56. Crystal SiDet
> > > OptPh=3D
> > > ot2
> > > USRBDX 200.
> &
> > > #endif
> > > #if 0
> > > USRBIN 10. OPTIPHOT -57. 1.5 1.5
> > > 11.5OptPh=3D
> > > 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=3D
> > > ll
> > > USRBIN -20. -20. -20. 100. 100.
> 100.&
> > > #endif
> > > #if 0
> > > USRYIELD 2401. OPTIPHOT -58. -1. -2.
> > > Scint=3D
> > > 1
> > > USRYIELD 4.13D-9 1.55D-9 100. 6.2831853 0.0
> 2703.&
> > > #endif
> > > #if 0
> > > USRYIELD 2501. OPTIPHOT -58. -1. -2.
> > > Scint=3D
> > > 2
> > > USRYIELD 4.13D-9 1.55D-9 100. 6.2831853 0.0
> 2703.&
> > > #endif
> > > #if 0
> > > USRTRACK 1. OPTIPHOT -59. Crystal
> > > 100.Scint=3D
> > > 1
> > > USRTRACK 4.13D-9 1.55D-9
> &
> > > #endif
> > > #if 0
> > > USRTRACK -1. OPTIPHOT -59. Crystal
> > > 100.Scint=3D
> > > 2
> > > USRTRACK 4.13D-9 1.55D-9
> &
> > > #endif
> > > USERDUMP 100. 21. 2.
> > > MGDRA=3D
> > > W
> > > RANDOMIZ 1.0
> > > START 10.
> > > STOP
> > >
> > > --=3D3D-fyAaG3HVQvBjw1SihXoi--
> > >
> > > --3D--
> >
> --
> =
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=

> ()".'."() Ercan Pilicer
> ( (T) ) Uludag University, Turkey
> (??) (??) 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=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=

>=20
>=20
Received on Sun Jan 25 2009 - 21:01:43 CET

This archive was generated by hypermail 2.2.0 : Sun Jan 25 2009 - 21:01:43 CET