! *** OPTICS=AD_ACCEL_HE-13MAR23 *** ! LCLS2sc matching routines ! ============================================================================== ! Modification History ! ------------------------------------------------------------------------------ ! 10-MAR-2021, M. Woodley ! * add LEI matching routines (search for "LEI matching") ! ------------------------------------------------------------------------------ ! 31-OCT-2018, Y. Nosochkov ! * update DASEL matching routines ! ------------------------------------------------------------------------------ ! 28-NOV-2016, M. Woodley ! * merge LCLS-II' version 28NOV16 with LCLS2sc version 04NOV16 ! * added definition of TBKYDAS0 parameter ... see SUBROUTINE MKICKDASEL ! 23-NOV-2016, Y. Nosochkov ! * add MKICKDASEL routine for setting DASEL kicker angle ! ------------------------------------------------------------------------------ ! 04-NOV-2016, Y. Nosochkov ! * add DASEL matching routines ! ------------------------------------------------------------------------------ ! 09-FEB-2016, M. Woodley ! * add 2nd-order dispersion matching for dogleg, SPDh, SPDs, and DL2B ! ------------------------------------------------------------------------------ ! 25-SEP-2015, Y. Nosochkov ! * update MDL2Mb routine, add new MDL2M2 routine ! 12-MAR-2015, Y. Nosochkov ! * add routines for A-line match and CUSXR match ! 03-APR-2014, M. Woodley ! * add missing "USE" command to MDLBXY routine ! ------------------------------------------------------------------------------ OPTION, -ECHO Dmdw : DRIF, L=1 TBKRDAS0 := 0 TBKRSPX0 := 0 MUXTCXI := 1.75 !TCAV-OTR phase advance (see MDGITCAV) KmaxQcm := 24.40/(Cb*Emax*Lcmq) !max |K1| @ Emax for a CM quadrupole KmaxQR := 47.31/(Cb*Emax*LQR) !max |K1| @ 165 A for 2Q10 KmaxQA := 150.00/(Cb*Emax*LQA) !max |K1| @ Emax for 1.26Q12 KmaxQM := 27.00/(Cb*Emax*LQM) !max |K1| @ Emax for 2Q4W KmaxQz := 63.80/(Cb*Emax*LQz) !max |K1| @ Emax for SLTR QF1545 KmaxQW := 16.50/(Cb*Emax*LQW) !max |K1| @ Emax for a QW KmaxQx := 20.20/(Cb*Emax*LQx) !max |K1| @ Emax for a 1.26Q3.5 KmaxQE := 100.00/(Cb*Emax*LQE) !max |K1| @ Emax for a QE KmaxQu := 40.00/(Cb*Emax*LQu) !max |K1| @ Emax for an undulator quadrupole KmaxQP := 38.58/(Cb*Emax*LQP) !max |K1| @ 150 A for 3.94Q17 KmaxQxi := 20.20/(Cb*0.15*LQx) !max |K1| @ 150 MeV for 1.26Q3.5 KmaxQMi := 6.75/(Cb*0.15*LQM) !max |K1| @ 150 MeV for 2Q4 KmaxQx1 := 20.20/(Cb*0.26*LQx) !max |K1| @ 260 MeV for 1.26Q3.5 KmaxQx2 := 20.20/(Cb*1.70*LQx) !max |K1| @ 1.7 GeV for 1.26Q3.5 KmaxSb := 870.00/(Cb*Emax*LSb) !max |K2| @ Emax for 1.38S3.00 ! ============================================================================== ! Injector/Linac matching (through bypass line) ! ------------------------------------------------------------------------------ ML0E : SUBROUTINE !no matching ... just check that final energy is 100 MeV BEAM, ENERGY=E0 USE, (GUN,L0) SAVEBETA, TWSSm, #E TWISS, BETA0=TWSS0, SAVE VALUE, TWSSm[ENERGY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, -SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML1E : SUBROUTINE VALUE, FracL1,GradL1 BEAM, ENERGY=Ei USE, L1 MATCH, BETX=1, BETY=1 VARY, FracL1, STEP=1.E-6 !VARY, GradL1, STEP=1.E-6 WEIGHT, ENERGY=1 CONSTR, #E, ENERGY=E1 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, E1,FracL1,GradL1,GradL1/FracL1,100*(1-FracL1) BEAM, ENERGY=Ei USE, L1 TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML2E : SUBROUTINE VALUE, FracL2,GradL2 BEAM, ENERGY=E1 USE, L2 MATCH, BETX=1, BETY=1 VARY, FracL2, STEP=1.E-6 !VARY, GradL2, STEP=1.E-6 WEIGHT, ENERGY=1 CONSTR, #E, ENERGY=E2 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, E2,FracL2,GradL2,GradL2/FracL2,100*(1-FracL2) BEAM, ENERGY=E1 USE, L2 TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.T., FILE="L2" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML3E : SUBROUTINE BEAM, ENERGY=E2 USE, L3 SAVEBETA, TWm, ENDL3B TWISS, BETX=1, BETY=1, SAVE SHOW, TWm VALUE, E3 PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.T., FILE="L3" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML4E : SUBROUTINE VALUE, FracL4,GradL4 BEAM, ENERGY=E3 USE, L4 MATCH, BETX=1, BETY=1 VARY, FracL4, STEP=1.E-6 !VARY, GradL4, STEP=1.E-6 WEIGHT, ENERGY=1 CONSTR, #E, ENERGY=Ef !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, Ef,FracL4,GradL4,GradL4/FracL4,100*(1-FracL4) BEAM, ENERGY=E3 USE, L4 TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.T., FILE="L4" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCMfodo1 : SUBROUTINE VALUE, KQCMfodo1 BEAM, ENERGY=1 USE, CMfodo1 CELL VARY, KQCMfodo1, STEP=1.E-6 CONSTR, #E, MUX=45/360, MUY=45/360 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQCMfodo1 !PRINT, FULL SAVEBETA, TWSSm, #S TWISS, SAVE SET, BXCMfodo1, TWSSm[BETX] SET, BYCMfodo1, TWSSm[BETY] VALUE, BXCMfodo1,BYCMfodo1 PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCMfodo2 : SUBROUTINE VALUE, KQCMfodo2 BEAM, ENERGY=1 USE, CMfodo2 CELL VARY, KQCMfodo2, STEP=1.E-6 CONSTR, #E, MUX=15/360, MUY=15/360 !30/30 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQCMfodo2 VALUE, KQCMfodo2*Lcmq !PRINT, FULL SAVEBETA, TWSSm, #S TWISS, SAVE SET, BXCMfodo2, TWSSm[BETX] SET, BYCMfodo2, TWSSm[BETY] VALUE, BXCMfodo2,BYCMfodo2 PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCMfodo3 : SUBROUTINE VALUE, KQCMfodo3 BEAM, ENERGY=1 USE, CMfodo3 CELL VARY, KQCMfodo3, STEP=1.E-6 CONSTR, #E, MUX=15/360, MUY=15/360 !30/30 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQCMfodo3 !PRINT, FULL SAVEBETA, TWSSm, #S TWISS, SAVE SET, BXCMfodo3, TWSSm[BETX] SET, BYCMfodo3, TWSSm[BETY] VALUE, BXCMfodo3,BYCMfodo3 PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MC0fodo : SUBROUTINE VALUE, KQColl0 BEAM, ENERGY=E1 USE, C0fodo CELL VARY, KQColl0, STEP=1.E-6 CONSTR, #E, MUX=45/360 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQColl0 USE, C0fodo SAVEBETA, TWSSm, #S TWISS, SAVE VALUE, TWSSm[BETX],TWSSm[BETY] USE, (3*C0fodo) !PRINT, FULL TWISS, BETA0=TWSSm, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MC1fodo : SUBROUTINE VALUE, KQColl1 BEAM, ENERGY=E1 USE, C1fodo CELL VARY, KQColl1, STEP=1.E-6 CONSTR, #E, MUX=45/360 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQColl1 USE, C1fodo SAVEBETA, TWSSm, #S TWISS, SAVE VALUE, TWSSm[BETX],TWSSm[BETY] USE, (3*C1fodo) !PRINT, FULL TWISS, BETA0=TWSSm, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MHEATER : SUBROUTINE !match with heater undulator OFF VALUE, ABh,ZhBBo BEAM, ENERGY=Ei USE, LSRHTR SET, K_und, 0 MATCH, BETX=1, BETY=1 VARY, ABh, STEP=1.E-6, LOWER=0 !VARY, ZhBBo, STEP=1.E-6 CONSTR, HTRUND, DX=0.075 !RMATRIX, #S/#E, RM(5,6)=R56HTR, WEIGHT(5,6)=1 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, ABh,ZhBBo,ABh*DEGRAD MATCH, BETX=1, BETY=1 RMATRIX, #S/#E, RM(5,6)=R56HTR, WEIGHT(5,6)=1 ENDMATCH !PRINT, FULL SURVEY TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MBC1 : SUBROUTINE VALUE, AB1 BEAM, ENERGY=E1 USE, BC1c MATCH, BETX=1, BETY=1, SURVEY VARY, AB1, STEP=1.E-6, UPPER=0 RMATRIX, #S/#E, RM(5,6)=R56BC1, WEIGHT(5,6)=1 !WEIGHT, XS=1 !CONSTR, BCX12B, XS=0.32 !maximum traverse (wall side) !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, AB1,AB1*DEGRAD !PRINT, FULL SURVEY TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MBC2 : SUBROUTINE VALUE, AB2 BEAM, ENERGY=E2 USE, BC2c MATCH, BETX=1, BETY=1, SURVEY VARY, AB2, STEP=1.E-6, LOWER=0 RMATRIX, #S/#E, RM(5,6)=R56BC2, WEIGHT(5,6)=1 !WEIGHT, XS=1 !CONSTR, BCX22B, XS=-0.630 !maximum traverse (aisle side) LMDIF, TOL=1.E-20, CALLS=10000 MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, AB2,AB2*DEGRAD !PRINT, FULL SURVEY TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCCDLU : SUBROUTINE VALUE, ABCXDLU0 BEAM, ENERGY=Ef USE, CCDLU SET, ABCXDLU, ABCXDLU0 MATCH, BETX=1, BETY=1 VARY, ABCXDLU, STEP=1.E-6, LOWER=0 RMATRIX, #S/#E, RM(5,6)=R56CCDLU, WEIGHT(5,6)=1 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, ABCXDLU !PRINT, FULL TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCCDLD : SUBROUTINE VALUE, ABCXDLD0 BEAM, ENERGY=Ef USE, CCDLD SET, ABCXDLD, ABCXDLD0 MATCH, BETX=1, BETY=1 VARY, ABCXDLD, STEP=1.E-6, LOWER=0 RMATRIX, #S/#E, RM(5,6)=R56CCDLD, WEIGHT(5,6)=1 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, ABCXDLD !PRINT, FULL TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MASTRA : SUBROUTINE VALUE, Lastra,DCAP0da[L],DCAP0db[L] BEAM, ENERGY=E0 USE, (GUN,L0) MATCH, SURVEY, ZS=Zgun VARY, Lastra, STEP=1.E-6 WEIGHT, ZS=1 CONSTR, ASTRA, ZS=Za !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, Lastra,DCAP0da[L],DCAP0db[L] BEAM, ENERGY=E0 USE, (GUN,L0) !PRINT, FULL SURVEY, X0=Xgun, Y0=Ygun, Z0=Zgun ENDSUBROUTINE ! ------------------------------------------------------------------------------ MTWSS0a : SUBROUTINE SET, QCM01[K1], 0 VALUE, BX0,AX0,BY0,AY0 BEAM, ENERGY=E0 USE, (GUN,L0) MATCH, BETA0=TWSS0 VARY, BX0, STEP=1.E-6, LOWER=0 VARY, AX0, STEP=1.E-6 CONSTR, ASTRA, BETX=BXa, ALFX=AXa !LMDIF, TOL=1.E-20, CALLS=10000 !SIMPLEX !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH SET, BY0, BX0 SET, AY0, AX0 BEAM, ENERGY=E0 USE, (GUN,L0) MATCH, BETA0=TWSS0 CONSTR, ASTRA, BETX=BXa, ALFX=AXa, BETY=BYa, ALFY=AYa ENDMATCH VALUE, BX0,AX0,BY0,AY0 BEAM, ENERGY=E0 USE, (GUN,L0) !PRINT, FULL SAVEBETA, TWSSm, BEAM0 !entrance to QCM01 TWISS, COUPLE, BETA0=TWSS0 VALUE, TWSSm[ENERGY],TWSSm[BETX],TWSSm[ALFX],TWSSm[BETY],TWSSm[ALFY] ENDSUBROUTINE ! ------------------------------------------------------------------------------ MTWSS0i : SUBROUTINE BEAM, ENERGY=E0 USE, (GUN,L0) SET, QCM01[K1], 0 VALUE, BX0,AX0,BY0,AY0 MATCH, BETA0=TWSS0 VARY, BX0, STEP=1.E-6, LOWER=0 VARY, AX0, STEP=1.E-6 VARY, BY0, STEP=1.E-6, LOWER=0 VARY, AY0, STEP=1.E-6 CONSTR, BEAM0, BETX=BXi, ALFX=AXi, BETY=BYi, ALFY=AYi !LMDIF, TOL=1.E-20, CALLS=10000 !SIMPLEX !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, BX0,AX0,BY0,AY0 BEAM, ENERGY=E0 USE, (GUN,L0) !PRINT, FULL SAVEBETA, TWSSi, BEAM0 SAVEBETA, TWSSm, ASTRA TWISS, COUPLE, BETA0=TWSS0 SHOW, TWSSi SHOW, TWSSm !VALUE, TWSSm[ENERGY],TWSSm[BETX],TWSSm[ALFX],TWSSm[BETY],TWSSm[ALFY] ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML0LH : SUBROUTINE VALUE, QCM01[K1],KQ0H01,KQ0H02,KQ0H03,KQ0H04 BEAM, ENERGY=Ei USE, (LCLS2scI,HTR) MATCH, BETA0=TWSSi VARY, QCM01[K1], STEP=1.E-6, LOWER=0 VARY, KQ0H01, STEP=1.E-6, UPPER=0 VARY, KQ0H02, STEP=1.E-6, LOWER=0 VARY, KQ0H03, STEP=1.E-6, LOWER=0 VARY, KQ0H04, STEP=1.E-6, UPPER=0 CONSTR, Q0H01[1], BETX10 CONSTR, OTRDG02, BETX=BMINDG0, ALFX=0, BETY=BMINDG0, ALFY=0 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQDG004,KQDG005,KQDG006,KQDG007,KQDG008,KQDG009 VALUE, KQDG010,KQDG011 MATCH, & BETX=TWSSm[BETX], ALFX=TWSSm[ALFX], BETY=TWSSm[BETY], ALFY=TWSSm[ALFY] VARY, KQDG010, STEP=1.E-6, LOWER=0 VARY, KQDG011, STEP=1.E-6, UPPER=0 CONSTR, OTRDG04, BETY<0.16, BETX<6 WEIGHT, MUX=1, BETX=0, ALFX=0, DX=0, DPX=0, & MUY=0, BETY=0, ALFY=0, DY=0, DPY=0, & T=0, PT=0 COUPLE, MTCXDG0/OTRDG04, MUX=270/360 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQDG010,KQDG011 !PRINT, FULL SURVEY, X0=Xgun, Y0=Ygun, Z0=Zdg0, TAPE="survey.tape" TWISS, COUPLE, SAVE, TAPE="twiss.tape", & BETX=TWSSm[BETX], ALFX=TWSSm[ALFX], BETY=TWSSm[BETY], ALFY=TWSSm[ALFY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=DX, VAXIS2=DY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MBPMDG000 : SUBROUTINE !match offset of kicked beam BEAM, ENERGY=Ei USE, (LCLS2scI,HTR) SAVEBETA, TWSSm, #E TWISS, COUPLE, BETA0=TWSSi BEAM, ENERGY=TWSSm[ENERGY] USE, DIAG0 VALUE, ZKSD0a MATCH, SURVEY VARY, ZKSD0a, STEP=1.E-6 WEIGHT, YS=1 CONSTR, BPMDG000, YS=0.005 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, ZKSD0a PRINT, FULL SURVEY, X0=Xgun, Y0=Ygun, Z0=Zdg0, TAPE="survey.tape" ENDSUBROUTINE ! ============================================================================== ! Bypass/LTU matching ! ------------------------------------------------------------------------------ ! set phase advance of FODO cell in bypass straight line TWSBYP : BETA0 MSTRT : SUBROUTINE USE, FODO VALUE, KQY CELL VARY, KQY, STEP=1.E-4, LOWER=0, UPPER=0.08*1.5 CONSTR, #E, MUX=45/360, MUY=45/360 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQY ! SAVEBETA, TWSBYP, #S ! TWISS ! VALUE, TWSBYP[BETX],TWSBYP[BETY],TWSBYP[ALFX],TWSBYP[ALFY] ENDSUBROUTINE ! ============================================================================== ! LEI matching ! ------------------------------------------------------------------------------ ML0LEIE : SUBROUTINE BEAM, ENERGY=E0LEI USE, (GUNLEI,L0LEI) VALUE, fG00 MATCH, BETX=1, BETY=1 VARY, fG00, STEP=1.E-6 WEIGHT, ENERGY=1 CONSTR, #E, ENERGY=Ei !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, fG00 BEAM, ENERGY=E0LEI USE, (GUNLEI,L0LEI) SAVEBETA, TWSSm, #E TWISS, BETA0=TWSS0LEI, SAVE VALUE, TWSSm[ENERGY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, -SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MTWSS0iLEI : SUBROUTINE !match LEI initial twiss BEAM, ENERGY=E0LEI USE, (GUNLEI,L0LEI) SET, QCM00[K1], 0 VALUE, BX0LEI,AX0LEI,BY0LEI,AY0LEI MATCH, BETA0=TWSS0LEI VARY, BX0LEI, STEP=1.E-6, LOWER=0 VARY, AX0LEI, STEP=1.E-6 VARY, BY0LEI, STEP=1.E-6, LOWER=0 VARY, AY0LEI, STEP=1.E-6 CONSTR, BEAM0LEI, BETX=BXiLEI, ALFX=AXiLEI, BETY=BYiLEI, ALFY=AYiLEI !SIMPLEX !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, BX0LEI,AX0LEI,BY0LEI,AY0LEI BEAM, ENERGY=E0LEI USE, (GUNLEI,L0LEI) !PRINT, FULL TWISS, COUPLE, BETA0=TWSS0LEI, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, -SPLINE, FILE="match", VMIN=0, VMAX=0.12 PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDLIX : SUBROUTINE !match DLI X-displacement BEAM, ENERGY=Ei USE, LEI_2 VALUE, LDLIB2B3 MATCH, BETA0=TWSSDLI, SURVEY VARY, LDLIB2B3, STEP=1.E-6 WEIGHT, XS=1 CONSTR, #E, XS=XDLI !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, LDLIB2B3 PRINT, FULL SURVEY ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDLIDX : SUBROUTINE !match DLI dispersion BEAM, ENERGY=Ei USE, LEI_2 VALUE, KQDLI1 MATCH, BETA0=TWSSDLI VARY, KQDLI1, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi CONSTR, CNTDLI, DX=0, DPX=0 CONSTR, BXDLI4b[1], DX=0, DPX=0 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQDLI1 TWISS, BETA0=TWSSDLI, SAVE, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE, FILE="match", RANGE=#S/#E !, VMIN=-0.5, VMAX=0.5 ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDLITW : SUBROUTINE !match DLI twiss BEAM, ENERGY=Ei USE, LEI_2 VALUE, BXDLI,AXDLI,BYDLI,AYDLI VALUE, KQDLI2,KQDLI3,KQDLI4 !VALUE, DDLI03[L],DDLI04[L],DDLI05[L],DDLI06[L] TWISS, BETA0=TWSSDLI, CHROM MATCH, BETA0=TWSSDLI VARY, BXDLI , STEP=1.E-3, LOWER=0.01, UPPER=100 VARY, BYDLI , STEP=1.E-3, LOWER=0.01, UPPER=100 VARY, AXDLI , STEP=1.E-4 VARY, AYDLI , STEP=1.E-4 VARY, KQDLI2, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDLI3, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDLI4, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi !VARY, DDLI04[L], STEP=1.E-4, LOWER=0.5 !VARY, DDLI05[L], STEP=1.E-4, LOWER=0.5 CONSTR, #S, BETX=BXDLI, BETY=BYDLI, ALFX=AXDLI, ALFY=AYDLI CONSTR, CNTDLI, ALFX=0, ALFY=0 WEIGHT, DDX=10 CONSTR, CNTDLI, DDX=0 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, #S , BETX<15, BETY<15 CONSTR, QDLI1[1], BETX<15, BETY<15 CONSTR, QDLI2[1], BETX<15, BETY<15 CONSTR, QDLI3[1], BETX<15, BETY<15 CONSTR, QDLI4[1], BETX<15, BETY<15 !LMDIF, TOL=1.E-20, CALLS=100000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, BXDLI,AXDLI,BYDLI,AYDLI VALUE, KQDLI2,KQDLI3,KQDLI4 !VALUE, DDLI03[L],DDLI04[L],DDLI05[L],DDLI06[L] PRINT, FULL TWISS, BETA0=TWSSDLI, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E!, VMIN=0, VMAX=12 PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E!, VMIN=-0.3, VMAX=0.3 TWISS, BETA0=TWSSDLI, CHROM, SAVE, TAPE="chrom.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DDX,DDY, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E!, VMIN=-2, VMAX=2 TWISS, BETA0=TWSSDLI, SAVE, DELTAP=0.01, TAPE="twiss_dp1.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=X, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E PLOT, TABLE=TWISS, HAXIS=S, VAXIS=PX, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDGIFODO : SUBROUTINE !match DIAGI 90 deg FODO BEAM, ENERGY=Ei USE, FODODGI VALUE, KQDGI05 CELL VARY, KQDGI05, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi CONSTR, #E, MUX=90/360, MUY=90/360 !LMDIF, TOL=1.E-20, CALLS=2000 ENDMATCH VALUE, KQDGI05 TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, FILE="match", RANGE=#S/#E ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDGITW : SUBROUTINE !match DIAGI beta BEAM, ENERGY=Ei USE, DIAGIS VALUE, KQDGI01,KQDGI02,KQDGI03,KQDGI04,KQDGI05,KQDGI08 MATCH, BETA0=TWSSDLI VARY, KQDGI01, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI02, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI03, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI04, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI08, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi CONSTR, QDGI05[1], ALFX=0, ALFY=0 CONSTR, QDGI06[1], ALFX=0, ALFY=0 CONSTR, QDGI07[1], ALFX=0, ALFY=0 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QDGI02[1], BETX<5, BETY<14 CONSTR, QDGI03[1], BETX<5, BETY<5 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, BYDGISa[1], BETX<5, BETY<5 CONSTR, BYDGISb[1], BETX<5, BETY<5 CONSTR, OTRDGI1[1], BETX<15, BETY<0.549 !LMDIF, TOL=1.E-20, CALLS=100000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQDGI01,KQDGI02,KQDGI03,KQDGI04,KQDGI05,KQDGI08 PRINT, FULL TWISS, BETA0=TWSSDLI, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, FILE="match", & STYLE=100, SPLINE, RANGE=#S/#E ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML0LEIDLI : SUBROUTINE !match twiss from LEI CM to LEI dogleg BEAM, ENERGY=Ei USE, LCLS2scI2 VALUE, QCM00[K1],KQ0D01,KQ0D02,KQ0D03,KQ0D04 VALUE, D0D01[L],D0D02[L],D0D03[L],D0D04[L] MATCH, BETA0=TWSSiLEI VARY, QCM00[K1], STEP=1.E-5 VARY, KQ0D01, STEP=1.E-5, LOWER=-KmaxQMi, UPPER=+KmaxQMi VARY, KQ0D02, STEP=1.E-5, LOWER=-KmaxQMi, UPPER=+KmaxQMi VARY, KQ0D03, STEP=1.E-5, LOWER=-KmaxQMi, UPPER=+KmaxQMi VARY, KQ0D04, STEP=1.E-5, LOWER=-KmaxQMi, UPPER=+KmaxQMi !VARY, D0D01[L], STEP=1.E-4, LOWER=0.427898, UPPER=3.5 !VARY, D0D02[L], STEP=1.E-4, LOWER=1.2, UPPER=3.5 !VARY, D0D04[L], STEP=1.E-4, LOWER=0.8, UPPER=3.5 CONSTR, #E, BETX=BXDLI, ALFX=AXDLI, BETY=BYDLI, ALFY=AYDLI WEIGHT, BETX=1.E-5, BETY=1.E-5 !for nominal match !WEIGHT, BETX=1.E-6, BETY=1.E-6 !for Bmag sensitivity study CONSTR, Q0D01[1], BETX<46, BETY<46 CONSTR, Q0D02[1], BETX<40, BETY<41 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, Q0D03[1], BETX<25, BETY<25 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH !VALUE, LD0D14 VALUE, QCM00[K1],KQ0D01,KQ0D02,KQ0D03,KQ0D04 VALUE, D0D01[L],D0D02[L],D0D03[L],D0D04[L] USE, (LCLS2scI2,LEI_2) PRINT, FULL TWISS, COUPLE, BETA0=TWSSiLEI, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, RANGE=#S/#E, FILE="match" !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & ! STYLE=100, SPLINE, RANGE=BEGLEI_2/#E, FILE="match", & ! TITLE="LEI Dogleg" !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & ! STYLE=100, SPLINE, RANGE=BEGLEI_2/#E, FILE="match", & ! TITLE="LEI Dogleg" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDLILH : SUBROUTINE !match from LEI dogleg to laser heater BEAM, ENERGY=Ei USE, (LCLS2scI2,LEI_2,HTR_2) SETK2scLEI VALUE, KQ0H03,KQ0H04,KQ0H05,KQ0H06,KQ0H07,KQ0H08 MATCH, BETA0=TWSSiLEI VARY, KQ0H03, STEP=1.E-6, LOWER=-KmaxQMi, UPPER=+KmaxQMi VARY, KQ0H04, STEP=1.E-6, LOWER=-KmaxQMi, UPPER=+KmaxQMi VARY, KQ0H05, STEP=1.E-6, LOWER=-KmaxQMi, UPPER=0 VARY, KQ0H06, STEP=1.E-6, UPPER=+KmaxQMi, LOWER=0 !VARY, KQ0H07, STEP=1.E-6, UPPER=+KmaxQMi, LOWER=0 !VARY, KQ0H08, STEP=1.E-6, LOWER=-KmaxQMi, UPPER=0 !CONSTR, Q0H06[1], BETX<40 CONSTR, HTRUNDB, BETX=BXlh, ALFX=AXlh, BETY=BYlh, ALFY=AYlh CONSTR, #E, X=0, PX=0, Y=0, PY=0 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQ0H03,KQ0H04,KQ0H05,KQ0H06,KQ0H07,KQ0H08 PRINT, FULL SURVEY, X0=Xgun2, Y0=Ygun2, Z0=-12.51795552, TAPE="survey.tape" TWISS, BETA0=TWSSiLEI, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, -SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MDGITCAV : SUBROUTINE !match phase advance from TCAV center to OTRDGI1 BEAM, ENERGY=Ei USE, (LCLS2scI2,DIAGIS) SETK2scDGI VALUE, KQDGI01,KQDGI02,KQDGI03,KQDGI04,KQDGI05,KQDGI08 MATCH, BETA0=TWSSiLEI VARY, KQDGI01, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI02, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI03, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI04, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI05, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi VARY, KQDGI08, STEP=1.E-4, LOWER=-KmaxQxi, UPPER=KmaxQxi CONSTR, QDGI05[1], ALFX=0, ALFY=0 CONSTR, QDGI06[1], ALFX=0, ALFY=0 CONSTR, QDGI07[1], ALFX=0, ALFY=0 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QDGI02[1], BETX<5, BETY<14 CONSTR, QDGI03[1], BETX<5, BETY<5 CONSTR, BYDGISa[1], BETX<5, BETY<5 CONSTR, BYDGISb[1], BETX<5, BETY<5 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, OTRDGI1[1], BETX>8.5, BETY<0.5 WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0, & X=0, PX=0, Y=0, PY=0, T=0, PT=0, MUY=0 COUPLE, MTCXDGI[1]/OTRDGI1[1], MUX=MUXTCXI WEIGHT, MUX=1.E-4 COUPLE, MTCXDGI[1]/OTRDGI2[1], MUX=MUXTCXI-0.5 !LMDIF, TOL=1.E-20, CALLS=100000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQDGI01,KQDGI02,KQDGI03,KQDGI04,KQDGI05,KQDGI08 PRINT, FULL TWISS, COUPLE, BETA0=TWSSiLEI, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, RANGE=#S/#E, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, RANGE=DIAGIS, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE, RANGE=DIAGIS, FILE="match" ENDSUBROUTINE ! ============================================================================== ! Downstream matching (from bypass line) ! ------------------------------------------------------------------------------ ! ------------------------------------------------------------------------------ ! match XYZ in HXR spreader ! a general note: ! using matching variables defined with formulas in .xsif (e.g. =...*SETKIKh) ! will invalidate the formula -- therefore, subsequent matching routines using ! a different parameter setting in the formula (SETKIKh in this example) may ! not work correctly -- hence, after matching this type of variable, comment ! it out in the matching routine (or do not call this routine) and then ! proceed with the next matching routine MXYZSPh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (FODOLa,SPRDh,BSYH1,BSYH2) VALUE, ABRSPh,TBRSPh,dABXSPh,AROSP1h,AROSP2h VALUE, ABCSP1h,ABCSP2h,dLDSP4h,dLDSP9h,dLDSP17h MATCH, BETX=1, BETY=1, SURVEY, XS=Xbyp, YS=Ybyp, ZS=Zbyp, & THETAS=0, PHIS=0, PSIS=0 VARY, ABRSPh , STEP=1.E-5 VARY, TBRSPh , STEP=1.E-4 VARY, dABXSPh, STEP=1.E-10 VARY, AROSP1h, STEP=1.E-8 VARY, AROSP2h, STEP=1.E-8 VARY, ABCSP1h, STEP=1.E-6 VARY, ABCSP2h, STEP=1.E-6 VARY, dLDSP4h, STEP=1.E-4 VARY, dLDSP9h, STEP=1.E-4 VARY, dLDSP17h, STEP=1.E-6 WEIGHT, XS=1, YS=1, ZS=1, THETAS=1, PHIS=1, PSIS=1 CONSTR, BYSP2hb, YS=Ybyp, PHIS=0, PSIS=0 CONSTR, QSP3h[1], XS=Xbyp/2, THETAS=0, PSIS=0 CONSTR, BSYEND, XS=0, YS=0, ZS=3224.022425, & THETAS=0, PHIS=0, PSIS=0 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, ABRSPh,TBRSPh,dABXSPh,AROSP1h,AROSP2h VALUE, ABCSP1h,ABCSP2h,dLDSP4h,dLDSP9h,dLDSP17h ! PRINT, SPRDh ; PRINT, BSYEND ! SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match XYZ in SXR spreader MXYZSPs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (FODOL,BYPM,BSYLTUsc) MATCH, BETX=1, BETY=1, SURVEY, XS=Xbyp, YS=Ybyp, ZS=Zbyp, & THETAS=0, PHIS=0, PSIS=0 VARY, dABXSPs , STEP=1.E-9 !!VARY, ABLXSPs0, STEP=1.E-5 VARY, ABCSP1s , STEP=1.E-6 VARY, ABCSP2s , STEP=1.E-6 VARY, AROSP1s , STEP=1.E-8 VARY, dLDSP13s, STEP=1.E-6 VARY, dLDSP2sa, STEP=1.E-4 VARY, dDCYd , STEP=1.E-4 WEIGHT, XS=1, YS=1, ZS=1, THETAS=1, PHIS=1, PSIS=1 CONSTR, BXSP1sb, XS=Xbyp-0.50, YS=Ybyp, THETAS=0, PHIS=0, PSIS=0 CONSTR, BSYENDb, XS=Xbyp, YS=Ybyp, ZS=3224.022425, & THETAS=0, PHIS=0, PSIS=0 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, dABXSPs,ABLXSPs0,ABCSP1s,ABCSP2s,AROSP1s VALUE, dLDSP13s,dLDSP2sa,dDCYd ! PRINT, full;!SPRDs ; PRINT, BSYENDb ! SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match XYZ in BSY dumpline MXYZSPd : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAOFF USE, (FODOLa,SPRDd) VALUE, ABYSPd,dLDSP6d MATCH, BETX=1, BETY=1, SURVEY, XS=Xbyp, YS=Ybyp, ZS=Zbyp, & THETAS=0, PHIS=0, PSIS=0 VARY, ABYSPd , STEP=1.E-5 VARY, dLDSP6d, STEP=1.E-4 WEIGHT, XS=1, YS=1, ZS=1, THETAS=1, PHIS=1, PSIS=1 CONSTR, MUWALLd, YS=0, ZS=3207.258425 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, ABYSPd,dLDSP6d ! PRINT, SPRDd ! SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! empirical match of initial XYZ at marker LTUSPLIT in BSY system MXYZBSY : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (SPRDs,FODOLb,FODOLc,BYPM) MATCH, BETX=1, BETY=1, SURVEY, XS=Xf, YS=Yf, ZS=Zf, & THETAS=THETAf, PHIS=PHIf, PSIS=PSIf VARY, Yf, STEP=1.E-4 VARY, Zf, STEP=1.E-4 WEIGHT, XS=1.E-6, YS=1, ZS=1, THETAS=1.E-6, PHIS=1.E-6, PSIS=1.E-6 CONSTR, S100B, XS=Xf, YS=YOFF*COS(2*AVB), ZS=-YOFF*SIN(2*AVB), & THETAS=THETAf, PHIS=PHIf, PSIS=PSIf !LMDIF, TOL=1.E-20 ENDMATCH VALUE, Xf,Yf,Zf,THETAf,PHIf,PSIf ! PRINT, FULL ! SURVEY, X0=Xf, Y0=Yf, Z0=Zf, THETA0=THETAf, PHI0=PHIf, PSI0=PSIf ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! fine tune XYZ in LTU/undulator of HXR MXYZDL2h : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, BSYLCLS2scH MATCH, BETX=1, BETY=1, SURVEY, XS=Xf, YS=Yf, ZS=Zf, & THETAS=THETAf, PHIS=PHIf, PSIS=PSIf VARY, AB3P, STEP=1.E-5 VARY, dLDU0h, STEP=1.E-4 VARY, dDUE1a, STEP=1.E-4 WEIGHT, XS=1, ZS=1 CONSTR, UMAHXh13[1], XS=-1.25, ZS=516.830984-0.009 CONSTR, ENDHXR, XS=-1.25, ZS=670.486492 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, AB3P/RADDEG,dLDU0h,dDUE1a ! PRINT, BEGHXR/ENDHXR ! SURVEY, X0=Xf, Y0=Yf, Z0=Zf, THETA0=THETAf, PHI0=PHIf, PSI0=PSIf ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! fine tune XYZ in LTU/undulator of SXR MXYZDL2s : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, BSYLCLS2scS MATCH, BETX=1, BETY=1, SURVEY, XS=Xf, YS=Yf, ZS=Zf, & THETAS=THETAf, PHIS=PHIf, PSIS=PSIf VARY, AB3Pb, STEP=1.E-5 VARY, DDVB, STEP=1.E-6 VARY, dDMM, STEP=1.E-6 VARY, dLDU0s, STEP=1.E-6 WEIGHT, XS=1, YS=1, ZS=1 CONSTR, DBMARK36b, YS=-0.8953042661, ZS=376.3663895 CONSTR, SXRTERM, XS=1.25, ZS=669.586492 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, AB3Pb,DDVB,dDMM,dLDU0s ! PRINT, SXRTERM ! SURVEY, X0=Xf, Y0=Yf, Z0=Zf, THETA0=THETAf, PHI0=PHIf, PSI0=PSIf ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match XYZ in the main SXR dumpline in BSY coordinate system ! (note: this match is applicable to HXR dumpline as well) MXYZDMP : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, BSYLCLS2scS MATCH, BETX=1, BETY=1, SURVEY, XS=Xf, YS=Yf, ZS=Zf, & THETAS=THETAf, PHIS=PHIf, PSIS=PSIf VARY, ABdm , STEP=1.E-5 VARY, ARODMP1S, STEP=1.E-5 VARY, ARODMP2S, STEP=1.E-5 VARY, dDSB0e , STEP=1.E-4 VARY, LDMP1 , STEP=1.E-4 VARY, dDWSDUMP, STEP=1.E-4 WEIGHT, XS=1, YS=1, ZS=1, PSIS=1 CONSTR, BYD2BA, ZS=689.6517 CONSTR, MQDMPb, ZS=704.3025 CONSTR, DUMPFACEb, XS=0.9399, YS=-2.654105, ZS=715.9266, PSIS=0 !LMDIF, TOL=1.E-20, CALLS=2000 ENDMATCH VALUE, ABdm,ARODMP1S,ARODMP2S,dDSB0e,LDMP1,dDWSDUMP ! PRINT, SXRTERM ; PRINT, BYD2BA ; PRINT, MQDMPb ; PRINT, DUMPFACEb ! SURVEY, X0=Xf, Y0=Yf, Z0=Zf, THETA0=THETAf, PHI0=PHIf, PSI0=PSIf ! USE, BSYLCLS2scH ! PRINT, HXRTERM ; PRINT, BYD2A ; PRINT, MQDMP ; PRINT, DUMPFACE ! SURVEY, X0=Xf, Y0=Yf, Z0=Zf, THETA0=THETAf, PHI0=PHIf, PSI0=PSIf ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! fine tune SXR drift length after QBP32 to match with end of CUSXR MZCUSXR : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (FODOL,BYPM) MATCH, BETX=1, BETY=1, SURVEY, XS=Xbyp, YS=Ybyp, ZS=Zbyp, & THETAS=0, PHIS=0, PSIS=0 VARY, zMRGCUS, STEP=1.E-6 WEIGHT, ZS=1 CONSTR, MRGCUSXR, ZS=ZCuS2 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, zMRGCUS ! PRINT, MRGCUSXR ! SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match HXR spreader line TBXSPM := 120. TBYSPM := 60. MSPRDh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, SPRDh VALUE, dLDSP3h,dLDSP5h,dLDSP6h VALUE, TBXSP,TBYSP,TAXSP,TAYSP VALUE, KQSP1h,KQSP2h,KQSP3h,KQSP4h,KQSP5h,KQSP6h,KQSP7h VALUE, KQSP10h,KQSP11h,KQSP12h,KQSP13h MATCH, BETA0=TWSSP !!VARY, dLDSP3h, STEP=1.E-3, LOWER=-2, UPPER=2 !!VARY, dLDSP5h, STEP=1.E-3, LOWER=-7, UPPER=7 !!VARY, dLDSP6h, STEP=1.E-3, LOWER=-7, UPPER=7 !VARY, TBXSP , STEP=1.E-3, LOWER=0.01 !uncomment beta/alpha for R56 tuning !VARY, TBYSP , STEP=1.E-3, LOWER=0.01 !VARY, TAXSP , STEP=1.E-4 !VARY, TAYSP , STEP=1.E-4 VARY, KQSP1h , STEP=1.E-4, LOWER=0, UPPER= 0.570 VARY, KQSP2h , STEP=1.E-4, UPPER=0, LOWER=-2.807 VARY, KQSP3h , STEP=1.E-4, LOWER=0, UPPER= 1.405 !10 GeV VARY, KQSP4h , STEP=1.E-4, UPPER=0, LOWER=-0.570 VARY, KQSP5h , STEP=1.E-4, LOWER=0, UPPER= 0.570 VARY, KQSP6h , STEP=1.E-4, UPPER=0, LOWER=-0.570 VARY, KQSP7h , STEP=1.E-4, LOWER=0, UPPER= 2.807 VARY, KQSP10h, STEP=1.E-4, UPPER=0, LOWER=-0.570 VARY, KQSP11h, STEP=1.E-4, LOWER=0, UPPER= 1.405 !10 GeV VARY, KQSP12h, STEP=1.E-4, UPPER=0, LOWER=-2.807 VARY, KQSP13h, STEP=1.E-4, LOWER=0, UPPER= 0.570 RMATRIX, #S/#E, RM(5,6)=R56SPh, WEIGHT(5,6)=100 CONSTR, #E, DX=0, DPX=0, DY=0, DPY=0 !WEIGHT, BETX=1.E-4, BETY=1 !!WEIGHT, BETX=3.E-6, BETY=3.E-6 !for R56 tuning only !CONSTR, #S, BETX-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP4h[1] , BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP5h[1] , BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP6h[1] , BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP7h[1] , BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP8h[1] , BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP9h[1] , BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP10h[1], BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 CONSTR, QSP11h[1], BETX<110, BETY<110, DX>-0.35,DX<0.35, DY>-0.15,DY<0.15 WEIGHT, BETX=1.E-5, BETY=1.E-5 WEIGHT, ALFX=0,ALFY=0,DX=0,DPX=0,DY=0,DPY=0,X=0,PX=0,Y=0,PY=0,T=0,PT=0 COUPLE, QSP1h[1]/QSP13h[1] COUPLE, QSP2h[1]/QSP12h[1] COUPLE, QSP3h[1]/QSP11h[1] !LMDIF, TOL=1.E-20, CALLS=20000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, dLDSP3h,dLDSP5h,dLDSP6h VALUE, TBXSP,TBYSP,TAXSP,TAYSP VALUE, KQSP1h,KQSP2h,KQSP3h,KQSP4h,KQSP5h,KQSP6h,KQSP7h VALUE, KQSP10h,KQSP11h,KQSP12h,KQSP13h TWISS, BETA0=TWSSP, MUX=0, MUY=0, SAVE, COUPLE !, TAPE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, STYLE=100, SPLINE, & FILE="LCLS2sc", RANGE=#S/#E PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, STYLE=100, SPLINE, & FILE="LCLS2sc", RANGE=#S/#E ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match HXR spreader line 2nd-order dispersion MSPRDh2 : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, SPRDh !SET, KSSP1h, 0 ; SET, TSSP1h, 0 !SET, KSSP2h, 0 ; SET, TSSP2h, 0 VALUE, KSSP1h,KSSP2h,TSSP1h,TSSP2h TWISS, CHROM, BETA0=TWSSP MATCH, BETA0=TWSSP VARY, KSSP1h, STEP=1.E-6 VARY, KSSP2h, STEP=1.E-6 VARY, TSSP1h, STEP=1.E-6, LOWER=-30, UPPER=+30 VARY, TSSP2h, STEP=1.E-6, LOWER=-30, UPPER=+30 WEIGHT, DDX=1, DDPX=10, DDY=1, DDPY=10 CONSTR, #E, DDX=0, DDPX=0, DDY=0, DDPY=0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KSSP1h,KSSP2h,TSSP1h,TSSP2h !MATCH, BETA0=TWSSP ! RMATRIX, #S/#E, RM(6,6)=1 ! TMATRIX, #S/#E, TM(6,6,6)=0 !ENDMATCH !PRINT, FULL TWISS, BETA0=TWSSP !, TAPE="twiss.tape" TWISS, CHROM, BETA0=TWSSP, SAVE=CHROM !, TAPE="chrom.tape" ! PLOT, TABLE=CHROM, HAXIS=S, VAXIS=BETX,BETY, & ! STYLE=100, SPLINE=.T., FILE="LCLS2sc" ! PLOT, TABLE=CHROM, HAXIS=S, VAXIS=DX,DY, & ! STYLE=100, SPLINE=.F., FILE="LCLS2sc" PLOT, TABLE=CHROM, HAXIS=S, VAXIS=DDX,DDY, & STYLE=100, SPLINE=.F., FILE="LCLS2sc" PLOT, TABLE=CHROM, HAXIS=S, VAXIS=DDPX,DDPY, & STYLE=100, SPLINE=.F., FILE="LCLS2sc" PLOT, TABLE=CHROM, HAXIS=S, VAXIS=WX,WY, & STYLE=100, SPLINE=.F., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match SXR spreader line MSPRDs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, SPRDs VALUE, KQSP1s,KQSP2s,KQSP3s,KQSP4s,KQSP5s,KQSP8s,KQSP9s MATCH, BETA0=TWSSP !!VARY, dLDSP2s, STEP=1.E-3, LOWER=-0.5, UPPER=5 !!VARY, dLDSP3s, STEP=1.E-3, LOWER=-0.5, UPPER=5 !!VARY, dLDSP5s, STEP=1.E-3, LOWER=-3, UPPER=5 !!VARY, dLDSP6s, STEP=1.E-3, LOWER=-2, UPPER=5 !VARY, KQSP1 , STEP=1.E-4, LOWER=0, UPPER= 0.570 !VARY, KQSP2 , STEP=1.E-4, UPPER=0, LOWER=-0.570 VARY, KQSP1s, STEP=1.E-4, LOWER=0, UPPER= 1.405 !10 GeV VARY, KQSP2s, STEP=1.E-4, UPPER=0, LOWER=-1.405 !10 GeV VARY, KQSP3s, STEP=1.E-4, LOWER=0, UPPER= 1.405 !10 GeV VARY, KQSP4s, STEP=1.E-4, UPPER=0, LOWER=-0.570 VARY, KQSP5s, STEP=1.E-4, LOWER=0, UPPER= 0.570 VARY, KQSP8s, STEP=1.E-4, UPPER=0, LOWER=-1.405 !10 GeV VARY, KQSP9s, STEP=1.E-4, LOWER=0, UPPER= 1.405 !10 GeV RMATRIX, #S/#E, RM(5,6)=R56SPs, WEIGHT(5,6)=100 CONSTR, QSP5s[1], ALFX=0, ALFY=0, DPX=0 CONSTR, #E, DX=0, DPX=0 WEIGHT, DY=1.E-4, DPY=1.E-4 CONSTR, #E, DY=0, DPY=0 WEIGHT, BETX=1.E-5, BETY=1 CONSTR, BKYSP1sa[1], BETX<100, BETY5000, BETY>5000 !LMDIF, TOL=1.E-20, CALLS=2000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQSP1d,KQSP2d TWISS, BETA0=TWSSP, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., RANGE=#S/#E, FILE="LCLS2sc" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.T., RANGE=#S/#E, FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match from bypass to LTUSPLIT MBYPSP : SUBROUTINE BEAM, ENERGY=Ei USE, (LCLS2scI,LCLS2scC,FODOLa) VALUE, KQBP25,KQBP26,KQBP27,KQBP35,KQBP28 MATCH, BETA0=TWSSi VARY, KQBP25, STEP=1.E-4, LOWER=-0.161, UPPER=0.161 VARY, KQBP26, STEP=1.E-4, LOWER=-0.161, UPPER=0.161 VARY, KQBP27, STEP=1.E-4, LOWER=-0.161, UPPER=0.161 VARY, KQBP35, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 VARY, KQBP28, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 CONSTR, #E, BETX=TBXSP, BETY=TBYSP, ALFX=TAXSP, ALFY=TAYSP WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QBP26[1], BETX<460, BETY<460 CONSTR, QBP27[1], BETX<460, BETY<460 CONSTR, QBP35[1], BETX<400, BETY<400 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQBP25,KQBP26,KQBP27,KQBP35,KQBP28 TWISS, BETA0=TWSSi, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=QBP23[1]/#E, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match DL2 dispersion in SXR MDL2Bb : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, DL2sc VALUE, KQDLb MATCH, BETA0=TWSSms VARY, KQDLb, STEP=1.E-4, LOWER=0, UPPER=0.570 !10 GeV CONSTR, #E, DX=0, DPX=0 !!RMATRIX, #S/#E, RM(1,6)=0, RM(2,6)=0, WEIGHT(1,6)=100, WEIGHT(2,6)=100 TMATRIX, #S/#E, TM(5,6,6)=0, WEIGHT(5,6,6)=0 !LMDIF, TOL=1.E-20, CALLS=1000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQDLb TWISS, SAVE, BETA0=TWSSms PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match DL2 2nd-order dispersion in SXR MDL2Bb2 : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, DL2sc VALUE, KSDL1,KSDL2 TWISS, CHROM, BETA0=TWSSms MATCH, BETA0=TWSSms VARY, KSDL1, STEP=1.E-6 VARY, KSDL2, STEP=1.E-6 WEIGHT, DDX=1, DDPX=10 CONSTR, #E, DDX=0, DDPX=0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KSDL1,KSDL2 !MATCH, BETA0=TWSSms ! RMATRIX, #S/#E, RM(6,6)=1 ! TMATRIX, #S/#E, TM(6,6,6)=0 !ENDMATCH !PRINT, FULL TWISS, BETA0=TWSSms !, TAPE="twiss.tape" TWISS, CHROM, BETA0=TWSSms, SAVE=CHROM !, TAPE="chrom.tape" ! PLOT, TABLE=CHROM, HAXIS=S, VAXIS=BETX,BETY, & ! STYLE=100, SPLINE=.T., FILE="LCLS2sc" ! PLOT, TABLE=CHROM, HAXIS=S, VAXIS=DX,DY, & ! STYLE=100, SPLINE=.F., FILE="LCLS2sc" PLOT, TABLE=CHROM, HAXIS=S, VAXIS=DDX,DDY, & STYLE=100, SPLINE=.F., FILE="LCLS2sc" PLOT, TABLE=CHROM, HAXIS=S, VAXIS=DDPX,DDPY, & STYLE=100, SPLINE=.F., FILE="LCLS2sc" PLOT, TABLE=CHROM, HAXIS=S, VAXIS=WX,WY, & STYLE=100, SPLINE=.F., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! find initial twiss for DL2 cells in SXR MDL2Mb : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, DL2sc VALUE, MBETXs,MALFXs,MBETYs,MALFYs MATCH, BETA0=TWSSms VARY, MBETXs, STEP=1.E-4, LOWER=0.01 VARY, MALFXs, STEP=1.E-5, LOWER=-20, UPPER=20 VARY, MBETYs, STEP=1.E-4, LOWER=0.01 VARY, MALFYs, STEP=1.E-5, LOWER=-20, UPPER=20 CONSTR, QDL12[1], ALFX=0, ALFY=0 CONSTR, QDL13[1], ALFX=0, ALFY=0 CONSTR, QDL14[1], ALFX=0, ALFY=0 CONSTR, QDL15[1], ALFX=0, ALFY=0 CONSTR, QDL16[1], ALFX=0, ALFY=0 CONSTR, QDL17[1], ALFX=0, ALFY=0 CONSTR, QDL18[1], ALFX=0, ALFY=0 CONSTR, QDL19[1], ALFX=0, ALFY=0 !LMDIF, TOL=1.E-20, CALLS=2000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, MBETXs,MALFXs,MBETYs,MALFYs TWISS, SAVE, BETA0=TWSSms PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match from bypass to horizontal dogleg bends in SXR LTU ! for nominal machine SC -> SXR with CUSXR beamline OFF MSPLTUs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (SPRDs,FODOLb,FODOLc,BYPM,BSYLTUsc) VALUE, KQBP36,KQBP30,KQBP31,KQBP32,KQBP33,KQBP34,KQDBL1,KQDBL2 MATCH, BETA0=TWSSP !!VARY, zQBP36, STEP=1.E-3, LOWER=-6, UPPER=6 !!VARY, zQBP30, STEP=1.E-3, LOWER=-6, UPPER=6 !!VARY, zQBP31, STEP=1.E-3, LOWER=-6, UPPER=6 !!VARY, zQBP32, STEP=1.E-3, LOWER=-6, UPPER=6 !!VARY, zQBP33, STEP=1.E-3, LOWER=-6, UPPER=6 !!VARY, zQBP34, STEP=1.E-3, LOWER=-6, UPPER=0.1 !comment VARY, KQBP36, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 VARY, KQBP30, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 VARY, KQBP31, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 VARY, KQBP32, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 !endcomment VARY, KQBP33, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 !10 GeV VARY, KQBP34, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 !10 GeV !VARY, KQDBL1, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !10 GeV !VARY, KQDBL2, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !10 GeV CONSTR, DBMARK34B, BETX=MBETXs, BETY=MBETYs, ALFX=MALFXs, ALFY=MALFYs !comment WEIGHT, BETX=1.E-5, BETY=1.E-5 !!CONSTR, QBP36[1], BETX<80 , BETY<80 CONSTR, QBP30[1], BETX<125, BETY<125 CONSTR, QBP31[1], BETX<155, BETY<125 CONSTR, QBP32[1], BETX<125, BETY<125 CONSTR, QBP33[1], BETX<125, BETY<125 CONSTR, QBP34[1], BETX<155, BETY<125 !CONSTR, QBP30[1], BETX<115, BETY<115 !QBP31 moves u/s !CONSTR, QBP31[1], BETX<155, BETY<115 !CONSTR, QBP32[1], BETX<115, BETY<115 !CONSTR, QBP33[1], BETX<115, BETY<115 !CONSTR, QBP34[1], BETX<155, BETY<115 !CONSTR, QBP30[1], BETX<100, BETY<100 !QBP31 moves d/s !CONSTR, QBP31[1], BETX<155, BETY<100 !CONSTR, QBP32[1], BETX<100, BETY<100 !CONSTR, QBP33[1], BETX<100, BETY<100 !CONSTR, QBP34[1], BETX<155, BETY<100 !SIMPLEX, TOL=1.E-20, CALLS=1000 !endcomment !LMDIF, TOL=1.E-20 !, CALLS=20000 !MIGRAD, TOL=1.E-20 ENDMATCH MATCH, BETA0=TWSSP VARY, KQBP33, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 !10 GeV VARY, KQBP34, STEP=1.E-4, LOWER=-0.650, UPPER=0.650 !10 GeV VARY, KQDBL1, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !10 GeV VARY, KQDBL2, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !10 GeV CONSTR, DBMARK34B, BETX=MBETXs, BETY=MBETYs, ALFX=MALFXs, ALFY=MALFYs !LMDIF, TOL=1.E-20 !, CALLS=20000 !MIGRAD, TOL=1.E-20 ENDMATCH !!VALUE, zQBP36,zQBP30,zQBP31,zQBP32,zQBP33,zQBP34 VALUE, KQBP36,KQBP30,KQBP31,KQBP32,KQBP33,KQBP34,KQDBL1,KQDBL2 PRINT, FULL TWISS, BETA0=TWSSP, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match", RANGE=#S/QDL20[2] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=CNTSP2S/DBLDL21, & STYLE=100, SPLINE=.T., FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=#S/DBLDL21, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match dispersion in VBENDb in SXR MVBENDb : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, VBENDb VALUE, KQVBb MATCH, BETX=10, BETY=10, ENERGY=Ef VARY, KQVBb, STEP=1.E-4, LOWER=-2.807, UPPER=0 RMATRIX, #S/#E, RM(3,6)=0, RM(4,6)=0 TMATRIX, #S/#E, TM(5,6,6)=0, WEIGHT(5,6,6)=0 !LMDIF, TOL=1.E-20, CALLS=2000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQVBb TWISS, BETX=10, BETY=10, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match from DL2sc to VBENDb in SXR MDL2VB : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, LTUsc VALUE, KQDL20,KQDL21,KQDL22,KQVM3b MATCH, BETA0=TWSSms VARY, KQDL20, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQDL21, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQDL22, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQVM3b, STEP=1.E-4, LOWER=-1.405, UPPER=1.405 CONSTR, BPMVB2B, ALFY=0, ALFX=0, BETY<10 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QDL21[1], BETX<45, BETY<45 CONSTR, QDL22[1], BETX<45, BETY<45 CONSTR, QVM3B[1], BETX<45, BETY<45 CONSTR, QVM4B[1], BETX<45, BETY<45 !LMDIF, TOL=1.E-20, CALLS=5000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQDL20,KQDL21,KQDL22,KQVM3b TWISS, BETA0=TWSSms, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=QDL19[1]/QVM4B[2], & STYLE=100, SPLINE=.T., FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, RANGE=QDL19[1]/QVM4B[2], & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match FODO quad strength for 45-deg wire-to-wire phase adv. in SXR MECELL : SUBROUTINE USE, ECELLb VALUE, KQED2b CELL VARY, KQED2b, STEP=1.E-4, LOWER=0, UPPER=0.555 CONSTR, #E, MUX=45/360 !LMDIF, TOL=1.E-20, CALLS=1000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQED2b TWISS, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=#S/#E, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! find periodic beta's at wire scanners in LTU emit-diag section in SXR MED2P : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, EDSYSb VALUE, BXEDs,AXEDs,BYEDs,AYEDs MATCH, ENERGY=Ef, BETX=BXEDs, ALFX=AXEDs, BETY=BYEDs, ALFY=AYEDs VARY, BXEDs, STEP=1.E-4, LOWER=0.01, UPPER=100 VARY, AXEDs, STEP=1.E-5, LOWER=-20 , UPPER=20 VARY, BYEDs, STEP=1.E-4, LOWER=0.01, UPPER=100 VARY, AYEDs, STEP=1.E-5, LOWER=-20 , UPPER=20 CONSTR, WS31b, BETX=BXEDs, BETY=BYEDs CONSTR, WS32b, BETX=BXEDs, BETY=BYEDs CONSTR, WS33b, BETX=BXEDs, BETY=BYEDs CONSTR, WS34b, BETX=BXEDs, BETY=BYEDs !LMDIF, TOL=1.E-20, CALLS=1000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, BXEDs,AXEDs,BYEDs,AYEDs TWISS, SAVE, BETX=BXEDs, ALFX=AXEDs, BETY=BYEDs, ALFY=AYEDs PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=#S/#E, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match from VBENDb to emit-diag section in SXR MED2DL2b : SUBROUTINE BEAM, ENERGY=Ei SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, SC_SXRi VALUE, KQVM3b,KQVM4b,KQEM1b,KQEM2b,KQEM3b,KQEM4b MATCH, BETA0=TWSSi VARY, KQVM3b, STEP=1.E-4, LOWER=-1.405, UPPER=1.405 VARY, KQVM4b, STEP=1.E-4, LOWER=-1.405, UPPER=1.405 VARY, KQEM1b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQEM2b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQEM3b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQEM4b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 CONSTR, WS31B, BETX=BXEDs, BETY=BYEDs, ALFX=AXEDs, ALFY=AYEDs CONSTR, QEM2B[1], BETX<65, BETY<65 CONSTR, QEM3B[1], BETX<70, BETY<70 !LMDIF, TOL=1.E-20 !, CALLS=20000 !MIGRAD, TOL=1.E-20 !, CALLS=2000 ENDMATCH MATCH, BETA0=TWSSi VARY, KQEM1b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQEM2b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQEM3b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 VARY, KQEM4b, STEP=1.E-4, LOWER=-2.807, UPPER=2.807 CONSTR, WS31B, BETX=BXEDs, BETY=BYEDs, ALFX=AXEDs, ALFY=AYEDs !LMDIF, TOL=1.E-20 !, CALLS=20000 !MIGRAD, TOL=1.E-20 !, CALLS=2000 ENDMATCH VALUE, KQVM3b,KQVM4b,KQEM1b,KQEM2b,KQEM3b,KQEM4b TWISS, BETA0=TWSSi, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=QDL19[1]/WS34B, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match beta through SXR undulator MUNDs : SUBROUTINE BEAM, ENERGY=Ei SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, SC_SXRi VALUE, KQUM1b,KQUM2b,KQUM3b,KQUM4b MATCH, BETA0=TWSSi VARY, KQUM1b, STEP=1.E-4, UPPER= 2.807, LOWER=-2.807*0 !10 GeV VARY, KQUM2b, STEP=1.E-4, LOWER=-2.807, UPPER= 2.807*0 VARY, KQUM3b, STEP=1.E-4, UPPER= 2.807, LOWER=-2.807*0 VARY, KQUM4b, STEP=1.E-4, LOWER=-2.807, UPPER= 2.807*0 WEIGHT, DX=0, DPX=0, DY=0, DPY=0, X=0, PX=0, Y=0, PY=0, T=0, PT=0 WEIGHT, BETX=1.E-5, BETY=1.E-5, ALFX=0, ALFY=0 COUPLE, MUQS[1]/MUQS[3] COUPLE, MUQS[2]/MUQS[4] COUPLE, MUQS[1]/MUQS[5] COUPLE, MUQS[2]/MUQS[6] COUPLE, MUQS[1]/MUQS[7] COUPLE, MUQS[2]/MUQS[8] COUPLE, MUQS[1]/MUQS[9] COUPLE, MUQS[2]/MUQS[10] COUPLE, MUQS[1]/MUQS[11] COUPLE, MUQS[2]/MUQS[12] COUPLE, MUQS[1]/MUQS[13] COUPLE, MUQS[2]/MUQS[14] COUPLE, MUQS[1]/MUQS[15] COUPLE, MUQS[2]/MUQS[16] COUPLE, MUQS[1]/MUQS[17] COUPLE, MUQS[2]/MUQS[18] COUPLE, MUQS[1]/MUQS[19] COUPLE, MUQS[2]/MUQS[20] COUPLE, MUQS[1]/MUQS[21] COUPLE, MUQS[2]/MUQS[22] COUPLE, MUQS[1]/MUQS[23] COUPLE, MUQS[2]/MUQS[24] COUPLE, MUQS[1]/MUQS[25] COUPLE, MUQS[2]/MUQS[26] COUPLE, MUQS[1]/MUQS[27] COUPLE, MUQS[2]/MUQS[28] COUPLE, MUQS[1]/MUQS[29] COUPLE, MUQS[2]/MUQS[30] COUPLE, MUQS[1]/MUQS[31] COUPLE, MUQS[2]/MUQS[32] !WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QUM2B[1] , BETX<65, BETY<65 CONSTR, QUM3B[1] , BETX<65, BETY<65 CONSTR, QUM4B[1] , BETX<65, BETY<65 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQUM1b,KQUM2b,KQUM3b,KQUM4b TWISS, BETA0=TWSSi, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=WS31B/QSXh47[2], & STYLE=100, SPLINE, FILE="match", TITLE="UNDS (SC beam)" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match dumpline in SXR COMMENT !old routine -- do not use MDUMPs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (LTUsc,SXRUND,DUMPLINEb) MATCH, BETA0=TWSSms VARY, kQDSXm, STEP=1.E-4, LOWER=-0.999, UPPER=0.999 !15 GeV VARY, KQUE1b, STEP=1.E-4, LOWER=-0.380, UPPER=0.380 !15 GeV VARY, KQUE2b, STEP=1.E-4, LOWER=-0.570, UPPER=0.570 !10 GeV VARY, KQDMPb, STEP=1.E-4, LOWER=-0.182, UPPER=0.182 !17 GeV CONSTR, OTRDMPB, DY>0.46, BETY<0.6, BETX>100 WEIGHT, ALFY=1.E-5 CONSTR, OTRDMPB, ALFY=0 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QUE1B[1] , BETY<150, BETX<350 CONSTR, QUE2B[1] , BETY<150, BETX<350 CONSTR, QDMP1B[1], BETY<150, BETX<350 CONSTR, QDMP2B[1], BETY<150, BETX<350 WEIGHT, DY=10, BETX=1.E-6 CONSTR, DUMPFACEB, DY<0.4 , BETX<390 !!do not use LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, kQDSXm,KQUE1b,KQUE2b,KQDMPb TWISS, BETA0=TWSSms, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QSX44[1]/DUMPFACEB, & TITLE="SXR undulator to Dump, w/o coupling", FILE="LCLS2sc" TWISS, BETA0=TWSSms, SAVE=TWISS, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QSX44[1]/DUMPFACEB, & TITLE="SXR undulator to Dump, with coupling", FILE="LCLS2sc" ENDSUBROUTINE ENDCOMMENT MDUMPs2 : SUBROUTINE !matching TCAV per Yuantao's requirements BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (LTUsc,SXRUND,DUMPLINEb) VALUE, kQDSXm,KQUE1b,KQUE2b,KQDMPb MATCH, BETA0=TWSSms !VARY, kQDSXm, STEP=1.E-4, LOWER=-KmaxQu, UPPER=+KmaxQu !10 GeV VARY, KQUE1b, STEP=1.E-4, LOWER= 0, UPPER=+KmaxQR !10 GeV VARY, KQUE2b, STEP=1.E-4, LOWER=-KmaxQR, UPPER= 0 !10 GeV VARY, KQDMPb, STEP=1.E-4, LOWER=-KmaxQP, UPPER= 0 !10 GeV WEIGHT, BETX=1.E-5 CONSTR, MTCXB, BETX>90 !need large betx (close to 100 m) WEIGHT, BETX=1.E-6 CONSTR, MTCXB, BETX<120 !!WEIGHT, DY=10, BETY=1.E-5, BETX=1 !relaxed bety !!CONSTR, OTRDMPB, DY>0.46, BETY<0.6, BETX>20 !need large betx WEIGHT, DY=10, BETY=1, BETX=1 !relaxed bety CONSTR, OTRDMPB, DY>0.46, BETY<1.2, BETX>20 !need large betx WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 WEIGHT, X=0, PX=0, Y=0, PY=0, T=0, PT=0 WEIGHT, MUX=10, MUY=0 COUPLE, MTCXB/OTRDMPB, MUX>70/360, MUX<110/360 WEIGHT, MUX=1.E-4 COUPLE, MTCXB/OTRDMPB, MUX=90/360 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QUE1B[1] , BETY<150, BETX<350 CONSTR, QUE2B[1] , BETY<150, BETX<350 CONSTR, QDMP1B[1], BETY<150, BETX<350 CONSTR, QDMP2B[1], BETY<150, BETX<350 WEIGHT, DY=10, BETX=1.E-6 CONSTR, DUMPFACEB, DY<0.4 , BETX<350 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, kQDSXm,KQUE1b,KQUE2b,KQDMPb PRINT, MTCXB ; PRINT, OTRDMPB TWISS, BETA0=TWSSms, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE, RANGE=QSXh44[1]/DUMPFACEB, & TITLE="SXR undulator to Dump, w/o coupling", FILE="match" TWISS, BETA0=TWSSms, SAVE=TWISS, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE, RANGE=QSXh44[1]/DUMPFACEB, & TITLE="SXR undulator to Dump, with coupling", FILE="match" ENDSUBROUTINE COMMENT MDUMPs4 : SUBROUTINE !for matching vs energy and undulator gap BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (LTUsc,SXRUND,DUMPLINEb) VALUE, kQDSXm,KQUE1b,KQUE2b,KQDMPb MATCH, BETA0=TWSSms VARY, kQDSXm, STEP=1.E-4, LOWER=-1.4276, UPPER=1.4276 !40 kG @ 10 GeV VARY, KQUE1 , STEP=1.E-4, LOWER=-0.149*17/Eu, UPPER=0.149*17/Eu VARY, KQUE2 , STEP=1.E-4, LOWER=-0.149*17/Eu, UPPER=0.149*17/Eu VARY, KQDMPB, STEP=1.E-4, LOWER=-0.182*17/Eu, UPPER=0.182*17/Eu CONSTR, OTRDMPB, BETX>10, DY<0.5 !betx should not be too small WEIGHT, BETX=1.E-4 CONSTR, OTRDMPB, BETX<100 !reasonably large betx is ok WEIGHT, BETY=1.E-2 CONSTR, OTRDMPB, BETY<10 !prefer small bety WEIGHT, BETY=1.E-5 CONSTR, OTRDMPB, BETY<2 !very small bety is ok WEIGHT, BETX=1.E-4 CONSTR, MTCXB, BETX>80 !need large betx WEIGHT, BETX=1.E-5 CONSTR, MTCXB, BETX<150 !but not too large due to small aperture WEIGHT, BETX=1.E-6 CONSTR, MTCXB, BETX=100 !try to keep it around 100 m WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, DUMPFACEB, BETX<350, BETY<350 WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 WEIGHT, X=0, PX=0, Y=0, PY=0, T=0, PT=0, MUY=0 WEIGHT, MUX=1.E-2 COUPLE, MTCXB/OTRDMPB, MUX=90/360 !need close to 90 deg phase advance WEIGHT, MUX=10 COUPLE, MTCXB/OTRDMPB, MUX=90/360 !need close to 90 deg phase advance !COUPLE, MTCXB/OTRDMPB, MUX>70/360, MUX<110/360 WEIGHT, BETX=1.E0 CONSTR, QUE2B[1], BETX<250 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QUE1B[1] , BETY<150, BETX<350 !CONSTR, QUE2B[1] , BETY<150, BETX<350 CONSTR, QDMP1B[1], BETY<150, BETX<350 LMDIF, TOL=1.E-20, CALLS=50000 MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, kQDSXm,KQUE1b,KQUE2b,KQDMPb PRINT, MTCXB ; PRINT, OTRDMPB ; PRINT, DUMPFACEB TWISS, BETA0=TWSSms, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., RANGE=QSX44[1]/DUMPFACEB, & TITLE="SXR undulator to Dump, w/o coupling", FILE="match", & VMIN=0, VMAX=400 PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QSX44[1]/DUMPFACEB, & TITLE="SXR undulator to Dump, w/o coupling", FILE="match" !TWISS, BETA0=TWSSms, SAVE=TWISS, COUPLE !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & ! STYLE=100, SPLINE=.T., RANGE=QSX44[1]/DUMPFACEB, & ! TITLE="SXR undulator to Dump, with coupling", FILE="match" ENDSUBROUTINE ENDCOMMENT ! ============================================================================== ! HXR matching ! ------------------------------------------------------------------------------ MDL2B : SUBROUTINE !match DL2 dogleg dispersion in HXR BEAM, ENERGY=Ef USE, DL23 VALUE, KQDL MATCH, BETX=1, BETY=1 VARY, KQDL, STEP=1.E-4, UPPER=0.827, LOWER=0 RMATRIX, #S/#E, RM(1,6)=0, RM(2,6)=0 TMATRIX, #S/#E, TM(5,6,6)=0, WEIGHT(5,6,6)=0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQDL TWISS, BETX=1, BETY=1, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! note: not exact match here is due to slight asymmetry of the two bend ! positions relative to the quad (per existing LCLS-I deck) ! and focusing due to R56 chicane bends COMMENT !old -- find symmetric beta through LTU bend cells MDL2M : SUBROUTINE BEAM, ENERGY=Ef USE, DL23 MATCH, BETA0=TWSSmh VARY, MBETXh, STEP=1.E-4, UPPER=100, LOWER= 0.01 VARY, MALFXh, STEP=1.E-5, UPPER= 20, LOWER=-20 VARY, MBETYh, STEP=1.E-4, UPPER=100, LOWER= 0.01 VARY, MALFYh, STEP=1.E-5, UPPER= 20, LOWER=-20 CONSTR, #E, BETX=MBETXh, ALFX=-MALFXh, BETY=MBETYh, ALFY=-MALFYh CONSTR, BPMDL3, BETX<5 !!do not use LMDIF, TOL=1.E-20 ENDMATCH VALUE, MBETXh,MALFXh,MBETYh,MALFYh TWISS, BETA0=TWSSmh, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., FILE="LCLS2sc" ENDSUBROUTINE ENDCOMMENT MDL2M2 : SUBROUTINE !new -- find symmetric beta through LTU bend cells BEAM, ENERGY=Ef USE, (MM1,DOGLG2A,DOGLG2B,MM2) VALUE, MBETXh,MALFXh,MBETYh,MALFYh MATCH, BETA0=TWSSmh VARY, MBETXh, STEP=1.E-4, UPPER=100, LOWER=0.01 VARY, MBETYh, STEP=1.E-4, UPPER=100, LOWER=0.01 VARY, MALFXh, STEP=1.E-5, UPPER=20 , LOWER=-20 VARY, MALFYh, STEP=1.E-5, UPPER=20 , LOWER=-20 WEIGHT, ALFX=5.E-2, ALFY=2.E-4 CONSTR, QDL31[1], ALFX=0, ALFY=0, BETX<5 CONSTR, QDL32[1], ALFX=0, ALFY=0, BETX<5 CONSTR, QDL33[1], ALFX=0, ALFY=0, BETX<5 CONSTR, QT12[1] , ALFX=0, ALFY=0 CONSTR, QT22[1] , ALFX=0, ALFY=0 CONSTR, QT32[1] , ALFX=0, ALFY=0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, MBETXh,MALFXh,MBETYh,MALFYh TWISS, BETA0=TWSSmh, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match from end of HXR spreader through BSY to LTU ! note two matching blocks in this routine MBSY : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (SPRDh,BSYH1,BSYH2,LTUH), RANGE=#S/MM1 VALUE, KQ50Q3,KQ4,KQ5,KQ6,KQA0,KQVM1,KQVM2,KQVM3,KQVM4 MATCH, BETA0=TWSSP !!VARY, zQ6 , STEP=1.E-4, LOWER=-5, UPPER=5 VARY, KQ50Q3, STEP=1.E-4 VARY, KQ4 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 VARY, KQ5 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 !17 GeV VARY, KQ6 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 VARY, KQA0 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 !VARY, KQVM1 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 !VARY, KQVM2 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 !VARY, KQVM3 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 !VARY, KQVM4 , STEP=1.E-4, LOWER=-0.949, UPPER=0.949 CONSTR, MM1, BETX=MBETXh, BETY=MBETYh, ALFX=MALFXh, ALFY=MALFYh WEIGHT, ALFX=2.E-5, ALFY=2.E-5 CONSTR, QVB2[1], ALFX=0, ALFY=0 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, Q4[1] , BETX<105, BETY<105 CONSTR, Q5[1] , BETX<120, BETY<105 CONSTR, Q6[1] , BETX<105, BETY<105 CONSTR, QA0[1] , BETX<105, BETY<105 CONSTR, QVM1[1], BETX<100, BETY<100 CONSTR, QVM2[1], BETX<100, BETY<100 CONSTR, QVB1[1], BETX<100, BETY<100 CONSTR, QVB2[1], BETX<100, BETY<100 CONSTR, QVB3[1], BETX<100, BETY<100 CONSTR, QVM3[1], BETX<100, BETY<100 SIMPLEX, TOL=1.E-20, CALLS=5000 !LMDIF, TOL=1.E-20, CALLS=20000 !MIGRAD, TOL=1.E-20, CALLS=5000 ENDMATCH MATCH, BETA0=TWSSP VARY, KQVM1, STEP=1.E-6, LOWER=-0.949, UPPER=0.949 VARY, KQVM2, STEP=1.E-6, LOWER=-0.949, UPPER=0.949 VARY, KQVM3, STEP=1.E-6, LOWER=-0.949, UPPER=0.949 VARY, KQVM4, STEP=1.E-6, LOWER=-0.949, UPPER=0.949 CONSTR, MM1, BETX=MBETXh, ALFX=MALFXh, BETY=MBETYh, ALFY=MALFYh !LMDIF, TOL=1.E-20 !, CALLS=5000 !MIGRAD, TOL=1.E-20 !, CALLS=5000 ENDMATCH !!VALUE, zQ6 VALUE, KQ50Q3,KQ4,KQ5,KQ6,KQA0,KQVM1,KQVM2,KQVM3,KQVM4 PRINT, FULL TWISS, BETA0=TWSSP, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=BSYH1/MM1, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=#S/MM1, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, RANGE=#S/MM1, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, STYLE=100, SPLINE=.T., & FILE="match", RANGE=LTU PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, STYLE=100, SPLINE=.T., & FILE="match", RANGE=LTU ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! not used any more COMMENT MBSYDL2 : SUBROUTINE !match BSY to DL2 dogleg BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (SPRDh,BSYH1,BSYH2,LTUH) MATCH, BETA0=TWSSP VARY, KQVM1, STEP=1.E-4, LOWER=-0.949, UPPER=0.949 VARY, KQVM2, STEP=1.E-4, LOWER=-0.949, UPPER=0.949 VARY, KQVM3, STEP=1.E-4, LOWER=-0.949, UPPER=0.949 VARY, KQVM4, STEP=1.E-4, LOWER=-0.949, UPPER=0.949 CONSTR, MM1, BETX=MBETXh, ALFX=MALFXh, BETY=MBETYh, ALFY=MALFYh !!do not use LMDIF, TOL=1.E-20, CALLS=5000 ENDMATCH VALUE, KQVM1,KQVM2,KQVM3,KQVM4 TWISS, BETA0=TWSSP, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, STYLE=100, SPLINE=.T., & FILE="LCLS2sc", RANGE=LTU PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, STYLE=100, SPLINE=.T., & FILE="LCLS2sc", RANGE=LTU ENDSUBROUTINE ENDCOMMENT ! ------------------------------------------------------------------------------ ! find periodic beta's at wire scanners in LTU emit-diag section in HXR MED2Ph : SUBROUTINE BEAM, ENERGY=Ef USE, EDSYS VALUE, BXEDh,AXEDh,BYEDh,AYEDh MATCH, BETX=BXEDh, ALFX=AXEDh, BETY=BYEDh, ALFY=AYEDh VARY, BXEDh, STEP=1.E-4, LOWER=0.01, UPPER=100 VARY, AXEDh, STEP=1.E-5, LOWER=-20 , UPPER=20 VARY, BYEDh, STEP=1.E-4, LOWER=0.01, UPPER=100 VARY, AYEDh, STEP=1.E-5, LOWER=-20 , UPPER=20 CONSTR, WS31, BETX=BXEDh, BETY=BYEDh CONSTR, WS32, BETX=BXEDh, BETY=BYEDh CONSTR, WS33, BETX=BXEDh, BETY=BYEDh CONSTR, WS34, BETX=BXEDh, BETY=BYEDh !LMDIF, TOL=1.E-20, CALLS=1000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, BXEDh,AXEDh,BYEDh,AYEDh TWISS, SAVE, BETX=BXEDh, ALFX=AXEDh, BETY=BYEDh, ALFY=AYEDh PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=#S/#E, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match from bend cells to emit-diag section MED2DL2 : SUBROUTINE BEAM, ENERGY=Ei SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, SC_HXRi VALUE, KQEM1,KQEM2,KQEM3,KQEM4 MATCH, BETA0=TWSSi VARY, KQEM1, STEP=1.E-6, UPPER=0.827, LOWER=-0.827 !17 GeV VARY, KQEM2, STEP=1.E-6, UPPER=0.827, LOWER=-0.827 VARY, KQEM3, STEP=1.E-6, UPPER=0.827, LOWER=-0.827 VARY, KQEM4, STEP=1.E-6, UPPER=0.827, LOWER=-0.827 CONSTR, WS31, BETX=BXEDh, ALFX=AXEDh, BETY=BYEDh, ALFY=AYEDh !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQEM1,KQEM2,KQEM3,KQEM4 TWISS, BETA0=TWSSi, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, RANGE=LTU, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match beta in the HXR undulator MUNDh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (LTUH,HXRUND), RANGE=MM1/#E VALUE, KQUM1,KQUM2,KQUM3,KQUM4 MATCH, BETA0=TWSSmh VARY, KQUM1, STEP=1.E-4, UPPER=0.827, LOWER=-0.827 !17 GeV VARY, KQUM2, STEP=1.E-4, UPPER=0.827, LOWER=-0.827 VARY, KQUM3, STEP=1.E-4, UPPER=0.827, LOWER=-0.827 VARY, KQUM4, STEP=1.E-4, UPPER=0.827, LOWER=-0.827 WEIGHT, DX=0, DPX=0, DY=0, DPY=0, X=0, PX=0, Y=0, PY=0, T=0, PT=0 WEIGHT, BETX=1.E-5, BETY=1.E-5, ALFX=0, ALFY=0 COUPLE, MUQh[1]/MUQh[3] COUPLE, MUQh[2]/MUQh[4] COUPLE, MUQh[1]/MUQh[5] COUPLE, MUQh[2]/MUQh[6] COUPLE, MUQh[1]/MUQh[7] COUPLE, MUQh[2]/MUQh[8] COUPLE, MUQh[1]/MUQh[9] COUPLE, MUQh[2]/MUQh[10] COUPLE, MUQh[1]/MUQh[11] COUPLE, MUQh[2]/MUQh[12] COUPLE, MUQh[1]/MUQh[13] COUPLE, MUQh[2]/MUQh[14] COUPLE, MUQh[1]/MUQh[15] COUPLE, MUQh[2]/MUQh[16] COUPLE, MUQh[1]/MUQh[17] COUPLE, MUQh[2]/MUQh[18] COUPLE, MUQh[1]/MUQh[19] COUPLE, MUQh[2]/MUQh[20] COUPLE, MUQh[1]/MUQh[21] COUPLE, MUQh[2]/MUQh[22] COUPLE, MUQh[1]/MUQh[23] COUPLE, MUQh[2]/MUQh[24] COUPLE, MUQh[1]/MUQh[25] COUPLE, MUQh[2]/MUQh[26] COUPLE, MUQh[1]/MUQh[27] COUPLE, MUQh[2]/MUQh[28] COUPLE, MUQh[1]/MUQh[29] COUPLE, MUQh[2]/MUQh[30] COUPLE, MUQh[1]/MUQh[31] COUPLE, MUQh[2]/MUQh[32] COUPLE, MUQh[1]/MUQh[33] WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QUM2[1], BETX<50, BETY<70 CONSTR, QUM3[1], BETX<50, BETY<50 CONSTR, QUM4[1], BETX<50, BETY<50 !LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQUM1,KQUM2,KQUM3,KQUM4 TWISS, BETA0=TWSSmh, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=QE35[1]/QHXh46[2], & STYLE=100, SPLINE=.T., FILE="match", TITLE="UNDH (SC beam)" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match dumpline in HXR COMMENT !old -- do not use MDUMPh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (LTUH,HXRUND,DUMPLINE), RANGE=MM1/#E MATCH, BETA0=TWSSmh VARY, kQHX46, STEP=1.E-4, LOWER=-0.999, UPPER=0.999 !15 GeV VARY, KQUE1 , STEP=1.E-4, LOWER=-0.149, UPPER=0.149 !17 GeV VARY, KQUE2 , STEP=1.E-4, LOWER=-0.169, UPPER=0.169 !15 geV VARY, KQDMP , STEP=1.E-4, LOWER=-0.182, UPPER=0.182 !17 GeV CONSTR, OTRDMP, DY>0.46, BETY<0.6, BETX>100 WEIGHT, ALFY=1.E-5 CONSTR, OTRDMP, ALFY=0 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QUE1[1] , BETY<150, BETX<350 CONSTR, QUE2[1] , BETY<150, BETX<350 CONSTR, QDMP1[1], BETY<150, BETX<350 CONSTR, QDMP2[1], BETY<150, BETX<350 WEIGHT, DY=10, BETX=1.E-6 CONSTR, DUMPFACE, DY<0.4 , BETX<450 !!do not use LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, kQHX46,KQUE1,KQUE2,KQDMP TWISS, BETA0=TWSSmh, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHXh44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, w/o coupling", FILE="LCLS2sc" TWISS, BETA0=TWSSmh, SAVE=TWISS, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHXh44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, with coupling", FILE="LCLS2sc" ENDSUBROUTINE ENDCOMMENT COMMENT !outdated due to new VPU undulator -- use MDUMPh3 instead MDUMPh2 : SUBROUTINE !matching TCAV per Yuantao's requirements BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (LTUH,HXRUND,DUMPLINE), RANGE=MM1/#E MATCH, BETA0=TWSSmh VARY, kQHX46, STEP=1.E-4, LOWER=-0.999, UPPER=0.999 !15 GeV VARY, KQUE1 , STEP=1.E-4, LOWER=-0.149, UPPER=0.149 !17 GeV VARY, KQUE2 , STEP=1.E-4, LOWER=-0.169, UPPER=0.169 !15 geV VARY, KQDMP , STEP=1.E-4, LOWER=-0.182, UPPER=0.182 !17 GeV WEIGHT, DY=10, BETY=1, BETX=1 !!CONSTR, OTRDMP, DY>0.46, BETY<9, BETX>20 !relaxed bety, large betx CONSTR, OTRDMP, DY>0.46, BETY<5, BETX>20 !relaxed bety, large betx !!WEIGHT, DY=1.E-3, BETY=1, BETX=1 !relaxed Dy !!WEIGHT, DY=1.E-2, BETY=1.E-4, BETX=1 !relaxed bety !!CONSTR, OTRDMP, DY>0.46, BETY<0.6, BETX>20 WEIGHT, BETX=1.E-4 !!CONSTR, MTCX, BETX>80 CONSTR, MTCX, BETX>56 !need large betx (up to 100 m) WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 WEIGHT, X=0, PX=0, Y=0, PY=0, T=0, PT=0 WEIGHT, MUX=10, MUY=0 COUPLE, MTCX/OTRDMP, MUX>70/360, MUX<110/360 WEIGHT, MUX=1.E-4 COUPLE, MTCX/OTRDMP, MUX=90/360 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QUE1[1] , BETY<150, BETX<350 CONSTR, QUE2[1] , BETY<150, BETX<350 CONSTR, QDMP1[1], BETY<150, BETX<350 CONSTR, QDMP2[1], BETY<150, BETX<350 WEIGHT, DY=10, BETX=1.E-6 CONSTR, DUMPFACE, DY<0.4 , BETX<350 !!do not use LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, kQHX46,KQUE1,KQUE2,KQDMP PRINT, MTCX ; PRINT, OTRDMP TWISS, BETA0=TWSSmh, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHXh44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, w/o coupling", FILE="LCLS2sc" TWISS, BETA0=TWSSmh, SAVE=TWISS, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHXh44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, with coupling", FILE="LCLS2sc" ENDSUBROUTINE ENDCOMMENT COMMENT !use Yuri's new matching routine MDUMPh4 MDUMPh3 : SUBROUTINE !matching TCAV per Yuantao's requirements BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (LTUH,HXRUND,DUMPLINE), RANGE=MM1/#E VALUE, kQHX46,KQUE1,KQUE2,KQDMP MATCH, BETA0=TWSSmh VARY, kQHX46, STEP=1.E-4, LOWER=-0.994, UPPER=0.994 !15 GeV VARY, KQUE1 , STEP=1.E-4, LOWER=-0.149, UPPER=0.149 !17 GeV VARY, KQUE2 , STEP=1.E-4, LOWER=-0.169*0.95, UPPER=0.169*0.95 !15 GeV VARY, KQDMP , STEP=1.E-4, LOWER=-0.182, UPPER=0.182 !17 GeV !VARY, kQHX46, STEP=1.E-4, LOWER=-3.730, UPPER=3.730 !4 GeV !VARY, KQUE1 , STEP=1.E-4, LOWER=-0.636, UPPER=0.636 !4 GeV !VARY, KQUE2 , STEP=1.E-4, LOWER=-0.636, UPPER=0.636 !4 GeV !VARY, KQDMP , STEP=1.E-4, LOWER=-0.777, UPPER=0.777 !4 GeV WEIGHT, BETX=1.E-5 CONSTR, MTCX, BETX>90 !need large betx (but <100 m) WEIGHT, BETX=1.E-6 CONSTR, MTCX, BETX<120 !!WEIGHT, DY=10, BETY=1.E-5, BETX=1 !relaxed bety !!CONSTR, OTRDMP, DY>0.46, BETY<0.6, BETX>20 !need large betx WEIGHT, DY=10, BETY=1, BETX=1 !relaxed bety CONSTR, OTRDMP, DY>0.46, BETY<2.0, BETX>20 !need large betx WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 WEIGHT, X=0, PX=0, Y=0, PY=0, T=0, PT=0 WEIGHT, MUX=10, MUY=0 COUPLE, MTCX/OTRDMP, MUX>70/360, MUX<110/360 WEIGHT, MUX=1.E-4 COUPLE, MTCX/OTRDMP, MUX=90/360 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QUE1[1] , BETY<150, BETX<350 CONSTR, QUE2[1] , BETY<150, BETX<350 CONSTR, QDMP1[1], BETY<150, BETX<350 CONSTR, QDMP2[1], BETY<150, BETX<350 WEIGHT, DY=10, BETX=1.E-6 CONSTR, DUMPFACE, DY<0.4 , BETX<350 !LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, kQHX46,KQUE1,KQUE2,KQDMP PRINT, MTCX ; PRINT, OTRDMP TWISS, BETA0=TWSSmh, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHXh44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, w/o coupling", FILE="LCLS2sc" TWISS, BETA0=TWSSmh, SAVE=TWISS, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHXh44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, with coupling", FILE="LCLS2sc" ENDSUBROUTINE ENDCOMMENT MDUMPh4 : SUBROUTINE !for matching vs energy and undulator gap BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (LTUH,HXRUND,DUMPLINE), RANGE=MM1/#E !VALUE, kQHX46,KQUE1,KQUE2,KQDMP MATCH, BETA0=TWSSmh !!! VARY, kQHX46, STEP=1.E-4, LOWER=-0.877*17/Eu, UPPER=0.877*17/Eu VARY, KQUE1 , STEP=1.E-4, LOWER=-0.149*17/Eu, UPPER=0.149*17/Eu VARY, KQUE2 , STEP=1.E-4, LOWER=-0.149*17/Eu, UPPER=0.149*17/Eu ! ------------------------------------------------------------------------------ ! match KQDMP for DY=0 at OTR ! also REMOVE constraint on OTR BETY below because it cannot be matched !!VARY, KQDMP, STEP=1.E-4, LOWER=-0.182*17/Eu, UPPER=0.182*17/Eu !!CONSTR, OTRDMP, DY=0 ! ------------------------------------------------------------------------------ ! for fixed KQDMP the DY is also fixed, so no need for constraint on DY CONSTR, OTRDMP, BETX>10 !betx should not be too small WEIGHT, BETX=1.E-4 CONSTR, OTRDMP, BETX<100 !reasonably large betx is ok WEIGHT, BETY=1.E-2 CONSTR, OTRDMP, BETY<10 !prefer small bety WEIGHT, BETY=1.E-5 CONSTR, OTRDMP, BETY<2 !very small bety is ok WEIGHT, BETX=1.E-4 CONSTR, MTCX, BETX>80 !need large betx WEIGHT, BETX=1.E-5 CONSTR, MTCX, BETX<150 !but not too large due to small aperture WEIGHT, BETX=1.E-6 CONSTR, MTCX, BETX=100 !try to keep it around 100 m WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, DUMPFACE, BETX<350, BETY<350 WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 WEIGHT, X=0, PX=0, Y=0, PY=0, T=0, PT=0, MUY=0 WEIGHT, MUX=1.E-2 COUPLE, MTCX/OTRDMP, MUX=90/360 !need close to 90 deg phase advance WEIGHT, MUX=10 COUPLE, MTCX/OTRDMP, MUX>70/360, MUX<110/360 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QUE1[1] , BETY<150, BETX<350 CONSTR, QUE2[1] , BETY<150, BETX<350 CONSTR, QDMP1[1], BETY<150, BETX<350 !LMDIF, TOL=1.E-20, CALLS=50000 ENDMATCH VALUE, KQUM1,KQUM2,KQUM3,KQUM4 VALUE, kQHX46,KQUE1,KQUE2 !!,KQDMP PRINT, MTCX ; PRINT, OTRDMP ; PRINT, DUMPFACE TWISS, BETA0=TWSSmh, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE=.T., RANGE=QHX44[1]/DUMPFACE, & TITLE="HXR undulator to Dump, w/o coupling", FILE="match" !TWISS, BETA0=TWSSmh, SAVE=TWISS, COUPLE !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & ! STYLE=100, SPLINE=.T., RANGE=QHX44[1]/DUMPFACE, & ! TITLE="HXR undulator to Dump, with coupling", FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match Twiss in A-line for SC beam ! match the A-line geometry in "LCLS2cu_match.mad8" before this match MALINEsc : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALON SETK2scA USE, BSYLCLS2scA, RANGE=#S/ALWALL VALUE, KQ10,KQ11,KQ19,KQ20,KQ27,KQ28,KQ30,KQ38 MATCH, BETA0=TWSSP VARY, KQ10, STEP=1.E-4 VARY, KQ11, STEP=1.E-4 VARY, KQ19, STEP=1.E-4 VARY, KQ20, STEP=1.E-4 VARY, KQ27, STEP=1.E-4 VARY, KQ28, STEP=1.E-4 VARY, KQ30, STEP=1.E-4 VARY, KQ38, STEP=1.E-4 CONSTR, ALWALL[1], BETX=7.834394963E+03, ALFX=4.252254496E+02, & BETY=4.018221983E+02, ALFY=2.208989467E+01, & DX=0, DPX=0 WEIGHT, ALFX=1.E-5, ALFY=1.E-5 CONSTR, Q19[1], ALFX=0, ALFY=0 WEIGHT, BETX=1.E-6 CONSTR, Q19[1], BETX<300 CONSTR, Q20[1], BETX<200 !LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQ10,KQ11,KQ19,KQ20,KQ27,KQ28,KQ30,KQ38 ! PRINT, FULL TWISS, COUPLE, SAVE, BETA0=TWSSP PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=rBETX,rBETY, SPLINE, & STYLE=100, TITLE="A-line", FILE="LCLS2sc", RANGE=BSYA PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, SPLINE, & STYLE=100, TITLE="A-line", FILE="LCLS2sc", RANGE=BSYA ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match XYZ in DASEL -- for optics with the septum moved d/s MKICKDASEL : SUBROUTINE BEAM, ENERGY=Ef SET, TBKRDAS0, TBKRDAS SET, TBKRDAS, 0.0 SET, SETSP, SETSPd ; SET, SETDA, DAON USE, DASEL, RANGE=#S/BLRDASb[1] MATCH, BETX=1, BETY=1, SURVEY, XS=0, YS=0, ZS=0, THETAS=0, PHIS=0, PSIS=0 VARY, ABKRDAS0, STEP=1.E-6 VARY, ZDDASBLXa, STEP=1.E-4 WEIGHT, YS=1 CONSTR, BLRDASa[1], YS=15.E-3 CONSTR, BPMDAS[1], YS=5.E-3 !LMDIF, TOL=1.E-20, CALLS=1000 ENDMATCH VALUE, ABKRDAS0,ZDDASBLXa !PRINT, RANGE=#S/BLRDASb[1] !SURVEY SET, TBKRDAS, TBKRDAS0 ENDSUBROUTINE MXYZDASEL : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 1.0 !for survey -- turn RC corrector roll angle ON SETK2scDA USE, (FODOLa,SPDASELA) VALUE, TBKRDAS,ABLRDAS0,ABRDAS1,TBRDAS1,ABRDAS2,TBRDAS2 VALUE, ARODAS1,ARODAS2,ARODAS3,dLDDAS5,dLDDAS,LDA04c,DDAS14aa[L] MATCH, BETX=1, BETY=1, SURVEY, XS=Xbyp, YS=Ybyp, ZS=Zbyp, & THETAS=0, PHIS=0, PSIS=0 VARY, TBKRDAS , STEP=1.E-6 VARY, ABLRDAS0, STEP=1.E-5 VARY, ABRDAS1 , STEP=1.E-8 VARY, TBRDAS1 , STEP=1.E-4 VARY, ABRDAS2 , STEP=1.E-5 VARY, TBRDAS2 , STEP=1.E-4 VARY, ARODAS1 , STEP=1.E-8 VARY, ARODAS2 , STEP=1.E-8 VARY, ARODAS3 , STEP=1.E-6 VARY, dLDDAS5 , STEP=1.E-4 VARY, dLDDAS , STEP=1.E-4 VARY, LDA04c , STEP=1.E-4 VARY, DDAS14aa[L], STEP=1.E-4 WEIGHT, XS=1, YS=1, ZS=1, THETAS=10, PHIS=10, PSIS=10 CONSTR, RODAS1, PHIS=0, PSIS=0 CONSTR, BRDAS1a, ZS=3010.8 -12.8 CONSTR, RODAS2, PSIS=0 CONSTR, QDAS11[1], XS=0.40, THETAS=0 CONSTR, WOODDOOR, ZS=3050.512 CONSTR, Q10[1], XS=6.690152099E-01, THETAS=8.737596263E-03, & YS=9.570135331E-03, PHIS =1.249881637E-04, & ZS=3.181405369E+03, PSIS =1.030684491E-02 WEIGHT, XS=1.E-4, YS=1.E-4, ZS=1.E-4, THETAS=1.E-3, PHIS=1.E-3, PSIS=1.E-3 CONSTR, MARC[1], XS=2.023616674E+01, THETAS=4.274979549E-01, & YS=2.146605227E-01, PHIS =4.305364423E-03, & ZS=3.278586348E+03, PSIS =9.365417912E-03 !LMDIF, TOL=1.E-20 ENDMATCH VALUE, TBKRDAS,ABLRDAS0,ABRDAS1,TBRDAS1,ABRDAS2,TBRDAS2 VALUE, ARODAS1,ARODAS2,ARODAS3,dLDDAS5,dLDDAS,LDA04c,DDAS14aa[L] !PRINT, FULL !PRINT, SPDASELA !SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 !, TAPE !--- !USE, SPDASELA !in BSY coordinate system !PRINT, FULL !SURVEY, X0=Xf, Y0=Yf, Z0=Zf, THETA0=THETAf, PHI0=PHIf, PSI0=PSIf ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match DASEL Twiss -- for optics with the septum moved d/s MDASEL : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 0.0 !for twiss -- turn RC corrector roll angle OFF SETK2scDA USE, SPDASELA, RANGE=#S/ALWALL MATCH, BETA0=TWSSP !!VARY, dLDDAS4, STEP=1.E-3, LOWER=-2.5, UPPER=5 !!VARY, dLDDAS11, STEP=1.E-3, LOWER=-5, UPPER=3 !!VARY, dLDDAS20, STEP=1.E-3, LOWER=-0.5, UPPER=10 VARY, KQDAS1 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !4 GeV VARY, KQDAS2 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS11, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS12, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS13, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS14, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !!VARY, KQDAS15, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !!VARY, KQDAS16, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS17, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !could be varied if needed VARY, KQDAS18, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS19, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQ10, STEP=1.E-4 VARY, KQ11, STEP=1.E-4 VARY, KQ19, STEP=1.E-4 VARY, KQ20, STEP=1.E-4 VARY, KQ27, STEP=1.E-4 VARY, KQ28, STEP=1.E-4 VARY, KQ30, STEP=1.E-4 VARY, KQ38, STEP=1.E-4 CONSTR, ALWALL[1], BETX=7.834394963E+03, ALFX=4.252254496E+02, & BETY=4.018221983E+02, ALFY=2.208989467E+01, & DX=0, DPX=0, DY=0, DPY=0 CONSTR, QDAS15[1], ALFX=0, ALFY=0 WEIGHT, ALFY=1.E-4, ALFX=1.E-2 CONSTR, Q19[1], ALFX=0, ALFY=0 WEIGHT, DX=1.E-3 CONSTR, QDAS16[1] , DX<0.55 CONSTR, QDAS18a[1], DX<0.55 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QDAS11[1], BETX<65, BETY<65 CONSTR, QDAS12[1], BETX<65, BETY<65 CONSTR, QDAS13[1], BETX<65, BETY<65 CONSTR, QDAS14[1], BETX<65, BETY<65 CONSTR, QDAS15[1], BETX<65, BETY<65 CONSTR, QDAS16[1], BETX<65, BETY<65 CONSTR, QDAS17[1], BETX<65, BETY<65 CONSTR, QDAS18a[1], BETX<65, BETY<65 CONSTR, QDAS19[1], BETX<65, BETY<65 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, BRDAS2b[1], BETX<65, BETY<65 RMATRIX, BEGDASEL/ENDDASEL, RM(5,6)=0, WEIGHT(5,6)=10 !LMDIF, TOL=1.E-20, CALLS=100000 ENDMATCH !VALUE, LDDAS,dLDDAS4,dLDDAS11,dLDDAS20,dLDDAS12 VALUE, KQDAS1,KQDAS2,KQDAS11,KQDAS12,KQDAS13,KQDAS14 VALUE, KQDAS17,KQDAS18,KQDAS19!,KQDAS15,KQDAS16 VALUE, KQ10,KQ11,KQ19,KQ20,KQ27,KQ28,KQ30,KQ38 !--- !PRINT, FULL !PRINT, CLASS=QUADRUPOLE TWISS, BETA0=TWSSP, SAVE, COUPLE !, TAPE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=100, SPLINE, FILE="match", RANGE=BEGDASEL/ENDDASEL, & VMIN=0,-4, VMAX=500,1 PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=rBETX,rBETY, VAXIS2=DX,DY, & STYLE=100, SPLINE, FILE="match", RANGE=BEGDASEL/#E, & VMIN=0,-10, VMAX=120,2 !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & ! STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=BEGDASEL/PRDAS14, & ! VMIN=0,-3, VMAX=500,1 !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DY, & ! STYLE=1, COLOUR=100, SPLINE=F, FILE="LCLS2sc", RANGE=BEGDASEL/DDAS5[1] !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, STYLE=100, SPLINE, & ! FILE="LCLS2sc", RANGE=#S/ENDDASEL !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, STYLE=100, SPLINE, & ! FILE="LCLS2sc", RANGE=#S/ENDDASEL !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & ! STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=#S/Q10[2], & ! VMIN=0,-2.5, VMAX=500,0.5 !PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=rBETX,rBETY, VAXIS2=DX,DY, & ! STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=ENDDASEL/#E, & ! VMIN=0,-8, VMAX=160,2 !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=rBETX,rBETY, STYLE=100, SPLINE, & ! FILE="LCLS2sc", RANGE=ENDDASEL/#E !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, STYLE=100, SPLINE, & ! FILE="LCLS2sc", RANGE=ENDDASEL/#E !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=rBETX,rBETY, COLOUR=100, SPLINE, & ! FILE="LCLS2sc", RANGE=#S/#E !PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, COLOUR=100, SPLINE, & ! FILE="LCLS2sc", RANGE=#S/#E !--- !TWISS, BETA0=TWSSP, COUPLE, SAVE=SIZE !PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCx,SIZEx, STYLE=1, SPLINE=F, & ! COLOUR=100, FILE="LCLS2sc", RANGE=BEGDASEL/#E, VMIN=0, VMAX=160 !PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCy,SIZEy, STYLE=1, SPLINE=F, & ! COLOUR=100, FILE="LCLS2sc", RANGE=BEGDASEL/#E, VMIN=0, VMAX=30 !TABLE, NAME=SIZE, columns=s,BSCx,BSCy,SIZEx,SIZEy !SELECT, OPTICS, FULL !SELECT, OPTICS, CLASS=QUADRUPOLE; SELECT, OPTICS, CLASS=SBEND !OPTICS, BETA0=TWSSP, COLUMNS=KEYWORD,NAME,MUX,MUY,BETX,BETY,DX,DY !OPTICS, BETA0=TWSSP, COLUMNS=NAME,L,K0L,K1L ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! match DASEL Twiss -- stage-1S optics with dump in S30 COMMENT MDASELS1S : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 0.0 !for twiss -- turn RC corrector roll angle OFF SETK2scDA SET, QDAS1b[K1], 0.0; SET, QDAS2b[K1], 0.0 !!SET, QDAS1a[K1], 0.0; SET, QDAS2a[K1], 0.0 !!SET, QDAS1a[K1], KQDAS1*2; SET, QDAS2a[K1], KQDAS2*2 SET, KQDAS1, -0.77; SET, KQDAS2, 0.75 SET, KQDAS11, 0.0 USE, SPDASEL1S MATCH, BETA0=TWSSP !VARY, KQDAS1 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !4 GeV !VARY, KQDAS2 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !LMDIF, TOL=1.E-20, CALLS=100000 ENDMATCH VALUE, KQDAS1,KQDAS2,KQDAS11 !--- PRINT, DASEL1S TWISS, BETA0=TWSSP, SAVE, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=BEGDASEL/#E, & VMIN=0,-4, VMAX=500,1 !--- TWISS, BETA0=TWSSP, COUPLE, SAVE=SIZE PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCx,SIZEx, STYLE=1, SPLINE=F, & COLOUR=100, FILE="LCLS2sc", RANGE=BEGDASEL/#E, VMIN=0, VMAX=25 PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCy,SIZEy, STYLE=1, SPLINE=F, & COLOUR=100, FILE="LCLS2sc", RANGE=BEGDASEL/#E, VMIN=0, VMAX=10 !TABLE, NAME=SIZE, columns=s,BSCx,BSCy,SIZEx,SIZEy !--- SET, ROLLRC, 1.0 !for survey -- turn RC corrector roll angle ON USE, (FODOLa,SPDASEL1S) PRINT, DASEL1S SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 !, TAPE ENDSUBROUTINE ENDCOMMENT ! ------------------------------------------------------------------------------ ! match DASEL Twiss -- stage-1B optics with dump in BSY COMMENT MDASELS1B : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 0.0 !for twiss -- turn RC corrector roll angle OFF SETK2scDA SET, QDAS1b[K1], 0.0; SET, QDAS2b[K1], 0.0 !!SET, QDAS1a[K1], 0.0; SET, QDAS2a[K1], 0.0 !!SET, QDAS1a[K1], KQDAS1*2; SET, QDAS2a[K1], KQDAS2*2 SET, KQDAS1, -0.77; SET, KQDAS2, 0.75 SET, KQDAS11, 0.0; SET, KQDAS12, 0.0; SET, KQDAS13, 0.0 USE, SPDASEL1B MATCH, BETA0=TWSSP !VARY, KQDAS1 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !4 GeV !VARY, KQDAS2 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS11, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS12, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !LMDIF, TOL=1.E-20, CALLS=100000 ENDMATCH VALUE, KQDAS1,KQDAS2,KQDAS11,KQDAS12,KQDAS13 !--- !PRINT, FULL !PRINT, CLASS=QUADRUPOLE TWISS, BETA0=TWSSP, SAVE, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=BEGDASEL/#E, & VMIN=0,-4, VMAX=500,1 !--- TWISS, BETA0=TWSSP, COUPLE, SAVE=SIZE PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCx,SIZEx, STYLE=1, SPLINE=F, & COLOUR=100, FILE="LCLS2sc", RANGE=BEGDASEL/#E, VMIN=0, VMAX=25 PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCy,SIZEy, STYLE=1, SPLINE=F, & COLOUR=100, FILE="LCLS2sc", RANGE=BEGDASEL/#E, VMIN=0, VMAX=10 !TABLE, NAME=SIZE, columns=s,BSCx,BSCy,SIZEx,SIZEy !--- SET, ROLLRC, 1.0 !for survey -- turn RC corrector roll angle ON USE, (FODOLa,SPDASEL1B) PRINT, DASEL1B SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 ENDSUBROUTINE ENDCOMMENT ! ------------------------------------------------------------------------------ ! match DASEL Twiss -- stage-1A relaxed optics to ESA COMMENT MDASELS1A : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 0.0 !for twiss -- turn RC corrector roll angle OFF !!SETK2scDA SET, QDAS1b[K1], 0.0; SET, QDAS2b[K1], 0.0; SET, QDAS18b[K1], 0.0 SET, KQDAS1, -0.963105987165 SET, KQDAS2, 0.908079265384 SET, KQDAS11, 0.0 SET, KQDAS12, 0.0 SET, KQDAS13, 0.0 SET, KQDAS14, -0.331257876143 SET, KQDAS15, 0.0 SET, KQDAS16, 0.296918853966 SET, KQDAS17, 0.0 SET, KQDAS18, -0.576427105477 SET, KQDAS19, 0.0 SET, KQ10, 0.039778774082 SET, KQ11, -0.036654167038 SET, KQ19, 0.037245357449 SET, KQ20, -0.01736283793 SET, KQ27, -0.072119322567 SET, KQ28, 0.014825636448 SET, KQ30, -0.028939590569 SET, KQ38, 0.034948816176 USE, SPDASELA, RANGE=#S/ALWALL MATCH, BETA0=TWSSP VARY, KQDAS1 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !4 GeV VARY, KQDAS2 , STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS11, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS12, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS13, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS14, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS15, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS16, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS17, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQDAS18, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 !VARY, KQDAS19, STEP=1.E-4, LOWER=-1.627, UPPER= 1.627 VARY, KQ10, STEP=1.E-4 VARY, KQ11, STEP=1.E-4 VARY, KQ19, STEP=1.E-4 VARY, KQ20, STEP=1.E-4 VARY, KQ27, STEP=1.E-4 VARY, KQ28, STEP=1.E-4 VARY, KQ30, STEP=1.E-4 VARY, KQ38, STEP=1.E-4 CONSTR, ALWALL[1], BETX=7.834394963E+03, ALFX=4.252254496E+02, & BETY=4.018221983E+02, ALFY=2.208989467E+01, & DX=0, DPX=0, DY=0, DPY=0 WEIGHT, ALFY=1.E-4, ALFX=1.E-4 CONSTR, Q19[1], ALFX=0, ALFY=0 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, QDAS11[1] , BETX<220, BETY<220 CONSTR, QDAS12[1] , BETX<220, BETY<220 CONSTR, QDAS13[1] , BETX<220, BETY<220 CONSTR, QDAS14[1] , BETX<220, BETY<220 CONSTR, QDAS15[1] , BETX<220, BETY<220 CONSTR, QDAS16[1] , BETX<220, BETY<220 CONSTR, QDAS17[1] , BETX<220, BETY<220 CONSTR, QDAS18a[1], BETX<220, BETY<220 CONSTR, QDAS19[1] , BETX<220, BETY<220 CONSTR, BRDAS2b[1], BETX<220, BETY<220 RMATRIX, BEGDASEL/ENDDASEL, RM(5,6)=0, WEIGHT(5,6)=1.E-20 !LMDIF, TOL=1.E-20, CALLS=100000 ENDMATCH VALUE, KQDAS1,KQDAS2,KQDAS11,KQDAS12,KQDAS13,KQDAS14 VALUE, KQDAS15,KQDAS16,KQDAS17,KQDAS18,KQDAS19 VALUE, KQ10,KQ11,KQ19,KQ20,KQ27,KQ28,KQ30,KQ38 !--- !PRINT, FULL !PRINT, CLASS=QUADRUPOLE TWISS, BETA0=TWSSP, SAVE, COUPLE PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX,DY, & STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=DASEL, & VMIN=0,-4, VMAX=500,1 PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=rBETX,rBETY, VAXIS2=DX,DY, & STYLE=1, COLOUR=100, SPLINE, FILE="LCLS2sc", RANGE=DASELA, & VMIN=0,-10, VMAX=120,2 !--- TWISS, BETA0=TWSSP, COUPLE, SAVE=SIZE PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCx,SIZEx, STYLE=1, SPLINE=F, & COLOUR=100, FILE="LCLS2sc", RANGE=DASELA, VMIN=0, VMAX=160 PLOT, TABLE=SIZE, HAXIS=S, VAXIS=BSCy,SIZEy, STYLE=1, SPLINE=F, & COLOUR=100, FILE="LCLS2sc", RANGE=DASELA, VMIN=0, VMAX=30 !TABLE, NAME=SIZE, columns=s,BSCx,BSCy,SIZEx,SIZEy !--- !SET, ROLLRC, 1.0 !for survey -- turn RC corrector roll angle ON !USE, (FODOLa,SPDASELA) !PRINT, DASELA !SURVEY, X0=Xbyp, Y0=Ybyp, Z0=Zbyp, THETA0=0, PHI0=0, PSI0=0 ENDSUBROUTINE ENDCOMMENT ! ============================================================================== ! R56 checks ! ------------------------------------------------------------------------------ TWSSEXT : BETA0, ENERGY=Ef BEAM, ENERGY=Ei USE, (LCLS2scI,LCLS2scC) SAVEBETA, TWSSEXT, BEGEXT TWISS, BETA0=TWSSi MR56BYPDOG : SUBROUTINE BEAM, ENERGY=Ef USE, (EXT,DLBM) MATCH, BETA0=TWSSEXT WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/#E, TM(6,6,6)=0 RMATRIX, #S/#E, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56SPs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, SPRDs MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/#E, TM(6,6,6)=0 RMATRIX, #S/#E, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56SPh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, SPRDh MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/#E, TM(6,6,6)=0 RMATRIX, #S/#E, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56VBENDs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, LTUsc MATCH, BETA0=TWSSms WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, VBinb/VBoutb, TM(6,6,6)=0 RMATRIX, VBinb/VBoutb, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56VBENDh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (SPRDh,BSYH1,BSYH2,LTUH) MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, VBin/VBout, TM(6,6,6)=0 RMATRIX, VBin/VBout, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56LTUs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, LTUsc MATCH, BETA0=TWSSms WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/#E, TM(6,6,6)=0 RMATRIX, #S/#E, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56LTUh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (SPRDh,BSYH1,BSYH2,LTUH) MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, BEGBSYH_1/#E, TM(6,6,6)=0 RMATRIX, BEGBSYH_1/#E, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56DOGLTUs : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPs ; SET, SETCUS, CUSOFF USE, (EXT,DLBM,LCLS2scS) MATCH, BETA0=TWSSEXT WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/ENDLTUS, TM(6,6,6)=0 RMATRIX, #S/ENDLTUS, RM(6,6)=1 ENDMATCH TWISS, CHROM, BETA0=TWSSEXT, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & COLOUR=100, SPLINE=.T., RANGE=#S/ENDSXR, FILE="LCLS2sc", & TITLE="Dispersion from EXT to end of SXR undulator" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DDX,DDY, & COLOUR=100, SPLINE=.T., RANGE=#S/ENDSXR, FILE="LCLS2sc", & TITLE="2nd order dispersion from EXT to end of SXR undulator" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=WX,WY, & COLOUR=100, SPLINE=.T., RANGE=#S/ENDSXR, FILE="LCLS2sc", & TITLE="W-function from EXT to end of SXR undulator" ENDSUBROUTINE MR56DOGLTUh : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALOFF USE, (EXT,DLBM,LCLS2scH) MATCH, BETA0=TWSSEXT WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/ENDLTUH, TM(6,6,6)=0 RMATRIX, #S/ENDLTUH, RM(6,6)=1 ENDMATCH TWISS, CHROM, BETA0=TWSSEXT, SAVE=TWISS PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & COLOUR=100, SPLINE=.T., RANGE=#S/ENDHXR, FILE="LCLS2sc", & TITLE="Dispersion from EXT to end of HXR undulator" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DDX,DDY, & COLOUR=100, SPLINE=.T., RANGE=#S/ENDHXR, FILE="LCLS2sc", & TITLE="2nd order dispersion from EXT to end of HXR undulator" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=WX,WY, & COLOUR=100, SPLINE=.T., RANGE=#S/ENDHXR, FILE="LCLS2sc", & TITLE="W-function from EXT to end of HXR undulator" ENDSUBROUTINE MR56ALINE : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPh ; SET, SETAL, ALON USE, BSYLCLS2scA, RANGE=#S/ALWALL MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, BEGBSYA_1/ENDB, TM(6,6,6)=0 RMATRIX, BEGBSYA_1/ENDB, RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56DASEL : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 0.0 SETK2scDA USE, SPDASELA, RANGE=#S/ALWALL MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, BEGDASEL[1]/ENDDASEL[1], TM(6,6,6)=0 RMATRIX, BEGDASEL[1]/ENDDASEL[1], RM(6,6)=1 ENDMATCH ENDSUBROUTINE MR56DASELA : SUBROUTINE BEAM, ENERGY=Ef SET, SETSP, SETSPd ; SET, SETDA, DAON SET, ROLLRC, 0.0 SETK2scDA USE, SPDASELA, RANGE=#S/ALWALL MATCH, BETA0=TWSSP WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, BEGDASEL[1]/ALWALL[1], TM(6,6,6)=0 RMATRIX, BEGDASEL[1]/ALWALL[1], RM(6,6)=1 ENDMATCH ENDSUBROUTINE ! ============================================================================== ! COMMANDs ! ------------------------------------------------------------------------------ OPTION, ECHO ! SC linac systems !ML0E !ML1E !ML2E !ML3E !ML4E !MCMfodo1 !MCMfodo2 !MCMfodo3 !MC0fodo !MC1fodo !MHEATER !MBC1 !MBC2 !MCCDLU !MCCDLD !MASTRA !MTWSS0a !MTWSS0i !ML0LH !MLHD0 !MD0Coll0 !MColl0Coll1 !MColl1L2 !ML2BC2 !MBC2EMIT2 !MEMIT2L3 ML3L4 !QCM34 is dead !MDLGEO !MDLDXY !MDLBXY !MDLDDXY !ML4DL !MBYPGEO !MDLBYP !MDIAG0 !MBPMDG000 ! LEI systems !ML0LEIE !MTWSS0iLEI !MDLIX !MDLIDX !MDLITW !MDGIFODO !MDGITW !ML0LEIDLI !MDLILH !MDGITCAV ! ------------------------------------------------------------------------------ !COMMENT !MSTRT !MXYZSPh !MXYZSPs !MXYZSPd !MXYZBSY !MXYZDL2h !MXYZDL2s !MXYZDMP !MZCUSXR !MSPRDh !MSPRDh2 !added by MDW !MSPRDs !MSPRDs2 !added by MDW !MSPRDd !MBYPSP !MDL2Bb !MDL2Bb2 !added by MDW !MDL2Mb !MSPLTUs !MVBENDb !MDL2VB !MECELL !MED2P !MED2DL2b !MUNDs !MDUMPs2 !MDL2B !MDL2M2 !MBSY !MED2Ph !MED2DL2 !MUNDh !MDUMPh4 !MALINEsc !insert matched K-values in SETK2scA !ENDCOMMENT ! ------------------------------------------------------------------------------ ! DASEL !EXDAS := 1.49E-9 !rms emittance from M. Guetg simulation !EYDAS := 5.33E-9 !EPDAS := 0.51E-2 !rms dp/p EXDAS := 25.E-6/(Ef/0.511E-3) !3.19E-9 @ 4GeV, emittance from parameter table EYDAS := 25.E-6/(Ef/0.511E-3) !3.19E-9 @ 4GeV, emittance from parameter table EPDAS := 1.0E-2 !rms dp/p, assuming full width = +/-0.02 -- too large for A-line ABSC := 65.E-9 !BSC acceptance of LCLS2 (w/o 2 mm orbit) ! SIZE = DASEL full beam size (=3 sigma) per DASEL emittance plus dp/p=+/-0.02 ! BSC = DASEL BSC based on LCLS2 BSC (with 2 mm orbit) plus dp/p=+/-0.02 ! betatron size and dispersion are added linearly SIZEx : string, "1.E3*(3*sqrt(EXDAS*betx)+2*EPDAS*ABS(dx))" SIZEy : string, "1.E3*(3*sqrt(EYDAS*bety)+2*EPDAS*ABS(dy))" BSCx : string, "1.E3*(sqrt(ABSC*betx)+2*EPDAS*ABS(dx))+2.0" BSCy : string, "1.E3*(sqrt(ABSC*bety)+2*EPDAS*ABS(dy))+2.0" !COMMENT !MKICKDASEL !MXYZDASEL !MDASEL !insert matched results into SETK2scDA !ENDCOMMENT COMMENT !for DASEL stage-1 only !MDASELS1S !stage-1S with dump in S30 !MDASELS1B !stage-1B with dump in BSY !MDASELS1A !stage-1A with dump at SL10 ENDCOMMENT COMMENT VALUE, Ef,Brhof VALUE, QDAS1a[K1]*LQM*Brhof,QDAS1b[K1]*LQM*Brhof VALUE, QDAS2a[K1]*LQM*Brhof,QDAS2b[K1]*LQM*Brhof VALUE, QDAS11[K1]*LQM*Brhof,QDAS12[K1]*LQM*Brhof,QDAS13[K1]*LQM*Brhof VALUE, QDAS14[K1]*LQM*Brhof,QDAS15[K1]*LQM*Brhof,QDAS16[K1]*LQM*Brhof VALUE, QDAS17[K1]*LQM*Brhof,QDAS18a[K1]*LQM*Brhof,QDAS18b[K1]*LQM*Brhof VALUE, QDAS19[K1]*LQM*Brhof VALUE, ABKRDAS1*Brhof,ABKRDAS2*Brhof,ABKRDAS3*Brhof VALUE, ABKRDAS4*Brhof,ABKRDAS5*Brhof,ABKRDAS6*Brhof VALUE, ABLRDAS*Brhof,ABRDAS1*Brhof,ABRDAS2*Brhof,ABRAM1*Brhof ENDCOMMENT ! ------------------------------------------------------------------------------ COMMENT !compute R56 MR56BYPDOG MR56SPs MR56SPh MR56VBENDs MR56VBENDh MR56LTUs MR56LTUh MR56DOGLTUs MR56DOGLTUh MR56ALINE MR56DASEL MR56DASELA ENDCOMMENT RETURN ! ============================================================================== ! working stuff ! ------------------------------------------------------------------------------ ! match geometry of DC spreader (HXR and SXR are the same) ! (NOTE: enable DC spreader definitions SPRD.xsif and BSYsc.xsif) MXYZSPdch : SUBROUTINE BEAM, ENERGY=Ef USE, SPRDkhDC VALUE, DSPDCH1[L] MATCH, SURVEY, XS=0, YS=0, ZS=0, THETAS=0, PHIS=0, PSIS=0 VARY, DSPDCH1[L], STEP=1.E-6 WEIGHT, YS=1 CONSTR, BLXSPHa, YS=0.015 LMDIF, TOL=1.E-20 MIGRAD, TOL=1.E-20 ENDMATCH VALUE, DSPDCH1[L] PRINT, FULL SURVEY ENDSUBROUTINE ! ------------------------------------------------------------------------------ COMMENT DO, TIMES=20 SET, R56SPH, R56SPH+0.05E-3 VALUE, R56SPH MSPRDH ENDDO RETURN ENDCOMMENT COMMENT DO, TIMES=30 SET, R56SPS, R56SPS-0.1E-3 VALUE, R56SPS MSPRDS ; MSPLTUS ENDDO RETURN ENDCOMMENT ! ------------------------------------------------------------------------------ RETURN