! New standard file using matlab file HBSRS_Ver03 oi july2025 
ap: maxamp,x_max=0.013,y_max=0.013
ap1: maxamp,x_max=0.010,y_max=0.0025
ap2: maxamp,x_max=0.013,y_max=0.010
ap3: maxamp,x_max=0.013,y_max=0.007
ap4: maxamp,x_max=0.007,y_max=0.007
! Elliptical aperture offset in horizontal plane
MM:MARKER
MD:MARKER
MA:MARKER

Loct1=0.09
Loct2=0.09

!  Drifts

D1a : EDRIFT, L= 1.3092
D1b : EDRIFT, L= 0.5
D1c : EDRIFT, L= 1.3092
D2 : EDRIFT, L= 0.11600
D3 : EDRIFT, L= 0.251000
D4 : EDRIFT, L= 0.089500 
D5 : EDRIFT, L= 0.089500
D6 : EDRIFT, L= 0.100700
D7a : EDRIFT, L= 0.0853
D7b : EDRIFT, L= 0.5
D7c : EDRIFT, L= 0.0853
D8 : EDRIFT, L= 0.099500

!D9A & D9B LENGTH = 0.145-Loct1/2
D9A : EDRIFT, L= 0.1
D9B : EDRIFT, L= 0.1 

D10 : EDRIFT, L= 0.0835 
D11 : EDRIFT, L= 0.494000 
D11a : EDRIFT, L= 0.1
D11b : EDRIFT, L= 0.040
D11c : EDRIFT, L= 0.3540
D12 : EDRIFT, L= 0.100700
D13 : EDRIFT, L= 0.089500 
D14 : EDRIFT, L= 0.352000
D15 : EDRIFT, L= 0.115500 
D16 : EDRIFT, L= 0.120800
D17 : EDRIFT, L= 0.322500
D18 : EDRIFT, L= 0.253500
!DOCT LENGTH = 0.04500-Loct2/2
DOCT : EDRIFT, L=0.0
! Dipoles

BM1LGB1   : CSBEND,L=0.360000, ANGLE = 0.0089870172, K1 = 0.0000, E1= 0.0158148039, E2= -0.0068296922, N_SLICES= 40 
BM1LGB2   : CSBEND,L=0.360000, ANGLE = 0.0071184443, K1 = 0.0000, E1= 0.0068296922, E2= 0.0002905702, N_SLICES= 40 
BM1LGB3   : CSBEND,L=0.360000, ANGLE = 0.0055516956, K1 = 0.0000, E1= -0.0002905702, E2= 0.0058425926, N_SLICES= 40 
BM1LGB4   : CSBEND,L=0.360000, ANGLE = 0.0052056249, K1 = 0.0000, E1= -0.0058425926, E2= 0.0110464226, N_SLICES= 40 
BM1LGB5   : CSBEND,L=0.360000, ANGLE = 0.0047686435, K1 = 0.0000, E1= -0.0110464226, E2= 0.0158148039, N_SLICES= 40 

BM2LGB1   : CSBEND,L=0.360000, ANGLE = 0.0034568582, K1 = 0.0000, E1= 0.0172652665, E2= -0.0138087713, N_SLICES= 40 
BM2LGB2   : CSBEND,L=0.360000, ANGLE = 0.0049638926, K1 = 0.0000, E1= 0.0138087713, E2= -0.0088438338, N_SLICES= 40 
BM2LGB3   : CSBEND,L=0.360000, ANGLE = 0.0070636541, K1 = 0.0000, E1= 0.0088438338, E2= -0.0017795585, N_SLICES= 40 
BM2LGB4   : CSBEND,L=0.360000, ANGLE = 0.0077549465, K1 = 0.0000, E1= 0.0017795585, E2= 0.0059731992, N_SLICES= 40 
BM2LGB5   : CSBEND,L=0.360000, ANGLE = 0.0112935045, K1 = 0.0000, E1= -0.0059731992, E2= 0.0172652665, N_SLICES= 40 

BM3   : CSBEND,L=1.0000, ANGLE = 0.0240558695, K1 = -1.8054710000000, E1= 0.0120279348, E2= 0.0120279348, N_SLICES= 40 
BM41   : CSBEND,L=0.4000, ANGLE = 0.0130201546, K1 = -1.0278610000000, E1= 0.0130201546, E2= 0.0000, N_SLICES= 40
BM42   : CSBEND,L=0.4000, ANGLE = 0.0130201546, K1 = -1.0278610000000, E1= 0.0000, E2= 0.0130201546, N_SLICES= 40 



!  Quadrupoles

Q1   : kquad, l = 0.30, k1 = 2.37256314, N_SLICES=20
Q2   : kquad, l = 0.250, k1 = -2.229688834, N_SLICES=20
Q3   : kquad, l = 0.20, k1 = -1.742713, N_SLICES=20
Q6   : kquad,  L= 0.250, K1= -2.279661, N_SLICES= 20

