Split-Pole Spectrograph: Difference between revisions
(Created page with "== SABRE == thumb|right|Installing particle shield on SABER (photo taken on May 5, 2022)") |
|||
(139 intermediate revisions by 7 users not shown) | |||
Line 1: | Line 1: | ||
[[File:Sideview of the split-pole spectrometer.jpg|none| 800px | Side view of the spilt-pole spectrometer]] | |||
[[File:SABER installing particle shield.png|thumb|right|Installing particle shield on | {| align="right" style="width:30%;" | ||
|- | |||
| [[File:SPS Pictue Annotated.png|thumb|Annotated picture of the SE-SPS, An plain picture is here : [[:File:SPS Magnet.png]]]] || [[File:SPS Picture ray.png|thumb| SE-SPS COSY simulation. An plain picture is here [[:File:SPS Sketch With Cosy.png]]]] | |||
|} | |||
The '''Super Enge Split-Pole Spectrograph''' | |||
<ref name="Enge1979">H.A. Enge, NIM '''162''', 161 (1979) https://doi.org/10.1016/0029-554X(79)90711-0</ref> | |||
<ref>H. A. Enge, NIM '''187''', 1 (1981) https://doi.org/10.1016/0029-554X(81)90465-1 </ref> | |||
is a magnet spectrometer to measure the spectrum of nuclear reactions. The concept and design were developed by [https://en.wikipedia.org/wiki/Harald_A._Enge Harald A. Enge]<ref>J. E. Spencer and H. A. Enge, NIM '''49''', 181 (1967) https://doi.org/10.1016/0029-554X(67)90684-2 </ref> at 1967, aimed to have a broad-momentum range spectrograph with <math> p_{max}/p_{min} \approx 2.8 </math> or <math> E_{max}/E_{min} \approx 8 </math>. The spectrometer was originally located at the Wright Nuclear Structure Laboratory (closed at 2013), at Yale University. It was moved to FSU in the fall of 2013. It consists of a reaction chamber, a '''split-pole magnetic spectrograph''', a '''position-sensitive ionization drift chamber''', and a '''plastic scintillator'''. It has an angular acceptance of 128 msr (vertical ±40 mrad, horizontal ±80 mrad). The maximum B-field is 1.63 T with a radius of curvature from 511 mm to 920 mm. The mean radius is 600 mm. The advantage of the split-pole instead of a single-pole magnet is the aberration (x|θ<sup>2</sup>) and (x|φ<sup>2</sup>) are almost zero <ref name="Enge1979" />. | |||
The '''Super''' Enge Split-Pole Spectrograph is an upgrade of the Yale Enge SPS. The major change is the redesign of the backward silicon detector array to the [[Split-Pole_Spectrograph#SABRE|SABRE]]. | |||
= SPS Experiment Guide = | |||
[[Media:SPS_Experiment_Guide.pdf]] | |||
= SPS Operating Procedures = | |||
I created this section as a place to store procedures for the chamber swaps, however, I expect there are other things we might want to document here. -p | |||
* [[Target Chamber Swaps]] | |||
= Online Simulator = | |||
https://fsunuc.physics.fsu.edu/SPS/ | |||
There was another simulator called [https://github.com/gwm17/websps websps] developed by Gordon. But I am unable to revive it after the server upgraded from Ubuntu 20.04 to Ubuntu 22.04, so I write a new one using javascript. The missing part is the energy loss calculation. | |||
= Magnet = | |||
[[File:Design of a Split-pole spectrograph.png|thumb|Design of a Split-pole spectrograph. Take from Ref. <ref name="Enge1979"/>]] | |||
The primary goal of a spectrograph is resolving momentum. A general discussion of magnetic spectrographs can be found at Ref. <ref name="Enge1979" /> and Ref. <ref> H. A. Enge, Physics Today '''20''', 65 (1967) https://doi.org/10.1063/1.3034401 </ref>. There are many designs from a simple single dipole to a combination of multiple dipoles and quadrupoles. | |||
The SPS magnet was designed for a large solid angle, large resolving power, and correction of kinematic broadening. Using two-directional focusing and second-order focusing spectrograph can achieve a large sold angle and resolving power. Second-order focusing means the second-order terms in the acceptance angles vanish, i.e. no aberration. | |||
The SPS contains 2 separate poles enveloped by a single coil. The split provides second-order double focusing over a broad range of momenta. The magnet can be rotated from 0 to 55 degrees in the lab. The magnetic field has an upper limit of 1.63 T (or 16.3 kG). | |||
{| class="wikitable" | |||
|+ This table is taken from B.P. Kay Ph.D thesis (2007) | |||
! Property !! Symbol !! Value | |||
|- | |||
| Orbital radius || <math> \rho </math> || 511 to 920 mm | |||
|- | |||
|Resolving power || <math> p/\Delta p </math> || 1st order of <math> (x|\theta) </math>4290 (at <math>\theta = \pm 80 </math> mrad) | |||
|- | |||
|rowspan="2"| Acceptance || Horizontal || 160 mrad | |||
|- | |||
| Vertical || 80 mrad | |||
|- | |||
| Dispersion || <math> D=(x|\delta) </math> || 1.96 | |||
|- | |||
|rowspan="2" | Magnification || <math> M_x = (x|x)</math> || 0.39 | |||
|- | |||
| <math> M_y = (y|y)</math> || 2.9 | |||
|- | |||
| Maximum field || <math> B </math> || 1.63 T | |||
|} | |||
== Transfer matrix and COSY INFINITY simulation == | |||
[[File:Schematics of ion-beam optical element.png|thumb|An illustration of the coordinate of an optical element. This is taken from H.A. Enge NIM 162, 161 (1979).]] | |||
{{Notice | need to fill up. Any 1st few orders transfer matrixes?}} | |||
The entrance coordinates of the beam are <math> x_1, y_1, \theta_1, \phi_1 </math> wiht momentum <math> \delta = p/p_0 </math>, and coordinate at exit are <math> x_2, y_2, \theta_2, \phi_2 </math>. The entrance and exit coordinates are related by | |||
<math> x_2 = f_x(x_1, y_1, \theta_1, \phi_1, \delta) </math> | |||
using Taylor expansion: | |||
<math> \frac{x_2}{\rho} = (x|x) \frac{x_1}{\rho} + (x|\theta) \theta_1 + (x|\delta) \delta + (x|\theta^2) \theta_1^2 + ... </math> | |||
In the above expansion, the term <math> (x|x) </math> is the '''magnification''' in the x-direction. <math> (x|\delta) </math> is the '''dispersion''', and <math> (x|\theta^2) </math> is '''aberration'''. The '''focal plane''' is the z-position that <math> (x|\theta) = 0 </math>, i.e. the exit <math> x_2 </math> does not depend on the entrance angle. | |||
=== Kinematic broadening === | |||
[[File:Kinematic correction of spectrometer.png|thumb|Kinematic correction of spectrometer. Taken from H. A. Enge NIM 162, 161 (1979)]] | |||
Kinematic broadening is the broadening of focus for the same reaction state. After a reaction, the angle and momentum of the recoil particle are related that the entrance angle <math> \theta_1 = f(\delta)</math> is a function of momentum. For each energy state, the relation between the angle and momentum is unique. For example, in a 2-body transfer reaction, the momentum vector is | |||
<math> ( p_x, p_z ) = ( k \sin(\theta), \gamma \beta \sqrt{m^2-k^2} + \gamma k \cos(\theta) ) </math> | |||
where <math> k, \theta</math> are the momentum and the scattering angle at the CM frame, <math> \gamma, \beta </math> are the Lorentz factor from Lab frame to CM frame, and <math> m </math> is the mass of the particle. All 5 coefficients are constant for a fixed energy state. And the different state is characterized by <math> k </math>. Defined the kinematic factor K: | |||
<math> K = \frac{1}{p} \frac{dp}{d\theta_1} = \frac{\beta E \sin(\theta)}{ k + \beta \sqrt{m^2 - k^2} \cos(\theta)} </math> | |||
The kinematic broadening can be corrected by shifting the focal plane by | |||
<math> \Delta z = - D M \rho K, D = (x|\delta), M = (x|x) </math> | |||
= Focal plane detector = | |||
[[File:Focal plane detector.png|400px|thumb|right|Front view of the opened camerabox. The SPS focal plane detector with the front window removed is at the bottom.]] | |||
[[File:Side Cross section view of the SPS focal plane detector.png|400px|thumb|Side Cross section view of the SPS focal plane detector. Taken from B.P. Kay Ph.D. thesis (2007).]] | |||
[[File:PID EDE annoteted.png|400px|thumb]] | |||
{{Notice | The drift ion chamber was repaired in summer 2018 }} | |||
The focal plane detector | |||
<ref> C. Marshal ''et. al'', IEEE Tran. Inst. and Meas. '''68''', 533 (2018) https://doi.org/10.1109/TIM.2018.2847938</ref> | |||
<ref name="Markham1975"> R. G. Markham and R. G. H. Robertson, NIM '''129''', 131 (1975) https://doi.org/10.1016/0029-554X(75)90122-6 </ref> | |||
consists of an ion drift chamber with a set of delay lines to detect the position of a particle along the focal plane and a plastic scintillator to detect the energy of the incoming particle. Using the energy loss of the particle through the ion chamber with the energy deposited in the scintillator, particles of different charges and masses can be identified. | |||
The typical pressure of the drift chamber is 70 to 300 Torr of isobutane gas [HC(CH3)3]. The pressure controls the density of the gas and affects the bias voltage, it further affects the drift velocity. | |||
{|class='wikitable' | |||
|+ Table of pressure and bias voltages. Data was taken from the Ph.D. thesis of Erin Good (2020) | |||
! Gas pressure (Torr) !! Anode bias (V) !! Cathode plate bias (V) | |||
|- | |||
| 70 || +1050 to +1035 || -550 to -500 | |||
|- | |||
| 80 || +1150 || -550 | |||
|- | |||
| 100 || +1250 || -600 | |||
|- | |||
| 110 || +1200 to 1320 || -620 to -600 | |||
|- | |||
| 125 || +1425 || -650 | |||
|- | |||
| 130 || +1360 || -725 | |||
|- | |||
| 150 || +1500 || -700 | |||
|- | |||
|} | |||
From bottom to top, the cathode plate, drift region (contains four biased field-shaping wire grids), Frisch grid (grounded), three anode wires, and pickup pads (which are with the delay lines). Electrons induced by any radiation will drift upward, pass the Frisch grid, are accelerated by the anodes, and hit the pickup pads. The pickup pads are strips with 45° against the anode wires, almost parallel to the particle trajectories<ref name="Markham1975" />. Each pickup strip is 0.09" (2.286 mm) wide and 1.4" (35.56 mm) long, and spaced 0.01" (0.245 mm). A total of 440 lead-coated copper strips with a 5 ns delay per strip results in a nominal total delay of 2.2 μs. Every 10 strips share a delay chip. The position of the hit position can then be determined by the time difference at the end of the delay line. | |||
There are two position-sensitive delay lines (separated by 42.8625 mm) in the focal plane detector. By reconstructing the particle trajectory using the position information of both delay lines, the resolution can be enhanced by correcting for the kinematic shift of the reaction. | |||
After passing the drift chamber, the particles will be stopped and detected in a plastic scintillator with a photomultiplier tube (PMT) at each end. Together with the energy loss, obtained by the cathode in the drift chamber, a ΔE-E particle identification can be done. | |||
== Outline of the algorithm == | |||
There are 9 readouts channels from the focal plane detector: | |||
{| class="wikitable" | |||
! readout !! type of signal | |||
|- | |||
| cathode || energy loss | |||
|- | |||
| Front delay line Left || timing | |||
|- | |||
| Front delay line Right || timing | |||
|- | |||
| Front anode || energy loss | |||
|- | |||
| Rear delay line Left || timing | |||
|- | |||
| Rear delay line Right || timing | |||
|- | |||
| Rear anode || energy loss | |||
|- | |||
| PMT Left || energy loss | |||
|- | |||
| PMT Right || energy loss | |||
|} | |||
The PID is usually using one of the PMT energy and either the cathode or anode energy. | |||
The coordinate at the Focal Plane is the conventional one, where z-axis is perpendicular to the focal plane detector, y-axis is the vertical, and x-axis is the z-axis cross y-axis. The positions of the front and Rear planes are constructed by the timestamp. Suppose the timestamp is in ns. | |||
<math> | |||
x_1 = \frac{t_{FL} - t_{FR}}{2} \frac{1}{2.1} ; ~~~ x_2 = \frac{t_{BL} - t_{BR}}{2} \frac{1}{1.98} | |||
</math> | |||
The position at the center of the focal plane is | |||
<math> | |||
x_{avg} = x_1 + x_2 | |||
</math> | |||
However, for different reactions, there is a z-offset, so that the focal plane is shifted | |||
<math> | |||
x_{avg} = \left(\frac{1}{2} - \frac{z_{o}}{D} \right) x_1 + \left(\frac{1}{2} + \frac{z_{o}}{D} \right) x_2 | |||
</math> | |||
where <math> D = 42.8625~[\textrm{mm}] </math> is the distance between the front and rear delay lines. | |||
=== Calculation of the z-offset === | |||
The z-offset depends on the reaction, the angle <math> \theta </math>, and the magnetic field of the spilt-pole. Suppose we know the KE <math> T_b </math> and momentum <math> P </math> of the ejectile or the interested particle that goes into the split-pole, The reaction is denoted as A(a,b)B, where a is the beam (projectile), A is the target, B is the heavy recoil (residual),and b is the recoil (ejectile). | |||
<math> P^2 = (m_b+T_b)^2 - m_b^2 = T_b(T_b + 2m_b)~[\textrm{MeV/c}^2]</math> | |||
<math> \rho = \frac{P}{c Z B}~[\textrm{m}] </math>, where <math> c = 299.792458 </math>, <math>Z</math> is the charge number, and <math>B</math> is the magnetic field in Tesla. | |||
<math> k = \frac{v \sin(\theta)}{ m_b + m_B - v \cos(\theta)},~~~ v = \sqrt{\frac{m_a m_B T_a}{T_b}}</math> | |||
<math> z_o = -\rho \delta_x M_x k [\textrm{m}]</math> | |||
where <math> \delta_x = 1.96 </math> is the x-dispersion and <math> M_x = 0.39</math> is the x-magnification. | |||
=== Simplification of the z-offset calculation (IW) === | |||
The Kinematic Factor k can be re-written as a function of the beam and ejectile momenta: | |||
<math> k = \frac{(p_a/p_b) \sin(\theta)}{ 1 + m_B/m_b - (p_a/p_b) \cos(\theta)} </math> | |||
Inserting that into the offset formula and using <math> p_b = B \rho q_b </math> gives | |||
<math> z_o = - \delta_x M_x \frac{(p_a/B q_b) \sin(\theta) }{ 1 + m_B/m_b - (p_a/(B \rho q_b)) \cos(\theta)} </math>, | |||
which only depends on the beam momentum and masses and the spectrograph B and <math> \theta </math>. | |||
The second order correction (tilt of focal plane) is included in the denominator through the <math> \rho </math> dependence. | |||
== Kinematic correction of the focal plane == | |||
As pointed out before, the [[Split-Pole_Spectrograph#Kinematic_broadening | kinematic broadening]] can be corrected. In SPS, the dispersion D is 1.96, magnification is 0.39. | |||
{|class='wikitable' | |||
| style="width: 400px;"| [[File:AnnotatedFocalPlaneRay.png | 400px|frameless| ]] | |||
| style="width: 400px;"| [[File:FPShift.gif|frame|]] | |||
|- | |||
|Simulated rays near the focal plane. || An animation on the shift of the focal panel. An optimum is reached at FP = -42 mm. | |||
|} | |||
[[File:XavgDiagram.png|thumb| construction of Xavg (X-average) on the virtual focal plan (a liner plane in this case). Need to redraw the picture, the Y-axis should be Z-axis, and it should be rotated 180 degree, so the particle is from bottom to top.]] | |||
A parallel shift of the focal plane maybe not be enough. Suppose the best focal plan is given by a function <math> z = f(x) </math>. The 2 positions extracted from the front and rear delay lines are <math> x_1, x_2 </math>, and the distance between the front and rear delay lines is <math> d</math>. The X-avg is the solution of the equation: | |||
<math> f(x) = \frac{x_2 - x_1}{d} \left( x - \frac{x_2 + x_1}{2} \right) </math> | |||
For a linear tilted plane <math> f(x) = m x + z_0 </math>, the X-avg is | |||
<math> X_{avg} = \frac{x_1^2 - x_2^2 - 2d z_0 }{2 (d m + x_1 - x_2) } </math> | |||
<!--[[File:TwoBodyKinematics.png|thumb]]--> | |||
= SABRE = | |||
[[File:SABER installing particle shield.png|thumb|right|Installing particle shield on SABRE (photo taken on May 5, 2022)]] | |||
SABRE is a '''S'''ilicon '''A'''rray for '''B'''ranching '''R'''atio '''E'''xperiments | |||
<ref> E. C. Good ''et. al'', NIM A '''1003''', 165299 (2021) https://www.sciencedirect.com/science/article/pii/S0168900221002837</ref> | |||
with the SPS. Its predecessor is the Yale Lamp Shade Array (YLSA). SABRE sits at backward angles from the target and covers roughly 30% of 4π. SABRE has both thick and thin dead-layer detectors, with the thin dead-layer detectors capable of reaching ~200 keV thresholds for protons and deuterons. | |||
= CeBrA = | |||
[[File:CeBrA array diagram.png|thumb| Solid Works drawing of the fully planned array, which will consist of 13 CeBr<sub>3</sub> detectors.]] | |||
The Cerium Bromide Array (CeBrA) <ref> A. L. Conley ''et al.'', NIM A '''1058''' (2023) 168827 https://doi.org/10.1016/j.nima.2023.168827 </ref> is a γ-ray detector array designed to be used in conjunction with the SE-SPS. Comprised of low-background CeBr<sub>3</sub> scintillators with Hamamatsu photomultipliers. There are currently 7 commissioned detectors in the array with varying size crystals (2-1x1 inch, 4-2x2 inch, and 1-3x4 inch crystal detectors; a schematic for the full array on is shown on the right). The goal of CeBrA is to establish coincident events with the light-ions detected in the focal plane detector of the SE-SPS and the corresponding γ-rays from the excited recoiling nucleus, which are called particle-gamma coincidences. The scattering chamber used for CeBrA differs from the usual sliding-seal chamber that is used with the SE-SPS. With a hemisphere shape, it sits at a fixed 35<sup>o</sup> angle relative to the SE-SPS and allows for a more detailed study of the electromagnetic transitions from nuclei excited in reactions using the SE-SPS. | |||
{|class='wikitable' | |||
| style="width: 400px;"| [[File:Current_array.jpg | 400px|frameless|]] | |||
|- | |||
|Current array setup for CeBrA as it was used in the Summer 2023 REU experiments studying the <sup>52</sup>Cr(d,pγ)<sup>53</sup>Cr and <sup>34</sup>S(d,pγ)<sup>35</sup>S reactions, which was a follow up from the previous REU from 2022. | |||
|} | |||
= Repositories = | |||
https://github.com/sesps | |||
= Contact = | |||
* Jeff Blackmon mailto:blackmon@lsu.edu | |||
* Mark-Christoph Spieker mailto:mspieker@fsu.edu | |||
* Ingo Wiedenhoever mailto:iwiedenhoever@fsu.edu | |||
* <span style=color:red">who should be contacted? </span> | |||
= References = |
Latest revision as of 15:52, 20 June 2024
![]() Annotated picture of the SE-SPS, An plain picture is here : File:SPS Magnet.png |
![]() SE-SPS COSY simulation. An plain picture is here File:SPS Sketch With Cosy.png |
The Super Enge Split-Pole Spectrograph [1] [2] is a magnet spectrometer to measure the spectrum of nuclear reactions. The concept and design were developed by Harald A. Enge[3] at 1967, aimed to have a broad-momentum range spectrograph with or . The spectrometer was originally located at the Wright Nuclear Structure Laboratory (closed at 2013), at Yale University. It was moved to FSU in the fall of 2013. It consists of a reaction chamber, a split-pole magnetic spectrograph, a position-sensitive ionization drift chamber, and a plastic scintillator. It has an angular acceptance of 128 msr (vertical ±40 mrad, horizontal ±80 mrad). The maximum B-field is 1.63 T with a radius of curvature from 511 mm to 920 mm. The mean radius is 600 mm. The advantage of the split-pole instead of a single-pole magnet is the aberration (x|θ2) and (x|φ2) are almost zero [1].
The Super Enge Split-Pole Spectrograph is an upgrade of the Yale Enge SPS. The major change is the redesign of the backward silicon detector array to the SABRE.
SPS Experiment Guide
Media:SPS_Experiment_Guide.pdf
SPS Operating Procedures
I created this section as a place to store procedures for the chamber swaps, however, I expect there are other things we might want to document here. -p
Online Simulator
https://fsunuc.physics.fsu.edu/SPS/
There was another simulator called websps developed by Gordon. But I am unable to revive it after the server upgraded from Ubuntu 20.04 to Ubuntu 22.04, so I write a new one using javascript. The missing part is the energy loss calculation.
Magnet
![](/wiki/images/thumb/8/84/Design_of_a_Split-pole_spectrograph.png/300px-Design_of_a_Split-pole_spectrograph.png)
The primary goal of a spectrograph is resolving momentum. A general discussion of magnetic spectrographs can be found at Ref. [1] and Ref. [4]. There are many designs from a simple single dipole to a combination of multiple dipoles and quadrupoles.
The SPS magnet was designed for a large solid angle, large resolving power, and correction of kinematic broadening. Using two-directional focusing and second-order focusing spectrograph can achieve a large sold angle and resolving power. Second-order focusing means the second-order terms in the acceptance angles vanish, i.e. no aberration.
The SPS contains 2 separate poles enveloped by a single coil. The split provides second-order double focusing over a broad range of momenta. The magnet can be rotated from 0 to 55 degrees in the lab. The magnetic field has an upper limit of 1.63 T (or 16.3 kG).
Property | Symbol | Value |
---|---|---|
Orbital radius | 511 to 920 mm | |
Resolving power | 1st order of 4290 (at mrad) | |
Acceptance | Horizontal | 160 mrad |
Vertical | 80 mrad | |
Dispersion | 1.96 | |
Magnification | 0.39 | |
2.9 | ||
Maximum field | 1.63 T |
Transfer matrix and COSY INFINITY simulation
![]() | need to fill up. Any 1st few orders transfer matrixes? |
The entrance coordinates of the beam are wiht momentum , and coordinate at exit are . The entrance and exit coordinates are related by
using Taylor expansion:
In the above expansion, the term is the magnification in the x-direction. is the dispersion, and is aberration. The focal plane is the z-position that , i.e. the exit does not depend on the entrance angle.
Kinematic broadening
Kinematic broadening is the broadening of focus for the same reaction state. After a reaction, the angle and momentum of the recoil particle are related that the entrance angle is a function of momentum. For each energy state, the relation between the angle and momentum is unique. For example, in a 2-body transfer reaction, the momentum vector is
where are the momentum and the scattering angle at the CM frame, are the Lorentz factor from Lab frame to CM frame, and is the mass of the particle. All 5 coefficients are constant for a fixed energy state. And the different state is characterized by . Defined the kinematic factor K:
The kinematic broadening can be corrected by shifting the focal plane by
Focal plane detector
![]() | The drift ion chamber was repaired in summer 2018 |
The focal plane detector [5] [6] consists of an ion drift chamber with a set of delay lines to detect the position of a particle along the focal plane and a plastic scintillator to detect the energy of the incoming particle. Using the energy loss of the particle through the ion chamber with the energy deposited in the scintillator, particles of different charges and masses can be identified.
The typical pressure of the drift chamber is 70 to 300 Torr of isobutane gas [HC(CH3)3]. The pressure controls the density of the gas and affects the bias voltage, it further affects the drift velocity.
Gas pressure (Torr) | Anode bias (V) | Cathode plate bias (V) |
---|---|---|
70 | +1050 to +1035 | -550 to -500 |
80 | +1150 | -550 |
100 | +1250 | -600 |
110 | +1200 to 1320 | -620 to -600 |
125 | +1425 | -650 |
130 | +1360 | -725 |
150 | +1500 | -700 |
From bottom to top, the cathode plate, drift region (contains four biased field-shaping wire grids), Frisch grid (grounded), three anode wires, and pickup pads (which are with the delay lines). Electrons induced by any radiation will drift upward, pass the Frisch grid, are accelerated by the anodes, and hit the pickup pads. The pickup pads are strips with 45° against the anode wires, almost parallel to the particle trajectories[6]. Each pickup strip is 0.09" (2.286 mm) wide and 1.4" (35.56 mm) long, and spaced 0.01" (0.245 mm). A total of 440 lead-coated copper strips with a 5 ns delay per strip results in a nominal total delay of 2.2 μs. Every 10 strips share a delay chip. The position of the hit position can then be determined by the time difference at the end of the delay line.
There are two position-sensitive delay lines (separated by 42.8625 mm) in the focal plane detector. By reconstructing the particle trajectory using the position information of both delay lines, the resolution can be enhanced by correcting for the kinematic shift of the reaction.
After passing the drift chamber, the particles will be stopped and detected in a plastic scintillator with a photomultiplier tube (PMT) at each end. Together with the energy loss, obtained by the cathode in the drift chamber, a ΔE-E particle identification can be done.
Outline of the algorithm
There are 9 readouts channels from the focal plane detector:
readout | type of signal |
---|---|
cathode | energy loss |
Front delay line Left | timing |
Front delay line Right | timing |
Front anode | energy loss |
Rear delay line Left | timing |
Rear delay line Right | timing |
Rear anode | energy loss |
PMT Left | energy loss |
PMT Right | energy loss |
The PID is usually using one of the PMT energy and either the cathode or anode energy.
The coordinate at the Focal Plane is the conventional one, where z-axis is perpendicular to the focal plane detector, y-axis is the vertical, and x-axis is the z-axis cross y-axis. The positions of the front and Rear planes are constructed by the timestamp. Suppose the timestamp is in ns.
The position at the center of the focal plane is
However, for different reactions, there is a z-offset, so that the focal plane is shifted
where is the distance between the front and rear delay lines.
Calculation of the z-offset
The z-offset depends on the reaction, the angle , and the magnetic field of the spilt-pole. Suppose we know the KE and momentum of the ejectile or the interested particle that goes into the split-pole, The reaction is denoted as A(a,b)B, where a is the beam (projectile), A is the target, B is the heavy recoil (residual),and b is the recoil (ejectile).
, where , is the charge number, and is the magnetic field in Tesla.
where is the x-dispersion and is the x-magnification.
Simplification of the z-offset calculation (IW)
The Kinematic Factor k can be re-written as a function of the beam and ejectile momenta:
Inserting that into the offset formula and using gives
,
which only depends on the beam momentum and masses and the spectrograph B and . The second order correction (tilt of focal plane) is included in the denominator through the dependence.
Kinematic correction of the focal plane
As pointed out before, the kinematic broadening can be corrected. In SPS, the dispersion D is 1.96, magnification is 0.39.
![]() |
|
Simulated rays near the focal plane. | An animation on the shift of the focal panel. An optimum is reached at FP = -42 mm. |
A parallel shift of the focal plane maybe not be enough. Suppose the best focal plan is given by a function . The 2 positions extracted from the front and rear delay lines are , and the distance between the front and rear delay lines is . The X-avg is the solution of the equation:
For a linear tilted plane , the X-avg is
SABRE
SABRE is a Silicon Array for Branching Ratio Experiments [7] with the SPS. Its predecessor is the Yale Lamp Shade Array (YLSA). SABRE sits at backward angles from the target and covers roughly 30% of 4π. SABRE has both thick and thin dead-layer detectors, with the thin dead-layer detectors capable of reaching ~200 keV thresholds for protons and deuterons.
CeBrA
The Cerium Bromide Array (CeBrA) [8] is a γ-ray detector array designed to be used in conjunction with the SE-SPS. Comprised of low-background CeBr3 scintillators with Hamamatsu photomultipliers. There are currently 7 commissioned detectors in the array with varying size crystals (2-1x1 inch, 4-2x2 inch, and 1-3x4 inch crystal detectors; a schematic for the full array on is shown on the right). The goal of CeBrA is to establish coincident events with the light-ions detected in the focal plane detector of the SE-SPS and the corresponding γ-rays from the excited recoiling nucleus, which are called particle-gamma coincidences. The scattering chamber used for CeBrA differs from the usual sliding-seal chamber that is used with the SE-SPS. With a hemisphere shape, it sits at a fixed 35o angle relative to the SE-SPS and allows for a more detailed study of the electromagnetic transitions from nuclei excited in reactions using the SE-SPS.
![]() |
Current array setup for CeBrA as it was used in the Summer 2023 REU experiments studying the 52Cr(d,pγ)53Cr and 34S(d,pγ)35S reactions, which was a follow up from the previous REU from 2022. |
Repositories
Contact
- Jeff Blackmon mailto:blackmon@lsu.edu
- Mark-Christoph Spieker mailto:mspieker@fsu.edu
- Ingo Wiedenhoever mailto:iwiedenhoever@fsu.edu
- who should be contacted?
References
- ↑ 1.0 1.1 1.2 1.3 H.A. Enge, NIM 162, 161 (1979) https://doi.org/10.1016/0029-554X(79)90711-0
- ↑ H. A. Enge, NIM 187, 1 (1981) https://doi.org/10.1016/0029-554X(81)90465-1
- ↑ J. E. Spencer and H. A. Enge, NIM 49, 181 (1967) https://doi.org/10.1016/0029-554X(67)90684-2
- ↑ H. A. Enge, Physics Today 20, 65 (1967) https://doi.org/10.1063/1.3034401
- ↑ C. Marshal et. al, IEEE Tran. Inst. and Meas. 68, 533 (2018) https://doi.org/10.1109/TIM.2018.2847938
- ↑ 6.0 6.1 R. G. Markham and R. G. H. Robertson, NIM 129, 131 (1975) https://doi.org/10.1016/0029-554X(75)90122-6
- ↑ E. C. Good et. al, NIM A 1003, 165299 (2021) https://www.sciencedirect.com/science/article/pii/S0168900221002837
- ↑ A. L. Conley et al., NIM A 1058 (2023) 168827 https://doi.org/10.1016/j.nima.2023.168827