#!/bin/bash

# Make a file to allow varying initial delta between -1% and +1%. Sort so delta=0 is first for
# proper fiducialization
sddssequence -pipe=out -define=deltaInput -seq=begin=-0.01,end=0.01,n=21 | sddsprocess -pipe -define=col,deltaInputAbs,"deltaInput abs" | 
sddssort -pipe=in -col=deltaInputAbs -col=deltaInput R56input.sdds

# Run matrix analysis for R56
elegant SF-Injector.ele
# Run tracking to give Cs vs *initial* delta
elegant SF-InjectorR56Tracking.ele

# Fit Cs vs *initial* delta at each s location
sddsprocess SF-InjectorR56Tracking.cen -pipe=out -clip=1,0 -process=Cdelta,first,%s0  \
    | sddsregroup -pipe -newparam=s -newcol=Cdelta0 \
    | sddspfit -pipe=in SF-InjectorR56Tracking.cenfit -col=Cdelta0,Cs -terms=3 -copyP

# Plot comparison.
\rm checkR56.png checkR56.eps
sddsplot3 checkR56 -graph=line,vary -param=s,Coefficient01 SF-InjectorR56Tracking.cenfit -graph=sym -legend=spec=Tracking \
    -col=s,R56 Aramis.mat -legend=spec=Matrix

