From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Wed Aug 10 2005 - 02:58:26 CEST
You geometry is incorrectly described, as you can see from the debugger:
**** Lookdb: Geometry error found ****
**** The point: 30.0000001 0. -99.9999999 ****
**** is not contained in any region ****
**** Lookdb: Geometry error found ****
**** The point: 35.0000001 0. -99.9999999 ****
**** is not contained in any region ****
**** Lookdb: Geometry error found ****
**** The point: 40.0000001 0. -99.9999999 ****
**** is not contained in any region ****
.......................
**** Lookdb: Geometry error found ****
**** The point: 145. 0. -69.9999999 ****
**** is not contained in any region ****
**** Lookdb: Geometry error found ****
**** The point: 150. 0. -69.9999999 ****
**** is not contained in any region ****
*** Error limit reached: geometry debugging aborted ***
Remember that + means "is inside" and - means "is outside".
Your region description should be:
* reg.....or.....or.....or.....or.....or.....or.....or.....or.....or.....
* The blackhole
1 5 +1 -2
* The vacuum (3 pieces: upstream, around and downstream of the target)
2 5OR +2 +4OR +2 -3OR +2 -5
* The concrete
3 5 +3 +5 -4
END
And please, don't play with the geometry accuracy parameter.
Leave it = 0.0 (the default): it is for very special cases,
certainly not for a beginner!
Alberto
On Tue, 9 Aug 2005, Khnkanos N. Sanosyan wrote:
> Dear FLUKA developers and experts,
>
> Thank you very much for your help, remarks and advices!
> I passed the problem mentioned in my previous message thanks to your
> advices. As it could be expected now I have a problem of geometry as following:
> 1. In .err file:
> .....
> Geofar: Particle in region 3 (cell # 0) in position
> 2.741794896E-01 -6.158208031E+00 0.000000000E+00
> is now causing trouble, requesting a step of 1.000000000E+08 cm
> to direction -3.143536261E-01 -2.319948821E-01 -9.205216850E-01, error
> count: 0
> R2: 6.164308602E+00 R3: 6.164308602E+00 cm
> X*U (2D): 1.342483429E+00 X*U (3D): 1.342483429E+00 cm
> X*UOLD(2D): 4.769731625E+00 X*UOLD(3D): 4.769731625E+00 cm
> Kloop: 36, Irsave: 4, Irsav2: 3, error code: -33
> Nfrom: 5000
> old direction 5.996473208E-01 -7.478345334E-01 2.848975277E-01,
> lagain, lstnew, lsense, lsnsct F F F F
> Particle index 7 total energy 1.572643271E-04 GeV Nsurf 4
> Try again to establish the current region moving the particle of a
> 1.000000000E-03 long step
> We succeeded in saving the particle: current region is n. 3 (cell #
> 0)
> Geofar: Particle in region 3 (cell # 0) in position
> 2.738651360E-01 -6.158440025E+00 -9.205216850E-04
> is now causing trouble, requesting a step of 1.000000000E+08 cm
> ...
> 2. in .log file:
> ENTER IN HIST INIT.
> HROPEN: ISTAT= 0
> Sono arrivato in BXDRAW 29 3 4 3 3.00000004 0.00051099906
> Sono arrivato in BXDRAW 29 4 3 7 0.000157264327 0.
>
> 3. in .out file:
> Exit being called from G1, next region not found.
>
> My input file is attached.
> Please, have you got any suggestion for me to leave geometry uncertainties
> in my task?
> Thank you for your help.
> Khnko
>
> Khnkanos N. SANOSYAN
>
> ASLS (Armenian Synchrotron Light Source)
> CANDLE (Center for the Advancement of Atjaryan St. 31, Avan,
> Natural Discoveries using Light Emission) Yerevan 375040 / ARMENIA
>
> Tel.: (+374-1) 612569 Private: (+374-1) 645865
> Fax: (+374-1) 629806
> Mail: khnko@asls.candle.am
> http: //www.candle.am/
>
>
-- Alberto Fassò SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025 Phone: (1 650) 926 4762 Fax: (1 650) 926 3569 fasso@slac.stanford.edu
This archive was generated by hypermail 2.1.6 : Wed Aug 10 2005 - 12:39:02 CEST