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