Dear Franz,
this is cause by beam is traveling exactly on the surface of the region
defined by the RPP.
The code tolerate a certain number of those cases, but in your case those
are 25% of your primaries and FLUKA aborts the run.
One could tune (relax a bit) the accuracy on the boundary identification in
this case using WHAT(2) of the GEOBEGIN card.
Best,
Vittorio
On Fri, May 6, 2016 at 1:53 PM, Franz Siegfried Englbrecht <
Franz.Englbrecht_at_physik.uni-muenchen.de> wrote:
>
> Dear colleagues,
>
> I encountered a problem in FLUKA when simulating an experiment of a proton
> beam interacting with a solid target in vacuum, which can be reduced to a
> simple simulation to reproduce the error as follows:
>
> A proton beam is shot from the source position at (0,0,0) in +z direction
> towards a solid RPP target in vacuum sitting at z=10cm with one face being
> parallel to an XY plane. The Extensions are Xmin = 0, Xmax = 5, YMin =
> -2.5, Ymax = 2.5, Zmin = 10, Zmax = 15.
>
> It seems that FLUKA has problems, when protons are hitting the face of the
> RPP parallel to one face, since the .err file reports:
>
> Geofar: Particle in region 3 (cell # 0) in position
> 0.000000000E+00 0.000000000E+00 1.000010264E+01
> is now causing trouble, requesting a step of 5.118263514E-04 cm
> to direction -9.417213493E-03 8.622442818E-04 9.999552853E-01, error
> count: 0
> R2: 0.000000000E+00 R3: 1.000010264E+01 cm
> X*U (2D): 0.000000000E+00 X*U (3D): 9.999655489E+00 cm
> X*UOLD(2D): -0.000000000E+00 X*UOLD(3D): 1.000006331E+01 cm
> Kloop: 88, Irsave: 3, Irsav2: 3, error code: -3
> Nfrom: 21
> old direction -6.589281883E-05 -2.803760406E-03 9.999960673E-01, lagain,
> lstnew, lsense, lsnsct F F F F
> Particle index 1 total energy 9.482567258E-01 GeV Nsurf 0
> We succeeded in saving the particle: current region is n. 2 (cell #
> 0)
>
> Also the .log file and .out file report:
>
> Abort called from GEOFAR reason TOO MANY ERRORS IN GEOMETRY: STOP Run
> stopped!
> STOP TOO MANY ERRORS IN GEOMETRY: STOP
>
>
> Please find my inputfile for this testcase attached.
>
> Sincerely,
>
> Franz Englbrecht
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue May 10 2016 - 00:09:53 CEST