Re: LATTICE and ROT-DEFIni

From: Francesco Cerutti <Francesco.Cerutti@cern.ch>
Date: Mon Apr 28 2008 - 11:20:14 CEST

Dear Markus,

in the LATTICE cards you have to associate CONTAINER regions (WHAT(1-3))
to transformations (SDUM) (and to lattice indexes (WHAT(4-6)), used
internally by the code and fully uncorrelated with the rest of the
input. For them in your case a logical choice would be 1.0 for the magnet
2 and 2.0 for the magnet 3, since these magnets - regions 5 and 6 -
represent the two lattice cells of your geometry. This is indeed the
default choice gotten by the cards that I put below. Keep in mind that
your magnet 1 is NOT a lattice cell, but the prototype itself, made by
regions 3 and 4 which prevent from defining a further container region
there).
Here the correct cards:
LATTICE 5.0 ROT#1
LATTICE 6.0 ROT#2
[In case the format is lost in the e-mail body, I put WHAT(1)=5.0/6.0 and
SDUM=ROT#1/ROT#2].

Since you want to use fully number-based input, you have to issue the card
GLOBAL 4.0
[WHAT(4)=4.0].

About ROT-DEFI cards, they must not lie inside the geometry description
(between GEOBEGIN and GEOEND) and they define - according to the
transformation matrices in the manual - a roto-translation carrying the
replica (x,y,z_old corresponding to the container coordinates in your
geometry) over the prototype (x,y,z_new corresponding to the prototype
coordinates in your geometry).
Thus:
GEOEND
ROT-DEFI 101.0 0.0 0.0 0.0 0.0 -100.0
ROT-DEFI 102.0 0.0 0.0 0.0 0.0 -200.0

Moreover, do not inject your beam just on a boundary (z=0), but e.g.
BEAMPOS 47.0 -0.1
[WHAT(3)=-0.1, i.e. 1mm upstream].

Cheers

Francesco

**************************************************
Francesco Cerutti
CERN-AB
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854
Received on Mon Apr 28 14:17:01 2008

This archive was generated by hypermail 2.1.8 : Mon Apr 28 2008 - 14:17:02 CEST