SwissFEL-Simulation/OpticsTools
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Optimizing SF+ matching

Browse Source
This commit is contained in:
reiche
2026-01-21 17:11:34 +01:00
parent 1bbf5c0ba5
commit aefed53812
4 changed files with 319 additions and 287 deletions
Download Patch File Download Diff File Expand all files Collapse all files

33
Scripts/Final/SFPlus/matchAthos.madx
View File

@@ -306,7 +306,7 @@ twiss,SEQUENCE=swissfel,range=satma01.mqua250$end/satun11$end,beta0=twissee;
plot,haxis=s,vaxis=betx,bety,range=#s/#e,colour=100;
scl=4;
scl=7;
satcb02.mqua230.k1=scl*0.2;
satcb02.mqua430.k1=-scl*0.2;
satcb03.mqua230.k1=scl*0.2;
@@ -340,6 +340,37 @@ ENDMATCH;
twiss,SEQUENCE=swissfel,range=s20sy02$start/satun11$end,beta0=twisssep;
plot,haxis=s,vaxis=betx,bety,range=satdi01$start/#e,colour=100;
Print,text='MATCH POINT: Athos XTCAV';
MATCH, SEQUENCE=swissfel,range=s20sy02$start/satbd01$end,beta0=twisssep;
VARY,NAME=satun26.mqua080.k1,STEP=0.0001,lower=-5.0,upper=5.0;
VARY,NAME=satun25.mqua080.k1,STEP=0.0001,lower=-5.0,upper=5.0;
VARY,NAME=satun24.mqua080.k1,STEP=0.0001,lower=-5.0,upper=5.0;
VARY,NAME=satun23.mqua080.k1,STEP=0.0001,lower=-5.0,upper=5.0;
CONSTRAINT,SEQUENCE=swissfel,range=satbd01.mqua010$start,betx=twissAT4->BETX;
CONSTRAINT,SEQUENCE=swissfel,range=satbd01.mqua010$start,bety=twissAT4->BETY;
CONSTRAINT,SEQUENCE=swissfel,range=satbd01.mqua010$start,alfx=twissAT4->ALFX;
CONSTRAINT,SEQUENCE=swissfel,range=satbd01.mqua010$start,alfy=twissAT4->ALFY;
!CONSTRAINT,SEQUENCE=swissfel,range=satun25.mqua080$start,betx < 80;
!CONSTRAINT,SEQUENCE=swissfel,range=satun25.mqua080$start,bety < 80;
LMDIF,CALLS=1000,TOLERANCE=1.e-21;
ENDMATCH;
! fixed values, given by Eduard
SATBD01.MQUA010.K1= 5.069789237265826e-01;
SATBD01.MQUA030.K1= -6.184665903305938e-01 ;
SATBD01.MQUA050.K1= 1.324854421128419e-01 ;
SATBD01.MQUA070.K1= 4.740280173415730e-01 ;
SATBD01.MQUA090.K1= -7.442648986155268e-01;
SATBD02.MQUA030.K1= -1.050806764409536e+00;
twiss,SEQUENCE=swissfel,range=s20sy02$start/satbd01$end,beta0=twisssep;
plot,haxis=s,vaxis=betx,bety,range=satdi01$start/#e,colour=100;
twiss,SEQUENCE=swissfel,range=#s/satun26$end,beta0=twiss0;
plot,haxis=s,vaxis=betx,bety,range=#s/#e,colour=100;