Re: Optical Photon Production by Scintillator Material, need help from Vahan Petrosyan on 2012-12-17 (fluka discuss archive)

From: Vahan Petrosyan <vpetrosyan_at_asls.candle.am>
Date: Mon, 17 Dec 2012 10:29:43 +0400

<201212121127.qBCBRKZp032216_at_smtp1.mi.infn.it> <CAKT5vuhfYOaz4Ai28MG7k=9Ruv2v2-odim5Y=tmpKjszLnz9zw_at_mail.gmail.com>
In-Reply-To: <CAKT5vuhfYOaz4Ai28MG7k=9Ruv2v2-odim5Y=tmpKjszLnz9zw_at_mail.gmail.com>
Content-Type: multipart/alternative;
boundary="------------080207000203040206030803"
X-Antivirus: avast! (VPS 121216-1, 16.12.2012), Outbound message
X-Antivirus-Status: Clean
X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.2.7 (smtp1.mi.infn.it [192.84.138.69]); Mon, 17 Dec 2012
07:29:58 +0100 (CET)
X-Virus-Scanned: clamav-milter 0.97.6 at smtp1.mi.infn.it
X-Virus-Status: Clean
Sender: owner-fluka-discuss_at_mi.infn.it

This is a multi-part message in MIME format.
--------------080207000203040206030803
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

Hi Nick. Thanks for your replay. It helped me much.

One question remains obscure for me.
1) How can one calculate WHAT(2) "fraction of deposited energy going
into i-th scin. photon emission" parameter?
2) Should it be specified in scintillator material parameters or it
can be calculated indirectly from other parameters too, for example
luminosity and emission peak wavelength?

I supposed formula Luminosity[ph/ev]*Energy(peak emission)[ev], am I right?

(I converted luminosity from [ph/Mev] to [ph/ev])

