From: Giuseppe Battistoni (Giuseppe.Battistoni@mi.infn.it)
Date: Wed Sep 10 2003 - 10:50:52 CEST
Date: Wed, 10 Sep 2003 10:47:16 +0200
From: Giuseppe Battistoni <giuseppe.battistoni@mi.infn.it>
Organization: INFN
User-Agent: Mozilla/5.0 (X11; U; Linux i686; en-US; rv:1.0.2) Gecko/20030708
X-Accept-Language: en-us, en
MIME-Version: 1.0
To: Christopher Knievel <Christopher.Knievel@lightsource.ca>
CC: "FLUKA Discuss (E-mail)" <fluka-discuss@mi.infn.it>
Subject: Re: Unknown error in geometry
References: <9DABE63F0EBE544C9C2313AA6AE117E2476B56@srv-mail.cls>
Content-Type: text/plain; charset=us-ascii; format=flowed
Content-Transfer-Encoding: 7bit
Sorry, we realize just now that a question of yours was still pending.
Just in case you never succeeded in solving your problem: in situations
like the one you are describing,
the use of geometry debugger is required. If you never did that, look at
the GEOEND directive in the
user reference guide (fluka2002.manual or the web page)
Best regards
Giuseppe Battistoni
Christopher Knievel wrote:
>Dear Alberto,
>
>Thank you very much for your reply. I understand what my problem was after you explained it to me and I went and made the changes to the input file. However it then gave me an error about tracking particles in Region 2. I also realized, after you pointed it out, about the problem of the undefined region left of plane 3 and right of plane 4 inside cylinder 1. I just went and redid my geometry now and made it a little different. Instead of enclosing the cylindrical beam stop in a cylindrical vacuum shell, I tried to enclose it in a vacuum sphere of radius 100 cm. Once again I received almost the exact same error message when I tried this approach. The error reads:
>
>Geofar: Particle in region 3 (cell # 0) in position -5.338737686E-01 3.151553314E+00 4.610356519E+01
> is now causing trouble, requesting a step of 1.808186627E+01 cm
> to direction -5.339776322E-01 5.616004021E-01 6.320386671E-01, error count: 0
> R2: 3.196452642E+00 R3: 4.621424058E+01 cm
> X*U (2D): 2.054990259E+00 X*U (3D): 3.119422615E+01 cm
> X*UOLD(2D): 2.108308487E+00 X*UOLD(3D): -2.700856933E-01 cm
> Kloop: 39, Irsave: 3, Irsav2: 3, error code: -33 Nfrom: 5000
> old direction 6.291678888E-01 7.755555675E-01 -5.158807504E-02, lagain, lstnew, lsense, lsnsct F F F F
> Particle index 7 total energy 7.251255011E-04 GeV Nsurf 0
> Try again to establish the current region moving the particle of a 1.000000000E-06 long step
> We succeeded in saving the particle: current region is n. 3 (cell # 0)
> *** Emfgeo: Ustep, Dnear, Nreg, Newreg, Iq(np), X, Y, Z, Ddnear 1.E-06
> 2.73727372 3 3 0 -0.533873769 3.15155331 46.1035652 0.
>
>Once again I will attach the input file which I currently using. Now the error says that the problem occurs in region 3. But I am pretty sure that I have described that region correctly.
>* Region 3 (beam stop)
> 3 5 +1 -2 +3
>This would correspond to points inside cylinder 1, to the right of plane 2 and to the left of plane 3. Therefore a point can be simulateneously >0 and <61. Once again I am at a loss and any help would be appreciated.
>
>Chris Knievel
>
>
>
>
>
>From: Alberto Fasso' [mailto:Alberto.Fasso@cern.ch]
>Sent: Tuesday, April 01, 2003 5:15 PM
>To: FLUKA Discuss (E-mail)
>Subject: Re: Unknown error in geometry
>
>
>
>Indeed, your geometry description is incorrect.
>I will try to list a few things which are obviously wrong or inconsistent, but the main message is that a geometry such as the one you are trying to describe cannot be done without the use of the OR operator.
>
>1) Region 1 (blackhole):
> 7 -2 4 -6
> that means "the space inside sphere 7, which is outside cylinder 2, on the
> left of plane 4 and on the right of plane 6"
> No point can be at the same time on the left of plane 4 and on the right of
> plane 6: it is the same as saying that its Z coordinate must be both < -5 and
> > 66. What you probably wanted was:
> OR 7 -2OR 7 4OR 7 -6
>
>2) Region 2 (vacuum at the beginning):
> 2 -1 3 -4 6 -5
> This would mean: "the space inside the wall of a cylindrical shell contained
> between cylinder 2 and cylinder 1, left of plane 3, right of plane 4, left
> of plane 6 and right of plane 5"
> Again, a point cannot be at the same time be left of plane 3 and right of
> plane 5. And it is not clear why you want vacuum only between 2 and 1 and not
> also inside 1. In any case, the space inside cyl. 1, left of plane 3 and
> right of plane 4 remains undefined. What you probably wanted was:
> OR 2 3 -4OR 2 6 -5OR 2 -1 6 -4
>
> Alberto
>
>
>
>On Tue, 1 Apr 2003, Christopher Knievel wrote:
>
>
>
>>Dear FLUKA users and authors:
>>
>>I am having a problem when I run my input with fluka. I am trying to
>>simulate an electron beam interacting with a cylindrical aluminum beam
>>stop. The error which I receive in the .err file is as follows:
>>
>> Geofar: Particle in region 3 (cell # 0) in position 5.730835694E+00 9.686117683E+00 3.982049539E+01
>> is now causing trouble, requesting a step of 9.549625271E+00 cm
>> to direction 1.929431495E-01 7.823632011E-01 5.921830482E-01, error count: 0
>> R2: 1.125448149E+01 R3: 4.138037224E+01 cm
>> X*U (2D): 8.683787525E+00 X*U (3D): 3.226480986E+01 cm
>> X*UOLD(2D): -2.805983023E+00 X*UOLD(3D): 2.332732094E+01 cm
>> Kloop: 49, Irsave: 3, Irsav2: 3, error code: -33 Nfrom: 5000
>> old direction 4.859413563E-01 -5.772006160E-01 6.562777210E-01, lagain, lstnew, lsense, lsnsct F F F F
>> Particle index 7 total energy 7.510855076E-04 GeV Nsurf 0
>> Try again to establish the current region moving the particle of a 1.000000000E-06 long step
>> We succeeded in saving the particle: current region is n. 3 (cell # 0)
>>
>>This error continues for many more cycles. In the .out file the error
>>reads as:
>>
>>0EXIT BEING CALLED FROM G1, NEXT REGION NOT FOUND.
>> IR= 3 XB= 0.57308359E+01 0.96861185E+01 0.39820496E+02 WB= 0.19294315E+00 0.78236320E+00 0.59218305E+00 DIST= 0.48005062E+01
>> IRLTGG= 0 SB= 0.57308359E+01 0.96861185E+01 0.39820496E+02 UB= 0.19294315E+00 0.78236320E+00 0.59218305E+00
>> PR 10 0.57308359E+01 0.96861185E+01 0.39820496E+02 0.19294315E+00 0.78236320E+00 0.59218305E+00 3
>> 3 3 1 -1 1 0 1 54 54 10
>> 93 101 1 93 95 4 0 105 10 10
>>-0.31547868E+02 0.48005062E+01 0.48005062E+01 0.48005062E+01
>>
>>I will attach my input file which I used. I understand that the error
>>means that the program is having problems tracking the particle from
>>point (x,y,z) to point (x',y',z') defined by the directional vectors
>>given above, I just do not know what to do to correct the problem. As
>>one can see from my input, the beam dump should be 61 cm long and it
>>is enclosed in a cylindrical vacuum shell that is 5 cm thick. At least this is what I tried to simulate.
>>Transport should also start at the front of the aluminum cylinder, which is
>>the coordinate (0,0,0). If anyone could assist me with my problem, it would
>>be greatly appreciated.
>>
>>Sincerely,
>>
>>Christopher Knievel
>>Health, Safety and Environment Assistant
>>Canadian Light Source
>>
>> <<dice.inp.txt>>
>>
>>
>>
>
>
>
>------------------------------------------------------------------------
>
>TITLE
> Test of Interaction with an Aluminum Beam Dump
>DEFAULTS NEW-DEFA
>*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
>BEAM -2.9 ELECTRON
>BEAMPOS 0.0 0.0 0.0
>GEOBEGIN 0.01 COMBINAT
> cylindrical beam stop
> ZCC 1 0.0 0.0 15.0
> XYP 2 0.0
> XYP 3 61.0
> SPH 4 0.0 0.0 30.0 100.0
> SPH 5 0.0 0.0 30.0 300.0
> END
>* black hole
> 1 5 +5 -4
>* vacuum surrounding
> 2 5OR +4 -1OR +4 +2OR +4 -3
>* beam stop
> 3 5 +1 -2 +3
> END
>GEOEND
>*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
>MATERIAL 13.0 26.982 2.70 10.0 ALUMINUM
>* Aluminum Beam Stop
>ASSIGNMAT 10.0 3.0
>* External Black Hole
>ASSIGNMAT 1.0 1.0
>* Vacuum
>ASSIGNMAT 2.0 2.0
>*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
>STEPSIZE 1.0E-09 0.000001 1.0 3.0
>SCORE 208.0 211.0
>OUTLEVEL 1.0 7.0
>* **** Scoring ****
>PHYSICS 2.0 EVAPORAT
>RESNUCLEI 3.0 11.0 3.0 1.0 DUMMY
>*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
>*USRTRACK 1.0 201.0 11.0 3.0 10.0 1.0 Trdetect
>*USRTRACK 2.9 0.01 &
>*USRBDX 1.0 201.0 11.0 3.0 2.0 1.0 bddetect
>*USRBDX 2.9 0.01 10.0 10.0 &
>*RANDOMIZE 1.
>*23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
>START 150.0 99999999. 0.0 1.0
>STOP
>
>
This archive was generated by hypermail 2.1.6 : Wed Sep 10 2003 - 10:51:04 CEST