Differences

This shows you the differences between two versions of the page.

--- group_secar2.0 [2018/09/11 16:47]
wagner
+++ group_secar2.0 [2018/09/13 16:48]
wagner [Reaction Kinematics Calculation]
@@ Line 2: / Line 2: @@
 -- by Alex, Nadeesha and Louis --
 ===== DAY1=====
@@ Line 11: / Line 12: @@
 {{ :dipole.png?300 |}}\\
 ====Understanding COSY Matrix output====
-Command: OV 1 1 0 ; {order 1, phase space dim 1 INCLUDING ENERGY, # of parameters 0}\\
+Command:
+<code>OV 1 1 0 ; {order 1, phase space dim 1 INCLUDING ENERGY, # of parameters 0}</code>\\
 Output:\\
 |1.0|0.0|0.0|0.0|-0.41|
@@ Line 28: / Line 30: @@
 {{ :reaction_kinematics.png?400 |}}
-Reaction summary for <sup>45</sup>V+p→γ+<sup>46</sup>Cr, E<sub>kin</sub>(<sup>45</sup>V)=135 MeV\\
+Reaction summary for <sup>45</sup>V+p → γ+<sup>46</sup>Cr, E<sub>kin</sub>(<sup>45</sup>V)=135 MeV\\
 The maximum γ energy is 8.481 MeV. The minimum γ energy is 7.249 MeV.\\
 The maximum <sup>46</sup>Cr energy is **132.634 MeV**. The minimum <sup>46</sup>Cr energy is **131.402 MeV**. The maximum <sup>46</sup>Cr angle is **0.13 degrees**.\\
@@ Line 35: / Line 37: @@
 With the spot size on the target, designed to be ±0.5 mm (1σ), the anglular spread then is ±15 mrad.\\
 Adding the calculated spread from the reaction with the target of 0.13 deg = ±2.26 mrad we get ±17.26 mrad of angular spread.\\
-⇒COSY beam input: SB 0.0005 **0.01726** 0 0.0005 **0.01726** 0 0 0.0047 0 0 0 ;\\
+⇒COSY beam input:
+<code>SB 0.0005 **0.01726** 0 0.0005 **0.01726** 0 0 0.0047 0 0 0 ;</code>\\
 Output for this beam passing the dipole:
 {{ :dipole46cr13.png?400 |}}
+----
 ===== DAY2 =====
@@ Line 48: / Line 52: @@
 === Quadrupole strenght variations for emittance measurment ===
-|B<sub>Q7</sub> [T]	|X<sub>max</sub>[m]	|ap [m]	|Bρ* [Tm]	|L<sub>Q7</sub> [m]	|K|
+|B<sub>Q7</sub> [T]	|X<sub>max</sub>[m]	|ap [m]	|Bρ* [Tm]	|L [m]	|K|
-|0.02	|2.91E-03	|0.13	|0.8	|0.34	|0.065384615|
+|0.02	|2.91E-03	|0.13	|0.8	|6.085	|0.065384615|
-|0.01	|2.27E-03	|0.13	|0.8	|0.34	|0.032692308|
+|0.01	|2.27E-03	|0.13	|0.8	|6.085	|0.032692308|
-|0	|1.62E-03	|0.13	|0.8	|0.34	|0|
+|0	|1.62E-03	|0.13	|0.8	|6.085	|0|
-|-0.01	|1.20E-03	|0.13	|0.8	|0.34	|-0.032692308|
+|-0.01	|1.20E-03	|0.13	|0.8	|6.085	|-0.032692308|
-|-0.02	|9.08E-04	|0.13	|0.8	|0.34	|-0.065384615|
+|-0.02	|9.08E-04	|0.13	|0.8	|6.085	|-0.065384615|
-|-0.03	|6.23E-04	|0.13	|0.8	|0.34	|-0.098076923|
+|-0.03	|6.23E-04	|0.13	|0.8	|6.085	|-0.098076923|
-|-0.04	|9.54E-04	|0.13	|0.8	|0.34	|-0.130769231|
+|-0.04	|9.54E-04	|0.13	|0.8	|6.085	|-0.130769231|
-|-0.05	|1.60E-03	|0.13	|0.8	|0.34	|-0.163461538|
+|-0.05	|1.60E-03	|0.13	|0.8	|6.085	|-0.163461538|
-|-0.06	|2.26E-03	|0.13	|0.8	|0.34	|-0.196153846|\\
+|-0.06	|2.26E-03	|0.13	|0.8	|6.085	|-0.196153846|\\
 *Bρ the magnetic rigidity was extracted from COSY by comand CONS(CHIM)\\
@@ Line 64: / Line 68: @@
 {{ :plot.png?400 |}}
+Fit result for parabola X<sub>max</sub><sup>2</sup> = σ<sub>11</sub> = a·K<sup>2</sup> - 2ab·K + ab<sup>2</sup> + c is:\\
 a = 0.000362572706 \\
 b = 0.0834414129893 \\
 c = 3.506973286e-07\\
-X<sub>max</sub><sup>2</sup> = σ<sub>11</sub> = aK<sup>2</sup> - 2abK + ab<sup>2</sup> + c
-The consistency check for beam size calculation in Q7 shows: X_calculated = 3.13mm, X_cosy_read_out = 3.16mm at the Q7 exit.\\
-Emittance = 3.045387898578198e-07  π m rad= 0.3 π mm mrad.\\
-Start beam emittance= XX * AX= 0.2 π mm mrad, so the emittance grows from the start to the focal point.\\
-This rmittance grow is possibly due to the optics aberrations. Input COSY file does a 4th-order-calculation. Also the measured points don't fit the parabola ideally, i.e. the optics is non-linear.\\
+The consistency check for beam size calculation in Q7 shows:\\
+X_calculated = 3.13mm, X_cosy_read_out = 3.16mm at the Q7 exit.\\
+Emittance calculated from fit parameters ε= √(ac) / L<sup>2</sup> = 3.05E-7 π m rad = **0.305 π mm mrad**.\\
+Start beam emittance= XX · AX= 0.2 π mm mrad, so the emittance grows from the start to the focal point.\\
+This emittance grow is possibly due to the optics aberrations. That is because the input COSY file does a 4th-order-calculation. Also the measured points don't fit the parabola ideally, i.e. the optics is non-linear.\\
+----
+===== DAY3=====
+==== Impact of Lenght Change of Quadrupole 2 ====
+The effective length of Q2 was reduced by 3% to 0.291m with the drift lenghts before and after increased by 0.0045m.\\
+This led to a resolving power change from 445 to 102 for the standard particle of the script.
+As an educated guess we changed the gradient of Q2 by 3% to compensate the length reduction. The resolving power recoverd to ~400.\\
+This Q2 gradient was then used as the starting value for the fit routine, that minimized the inverse of the resolving power. Additionally the gradients of Q3 and Q4 were set as free parameters of the fit, too. The best fit value after 38 steps was 463, wich is even **5% better!** then the original resolving power.\\
+Off course then one has to check if the max angle and energy acceptance is still met. The fit values for Q2,Q3 and Q4 stayed close to the starting values so one has to keep in mind that the fit routine might not have discovered a global minimum but only a local one.
+{{ :fitted.png?400 |}}
+----
+===== DAY4=====
+==== Optics misalignments ====
+**Beam position. Longitudinal offset.**
+We populated "X-θ" phase-space with rays distributed on the border of ellipse by the loop:
+<code>
+LOOP PHI 0.0 360 5;
+SRXX:= XX*COS(PHI*3.14159/180);
+SRAX:= AX*SIN(PHI*3.14159/180);
+ENDLOOP;
+</code>
+Rays in another planes are on-axis:
+<code>
+SRYY:= 0.0;
+SRAY:= 0.0;
+SRDE:= 0.0;
+</code>
+The rays now look like this (not all the rays are displayed):
+{{ :rays.png?400 |}}
+Created loop for target position. Observed piece-wise curve for resolution vs target offset (meters) (see figure below) in case of 4th order calculation. If order reduced to the 1st, than curve is smooth, but maximum resolution is shifted.
+{{ :target_offset.png?450 |}}
+The resolution drops by 5% at the offsets (-5 mm, +7 mm).\\
+**Beam position. Transverse offset.**\\
+{{ :transv.png?450 |}}
+The resolution drops by 5% at the beam offsets (-0.8 mm, +0.5 mm). The offset is of the beam half-size order. \\
+**Beam size.**\\
+Change beam size with SCALING_FACTOR = (0.5..1.5):
+<code>
+SRXX:= XX*COS(PHI*3.14159/180)*SCALING_FACTOR
+</code>
+Resolution reasonably drops with larger beam size.
+{{ :beam_size.png?400 |}}
+**Dipole 4 position. Transvers offset**\\
+{{ :transvers_shift_b4.jpg?450 |}}\\
+**Dipole 4 position. Longitudinal offset**\\
+{{ :longitudinal_shift_b4.jpg?450 |}}\\
+**Dipole 4 position. Pitch**\\
+{{ :pitch_b4.jpg?400 |}}
+|Element|Variable |Range of ±5%|
+|Beam|transverse|-5..+7 mm|
+|Beam|longitudinal|-0.8..+0.5 mm|
+|Beam|size|-6%..+6%|
+|B4|transverse|-4..+1 mm|
+|B4|longitudinal|-10..+55 mm|
+|B4|pitch|-0.11..0.024 rad|
+==== Recoil tracking in SECAR ====
+**<sup>45</sup>V+p → γ+<sup>46</sup>Cr**
+Beam is set to 132MeV <sup>46</sup>Cr<sup>14+</sup>\\
+<code>
+RP 132.0 46*PARA(1) 14*PARA(2);
+{Design acceptances, depends on kinematics of a specific reaction!}
+XX:=0.0005;
+AX:=0.01726;
+YY:=0.0005;
+AY:=0.01726;
+DE:=0.0047;
+</code>
+Recoil beam fits the SECAR very well.
+{{ :46cr_recoil_in_secar.png?300 |}}
+** ±1 charge state **\\
+Added extra rays:
+<code>
+{Define special rays for different charge states}
+SR 0 -AX 0 0 0 0.0 0 1/14 2;
+SR 0 0 0 0 0 0.0 0  1/14 2;
+SR 0 AX 0 0 0 0.0 0 1/14 2;
+SR 0 -AX 0 0 0 0.0 0 -1/14 3;
+SR 0 0 0 0 0 0.0 0 -1/14 3;
+SR 0 AX 0 0 0 0.0 0 -1/14 3;
+{Beam 45V14+}
+SR 0 -AX 0 0 0 3/132 -1/46 0/14 5;
+SR 0 0 0 0 0 3/132 -1/46 0/14 5;
+SR 0 AX 0 0 0 3/132 -1/46 0/14 5;
+</code>
+{{ :different_beams.png?400 |}}
+One can see how Wien filter cleans the recoils from the projectile beam at FP2.
+----
+THANK YOU FOR VISITING OUR PAGE

JINA-CEE ION OPTICS SUMMER SCHOOL: JIOSS 2018

User Tools

Site Tools

Differences

Page Tools