On 12.12.2012 19:58, Nicholas Bolibruch wrote:
> Hi Vahan, sorry I should have been more clear.
>
> >5) You need to set an appropriate upper and lower energy in the
> USRBDX card >for binning the optical wavelengths. I usually use 1E-9
> to 1E-8 initially.
>
> The USRBDX card for OPTIPHO needs to have the binning range set to
> appropriate energies.
>
>
> >6) You may want to set an OPT-PROP (sdum blank) card to set an index
> of >refraction, absorption and diffusion coefficients. This will
> increase the >processing time per particle however.
>
> This should read "You want to set an OPT-PROP ...". From this card you
> need to set _at least_ the absorption coefficient and the material
> you're assigning the properties to (ie. YAG). Setting diffusion
> coefficients will increase the time per particle.
>
> Attached is your last INP file modified to work. I include comments
> on cards I've changed.
>
> --Nick
>
>
> On Wed, Dec 12, 2012 at 5:23 AM, Vahan Petrosyan
> <vpetrosyan_at_asls.candle.am <mailto:vpetrosyan_at_asls.candle.am>> wrote:
>
>
> Many thanks for your help Nick!
>
> I did corrections that you mention and was able to get nice USRBDX
> result for photons but not for optical ones. It still zeroes.....(:
> I am attaching new corrected input file.
>
> As I understood from your explanation card OPT-PROP (sdum
> WV-LIMITS) is
> optional and card OPT-PROP (sdum SENSITIV ) is used only for reducing
> CPU load not for getting more precise results....?
> I switched on also OPT-PROD (sdum SCINT-WV and WV-LIMITS) but
> without
> luck!
> I am little confused about WHAT(2) ( "fraction") parameter of
> OPT-PROD(
> sdum SCINT-WV ), is it the same as /Luminosity [photons/Mev]???? /
>
> How can I get USRBDX distribution of optical photons?
>
> With regards V. Petrosyan!
>
> On 10.12.2012 21:44, Nicholas Bolibruch wrote:
> > Hi Vahan,
> >
> > Here are the issues with your input:
> >
> > 1) your source is at 20 MeV, and is isotropic. If you test it
> with a
> > beam, make sure its in the positive direction, and set the primaries
> > lower as the processing time is ~2 hrs on an i7 if it's made
> > directional with your current settings. Using the SENSITIV sdum on
> > OPT-PROP you can probably lower the OPT-PHOT production (and reduce
> > CPU load) if you are trying to make other measurements at the same
> > time in your simulation.
> >
> > 2) when entering in the wavelengths for the OPT-PROP (sdum
> WV-LIMITS)
> > cards you want them in cm, they were all entered in as meters, you
> > technically don't need them to run the simulation
> >
> > 3) same issue with OPT-PROD (sdum SCINT-WV), it was also in meters.
> > This card is very much needed.
> >
> > 4) you explained you want an electron beam. But you have the source
> > card set to isotropic at 20 MeV. :-)
> >
> > 5) You need to set an appropriate upper and lower energy in the
> USRBDX
> > card for binning the optical wavelengths. I usually use 1E-9 to
> 1E-8
> > initially.
> >
> > 6) You may want to set an OPT-PROP (sdum blank) card to set an index
> > of refraction, absorption and diffusion coefficients. This will
> > increase the processing time per particle however.
> >
> > Hope this helps.
> >
> > --Nick
> >
> > On Mon, Dec 10, 2012 at 7:11 AM, Vahan Petrosyan
> > <vpetrosyan_at_asls.candle.am <mailto:vpetrosyan_at_asls.candle.am>
> <mailto:vpetrosyan_at_asls.candle.am
> <mailto:vpetrosyan_at_asls.candle.am>>> wrote:
> >
> > I want to get optical photons distribution (angular, energy)
> that
> > irradiated by Yag screen ( scintillator material Yag(Ce) yttrium
> > aluminum garnet ) via USRBDX card when electron beam with
> 5Mev energy
> > hits it.
> >
> http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf
> >
> > I want to define scintillator material properties with fluka
> using
> > OPT-PROD card with SCINT-WV option but fluka gives an error.
> > I have never done such kind of simulation and need experts
> help in
> > order
> > to do it correctly. I am attached my input file and ask you to
> > look and
> > say what I did wrong.
> >
> > btw. I played with parameters of my input file a little and
> > switched on
> > Cerenkov and Transition(TR-Radia) radiation and disable
> OPT-PROD card
> > with SCINT-WV in order to get some results but USRBDX plot
> with option
> > OPTIPHOT gives me an empty data range.
> >
> > can someone help me to define correct parameters for
> scintillator
> > material Yag(Ce) and get correct distribution (angular,
> energy) of
> > optical photons via USRBDX card?
> > http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet
> >
> > Thanks!
> >
> >
> >
>
>
> --------------010902060600070909020501
> Content-Type: text/html; charset=ISO-8859-1
> Content-Transfer-Encoding: 7bit
>
> <html>
> <head>
> <meta content="text/html; charset=ISO-8859-1"
> http-equiv="Content-Type">
> </head>
> <body text="#000000" bgcolor="#FFFFFF">
> <div class="moz-cite-prefix">Many thanks for your help Nick! 
> 
> I did corrections that you mention and was able to get nice
> USRBDX
> result for photons but not for optical ones.  It still
> zeroes.....(: 
> I am attaching new corrected input file. 
> 
> As I understood from your explanation card OPT-PROP (sdum
> WV-LIMITS) is optional and card OPT-PROP (sdum SENSITIV )
> is used
> only for reducing CPU load not for getting more precise
> results....? 
> I switched on also  OPT-PROD (sdum SCINT-WV and
> WV-LIMITS)  but
> without luck! 
> I am little confused about WHAT(2) ( "fraction") parameter of
> OPT-PROD( sdum SCINT-WV ), is it the same as <em
> style="font-weight: bold; font-style: normal; color:
> rgb(0, 0,
> 0); font-family: arial, sans-serif; font-size: small;
> font-variant: normal; letter-spacing: normal;
> line-height: 16px;
> orphans: 2; text-align: left; text-indent: 0px;
> text-transform:
> none; white-space: normal; widows: 2; word-spacing: 0px;
> -webkit-text-size-adjust: auto;
> -webkit-text-stroke-width: 0px;
> background-color: rgb(255, 255, 255);">Luminosity
> [photons/Mev]???? 
> 
> How can I get  USRBDX distribution of optical photons? 
> 
> With regards V. Petrosyan! 
> 
> On 10.12.2012 21:44, Nicholas Bolibruch wrote: 
> </div>
> <blockquote
> cite="mid:CAKT5vujSqavQWEx0TPcXoRaDDZV-3ewEgNw664fbA+DoWAHX8Q_at_mail.gmail.com
> <mailto:mid%3ACAKT5vujSqavQWEx0TPcXoRaDDZV-3ewEgNw664fbA%2BDoWAHX8Q_at_mail.gmail.com>"
> type="cite">Hi Vahan,  
> 
> Here are the issues with your input: 
> 
> 1) your source is at 20 MeV, and is isotropic.  If you
> test it
> with a beam, make sure its in the positive direction, and
> set the
> primaries lower as the processing time is ~2 hrs on an i7
> if it's
> made directional with your current settings.  Using
> the SENSITIV
> sdum on OPT-PROP you can probably lower the OPT-PHOT production
> (and reduce CPU load) if you are trying to make other
> measurements
> at the same time in your simulation. 
> 
> 2) when entering in the wavelengths for the OPT-PROP (sdum
> WV-LIMITS) cards you want them in cm, they were all entered
> in as
> meters, you technically don't need them to run the
> simulation 
> 
> 3) same issue with OPT-PROD (sdum SCINT-WV), it was also in
> meters.  This card is very much needed. 
> 
> 4) you explained you want an electron beam. But you have the
> source card set to isotropic at 20 MeV. :-) 
> 
> 5) You need to set an appropriate upper and lower energy in the
> USRBDX card for binning the optical wavelengths.  I
> usually use
> 1E-9 to 1E-8 initially. 
> 
> 6) You may want to set an OPT-PROP (sdum blank) card to set an
> index of refraction, absorption and diffusion coefficients.
> This
> will increase the processing time per particle however. 
> 
> Hope this helps. 
> 
> --Nick 
> 
> <div class="gmail_quote">
> On Mon, Dec 10, 2012 at 7:11 AM, Vahan Petrosyan <span
> dir="ltr"><<a moz-do-not-send="true"
> href="mailto:vpetrosyan_at_asls.candle.am
> <mailto:vpetrosyan_at_asls.candle.am>"
> target="_blank">vpetrosyan_at_asls.candle.am
> <mailto:vpetrosyan_at_asls.candle.am></a>>
> wrote: 
> <blockquote class="gmail_quote" style="margin:0 0 0
> .8ex;border-left:1px #ccc solid;padding-left:1ex">
> I want to get optical photons distribution (angular,
> energy)
> that 
> irradiated by Yag screen ( scintillator material Yag(Ce)
> yttrium 
> aluminum garnet )  via USRBDX card when electron
> beam with
> 5Mev energy 
> hits it. 
>      <a moz-do-not-send="true"
> href="http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf"
>
> target="_blank">http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf</a> 
> 
> I want to define scintillator material properties with
> fluka
> using 
> OPT-PROD card with SCINT-WV option but fluka gives an
> error. 
> I have never done such kind of simulation and need experts
> help in order 
> to do it correctly. I am attached my input file and ask
> you to
> look and 
> say what I did wrong. 
> 
> btw. I played with parameters of my input file a little and
> switched on 
> Cerenkov and Transition(TR-Radia) radiation and disable
>  OPT-PROD card 
> with SCINT-WV in order to get some results but USRBDX plot
> with option 
> OPTIPHOT gives me an empty data range. 
> 
> can someone help me to define correct parameters for
> scintillator 
> material  Yag(Ce) and get correct
>  distribution (angular,
> energy) of 
> optical photons via USRBDX card? 
> <a moz-do-not-send="true"
>
> href="http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet"
>
> target="_blank">http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet</a> 
> 
> Thanks! 
> 
> </blockquote>
> </div>
> 
> 
> </blockquote>
> 
> </body>
> </html>
>
> --------------010902060600070909020501--
>
> --------------070706020302040203090104
> Content-Type: application/x-fluka;
> name="Yag.inp"
> Content-Transfer-Encoding: base64
> Content-Disposition: attachment;
> filename="Yag.inp"
>
> KiBEOlxmbHVrYVxmbHVwaXhcWWFnLmRhdAoqIENyZWF0ZWQ6IDUuMTIuMjAxMgoqIEF0OiAy
> MTo0MToxNwpUSVRMRQpZYWcgU2NyZWVuCkdMT0JBTCAgICAgICAgICAgICAgICAgICAgIDAu
> MCAgICAgICAwLjAgICAgICAgIDEuICAgICAgICAxLgpERUZBVUxUUyAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgUFJFQ0lT
> SU8KQkVBTSAgICAgICAgICAgLTAuMDIgICAgLTZFLTA1ICAgIC0wLjA2OCAgICAgIC0wLjIg
> ICAgICAtMC4yICAgICAgICAgIEVMRUNUUk9OCkJFQU1QT1MgICAgICAgICAgMC4wICAgICAg
> IDAuMCAgICAgICAxMC4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAgICBORUdBVElWRQpH
> RU9CRUdJTiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgQ09NQk5BTUUKICAgIDAgICAgMCAgICAgICAgICAgICAgICAgICBN
> Qy1DQUQKKiBDQ0QKUkNDIENDRCAgICAgICAgMC4wMCAtMS4wMCAwLjAwIDAuMDAgLTAuMjAg
> MC4wMCAwLjIwCiogU3BCClNQSCBTcEIgICAgICAgIDAuMDAgMC4wMCAwLjAwIDEwMDAuMDAK
> KiBTcFYKU1BIIFNwViAgICAgICAgMC4wMCAwLjAwIDAuMDAgOTUwLjAwCiogWWFnCkJPWCBZ
> YWcgICAgICAgIC0wLjUwIC0wLjM1IC0wLjM1IDEuMDAgMC4wMCAwLjAwIDAuMDAgMC43MSAw
> LjcxIDAuMDAgLTAuMDEgMC4wMQpFTkQKKiBSZWcgIyAxCiogQmxhY2tIOyBhc3NpZ25lZCBt
> YXRlcmlhbDogQmxhY2tob2xlOyBtYXQgIyAoMSkKQmxhY2tIICAgICAgIDUgK1NwQiAgLVNw
> VgoqIFJlZyAjIDIKKiBZYWc7IGFzc2lnbmVkIG1hdGVyaWFsOiBCbGFja2hvbGU7IG1hdCAj
> ICgxKQpZYWcgICAgICAgICAgNSArWWFnCiogUmVnICMgMwoqIENDRDsgYXNzaWduZWQgbWF0
> ZXJpYWw6IEJsYWNraG9sZTsgbWF0ICMgKDEpCkNDRCAgICAgICAgICA1ICtDQ0QKKiBSZWcg
> IyA0CiogVmFjdXVtOyBhc3NpZ25lZCBtYXRlcmlhbDogQmxhY2tob2xlOyBtYXQgIyAoMSkK
> VmFjdXVtICAgICAgIDUgK1NwViAgLVlhZyAgLUNDRApFTkQKR0VPRU5ECk1BVEVSSUFMICAg
> ICAgICAgMzMuICAgNzQuOTIxNiAgICAgIDUuNzMgICAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICBBUlNFTklDCiogR2xhc3MgbGVhZAoqCk1BVEVSSUFMICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgIDYuMjIgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBMRUFELUdMQQpD
> T01QT1VORCAgIC0wLjE1NjQ1MyAgICBPWFlHRU4gLTAuMDgwODY2ICAgU0lMSUNPTiAtMC4w
> MDgwOTIgIFRJVEFOSVVNTEVBRC1HTEEKQ09NUE9VTkQgICAtMC4wMDI2NTEgICBBUlNFTklD
> IC0wLjc1MTkzOCAgICAgIExFQUQgICAgICAgICAgICAgICAgICAgIExFQUQtR0xBCk1BVEVS
> SUFMICAgICAgICAgMzkuICA4OC45MDU4NSAgICAgNC40NjkgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICBZVFRSSVVNCk1BVEVSSUFMICAgICAgICAgNTguICAgMTQwLjExNSAgICAg
> Ni42NTcgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBDRVJJVU0KKiBZYWcgc2NpbnRp
> bGxhdG9yIHNjcmVlbgoqCk1BVEVSSUFMICAgICAgICAgICAgICAgICAgICAgICAgICAgIDQu
> NTUgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBZQUcKQ09NUE9VTkQgICAgICAgICAz
> LjAgICBZVFRSSVVNICAgICAgIDUuMCAgQUxVTUlOVU0gICAgICAxMi4wICAgIE9YWUdFTllB
> RwpDT01QT1VORCAgICAgICAgIDEuMCAgICBDRVJJVU0gICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgWUFHCkFTU0lHTk1BICAgIEJMQ0tIT0xFICAgIEJsYWNrSApB
> U1NJR05NQSAgICAgICAgIFlBRyAgICAgICBZYWcKQVNTSUdOTUEgICAgICBWQUNVVU0gICAg
> ICAgQ0NECkFTU0lHTk1BICAgICAgVkFDVVVNICAgIFZhY3V1bQpSQU5ET01JWiAgICAgICAg
> ICAxLgpQSE9UT05VQyAgICAgICAgICAxLiAgICAgICAgICAgICAgICAgICAgICAgIFZBQ1VV
> TSAgQExBU1RNQVQKT1BULVBST0QgICAgMi4wN0UtMTEgICAgICAxRS01ICAgICAgICAgICAg
> ICAgICBZQUcgICAgICAgWUFHICAgICAgICAgIENFUkVOS09WCk9QVC1QUk9EICAgICAwLjAx
> RTE1ICAgICAgMUUxNyAgICAgICAgICAgICAgICAgWUFHICAgICAgIFlBRyAgICAgICAgICBU
> Ui1SQURJQQpPUFQtUFJPRCAgICAgIDUuNUUtNSAgICAxNjcwMC4gICAgICAgICAgICAgICAg
> IFlBRyAgICAgICBZQUcgICAgICAgICAgU0NJTlQtV1YKT1BULVBST1AgICAgICAgIDJFLTUg
> ICAgNS41RS01ICAgIDEuMUUtNCAgICAgICBZQUcgICAgICAgWUFHICAgICAgICAgIFdWLUxJ
> TUlUCipPUFQtUFJPUCAgICAgICAgMS44MiAgICAgICAgICAgICAgICAgICAgICAgICAgIFlB
> RyAgICAgICBZQUcgICAgICAgICAgV1YtTElNSVQKVVNSQkRYICAgICAgICAgIDEwMS4gICAg
> UEhPVE9OICAgICAgLTU3LiAgICAgICBZYWcgICAgVmFjdXVtICAgICAgMi4wOFVzeFBoClVT
> UkJEWCAgICAgICAgICAxRS01ICAgICAgMUUtOSAgICAgIDEwMC4gICAgICAgICAgICAgICAg
> ICAgICAgICAgIDEwMC4mClVTUkJEWCAgICAgICAgICAgIDEuICAgIFBIT1RPTiAgICAgIC01
> Ny4gICAgICAgWWFnICAgIFZhY3V1bSAgICAgIDIuMDhVc3hDUGgKVVNSQkRYICAgICAgICAg
> ICAgICAgICAgICAgICAgICAgICAgMTAwLiAgICAgICAgICAgICAgICAgICAgICAgICAgMTAw
> LiYKVVNSQkRYICAgICAgICAgIDEwMS4gIE9QVElQSE9UICAgICAgLTU3LiAgICAgICBZYWcg
> ICAgVmFjdXVtICAgICAgMi4wOFVzeE9QClVTUkJEWCAgICAgICAgICAgICAgICAgICAgICAg
> ICAgICAgIDEwMC4gICAgICAgICAgICAgICAgICAgICAgICAgIDEwMC4mClNUQVJUICAgICAg
> ICAxMDAwMDAuClNUT1AK
> --------------070706020302040203090104--
>
>
>
>
>
> --
> --
> Nicholas Bolibruch
>
>
>

