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