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/14 00:33] gastis |
tripleg_group [2018/09/14 09:32] gao |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | This is the tutorial written by the ' | + | This is the tutorial written by the ' |
In this tutorial, we will explain step by step how we work on the project [[https:// | In this tutorial, we will explain step by step how we work on the project [[https:// | ||
- | |||
**Q3.** | **Q3.** | ||
We use the example from the wiki [[https:// | 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 40: | Line 41: | ||
**Q4.** | **Q4.** | ||
- | In step 4, we calculate the reaction kinematics for the 15O(a, | + | In step 4, we calculate the reaction kinematics for the 15O(a, |
{{ : | {{ : | ||
- | The energy spread in this case is about +-2% (maximum energy acceptance of SECAR is +-3.1%), and the angular spread is ~ +-10mrad (maximum angular acceptance of SECAR ~ +- 25mrad). At this energy the reaction products fit in our systems acceptance | + | 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 |
+ | |||
+ | The maximum | ||
+ | For reaching the maximum angular acceptance we need energies above 20MeV/u. | ||
+ | |||
+ | |||
+ | |||
+ | * all energies are in the lab system. | ||
+ | |||
+ | |||
+ | {{ : | ||
-------------- | -------------- | ||
Line 66: | Line 77: | ||
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). These equations were dirived | + | where P1, P2, P3 are the fit parameters that we found (see the image above). These equations were derived |
- | Using these equations and numbers we get: | + | Using these equations and 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**: | **Correction**: | ||
the equation s22 needs to be replaced by s22 = P3 -s11 -2*Ls12)/ | 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/ | 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.45mm*mrad | + | |
- | By increasing the quality of the fit the number should go closer to 0.3mm*mrad | + | By repeating the process using the above corrections we got: epsilon=0.47mm*mrad |
- | The nominal value of epsilon according to the .fox file is: epsilon=XX*AX= 0.001*0.0002 = 0.2 mm*mrad | + | |
- | You can perform these calculations using this spreadsheet: | + | By using more precise |
- | {{ : | + | |
+ | The nominal value of epsilon according to the ".fox" | ||
+ | |||
+ | You can perform these calculations using this spreadsheet: | ||