--------------080207000203040206030803
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit

<html>
 <head>
 <meta content="text/html; charset=ISO-8859-1"
 http-equiv="Content-Type">
 </head>
 <body text="#000000" bgcolor="#FFFFFF">
 <div class="moz-cite-prefix">Hi Nick. Thanks for your replay. It
 helped me much. 
 
 One question remains obscure for me. 
  1) How can one calculate WHAT(2) "fraction  of deposited energy
 going into i-th scin. photon emission"  parameter? 
  2) Should it be specified in scintillator material parameters or
 it can be calculated indirectly from other parameters too, for
 example luminosity and emission peak wavelength?  
 
 I supposed formula  Luminosity[ph/ev]*Energy(peak emission)[ev],
 am I right? 
 
 (I converted luminosity from [ph/Mev] to [ph/ev]) 
 
 
 On 12.12.2012 19:58, Nicholas Bolibruch wrote: 
 </div>
 <blockquote
cite="mid:CAKT5vuhfYOaz4Ai28MG7k=9Ruv2v2-odim5Y=tmpKjszLnz9zw_at_mail.gmail.com"
 type="cite">Hi Vahan, sorry I should have been more clear.  
 
 >5) You need to set an appropriate upper and lower energy in
 the USRBDX card >for binning the optical wavelengths.  I
 usually use 1E-9 to 1E-8 initially. 
 
 The USRBDX card for OPTIPHO needs to have the binning range set to
 appropriate energies. 
 
 
 >6) You may want to set an OPT-PROP (sdum blank) card to set an
 index of >refraction, absorption and diffusion coefficients.
 This will increase the >processing time per particle however. 
 
 This should read "You want to set an OPT-PROP ...". From this card
 you need to set _at least_ the absorption coefficient and the
 material you're assigning the properties to (ie. YAG).  Setting
 diffusion coefficients will increase the time per particle. 
 
 Attached is your last INP file modified to work.  I include
 comments on cards I've changed. 
 
 --Nick 
 
 
 <div class="gmail_quote">On Wed, Dec 12, 2012 at 5:23 AM, Vahan
 Petrosyan <<a moz-do-not-send="true"
 href="mailto:vpetrosyan_at_asls.candle.am" target="_blank">vpetrosyan_at_asls.candle.am</a>>
 wrote: 
 <blockquote class="gmail_quote" style="margin:0 0 0
 .8ex;border-left:1px #ccc solid;padding-left:1ex"> 
 <div class="im">
 Many thanks for your help Nick! 
 
 I did corrections that you mention and was able to get nice
 USRBDX 
 result for photons but not for optical ones.  It still
 zeroes.....(: 
 I am attaching new corrected input file. 
 
 As I understood from your explanation card OPT-PROP (sdum
 WV-LIMITS) is 
 optional and card OPT-PROP (sdum SENSITIV ) is used only for
 reducing 
 CPU load not for getting more precise results....? 
 I switched on also  OPT-PROD (sdum SCINT-WV and WV-LIMITS)
  but without 
 luck! 
 I am little confused about WHAT(2) ( "fraction") parameter
 of OPT-PROD( 
 </div>
 sdum SCINT-WV ), is it the same as /Luminosity
 [photons/Mev]???? / 
 <div>
 <div class="h5"> 
 How can I get  USRBDX distribution of optical photons? 
 
 With regards V. Petrosyan! 
 
 On 10.12.2012 21:44, Nicholas Bolibruch wrote: 
 > Hi Vahan, 
 > 
 > Here are the issues with your input: 
 > 
 > 1) your source is at 20 MeV, and is isotropic.  If
 you test it with a 
 > beam, make sure its in the positive direction, and
 set the primaries 
 > lower as the processing time is ~2 hrs on an i7 if
 it's made 
 > directional with your current settings.  Using the
 SENSITIV sdum on 
 > OPT-PROP you can probably lower the OPT-PHOT
 production (and reduce 
 > CPU load) if you are trying to make other
 measurements at the same 
 > time in your simulation. 
 > 
 > 2) when entering in the wavelengths for the OPT-PROP
 (sdum WV-LIMITS) 
 > cards you want them in cm, they were all entered in
 as meters, you 
 > technically don't need them to run the simulation 
 > 
 > 3) same issue with OPT-PROD (sdum SCINT-WV), it was
 also in meters. 
 > This card is very much needed. 
 > 
 > 4) you explained you want an electron beam. But you
 have the source 
 > card set to isotropic at 20 MeV. :-) 
 > 
 > 5) You need to set an appropriate upper and lower
 energy in the USRBDX 
 > card for binning the optical wavelengths.  I usually
 use 1E-9 to 1E-8 
 > initially. 
 > 
 > 6) You may want to set an OPT-PROP (sdum blank) card
 to set an index 
 > of refraction, absorption and diffusion coefficients.
 This will 
 > increase the processing time per particle however. 
 > 
 > Hope this helps. 
 > 
 > --Nick 
 > 
 > On Mon, Dec 10, 2012 at 7:11 AM, Vahan Petrosyan 
 </div>
 </div>
 <div>
 <div class="h5">> <<a moz-do-not-send="true"
 href="mailto:vpetrosyan_at_asls.candle.am">vpetrosyan_at_asls.candle.am</a>
 <mailto:<a moz-do-not-send="true"
 href="mailto:vpetrosyan_at_asls.candle.am">vpetrosyan_at_asls.candle.am</a>>>
 wrote: 
 > 
 >     I want to get optical photons distribution
 (angular, energy) that 
 >     irradiated by Yag screen ( scintillator material
 Yag(Ce) yttrium 
 >     aluminum garnet )  via USRBDX card when electron
 beam with 5Mev energy 
 >     hits it. 
 >     <a moz-do-not-send="true"
