3, 5, 190, 194, 17.8, 28 Number of energy <strong>and</strong> theta meshpoints, energy range, theta range of detectedparticle, in this case target.190, 192, 194 Energy meshpoints−3, −3.5, −4., −4.5, −5. Projectile scattering angle meshpoints (appended minus sign when target wasactually detected)3, 5, 190, 194, 28, 32.8 Number of energy <strong>and</strong> theta meshpoints, etc. for experiemnt 2..190, 192, 194−5., −5.7, −6.4, −7.1, −7.71 As above, for experiment 2.3 Number of meshpoints to define dE/dx190, 192, 194 Meshpoints10, 10, 10 dE/dx values20, 20 Number of subdivisions for numerical integration0 Use the same dE/dx table for experiment 220, 20 Number of subdivisions for experiment 2.In the present case input to OP,INTG is very simple — we defined axially symmetric ring detectors, so itis not necessary to define theta-phi dependence to reproduce arbitrary shape. Also the geometry is identicalfor both experiments. This corresponds to usual MINIBALL situation. OP,INTG offers the full descriptionof an experiment, not only-as OP,POIN does-limited to mean values of bombarding energy <strong>and</strong> scatteringangle. Input to OP,INTG is identical to that of GOSIA.OP, CORR Converts the experimental yields (unit 3) to the values corresponding to „point yields”(unit 4). Corrected yields thus replace „original” yields <strong>and</strong> are subsequently used in the fittingprocedure. Should immediately follow OP,INTGOP, RESTcode.Overwrites the set of matrix elements defined by ME. Should be used when re-running the1,.5 Matrix element #1 set to .5 by h<strong>and</strong>. Of course remove this when running fitting again, otherwiseonce again the value will be reset.0, 0 Ends the input.A sequenceOP, REST0, 0means read the matrix elements as stored with no h<strong>and</strong>-made changesOP, MINI2100, 5,.0001,.0001, 1.01, 0, 0, 1, 0, 0 Fitting parameters as in GOSIA, but the way GOSIA2 treatsfitting is different. GOSIA2, along with fitting the matrix elements fits normalization constants insuch a way that the constants are the same for corresponding <strong>data</strong>sets from both projectile <strong>and</strong> target<strong>excitation</strong>. The parameters given here are recommended ones. Before OP,MINI is executed OP,MAPmust be executed (only once for each case (projectile or target <strong>excitation</strong>).144
OP, ERRO Estimation of errors of fitted matrix elements <strong>and</strong> optionally creation of SIGMA input file.Executed separately for both target <strong>and</strong> projectile exactly like in GOSIA.Parallel input for target — very similar to that of projectile. Below there is a list of optionsjust pointing out what are the differences (easy to guess):2 –––––––-input must be labelled 2OP,FILE –––– change kr88 to c12 where appropriate22,3,1c12.out25,3,1kr88.inp26,3,1c12.inp14,3,1dum.1417,3,1dum.1715,3,1dum.1512,3,1kr88.me32,3,1c12.me9,3,1det.inp11,3,2crf.dat3,3,1c12.org27,3,1c12.map4,3,1c12.cor7,3,1kr88.map18,3,1dum.1833,3,1c12.smr13,3,1cnor.dat23,3,1c12.raw0,0,0OP,TITLTarget excited –––––––––––—or any other titleOP,GOSILEVE –––––––––––––––—level scheme of c12 instead of kr881,1,0,02,1,2,.83,1,4,2.4,1,2,1.0,0,0,0ME ––––––––––––––––––––––-c12 matrix elements145
- 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: 0 Correction for in-flight decay ch
- 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 and 174: 90145901459014590145901459014590145
- Page 175 and 176: .10.028921.10.026031.10.023431.10.0
- 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