! Physica command file.
! Calculates x position residuals and z position residuals.
! Former "senw_stat.pcm" modified for the TWIST database.
!
! Frame data are supplied in a file with 2 columns of numbers: x u
! Comments are to be stripped by database.
! Guard wire lines can be stripped as well if only sense wire 
! statistics is to be shown.

cl
orien land
disable journal


file=`191PC0P.txt'
inquire `enter filename >' file

title_line = `global title'
inquire `enter global title >' title_line

read file x u
NW=last(x)          ! number of wires in the plane

window 0
 set
  %xloc   50
  %yloc   97.5
  %txthit  2
   txtang  0
   cursor -2

 text title_line
 set cursor 1

generate n 4 1 NW+3

z=x-u
! list n x u z

! fit data and plot residuals in x direction
scalar\vary a b
fit x=a*n+b
fit\update xfit
! list x xfit
resid=x-xfit
! list n x resid

WINDOW 5

 set
  %xloc   50
  %yloc   91
  %txthit  5
   txtang  0
   cursor -2

 text `x'

SET nsxinc 5
SET nsyinc 2
SET %ytics -0.5
SET %yticl -1
SET %xtics -0.5
SET %xticl -1
SET %xnumsz 2.0
SET %ynumsz 2.0

!goto Zpos

!set histyp 2
set cursor 1
set pchar 1
label\xaxis `wire number'
label\yaxis `residuals, <mu>m'
graph n 1000*resid
statistics resid xrms\rms xsig\sdev xmin\min xi\imin xmax\max xm\imax
xi=xi+2
xm=xm+2
! show xmin xi xmax xm xsig
bin\nbins resid residb residcount 400 -.20 .20

a=abs(a)
xloc=50
yloc=20
 set
   %xloc xloc
   %yloc yloc
   %txthit  3
    cursor -1
          
  text `s='//rchar(a,`f6.4')//`<Pm>'//rchar(xsig,`f6.4')//` mm'
          
          
!goto END

WINDOW 7

!GOTO Zpos

set histyp 2
set pchar 0
label\xaxis `residual, <mu>m'
label\yaxis `number of wires'

a1=2000*xmin
a2=2000*xmax
! show a1 a2

scales a1 a2 0 0 60 0

set autoscale yaxis on

!SET autoscale on

graph\histogram 1000*residb residcount
xloc=65
yloc=85
 set
  %xloc xloc
  %yloc yloc
  %txthit  3
   cursor -1


! text `rms='//rchar(1000*xrms,`f5.2')//`<mu>m'
!yloc=yloc-8
! set
!  %yloc yloc

  text `<sigma>='//rchar(1000*xsig,`f5.2')//`<mu>m'
yloc=yloc-8
 set
  %yloc yloc

  text `min='//rchar(1000*xmin,`f5.1')//`<mu>m'
yloc=yloc-8
 set
  %yloc yloc

  text `wire '//rchar(xi,`f4.0')
yloc=yloc-8
 set
  %yloc=yloc

  text `max='//rchar(1000*xmax,`f5.1')//`<mu>m'
yloc=yloc-8
 set
  %yloc=yloc

 text `wire '//rchar(xm,`f4.0')
  set cursor 1

Zpos:

! fit data and plot residuals in z direction
scalar\vary c d
fit z=c*n+d
fit\update zfit
residz=z-zfit
! list n z residz

WINDOW 6

 set
  %xloc   50
  %yloc   93
  %txthit  5
   txtang  0
   cursor -0

 text `z'

SET autoscale on

SET histyp 0
 set cursor 1
set pchar 1
label\xaxis `wire number'
label\yaxis `residual, <mu>m'
graph n 1000*residz
statistics residz zmin\min zi\imin zmax\max zm\imax zsig\sdev
zi=zi+2
zm=zm+2
! show zmin zi zmax zm zsig
bin\nbins residz residzb residzcount 400 -.2 .2


!goto END

WINDOW 8

set histyp 2
set pchar 0

b1=2000*zmin
b2=2000*zmax
! show b1 b2

scales b1 b2 0 0 60 0

set autoscale yaxis on

label\xaxis `residuals, <mu>m'
label\yaxis `number of wires'
graph 1000*residzb residzcount
xloc=65
yloc=85
 set
  %xloc xloc
  %yloc yloc
  %txthit  3
  cursor -1

 text `<sigma>='//rchar(1000*zsig,`f5.2')//`<mu>m'
yloc=yloc-8
 set
  %yloc yloc

 text `min='//rchar(1000*zmin,`f6.1')//`<mu>m'
yloc=yloc-8
 set
  %yloc yloc

 text `wire '//rchar(zi,`f4.0')
yloc=yloc-8
 set
  %yloc yloc

 text `max='//rchar(1000*zmax,`f6.1')//`<mu>m'
yloc=yloc-8
 set
  %yloc=yloc

 text `wire '//rchar(zm,`f4.0')
 set cursor 1

END:
default