href="http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf"

target="_blank">http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf</a> 
 > 
 >     I want to define scintillator material properties
 with fluka using 
 >     OPT-PROD card with SCINT-WV option but fluka
 gives an error. 
 >     I have never done such kind of simulation and
 need experts help in 
 >     order 
 >     to do it correctly. I am attached my input file
 and ask you to 
 >     look and 
 >     say what I did wrong. 
 > 
 >     btw. I played with parameters of my input file a
 little and 
 >     switched on 
 >     Cerenkov and Transition(TR-Radia) radiation and
 disable  OPT-PROD card 
 >     with SCINT-WV in order to get some results but
 USRBDX plot with option 
 >     OPTIPHOT gives me an empty data range. 
 > 
 >     can someone help me to define correct parameters
 for scintillator 
 >     material  Yag(Ce) and get correct  distribution
 (angular, energy) of 
 >     optical photons via USRBDX card? 
 >     <a moz-do-not-send="true"
 href="http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet"
 target="_blank">http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet</a> 
 > 
 >     Thanks! 
 > 
 > 
 > 
 
 
 </div>
 </div>
 --------------010902060600070909020501 
 Content-Type: text/html; charset=ISO-8859-1 
 Content-Transfer-Encoding: 7bit 
 
 <html> 
    <head> 
      <meta content="text/html; charset=ISO-8859-1" 
        http-equiv="Content-Type"> 
    </head> 
    <body text="#000000" bgcolor="#FFFFFF"> 
      <div class="moz-cite-prefix">Many thanks for your
 help Nick! 
         
 <div class="im">       I did corrections that you mention and
 was able to get nice USRBDX 
 </div>
        result for photons but not for optical ones.&nbsp;
 It still 
        zeroes.....(: 
        I am attaching new corrected input file. 
         
 <div class="im">       As I understood from your explanation
 card OPT-PROP (sdum 
        WV-LIMITS) is optional and card OPT-PROP (sdum
 SENSITIV ) is used 
        only for reducing CPU load not for getting more
 precise 
 </div>
        results....? 
        I switched on also&nbsp; OPT-PROD (sdum SCINT-WV
 and WV-LIMITS)&nbsp; but 
        without luck! 
 <div class="im">       I am little confused about WHAT(2) (
 "fraction") parameter of 
 </div>
        OPT-PROD( sdum SCINT-WV ), is it the same as <em 
          style="font-weight: bold; font-style: normal; color:
 rgb(0, 0, 
          0); font-family: arial, sans-serif; font-size: small; 
          font-variant: normal; letter-spacing: normal;
 line-height: 16px; 
          orphans: 2; text-align: left; text-indent: 0px;
 text-transform: 
          none; white-space: normal; widows: 2; word-spacing:
 0px; 
          -webkit-text-size-adjust: auto;
 -webkit-text-stroke-width: 0px; 
          background-color: rgb(255, 255, 255);">Luminosity 
          [photons/Mev]???? 
         
        How can I get&nbsp; USRBDX distribution of optical
 photons? 
         
        With regards V. Petrosyan! 
         
        On 10.12.2012 21:44, Nicholas Bolibruch
 wrote: 
      </div> 
      <blockquote 
 cite="<a moz-do-not-send="true"
href="mailto:mid%3ACAKT5vujSqavQWEx0TPcXoRaDDZV-3ewEgNw664fbA%2BDoWAHX8Q_at_mail.gmail.com">mid:CAKT5vujSqavQWEx0TPcXoRaDDZV-3ewEgNw664fbA+DoWAHX8Q_at_mail.gmail.com</a>" 
        type="cite">Hi Vahan,&nbsp; 
         
        Here are the issues with your input: 
         
        1) your source is at 20 MeV, and is
 isotropic.&nbsp; If you test it 
 <div class="im">       with a beam, make sure its in the
 positive direction, and set the 
        primaries lower as the processing time is ~2 hrs on
 an i7 if it's 
 </div>
        made directional with your current settings.&nbsp;
 Using the SENSITIV 
 <div class="im">       sdum on OPT-PROP you can probably lower
 the OPT-PHOT production 
        (and reduce CPU load) if you are trying to make other
 measurements 
 </div>
        at the same time in your simulation. 
         
 <div class="im">       2) when entering in the wavelengths for
 the OPT-PROP (sdum 
        WV-LIMITS) cards you want them in cm, they were all
 entered in as 
 </div>
        meters, you technically don't need them to run the
 simulation 
         
 <div class="im">       3) same issue with OPT-PROD (sdum
 SCINT-WV), it was also in 
 </div>
        meters.&nbsp; This card is very much needed.
 
         
 <div class="im">       4) you explained you want an electron
 beam. But you have the 
 </div>
        source card set to isotropic at 20 MeV. :-) 
         
 <div class="im">       5) You need to set an appropriate upper
 and lower energy in the 
 </div>
        USRBDX card for binning the optical
 wavelengths.&nbsp; I usually use 
        1E-9 to 1E-8 initially. 
         
 <div class="im">       6) You may want to set an OPT-PROP
 (sdum blank) card to set an 
        index of refraction, absorption and diffusion
 coefficients. This 
 </div>
        will increase the processing time per particle
 however. 
         
        Hope this helps. 
         
        --Nick 
         
        <div class="gmail_quote"> 
          On Mon, Dec 10, 2012 at 7:11 AM, Vahan Petrosyan
 <span 
            dir="ltr">&lt;<a moz-do-not-send="true" 
              href="mailto:<a moz-do-not-send="true"
 href="mailto:vpetrosyan_at_asls.candle.am">vpetrosyan_at_asls.candle.am</a>"
 target="_blank"><a moz-do-not-send="true"
 href="mailto:vpetrosyan_at_asls.candle.am">vpetrosyan_at_asls.candle.am</a></a>&gt; 
          wrote: 
          <blockquote class="gmail_quote" style="margin:0 0
 0 
            .8ex;border-left:1px #ccc
 solid;padding-left:1ex"> 
 <div class="im">           I want to get optical photons
 distribution (angular, energy) 
 </div>
            that 
 <div class="im">           irradiated by Yag screen (
 scintillator material Yag(Ce) 
 </div>
            yttrium 
            aluminum garnet ) &nbsp;via USRBDX card when
 electron beam with 
            5Mev energy 
            hits it. 
            &nbsp; &nbsp; &nbsp;<a
 moz-do-not-send="true" 
 href="<a moz-do-not-send="true"
