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 b7cdc435d0
commit 13873bbd69
4 changed files with 319 additions and 287 deletions
Download Patch File Download Diff File Expand all files Collapse all files

38
Scripts/Final/SFPlus/matchInjector.madx
View File

@@ -3,9 +3,6 @@
L = 8.72;
alpha = 1;
beta = alpha*L;
Print,text='MATCH POINT: Laser Heater';
@@ -22,12 +19,15 @@ CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua230$end, alfx = TwissM1->ALFX;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua230$end, alfy = TwissM1->ALFY;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua230$end, betx = TwissM1->BETX;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua230$end, bety = TwissM1->BETY;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua130$end, betx < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua130$end, bety < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinlh01.mqua070$end, betx < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinlh01.mqua070$end, bety < 100;
LMDIF,CALLS=1000,TOLERANCE=1.e-21;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua130$end, betx < 65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb04.mqua130$end, bety < 65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb03.mqua230$end, betx < 65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb03.mqua230$end, bety < 65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb03.mqua130$end, betx < 65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinsb03.mqua130$end, bety < 65;
SIMPLEX,CALLS=1000,TOLERANCE=1.e-21;
LMDIF,CALLS=1000,TOLERANCE=1.e-21;
ENDMATCH;
@@ -42,10 +42,10 @@ VARY,NAME=sinbc01.MQUA050.k1,STEP=0.0001;
VARY,NAME=sinbc01.MQUA070.k1,STEP=0.0001;
VARY,NAME=sinbc01.MQUA090.k1,STEP=0.0001;
VARY,NAME=sinbc01.MQUA110.k1,STEP=0.0001;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc02$end, betx < 5;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc02$end, betx = 3;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc02$end, alfx = 0;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc01.mqua070$end, bety < 32;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc02$end, bety < 59;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc02$end, bety = 25;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinbc01$end, x = 0;
SIMPLEX,CALLS=1000,TOLERANCE=1.e-21;
ENDMATCH;
@@ -53,7 +53,7 @@ ENDMATCH;
twiss,SEQUENCE=swissfel,range=#s/sinbc02$end,beta0 = Twiss0;
plot,haxis=s,vaxis=betx,bety,range=#s/#e,colour=100;
sinma01.mqua020.k1=-0.1;
Print,text='MATCH POINT: Injector-TDS';
MATCH,SEQUENCE=SwissFEL,range=#s/sindi02$end,beta0 = Twiss0;
@@ -68,10 +68,10 @@ CONSTRAINT,SEQUENCE=SwissFEL,range=sindi02.mqua020$start,bety=TwissM2->BETY;
CONSTRAINT,SEQUENCE=SwissFEL,range=sindi02.mqua020$start,alfy=TwissM2->ALFY;
CONSTRAINT,SEQUENCE=SwissFEL,range=sindi02.mqua020$start,betx=TwissM2->BETX;
CONSTRAINT,SEQUENCE=SwissFEL,range=sindi02.mqua020$start,alfx=TwissM2->ALFX;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinma01$end, betx < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sinma01$end, bety < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sincb01$end, betx < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sincb01$end, bety < 100;
CONSTRAINT,SEQUENCE=SwissFEL,range=sincb01.mqua230$end, betx <65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sincb01.mqua230$end, bety < 65;
CONSTRAINT,SEQUENCE=SwissFEL,range=sincb01$end, betx < 50;
CONSTRAINT,SEQUENCE=SwissFEL,range=sincb01$end, bety < 50;
LMDIF,CALLS=1000,TOLERANCE=1.e-21;
SIMPLEX,CALLS=1000,TOLERANCE=1.e-21;
ENDMATCH;
@@ -81,10 +81,10 @@ twiss,SEQUENCE=swissfel,range=#s/sindi02$end,beta0 = Twiss0;
plot,haxis=s,vaxis=betx,bety,range=sinbc02$start/#e,colour=100;
s10cb02.mqua230.k1=0.2;
s10cb02.mqua430.k1=-0.2;
s10cb02.mqua230.k1=-0.2;
s10cb02.mqua430.k1=0.2;
muLin1=0.2;
muLin1=0.13;
Print,text='MATCH POINT: FODO - Linac 1';
match,SEQUENCE=swissfel,range=S10cb02$start/s10cb02$end;
Vary,name=S10CB02.MQUA230.k1,step=0.0001;
@@ -113,7 +113,7 @@ VARY,NAME=sindi02.mqua050.k1,STEP=0.0001,lower=-5.5,upper=5.5;
VARY,NAME=sindi02.mqua060.k1,STEP=0.0001,lower=-5.5,upper=5.5;
VARY,NAME=sindi02.mqua090.k1,STEP=0.0001,lower=-5.5,upper=5.5;
VARY,NAME=s10cb01.mqua230.k1,STEP=0.0001,lower=-5.5,upper=5.5;
VARY,NAME=s10cb01.mqua430.k1,STEP=0.0001,lower=-5.5,upper=5.5;
!VARY,NAME=s10cb01.mqua430.k1,STEP=0.0001,lower=-5.5,upper=5.5;
CONSTRAINT,SEQUENCE=swissfel,range=s10cb02$start,betx=nbetx;
CONSTRAINT,SEQUENCE=swissfel,range=s10cb02$start,alfx=nalfx;