Re: residual nucleus infomation

From: Shirley Li <shirleyli.phys_at_gmail.com>
Date: Mon, 28 Jan 2013 16:40:02 -0500

Hi Mary,

I just noticed that it seems like all the residual nuclei in USDRAW from
reactions ICODE = 101 or 300 are scored in USRRNC as well. But, there are
residual nuclei from reactions ICODE = 102, 103, 208 etc, and they are not
scored in USRRNC. Are there any particular reasons for this?

Thank you very much!

Shirley


On Mon, Jan 28, 2013 at 2:53 PM, Mary Chin <mary.chin_at_cern.ch> wrote:

> Hi Shirley,
>
> The two will not match when multigroup transport is in operation because
> it is non-analog. Outside multigroup transport, the correspondence should
> be there, but not necessarily at a call-by-call level. USDRAW is called per
> collision, whereas USRRNC is called when a resnuc stops. Following a
> collision a product may be transported, and will have its USRRNC call only
> when it finally stops.
>
> :) mary
>
>
>
> On Sun, 27 Jan 2013, Shirley Li wrote:
>
> Hi Mary,
>> Thanks! I added looking over NP and NPHEAV to get all the products in
>> mgdraw. For isotopes like heavier
>> than He, I count the number from either (ICRES, IBRES) or (ICHEAV,
>> IBHEAV) from USDRAW. Now most of the
>> isotopes number match between the two routines! :) However, some of the
>> isotopes like 9Be, 11B numbers
>> still mismatch by several percent. And, I get much more 13C from mgdraw
>> than from usrrnc.
>> I attached my inputs. Did I miss something in the loop in mgdraw?
>>
>> Thank you so much!
>>
>> Shirley
>>
>> On Sat, Jan 26, 2013 at 5:34 PM, Mary Chin <mary.chin_at_cern.ch> wrote:
>> Hi Shirley,
>>
>> If we take the following collision as an example:
>> p + 12C -> 7Be + n + p + a
>>
>> 7Be, 1H and 4He are going to get scored as RESNUCLE i.e. we get
>> three out of a single
>> collision. USDRAW is called once. USRRNC is called three times.
>>
>> If we try to print out from a USDRAW call, without any looping
>> (over an array or something),
>> we can't possibly get more than one output.
>>
>> To do the cross-check you intend to, I would loop through NP and
>> NPHEAV in USDRAW to get the
>> full suite of products following a collision. (I have done that
>> before, they matched all
>> right!)
>>
>> :) mary
>>
>>
>> On Sat, 26 Jan 2013, Shirley Li wrote:
>>
>> Hi Mary,
>> Thank you very much for your explanation!
>>
>> So, if I check ICRES from mgdraw.f and IZ from usrrnc, should
>> I expect exactly
>> the same number? I'm
>> running 1000 muons through oil (CH2) to get isotope
>> production. I print out
>> ICRES, IBRES and WTRACK
>> from mgdraw.f and IZ, IA, WEE from usrrnc.f. However, the
>> numbers don't always
>> match. E.g., I got 20
>> 7Be from mgdraw.f and 33 7Be from usrrnc.f. I understand that
>> I can get isotope
>> production from
>> RESNUCLEi scoring, but I just want to make sure that
>> everything is consistent.
>> Thanks a lot!
>>
>> Shirley
>>
>>
>> On Fri, Jan 25, 2013 at 5:02 PM, Mary Chin <mary.chin_at_cern.ch>
>> wrote:
>> Dear Shirley,
>>
>> ICRES is a variable in RESNUC, which is an include file.
>> IZ is an argument for USRRNC, which is a subroutine.
>>
>> ICRES will remain undefined unless the user routine
>> contains the line:
>> INCLUDE ‘(RESNUC)’
>> Incidently, in the usrrnc.f template distributed, this
>> line is absent.
>>
>> If IZ is already available as an argument in usrrnc.f,
>> why would one
>> attempt/wish to use
>> ICRES instead? One would use ICRES, for example, in
>> mgdraw.f, where the
>> atomic number is
>> not readily available as an argument.
>>
>> :) mary
>>
>> On Thu, 24 Jan 2013, Shirley Li wrote:
>>
>> Dear FLUKA users,
>> I'm not sure what's the difference between ICRES
>> (residual nucleus
>> atomic
>> number) in RESNUC and IZ
>> (also residual nucleus atomic number) in
>> usrrnc.f. Could anyone
>> explain this
>> to me?
>>
>> Thank you very much!
>>
>> Best,
>> Shirley
>>
>>
>>
>>
>>
>>
Received on Tue Jan 29 2013 - 09:49:50 CET