href="http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf"

target="_blank">http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf</a>" 
 
 target="_blank"><a moz-do-not-send="true"
href="http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf"

target="_blank">http://www.detectors.saint-gobain.com/uploadedFiles/SGdetectors/Documents/Product_Data_Sheets/YAG-Data-Sheet.pdf</a></a> 
             
 <div class="im">           I want to define scintillator
 material properties with fluka 
 </div>
            using 
            OPT-PROD card with SCINT-WV option but fluka gives
 an error. 
 <div class="im">           I have never done such kind of
 simulation and need experts 
 </div>
            help in order 
 <div class="im">           to do it correctly. I am attached
 my input file and ask you to 
 </div>
            look and 
            say what I did wrong. 
             
 <div class="im">           btw. I played with parameters of my
 input file a little and 
 </div>
            switched on 
 <div class="im">           Cerenkov and Transition(TR-Radia)
 radiation and disable 
 </div>
            &nbsp;OPT-PROD card 
 <div class="im">           with SCINT-WV in order to get some
 results but USRBDX plot 
 </div>
            with option 
            OPTIPHOT gives me an empty data range. 
             
 <div class="im">           can someone help me to define
 correct parameters for 
 </div>
            scintillator 
            material &nbsp;Yag(Ce) and get correct
 &nbsp;distribution (angular, 
            energy) of 
            optical photons via USRBDX card? 
            <a moz-do-not-send="true" 
              href="<a moz-do-not-send="true"
 href="http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet"
 target="_blank">http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet</a>" 
              target="_blank"><a moz-do-not-send="true"
 href="http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet"
 target="_blank">http://en.wikipedia.org/wiki/Yttrium_aluminum_garnet</a></a> 
             
            Thanks! 
             
          </blockquote> 
        </div> 
         
         
      </blockquote> 
       
    </body> 
 </html> 
 
 --------------010902060600070909020501-- 
 
 --------------070706020302040203090104 
 Content-Type: application/x-fluka; 
   name="Yag.inp" 
 Content-Transfer-Encoding: base64 
 Content-Disposition: attachment; 
   filename="Yag.inp" 
 
