From: Giuseppe Battistoni (Giuseppe.Battistoni@mi.infn.it)
Date: Tue Feb 06 2007 - 10:34:03 CET
Hello,
I gave a look at your input file.
You have a thick target and you are requiring the yield of particles
(neutrons) emerging from inelastic interactions. In this way you are
counting not only neutrons produced by muons, but also those from other
secondary interactions.
I cannot say if this is a meaningful answer to your question, since it is
not clear to me if 27.1 x10-5 neutrons per muon per
g/cm2 includes reinteractions or not.
If you need to isolate neutron production by muons only there could be
some solution. For instance you could
use a thin target and at the same time you can bias the inelastic cross
sections of primary muons (see LAM-BIAS with INEPRI sdum)
Alternatively you can select particles of interest by means of a fluscw.f
routine (activated by USERWEIG). This requires also the use of
stuprf.f to tag with the the available user flags neutrons produced by
muons. This could not be an easy implementation...
Giuseppe Battistoni
On Fri, 2 Feb 2007, Lindley Winslow wrote:
> Hello all,
>
> I am trying to reproduce neutron yield results in C10H22 with the
> attached setup and scoring using USRYIELD. The paper (PhysRevD.
> 64.013012 this may be familiar to some) predicts 27.1 x10-5 neutrons
> per muon per g/cm2 where I am getting 37.2x10-5 neutrons per muon per
> g/cm2 at Emu=285.0GeV. My numbers from integrating dN/dOmega/dE and
> dN/dE/dOmega are consistent so I was thinking I had another issue.
>
> I was wondering if by scoring neutrons emerging from inelastic
> interactions I will be double counting neutrons from neutron
> scattering or neutron spallation.
>
> Thank you,
> Lindley
>
>
This archive was generated by hypermail 2.1.6 : Tue Feb 06 2007 - 10:43:19 CET