RE: [fluka-discuss]: Multiple sources

From: Eleftherios Skordis <eleftherios.skordis_at_cern.ch>
Date: Sat, 11 May 2019 09:13:15 +0000

Dear Jyoti,

at the moment it is not possible to define so many different sources with cards.
The source routine gives you the freedom to define both the energy, the position and the direction of the particles that are generated.

Kind regards

Lefteris


-----------------------

Eleftherios Skordis
Dep. EN/STI, CERN
CH-1211 GENEVA 23
SWITZERLAND

OFFICE: +41-22-7679541<tel:%2B41-22-7675461>
________________________________
From: jyoti.grg_at_britatom.gov.in [jyoti.grg_at_britatom.gov.in]
Sent: 10 May 2019 13:34
To: Eleftherios Skordis
Cc: fluka-discuss_at_fluka.org
Subject: RE: [fluka-discuss]: Multiple sources

Dear Lefteris,

Thank you. It worked.
Can we do the same without using subroutine? i.e., defining multiple sources without using source file.

With Regards,
Jyoti Garg


On 05/08/19 11:17 PM, Eleftherios Skordis <eleftherios.skordis_at_cern.ch> wrote:

Dear Jyoti,

defining more than one source, basically means that you are just sampling primary particles from another set of coordinates and/or energies.
That means that the primary particle will originate from either of the sources with a chance that you set.

There are many implementations you can do.
e.g.
RNDNUM = FLRNDM(DUMMY)
IF ( RNDNUM .GT. 0.5 ) THEN
XORIG = 10
YORIG = 10
RMIN = 0
RMAX = 4
ZMIN = 10
ZMAX = 10
ELSE
           XORIG = 0
YORIG = 0
RMIN = 0
RMAX = 5
ZMIN = -10
ZMAX = +10
END IF
R=SQRT(FLRNDM(DUMMY)) * (RMAX-RMIN) + RMIN
THETA = FLRNDM(DUMMY) * TWOPIP
    x = R * cos(THETA) + XORIG
y = R * sin(THETA) + YORIG
z = ZMIN + (ZMAX - ZMIN) * FLRNDM(DUMMY)
      XFLK (NPFLKA) =x
      YFLK (NPFLKA) =y
      ZFLK (NPFLKA) =z
WRITE(*,*) x,y,z


Unless you are interested in sampling more than one particles in the primary particle stack but that should be of interest only in the case of event by event scoring.

Kind regards

Lefteris




-----------------------

Eleftherios Skordis
Dep. EN/STI, CERN
CH-1211 GENEVA 23
SWITZERLAND

OFFICE: +41-22-7679541<tel:%2B41-22-7675461>
________________________________
From: jyoti.grg_at_britatom.gov.in [jyoti.grg_at_britatom.gov.in]
Sent: 08 May 2019 13:11
To: Eleftherios Skordis; fluka-discuss_at_fluka.org
Subject: RE: [fluka-discuss]: Multiple sources


Dear Lefteris
Can you please help me to describe two sources in subroutine. Then I will write it for nine.
I am facing problem to define more than one source.
I have attached the subroutine file for a single source. Kindly find the attached file.

Thanks and Regards,
Jyoti Garg


On 05/08/19 03:50 PM, Eleftherios Skordis <eleftherios.skordis_at_cern.ch> wrote:

Dear Jyoti

The easiest way I can imagine is that you can run 9 different simulations with the 9 sources and then add/normalize the dose in post processing.

Otherwise, you can describe the 9 sources in your subroutine and then just use a random number generator to sample (depending on your sources) each different source.

The results are equivalent if normalized correctly.

Hope this helps

Kind regards

Lefteris








-----------------------

Eleftherios Skordis
Dep. EN/STI, CERN
CH-1211 GENEVA 23
SWITZERLAND

OFFICE: +41-22-7679541<tel:%2B41-22-7675461>
________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Jyoti [jyoti.grg_at_britatom.gov.in]
Sent: 08 May 2019 10:16
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: Multiple sources


Dear FLUKA experts,
I have nine Co-60 sources (cylinderical) and I need to calculate dose at a point from these sources using FLUKA.
I am able to design a single cylinderical source using source subroutine.
How to design multiple sources? Please guide.

With Regards,
Jyoti Garg




__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Mon May 13 2019 - 23:35:17 CEST

This archive was generated by hypermail 2.3.0 : Tue May 14 2019 - 00:56:49 CEST