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