!------------------------------------------------------------------------------------------------------
! ELEGANT lattice for Bubble Chamber Experiment
! August 3, 2015
! C. Tennant
!------------------------------------------------------------------------------------------------------
!
! - using songsheet ACC2008000-1100.pdf

"cle" ! clear whole RPN stack for safety

% 1 atan 4 * sto pi
% pi 180 / sto cdtor 
% 180 pi / sto crtod

TWI: TWISS, BETAX=22.96, ALPHAX=10.14, BETAY=29.69, ALPHAY=-4.56, FROM_BEAM=1
EM:  EMITTANCE, EMITNX=9.55e-007, EMITNY=1.00e-007

W0: WATCH, FILENAME="%s.W0", MODE=COORD
W1: WATCH, FILENAME="%s.W1", MODE=COORD
W2: WATCH, FILENAME="%s.W2", MODE=COORD

!---------------------------------------------------------------
!     Region 0: EXIT OF CRYOUNIT TO ENTRANCE OF 5 MEV DIPOLE
!---------------------------------------------------------------
! Assume nothing has changed in this region since ~2013 (check with Riad)

! DRIFTS

D0:        DRIFT, L=0.305322
D1:        DRIFT, L=0.224795
D2:        DRIFT, L=0.294893
D3:        DRIFT, L=0
D4:        DRIFT, L=0.1068
D5:        DRIFT, L=1.205
D6:        DRIFT, L=0.6788
D7:        DRIFT, L=0.0743
D8:        DRIFT, L=0.8521
D9:        DRIFT, L=0.191895
D10:       DRIFT, L=0.132852
D11:       DRIFT, L=0.0550586
D12:       DRIFT, L=1.01011

! VIEWERS/BPMS/HARPS

ITV0L01:   WATCH, FILENAME="%s.TV0L01", MODE=COORD
ITV0L02:   WATCH, FILENAME="%s.TV0L02", MODE=COORD
IPM0L01:   MONITOR, L=0
IPM0L02:   MONITOR, L=0

! CORRECTORS

MBH0L01H:  HKICK, L=0, KICK=0, TILT=0
MBH0L01V:  VKICK, L=0, KICK=0, TILT=0
MHE0L01V:  VKICK, L=0, KICK=0, TILT=0
MHB0L01AH: HKICK, L=0, KICK=0, TILT=0
MHB0L01AV: VKICK, L=0, KICK=0, TILT=0
MBH0L02H:  HKICK, L=0, KICK=0, TILT=0
MBH0L02V:  VKICK, L=0, KICK=0, TILT=0
MHE0L02H:  HKICK, L=0, KICK=0, TILT=0 ! does not appear in any line definition

! QUADRUPOLES

MQS0L01:   KQUAD, L=0.08, K1=0.0, TILT=0.785398163397448
MQS0L01A:  KQUAD, L=0.08, K1=0.0, TILT=0.785398163397448
MQS0L02:   KQUAD, L=0.08, K1=0.0, TILT=0.785398163397448

! taken from 10 MeV/c to 5D line all-save: 01=-0.5514, 02=-5.3872, 02A=5.5193
MQJ0L01:  KQUAD,L=0.15,K1=  2.9248
MQJ0L02:  KQUAD,L=0.15,K1=-11.13273193491806
MQJ0L02A: KQUAD,L=0.15,K1=  7.67444177359485

! TYPICALLY MEASURE TWISS PARAMETERS AT THE ENTRANCE TO MQJ0L02 

!BUBBLE_A:  LINE=(D0,ITV0L01,D1,MBH0L01H,MBH0L01V,D2,IPM0L01,D3,MQS0L01,D4,MQJ0L01,D5,MHE0L01V,D6,MQS0L01A,D7,      &
!                 MHB0L01AH,MHB0L01AV,D8,ITV0L02,MQJ0L02,D9,IPM0L02,D3,MQS0L02,D10,MBH0L02H,MBH0L02V,D11,MQJ0L02A,D12)

