-- Stefan E. Mueller Department of Information Services and Computing - Computational Science and Institute of Radiation Physics Helmholtz-Zentrum Dresden-Rossendorf Tel: +49 (0351) 260 3847 Stefan.Mueller_at_hzdr.de http://www.hzdr.de Vorstand: Prof. Dr. Sebastian M. Schmidt, Dr. Diana Stiller Vereinsregister: VR 1693 beim Amtsgericht Dresden On Fri, 22 Jul 2022, emilio wrote: > Dear Paola, > > thank you very much for your answer. > > I tried using USRBDX, and it works. > > However... I cannot see any particle coming out.... > > My guess is that the radioactive decay never take place: I tried to replace ISOTOPES with ELECTRON, and I can actually see the particles, so it's not a problem of USRBDX... > > I only included the RADDECAY card, selecting the "semi-analogue" mode and leaving everything else empty; from my understanding this should be sufficient to have radioactive decay (without cards such as IRRPROFI or DCYTIMES, etc...), is it correct? > > RADDECAY 2. 0.0 > > In the attached files you can find the input card and an example of the output > > Thank you again for your help > Best Regards, > Emilio > > > > -----Original Messages----- > > From: "Paola Sala" <paola.sala_at_mi.infn.it> > > Sent Time: 2022-07-21 17:22:37 (Thursday) > > To: emilio <emilio_at_impcas.ac.cn>, fluka-discuss_at_fluka.org > > Cc: > > Subject: Re: [fluka-discuss]: Problem simulating radioactive decay of isotopes at rest > > > > Dear Emilio, > > > > The problem is in the combination of an ISOTOPE beam with the USRYIELD > > cards. USRYIELD calculates angular quantities with respect to the > > direction of the beam, and since some of the quantities can be > > calculated with respect to the center-of-mass frame, the code tries to > > initialize a center-of-mass momentum, which fails if the beam momentum > > is zero. > > > > I have a question and a solution, besides that this specific situation > > will be > > protected in a next respin > > > > Question: are you sure that you need USRYIELD instead of USRBDX? I > > mean, what is the purpose of calculating angles with respect to a > > fixed (beam) direction when source and geometry are perfectly > > isotropic? Why not to use USRBDX that calculates angles with respect > > to the normal to the boundary surface? (see note 1 to the USRYIELD > > description in the manual) > > > > Solution, in case you really need USRYIELD: add another USRYIELD card > > with SDUM=BEAMDEF (see again note 1) , in order to define a dummy beam > > just to initialize the scoring. For instance > > USRYIELD PROTON 1.0 1. 0. 0. 1.0 BEAMDEF > > > > Again from the fluka manual, USRYIELD description: > > For SDUM = BEAMDEF: > > WHAT(1) = projectile particle index, or corresponding name > > Default= IJBEAM (beam particle) > > WHAT(2)= target particle index, or corresponding name (used by the > > code to define the c.m.s. frame) > > Default : 1.0 (proton) > > WHAT(3) = projectile momentum > > Default= PBEAM (beam momentum) > > WHAT(4,5,6) = projectile direction cosines > > Default = UBEAM, VBEAM, WBEAM (beam direction cosines) > > > > Hope this helps > > Paola > > > > On 7/21/22 06:44, emilio wrote: > > > > > > Dear Fluka experts, > > > > > > > > > I am trying to simulate the decay of radioactive isotopes, however > > > when I try to run the got (via Flair), i got a time-out error and no > > > output is produced. > > > > > > > > > In the attached file you can find the input card; I do believe however > > > that the issue is in the primary particle generation. > > > > > > The radioactive isotopes are created with the following BEAM card > > > (which is similar to the one I found here: > > > https://indico.cern.ch/event/604018/contributions/2435668/attachments/1395470/2242783/19_Activation_2017.pdf) > > > > > > > > > BEAM ISOTOPE > > > HI-PROPE 34. 83. > > > > > > > > > The RADDECAY card is present as well (semi-analogue mode) > > > > > > > > > At first I thought the problem was that the isotopes are created at > > > rest, so I changed the energy to 1 keV, but I got the same result. > > > However if I try to change particle (for example, using electrons) the > > > code runs without problem. > > > > > > > > > Any idea where could be the issue? > > > > > > Thank you very much for your help > > > > > > Best Regards > > > > > > Emilio Ciuffoli > > > > > > > > > > > > > > > > __________________________________________________________________________ > > You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info > </emilio_at_impcas.ac.cn></paola.sala_at_mi.infn.it>
This archive was generated by hypermail 2.3.0 : Fri Jul 22 2022 - 19:17:35 CEST