622 |
enddo |
enddo |
623 |
if(DEBUG.EQ.1) |
if(DEBUG.EQ.1) |
624 |
$ print*,'pfaeta2 *** warning *** angle out of range: ',angle |
$ print*,'pfaeta2 *** warning *** angle out of range: ',angle |
625 |
if(angle.lt.angL(1))iang=1 |
if(angle.le.angL(1))iang=1 |
626 |
if(angle.gt.angR(nangbin))iang=nangbin |
if(angle.ge.angR(nangbin))iang=nangbin |
627 |
98 continue !jump here if ok |
98 continue !jump here if ok |
628 |
|
|
629 |
|
|
740 |
enddo |
enddo |
741 |
if(DEBUG.EQ.1) |
if(DEBUG.EQ.1) |
742 |
$ print*,'pfaeta3 *** warning *** angle out of range: ',angle |
$ print*,'pfaeta3 *** warning *** angle out of range: ',angle |
743 |
if(angle.lt.angL(1))iang=1 |
if(angle.le.angL(1))iang=1 |
744 |
if(angle.gt.angR(nangbin))iang=nangbin |
if(angle.ge.angR(nangbin))iang=nangbin |
745 |
98 continue !jump here if ok |
98 continue !jump here if ok |
746 |
|
|
747 |
|
|
856 |
enddo |
enddo |
857 |
if(DEBUG.EQ.1) |
if(DEBUG.EQ.1) |
858 |
$ print*,'pfaeta4 *** warning *** angle out of range: ',angle |
$ print*,'pfaeta4 *** warning *** angle out of range: ',angle |
859 |
if(angle.lt.angL(1))iang=1 |
if(angle.le.angL(1))iang=1 |
860 |
if(angle.gt.angR(nangbin))iang=nangbin |
if(angle.ge.angR(nangbin))iang=nangbin |
861 |
98 continue !jump here if ok |
98 continue !jump here if ok |
862 |
|
|
863 |
|
|
1961 |
enddo |
enddo |
1962 |
if(DEBUG.eq.1) |
if(DEBUG.eq.1) |
1963 |
$ print*,'pfacorr *** warning *** angle out of range: ',angle |
$ print*,'pfacorr *** warning *** angle out of range: ',angle |
1964 |
if(angle.lt.angL(1))iang=1 |
if(angle.le.angL(1))iang=1 |
1965 |
if(angle.gt.angR(nangbin))iang=nangbin |
if(angle.ge.angR(nangbin))iang=nangbin |
1966 |
98 continue !jump here if ok |
98 continue !jump here if ok |
1967 |
|
|
1968 |
pfacorr = fcorr(iview,lad,iang) |
pfacorr = fcorr(iview,lad,iang) |