TrackerLevel2/doc/libTrkLevel2_F77.txt

                                  ============
                                  LibTrkLevel2
                                     V0.01
                                  ============

List of F77 subroutines included in the library:

- DOTRACK
- DOTRACK2

- GRKUTA
- GRKUTA2
- GUFLD

- READB
- READBINNER
- READBOUTER

- INTER_B
- INTER_B_INNER
- INTER_B_OUTER


Usage inside F77 code:

--------------------------------------------------------------------------------
CALL READ_B(*PATH)
--------------------------------------------------------------------------------

PATH          (CHARACTER) Path of the directory that contains field-map files:
                 - External_top_map_n4_150402.rz  
                 - measure_n3_290302.rz  
                 - measure_n4_110402_corrected.rz

Read the magnetic field maps (rz files) and fill related common blocks.
===> To be called before using TRACK and INTER_B
It calls READBINNER and READBOUTER.

--------------------------------------------------------------------------------
CALL DOTRACK(*NPOINT,*ZIN,XOUT*,YOUT*,*AL_P,IFAIL*)
--------------------------------------------------------------------------------

NPOINT        (INTEGER) Num. points 
ZIN(NPOINT)   (DOUBLE)  Input z-coordinates (in cm)
XOUT(NPOINT)  (DOUBLE)  Output x-coordinates (in cm)
YOUT(NPOINT)  (DOUBLE)  Output y-coordinates (in cm)
AL_P(6)       (DOUBLE)  Track state vector
IFAIL         (INTEGER) Error flag (ok if IFLAG=0)

Calculate particle-trajectory coordinates, for the state vector AL_P, 
at z-coordinates set by the user. 
Make use of GRKUTA, routine to solve kinematic equations in a magnetic field
using Runge-Kutta method, adapted from GEANT3 library.

--------------------------------------------------------------------------------
CALL DOTRACK2(*NPOINT,*ZIN,XOUT*,YOUT*,THXOUT*,THYOUT*,TLOUT*,*AL_P,IFAIL*)
--------------------------------------------------------------------------------

NPOINT         (INTEGER) Num. points 
ZIN(NPOINT)    (DOUBLE)  Input z-coordinates (in cm)
XOUT(NPOINT)   (DOUBLE)  Output x-coordinates (in cm)
YOUT(NPOINT)   (DOUBLE)  Output y-coordinates (in cm)
THXOUT(NPOINT) (DOUBLE)  Output projected angle (in deg) on XZ plane
THYOUT(NPOINT) (DOUBLE)  Output projected angle (in deg) on YZ plane
TLOUT(NPOINT)  (DOUBLE)  Output track length (in cm) 
AL_P(6)        (DOUBLE)  Track state vector
IFAIL          (INTEGER) Error flag (ok if IFLAG=0)

Calculate particle-trajectory coordinates, for the state vector AL_P, 
at z-coordinates set by the user. 
Make use of GRKUTA2, routine to solve kinematic equations in a magnetic field
using Runge-Kutta method, adapted from GEANT3 library.
Same as DOTRACK and GRKUTA, modified in order to give as output the projected
angles and the track length.
The track length assigned to each track point i is evaluated between the i-th 
plane and the closest lower(upper) one, reference plane included, if the plane 
is above(below) the reference plane.

--------------------------------------------------------------------------------
CALL GUFLD(*V,*F)
--------------------------------------------------------------------------------

V(3)           (REAL) x,y,z spatial coordinates (in m)
F(3)           (REAL) Magnetic-field components (in kGauss)

Evaluate the three components of the magnetic field, by interpolating the
measured magnetic field map (used by GRKUTA/GRKUTA2). 
It calls INTER_B.

--------------------------------------------------------------------------------
CALL INTER_B(*X,*Y,*Z,B*)
--------------------------------------------------------------------------------

X,Y,Z         (DOUBLE) Spatial coordinates (in m)
B(3)          (DOUBLE) Magnetic-field components (in T)

