--- DarthVader/TrackerLevel2/src/F77/readallparam.f 2006/05/30 16:30:37 1.2 +++ DarthVader/TrackerLevel2/src/F77/readallparam.f 2007/08/17 14:36:06 1.3 @@ -149,7 +149,9 @@ c2f_error = 0 -******retrieve ANGULAR BINNING info +* ====================================================== +* read ANGULAR BINNING info +* ====================================================== fname_binning='binning.dat' c if(C2F_DEBUG.eq.1) if(VERBOSE) @@ -188,8 +190,10 @@ +* ====================================================== +* read eta2-eta3-eta4 parameters +* ====================================================== do ieta=2,4 !loop on eta 2,3,4 -******retrieve correction parameters 200 format(' Opening eta',i1,' files...') c if(C2F_DEBUG.eq.1)write(*,200)ieta if(VERBOSE)write(*,200)ieta @@ -214,7 +218,8 @@ $ ,c2f_path(1:c2f_pathlen) $ //fname_param(1:LNBLNK(fname_param)) c2f_error = 1 - return +c return + goto 2001 endif netaval=0 do ival=1,netavalmax @@ -248,12 +253,65 @@ netaval=netaval+1 enddo 2000 close(10) + 2001 continue * print*,'... done' enddo enddo enddo !end loop on eta 2,3,4 +* ====================================================== +* read landi correction +* ====================================================== + if(VERBOSE)print*,' Opening landi-correction files...' + + 301 format('lcorr-l',i1,'.dat') + do ilad=1,nladders_view + + write(fname_param,301)ilad + open(10, + $ FILE=c2f_path(1:c2f_pathlen) + $ //fname_param(1:LNBLNK(fname_param)) + $ ,STATUS='OLD' + $ ,IOSTAT=iostat + $ ) + if(iostat.ne.0)then + if(.true.) + $ print*,'READETAPARAM: ** Error opening file ** ' + $ ,c2f_path(1:c2f_pathlen) + $ //fname_param(1:LNBLNK(fname_param)) + c2f_error = 1 +c return + goto 3001 + endif + + do iang=1,nangbin + + read(10,*,IOSTAT=iostat) + $ ii,aal,aar,(fcorr(iv,ilad,iang),iv=1,nviews) + + write(*,*) + $ ii,aal,aar,(fcorr(iv,ilad,iang),iv=1,nviews) + + if(ii.ne.iang)print*,'parametri eta incasinati' + if(aal.ne.angl(iang))print*,'parametri eta incasinati' + if(aar.ne.angr(iang))print*,'parametri eta incasinati' + + if(iostat.ne.0)then + if(.true.) + $ print*,'READETAPARAM: ' + $ //'*** Error reading file *** ' + $ ,c2f_path(1:c2f_pathlen) + $ //fname_param(1:LNBLNK(fname_param)) + c2f_error = 1 + goto 3000 + endif + + enddo !end loop on angular bins + 3000 close(10) + 3001 continue + + enddo ! end loop on ladders return end