Rachel simulated yield output - Revision history

Hayes: /* Accuracy checking */

2012-01-28T16:02:10Z

Accuracy checking

← Older revision		Revision as of 16:02, 28 January 2012
Line 28:		Line 28:
	===Accuracy checking===		===Accuracy checking===

-	The predicted absolute counts (cross sections) for p-<math>\gamma</math> events were checked against perturbation calculations for weakly-coupled two-state systems where the perturbation method is valid. Both the total count rates and solid angle factors have been re-checked in version 2.0.8.beta and agree with the perturbation calculations where they are applicable.	+	The predicted absolute counts (cross sections) for p-<math>\gamma</math> events were checked against perturbation calculations for weakly-coupled two-state systems where the perturbation method is valid. Both the total count rates and solid angle factors have been re-checked in version 2.0.8.beta and agree with the perturbation calculations where they are applicable. Many tests of fits to simulated yield data without the optional random error scatter resulted in <math>\chi^2</math> values very close to zero, as expected.

-	The ''intrinsic'' efficiency has been estimated more accurately in version 2.0.8.beta. This will make a noticeable change in the "Raw counts" column.	+	The ''intrinsic'' efficiency has been estimated more accurately in version 2.0.8.beta. This will make a noticeable change in the "Raw counts" column.

	===Simulated data passed to Gosia===		===Simulated data passed to Gosia===

	Note that Rachel operates Gosia such that efficiency-corrected "incident gamma" counts are passed to Gosia, but first they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected, which usually did not result in accuracy errors.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.		Note that Rachel operates Gosia such that efficiency-corrected "incident gamma" counts are passed to Gosia, but first they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected, which usually did not result in accuracy errors.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.

Hayes: Major changes reflect the output and operation of version 1.

2012-01-28T16:00:05Z

Major changes reflect the output and operation of version 1.

← Older revision		Revision as of 16:00, 28 January 2012
Line 1:		Line 1:
-	''This page uses the output of Rachel version ~~2.0.8~~ beta.''	+	''This page uses the output of Rachel version 1, which is more intuitive than that of the beta versions.''

	The output to the terminal from a "simulated yields" calculation has the following format:		The output to the terminal from a "simulated yields" calculation has the following format:
		+
		+	The simulation output now uses the standard definition of efficiency, and the default efficiency curve for any detector can be adjusted by the user. The user will find prompts for the details of the planned beam run, and estimated absolute measured counts will be calculated and displayed. The display format looks like the following:

	<pre>		<pre>
-	Experiment 1 Detector 2	+	Experiment 1
-	~~Raw Count~~ Gamma **~~Est. intr. "Corrected"~~	+	Target excitation by Z,A = 54,136 at 517 MeV. Mean scattering angle = 45.0 deg.
-	Band Spin Band Spin ~~Counts~~ Error ~~Rate(Hz) E~~(keV) ~~effic.~~(1) **~~Counts~~ Error	+	Detector 1
-	----------------------------------------------------------------------------------------------------------------	+	Single crystal at theta, phi = 45.0, 45.0 deg. Solid angle = 0.06 sr.
-	gsb 18.0 ~~gsb 16~~.~~0 864 29~~ 0.~~002001 686~~.0 0.~~271 3266 198~~	+	This simulation does not include random gaussian scatter. Counts represent the cross sections calculated by Gosia.
-	~~gsb 16.~~0 gsb 14.0 ~~21056 145~~ 0.~~048741 658~~.7 0.~~277 78074 3941~~	+	Transition \| Gammas incident \| Fraction of \| Observed Counts
-	~~gsb 14~~.0 ~~gsb 12.~~0 ~~400895 633~~ 0.~~927999 626~~.6 0.~~284~~ ~~1360201 68044~~	+	\| on Ge detectors \| Gamma Incident \|
-	gsb 12.0 ~~gsb 10~~.~~0 5914064 2432 13~~.~~689963 579~~.7 0.~~296~~ ~~21145020 1057287~~	+	Initial Final \| \| Energy **Detector Gammas \|
-	~~gsb 10.~~0 ~~gsb 8~~.~~0 65755367 8109 152.21149 512~~.4 0.~~315 211214415 10560753~~	+	Band Spin Band Spin \| Counts Error \| (keV) Efficiency Detected \| Rate(Hz) Counts **Error
-	gsb 8.0 ~~gsb 6~~.0 ~~527950313 22977 1222~~.~~1072 426~~.3 0.~~343 1470121661~~ ~~73506111~~	+	-------------------------------------------------------------------------------------------------------------------------
-	~~gsb 6.~~0 ~~gsb 4~~.~~0 2866371774 53539 6635~~.~~1198 325~~.~~6 0~~.~~379 7435917854 371795919~~	+	gsb 1.5 gsb 0.5 4.662e+06 2.332e+05 100.0 0.00049 0.09962 1.075e+00 4.644e+05 6.815e+02
-	~~gsb 4.0 gsb~~ 2.~~0 8174712130~~ ~~90414 18922.944 213~~.5 0.~~382~~ ~~21045773464 1052288699~~	+	gsb 2.5 gsb 0.5 1.570e+07 7.848e+05 200.0 0.00179 0.36853 1.339e+01 5.784e+06 2.405e+03
-	~~gsb~~ 2.~~0 gsb 0~~.0 ~~706191594 26574 1634~~.~~7027~~ 93.3 0.~~078 8631491935 431574719~~	+	a 2.5 gsb 0.5 1.079e+04 5.647e+02 250.0 0.00189 0.38907 9.719e-03 4.199e+03 6.480e+01
		+	a 2.5 gsb 1.5 9.816e+02 7.705e+01 150.0 0.00135 0.27823 6.322e-04 2.731e+02 1.653e+01
		+	a 2.5 b 3.5 4.689e+05 2.376e+04 75.0 0.00015 0.0312 3.387e-02 1.463e+04 1.210e+02
	</pre>		</pre>

		+	The "Observed Counts" are calculated using the efficiency curve from the detector library or a detector created by the user. The "Detector Efficiency" <math>\epsilon</math> is the usual definition. Possion counting errors are based on the observed counts, while additional error on the ''total'' efficiency-corrected counts may be added by the user.

