A weight Problem in uniform sampling

From: SEMON Laurent (laurent.semon@synchrotron-soleil.fr)
Date: Thu Jun 21 2007 - 09:17:50 CEST

  • Next message: SEMON Laurent: "Option EMF-BIAS"

    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&#8217;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'>&nbsp;</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'>&nbsp;</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'>&nbsp;</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'>&nbsp;</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 =
    &nbsp;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'>&nbsp;</span></font></p>

    <p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
    style=3D'font-size:
    10.0pt;font-family:Arial'>&nbsp;</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 =
    &#8230;</span></font></p>

    <p class=3DMsoNormal><font size=3D2 face=3DArial><span lang=3DEN-GB =
    style=3D'font-size:
    10.0pt;font-family:Arial'>&nbsp;</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'>&nbsp;</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--


  • Next message: SEMON Laurent: "Option EMF-BIAS"

    This archive was generated by hypermail 2.1.6 : Wed Jun 27 2007 - 14:32:34 CEST