Evaluate the three components of the magnetic field, by interpolating the
measured magnetic field map. 
It calls INTER_B_INNER and INTER_B_OUTER.
1	pam-fi	1.1	============
2			LibTrkLevel2
3			V0.01
4			============
5
6			List of F77 subroutines included in the library:
7
8			- DOTRACK
9			- DOTRACK2
10
11			- GRKUTA
12			- GRKUTA2
13			- GUFLD
14
15			- READB
16			- READBINNER
17			- READBOUTER
18
19			- INTER_B
20			- INTER_B_INNER
21			- INTER_B_OUTER
22
23
24			Usage inside F77 code:
25
26			--------------------------------------------------------------------------------
27			CALL READ_B(*PATH)
28			--------------------------------------------------------------------------------
29
30			PATH (CHARACTER) Path of the directory that contains field-map files:
31			- External_top_map_n4_150402.rz
32			- measure_n3_290302.rz
33			- measure_n4_110402_corrected.rz
34
35			Read the magnetic field maps (rz files) and fill related common blocks.
36			===> To be called before using TRACK and INTER_B
37			It calls READBINNER and READBOUTER.
38
39			--------------------------------------------------------------------------------
40			CALL DOTRACK(NPOINT,ZIN,XOUT,YOUT,AL_P,IFAIL)
41			--------------------------------------------------------------------------------
42
43			NPOINT (INTEGER) Num. points
44			ZIN(NPOINT) (DOUBLE) Input z-coordinates (in cm)
45			XOUT(NPOINT) (DOUBLE) Output x-coordinates (in cm)
46			YOUT(NPOINT) (DOUBLE) Output y-coordinates (in cm)
47			AL_P(6) (DOUBLE) Track state vector
48			IFAIL (INTEGER) Error flag (ok if IFLAG=0)
49
50			Calculate particle-trajectory coordinates, for the state vector AL_P,
51			at z-coordinates set by the user.
52			Make use of GRKUTA, routine to solve kinematic equations in a magnetic field
53			using Runge-Kutta method, adapted from GEANT3 library.
54
55			--------------------------------------------------------------------------------
56			CALL DOTRACK2(NPOINT,ZIN,XOUT,YOUT,THXOUT,THYOUT,TLOUT,AL_P,IFAIL*)
57			--------------------------------------------------------------------------------
58
59			NPOINT (INTEGER) Num. points
60			ZIN(NPOINT) (DOUBLE) Input z-coordinates (in cm)
61			XOUT(NPOINT) (DOUBLE) Output x-coordinates (in cm)
62			YOUT(NPOINT) (DOUBLE) Output y-coordinates (in cm)
63			THXOUT(NPOINT) (DOUBLE) Output projected angle (in deg) on XZ plane
64			THYOUT(NPOINT) (DOUBLE) Output projected angle (in deg) on YZ plane
65			TLOUT(NPOINT) (DOUBLE) Output track length (in cm)
66			AL_P(6) (DOUBLE) Track state vector
67			IFAIL (INTEGER) Error flag (ok if IFLAG=0)
68
69			Calculate particle-trajectory coordinates, for the state vector AL_P,
70			at z-coordinates set by the user.
71			Make use of GRKUTA2, routine to solve kinematic equations in a magnetic field
72			using Runge-Kutta method, adapted from GEANT3 library.
73			Same as DOTRACK and GRKUTA, modified in order to give as output the projected
74			angles and the track length.
75			The track length assigned to each track point i is evaluated between the i-th
76			plane and the closest lower(upper) one, reference plane included, if the plane
77			is above(below) the reference plane.
78
79			--------------------------------------------------------------------------------
80			CALL GUFLD(V,F)
81			--------------------------------------------------------------------------------
82
83			V(3) (REAL) x,y,z spatial coordinates (in m)
84			F(3) (REAL) Magnetic-field components (in kGauss)
85
86			Evaluate the three components of the magnetic field, by interpolating the
87			measured magnetic field map (used by GRKUTA/GRKUTA2).
88			It calls INTER_B.
89
90			--------------------------------------------------------------------------------
91			CALL INTER_B(X,Y,Z,B)
92			--------------------------------------------------------------------------------
93
94			X,Y,Z (DOUBLE) Spatial coordinates (in m)
95			B(3) (DOUBLE) Magnetic-field components (in T)
96
97			Evaluate the three components of the magnetic field, by interpolating the
98			measured magnetic field map.
99			It calls INTER_B_INNER and INTER_B_OUTER.