#!/bin/csh
set nonomatch

sddssequence -pipe=out -define=delta -sequence=begin=-1e-6,end=1e-6,n=3 \
    | sddsprocess -pipe -define=col,absDelta,"delta abs" \
    | sddssort -pipe=in delta.sdds -col=absDelta,incr 

elegant nls_2pass_lattice_work_mal.ele

sddsxref nls_inj_fel.twi nls_inj_fel.fin -pipe=out -leave=* -transfer=parameter,MAL.DP -rename=parameter,MAL.DP=delta \
    | sddsprocess -print=column,Rootname,%s\#%ld,ElementName,ElementOccurence -pipe=in nls_inj_fel.tmp

sddsvslopes nls_inj_fel.tmp -pipe=out -columns=betax,betay,alphax,alphay,etax,etaxp -independentVariable=delta \
| sddsxref -pipe nls_inj_fel.tmp -take=s \
| sddsprocess -pipe=in nls_inj_fel.fit \
    "-define=col,bx,betaxSlope betaxIntercept /" \
    "-define=col,by,betaySlope betayIntercept /" \
    "-define=col,ax,alphaxSlope alphaxIntercept bx * -" \
    "-define=col,ay,alphaySlope alphayIntercept by * -" \
    "-define=col,Wx,ax sqr bx sqr + sqrt" \
    "-define=col,Wy,ay sqr by sqr + sqrt" 

\rm nls_inj_fel.tmp

sddsplot -column=s,W? -sep nls_inj_fel.fit -title=





  
     
