1 |
============ |
2 |
LIBTRK |
3 |
(R-1.01) |
4 |
============ |
5 |
|
6 |
LIBTRK is a library containing fortran subroutines of general utility. |
7 |
The library is stored in the directory trk-software-ground/bin. Both a static |
8 |
(.a) and dynamic (.so) libraries have been created. |
9 |
|
10 |
List of subroutines |
11 |
------------------- |
12 |
|
13 |
- TRACK |
14 |
- READ_B |
15 |
- INTER_B |
16 |
|
17 |
|
18 |
Usage: |
19 |
|
20 |
-------------------------------------------------------------------------------- |
21 |
CALL READ_B |
22 |
-------------------------------------------------------------------------------- |
23 |
|
24 |
Read the magnetic field maps (rz files) and fill related common blocks. |
25 |
===> To be called once, before using TRACK and INTER_B |
26 |
|
27 |
-------------------------------------------------------------------------------- |
28 |
CALL TRACK(*NPTOF,*ZIN,XOUT*,YOUT*,*AL_P,IFAIL*) |
29 |
-------------------------------------------------------------------------------- |
30 |
|
31 |
NPTOF (INTEGER) Num. points |
32 |
ZIN(NPTOF) (DOUBLE) Input z-coordinates (in cm) |
33 |
XOUT(NPTOF) (DOUBLE) Output x-coordinates (in cm) |
34 |
YOUT(NPTOF) (DOUBLE) Output y-coordinates (in cm) |
35 |
AL_P(6) (DOUBLE) Track state vector |
36 |
IFAIL (INTEGER) Error flag (ok if IFLAG=0) |
37 |
|
38 |
Calculate particle-trajectory coordinates, for the state vector AL_P, |
39 |
at z-coordinates set by the user. |
40 |
(Make use of GRKUTA, routine to solve kinematic equations in a magnetic field |
41 |
using Runge-Kutta method, adapted from GEANT3 library) |
42 |
|
43 |
-------------------------------------------------------------------------------- |
44 |
CALL INTER_B(*X,*Y,*Z,B*) |
45 |
-------------------------------------------------------------------------------- |
46 |
|
47 |
X,Y,Z (DOUBLE) Coordinates (in m) |
48 |
B(3) (DOUBLE) Magnetic field (in T) |
49 |
|
50 |
Evaluate the three components of the magnetic field, by interpolating the |
51 |
measured magnetic field map. |
52 |
(Used by TRACK) |