RE: [fluka-discuss]: about annular cylindrical source

From: Joachim Vollaire <joachim.vollaire_at_cern.ch>
Date: Wed, 24 May 2017 06:06:26 +0000

Dear Arghya,

Just to complement the answer from Luigi, I would like to draw to your attention on the fact that FLUKA allows you to simulate the decay of a radionuclide (Ir-192 in your case) as a "source/primary particle". And again this is treated in the advanced course material in the presentation dealing with radioactivity / activation. You have an example for Co-60 in the course but you just need to change it to Ir 192 (corresponding Z and A)

https://indico.cern.ch/event/489973/contributions/2000420/attachments/1272671/1886660/14_AdvancedActivation2016.pdf

Note that you need the RADDECAY/DCYSCORE card activated in semi-analogue mode to call for the decay of the isotope and request the detectors to be associated to the decay... In that case the normalization will be per Bq. If you have already implemented your source routine maybe it is worth comparing the two methods.

Greetings
Joachim
 


-----Original Message-----
From: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of Luigi Salvatore Esposito
Sent: 22 May 2017 21:22
To: arghyac_at_barc.gov.in
Cc: fluka <fluka-discuss_at_fluka.org>; owner-fluka-discuss_at_mi.infn.it
Subject: Re: [fluka-discuss]: about annular cylindrical source

Dear Arghya,
To sample from an energy spectrum, you can refer to the same slides of my previous email: there is an example from page 36 and following.
About the position/direction cosines, you can adapt the following pseudocode (assuming that the cylinder axis is along the z-axis), to be inserted in the source.f.

* Cosines (tx,ty,tz)
      TXFLK (NPFLKA) = TXX
      TYFLK (NPFLKA) = TYY
      TZFLK (NPFLKA) = TZZ
* Particle coordinates
      R = SQRT(FLRNDM(DUMMY)) * (RMAX-RMIN) + RMIN
      theta = FLRNDM(DUMMY) * TWOPIP
      XFLK (NPFLKA) = R * cos( theta )
      YFLK (NPFLKA) = R * sin( theta )
      ZFLK (NPFLKA) = ZMIN + (ZMAX - ZMIN) * FLRNDM(DUMMY)

In the same slides, some useful auxiliary routines are documented at page 31 and following.
Hope this could help.
Best regards, luigi

On 22.05.17, 13:40, "arghyac_at_barc.gov.in" <arghyac_at_barc.gov.in> wrote:

      Thank you very much for your quick response.
    
        But in my case I cannot use Beampos card as I have to feed energy
    spectrum also. So I need to sample both annular cylindrical source
    geometry and energy spectrum using source.f subroutine only.
      Please suggest me how to sample both using source.f subroutine.
    
    Thanking you in advance.
    
    
    
    On 2017-05-22 15:37, Luigi Salvatore Esposito wrote:
> Dear Arghya,
> You can use BEAMPOS card with SDUM = CYLI-VOL.
> More details in the manual and in the example given in these slides
> (from page 8 to 11)
> https://indico.cern.ch/event/489973/contributions/2148893/attachments/1271819/1884943/11_AdvancedSources2016.pdf
> Best regards, luigi
>
> On 22.05.17, 07:44, "owner-fluka-discuss_at_mi.infn.it on behalf of
> arghyac_at_barc.gov.in" <owner-fluka-discuss_at_mi.infn.it on behalf of
> arghyac_at_barc.gov.in> wrote:
>
> Dear experts,
>
> I need to design one annular cylindrical source (Rmin = 2 cm
> and
> Rmax = 5 cm) in which source material will be within Rmin and Rmax
> i.e
> annular region. Source.f card will be for Ir-192 spectrum.
>
> How to sample this kind of Cylindrical source in source.f card.
> Please give me some suggestions on this.
>
> Thank you in advance.
>
> With regards,
> Arghya Chattaraj
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id=acc_info
>
>
>
> ��������������������������������������������������������(�Ƨ��ځ맲��r���*'~�&�?�����x%��Z�?m����?0��n���?ߖ�?��i�'i���w�
    
    

b��jy����{.n�+��b�w�ka�ɚ�X��X���a��l��0��n��+���F����i�"��
(Ƨځ맲r*'~&x%Zm 0nߖi'iw
Received on Wed May 24 2017 - 09:05:36 CEST

This archive was generated by hypermail 2.3.0 : Wed May 24 2017 - 09:05:38 CEST