--- DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/05/14 11:03:06 1.10 +++ DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/05/15 16:22:18 1.11 @@ -164,7 +164,8 @@ end *** * * * *** * * * *** * * * *** * * * *** * * * *** * * * *** - real function ris_eta(ic,angle) +c real function riseta(ic,angle) + real function riseta(iview,angle) *-------------------------------------------------------------- * this function returns the average spatial resolution * (in cm) for the ETA algorithm (function pfaeta(ic,angle)) @@ -179,35 +180,38 @@ include 'level1.f' include 'calib.f' - ris_eta = 0 + riseta = 0 - if(mod(int(VIEW(ic)),2).eq.1)then !Y-view +c if(mod(int(VIEW(ic)),2).eq.1)then !Y-view + if(mod(iview,2).eq.1)then !Y-view if( abs(angle).ge.e2fay.and.abs(angle).le.e2tay )then - ris_eta = risy_eta2(angle) + riseta = risy_eta2(angle) elseif( abs(angle).ge.e3fay.and.abs(angle).le.e3tay )then - ris_eta = risy_cog(angle) !ATTENZIONE!! + riseta = risy_cog(angle) !ATTENZIONE!! elseif( abs(angle).ge.e4fay.and.abs(angle).le.e4tay )then - ris_eta = risy_cog(angle) !ATTENZIONE!! + riseta = risy_cog(angle) !ATTENZIONE!! else - ris_eta = risy_cog(angle) + riseta = risy_cog(angle) endif else !X-view if( abs(angle).ge.e2fax.and.abs(angle).le.e2tax )then - ris_eta = risx_eta2(angle) + riseta = risx_eta2(angle) elseif( abs(angle).ge.e3fax.and.abs(angle).le.e3tax )then - ris_eta = risx_eta3(angle) + riseta = risx_eta3(angle) elseif( abs(angle).ge.e4fax.and.abs(angle).le.e4tax )then - ris_eta = risx_eta4(angle) + riseta = risx_eta4(angle) else - ris_eta = risx_cog(angle) + riseta = risx_cog(angle) endif endif + print*,'---- ',riseta,iview,angle + 100 return end