RE: optical photons by scintillation/cherenkov process

Hi Ercan,

the conversion from energy to wavelength makes sense, as actually does
the result. Contrary to Cherenkov, scintillation light in FLUKA is only
generated for monochromatic photons (maximum of 3 different lines
possible by repeating OPT-PROD with SCINTILL, c.f., manual 12.2.3). For
typical problems this is sufficient as one is interested in few
(sometimes only one, e.g., 128nm for liquid Ar) characteristic lines.
One could possibly think of a dedicated user source routine, however
when checking yesterday we couldn't find a readily prepared interface.

cheers

Markus

> -----Original Message-----
> From: ercan pilicer [mailto:ercanpilicer_at_gmail.com]
> Sent: Montag, 26. J=E4nner 2009 08:55
> To: Markus Brugger
> Cc: fluka-discuss_at_fluka.org
> Subject: RE: optical photons by scintillation/cherenkov process
>=20
>=20
> here i attached my input and mgdraw.f routine
>=20
>=20
> On Sun, 2009-01-25 at 20:35 +0100, Markus Brugger wrote:
> > Hi Ecran,
> >
> > 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
> > >
> > > hi again
> > >
> > > 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 ;
> > >
> > > * 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 =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=

> > >
> > >
> >
> --
> =
=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
> (=F6) (=F6) 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=
Received on Tue Jan 27 2009 - 13:35:45 CET