-	* The "Corrected Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.
-	* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute ''intrinsic'' efficiency for Gammasphere.
-	* The "Est. intr. effic." column gives the ''assumed'' absolute ''intrinsic'' efficiency <math>(\epsilon_{intr})</math> of the Ge detector used to calculate the value in the "Raw counts" column. For example, at 1333 keV, Gammasphere typically has an absolute efficiency of 9% for 100 installed detectors. Since the solid angle coverage of these 100 detectors is ~6.1 sr, the ''intrinsic'' efficiency is approximately <math>9% * 4\pi / 6.1 = 18.5%</math>.

	===Accuracy checking===		===Accuracy checking===
Line 29:		Line 31:

	The ''intrinsic'' efficiency has been estimated more accurately in version 2.0.8.beta. This will make a noticeable change in the "Raw counts" column.		The ''intrinsic'' efficiency has been estimated more accurately in version 2.0.8.beta. This will make a noticeable change in the "Raw counts" column.
-
-	~~===Estimating count rates in a planned experiment===~~
-
-	Since the "Corrected counts" column gives the absolute emitted particle-gamma counts into the particle- and gamma-detector solid angles (for 100% efficiency), these values can be combined with better known values of the particle and gamma efficiencies to predict more accurate "raw" count rates, if the approximate predicted raw counts are not accurate enough.

	===Simulated data passed to Gosia===		===Simulated data passed to Gosia===

-	Note that ~~the~~ "~~Corrected Counts~~" ~~are not sent to Gosia as they appear in the table.~~	+	Note that Rachel operates Gosia such that efficiency-corrected "incident gamma" counts are passed to Gosia, but first they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected, which usually did not result in accuracy errors.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.
-		+
-	~~Before the corrected~~ counts are ~~sent~~ to Gosia they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.	+

Hayes at 19:14, 29 June 2011

2011-06-29T19:14:53Z

← Older revision		Revision as of 19:14, 29 June 2011
Line 21:		Line 21:

	* The "Corrected Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.		* The "Corrected Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.
-	* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute efficiency for Gammasphere.	+	* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute ''intrinsic'' efficiency for Gammasphere.
	* The "Est. intr. effic." column gives the ''assumed'' absolute ''intrinsic'' efficiency <math>(\epsilon_{intr})</math> of the Ge detector used to calculate the value in the "Raw counts" column. For example, at 1333 keV, Gammasphere typically has an absolute efficiency of 9% for 100 installed detectors. Since the solid angle coverage of these 100 detectors is ~6.1 sr, the ''intrinsic'' efficiency is approximately <math>9% * 4\pi / 6.1 = 18.5%</math>.		* The "Est. intr. effic." column gives the ''assumed'' absolute ''intrinsic'' efficiency <math>(\epsilon_{intr})</math> of the Ge detector used to calculate the value in the "Raw counts" column. For example, at 1333 keV, Gammasphere typically has an absolute efficiency of 9% for 100 installed detectors. Since the solid angle coverage of these 100 detectors is ~6.1 sr, the ''intrinsic'' efficiency is approximately <math>9% * 4\pi / 6.1 = 18.5%</math>.

