1 |
pam-fi |
1.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) |