I guess,
thgere are not too many Fortran guys around, so I took upon my self to
change the gbrahms to do just this.
It only took about 1/2 hour so it was not too bad.
The control/spot command is implemented xommandwise exaclt as Christian
suggested i.e.
> user/control/spot <meanx> <varx> <meany> <vary> <meanz> <varz>
> <absolut maxx> <absolut maxy> <absolut maxz>
The limits is ONLY imposed for gaussian distributions, and the condition is
abs(mean+randomgauss*var) .lt. <Maxdist>.
Note it is possible but not likelt to create near infinite loops e.g.
cont/spot 0 0 0 0 400 0.1 0 0 10
will put the center at 4 meter and a width of 1 mm and req the vert to be
within +-10 cm.
Yuo calculate the number of thros to get there!
Modified code
guxint.F
guinti.cdf
gxuser.F
gcuser.inc
> > u/c/spot 0. 0. 0. 0. -0.95 77.16
>
> There seem to be some problem with the the sign of the mean also,
> since the vertex distribution that Erik get out has a mean of 0.95.
>
By the way there is NO sign problem in gbrahms - it is easy to see in the
code secondly one can try e.g. by
the following test case
geodef
geofin
cards ukin 1
kine 1 1 1 90 90 0 0 8
cont/spot 0 0 0 0 5 0.1 0 0 0
trig;cont/print vert 1
gives you Vz distributed around 5.
If you get otherwise, please check the analysis code.
Enjoy
Flemming
------------------------------------------------------
Flemming Videbaek
Physics Department
Brookhaven National Laboratory
tlf: 631-344-4106
fax 631-344-1334
e-mail: videbaek@bnl.gov
This archive was generated by hypermail 2b29 : Fri Jun 01 2001 - 17:17:01 EDT