-	===~~Estimating count rates in a planned experiment~~===	+	===Accuracy checking===

-	~~Since the "Corrected counts" column gives the~~ absolute ~~emitted particle-gamma~~ counts ~~into the particle~~- ~~and~~ gamma-~~detector~~ solid ~~angles (for 100% efficiency), these values can be combined~~ with ~~better known values of~~ the ~~particle and gamma efficiencies to predict~~ more ~~accurate~~ "~~raw~~" count rates.	+	The predicted absolute counts (cross sections) for p-<math>\gamma</math> events were checked against perturbation calculations for weakly-coupled two-state systems where the perturbation method is valid. Both the total count rates and solid angle factors have been re-checked in version 2.0.8.beta and agree with the perturbation calculations where they are applicable.
		+
		+	The ''intrinsic'' efficiency has been estimated more accurately in version 2.0.8.beta. This will make a noticeable change in the "Raw counts" column.
		+
		+	===Estimating count rates in a planned experiment===

		+	Since the "Corrected counts" column gives the absolute emitted particle-gamma counts into the particle- and gamma-detector solid angles (for 100% efficiency), these values can be combined with better known values of the particle and gamma efficiencies to predict more accurate "raw" count rates, if the approximate predicted raw counts are not accurate enough.

	===Simulated data passed to Gosia===		===Simulated data passed to Gosia===

Hayes: fixed math formatting

2011-06-29T19:03:05Z

fixed math formatting

← Older revision		Revision as of 19:03, 29 June 2011
Line 22:		Line 22:
	* The "Corrected Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.		* The "Corrected Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.
	* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute efficiency for Gammasphere.		* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute efficiency for Gammasphere.
-	* The "Est. intr. effic." column gives the ''assumed'' absolute ''intrinsic'' efficiency <math>(\epsilon_{intr})</math> of the Ge detector used to calculate the value in the "Raw counts" column. For example, at 1333 keV, Gammasphere typically has an absolute efficiency of 9% for 100 installed detectors. Since the solid angle coverage of these 100 detectors is ~6.1 sr, the ''intrinsic'' efficiency is <math>~9% * 4 * \pi / 6.1 = 18.5%</math>.	+	* The "Est. intr. effic." column gives the ''assumed'' absolute ''intrinsic'' efficiency <math>(\epsilon_{intr})</math> of the Ge detector used to calculate the value in the "Raw counts" column. For example, at 1333 keV, Gammasphere typically has an absolute efficiency of 9% for 100 installed detectors. Since the solid angle coverage of these 100 detectors is ~6.1 sr, the ''intrinsic'' efficiency is approximately <math>9% * 4\pi / 6.1 = 18.5%</math>.

	===Estimating count rates in a planned experiment===		===Estimating count rates in a planned experiment===

Hayes at 18:50, 29 June 2011

2011-06-29T18:50:21Z

← Older revision		Revision as of 18:50, 29 June 2011
Line 1:		Line 1:
-	''This page ~~stub will be updated when a new version~~ of Rachel ~~is released~~.''	+	''This page uses the output of Rachel version 2.0.8 beta.''

-	The output to the terminal ~~following~~ a "simulated yields" calculation ~~will be similar to~~ the following ~~until the next release version where the formatting will be improved.~~	+	The output to the terminal from a "simulated yields" calculation has the following format:

	<pre>		<pre>
-	Experiment 1 Detector 4	+	Experiment 1 Detector 2
-	Raw Gamma ~~Sim~~. "~~Measured~~"	+	Raw Count Gamma **Est. intr. "Corrected"
-	Band Spin Band Spin Counts Error E(keV) ~~eps~~ Counts Error	+	Band Spin Band Spin Counts Error Rate(Hz) E(keV) effic.(1) **Counts Error
-	--------------------------------------------------------------------------------------------------------	+	----------------------------------------------------------------------------------------------------------------
-	gsb 6.0 gsb 4.0 ~~32055690~~ ~~5662~~ 325.6 0.~~251~~ ~~127940090 6397044~~	+	gsb 18.0 gsb 16.0 864 29 0.002001 686.0 0.271 3266 198
-	gsb 4.0 gsb 2.0 ~~27961778 5288~~ 213.5 0.~~204 137318753 6865987~~	+	gsb 16.0 gsb 14.0 21056 145 0.048741 658.7 0.277 78074 3941
-	gsb 2.0 gsb 0.0 ~~2352850 1534~~ 93.3 0.~~065 36044668 1802387~~	+	gsb 14.0 gsb 12.0 400895 633 0.927999 626.6 0.284 1360201 68044
		+	gsb 12.0 gsb 10.0 5914064 2432 13.689963 579.7 0.296 21145020 1057287
		+	gsb 10.0 gsb 8.0 65755367 8109 152.21149 512.4 0.315 211214415 10560753
		+	gsb 8.0 gsb 6.0 527950313 22977 1222.1072 426.3 0.343 1470121661 73506111
		+	gsb 6.0 gsb 4.0 2866371774 53539 6635.1198 325.6 0.379 7435917854 371795919
		+	gsb 4.0 gsb 2.0 8174712130 90414 18922.944 213.5 0.382 21045773464 1052288699
		+	gsb 2.0 gsb 0.0 706191594 26574 1634.7027 93.3 0.078 8631491935 431574719
	</pre>		</pre>

