1 |
pam-fi |
1.1 |
subroutine readalignparam |
2 |
|
|
|
3 |
|
|
include '../common/commontracker.f' |
4 |
|
|
include '../common/common_align.f' |
5 |
|
|
|
6 |
|
|
character*50 fname_param |
7 |
|
|
integer sensor |
8 |
|
|
c character*120 cmd1 |
9 |
|
|
c character*120 cmd2 |
10 |
|
|
|
11 |
|
|
c 100 format('../common/align_param/parameters_l',i1,'s',i1,'.dat') |
12 |
|
|
100 format('parameters_l',i1,'s',i1,'.dat') |
13 |
|
|
|
14 |
|
|
do ilad=1,nladders_view |
15 |
|
|
do is=1,2 |
16 |
|
|
|
17 |
|
|
write(fname_param,100)ilad,is |
18 |
|
|
c$$$ cmd1='cp $TRK_GRND/source/common/align_param/' |
19 |
|
|
c$$$ $ //fname_param(1:LNBLNK(fname_param))//' .' |
20 |
|
|
c$$$ call system(cmd1) |
21 |
|
|
print *,'Opening file: ',fname_param |
22 |
|
|
|
23 |
|
|
open(10, |
24 |
|
|
$ FILE='./bin-aux/'//fname_param(1:LNBLNK(fname_param)) |
25 |
|
|
$ ,STATUS='UNKNOWN' |
26 |
|
|
$ ) |
27 |
|
|
|
28 |
|
|
do ip=1,nplanes |
29 |
|
|
* |
30 |
|
|
* NB! NB! NB! NB! NB! |
31 |
|
|
* The file labelled for example "l1s1" include the alignment |
32 |
|
|
* parameters of the column of sensors: |
33 |
|
|
* |
34 |
|
|
* PLANE LADDER SENSOR |
35 |
|
|
* ------------------- |
36 |
|
|
* 1 1 2 |
37 |
|
|
* 2 1 1 |
38 |
|
|
* 3 1 1 |
39 |
|
|
* 4 1 1 |
40 |
|
|
* 5 1 1 |
41 |
|
|
* 6 1 1 |
42 |
|
|
* |
43 |
|
|
* This is becouse the plane 1 (bottom plane) is up-side-down |
44 |
|
|
* |
45 |
|
|
sensor=is |
46 |
|
|
if(ip.eq.1)sensor=mod(is,2)+1 |
47 |
|
|
read(10,*)omega(ip,ilad,sensor) |
48 |
|
|
read(10,*)beta(ip,ilad,sensor) |
49 |
|
|
read(10,*)gamma(ip,ilad,sensor) |
50 |
|
|
|
51 |
|
|
c N.B. I convert angles from microradiants to radiant |
52 |
|
|
omega(ip,ilad,sensor)=omega(ip,ilad,sensor)/1.d6 |
53 |
|
|
beta(ip,ilad,sensor)=beta(ip,ilad,sensor)/1.d6 |
54 |
|
|
gamma(ip,ilad,sensor)=gamma(ip,ilad,sensor)/1.d6 |
55 |
|
|
|
56 |
|
|
read(10,*)dx(ip,ilad,sensor) |
57 |
|
|
read(10,*)dy(ip,ilad,sensor) |
58 |
|
|
read(10,*)dz(ip,ilad,sensor) |
59 |
|
|
enddo |
60 |
|
|
|
61 |
|
|
close(10) |
62 |
|
|
c$$$ cmd2='rm -f ' |
63 |
|
|
c$$$ $ //fname_param(1:LNBLNK(fname_param)) |
64 |
|
|
c$$$ call system(cmd2) |
65 |
|
|
enddo |
66 |
|
|
enddo |
67 |
|
|
|
68 |
|
|
|
69 |
|
|
|
70 |
|
|
end |