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	============
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.