Q4   : CSBEND, L=0.20, ANGLE = -0.0008726646, K1 = 2.598000, E1= 0.0000, E2= 0.0000, N_SLICES= 20
Q5   : CSBEND, L=0.20, ANGLE = -0.0008726646, K1 = 2.344433, E1= 0.0000, E2= 0.0000, N_SLICES= 20
Q7   : CSBEND,  L=0.450, ANGLE = -0.0015707963, K1 = 3.718301, E1= 0.0000, E2= 0.0000, N_SLICES= 40
Q8   : CSBEND, L=0.50, ANGLE = -0.0017453293, K1 = 3.722488, E1= 0.0000, E2= 0.0000, N_SLICES= 40

! Sextupoles

S1   : KSEXT,  L= 0.20, K2= -136.177172, N_SLICES= 40
S2   : KSEXT,  L= 0.20, K2= 151.680556, N_SLICES= 40 
S3   : KSEXT,  L= 0.20, K2= -118.92727, N_SLICES= 40 

FCHV1   : KSEXT,  L= 0.10, K2= 0.0000000000000000, N_SLICES= 40 
FCHV2   : KSEXT,  L= 0.10, K2= 0.0000000000000000, N_SLICES= 40 

!  Octopoles

!OXY3   : MULT,  L= 0.090000000000, ORDER= 3, KNL= -1429.08, N_SLICES= 40
!OXY4   : MULT,  L= 0.090000000000, ORDER= 3, KNL= 214.38, N_SLICES= 40
OXY3   : koct, l=0.09, k3 = -15878.66666666666666666667, N_SLICES=40
OXY4   : koct, l=0.09, k3 = 2382.0000, N_SLICES=40



rf: rfca,freq="9.598208e+02 6000 mev / beta.p c_mks * / rec 1600 *",volt=7.00e6,phase="180 2.76541554 7.00 / dasin -"

BPM: MONI
HCOR: HKICK,L=0.0
VCOR: VKICK,L=0.0
BM1LGB : line = ( BM1LGB1, BM1LGB2, BM1LGB3, BM1LGB4, BM1LGB5)
BM2LGB : line = ( BM2LGB1, BM2LGB2, BM2LGB3, BM2LGB4, BM2LGB5)

ICELLI   : line = ( ap, D1a, D1b, D1c, Q1, D2, FCHV1, D3, Q2, D4, BM1LGB1, &
                  BM1LGB2, BM1LGB3, BM1LGB4, BM1LGB5, D5, Q3, D6, S1, &
                  D7a, D7b, D7c, Q4, D8, S2, D9A, OXY3, D9B, Q5, D11a, D11b, D11c, DOCT, &
                  OXY4, DOCT, D10, S3, D12, Q6, D13, BM2LGB1, BM2LGB2, BM2LGB3, &
                  BM2LGB4, BM2LGB5, ap, ap4, D14, Q7, D15, BM3, D16, Q8, &
                  D17, FCHV2, D18, BM41, BM42, D18, FCHV2, D17, Q8, &
                  D16, BM3, D15, Q7, D14, ap, BM2LGB5, BM2LGB4, BM2LGB3, &
                  BM2LGB2, BM2LGB1, D13, Q6, D12, S3, D10, DOCT, &
                  OXY4, DOCT, D11, Q5, D9B, OXY3, D9A, S2, D8, Q4, D7c, D7b, D7a, S1, &
                  D6, Q3, D5, BM1LGB5, BM1LGB4, BM1LGB3, BM1LGB2, BM1LGB1, &
                  D4, Q2, D3, FCHV1, D2, Q1, D1c, D1b, D1a )

ICELL   : LINE = ( ap, D1a, D1b, D1c, Q1, D2, FCHV1, D3, Q2, D4, BM1LGB1, &
                  BM1LGB2, BM1LGB3, BM1LGB4, BM1LGB5, D5, Q3, D6, S1, &
                  D7a, D7b, D7c, Q4, D8, S2, D9A, OXY3, D9B, Q5, D11, DOCT, &
                  OXY4, DOCT, D10, S3, D12, Q6, D13, BM2LGB1, BM2LGB2, BM2LGB3, &
                  BM2LGB4, BM2LGB5, ap, ap4, D14, Q7, D15, BM3, D16, Q8, &
                  D17, FCHV2, D18, BM41, BM42, D18, FCHV2, D17, Q8, &
                  D16, BM3, D15, Q7, D14, ap, BM2LGB5, BM2LGB4, BM2LGB3, &
                  BM2LGB2, BM2LGB1, D13, Q6, D12, S3, D10, DOCT, &
                  OXY4, DOCT, D11, Q5, D9B, OXY3, D9A, S2, D8, Q4, D7c, D7b, D7a, S1, &
                  D6, Q3, D5, BM1LGB5, BM1LGB4, BM1LGB3, BM1LGB2, BM1LGB1, &
                  D4, Q2, D3, FCHV1, D2, Q1, D1c, D1b, D1a )



W1: watch, filename="%s.w1", mode="coord", interval=1
W3: watch, filename="%s.w3", mode="centroid", interval=1
W2: watch, filename="%s.w2", mode="parameters", interval=1

malin: malign,on_pass=0
sr: sreffects,qexcitation=0
ring             : line = (ICELLi,31*ICELL)
ringw            : line = (malin,w1,ring)
ringrf           : line=(ring,sr,rf)

RETURN