KiBEOlxmbHVrYVxmbHVwaXhcWWFnLmRhdAoqIENyZWF0ZWQ6IDUuMTIuMjAxMgoqIEF0OiAy 
MTo0MToxNwpUSVRMRQpZYWcgU2NyZWVuCkdMT0JBTCAgICAgICAgICAgICAgICAgICAgIDAu 
MCAgICAgICAwLjAgICAgICAgIDEuICAgICAgICAxLgpERUZBVUxUUyAgICAgICAgICAgICAg 
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgUFJFQ0lT 
SU8KQkVBTSAgICAgICAgICAgLTAuMDIgICAgLTZFLTA1ICAgIC0wLjA2OCAgICAgIC0wLjIg 
ICAgICAtMC4yICAgICAgICAgIEVMRUNUUk9OCkJFQU1QT1MgICAgICAgICAgMC4wICAgICAg 
IDAuMCAgICAgICAxMC4gICAgICAgMC4wICAgICAgIDAuMCAgICAgICAgICBORUdBVElWRQpH 
RU9CRUdJTiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg 
ICAgICAgICAgICAgICAgQ09NQk5BTUUKICAgIDAgICAgMCAgICAgICAgICAgICAgICAgICBN 
Qy1DQUQKKiBDQ0QKUkNDIENDRCAgICAgICAgMC4wMCAtMS4wMCAwLjAwIDAuMDAgLTAuMjAg 
MC4wMCAwLjIwCiogU3BCClNQSCBTcEIgICAgICAgIDAuMDAgMC4wMCAwLjAwIDEwMDAuMDAK 
KiBTcFYKU1BIIFNwViAgICAgICAgMC4wMCAwLjAwIDAuMDAgOTUwLjAwCiogWWFnCkJPWCBZ 
YWcgICAgICAgIC0wLjUwIC0wLjM1IC0wLjM1IDEuMDAgMC4wMCAwLjAwIDAuMDAgMC43MSAw 
LjcxIDAuMDAgLTAuMDEgMC4wMQpFTkQKKiBSZWcgIyAxCiogQmxhY2tIOyBhc3NpZ25lZCBt 
YXRlcmlhbDogQmxhY2tob2xlOyBtYXQgIyAoMSkKQmxhY2tIICAgICAgIDUgK1NwQiAgLVNw 
VgoqIFJlZyAjIDIKKiBZYWc7IGFzc2lnbmVkIG1hdGVyaWFsOiBCbGFja2hvbGU7IG1hdCAj 
ICgxKQpZYWcgICAgICAgICAgNSArWWFnCiogUmVnICMgMwoqIENDRDsgYXNzaWduZWQgbWF0 
ZXJpYWw6IEJsYWNraG9sZTsgbWF0ICMgKDEpCkNDRCAgICAgICAgICA1ICtDQ0QKKiBSZWcg 
IyA0CiogVmFjdXVtOyBhc3NpZ25lZCBtYXRlcmlhbDogQmxhY2tob2xlOyBtYXQgIyAoMSkK 
VmFjdXVtICAgICAgIDUgK1NwViAgLVlhZyAgLUNDRApFTkQKR0VPRU5ECk1BVEVSSUFMICAg 
ICAgICAgMzMuICAgNzQuOTIxNiAgICAgIDUuNzMgICAgICAgICAgICAgICAgICAgICAgICAg 
ICAgICBBUlNFTklDCiogR2xhc3MgbGVhZAoqCk1BVEVSSUFMICAgICAgICAgICAgICAgICAg 
ICAgICAgICAgIDYuMjIgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBMRUFELUdMQQpD 
T01QT1VORCAgIC0wLjE1NjQ1MyAgICBPWFlHRU4gLTAuMDgwODY2ICAgU0lMSUNPTiAtMC4w 
MDgwOTIgIFRJVEFOSVVNTEVBRC1HTEEKQ09NUE9VTkQgICAtMC4wMDI2NTEgICBBUlNFTklD 
IC0wLjc1MTkzOCAgICAgIExFQUQgICAgICAgICAgICAgICAgICAgIExFQUQtR0xBCk1BVEVS 
SUFMICAgICAgICAgMzkuICA4OC45MDU4NSAgICAgNC40NjkgICAgICAgICAgICAgICAgICAg 
ICAgICAgICAgICBZVFRSSVVNCk1BVEVSSUFMICAgICAgICAgNTguICAgMTQwLjExNSAgICAg 
Ni42NTcgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBDRVJJVU0KKiBZYWcgc2NpbnRp 
bGxhdG9yIHNjcmVlbgoqCk1BVEVSSUFMICAgICAgICAgICAgICAgICAgICAgICAgICAgIDQu 
NTUgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBZQUcKQ09NUE9VTkQgICAgICAgICAz 
LjAgICBZVFRSSVVNICAgICAgIDUuMCAgQUxVTUlOVU0gICAgICAxMi4wICAgIE9YWUdFTllB 
RwpDT01QT1VORCAgICAgICAgIDEuMCAgICBDRVJJVU0gICAgICAgICAgICAgICAgICAgICAg 
ICAgICAgICAgICAgICAgICAgWUFHCkFTU0lHTk1BICAgIEJMQ0tIT0xFICAgIEJsYWNrSApB 
U1NJR05NQSAgICAgICAgIFlBRyAgICAgICBZYWcKQVNTSUdOTUEgICAgICBWQUNVVU0gICAg 
ICAgQ0NECkFTU0lHTk1BICAgICAgVkFDVVVNICAgIFZhY3V1bQpSQU5ET01JWiAgICAgICAg 
ICAxLgpQSE9UT05VQyAgICAgICAgICAxLiAgICAgICAgICAgICAgICAgICAgICAgIFZBQ1VV 
TSAgQExBU1RNQVQKT1BULVBST0QgICAgMi4wN0UtMTEgICAgICAxRS01ICAgICAgICAgICAg 
ICAgICBZQUcgICAgICAgWUFHICAgICAgICAgIENFUkVOS09WCk9QVC1QUk9EICAgICAwLjAx 
RTE1ICAgICAgMUUxNyAgICAgICAgICAgICAgICAgWUFHICAgICAgIFlBRyAgICAgICAgICBU 
Ui1SQURJQQpPUFQtUFJPRCAgICAgIDUuNUUtNSAgICAxNjcwMC4gICAgICAgICAgICAgICAg 
IFlBRyAgICAgICBZQUcgICAgICAgICAgU0NJTlQtV1YKT1BULVBST1AgICAgICAgIDJFLTUg 
ICAgNS41RS01ICAgIDEuMUUtNCAgICAgICBZQUcgICAgICAgWUFHICAgICAgICAgIFdWLUxJ 
TUlUCipPUFQtUFJPUCAgICAgICAgMS44MiAgICAgICAgICAgICAgICAgICAgICAgICAgIFlB 
RyAgICAgICBZQUcgICAgICAgICAgV1YtTElNSVQKVVNSQkRYICAgICAgICAgIDEwMS4gICAg 
UEhPVE9OICAgICAgLTU3LiAgICAgICBZYWcgICAgVmFjdXVtICAgICAgMi4wOFVzeFBoClVT 
UkJEWCAgICAgICAgICAxRS01ICAgICAgMUUtOSAgICAgIDEwMC4gICAgICAgICAgICAgICAg 
ICAgICAgICAgIDEwMC4mClVTUkJEWCAgICAgICAgICAgIDEuICAgIFBIT1RPTiAgICAgIC01 
Ny4gICAgICAgWWFnICAgIFZhY3V1bSAgICAgIDIuMDhVc3hDUGgKVVNSQkRYICAgICAgICAg 
ICAgICAgICAgICAgICAgICAgICAgMTAwLiAgICAgICAgICAgICAgICAgICAgICAgICAgMTAw 
LiYKVVNSQkRYICAgICAgICAgIDEwMS4gIE9QVElQSE9UICAgICAgLTU3LiAgICAgICBZYWcg 
ICAgVmFjdXVtICAgICAgMi4wOFVzeE9QClVTUkJEWCAgICAgICAgICAgICAgICAgICAgICAg 
ICAgICAgIDEwMC4gICAgICAgICAgICAgICAgICAgICAgICAgIDEwMC4mClNUQVJUICAgICAg 
 ICAxMDAwMDAuClNUT1AK 
 --------------070706020302040203090104-- 
 
 
 </blockquote>
 </div>
 
 
 
 -- 
 -- 
 Nicholas Bolibruch 
 
 
 
 </blockquote>
 
 </body>
</html>

--------------080207000203040206030803--
Received on Mon Dec 17 2012 - 09:45:09 CET

This archive was generated by hypermail 2.2.0 : Mon Dec 17 2012 - 09:45:23 CET