From: SEMON Laurent (laurent.semon@synchrotron-soleil.fr)
Date: Thu Jun 21 2007 - 09:17:50 CEST
This is a multi-part message in MIME format.
------_=_NextPart_001_01C7B3D4.46B4C163
Content-Type: text/plain;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable
Hello dear's users of Fluka,
=20
I try to sample a synchrotron source for Fluka calculations and i use
for that source.f . I perform an uniform sampling over my energy range
[1:500] keV and this sampling is very efficient (I manage to sample all
the energies).=20
My problem comes from the weight calculation. I use the formula w =3D
N(E)/A*(Emax-Emin) where A is the integral of the spectra over the
energy interval [Emin:Emax] and N(E) the number of photons per second
per keV as function of energy. For the lowest energies (for example 2
keV) N(E)=3D10E+20 and w=3D145 but for the highest energies (500 keV for
example) N(E)=3D10E-10 and w=3D10E-24!!!
With this huge weight interval Fluka calculations never converge.
=20
My questions are:
=20
1) Is it possible to perform good calculations with this huge weight
interval even if the spectra sampling is good?
=20
2) If an uniform sampling is not desirable for this type of spectra
(because of the huge weight interval), do you recommends a better
function that 1/N(E) to biase the spectra ?
=20
=20
Thank you very much for all informations ...
=20
Best regards
=20
Laurent
------_=_NextPart_001_01C7B3D4.46B4C163
Content-Type: text/html;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable
<html>
<head>
<meta http-equiv=3DContent-Type content=3D"text/html; =
charset=3Dus-ascii">
<meta name=3DGenerator content=3D"Microsoft Word 10 (filtered)">
<style>
<!--
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline;}
span.StyleCourrierlectronique17
{font-family:Arial;
color:windowtext;}
@page Section1
{size:595.3pt 841.9pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.Section1
{page:Section1;}
-->
</style>
</head>
<body lang=3DFR link=3Dblue vlink=3Dpurple>
<div class=3DSection1>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>Hello dear’s users of =
Fluka,</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>I try to sample a synchrotron source for Fluka
calculations and i use for that source.f . I perform an uniform sampling =
over
my energy range [1:500] keV and this sampling is very efficient (I =
manage to
sample all the energies). </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>My problem comes from the weight calculation. =
I use
the formula w =3D N(E)/A*(Emax-Emin) where A is the integral of the =
spectra over the
energy interval [Emin:Emax] and N(E) the number of photons per second =
per keV
as function of energy. For the lowest energies (for example 2 keV) =
N(E)=3D10E+20
and w=3D145 but for the highest energies (500 keV for example) =
N(E)=3D10E-10 and
w=3D10E-24!!!</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>With this huge weight interval Fluka =
calculations never
converge.</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>My questions are:</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>1) Is it possible to perform good calculations =
with
this huge weight interval even if the spectra sampling is =
good?</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>2) If an uniform sampling is not desirable for =
this
type of spectra (because of the huge weight interval), do you recommends =
a
better function that 1/N(E) to biase the spectra ?</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>Thank you very much for all informations =
…</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>Best regards</span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'> </span></font></p>
<p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
style=3D'font-size:
10.0pt;font-family:Arial'>Laurent</span></font></p>
</div>
</body>
</html>
------_=_NextPart_001_01C7B3D4.46B4C163--
This archive was generated by hypermail 2.1.6 : Wed Jun 27 2007 - 14:32:34 CEST