9 |
|
|
10 |
c2f_error = 0 |
c2f_error = 0 |
11 |
|
|
12 |
if(C2F_DEBUG.eq.1)print * |
c if(C2F_DEBUG.eq.1)print * |
13 |
|
if(verbose)print * |
14 |
$ ,c2f_path(1:c2f_pathlen)//'trk-maskvk.dat' |
$ ,c2f_path(1:c2f_pathlen)//'trk-maskvk.dat' |
15 |
open(10, |
open(10, |
16 |
$ FILE=c2f_path(1:c2f_pathlen)//'trk-maskvk.dat' !<<<< C2F |
$ FILE=c2f_path(1:c2f_pathlen)//'trk-maskvk.dat' !<<<< C2F |
18 |
$ ,IOSTAT=iostat |
$ ,IOSTAT=iostat |
19 |
$ ) |
$ ) |
20 |
if(iostat.ne.0)then |
if(iostat.ne.0)then |
21 |
if(C2F_DEBUG.eq.1) |
c if(C2F_DEBUG.eq.1) |
22 |
|
if(debug) |
23 |
$ print*,'READVKMASK: *** Error opening file ***' |
$ print*,'READVKMASK: *** Error opening file ***' |
24 |
c2f_error = 1 |
c2f_error = 1 |
25 |
return |
return |
30 |
$ ) |
$ ) |
31 |
$ (mask_vk(iv,i),i=1,24) |
$ (mask_vk(iv,i),i=1,24) |
32 |
if(iostat.ne.0)then |
if(iostat.ne.0)then |
33 |
if(C2F_DEBUG.eq.1) |
c if(C2F_DEBUG.eq.1) |
34 |
|
if(debug) |
35 |
$ print*,'READVKMASK: *** Error reading file ***' |
$ print*,'READVKMASK: *** Error reading file ***' |
36 |
c2f_error = 1 |
c2f_error = 1 |
37 |
goto 1000 |
goto 1000 |
56 |
201 format('trk-LADDER',i1,'-mip.dat') |
201 format('trk-LADDER',i1,'-mip.dat') |
57 |
do ilad=1,nladders_view |
do ilad=1,nladders_view |
58 |
write(fname_param,201)ilad |
write(fname_param,201)ilad |
59 |
if(C2F_DEBUG.eq.1)print * |
c if(C2F_DEBUG.eq.1)print * |
60 |
|
if(VERBOSE)print * |
61 |
$ ,c2f_path(1:c2f_pathlen)!<<<< C2F |
$ ,c2f_path(1:c2f_pathlen)!<<<< C2F |
62 |
$ //fname_param(1:LNBLNK(fname_param)) |
$ //fname_param(1:LNBLNK(fname_param)) |
63 |
open(10, |
open(10, |
67 |
$ ,IOSTAT=iostat |
$ ,IOSTAT=iostat |
68 |
$ ) |
$ ) |
69 |
if(iostat.ne.0)then |
if(iostat.ne.0)then |
70 |
if(C2F_DEBUG.eq.1)print* |
c if(C2F_DEBUG.eq.1)print* |
71 |
|
if(DEBUG)print* |
72 |
$ ,'READMIPPARAM: *** Error opening file *** ' |
$ ,'READMIPPARAM: *** Error opening file *** ' |
73 |
c2f_error = 1 |
c2f_error = 1 |
74 |
return |
return |
100 |
201 format('charge-l',i1,'.dat') |
201 format('charge-l',i1,'.dat') |
101 |
do ilad=1,nladders_view |
do ilad=1,nladders_view |
102 |
write(fname_param,201)ilad |
write(fname_param,201)ilad |
103 |
if(C2F_DEBUG.eq.1)print * |
c if(C2F_DEBUG.eq.1)print * |
104 |
|
if(VERBOSE)print * |
105 |
$ ,c2f_path(1:c2f_pathlen)!<<<< C2F |
$ ,c2f_path(1:c2f_pathlen)!<<<< C2F |
106 |
$ //fname_param(1:LNBLNK(fname_param)) |
$ //fname_param(1:LNBLNK(fname_param)) |
107 |
open(10, |
open(10, |
151 |
|
|
152 |
******retrieve ANGULAR BINNING info |
******retrieve ANGULAR BINNING info |
153 |
fname_binning='binning.dat' |
fname_binning='binning.dat' |
154 |
if(C2F_DEBUG.eq.1) |
c if(C2F_DEBUG.eq.1) |
155 |
|
if(VERBOSE) |
156 |
$ print * |
$ print * |
157 |
$ ,c2f_path(1:c2f_pathlen) !<<<< C2F |
$ ,c2f_path(1:c2f_pathlen) !<<<< C2F |
158 |
$ //fname_binning(1:LNBLNK(fname_binning)) |
$ //fname_binning(1:LNBLNK(fname_binning)) |
163 |
$ ,IOSTAT=iostat |
$ ,IOSTAT=iostat |
164 |
$ ) |
$ ) |
165 |
if(iostat.ne.0)then |
if(iostat.ne.0)then |
166 |
if(C2F_DEBUG.eq.1) |
c if(C2F_DEBUG.eq.1) |
167 |
|
if(DEBUG) |
168 |
$ print*,'READETAPARAM: *** Error in opening file *** ' |
$ print*,'READETAPARAM: *** Error in opening file *** ' |
169 |
c2f_error = 1 |
c2f_error = 1 |
170 |
return |
return |
191 |
do ieta=2,4 !loop on eta 2,3,4 |
do ieta=2,4 !loop on eta 2,3,4 |
192 |
******retrieve correction parameters |
******retrieve correction parameters |
193 |
200 format(' Opening eta',i1,' files...') |
200 format(' Opening eta',i1,' files...') |
194 |
if(C2F_DEBUG.eq.1)write(*,200)ieta |
c if(C2F_DEBUG.eq.1)write(*,200)ieta |
195 |
|
if(VERBOSE)write(*,200)ieta |
196 |
|
|
197 |
201 format('eta',i1,'-bin',i1,'-l',i1,'.dat') |
201 format('eta',i1,'-bin',i1,'-l',i1,'.dat') |
198 |
202 format('eta',i1,'-bin',i2,'-l',i1,'.dat') |
202 format('eta',i1,'-bin',i2,'-l',i1,'.dat') |
208 |
$ ,IOSTAT=iostat |
$ ,IOSTAT=iostat |
209 |
$ ) |
$ ) |
210 |
if(iostat.ne.0)then |
if(iostat.ne.0)then |
211 |
if(C2F_DEBUG.eq.1) |
c if(C2F_DEBUG.eq.1) |
212 |
|
if(DEBUG) |
213 |
$ print*,'READETAPARAM: ** Error opening file ** ' |
$ print*,'READETAPARAM: ** Error opening file ** ' |
214 |
$ ,c2f_path(1:c2f_pathlen) |
$ ,c2f_path(1:c2f_pathlen) |
215 |
$ //fname_param(1:LNBLNK(fname_param)) |
$ //fname_param(1:LNBLNK(fname_param)) |
235 |
$ (feta4(ival,iv,ilad,iang),iv=1,nviews) |
$ (feta4(ival,iv,ilad,iang),iv=1,nviews) |
236 |
if(iostat.ne.0)then |
if(iostat.ne.0)then |
237 |
if(netaval.eq.0)then |
if(netaval.eq.0)then |
238 |
if(C2F_DEBUG.eq.1)print*,'READETAPARAM: ' |
c if(C2F_DEBUG.eq.1)print*,'READETAPARAM: ' |
239 |
|
if(DEBUG)print*,'READETAPARAM: ' |
240 |
$ //'*** Error reading file *** ' |
$ //'*** Error reading file *** ' |
241 |
$ ,c2f_path(1:c2f_pathlen) |
$ ,c2f_path(1:c2f_pathlen) |
242 |
$ //fname_param(1:LNBLNK(fname_param)) |
$ //fname_param(1:LNBLNK(fname_param)) |
286 |
do is=1,2 |
do is=1,2 |
287 |
|
|
288 |
write(fname_param,100)ilad,is |
write(fname_param,100)ilad,is |
289 |
if(C2F_DEBUG.eq.1)print * |
c if(C2F_DEBUG.eq.1)print * |
290 |
|
if(VERBOSE)print * |
291 |
$ ,c2f_path(1:c2f_pathlen) |
$ ,c2f_path(1:c2f_pathlen) |
292 |
$ //fname_param |
$ //fname_param |
293 |
|
|
363 |
|
|
364 |
c2f_error = 0 |
c2f_error = 0 |
365 |
|
|
366 |
if(C2F_DEBUG.eq.1)print * |
c if(C2F_DEBUG.eq.1)print * |
367 |
|
if(VERBOSE)print * |
368 |
$ ,c2f_path(1:c2f_pathlen)//'mech_pos.dat' |
$ ,c2f_path(1:c2f_pathlen)//'mech_pos.dat' |
369 |
open(10 |
open(10 |
370 |
$ ,FILE= |
$ ,FILE= |