Integrated yields

From GOSIA

Jump to: navigation, search
(more clean-up)
(separated p-gamma event calculation into a new subsection)
 
Line 14: Line 14:
After an integration, the cross section in mb is obtained by multiplying the yield by
After an integration, the cross section in mb is obtained by multiplying the yield by
-
the Ge solid angle in <math>sr</math> and dividing by the target thickness in <math>mg/cm^2</math>.  You can also calculate the absolute
+
the Ge solid angle in <math>sr</math> and dividing by the target thickness in <math>mg/cm^2</math>.   
-
p-gamma counts expected using the equation following equation 6.44b in the manual:
+
-
Detected p-gamma coincident events  
+
===Detected p-gamma coincident events===
-
<math>N_i = 10^{-30} [Q / qe] [N_A / A] Y(I_i-->I_f) \epsilon_p \epsilon_\gamma \Delta \Omega_\gamma</math>
+
The absolute p-gamma counts expected can be calculated using the equation following equation 6.44b in the manual:
-
If the absolute efficiency is known well, then it is possible to retrieve
+
<math>N_i = 10^{-30} [Q / qe] [N_A / A] Y(I_i-->I_f) \epsilon_p \epsilon_\gamma \Delta \Omega_\gamma</math>, (1)
 +
 
 +
where <math>Q</math> is the integrated beam current incident on the target during the experiment, <math>q</math> is the average charge state of the beam, <math>e</math> is the electron charge,
 +
 
 +
If the absolute efficiency is known well, then it is possible to reproduce
the actual counts measured in the photopeak by putting the efficiency into this
the actual counts measured in the photopeak by putting the efficiency into this
equation.
equation.
Line 31: Line 34:
<math>\epsilon_\gamma</math> is 0.1 to 0.15.  If the solid angle subtended by the crystal is  ~0.1 sr, then for one crystal <math>\epsilon_\gamma \Delta\Omega_\gamma ~ 0.01</math>.
<math>\epsilon_\gamma</math> is 0.1 to 0.15.  If the solid angle subtended by the crystal is  ~0.1 sr, then for one crystal <math>\epsilon_\gamma \Delta\Omega_\gamma ~ 0.01</math>.
-
If the laboratory setup and the EM matrix are accurately represented in the Gosia input, then the absolute counts can be obtained:
+
If the laboratory setup and the EM matrix are accurately represented in the Gosia input, then the absolute counts can be reproduced.  By appropriately combining the terms in equation (1) above, it can be rewritten
-
<pre>
+
<math> N_\gamma = 10^{-30} N_p [N_A / A] Y_\gamma \epsilon_p \epsilon_\gamma \Delta\Omega_\gamma </math> (2),
-
N...
+
 
-
</pre>
+
where <math>N_p</math> is the total number of ''incident'' particles, <math>N_A</math> is Avogadro's number, <math>A</math> is the mass number of the target species, <math>Y_\gamma</math> the yield given by Gosia (using OP,INTI), <math>\epsilon_p</math> is the particle-detection efficiency, <math>\epsilon_\gamma</math> is the <math>\gamma</math>-detector absolute photopeak efficiency, and <math>\Delta\Omega_\gamma</math> is the solid angle subtended by the Ge detector.
Refer to the page on the [[particle_singles | particle singles]], which is a cross section given in the same output with the integrated yields.
Refer to the page on the [[particle_singles | particle singles]], which is a cross section given in the same output with the integrated yields.
 +
==Representing a <math>4\pi</math> Array in Gosia==
-
 
+
A <math>4\pi</math> array can be represented in Gosia in a single calculation without summing the output of a large number of detectors.  See the page on [[four_pi_arrays | 4pi arrays]].
-
==Quickly representing a summed 4pi array==
+
-
 
+
-
Gosia can quickly integrate the p-gamma events over a <math>4\pi</math> Ge array without adding a detector at every laboratory position.  In order to do this, the output file 9 can be modified to make the first two attenuation coefficients 0 for all orders.  (File 9 is the output of OP,GDET, commonly called the *.gdt file.  For Gammasphere the *.gdt file would look something like
+
-
 
+
-
<pre>
+
-
  1
+
-
  25.0
+
-
  5.000E-02
+
-
0.  0.  0.9951
+
-
0.  0.  0.9854
+
-
0.  0.  0.9711
+
-
0.  0.  0.9521
+
-
0.  0.  0.9379
+
-
0.  0.  0.9013
+
-
0.  0.  0.8699
+
-
0.  0.  0.8349
+
-
</pre>
+
-
 
+
-
where the first two columns have been set to 0. to simulate a 4pi arrayIn this case the <math>\epsilon_\gamma</math> value would
+
-
still be the absolute photopeak efficiency, but <math>\Delta\Omega_\gamma</math> would be 4*pi.
+
-
 
+
-
See the manual entries on OP,GDET and "Gamma Detector Solid Angle Attenuation Factors" for more information on these attenuation coefficients.
+

Latest revision as of 13:40, 15 June 2011

Personal tools