--- DarthVader/TrackerLevel2/src/F77/analysissubroutines.f 2007/08/20 16:07:16 1.28 +++ DarthVader/TrackerLevel2/src/F77/analysissubroutines.f 2007/08/28 13:25:50 1.30 @@ -859,7 +859,7 @@ if( nsatstrips(icx).gt.0 )then stripx = stripx + cog(4,icx) resxPAM = pitchX*1e-4/sqrt(12.) - cc=cog(4,icx) +cc cc=cog(4,icx) c$$$ print*,icx,' *** ',cc c$$$ print*,icx,' *** ',resxPAM endif @@ -999,7 +999,7 @@ if( nsatstrips(icy).gt.0 )then stripy = stripy + cog(4,icy) resyPAM = pitchY*1e-4/sqrt(12.) - cc=cog(4,icy) +cc cc=cog(4,icy) c$$$ print*,icy,' *** ',cc c$$$ print*,icy,' *** ',resyPAM endif @@ -1253,6 +1253,7 @@ c$$$ PFAx = 'COG4'!PFA c$$$ PFAy = 'COG4'!PFA + if(icx.gt.nclstr1.or.icy.gt.nclstr1)then print*,'xyzpam: ***WARNING*** clusters ',icx,icy $ ,' does not exists (nclstr1=',nclstr1,')' @@ -1265,6 +1266,8 @@ call idtoc(pfaid,PFAx) call idtoc(pfaid,PFAy) +cc print*,PFAx,PFAy + c$$$ call xyz_PAM(icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy) c$$$ print*,icx,icy,sensor,PFAx,PFAy,ax,ay,bfx,bfy @@ -2240,6 +2243,9 @@ ym3=yPAM zm3=zPAM * find the circle passing through the three points +c$$$ call tricircle(3,xp,zp,angp,resp,chi +c$$$ $ ,xc,zc,radius,iflag) + iflag = DEBUG call tricircle(3,xp,zp,angp,resp,chi $ ,xc,zc,radius,iflag) c print*,xc,zc,radius