From: Giuseppe Battistoni (Giuseppe.Battistoni@mi.infn.it)
Date: Fri Aug 06 2004 - 11:51:43 CEST
We always give cosines in the source and it works (with a difference, see
below). Therefore there must be an error somewhere
in your source.
Are you sure that you are normalizing the direction cosines
properly (i.e. within the rounding of the machine)?
FLUKA makes a check on that.
You should define
TZ(LSTACK) = SQRT( ONEONE - TX(LSTACK)**2 - TY(LSTACK)**2)
and then adjust the sign of TZ according to your needs.
Giuseppe Battistoni
On Thu, 5 Aug 2004, LEE, KERRY T. (JSC-SR) (UHCL) wrote:
> Date: Thu, 5 Aug 2004 14:21:13 -0500
> From: "LEE, KERRY T. (JSC-SR) (UHCL)" <kerry.t.lee1@jsc.nasa.gov>
> To: "'fluka-discuss@fluka.org'" <fluka-discuss@fluka.org>
> Subject: USER SOURCE
>
>
> Hi,
>
> I am using the user SOURCE routine and I am trying to set the direction
> cosines of the beam in my source.f routine. I have successfully run
> changing both cosX and cosY, but when I try and change cosZ nothing happens.
> When I set coxZ to 1 or -1 the code gets hung in an infinite loop, and the
> error output shows that cosX=0, cosY=0, and cosZ=0, therefore, it thinks
> there is no movement.
>
> I use the following lines to set the direction cosines
>
> TX (LSTACK) = cosx_
> TY (LSTACK) = cosy_
> TZ (LSTACK) = cosz_
>
> Is there something additional that must be done to change cosZ properly?
>
> When I use the BEAMPOS card without SOURCE everything works as expected.
>
> Thanks
> Kerry Lee
>
This archive was generated by hypermail 2.1.6 : Fri Aug 06 2004 - 12:09:43 CEST