-	~~Three columns have ambiguous names:~~

-	* The "~~Measured~~ Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.	+	* The "Corrected Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.
	* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute efficiency for Gammasphere.		* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute efficiency for Gammasphere.
-	* The "~~Sim~~. ~~eps~~" column gives the assumed absolute efficiency <math>(\epsilon_{~~abs~~})</math> of the Ge detector used to calculate the value in the "Raw counts" column.	+	* The "Est. intr. effic." column gives the ''assumed'' absolute ''intrinsic'' efficiency <math>(\epsilon_{intr})</math> of the Ge detector used to calculate the value in the "Raw counts" column. For example, at 1333 keV, Gammasphere typically has an absolute efficiency of 9% for 100 installed detectors. Since the solid angle coverage of these 100 detectors is ~6.1 sr, the ''intrinsic'' efficiency is <math>~9% * 4 * \pi / 6.1 = 18.5%</math>.

-	~~Note that the "Measured Counts" are not sent to Gosia as they appear~~ in ~~the table.~~	+	===Estimating count rates in a planned experiment===

-	Before the ~~measured~~ counts are sent to Gosia they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.	+	Since the "Corrected counts" column gives the absolute emitted particle-gamma counts into the particle- and gamma-detector solid angles (for 100% efficiency), these values can be combined with better known values of the particle and gamma efficiencies to predict more accurate "raw" count rates.
		+
		+
		+	===Simulated data passed to Gosia===
		+
		+	Note that the "Corrected Counts" are not sent to Gosia as they appear in the table.
		+
		+	Before the corrected counts are sent to Gosia they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.

Hayes: added link to page on normalizations

2011-06-14T18:20:17Z

added link to page on normalizations

← Older revision		Revision as of 18:20, 14 June 2011
Line 21:		Line 21:
	Note that the "Measured Counts" are not sent to Gosia as they appear in the table.		Note that the "Measured Counts" are not sent to Gosia as they appear in the table.

-	Before the measured counts are sent to Gosia they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[]] for ~~details~~.	+	Before the measured counts are sent to Gosia they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[Normalization_of_yield_data \| Normalization of yield data]] for more detail.

Hayes: page creation

2011-06-14T18:18:47Z

page creation

New page

''This page stub will be updated when a new version of Rachel is released.''

The output to the terminal following a "simulated yields" calculation will be similar to the following until the next release version where the formatting will be improved.

<pre>
Experiment 1 Detector 4
Raw Gamma Sim. "Measured"
Band Spin Band Spin Counts Error E(keV) eps Counts Error
--------------------------------------------------------------------------------------------------------
gsb 6.0 gsb 4.0 32055690 5662 325.6 0.251 127940090 6397044
gsb 4.0 gsb 2.0 27961778 5288 213.5 0.204 137318753 6865987
gsb 2.0 gsb 0.0 2352850 1534 93.3 0.065 36044668 1802387
</pre>

Three columns have ambiguous names:

* The "Measured Counts" column gives the true particle-gamma count emitted into the solid angles of the particle and Ge detectors. This is equivalent to the total counts for the beam run if the detectors were blackbodies (had 100% efficiency), or if the raw counts were corrected for the absolute p-gamma efficiency.
* The "Raw Counts" column is only used in this output and is not passed to Gosia. The raw counts are the estimated p-gamma counts in the raw data, assuming an approximate absolute efficiency for Gammasphere.
* The "Sim. eps" column gives the assumed absolute efficiency <math>(\epsilon_{abs})</math> of the Ge detector used to calculate the value in the "Raw counts" column.

Note that the "Measured Counts" are not sent to Gosia as they appear in the table.

Before the measured counts are sent to Gosia they are divided by the solid angle of each Ge detector to match Gosia's calculation of yield per unit solid angle. (This is true in Rachel versions >= 2.0.7.beta; in earlier versions this solid angle factor was neglected.) All yields sent to Gosia are then multiplied by <math>4\pi</math> so that weak transitions have counts >=1. This overall multiplicative constant of <math>4\pi</math> has no effect on fitting. Refer to the page [[]] for details.