Re: fluka + pythia

Dear Ercan,

the error is due to the fact that some of the particles you read from the
file have an energy higher than the 5 GeV you have indicated by the BEAM
command in your input file. That value is used to set the higher limit of
cross section tabulations, so it should be higher than the energy of any
particle you read with your source.

In the manual, for each particle that can be transported with FLUKA, the
corresponding PDG number is also reported. PDG numbers non reported there
represent particles with a half-life too short to be transported.
You should run Pythia letting those particles decay, so that the file
would contain only particles that FLUKA can transport.

As an alternative solution, not very elegant (but perhaps sufficient for
your purposes?) you can replace each of those particles by a pion with
the same charge and the same total energy (mass+kinetic). At least, this
trick will conserve energy. But of course, it will ignore all particles lighter
than the pion that can be produced in the decay (electrons, photons,
neutrinos). So, although energy is conserved, neutrino energy which would
escape from the system will remain in the system. I would much prefer
having Pythia doing the decay.

Alberto

On Sat, 18 Jun 2011, Ercan Pilicer wrote:

> dear Fluka users
>
> i want to initiate a simulation with the primaries generated by Pythia.
> in the attachment there are input and source files.
>
> i got the error attached, but i am not sure whether source.f is
> correct. what might be the reason?
>
> another question is, since Fluka particle codes are differing from PDG
> codes, what number should be used
> for the particles not listed in Fluka particle names (Table5.1)?
>
> or mass number given in Pyhtia can be used to evaluate kinetic energy
> without looking for the particle number?
>
> thanks in advance
> e.
>
>
> =====================================
>
> Ercan Pilicer
> Uludag University
> High Energy Physics Department
> 16059 Bursa, TURKEY
>
> =====================================
Received on Mon Jun 20 2011 - 09:35:56 CEST