This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
tripleg_group [2018/09/13 23:19] gastis |
tripleg_group [2018/09/14 09:32] (current) gao |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | This is the tutorial written by the ' | + | This is the tutorial written by the ' |
- | In step 2, we chose the 15O+4He radiative capture reaction to study. We use the example from the wiki [[https:// | + | |
+ | In this tutorial, we will explain step by step how we work on the project [[https:// | ||
+ | |||
+ | |||
+ | **Q3.** | ||
+ | We use the example from the wiki [[https:// | ||
1, Find the lines as shown in the picture below: | 1, Find the lines as shown in the picture below: | ||
{{: | {{: | ||
+ | |||
In the picture above, each of the commands (DL, DP, MQ ...) represents an ion optics element. (For a complete reference of the commands used in the cosy script, one has to refer to the manual). Here in the step 3 of [[https:// | In the picture above, each of the commands (DL, DP, MQ ...) represents an ion optics element. (For a complete reference of the commands used in the cosy script, one has to refer to the manual). Here in the step 3 of [[https:// | ||
radius and a drift. So we modify this part of the script as shown in the picture below: | radius and a drift. So we modify this part of the script as shown in the picture below: | ||
Line 23: | Line 30: | ||
SB 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.025 0.0 0.0 0.0 | SB 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.025 0.0 0.0 0.0 | ||
| | ||
- | 4, After changing the beam energy spread, run the script again, we get the following result: | + | After changing the beam energy spread, run the script again, we get the following result: |
{{:: | {{:: | ||
Again, the result looks reasonable. | Again, the result looks reasonable. | ||
5, In step 3c, we also need to calculate with the first order matrix elements the position of ions using. This can be done by the equation that we mentioned above. We write it down here again: | 5, In step 3c, we also need to calculate with the first order matrix elements the position of ions using. This can be done by the equation that we mentioned above. We write it down here again: | ||
x1 = (x|x)x0 + (x|a)a0 + (x|y)y0 + (x|b|b0 + (x|l)l0 + (x|δk)δk0 | x1 = (x|x)x0 + (x|a)a0 + (x|y)y0 + (x|b|b0 + (x|l)l0 + (x|δk)δk0 | ||
- | 6, In step 4, we need to calculate the reaction kinematics. | + | |
- | some pictures to be added. | + | |
+ | |||
+ | |||
+ | **Q4.** | ||
+ | |||
+ | In step 4, we calculate the reaction kinematics | ||
+ | {{ : | ||
+ | |||
+ | The energy spread in this case is about +-2% (maximum energy acceptance of SECAR: +-3.1%), and the angular spread is ~ +-10mrad (maximum angular acceptance of SECAR: +- 25mrad). At this energy the reaction products fit in our system without any problem. | ||
+ | |||
+ | The maximum energy that we can go until we reach the maximum energy acceptance of SECAR is ~8.2 MeV/u. | ||
+ | For reaching the maximum angular acceptance we need energies above 20MeV/u. | ||
+ | |||
+ | |||
+ | |||
+ | * all energies are in the lab system. | ||
+ | |||
+ | |||
+ | {{ :gammas19ne.png? | ||
-------------- | -------------- | ||
| | ||
- | Tuesday, 11-Sep | ||
- | |||
| | ||
- | We completed part 6. We created | + | **Q5.** |
+ | |||
+ | Bibliography for the 15O(p, | ||
| | ||
+ | **Q6.** | ||
+ | | ||
{{ : | {{ : | ||
| | ||
| | ||
- | | + | We got our data by varying the field of the quadrupole Q7 in the simplified SECAR model " |
- | + | For getting | |
- | For getting | + | WRITE 6 ' |
- | WRITE 6 ' | + | |
- | right before the command " | + | |
- | + | ||
- | For calculating the emittance using the quadrupole variation method, we used the following equations: | + | |
| | ||
+ | For calculating the emittance we used the following equations: | ||
s11 = P1/ | s11 = P1/ | ||
s12 = (P2 - 2*P1/L)/ 2*L^2 | s12 = (P2 - 2*P1/L)/ 2*L^2 | ||
- | s22 = P3/L^2 - P1/ | + | s22 = P3/L^2 - P1/ |
- | | + | where P1, P2, P3 are the fit parameters that we found (see the image above). |
- | where P1, P2, P3 are the fit parameters that we found (see the image above). | + | |
- | + | ||
- | + | Using these equations and numbers we get: | |
- | By using these numbers we get: | + | |
- | + | ||
epsilon = sqrt (s11*s22 - s12*s12) = 2.19e-7 m*rad = 0.219 mm*mrad | epsilon = sqrt (s11*s22 - s12*s12) = 2.19e-7 m*rad = 0.219 mm*mrad | ||
| | ||
- | | + | **Correction**: |
- | | + | the equation s22 needs to be replaced by s22 = P3 -s11 -2*Ls12)/ |
- | | + | |
- | {{ : | + | For getting the correct parabola one should plot the quantity K=Gradient*Leff/ |
+ | |||
+ | By repeating the process using the above corrections we got: epsilon=0.47mm*mrad | ||
+ | |||
+ | By using more precise fit parameters in the calculation (excel returns rounded values) the number should go closer to 0.3mm*mrad | ||
+ | |||
+ | The nominal value of epsilon according to the ".fox" | ||
+ | |||
+ | You can perform these calculations using this spreadsheet: | ||
- | | ||
- | | + | **Q7.** |
1, Change the effective length of Q2, the resolving power degraded slightly. | 1, Change the effective length of Q2, the resolving power degraded slightly. | ||
2, To restore the optimum resolving power, change the field strength of Q2 accordingly. We don't do it by hand, we use the FIT tookit to do that for us. | 2, To restore the optimum resolving power, change the field strength of Q2 accordingly. We don't do it by hand, we use the FIT tookit to do that for us. |