Dear Ogawa,
You have defined HEAVYION as primary particle. Thus IJBEAM=-2 and the
source will never enter the branch where you set
ILOFLK (NPFLKA) = jpart(sample)
Furthermore, you have modified in the default source the index of AM in
your source:
PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
& + TWOTWO * AM (ILOFLK(NPFLKA)) ) )
default:
PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
& + TWOTWO * AM (IONID) ) )
In case of a HEAVYION primary ILOFLK(NPFLKA) contains the identity of the
ion, in your case carbon (601230) and IONID=-2. Thus, using ILOFLK(NPFLKA)
to index AM causes a crash, while IONID would not.
I suggest to set your beam particle to a hadron of the highest energy
of any particle in your data file for your simulation.
Best regards
Stefan
On Mon, 11 Jan 2010, Ogawa Tatsuhiko wrote:
> Dear FLUKA experst,
>
> Please find the attached input file. This is the input file I used actually. I
> also attached the source.f file again for convenience.
>
> I would appreciate if you could point out what is wrong in this input.
>
> Sincerely,
> Ogawa
>
>
-- ___________________________________ Stefan Roesler CERN, DGS-RP CH-1211 Geneva 23 Switzerland Phone: +41-22-7679891 Fax: +41-22-7669639 E-mail: Stefan.Roesler_at_cern.chReceived on Mon Jan 11 2010 - 11:12:22 CET
This archive was generated by hypermail 2.2.0 : Mon Jan 11 2010 - 11:12:23 CET