Re: Optical Photon Production by Scintillator Material, need help

From: Vahan Petrosyan <vpetrosyan_at_asls.candle.am>
Date: Wed, 12 Dec 2012 15:23:18 +0400

<CAKT5vujSqavQWEx0TPcXoRaDDZV-3ewEgNw664fbA+DoWAHX8Q_at_mail.gmail.com>
In-Reply-To: <CAKT5vujSqavQWEx0TPcXoRaDDZV-3ewEgNw664fbA+DoWAHX8Q_at_mail.gmail.com>
Content-Type: multipart/mixed;
  boundary="------------070706020302040203090104"
X-Antivirus: avast! (VPS 121211-2, 12.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]); Wed, 12 Dec 2012
12:23:37 +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.
--------------070706020302040203090104
Content-Type: multipart/alternative;
  boundary="------------010902060600070909020501"

--------------010902060600070909020501
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

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>> 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.&nbsp; 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&nbsp; OPT-PROD (sdum SCINT-WV and WV-LIMITS)&nbsp; 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&nbsp; 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"
       type="cite">Hi Vahan,&nbsp; <br>
       <br>
       Here are the issues with your input:<br>
       <br>
       1) your source is at 20 MeV, and is isotropic.&nbsp; 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.&nbsp; 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.&nbsp; 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.&nbsp; 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">&lt;<a moz-do-not-send="true"
             href="mailto:vpetrosyan_at_asls.candle.am" target="_blank">vpetrosyan_at_asls.candle.am</a>&gt;</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 ) &nbsp;via USRBDX card when electron beam with
           5Mev energy<br>
           hits it.<br>
           &nbsp; &nbsp; &nbsp;<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
           &nbsp;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 &nbsp;Yag(Ce) and get correct &nbsp;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--
Received on Wed Dec 12 2012 - 13:10:17 CET

This archive was generated by hypermail 2.2.0 : Wed Dec 12 2012 - 13:10:32 CET