! *** OPTICS=AD_ACCEL-14APR2025s *** ! LCLS2cu matching routines OPTION, -ECHO ! ------------------------------------------------------------------------------ TCUSXR0 := TCUSXR TBSY : BETA0 TWSa : BETA0 TWSb : BETA0 TWSm : BETA0 TWm : BETA0 TWm26 : BETA0 GLmaxQUM := 100 !BMAX = 100 kG @ 120 A (bipolar) for 1.26Q12 KmaxQUM15 := GLmaxQUM/LQA/(Cb*15) !max |K1| @ 15 GeV for QUM1-4 KmaxQUM10 := GLmaxQUM/LQA/(Cb*10) !max |K1| @ 10 GeV for QUM1-4 GLmaxQUMb := 67 !BMAX = 67 kG @ 120 A (bipolar) for 1.085Q4.31 KmaxQUMb15 := GLmaxQUMb/LQE/(Cb*15) !max |K1| @ 15 GeV for QUM1B-4B KmaxQUMb10 := GLmaxQUMb/LQE/(Cb*10) !max |K1| @ 10 GeV for QUM1B-4B GLmaxQu := 40 !BMAX = 40 kG @ 6 A (bipolar) for 0.433Q3.1 KmaxQu15 := GLmaxQu/LQu/(Cb*15) !max |K1| @ 15 GeV for HXR undulator quads KmaxQu10 := GLmaxQu/LQu/(Cb*10) !max |K1| @ 10 GeV for SXR undulator quads !VALUE, KmaxQUM15,KmaxQUMb15,KmaxQu15 !VALUE, KmaxQUM10,KmaxQUMb10,KmaxQu10 KminQUM1b := -1 ; KmaxQUM1b := +1 !bipolar KminQUM2b := -1 ; KmaxQUM2b := +1 !bipolar KminQUM3b := -1 ; KmaxQUM3b := +1 !bipolar KminQUM4b := -1 ; KmaxQUM4b := +1 !bipolar KminQSX16 := -1 ; KmaxQSX16 := +1 !bipolar KminQSX47 := -1 ; KmaxQSX47 := +1 !bipolar KminQUE1b := 0 ; KmaxQUE1b := +1 !unipolar QF KminQUE2b := -1 ; KmaxQUE2b := 0 !unipolar QD KminQDMPb := -1 ; KmaxQDMPb := 0 !unipolar QD MBETX := 1 MALFX := 0 MBETY := 1 MALFY := 0 ! ============================================================================== ! Original LCLS matching SUBROUTINEs ! ------------------------------------------------------------------------------ ML0E : SUBROUTINE !check final energy BEAM, ENERGY=E00 USE, (GUNL0a,L0aL0b) SAVEBETA, TWSa, L0BBEG SAVEBETA, TWSb, #E PRINT, FULL TWISS, COUPLE, BETA0=TWSSC, SAVE, TAPE="twiss.tape" VALUE, E0i,TWSa[ENERGY] VALUE, Ei,TWSb[ENERGY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML1E : SUBROUTINE !check BC1 energy BEAM, ENERGY=Ei USE, (L1,BC1) SAVEBETA, TWSm, XEND PRINT, FULL TWISS, BETX=10, BETY=10, SAVE, TAPE="twiss.tape" VALUE, EBC1,TWSm[ENERGY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML2E : SUBROUTINE !check BC2 energy BEAM, ENERGY=EBC1 USE, L2 SAVEBETA, TWSm, #E PRINT, FULL TWISS, BETX=10, BETY=10, SAVE, TAPE="twiss.tape" VALUE, EBC2,TWSm[ENERGY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML3E : SUBROUTINE !match final energy (LI25+LI26) BEAM, ENERGY=EBC2 USE, L3 VALUE, GradL3 MATCH, BETX=10, BETY=10 VARY, GradL3, STEP=1.E-6 WEIGHT, ENERGY=1 CONSTR, #E, ENERGY=Ef !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, GradL3 BEAM, ENERGY=EBC2 USE, L3 SAVEBETA, TWSm, #E PRINT, FULL TWISS, BETX=10, BETY=10, SAVE, TAPE="twiss.tape" VALUE, Ef,TWSm[ENERGY] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=ENERGY, & STYLE=100, SPLINE=.F., FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML1c : SUBROUTINE !match L1 phase advance per cell (coasting) VALUE, KQL1 USE, L1c CELL VARY, KQL1, STEP=1.E-5, LOWER=0, UPPER=100 CONSTR, #E, MUX=MU_L1, MUY=MU_L1 !LMDIF !MIGRAD ENDMATCH VALUE, 1.448192E-25 !expected penalty function with fitting OFF VALUE, KQL1 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML2c : SUBROUTINE !match L2 phase advance per cell (coasting) VALUE, KQL2 USE, L2c CELL VARY, KQL2, STEP=1.E-5, LOWER=0, UPPER=100 CONSTR, #E, MUX=MU_L2, MUY=MU_L2 !LMDIF !MIGRAD ENDMATCH VALUE, 2.277840E-24 !expected penalty function with fitting OFF VALUE, KQL2 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML3c : SUBROUTINE !match L3 phase advance per cell (coasting) VALUE, KQFL3,KQDL3 USE, L3c CELL VARY, KQFL3, STEP=1.E-5, LOWER=0 , UPPER=100 VARY, KQDL3, STEP=1.E-5, LOWER=-100, UPPER=0 CONSTR, #E, MUX=MUX_L3, MUY=MUY_L3 !LMDIF !MIGRAD ENDMATCH VALUE, 1.809822E-13 !expected penalty function with fitting OFF VALUE, KQFL3,KQDL3 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! (NOTE: optics between cathode and DL1beg is meaningless ... used only to ! initialize sumL and show beamline layout) MTWSSC : SUBROUTINE !find Twiss at cathode VALUE, CBETX,CALFX,CBETY,CALFY USE, GUNL0a MATCH, ENERGY=E00, BETX=CBETX, ALFX=CALFX, BETY=CBETY, ALFY=CALFY VARY, CBETX, STEP=1.E-4, LOWER=0.1, UPPER=100 VARY, CALFX, STEP=1.E-4, LOWER=-10, UPPER=10 VARY, CBETY, STEP=1.E-4, LOWER=0.1, UPPER=100 VARY, CALFY, STEP=1.E-4, LOWER=-10, UPPER=10 CONSTR, #E, BETX=TBETX, ALFX=TALFX, BETY=TBETY, ALFY=TALFY !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.434054E-15 !expected penalty function with fitting OFF VALUE, CBETX,CALFX,CBETY,CALFY ENDSUBROUTINE ! ------------------------------------------------------------------------------ MGSPEC : SUBROUTINE !match gun spectrometer VALUE, QG02[K1],QG03[K1] USE, (GUNBXG,GSPEC) MATCH, BETA0=TWSSC VARY, QG02[K1], STEP=1.E-4, LOWER=-100, UPPER=+100 VARY, QG03[K1], STEP=1.E-4, LOWER=-100, UPPER=+100 RMATRIX, YAG01/YAGG1, RM(1,2)=0, WEIGHT(1,2)=6, & RM(3,4)=0, WEIGHT(3,4)=6 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 3.013681E-13 !expected penalty function with fitting OFF VALUE, QG02[K1],QG03[K1] ENDSUBROUTINE ! ------------------------------------------------------------------------------ MHTRR56 : SUBROUTINE !set Heater-chicane R56 BEAM, ENERGY=Ei USE, LSRHTR SET, K_und, 0 !turn off undulator VALUE, ABh1,BBh1 MATCH, BETX=10, BETY=10 !VARY, BBh1, STEP=1.E-6, UPPER=0 VARY, ABh1, STEP=1.E-6, UPPER=0 RMATRIX, #S/#E, RM(5,6)=R56HTR, WEIGHT(5,6)=1 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.235206E-27 !expected penalty function with fitting OFF VALUE, ABh1,BBh1 PRINT, FULL SURVEY ENDSUBROUTINE ! ------------------------------------------------------------------------------ MBSYS : SUBROUTINE !match DL1 bend system dispersion VALUE, KQB USE, (L0aL0b,DL1) MATCH, BETA0=TWSS0 VARY, KQB, STEP=1.E-5, LOWER=-39, UPPER=+39 RMATRIX, BEGDL1_2/BX02B, RM(1,6)=0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 2.099758E-22 !expected penalty function with fitting OFF VALUE, KQB ENDSUBROUTINE ! ------------------------------------------------------------------------------ MED0 : SUBROUTINE !match from exit of L0a into ED0 BEAM, ENERGY=E0i VALUE, KQA01,KQA02,KQE01,KQE02,KQE03,KQE04 USE, (L0aL0b,DL1_1) MATCH, BETA0=TWSS0 VARY, KQA01, STEP=1.E-5, LOWER=-39, UPPER=+39 VARY, KQA02, STEP=1.E-5, LOWER=-39, UPPER=+39 VARY, KQE01, STEP=1.E-5, LOWER=-39, UPPER=+39 VARY, KQE02, STEP=1.E-5, LOWER=-39, UPPER=+39 VARY, KQE03, STEP=1.E-5, LOWER=-39, UPPER=+39 VARY, KQE04, STEP=1.E-5, LOWER=-39, UPPER=+39 !CONSTR, HTRUND, BETX<12,BETX>8,BETY<12,BETY>8 CONSTR, HTRUND, BETX=12,BETY=12 CONSTR, MRK0, BETX=BMIN0,ALFX=0.,BETY=BMIN0,ALFY=0. !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 7.202187E-23 !expected penalty function with fitting OFF (htr fully ON) VALUE, KQA01,KQA02,KQE01,KQE02,KQE03,KQE04 BEAM, ENERGY=E0i USE, (L0aL0b,DL1_1) PRINT, FULL TWISS, COUPLE, BETA0=TWSS0, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MED2L1 : SUBROUTINE !match from ED0 through injection bend system to L1 BEAM, ENERGY=E0i USE, (L0aL0b,DL1,L1) VALUE, KQM01,KQM02,KQM03,KQM04 MATCH, BETA0=TWSS0 VARY, KQM01, STEP=1.E-5, LOWER=-39,UPPER=+39 VARY, KQM02, STEP=1.E-5, LOWER=-39,UPPER=+39 VARY, KQM03, STEP=1.E-5, LOWER=-39,UPPER=+39 VARY, KQM04, STEP=1.E-5, LOWER=-39,UPPER=+39 CONSTR, QA11[1], BETY=BET11, BETX=BET12 CONSTR, QA12[1], BETY=BET12, BETX=BET11 CONSTR, QM03[1], BETY<10 CONSTR, OTR4, BETX<1.7, BETY>8.0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 8.576116E-02 !expected penalty function with fitting OFF VALUE, KQM01,KQM02,KQM03,KQM04 BEAM, ENERGY=E0i USE, (L0aL0b,DL1_1,L1) PRINT, FULL TWISS, COUPLE, BETA0=TWSS0, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MSPEC : SUBROUTINE !match 135-MeV spectrometer BEAM, ENERGY=E0i USE, (L0aL0b,DL1_1,SPECBL) VALUE, KQS01,KQS02,KQM01,KQM02 MATCH, BETA0=TWSS0 VARY, KQS01, STEP=1.E-5, LOWER=-40, UPPER=+40 ! nominal fit at YAGS2 VARY, KQS02, STEP=1.E-5, LOWER=-40, UPPER=+40 ! nominal fit at YAGS2 !!VARY, KQM01, STEP=1.E-5, LOWER=-40, UPPER=+40 ! only use when 3-keV resolution on OTR required !!VARY, KQM02, STEP=1.E-5, LOWER=-40, UPPER=+40 ! only use when 3-keV resolution on OTR required WEIGHT, BETX=10 ! all fits CONSTR, YAGS2, BETX<0.120,BETY>5,BETY<12 ! nominal fit at YAGS2 !!CONSTR, OTRS1, BETX<0.07,BETY>5,BETY<12 ! only use when 3-keV resolution on OTR required CONSTR, YAGS2, DX<-0.9*E0i/Ei ! nominal fit at YAGS2 !!CONSTR, OTRS1, DX=-1.0*E0i/Ei ! only use when 3-keV resolution on OTR required !!CONSTR, DS2, BETX<35 ! only use when 3-keV resolution on OTR required !!CONSTR, BPM12, BETY<35 ! only use when 3-keV resolution on OTR required !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.377848E-02 !expected penalty function with fitting OFF VALUE, KQS01,KQS02,KQM01,KQM02 BEAM, ENERGY=E0i USE, (L0aL0b,DL1_1,SPECBL) PRINT, FULL TWISS, COUPLE, BETA0=TWSS0, 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 ! ------------------------------------------------------------------------------ MBC1R56 : SUBROUTINE !match BC1 R56 BEAM, ENERGY=EBC1 VALUE, AB11,AB11S,LB11S,AB11L,LB11L,BB11 USE, BC1c MATCH, BETX=10, BETY=10 VARY, AB11, STEP=1.E-6, UPPER=0 RMATRIX, #S/#E, RM(5,6)=R56BC1 !0.0455 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 2.068490E-26 !expected penalty function with fitting OFF VALUE, AB11,AB11S,LB11S,AB11L,LB11L,BB11 PRINT, FULL SURVEY ENDSUBROUTINE ! ------------------------------------------------------------------------------ MBC1e : SUBROUTINE !match through BC1 chicane into emittance measurement system BEAM, ENERGY=E0i USE, (L0aL0b,DL1,L1,BC1) VALUE, KQA12,KQ21201,KQM11,KQM12,KQM13 MATCH, BETA0=TWSS0 VARY, KQA12, STEP=1.E-6, LOWER=-100, UPPER=+100 VARY, KQ21201, STEP=1.E-6, LOWER=-100, UPPER= 0 VARY, KQM11, STEP=1.E-6, LOWER=-100, UPPER=+100 VARY, KQM12, STEP=1.E-6, LOWER=-100, UPPER=+100 VARY, KQM13, STEP=1.E-6, LOWER=-100, UPPER=+100 CONSTR, BX14B, BETX<4, ALFX<0 CONSTR, WS12, BETX=BMIN1, ALFX=0, BETY=BMIN1, ALFY=0 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.675593E-22 !expected penalty function with fitting OFF VALUE, KQA12,KQ21201,KQM11,KQM12,KQM13 BEAM, ENERGY=E0i USE, (L0aL0b,DL1,L1,BC1) PRINT, FULL TWISS, COUPLE, BETA0=TWSS0, 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 ! ------------------------------------------------------------------------------ MBC1L2 : SUBROUTINE !match BC1 to L2 BEAM, ENERGY=E0i VALUE, KQM14,KQM15,KQ21401,KQ21501 USE, (L0aL0b,DL1,L1,BC1,LI21) MATCH, BETA0=TWSS0 !!VARY, KQ21301, STEP=1.E-5, LOWER=-10, UPPER= 0 !leave this quad OFF - unhook cables? VARY, KQM14, STEP=1.E-5, LOWER= 0, UPPER=10 VARY, KQM15, STEP=1.E-5, LOWER=-10, UPPER= 0 VARY, KQ21401, STEP=1.E-5, LOWER= 0, UPPER=10 VARY, KQ21501, STEP=1.E-5, LOWER=-10, UPPER= 0 CONSTR, QM15[1], BETY10 CONSTR, Q24701A[2], BETY<50 CONSTR, Q24701A[2], BETX<50 CONSTR, Q24901A[2], BETY<75 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 7.347058E-13 !expected penalty function with fitting OFF VALUE, KQ24501,KQ24601,KQ24701,KQM21 ENDSUBROUTINE ! ------------------------------------------------------------------------------ MBC2L3 : SUBROUTINE !match BC2 to L3 BEAM, ENERGY=E0i VALUE, KQM22,KQ24901,KQ25201,KQ25301,KQ25401,KQ25501 USE, (L0aL0b,DL1,L1,BC1,L2,BC2,LI25) MATCH, BETA0=TWSS0 VARY, KQM22, STEP=1.E-6, LOWER= -0.65, UPPER= 0 VARY, KQ24901, STEP=1.E-6 !, LOWER=-10, UPPER= 0 VARY, KQ25201, STEP=1.E-6 !, LOWER= 0, UPPER=10 VARY, KQ25301, STEP=1.E-6 !, LOWER=-10, UPPER= 0 VARY, KQ25401, STEP=1.E-6 !, LOWER= 0, UPPER=10 VARY, KQ25501, STEP=1.E-6 !, LOWER=-10, UPPER= 0 CONSTR, Q25201[1], BETX<60 CONSTR, Q25301[1], BETY<64 CONSTR, Q25401[1], BETX=BET31, BETY=BET32 CONSTR, Q25501[1], BETX=BET33, BETY=BET34 CONSTR, Q25601[1], BETX=BET31, BETY=BET32 CONSTR, Q25701[1], BETX=BET33, BETY=BET34 CONSTR, Q25801[1], BETX=BET31, BETY=BET32 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 4.578278E-01 !expected penalty function with fitting OFF VALUE, KQM22,KQ24901,KQ25201,KQ25301,KQ25401,KQ25501 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML26 : SUBROUTINE !match LI26 betas VALUE, KQ25801,KQ25901,KQ26201,KQ26301 USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3) MATCH, BETA0=TWSS0 VARY, KQ25801, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ25901, STEP=1.E-5 !, LOWER=-100, UPPER= 0 VARY, KQ26201, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ26301, STEP=1.E-5 !, LOWER=-100, UPPER= 0 CONSTR, Q26301[1], BETX=BET33, BETY=BET34 CONSTR, Q26401[1], BETX=BET31, BETY=BET32 CONSTR, Q26501[1], BETX=BET33, BETY=BET34 CONSTR, Q26601[1], BETX=BET31, BETY=BET32 CONSTR, Q26701[1], BETX=BET33, BETY=BET34 CONSTR, Q26801[1], BETX=BET31, BETY=BET32 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.522531E-02 !expected penalty function with fitting OFF VALUE, KQ25801,KQ25901,KQ26201,KQ26301 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML27 : SUBROUTINE !match LI27 betas VALUE, KQ26801,KQ26901,KQ27201,KQ27301 USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3) MATCH, BETA0=TWSS0 VARY, KQ26801, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ26901, STEP=1.E-5 !, LOWER=-100, UPPER= 0 VARY, KQ27201, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ27301, STEP=1.E-5 !, LOWER=-100, UPPER= 0 CONSTR, Q27301[1], BETX=BET33, BETY=BET34 CONSTR, Q27401[1], BETX=BET31, BETY=BET32 CONSTR, Q27501[1], BETX=BET33, BETY=BET34 CONSTR, Q27601[1], BETX=BET31, BETY=BET32 CONSTR, Q27701[1], BETX=BET33, BETY=BET34 CONSTR, Q27801[1], BETX=BET31, BETY=BET32 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 5.351428E-03 !expected penalty function with fitting OFF VALUE, KQ26801,KQ26901,KQ27201,KQ27301 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML28 : SUBROUTINE !match LI28 betas VALUE, KQ27801,KQ27901,KQ28201,KQ28301 USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3) MATCH, BETA0=TWSS0 VARY, KQ27801, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ27901, STEP=1.E-5 !, LOWER=-100, UPPER= 0 VARY, KQ28201, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ28301, STEP=1.E-5 !, LOWER=-100, UPPER= 0 CONSTR, Q28301[1], BETX=BET33, BETY=BET34 CONSTR, Q28401[1], BETX=BET31, BETY=BET32 CONSTR, Q28501[1], BETX=BET33, BETY=BET34 CONSTR, Q28601[1], BETX=BET31, BETY=BET32 CONSTR, Q28701[1], BETX=BET33, BETY=BET34 CONSTR, Q28801[1], BETX=BET31, BETY=BET32 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 2.676795E-03 !expected penalty function with fitting OFF VALUE, KQ27801,KQ27901,KQ28201,KQ28301 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML29 : SUBROUTINE !match LI29 betas VALUE, KQ28801,KQ28901,KQ29201,KQ29301 USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3) MATCH, BETA0=TWSS0 VARY, KQ28801, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ28901, STEP=1.E-5 !, LOWER=-100, UPPER= 0 VARY, KQ29201, STEP=1.E-5 !, LOWER= 0, UPPER=100 VARY, KQ29301, STEP=1.E-5 !, LOWER=-100, UPPER= 0 CONSTR, Q29301[1], BETX=BET33, BETY=BET34 CONSTR, Q29401[1], BETX=BET31, BETY=BET32 CONSTR, Q29501[1], BETX=BET33, BETY=BET34 CONSTR, Q29601[1], BETX=BET31, BETY=BET32 CONSTR, Q29701[1], BETX=BET33, BETY=BET34 CONSTR, Q29801[1], BETX=BET31, BETY=BET32 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.515926E-03 !expected penalty function with fitting OFF VALUE, KQ28801,KQ28901,KQ29201,KQ29301 ENDSUBROUTINE ! ------------------------------------------------------------------------------ ML30 : SUBROUTINE !match "periodic" (FODO cell) LI30 betas BEAM, ENERGY=E0i SET, KQ30401, KQFL3 SET, KQ30501, KQDL3 SET, KQ30601, KQFL3 SET, KQ30701, KQDL3 SET, KQ30801, KQFL3 VALUE, KQ29801,KQ29901,KQ30201,KQ30301 USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3) MATCH, BETA0=TWSS0 VARY, KQ29801, STEP=1.E-6, LOWER=0 VARY, KQ29901, STEP=1.E-6, UPPER=0 VARY, KQ30201, STEP=1.E-6, LOWER=0 VARY, KQ30301, STEP=1.E-6, UPPER=0 CONSTR, Q30301[1], BETX=BET33, BETY=BET34 CONSTR, Q30401[1], BETX=BET31, BETY=BET32 CONSTR, Q30501[1], BETX=BET33, BETY=BET34 CONSTR, Q30601[1], BETX=BET31, BETY=BET32 CONSTR, Q30701[1], BETX=BET33, BETY=BET34 CONSTR, Q30801[1], BETX=BET31, BETY=BET32 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 1.694700E-01 !expected penalty function with fitting OFF VALUE, KQ29801,KQ29901,KQ30201,KQ30301 ENDSUBROUTINE ! ============================================================================== ! LCLS-II SXR matching SUBROUTINEs ! ------------------------------------------------------------------------------ MKICKCUSXR : SUBROUTINE !match geometry and Twiss of CUSXR line (1) BEAM, ENERGY=Ef SET, TCUSXR0, TCUSXR SET, TCUSXR, 0.0 SET, SETCUS, CUSON USE, CUSXR, RANGE=#S/BLRCUSb[1] VALUE, ABKRCUS0,LDCUSBLRa*COS(ABKRCUS) MATCH, BETA0=TWSCuS1, SURVEY, XS=0, YS=0, ZS=0, THETAS=0, PHIS=0, PSIS=0 VARY, ABKRCUS0, STEP=1.E-6 VARY, LDCUSBLRa, STEP=1.E-4 WEIGHT, YS=1 CONSTR, BLRCUSa[1], YS=15.E-3 CONSTR, BPMCUS[1], YS=5.E-3 !LMDIF, TOL=1.E-20, CALLS=1000 ENDMATCH VALUE, ABKRCUS0,LDCUSBLRa*COS(ABKRCUS) PRINT, RANGE=#S/BLRCUSb[1] SURVEY SET, TCUSXR, TCUSXR0 ENDSUBROUTINE MXYZCUSXR : SUBROUTINE !match geometry and Twiss of CUSXR line (2) BEAM, ENERGY=Ef SET, SETCUS, CUSON USE, CUSXR VALUE, TCUSXR,TBRCUS,ABLRCUS0,dABRCUS,dLDCUS0,AROCUS1,AROCUS2,AROCUS3 MATCH, BETA0=TWSCuS1, SURVEY, XS=XCuS1, YS=YCuS1, ZS=ZCuS1, & THETAS=THCuS1, PHIS=PHCuS1, PSIS=PSCuS1 VARY, TCUSXR, STEP=1.E-5 VARY, TBRCUS, STEP=1.E-5 VARY, ABLRCUS0, STEP=1.E-5 VARY, dABRCUS, STEP=1.E-5 VARY, dLDCUS0, STEP=1.E-4 VARY, AROCUS1, STEP=1.E-6 VARY, AROCUS2, STEP=1.E-6 VARY, AROCUS3, STEP=1.E-6 WEIGHT, XS=1, YS=1, ZS=1, THETAS=1, PHIS=1, PSIS=1 CONSTR, #E, XS=XCuS2, YS=YCuS2, ZS=ZCuS2, & THETAS=THCuS2, PHIS=PHCuS2, PSIS=PSCuS2 CONSTR, QCUS1[1] , PSIS=0 CONSTR, MIDCUSXR , PSIS=0 CONSTR, QCuS10[1], PSIS=0 !LMDIF, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, TCUSXR,TBRCUS,ABLRCUS0,dABRCUS,dLDCUS0,AROCUS1,AROCUS2,AROCUS3 PRINT, FULL SURVEY, X0=XCuS1, Y0=YCuS1, Z0=ZCuS1, & THETA0=THCuS1, PHI0=PHCuS1, PSI0=PSCuS1 ENDSUBROUTINE MCUSXR : SUBROUTINE !match geometry and Twiss of CUSXR line (3) COMMENT !use this COMMENT/UNCOMMENT block to set-up BEAM, ENERGY=E00 !TBXCUS1 (copy result to CUSXR.xsif), then SET, SETCUS, CUSON !turn it off, otherwise it leads to a larger USE, CU_SXR !penalty value in MCUSXR, when verifying the SAVEBETA, TWSCUS1, BEGCLTS !match (LMDIF off) with prior MXYZCUSXR routine TWISS, BETA0=TWSSC SET, TBXCUS1, TWSCUS1[BETX]; SET, TBYCUS1, TWSCUS1[BETY] SET, TAXCUS1, TWSCUS1[ALFX]; SET, TAYCUS1, TWSCUS1[ALFY] ENDCOMMENT BEAM, ENERGY=Ef SET, SETCUS, CUSON USE, CUSXR VALUE, KQCuS1,KQCuS2,KQCuS3,KQCuS4,KQCuS5,KQCuS9,KQCuS10 MATCH, BETA0=TWSCuS1, ENERGY=Ef VARY, KQCuS1 , STEP=1.E-4, LOWER=-0.570, UPPER=0.570 !10 GeV VARY, KQCuS2 , STEP=1.E-4, LOWER=-0.570, UPPER=0.570 VARY, KQCuS3 , STEP=1.E-4, LOWER=-0.570, UPPER=0.570 VARY, KQCuS4 , STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !10 GeV VARY, KQCuS5 , STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !!VARY, KQCuS6 , STEP=1.E-4, LOWER=-2.807, UPPER=2.807 !!VARY, KQCuS7 , STEP=1.E-4, LOWER=-0.570, UPPER=0.570 !!VARY, KQCuS8 , STEP=1.E-4, LOWER=-0.570, UPPER=0.570 VARY, KQCuS9 , STEP=1.E-4, LOWER=-0.570, UPPER=0.570 !10 GeV VARY, KQCuS10, STEP=1.E-4, LOWER=-0.570, UPPER=0.570 !!VARY, dLDCuS2 , STEP=1.E-4, LOWER=-2.0, UPPER=2.0 !!VARY, dLDCuS3 , STEP=1.E-4, LOWER=-2.0, UPPER=2.0 !!VARY, dLDCuS4 , STEP=1.E-4, LOWER=-2.0, UPPER=1.0 !!VARY, ABYCUS , STEP=1.E-5 RMATRIX, #S/#E, RM(5,6)=0, WEIGHT(5,6)=100 WEIGHT, DX=10, DY=10, DPX=100, DPY=100 CONSTR, #E, DX=0, DY=0, DPX=0, DPY=0 !!CONSTR, MIDCUSXR, DPX<0, DPY>0 CONSTR, MIDCUSXR, ALFX=0, ALFY=0 WEIGHT, BETX=1.E-6, BETY=1.E-6 CONSTR, QCuS2[1], BETX<30, BETY<155 CONSTR, QCuS3[1], BETX<30, BETY<130 CONSTR, QCuS4[1], BETX<30, BETY<30 CONSTR, QCuS5[1], BETX<30, BETY<30 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-6, BETY=1.E-6, ALFX=1.E-5, ALFY=1.E-5 WEIGHT, ALFX=0, ALFY=0 COUPLE, BLRCUSa[1]/BRCUS1a[1] COUPLE, QCUS1[1]/QCUS10[1] COUPLE, QCUS2[1]/QCUS9[1] COUPLE, QCUS3[1]/QCUS8[1] COUPLE, QCUS4[1]/QCUS7[1] !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, TBXCUS1,TBYCUS1,TAXCUS1,TAYCUS1 VALUE, KQCuS1,KQCuS2,KQCuS3,KQCuS4,KQCuS5,KQCuS9,KQCuS10 !!VALUE, KQCuS8!,KQCuS7,KQCuS6 !!VALUE, dLDCuS2,dLDCuS3,dLDCuS4 !!VALUE, ABYCUS PRINT, MIDCUSXR TWISS, BETA0=TWSCuS1, SAVE=TWISS !, TAPE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, TITLE="Cu-linac to SXR", FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE, TITLE="Cu-linac to SXR", FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCLTS : SUBROUTINE !match into CLTS (input Twiss defined in CUSXR.xsif) BEAM, ENERGY=E0i USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3,CLTH0,CLTH1) VALUE, KQ30701,KQ30801,KQ50Q1,KQ50Q2 MATCH, BETA0=TWSS0 VARY, KQ30701, STEP=1.E-6, UPPER=0 VARY, KQ30801, STEP=1.E-6, LOWER=0 VARY, KQ50Q1, STEP=1.E-6, UPPER=0 VARY, KQ50Q2, STEP=1.E-6, LOWER=0 CONSTR, #E, BETX=TBXCUS1, ALFX=TAXCUS1, BETY=TBYCUS1, ALFY=TAYCUS1 !LMDIF, TOL=1.E-20 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, 2.435538E-19 !expected penalty function with fitting OFF VALUE, KQ30701,KQ30801,KQ50Q1,KQ50Q2 BEAM, ENERGY=E0i USE, (L0aL0b,DL1,L1,BC1,L2,BC2,L3,CLTH0,CLTH1,CLTS) PRINT, FULL TWISS, BETA0=TWSS0, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, RANGE=LI30BEG/#E, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ MCULTUs : SUBROUTINE ! match from CUSXR to horizontal dogleg bends in SXR LTU BEAM, ENERGY=E00 SET, SETCUS, CUSON SETK2cuS USE, CU_SXR VALUE, KQBP33,KQBP34,KQDBL1,KQDBL2,KQDL11,Ef MATCH, ENERGY=E00, BETA0=TWSSC 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 VARY, KQDBL1, STEP=1.E-4, LOWER=-1.871, UPPER=1.871 !15 GeV VARY, KQDBL2, STEP=1.E-4, LOWER=-1.871, UPPER=1.871 VARY, KQDL11, STEP=1.E-4, LOWER= 0 , UPPER=0.570 !10 GeV CONSTR, QBP34[1], BETX<100 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 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 WEIGHT, BETX=1.E-6, BETY=1.E-6 COUPLE, QDL11[1]/QDL13[1] !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQBP33,KQBP34,KQDBL1,KQDBL2,KQDL11 PRINT, FULL TWISS, SAVE, BETA0=TWSSC PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE=.T., FILE="match", RANGE=BEGCLTS/QDL20[2] PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, & STYLE=100, SPLINE=.T., FILE="match", RANGE=BEGCLTS/QDL20[2] ENDSUBROUTINE ! ------------------------------------------------------------------------------ MSXR : SUBROUTINE !match into SXR undulator BEAM, ENERGY=E00 SET, SETCUS, CUSON SETK2cuS USE, CU_SXR SAVEBETA, TWm, ENDL3 TWISS, BETA0=TWSSC BEAM, ENERGY=Ef USE, (SXCEL26,SXCEL27,SXCEL28,SXCEL29, & SXCEL30,SXCEL31,SXCEL32,SXCEL33) SAVEBETA, TWm26, #S TWISS SET, MBETX, TWm26[BETX] SET, MALFX, TWm26[ALFX] SET, MBETY, TWm26[BETY] SET, MALFY, TWm26[ALFY] SET, KminQUM1b, -4.272301476446 ; SET, KmaxQUM1b, 4.262158217579 SET, KminQUM2b, -4.277101524229 ; SET, KmaxQUM2b, 4.263591916062 SET, KminQUM3b, -4.236370592453 ; SET, KmaxQUM3b, 4.285204762788 SET, KminQUM4b, -4.258539234185 ; SET, KmaxQUM4b, 4.243501189751 SET, KminQSX16, -3.21206205 ; SET, KmaxQSX16, 3.21206205 BEAM, ENERGY=Ef USE, BSYLTUS VALUE, KQUM1b,KQUM2b,KQUM3b,KQUM4b,kQSX16 MATCH, BETX=TWm[BETX], ALFX=TWm[ALFX], BETY=TWm[BETY], ALFY=TWm[ALFY] VARY, kQUM1b, STEP=1.E-6, LOWER=KminQUM1b, UPPER=KmaxQUM1b VARY, kQUM2b, STEP=1.E-6, LOWER=KminQUM2b, UPPER=KmaxQUM2b VARY, kQUM3b, STEP=1.E-6, LOWER=KminQUM3b, UPPER=KmaxQUM3b VARY, kQUM4b, STEP=1.E-6, LOWER=KminQUM4b, UPPER=KmaxQUM4b VARY, kQSX16, STEP=1.E-6, LOWER=KminQSX16, UPPER=KmaxQSX16 CONSTR, SXR26BEG, BETX=MBETX, ALFX=MALFX, BETY=MBETY, ALFY=MALFY WEIGHT, BETX=1.E-03,BETY=1.E-03 CONSTR, QUM1B[1], BETX<76, BETY<76 CONSTR, QUM2B[1], BETX<76, BETY<76 CONSTR, QUM3B[1], BETX<76, BETY<76 CONSTR, QUM4B[1], BETX<76, BETY<76 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQUM1b,KQUM2b,KQUM3b,KQUM4b,kQSX16 PRINT, FULL TWISS, BETX=TWm[BETX], ALFX=TWm[ALFX], BETY=TWm[BETY], ALFY=TWm[ALFY], & COUPLE, SAVE, TAPE="twiss.tape", RTAPE="rmat.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, VMIN=0, VMAX=80, & STYLE=100, SPLINE, RANGE=WS33B/SXR46END, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=SXR26BEG/SXR33END, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ comment !use MSXR instead MUNDs : SUBROUTINE !match SXR undulator (no QSX16, normal SXR/DMPS polarity) BEAM, ENERGY=E00 SET, SETCUS, CUSON SETK2cuS USE, CU_SXR SAVEBETA, TWm, ENDL3 TWISS, BETA0=TWSSC SHOW, TWm VALUE, QSXh17[K1],QSXh18[K1],QSXh19[K1],QSXh20[K1],QSXh21[K1],& QSXh22[K1],QSXh23[K1],QSXh24[K1],QSXh25[K1] VALUE, QSXh26[K1],QSXh27[K1],QSXh28[K1],QSXh29[K1],QSXh30[K1],& QSXh31[K1],QSXh32[K1],QSXh33[K1],QSXh34[K1],QSXh35[K1] VALUE, QSXh36[K1],QSXh37[K1],QSXh38[K1],QSXh39[K1],QSXh40[K1],& QSXh41[K1],QSXh42[K1],QSXh43[K1],QSXh44[K1],QSXh45[K1] VALUE, QSXh46[K1],QSXh47[K1] BEAM, ENERGY=TWm[ENERGY] USE, BSYLTUS VALUE, KQUM1b,KQUM2b,KQUM3b,KQUM4b MATCH, BETX=TWm[BETX], ALFX=TWm[ALFX], BETY=TWm[BETY], ALFY=TWm[ALFY] VARY, KQUM1b, STEP=1.E-6, LOWER=0, UPPER=+KmaxQUMb10 VARY, KQUM2b, STEP=1.E-6, UPPER=0, LOWER=-KmaxQUMb10 VARY, KQUM3b, STEP=1.E-6, UPPER=0, LOWER=-KmaxQUMb10 VARY, KQUM4b, STEP=1.E-6, LOWER=0, UPPER=+KmaxQUMb10 WEIGHT, X=0, PX=0, BETX=1, ALFX=0, DX=0, DPX=0 WEIGHT, Y=0, PY=0, BETY=1, ALFY=0, DY=0, DPY=0 WEIGHT, T=0, PT=0 COUPLE, QSXh26[1]/QSXh28[1] COUPLE, QSXh27[1]/QSXh29[1] COUPLE, QSXh26[1]/QSXh30[1] COUPLE, QSXh27[1]/QSXh31[1] COUPLE, QSXh26[1]/QSXh32[1] COUPLE, QSXh27[1]/QSXh33[1] COUPLE, QSXh26[1]/QSXh34[1] COUPLE, QSXh27[1]/QSXh35[1] COUPLE, QSXh26[1]/QSXh36[1] COUPLE, QSXh27[1]/QSXh37[1] COUPLE, QSXh26[1]/QSXh38[1] COUPLE, QSXh27[1]/QSXh39[1] COUPLE, QSXh26[1]/QSXh40[1] COUPLE, QSXh27[1]/QSXh41[1] COUPLE, QSXh26[1]/QSXh42[1] COUPLE, QSXh27[1]/QSXh43[1] COUPLE, QSXh26[1]/QSXh44[1] COUPLE, QSXh27[1]/QSXh45[1] COUPLE, QSXh26[1]/QSXh46[1] !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQUM1b,KQUM2b,KQUM3b,KQUM4b PRINT, FULL TWISS, BETX=TWm[BETX], ALFX=TWm[ALFX], BETY=TWm[BETY], ALFY=TWm[ALFY], & SAVE, TAPE="twiss.tape", RTAPE="rmat.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=WS33B/SXR47END, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=BEGBSYS/SXR47END, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=DX,DY, RANGE=BEGBSYS/SXR47END, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE endcomment ! ------------------------------------------------------------------------------ ! matching TCAV per Yuantao's requirements MDUMPs : SUBROUTINE !match dumpline in SXR BEAM, ENERGY=Ef SET, SETCUS, CUSON SETK2cuS USE, (CLTS,BSYS,LTUS,SXRUND,DUMPLINEb) VALUE, kQSX47,KQUE1b,KQUE2b,KQDMPb MATCH, BETA0=TWSCuS1 !VARY, kQSX47, STEP=1.E-4, LOWER=-0.994, UPPER= 0 ! 0.994 !15 GeV VARY, KQUE1b, STEP=1.E-4, UPPER= 0.380, LOWER=-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 ! 0.182 !17 GeV !WEIGHT, BETX=1.E-5 !CONSTR, MTCXB, BETX>85 !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.8, BETX>20 !need large betx WEIGHT, X=0, PX=0, BETX=0, ALFX=0, DX=0, DPX=0, & Y=0, PY=0, BETY=0, ALFY=0, DY=0, DPY=0, & T=0, PT=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, kQSX47,KQUE1b,KQUE2b,KQDMPb PRINT, MTCXB ; PRINT, OTRDMPB TWISS, BETA0=TWSCuS1, 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=TWSCuS1, 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 ! ============================================================================== ! LCLS-II HXR matching SUBROUTINEs ! ------------------------------------------------------------------------------ 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 MDL2M2 : SUBROUTINE !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 ! ------------------------------------------------------------------------------ MBSY : SUBROUTINE !match BSY to HXR LTU BEAM, ENERGY=ENBEGL3 SET, SETCUS, CUSOFF SETK2cuH USE, (L3,CLTH0,CLTH1,CLTH2,BSYH1,BSYH2,LTUH) VALUE, KQ50Q3,KQ4,KQ5,KQ6,KQA0 MATCH, BETA0=TWBEGL3 !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 CONSTR, MM1, BETX=MBETXh, BETY=MBETYh, ALFX=MALFXh, ALFY=MALFYh comment WEIGHT, ALFX=2.E-5, ALFY=2.E-5 CONSTR, QVB2[1], ALFX=0, ALFY=0 WEIGHT, BETX=1.E-5, BETY=1.E-5 CONSTR, Q50Q3[1], BETX<140, BETY<105 CONSTR, Q4[1] , BETX<105, BETY<140 CONSTR, Q5[1] , BETX<105, BETY<105 CONSTR, Q6[1] , BETX<105, BETY<115 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 endcomment !LMDIF, TOL=1.E-20 !, CALLS=20000 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQ50Q3,KQ4,KQ5,KQ6,KQA0 TWISS, BETA0=TWBEGL3, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=Q30401[1]/MUWALL, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=Q30401[1]/MM2, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! NOTE: keep SC values for QEM1-4 ... don't rematch MED2DL2 : SUBROUTINE !match from bend cells to emit-diag section BEAM, ENERGY=Ef USE, LTUM VALUE, KQEM1,KQEM2,KQEM3,KQEM4 MATCH, BETA0=TWSSmh 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=TWSSmh, SAVE PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------ ! NOTE: keep Yuri's energy-scalable values for QUM1-4 ... don't rematch MUNDh : SUBROUTINE !match through the HXR undulator BEAM, ENERGY=Ef SET, SETCUS, CUSOFF SETK2cuH USE, (LTUH,HXRUND), RANGE=MM1/#E VALUE, KQUM1,KQUM2,KQUM3,KQUM4 MATCH, BETA0=TWSSmh VARY, KQUM1, STEP=1.E-6, LOWER=-KmaxQUM15, UPPER=+KmaxQUM15 !15 GeV VARY, KQUM2, STEP=1.E-6, LOWER=-KmaxQUM15, UPPER=+KmaxQUM15 VARY, KQUM3, STEP=1.E-6, LOWER=-KmaxQUM15, UPPER=+KmaxQUM15 VARY, KQUM4, STEP=1.E-6, LOWER=-KmaxQUM15, UPPER=+KmaxQUM15 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, QHXh13[1]/QHXh15[1] COUPLE, QHXh14[1]/QHXh16[1] COUPLE, QHXh13[1]/QHXh17[1] COUPLE, QHXh14[1]/QHXh18[1] COUPLE, QHXh13[1]/QHXh19[1] COUPLE, QHXh14[1]/QHXh20[1] COUPLE, QHXh13[1]/QHXh21[1] COUPLE, QHXh14[1]/QHXh22[1] COUPLE, QHXh13[1]/QHXh23[1] COUPLE, QHXh14[1]/QHXh24[1] COUPLE, QHXh13[1]/QHXh25[1] COUPLE, QHXh14[1]/QHXh26[1] COUPLE, QHXh13[1]/QHXh27[1] COUPLE, QHXh14[1]/QHXh28[1] COUPLE, QHXh13[1]/QHXh29[1] COUPLE, QHXh14[1]/QHXh30[1] COUPLE, QHXh13[1]/QHXh31[1] COUPLE, QHXh14[1]/QHXh32[1] COUPLE, QHXh13[1]/QHXh33[1] COUPLE, QHXh14[1]/QHXh34[1] COUPLE, QHXh13[1]/QHXh35[1] COUPLE, QHXh14[1]/QHXh36[1] COUPLE, QHXh13[1]/QHXh37[1] COUPLE, QHXh14[1]/QHXh38[1] COUPLE, QHXh13[1]/QHXh39[1] COUPLE, QHXh14[1]/QHXh40[1] COUPLE, QHXh13[1]/QHXh41[1] COUPLE, QHXh14[1]/QHXh42[1] COUPLE, QHXh13[1]/QHXh43[1] COUPLE, QHXh14[1]/QHXh44[1] COUPLE, QHXh13[1]/QHXh45[1] COUPLE, QHXh14[1]/QHXh46[1] 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 !MIGRAD, TOL=1.E-20 ENDMATCH VALUE, KQUM1,KQUM2,KQUM3,KQUM4 !,KQHX13 TWISS, BETA0=TWSSmh, SAVE, TAPE="twiss.tape" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=QE35[1]/HXR46END, & STYLE=100, SPLINE, FILE="match" PLOT, TABLE=TWISS, HAXIS=S, VAXIS=BETX,BETY, RANGE=HXR13BEG/HXR46END, & STYLE=100, SPLINE, FILE="match" ENDSUBROUTINE ! ------------------------------------------------------------------------------! ------------------------------------------------------------------------------ ! NOTE: keep Yuri's energy-scalable values for these quads ... don't rematch ! matching TCAV per Yuantao's requirements MDUMPh : SUBROUTINE !match dumpline in HXR BEAM, ENERGY=Ef SET, SETCUS, CUSOFF SETK2cuH USE, (LTUH,HXRUND,DUMPLINE), RANGE=MM1/#E VALUE, KQHX46,KQUE1,KQUE2,KQDMP MATCH, BETA0=TWSSmh VARY, KQHX46, STEP=1.E-6, LOWER=-0.9874*(15/Ef), UPPER= 0.9874*(15/Ef) ! 41.5 kG @ 15 GeV VARY, KQUE1, STEP=1.E-6, LOWER=-0.1464*(15/Ef), UPPER= 0.1464*(15/Ef) ! 40.3 kG @ 15 GeV VARY, KQUE2, STEP=1.E-4, LOWER=-0.1464*(15/Ef), UPPER= 0.1464*(15/Ef) !-40.3 kG @ 15 GeV VARY, KQDMP, STEP=1.E-4, LOWER=-0.1789*(15/Ef), UPPER= 0.1789*(15/Ef) !-38.5 kG @ 15 GeV CONSTR, QUE1[1], BETX<200 CONSTR, QUE2[1], BETY<285 CONSTR, DUMPFACE, DY<1.0 WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 COUPLE, MTCX/OTRDMP, MUX=90/360 !LMDIF, TOL=1.E-20, CALLS=10000 !MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQHX46,KQUE1,KQUE2,KQDMP PRINT, FULL TWISS, BETA0=TWSSmh, SAVE=TWISS, TAPE="twiss.tape" 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="match" 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="match" ENDSUBROUTINE MDUMPh_alt : SUBROUTINE !match dumpline in HXR BEAM, ENERGY=Ef SET, SETCUS, CUSOFF SETK2cuH USE, (LTUH,HXRUND,DUMPLINE), RANGE=MM1/#E VALUE, KQHX46,KQHX47,KQUE1,KQUE2,KQDMP MATCH, BETA0=TWSSmh VARY, KQHX47, STEP=1.E-6, LOWER=-0.9874*(15/Ef), UPPER= 0.9874*(15/Ef) ! 41.5 kG @ 15 GeV VARY, KQUE1, STEP=1.E-6, LOWER=-0.1464*(15/Ef), UPPER= 0.1464*(15/Ef) ! 40.3 kG @ 15 GeV VARY, KQUE2, STEP=1.E-4, LOWER=-0.1464*(15/Ef), UPPER= 0.1464*(15/Ef) !-40.3 kG @ 15 GeV !VARY, KQDMP, STEP=1.E-4, LOWER=-0.1789*(15/Ef), UPPER= 0.1789*(15/Ef) !-38.5 kG @ 15 GeV CONSTR, ENDDMPH_1, BETX=39.414, ALFX=1.653, BETY=275.158, ALFY=1.473 !CONSTR, QUE1[1], BETX>150, BETX<200 !CONSTR, QUE2[1], BETY>250, BETY<285 !CONSTR, DUMPFACE, DY<1.0 WEIGHT, BETX=0, BETY=0, ALFX=0, ALFY=0, DX=0, DPX=0, DY=0, DPY=0 COUPLE, MTCX/OTRDMP, MUX=90/360 !!LMDIF, TOL=1.E-20, CALLS=10000 !!MIGRAD, TOL=1.E-20, CALLS=10000 ENDMATCH VALUE, KQHX47,KQUE1,KQUE2,KQDMP PRINT, FULL TWISS, BETA0=TWSSmh, SAVE=TWISS, TAPE="twiss.tape" 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="match" 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="match" ENDSUBROUTINE ! ============================================================================== ! R56 checks ! ------------------------------------------------------------------------------ MR56CUSXR : SUBROUTINE BEAM, ENERGY=Ef SET, SETCUS, CUSON USE, CUSXR MATCH, BETA0=TWSCuS1 WEIGHT, I2=1, I5=1, I5I2=1 TMATRIX, #S/#E, TM(6,6,6)=0 RMATRIX, #S/#E, RM(6,6)=1 ENDMATCH ENDSUBROUTINE ! ============================================================================== ! Cu linac quad strengths ! ------------------------------------------------------------------------------ L2K1in : SUBROUTINE !initialize L2 quad K1s to FODO cell values SET, KQ21401, KQL2 SET, KQ21501, -KQL2 SET, KQ21601, KQL2 SET, KQ21701, -KQL2 SET, KQ21801, KQL2 SET, KQ21901, -KQL2 SET, KQ22201, KQL2 SET, KQ22301, -KQL2 SET, KQ22401, KQL2 SET, KQ22501, -KQL2 SET, KQ22601, KQL2 SET, KQ22701, -KQL2 SET, KQ22801, KQL2 SET, KQ22901, -KQL2 SET, KQ23201, KQL2 SET, KQ23301, -KQL2 SET, KQ23401, KQL2 SET, KQ23501, -KQL2 SET, KQ23601, KQL2 SET, KQ23701, -KQL2 SET, KQ23801, KQL2 SET, KQ23901, -KQL2 SET, KQ24201, KQL2 SET, KQ24301, -KQL2 SET, KQ24401, KQL2 SET, KQ24501, -KQL2 SET, KQ24601, KQL2 ENDSUBROUTINE ! ------------------------------------------------------------------------------ L2K1out : SUBROUTINE !output L2 quad K1s VALUE, KQ21401 VALUE, KQ21501 VALUE, KQ21601 VALUE, KQ21701 VALUE, KQ21801 VALUE, KQ21901 VALUE, KQ22201 VALUE, KQ22301 VALUE, KQ22401 VALUE, KQ22501 VALUE, KQ22601 VALUE, KQ22701 VALUE, KQ22801 VALUE, KQ22901 VALUE, KQ23201 VALUE, KQ23301 VALUE, KQ23401 VALUE, KQ23501 VALUE, KQ23601 VALUE, KQ23701 VALUE, KQ23801 VALUE, KQ23901 VALUE, KQ24201 VALUE, KQ24301 VALUE, KQ24401 VALUE, KQ24501 VALUE, KQ24601 ENDSUBROUTINE ! ------------------------------------------------------------------------------ L3K1in : SUBROUTINE !initialize L3 quad K1s to FODO cell values SET, KQ25201, KQFL3 SET, KQ25301, KQDL3 SET, KQ25401, KQFL3 SET, KQ25501, KQDL3 SET, KQ25601, KQFL3 SET, KQ25701, KQDL3 SET, KQ25801, KQFL3 SET, KQ25901, KQDL3 SET, KQ26201, KQFL3 SET, KQ26301, KQDL3 SET, KQ26401, KQFL3 SET, KQ26501, KQDL3 SET, KQ26601, KQFL3 SET, KQ26701, KQDL3 SET, KQ26801, KQFL3 SET, KQ26901, KQDL3 SET, KQ27201, KQFL3 SET, KQ27301, KQDL3 SET, KQ27401, KQFL3 SET, KQ27501, KQDL3 SET, KQ27601, KQFL3 SET, KQ27701, KQDL3 SET, KQ27801, KQFL3 SET, KQ27901, KQDL3 SET, KQ28201, KQFL3 SET, KQ28301, KQDL3 SET, KQ28401, KQFL3 SET, KQ28501, KQDL3 SET, KQ28601, KQFL3 SET, KQ28701, KQDL3 SET, KQ28801, KQFL3 SET, KQ28901, KQDL3 SET, KQ29201, KQFL3 SET, KQ29301, KQDL3 SET, KQ29401, KQFL3 SET, KQ29501, KQDL3 SET, KQ29601, KQFL3 SET, KQ29701, KQDL3 SET, KQ29801, KQFL3 SET, KQ29901, KQDL3 SET, KQ30201, KQFL3 SET, KQ30301, KQDL3 SET, KQ30401, KQFL3 SET, KQ30501, KQDL3 SET, KQ30601, KQFL3 SET, KQ30701, KQDL3 SET, KQ30801, KQFL3 ENDSUBROUTINE ! ------------------------------------------------------------------------------ L3K1out : SUBROUTINE !output L3 quad K1s VALUE, KQ25201 VALUE, KQ25301 VALUE, KQ25401 VALUE, KQ25501 VALUE, KQ25601 VALUE, KQ25701 VALUE, KQ25801 VALUE, KQ25901 VALUE, KQ26201 VALUE, KQ26301 VALUE, KQ26401 VALUE, KQ26501 VALUE, KQ26601 VALUE, KQ26701 VALUE, KQ26801 VALUE, KQ26901 VALUE, KQ27201 VALUE, KQ27301 VALUE, KQ27401 VALUE, KQ27501 VALUE, KQ27601 VALUE, KQ27701 VALUE, KQ27801 VALUE, KQ27901 VALUE, KQ28201 VALUE, KQ28301 VALUE, KQ28401 VALUE, KQ28501 VALUE, KQ28601 VALUE, KQ28701 VALUE, KQ28801 VALUE, KQ28901 VALUE, KQ29201 VALUE, KQ29301 VALUE, KQ29401 VALUE, KQ29501 VALUE, KQ29601 VALUE, KQ29701 VALUE, KQ29801 VALUE, KQ29901 VALUE, KQ30201 VALUE, KQ30301 VALUE, KQ30401 VALUE, KQ30501 VALUE, KQ30601 VALUE, KQ30701 VALUE, KQ30801 ENDSUBROUTINE ! ============================================================================== ! COMMANDs ! ------------------------------------------------------------------------------ COMMENT SET, MU_L1, 75/360 ML1c USE, L1c SAVEBETA, TWSSL1, QFL1[1] TWISS SET, BET11, TWSSL1[BETX] SET, BET12, TWSSL1[BETY] VALUE, BET11,BET12 ENDCOMMENT COMMENT SET, MU_L2, 55.500/360 ML2c USE, L2c SAVEBETA, TWSSL2, QFL2[1] TWISS SET, BET21, TWSSL2[BETX] SET, BET22, TWSSL2[BETY] VALUE, BET21,BET22 ENDCOMMENT COMMENT !SET, MUX_L3, 36.175/360 !July 13, 2005 - set for n*360 deg BX24A to BX31A (TCAV3 -> 25-2d) SET, MUX_L3, 30.175/360 !July 13, 2008 - set for best WS28 45-deg mux phase advances SET, MUY_L3, 30.130/360 !July 13, 2005 - set for 3*90 deg TCAV3[1] to OTR30 (TCAV3 -> 25-2d) ML3c USE, L3c SAVEBETA, TWSSFL3, QFL3[1] SAVEBETA, TWSSDL3, QDL3[1] TWISS SET, BET31, TWSSFL3[BETX] SET, BET32, TWSSFL3[BETY] SET, BET33, TWSSDL3[BETX] SET, BET34, TWSSDL3[BETY] VALUE, BET31,BET32,BET33,BET34 ENDCOMMENT !ML0E !ML1E !ML2E !ML3E COMMENT !Cu linac matching routines !MTWSSC !MGSPEC !MHTRR56 !MBSYS !MED0 !MED2L1 !MSPEC !MBC1R56 !MBC1e !**!L2K1in !use ONLY when you're changing L2 phase advance per cell !MBC1L2 !ML22 !ML23 !ML24 !!L2K1out !MBC2R56 !MBC2i !**!L3K1in !use ONLY when you're changing L3 phase advance per cell !MBC2L3 !ML26 !ML27 !ML28 !ML29 !ML30 ENDCOMMENT ! L3/DMPS match !MCLTS !MCUSXR !MCULTUs !MSXR !MDUMPs ! CLTH_1/HXR match !MDL2B !MDL2M2 !MBSY !MBSY2 !MED2DL2 !MUNDh !MDUMPh !MDUMPh_alt ! ------------------------------------------------------------------------------ !MR56CUSXR !check R56 ! ------------------------------------------------------------------------------ RETURN