.10.626331.10.417631.10.291841.10.211341.10.157521.10.1<strong>2007</strong>1.10.098891.10.087491.10.077611.10.069081.10.061631.10.055101.10.049361.10.044271.10.039741.10.035731.10.032131These three records, NDST, UPL <strong>and</strong> YNRM, determine thefollowing for experiment #3. There is one <strong>data</strong> set (usually,but not always equal to the number of gamma detectors enteredabove). Any unobserved yield calculated to be greater than0.1×Y 5→4 will enter the chi-squared statistic. Thenormalization constant with respect to experiment #1 is0.41763. These records continue for each experiment with apositive value of NANG.174
.10.028921.10.026031.10.023431.10.021081.10.018961.10.017051.10.015331.10.013771.10.012373 !NTAP0,01,1.03,112.,1.0,00,0OP,INTG5,7,634,650,22.000,26.000634,638,642,646,65021.,22.,23.,24.,25.,26.,27.5,4,634,650,26.000,30.000634,638,642,646,65025.,27.,29.,315,4,634,650,30.000,34.000634,638,642,646,65029.,31.,33.,355,4,634,650,34.000,38.000634,638,642,646,65033.,35.,37.,395,4,634,650,38.000,42.000NTAP=3. Experimental gamma-rayyields <strong>data</strong> will be read from tape 3defined as the file g1demo.yld in OP,FILEabove. This will be changed to 0 (no <strong>data</strong>read) or 4 for different operationsdescribed in the tutorial in this chapter.(Note that comments such as “!NTAP”are helpful, but not legal on every line.When in doubt, do not use a comment.)634,638,642,646,65037.,39.,41.,435,4,634,650,42.000,46.000634,638,642,646,65041.,43.,45.,47Number (NBRA=0) <strong>and</strong> weight (WBRA=0) ofbranching ratio <strong>data</strong> to be included in chi-squared.One experimental lifetime (NL=1) will be included inthe chi-squared statistic with weight WL=1.0.The experimental lifetime of state INDEX=3 is 112(1)ps.No experimental values of mixing ratios or matrixelements will be entered.Integration limits of energy <strong>and</strong>scattering angle for experiment #1.Because azimuthal symmetry <strong>and</strong>360-degree coverage of the particledetector were specified in EXPT(above), only the polar angleinformation is entered. In the presentexample, each experiment requiresonly three lines of input. Otherchoices of particle-detector geometrymay require more.175
- Page 1:
COULOMB EXCITATION DATA ANALYSIS CO
- Page 4 and 5:
10 MINIMIZATION BY SIMULATED ANNEAL
- Page 6 and 7:
1 INTRODUCTION1.1 Gosia suite of Co
- Page 8 and 9:
104 Ru, 110 Pd, 165 Ho, 166 Er, 186
- Page 13 and 14:
Figure 1: Coordinate system used to
- Page 15 and 16:
Cλ E =1.116547 · (13.889122) λ (
- Page 17 and 18:
Figure 2: The orbital integrals R 2
- Page 19 and 20:
2.2 Gamma Decay Following Electroma
- Page 21 and 22:
where :d 2 σ= σ R (θ p ) X R kχ
- Page 23 and 24:
Formula 2.49 is valid only for t mu
- Page 25 and 26:
à XK(α) =exp−iτ i (E γ )x i (
- Page 27 and 28:
important to have an accurate knowl
- Page 29 and 30:
3 APPROXIMATE EVALUATION OF EXCITAT
- Page 31 and 32:
with the reduced matrix element M c
- Page 33 and 34:
q (20)s (0 + → 2 + ) · M 1 ζ (2
- Page 35 and 36:
esults of minimization and error ru
- Page 37 and 38:
adjustment of the stepsize accordin
- Page 39 and 40:
approximation reliability improves
- Page 41 and 42:
Zd 2 σ(I → I f )Y (I → I f )=s
- Page 43 and 44:
4.5 MinimizationThe minimization, i
- Page 45 and 46:
X(CC k Yk c − Yk e ) 2 /σ 2 k =m
- Page 47 and 48:
However, estimation of the stepsize
- Page 49 and 50:
It can be shown that as long as the
- Page 51 and 52:
een exceeded; third, the user-given
- Page 53 and 54:
where f k stands for the functional
- Page 55 and 56:
x i + δx i Rx iexp ¡ − 1 2 χ2
- Page 57 and 58:
method used for the minimization, i
- Page 59 and 60:
OP,ERRO (ERRORS) (5.6):Activates th
- Page 61 and 62:
-----OP,SIXJ (SIX-j SYMBOL) (5.25):
- Page 63 and 64:
5.3 CONT (CONTROL)This suboption of
- Page 65 and 66:
I,I1 Ranges of matrix elements to b
- Page 67 and 68:
CODE DEFAULT OTHER CONSEQUENCES OF
- Page 69 and 70:
5.4 OP,CORR (CORRECT )This executio
- Page 71 and 72:
5.6 OP,ERRO (ERRORS)ThemoduleofGOSI
- Page 73 and 74:
5.7 OP,EXIT (EXIT)This option cause
- Page 75 and 76:
M AControls the number of magnetic
- Page 77 and 78:
5.10 OP,GDET (GE DETECTORS)This opt
- Page 79 and 80:
5.12 OP,INTG (INTEGRATE)This comman
- Page 81 and 82:
¡ dE¢dx1 ..¡ dEdx¢Stopping powe
- Page 83 and 84:
NI1, NI2 Number of subdivisions of
- Page 85 and 86:
5.13 LEVE (LEVELS)Mandatory subopti
- Page 87 and 88:
5.15 ME (OP,COUL)Mandatory suboptio
- Page 89 and 90:
Figure 10: Model system having 4 st
- Page 91 and 92:
ME =< INDEX2||E(M)λ||INDEX1 > The
- Page 93 and 94:
When entering matrix elements in th
- Page 95 and 96:
There are no restrictions concernin
- Page 97 and 98:
5.18 OP,POIN (POINT CALCULATION)Thi
- Page 99 and 100:
5.20 OP,RAW (RAW UNCORRECTED γ YIE
- Page 101 and 102:
5.21 OP,RE,A (RELEASE,A)This option
- Page 103 and 104:
5.25 OP,SIXJ (SIXJ SYMBOL)This stan
- Page 105 and 106:
5.27 OP,THEO (COLLECTIVE MODEL ME)C
- Page 107 and 108:
2,5,1,-2,23,5,1,-2,23,6,1,-2,2Matri
- Page 109 and 110:
5.29 OP,TROU (TROUBLE)This troubles
- Page 111 and 112:
to that of the previous experiment,
- Page 113 and 114:
To reduce the unnecessary input, on
- Page 115 and 116:
OP,STAR or OP,POIN under OP,GOSI. N
- Page 117 and 118:
5.31 INPUT OF EXPERIMENTAL γ-RAY Y
- Page 119 and 120:
6 QUADRUPOLE ROTATION INVARIANTS -
- Page 121 and 122:
*½P 5 (J) = s(E2 × E2) J ׯh¾
- Page 123 and 124: The expectation value of cos3δ can
- Page 125 and 126: where ē is an arbitratry vector. D
- Page 127 and 128: achieved using “mixed“ calculat
- Page 129 and 130: TAPE9 Contains the parameters neede
- Page 131 and 132: TAPE18 Input file, containing the i
- Page 133 and 134: 7.4.4 CALCULATION OF THE INTEGRATED
- Page 135 and 136: OP,EXITInput: TAPE4,TAPE7,TAPE9Outp
- Page 137 and 138: OP,ERRO0,MS,MEND,1,0,RMAXand the fi
- Page 139 and 140: 8 SIMULTANEOUS COULOMB EXCITATION:
- Page 141 and 142: 4, 3, 1kr88.corKr corrected yields
- Page 143 and 144: 0 Correction for in-flight decay ch
- Page 145 and 146: OP, ERRO Estimation of errors of fi
- Page 147 and 148: 9 COULOMB EXCITATION OF ISOMERIC ST
- Page 149 and 150: configurations with a probability e
- Page 151 and 152: The average range covered by each m
- Page 153 and 154: SFX,NTOTI1(1),I2(1),RSIGN(1)I1(2),I
- Page 155 and 156: 11.2 LearningtoWriteGosiaInputsThe
- Page 157 and 158: (1.6 MeV)1.1 MeV0.75 MeV0.4 MeV0.08
- Page 159 and 160: Define the germaniumdetector geomet
- Page 161 and 162: Figure 15: Flow diagram for Gosia m
- Page 163 and 164: gosia < 2-make-correction-factors.i
- Page 165 and 166: Issue the commandgosia < 9-diag-err
- Page 167 and 168: At this point, it is suggested to c
- Page 169 and 170: calculation.) In this case, a copy
- Page 171 and 172: 4,-4, -3.705, 3,44,5, 4.626, 3.,7.5
- Page 173: 90145901459014590145901459014590145
- Page 177 and 178: 5,5,634,650,82.000,84.000634,638,64
- Page 179 and 180: ***********************************
- Page 181 and 182: *** CHISQ= 0.134003E+01 ***MATRIX E
- Page 183 and 184: CALCULATED AND EXPERIMENTAL YIELDS
- Page 185 and 186: 11.7 Annotated excerpt from a Coulo
- Page 187 and 188: 11.8 Accuracy and speed of calculat
- Page 189 and 190: 18,10.056,0.068,0.082,0.1,0.12,0.15
- Page 191 and 192: line 152 Eu 182 Tanumber (keV) (keV
- Page 193 and 194: 1.6 Normalization between data sets
- Page 195 and 196: 13 GOSIA 2007 RELEASE NOTESThese no
- Page 197 and 198: Matrix elements 500(April 1990, T.
- Page 199 and 200: 14 GOSIA Manual UpdatesDATE UPDATE2
- Page 201 and 202: [KIB08]T.Kibédi,T.W.Burrows,M.B.Tr