Integrated yields

From GOSIA

Jump to: navigation, search
(partial 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>
+
-
N
+
-
<\pre>
+
-
 
+
-
Refer to the page on the [particle_singles], which is a cross section given in the same output with the integrated yields.
+
-
 
+
-
 
+
-
Maybe the answer to your last question is obvious now, but...  I would not try
+
<math> N_\gamma = 10^{-30} N_p [N_A / A] Y_\gamma \epsilon_p \epsilon_\gamma \Delta\Omega_\gamma </math> (2),
-
to think of this as correcting for 4pi, since the gamma-ray angular
+
-
distribution has to be calculated as a function of theta, phi for 4pi in order
+
-
to make that correction.  If you really want Gosia to tell you what the total
+
-
yield would be for a 4pi array, you can change the *.gdt file entries after
+
-
running OP,GDET.  For Gammasphere the *.gdt file would look something like
+
-
<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.
-
  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 array--no
+
Refer to the page on the [[particle_singles | particle singles]], which is a cross section given in the same output with the integrated yields.
-
angular attenuation.  Then you would have in the YIELD column of the output the
+
-
p-gamma events where the particle hit the detector as you defined it and the
+
-
gamma ray is measured at all angles.  In this case the epsilon_gamma would
+
-
still be the absolute photopeak efficiency, but DeltaOmega_gamma would be 4*pi.
+
-
See pages 48 and 117 in the newest manual version, if you want to do this.  You
+
==Representing a <math>4\pi</math> Array in Gosia==
-
should be able to figure out that the first two entries in the file above
+
-
should be zero to represent a perfect 4pi array.
+
-
I hope that is all clear.  I can't ever figure out a brief way to explain
+
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]].
-
things.  Maybe after you figure this out, you can put it on the Wiki.  Your
+
-
questions come up often.
+

Latest revision as of 13:40, 15 June 2011

Retrieved from "https://www-user.pas.rochester.edu/~gosia/mediawiki/index.php/Integrated_yields"
Personal tools