!!!!!!!!!! Start is 20cm after photocathode !!!!!!!!!!!!

PARMMATCH: MARK
D1:   DRIFT, L=0.5
SOL1:   SOLENOID, L=0.05002, KS=+14.0
D2:   DRIFT, L=0.75
DS:   CSBEND, L=0.076, ANGLE=-0.261799388, K1=-5.93309223683602 &
	, TILT=0 &
	, E1=-0.2618, HGAP=0.0641927050701428, FINT=0.5, &
	  E2=0, EDGE_ORDER=1
D3:   DRIFT, L=0.75
SOL2:   SOLENOID, L=0.04735, KS=+15.0
D4:   DRIFT, L=0.595
WIEN: MATR, L=0.3095, filename="/a/itsuser/grames/elegant/wienfields/w_h_p00.dat", order=1
D5:   DRIFT, L=0.595
SOL3A:  SOLENOID, L=0.03277, KS=-12.0
SOL3B:  SOLENOID, L=0.03277, KS=+12.0
D6:   DRIFT, L=0.674
BUNCHER:  RFCA, L=0.152, VOLT=0, PHASE=180, FREQ=1497, BODY_FOCUS_MODEL="SRS", END1_FOCUS=1, END2_FOCUS=1, CHANGE_P0=1
D7:   DRIFT, L=0.674
A1:   MONITOR
SOL4A:  SOLENOID, L=0.03277, KS=-12.0
SOL4B:  SOLENOID, L=0.03277, KS=+12.0

!!!!!!!!!! Chopper Start !!!!!!!!!!!!

  D8C1: DRIFT, L=0.08225
D8C1S5: DRIFT, L=0.75
  SOL5: SOLENOID, L=0.0889, KS=12
D8C5MS: DRIFT, L=0.05
   AMS: MONITOR
D8MSC6: DRIFT, L=0.05
  SOL6: SOLENOID, L=0.0889, KS=12
D8S6C2: DRIFT, L=0.75
  D8C2: DRIFT, L=0.08225

SOL7A:  SOLENOID, L=0.03277, KS=-10.0
SOL7B:  SOLENOID, L=0.03277, KS=+10.0
D9:   DRIFT, L=0.5
A2:   MONITOR
D10:  DRIFT, L=0.5
QS1:  KQUAD, L=0.07, K1=0, TILT=0.785398163397448

!!!!!!!!!! Cryounit Start !!!!!!!!!!!!

D11:  DRIFT, L=0.7844
SRF1: RFCA, L=0.5, VOLT=3500000, PHASE=90, FREQ=1.497E+09 &
        , BODY_FOCUS_MODEL="SRS", END1_FOCUS=1, END2_FOCUS=1, CHANGE_P0=1
D12:  DRIFT, L=0.25
SRF2: RFCA, L=0.5, VOLT=3500000, PHASE=90, FREQ=1.497E+09 &
        , BODY_FOCUS_MODEL="SRS", END1_FOCUS=1, END2_FOCUS=1, CHANGE_P0=1
D13:  DRIFT, L=0.23439
D14:  drift, L=2.0

!!!!!!!!!! MeV Start !!!!!!!!!!!!

b1: CSBEND, L=0.15, ANGLE=+0.436332313, E1=0.218166157, E2=0.218166157, TILT=-1.570796327, N_KICKS=20,
b2: CSBEND, L=0.15, ANGLE=+0.436332313, E1=0.218166157, E2=0.218166157, TILT=+1.570796327, N_KICKS=20,

qa: quadrupole, l=0.150, k1=0
qb: quadrupole, l=0.075, k1=0

mda: drift, l=0.501120203

MKMATCH: mark,fitpoint=1
MKFINAL: mark,fitpoint=1

q1: quadrupole, l=0.150, k1= 1
q2: quadrupole, l=0.150, k1=-1
q3: quadrupole, l=0.150, k1= 1
q4: quadrupole, l=0.150, k1=-1

qw: quadrupole, l=0.150, k1=-1
qx: quadrupole, l=0.150, k1= 1
qy: quadrupole, l=0.150, k1=-1
qz: quadrupole, l=0.150, k1= 1

md1: drift, l=1.0
md2: drift, l=1.0

! updated raster 9/22 for 6GeV length (coil=25cm, pipe=40cm)
    rasterx: drift, l=0.4
    rastery: drift, l=0.4
    hd_diag: drift, l=2.87826

  ibc_start: monitor
     ibc_d1: drift, l=2.12174
  ibc_target: monitor
! ibc_target: watch, filename="ibc_target.watch", mode=coordinate
     ibc_d2: drift, l=0.34085
    ibc_sol: solenoid, l=0.40, b=0.00
   ibc_exit: monitor
     d2dump: drift, l=1.00000

!!!!!!!!!! gun2cryounit !!!!!!!!!!!!

gun2chopper: line=(PARMMATCH, D1, SOL1, D2, DS, D3, SOL2, D4, &
WIEN, D5, SOL3A, SOL3B, D6, BUNCHER, D7, A1, SOL4A, SOL4B)

chopper2unit: line=(D8C1, D8C1S5, SOL5, D8C5MS, AMS, D8MSC6, &
SOL6, D8S6C2, D8C2, SOL7A, SOL7B)

unit2mev: line=(D9, A2, D10, QS1, D11, SRF1, D12, SRF2, D13, D14)

anne: line=(gun2chopper,chopper2unit,unit2mev)

!!!!!!!!!! cryounit2hdicedump !!!!!!!!!!!!

tele: line=(md1, q1, md2, q2, md2, q3, md2, q4, md1)
chic: line=(b1, mda, qa, mda, qb, MKMATCH, qb, mda, qa, mda, b2)
elet: line=(md1, qw, md2, qx, md2, qy, md2, qz, md1)
hdic: line=(rasterx, rastery, hd_diag, &
ibc_start, ibc_d1, MKFINAL, ibc_target, ibc_sol, ibc_exit, d2dump)

george: line=(tele, chic, elet, hdic)