User Tools

Site Tools


group_secar2.0

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Last revision Both sides next revision
group_secar2.0 [2018/09/12 15:42]
wagner [Impact of Lenght Change in Quadrupole 2]
group_secar2.0 [2018/09/13 16:48]
plastun
Line 2: Line 2:
  
 -- by Alex, Nadeesha and Louis -- -- by Alex, Nadeesha and Louis --
 +
  
 ===== DAY1===== ===== DAY1=====
Line 11: Line 12:
 {{ :dipole.png?300 |}}\\ {{ :dipole.png?300 |}}\\
 ====Understanding COSY Matrix output==== ====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:\\ Output:\\
 |1.0|0.0|0.0|0.0|-0.41| |1.0|0.0|0.0|0.0|-0.41|
Line 28: Line 30:
 {{ :reaction_kinematics.png?400 |}} {{ :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 γ 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**.\\ 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.\\ 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.\\ 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: Output for this beam passing the dipole:
 {{ :dipole46cr13.png?400 |}} {{ :dipole46cr13.png?400 |}}
  
 +----
 ===== DAY2 ===== ===== DAY2 =====
  
Line 77: Line 81:
 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.\\ 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===== ===== DAY3=====
-==== Impact of Lenght Change in Quadrupole 2 ====+==== 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.\\ The effective length of Q2 was reduced by 3% to 0.291m with the drift lenghts before and after increased by 0.0045m.\\
Line 84: Line 89:
 As an educated guess we changed the gradient of Q2 by 3% to compensate the length reduction. The resolving power recoverd to ~400.\\ 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 Q",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.+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 |}} {{ :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 ======
  
/srv/thewikis/JIOSS/data/pages/group_secar2.0.txt · Last modified: 2018/09/13 16:49 by plastun