RE: tracking or geometry problem?

From: Vasilis Vlachoudis <>
Date: Tue, 31 Aug 2010 12:11:13 +0200

Dear Makis,

You can always break the ROT-DEFI cards into two or more cards
One with only the translation and a subsequent one with only the
rotation. Of course all of them should have the same SDUM or WHAT(1)
You have to remember that the order of the transformations is important.


-----Original Message-----
[] On Behalf Of Francesco Cerutti
Sent: Monday, 30 August 2010 10:25
To: Chrysostomos Valderanis
Subject: RE: tracking or geometry problem?

Dear Makis,

indeed you did not do anything stupid. The unfortunate problem is that
also free format cards turn out to be (unintentionally) limited to 80
characters. So your ROT-DEFI cards are not interpreted correctly.
Then, for the time being, the increaded precision allowed by the free
format has to be reserved for the WHATs actually needing it (especially
the rotation angle), without breaking the 80 character card limit.



Francesco Cerutti
CH-1211 Geneva 23
tel. ++41 22 7678962
fax ++41 22 7668854

On Sat, 28 Aug 2010, Chrysostomos Valderanis wrote:

> Dear Francesco,
> =20
> Thank you very much for your answer! You remark about the increased
> through the use of the FREE card for the ROT-DEFI was really
> Unfortunately I can not implement it correctly. In order to write the
> in free format I filled all the values which where asumed in fixed
> input. Flair seems to interpret correctly the cards. But whem I look
at the
> Fluka output file, under the coordinate transformations section, I see
> different transformations. The net outcome is that in the free format
> version of my input nothing works due to geometry errors.
> I guess I did something stupid in the way I write the ROT-DEFI cards
in free
> format, but I cannot locate it.
> =20
> Any hint?
> Makis
> =20
> =20
> From: Francesco Cerutti
> Sent: Thu 26/08/2010 18:55
> To: Chrysostomos Valderanis
> Cc:
> Subject: Re: tracking or geometry problem?

> Dear Makis,
> > 1. Since only tracking is required in my code and no fluctuations
> (ie
> > pencil beam), I though that the boundary crossings of two events are
> > to be *exactly* the same. As you can see from the logfile this is
not the
> > case. What am I missing?
> actually the initial step is randomized, in order to avoid unnatural
> effects for transport in a material. We are now going to remove this
> randomization for transport in vacuum. Anyway, you can already get
> identical trajectories by setting the tracking precision through the
> MGNFIELD card (WHAT(1-2-3)), and the STEPSIZE card too. See
> ansport1009.pdf
> (from slide 29 onwards)
> > 2.The proton is not crossing all the boundaries I was hoping to. As
> can
> > see from the logfile there is a region succesion that is aborted at
> > 265. the proton is exiting a region I was hoping to enter first. The
> > result of this is that the proton is not seeing the full magnetic
> and
> > gets a wrong direction. I checked with the geometry and didn't find
> anything
> > unusual. Is there anything else to check?
> There are in fact some problems in your geometry. In particular, the
> transformations do not bring the Q1En container exactly upon the
> prototype. In order to do that, the z-coordinate of the Q1En vertex
> be 3993.96549586196, and not 3994.012495861956, being an accuracy
> than 1e-5 totally unacceptable (indeed you need much better). With
> correction, you will see the same problem occurring later along your
> magnet string, for analogous reasons.
> Otherwise, note that you can increase the accuracy of the
> definitions by using the free format (only) for the ROT-DEFI cards
(for this
> purpose, insert a FREE card and a FIXED card just before and after -
> respectively - the ROT-DEFI card list).
> Ciao
> Francesco
Received on Tue Aug 31 2010 - 14:22:46 CEST

This archive was generated by hypermail 2.2.0 : Tue Aug 31 2010 - 14:23:16 CEST