Re: Use of Resnucle and Part-thr

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
>>> **************************************************
>>>
>>>
>>
>
>
> --
> **************************************************
> 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 Thu Feb 11 2010 - 12:40:06 CET