| 1 | beam, particle='proton', e=0.9468, ex=4.08e-8, ey=4.08e-8;
|
|---|
| 2 |
|
|---|
| 3 | //all lengths in m
|
|---|
| 4 | //drift at the target-capture interface
|
|---|
| 5 | o0: drift, l=0.1;
|
|---|
| 6 | //length representing the distance from the physical start of the Gabor lens to the start of the hard edge focusing area
|
|---|
| 7 | o1: drift, l=0.15;
|
|---|
| 8 |
|
|---|
| 9 | //gl1: solenoid, l=0.857, b=1.2868;
|
|---|
| 10 | !gl1: solenoid, l=0.857, ks=2.290287;
|
|---|
| 11 | //gl2: solenoid, l=0.857, b=0.6671;
|
|---|
| 12 | !gl2: solenoid, l=0.857, ks=1.187325;
|
|---|
| 13 | //drift for the first RF cavity
|
|---|
| 14 | !o2: drift, l=0.5;
|
|---|
| 15 | //gl3: solenoid, l=0.857, b=0.8139;
|
|---|
| 16 | !gl3: solenoid, l=0.857, ks=1.448605;
|
|---|
| 17 | //drift for the energy selection;
|
|---|
| 18 | !o3: drift, l=2.55;
|
|---|
| 19 | //gl4: solenoid, l=0.857, b=1.3555;
|
|---|
| 20 | !gl4: solenoid, l=0.857, ks=2.412561;
|
|---|
| 21 | !o4: drift, l=2.9;
|
|---|
| 22 | //gl5: solenoid, l=0.857, b=1.3351;
|
|---|
| 23 | !gl5: solenoid, l=0.857, ks=2.376253;
|
|---|
| 24 |
|
|---|
| 25 |
|
|---|
| 26 | o0: drift, l=0.10;
|
|---|
| 27 | o1: drift, l=0.15;
|
|---|
| 28 | gl1: solenoid, l=0.857, ks=2.290287; ! anode radius of 3.65cm for all lenses
|
|---|
| 29 | o2: drift, l=0.15;
|
|---|
| 30 | o2a: drift, l=0.15;
|
|---|
| 31 | gl2: solenoid, l=0.857, ks=1.187325;
|
|---|
| 32 | o3: drift, l=0.15;
|
|---|
| 33 | //RF cavity 1
|
|---|
| 34 | s1rf1: drift, l=0.5;
|
|---|
| 35 | o5: drift, l=0.15;
|
|---|
| 36 | gl3: solenoid, l=0.857, ks=1.30567;
|
|---|
| 37 | o6: drift, l=0.15;
|
|---|
| 38 | o7: drift, l=2.55;
|
|---|
| 39 | s1rf2: drift, l=0.5;
|
|---|
| 40 | o8: drift, l=0.15;
|
|---|
| 41 | gl4: solenoid, l=0.857, ks=2.50424;
|
|---|
| 42 | o9: drift, l=0.15;
|
|---|
| 43 | o10: drift, l=0.1;
|
|---|
| 44 | o11: drift, l=0.15;
|
|---|
| 45 | gl5: solenoid, l=0.857, ks=0.411532;
|
|---|
| 46 | o12: drift, l=0.15;
|
|---|
| 47 | //place for octupole now mounted on the flange of gl5
|
|---|
| 48 | oct: drift, l=0.1;
|
|---|
| 49 | //RF cavity 2
|
|---|
| 50 | s1rf2: drift, l=0.5;
|
|---|
| 51 | //the extraction dipole is positioned after the RF2
|
|---|
| 52 |
|
|---|
| 53 | capture: line=(o0,o1,gl1,o2,o2a,gl2,o3);
|
|---|
| 54 | energyselection: line=(s1rf1,o5,gl3,o6,o7);
|
|---|
| 55 | matching: line=(o8,gl4,o9,o10,o11,gl5,o12, oct, s1rf2);
|
|---|
| 56 | ccapv4: line=(capture, energyselection, matching);
|
|---|
| 57 |
|
|---|
| 58 | // Assumed Initial Parameters
|
|---|
| 59 | //init_betx := 1.36e-4;
|
|---|
| 60 | //init_bety := 1.36e-4;
|
|---|
| 61 |
|
|---|
| 62 | //beam, particle=proton, sequence=ccapv6, energy = 0.95327231, NPART=1.05E11, sige=4.5e-4, ex=4.08e-8, ey=4.08e-8;
|
|---|
| 63 | Beam, particle = proton, sequence = injection, energy = 450,
|
|---|
| 64 | NPART=1.05E11, sige=4.5e-4 ;
|
|---|
| 65 | use, sequence = ccapv4;
|
|---|
| 66 |
|
|---|
| 67 | //select, flag = twiss;
|
|---|
| 68 | //twiss, sequence=ccapv4, betx=init_betx,bety=init_bety, save, file=full.twiss;
|
|---|
| 69 |
|
|---|
| 70 | select, flag = twiss, column=name,s,betx,alfx,mux,
|
|---|
| 71 | bety,alfy, muy,dx,dpx,dy;
|
|---|
| 72 | twiss, betx=0.002468,alfx=0.022969,bety=0.002468, alfy=0.022969, Dx=0.0, DPx=0.0, save, centre, file=twissnewi0.out;
|
|---|
| 73 |
|
|---|
| 74 | plot,file=injmatch0, haxis=s, vaxis=betx,bety, colour=100;
|
|---|
| 75 | //plot,file=s1optics, haxis=s, vaxis=alfx,alfy, colour=100;
|
|---|
| 76 |
|
|---|
