Re: Use of Resnucle and Part-thr

> Unexpectedly scored particles are neutron, proton, alpha.

Running your input with 2 million primaries, I see neutrons, protons, alphas,
and pions only on the front face of the target (as expected due to interactions
in air) and also, in quite reduced amount, at the following boundary, 1.5 mm
downstream. I guess that the reason for the latter ones is the fact that, due
to your sharp condition in usrmed (Zz .ge. ZERZER) and to the - not zero! -
accuracy in boundary identification, they are not killed on the front face and
do a further step

Cheers

Francesco

On Thu, 11 Feb 2010, t_ogawa wrote:

> Dear Francesco,
> CC: FLUKA discuss,
>
> I am sorry for poor explanation.
> It's true that I see secondaries on the front surface and they are expected.
> The problem is I saw particles in the other tallies located inside iron
> target. Could you tell me why secondary particles are scored inside the
> target?
>
> Sincerely,
> Ogawa
>
> Francesco Cerutti wrote:
>>
>> Dear Ogawa,
>>
>> the secondaries you see on the front face of your target (according to your
>> scoring I see them only there) are created in air, where you do not kill
>> them
>>
>> Ciao
>>
>> Francesco
>>
>> **************************************************
>> Francesco Cerutti
>> CERN-EN/STI
>> CH-1211 Geneva 23
>> Switzerland
>> tel. ++41 22 7678962
>> fax ++41 22 7668854
>>
>> On Wed, 10 Feb 2010, t_ogawa wrote:
>>
>>> Dear Francesco,
>>> CC FLUKA-discuss,
>>>
>>> Thank you very mush for your advice.
>>> I would like one more piece of advice to complete this calculation.
>>>
>>> I ran FLUKA with executable including usrmed.f (attached) intending to
>>> transport only primary particles inside iron but I still saw some
>>> particles without being killed.
>>> Could you help me to find the error in the subroutine or .inp file to kill
>>> all particles but original beam?
>>>
>>> Sincerely,
>>> Ogawa
>>>
>>> Francesco Cerutti wrote:
>>>>
>>>> Dear Ogawa,
>>>>
>>>> indeed you uncovered a bug since the *proton* threshold in fact
>>>> interferes even with ion transport (though it was not intended to). If
>>>> the ion kinetic energy per nucleon is below the proton threshold, only
>>>> dE/dx and possibly delta production apply (and not nuclear interactions).
>>>> Anyway, as an available solution, you can kill every particle but ions by
>>>> resetting to zero the weight WEE=ZERZER in the usrmed user routine as a
>>>> function of the particle identity IJ. Note that heavy ions, in addition
>>>> to IJ=-2, carry also IJ values less than -6, in case they are reaction
>>>> products. As you know, light ions (deuterons, tritons, 3He's and alphas)
>>>> cover the range between -3 and -6. To activate usrmed, you need a
>>>> MAT-PROP card with SDUM=USERDIRE applying to materials.
>>>>
>>>> Ciao
>>>>
>>>> Francesco
>>>>
>>>> **************************************************
>>>> Francesco Cerutti
>>>> CERN-EN/STI
>>>> CH-1211 Geneva 23
>>>> Switzerland
>>>> tel. ++41 22 7678962
>>>> fax ++41 22 7668854
>>>>
>>>> On Mon, 1 Feb 2010, t_ogawa wrote:
>>>>
>>>>> Dear FLUKA experts,
>>>>>
>>>>> Could someone tell me how to score residual nuclides induced by ions
>>>>> without contribution from secondaries.
>>>>> I am trying to score residual nuclide in a target made of iron induced
>>>>> by 400MeV/u carbon ion using resnucle card in a following way
>>>>>
>>>>> RESNUCLE 3. -30. TARCYL01
>>>>> 1.88496res1
>>>>>
>>>>> In order to disregard contribution from secondary particles, I put
>>>>>
>>>>> PART-THR -10.0 1.0 @LASTPAR
>>>>>
>>>>> Then I found no residual nuclides scored in the unit -30. When I turn
>>>>> off PART-THR card, scoring was successful.
>>>>> Many thanks in advance.
>>>>>
>>>>> Sincerely,
>>>>> Ogawa
>>>>>
>>>>> --
>>>>> **************************************************
>>>>> ATLAS Detector Operation Group
>>>>> Department of Physics
>>>>> European Organization for Nuclear Research(CERN)
>>>>> --------------------------------------------------
>>>>> Department of Nuclear Engineering and Management
>>>>> Graduate School of Engineering
>>>>> The University of Tokyo
>>>>>
>>>>> Tatsuhiko Ogawa
>>>>> E-mail: ogawa_at_n.t.u-tokyo.ac.jp
>>>>> togawa_at_cern.ch
>>>>> Tel: +81-3-5841-2905
>>>>> Facsimile : +81-3-3813-2010
>>>>> Address: 2-11-16 Engineering Bld.12th Room421
>>>>> Yayoi Bunkyo-ku, Tokyo 113-0032
>>>>> **************************************************
Received on Sat Feb 13 2010 - 19:16:26 CET