BUBBLE_A:  LINE=(MQJ0L02,D9,IPM0L02,D3,MQS0L02,D10,MBH0L02H,MBH0L02V,D11,MQJ0L02A,D12)

!---------------------------------------------------------------
!     Region 1: 5D LINE
!---------------------------------------------------------------

! Used drawing "ACC2008000-1100.pdf"
! S = 0 starts from starting point of 5 MeV dipole 
! 5 MeV bend, 25 degrees into 5D line

! 5 MeV DIPOLE
! bend radius of 0.2410 m and 25 degrees

MDL0L02:    CSBEND, L=0.131947, ANGLE="25.0 cdtor * ", E1=" 0.0 cdtor * ", E2="25.0 cdtor * ", &
            EDGE_ORDER=2, HGAP=0.013564, FINT=0.5, NONLINEAR=1, N_KICKS=25, INTEGRATION_ORDER=4

! DRIFTS

D5D00:      DRIFT, L=0.1747 ! subtract off length of straight line from survey center to dipole exit
D5D01:      DRIFT, L=0.1747
D5D02:      DRIFT, L=0.1907
D5D03:      DRIFT, L=0.1305
D5D04:      DRIFT, L=1.1165
D5D05:      DRIFT, L=0.1185
D5D06:      DRIFT, L=0.3272
D5D07:      DRIFT, L=0.0000
D5D08:      DRIFT, L=0.6668
D5D09:      DRIFT, L=0.1138
D5D10:      DRIFT, L=0.2639
D5D11:      DRIFT, L=0.7468

! QUADRUPOLES (K1 > 0 --> focus horizontally)

MQD5D00: KQUAD,L=0.15,K1=-0.38350060247329
MQD5D01: KQUAD,L=0.15,K1= 4.31511965815579

! 1.00 mm round spot
! 5D00=-2.94, 5D01=0.85
! 0.75 mm round spot
! 5D00= 2.61, 5D01=0.38
! 0.50 mm round spot
! 5D00=12.00, 5D01= 6.40

! VIEWERS/SUPERHARPS/BPMS/FARADAY CP

ITV5D00:    WATCH, FILENAME="%s.TV5D00", MODE=COORD
ITV5D01:    WATCH, FILENAME="%s.TV5D01", MODE=COORD
IPM5D00:    MONITOR, L=0
IPM5D01:    MONITOR, L=0

RADIATOR:   MARK, FITPOINT=1
COLLIMATOR: MARK, FITPOINT=1
CHAMBER:    MARK, FITPOINT=1
DUMP:       MARK, FITPOINT=1

! CORRECTORS

MBH5D00H:   HKICK, L=0, KICK=0, TILT=0
MBH5D00V:   VKICK, L=0, KICK=0, TILT=0
MBH5D00AH:  HKICK, L=0, KICK=0, TILT=0
MBH5D00AV:  VKICK, L=0, KICK=0, TILT=0
MBH5D01H:   HKICK, L=0, KICK=0, TILT=0
MBH5D01V:   VKICK, L=0, KICK=0, TILT=0
MBH5D01AH:  HKICK, L=0, KICK=0, TILT=0
MBH5D01AV:  VKICK, L=0, KICK=0, TILT=0

BUBBLE_B:   LINE=(MDL0L02,D5D00,MBH5D00H,MBH5D00V,D5D01,ITV5D00,D5D02,MQD5D00,D5D03,MBH5D00AH,MBH5D00AV,IPM5D00,D5D04,MQD5D01,D5D05,&
                  MBH5D01H,MBH5D01V,IPM5D01,D5D06,ITV5D01,D5D07,MBH5D01AH,MBH5D01AV,D5D08,RADIATOR,D5D09,COLLIMATOR,D5D10,CHAMBER,  &
				  D5D11,DUMP)

!---------------------------------------------------------------
!     LINE DEFINITIONS
!---------------------------------------------------------------

BUBBLE_CHAMBER: LINE=(BUBBLE_A, BUBBLE_B)