From: Francesco Cerutti (Francesco.Cerutti@cern.ch)
Date: Wed Nov 01 2006 - 10:18:07 CET
Hi Joseph
I think that in fact your problem is not concerning your source, but your
geometry instead. If you try to debug it just around the point where you
place your primary particle, using for instance
GEOEND 0.0 0.0 10.0 -0.5 -0.1 5.0DEBUG
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+..
(instead of GEOEND alone -> see the manual)
you can immediately see that there is a huge amount of points not
belonging to any region.
The reason is that the operator "+" means intersection and not union. So
the region R2air (where I guess that your primary particle was supposed to
start) is in fact empty, due to the intersection of bodies not overlapping
at all.
(Moreover, following the recent suggestion of Alberto, in the
source.f statement
upart=sqrt(1.0-vpart*vpart-upart*upart)
replace 1.0 by 1.0d0 or oneone (double precision required)).
Greetings
Francesco
**************************************************
Francesco Cerutti
CERN-AB
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
This archive was generated by hypermail 2.1.6 : Thu Nov 02 2006 - 10:10:56 CET