From: Paolo MAESTRO (paolo.maestro@pi.infn.it)
Date: Wed May 25 2005 - 22:03:17 CEST
Your guess is right.
I have modified the source routine so that the energy of the
primary particle is generated according to a power-law spectrum.
As you explained, the error occurred any time the randomly-extracted
energy overcame the BEAM card value.
I solved the problem simply writing in the card a value greater than
the upper limit of the spectrum.
Thanks for your help
Paolo
On Wed, 25 May 2005, Giuseppe Battistoni wrote:
> That error message means that you are sampling an energy greater than
> the value declared in the BEAM card.
> Are you using a source user routine? Are you using FLUKA with rqmd and
> dpmjet?
> The valued in the BEAM card is used by FLUKA as the maximum energy for
> initial tabulation of cross sections , dE/dx etc.
> If by chance you overcome that value the program goes crazy.
>
> Giuseppe
>
> On Thu, 19 May 2005, maestro wrote:
>
> > Date: Thu, 19 May 2005 18:17:53 +0200
> > From: maestro <paolo.maestro@pi.infn.it>
> > To: fluka-discuss@mi.infn.it
> > Subject: Error message explanation needed
> >
> >
> > Dear all,
> > I am simulating relativistic ions interaction in a calorimeter .
> > At a given point, my job gets into a never-ending loop
> > during which the following message is written and repeated
> > continuously in the myjob.err file:
> > ------------------------------------------------
> > Stepop, dp/dx=<0!! ij, mmat, ekin, dpdx0 -2 30 7621.6471 -0.000335210437
> > *** Stepop: Trange < 0, Ij,Pla,Pthr -2 7632.81378 0.516069975
> > Kaskad: Dool<0!!, DOOL,IJ,PLA,PLAOLD,EKIN,TUSTEP,TVSTEP,MMAT
> > -1.67605039E-19 -2
> > 7632.81378 7632.81378 7621.6471 5.E-16 5.E-16 30
> > 0 GEODEN CALLED WITH AN ECONTR VALUE OF: -0.168E-18
> > ICALL = 131 IJ = 208 MREG,NEWREG = 48 48 IEST = 0
> > --------------------------------------------------
> > Can you explain it ? I guess that the problem concerns material n. 30
> > (C2H4 i.e. light polyetylene with density = 0.035 g/cm3)
> > but I don't understand why.
> > Maybe could it be related to a wrong mat.pemf card compilation ?
> >
> > Many thanks for your help.
> > Cheers
> > Paolo
> >
> >
> >
> > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
> > Dr. Paolo Maestro
> > Physics Department & INFN
> > University of Siena
> > via Roma 56, 53100 Siena (Italy)
> > Phone +39 050 2214356
> > Fax +39 050 2214317
> >
> >
>
>
This archive was generated by hypermail 2.1.6 : Thu May 26 2005 - 09:31:16 CEST