Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/util/src/G4NuclearRadii.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/util/src/G4NuclearRadii.cc (Version 11.3.0) and /processes/hadronic/util/src/G4NuclearRadii.cc (Version 11.0.p2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27 // Geant4 class G4NuclearRadii                     27 // Geant4 class G4NuclearRadii
 28 //                                                 28 //
 29 // Author V.Ivanchenko 27.05.2019                  29 // Author V.Ivanchenko 27.05.2019
 30 //                                                 30 //
 31 //                                                 31 //
 32                                                    32 
 33 #include "G4NuclearRadii.hh"                       33 #include "G4NuclearRadii.hh"
 34 #include "G4Pow.hh"                                34 #include "G4Pow.hh"
 35 #include "G4PhysicalConstants.hh"                  35 #include "G4PhysicalConstants.hh"
 36 #include "G4ParticleDefinition.hh"                 36 #include "G4ParticleDefinition.hh"
 37 #include "G4NucleiProperties.hh"                   37 #include "G4NucleiProperties.hh"
 38 #include "G4IsotopeList.hh"                    << 
 39 #include "G4Log.hh"                            << 
 40 #include "G4Exp.hh"                            << 
 41                                                    38 
 42 G4Pow* G4NuclearRadii::fG4pow = G4Pow::GetInst     39 G4Pow* G4NuclearRadii::fG4pow = G4Pow::GetInstance();
 43                                                <<  40 const G4double fAlpha = 0.5*CLHEP::fine_structure_const*CLHEP::hbarc;
 44 namespace                                      <<  41 const G4double fInvep = 1.0/CLHEP::eplus;
 45 {                                              << 
 46   const G4double llog10 = G4Log(10.);          << 
 47   const G4double fAlpha = 0.5*CLHEP::fine_stru << 
 48   const G4double fInvep = 1.0/CLHEP::eplus;    << 
 49 }                                              << 
 50                                                    42 
 51 G4double G4NuclearRadii::ExplicitRadius(G4int      43 G4double G4NuclearRadii::ExplicitRadius(G4int Z, G4int A)
 52 {                                                  44 {
 53   G4double R = 0.0;                                45   G4double R = 0.0;
 54   // Special rms radii for light nucleii           46   // Special rms radii for light nucleii
 55   if(Z <= 4) {                                     47   if(Z <= 4) {
 56     if(A == 1)                { R = 0.895*CLHE     48     if(A == 1)                { R = 0.895*CLHEP::fermi; }// p
 57     else if(A == 2)           { R = 2.13*CLHEP     49     else if(A == 2)           { R = 2.13*CLHEP::fermi; }// d
 58     else if(Z == 1 && A == 3) { R = 1.80*CLHEP     50     else if(Z == 1 && A == 3) { R = 1.80*CLHEP::fermi; }// t
 59     else if(Z == 2 && A == 3) { R = 1.96*CLHEP     51     else if(Z == 2 && A == 3) { R = 1.96*CLHEP::fermi; }// He3
 60     else if(Z == 2 && A == 4) { R = 1.68*CLHEP     52     else if(Z == 2 && A == 4) { R = 1.68*CLHEP::fermi; }// He4
 61     else if(Z == 3)           { R = 2.40*CLHEP     53     else if(Z == 3)           { R = 2.40*CLHEP::fermi; }// Li7
 62     else if(Z == 4)           { R = 2.51*CLHEP     54     else if(Z == 4)           { R = 2.51*CLHEP::fermi; }// Be9
 63   }                                                55   }
 64   return R;                                        56   return R;
 65 }                                                  57 }
 66                                                    58 
 67 G4double G4NuclearRadii::Radius(G4int Z, G4int     59 G4double G4NuclearRadii::Radius(G4int Z, G4int A)
 68 {                                                  60 {
 69   G4double R = ExplicitRadius(Z, A);               61   G4double R = ExplicitRadius(Z, A);
 70   if(0.0 == R) {                                   62   if(0.0 == R) {
 71     if (A <= 50) {                                 63     if (A <= 50) {
 72       G4double y = 1.1;                            64       G4double y = 1.1;
 73       if( A <= 15)      { y = 1.26; }              65       if( A <= 15)      { y = 1.26; }
 74       else if( A <= 20) { y = 1.19; }              66       else if( A <= 20) { y = 1.19; }
 75       else if( A <= 30) { y = 1.12; }              67       else if( A <= 30) { y = 1.12; }
 76       G4double x = fG4pow->Z13(A);                 68       G4double x = fG4pow->Z13(A);
 77       R = y*(x - 1./x);                            69       R = y*(x - 1./x);
 78     } else {                                       70     } else {
 79       R = fG4pow->powZ(A, 0.27);                   71       R = fG4pow->powZ(A, 0.27);
 80     }                                              72     }
 81     R *= CLHEP::fermi;                             73     R *= CLHEP::fermi;
 82   }                                                74   }
 83   return R;                                        75   return R;
 84 }                                                  76 }
 85                                                    77 
 86 G4double G4NuclearRadii::RadiusRMS(G4int Z, G4     78 G4double G4NuclearRadii::RadiusRMS(G4int Z, G4int A)
 87 {                                                  79 {
 88   G4double R = ExplicitRadius(Z, A);               80   G4double R = ExplicitRadius(Z, A);
 89   if(0.0 == R) {                                   81   if(0.0 == R) {
 90     R = 1.24*fG4pow->powZ(A, 0.28)*CLHEP::ferm     82     R = 1.24*fG4pow->powZ(A, 0.28)*CLHEP::fermi;
 91   }                                                83   }
 92   return R;                                        84   return R;
 93 }                                                  85 }
 94                                                    86 
 95 G4double G4NuclearRadii::RadiusNNGG(G4int Z, G     87 G4double G4NuclearRadii::RadiusNNGG(G4int Z, G4int A)
 96 {                                                  88 {
 97   G4double R = ExplicitRadius(Z, A);               89   G4double R = ExplicitRadius(Z, A);
 98   if(0.0 == R) {                                   90   if(0.0 == R) {
 99     if(A > 20) {                                   91     if(A > 20) {
100       R = 1.08*fG4pow->Z13(A)*(0.85 + 0.15*G4E     92       R = 1.08*fG4pow->Z13(A)*(0.85 + 0.15*G4Exp(-(G4double)(A - 21)/40.));
101     } else {                                       93     } else {
102       R = 1.08*fG4pow->Z13(A)*(1.0 + 0.3*G4Exp     94       R = 1.08*fG4pow->Z13(A)*(1.0 + 0.3*G4Exp(-(G4double)(A - 21)/10.));
103     }                                              95     }
104     R *= CLHEP::fermi;                             96     R *= CLHEP::fermi;
105   }                                                97   }
106   return R;                                        98   return R;
107 }                                                  99 }
108                                                   100 
109 G4double G4NuclearRadii::RadiusECS(G4int Z, G4    101 G4double G4NuclearRadii::RadiusECS(G4int Z, G4int A)
110 {                                                 102 {
111   G4double R=0.;                                  103   G4double R=0.;
112   const G4double c[3]={0.77329745, 1.38206072,    104   const G4double c[3]={0.77329745, 1.38206072, 30.28295235};
113   const G4double c1=c[0];                         105   const G4double c1=c[0];
114   const G4double c2=c[1];                         106   const G4double c2=c[1];
115   const G4double c3=c[2];                         107   const G4double c3=c[2];
116                                                   108 
117   // Special rms radii for light nuclei           109   // Special rms radii for light nuclei
118   if (A <= 30) {                                  110   if (A <= 30) {
119     G4double vn = 0.5*A + fG4pow->powN(0.028*A    111     G4double vn = 0.5*A + fG4pow->powN(0.028*A,2) - fG4pow->powN(0.011*A,3);
120     G4double dev = vn - (A-Z);                    112     G4double dev = vn - (A-Z);
121     R = c1*fG4pow->Z13(A) + c2/fG4pow->Z13(A)     113     R = c1*fG4pow->Z13(A) + c2/fG4pow->Z13(A) + c3*dev*dev/(A*A);
122   } else if (A<=50){                              114   } else if (A<=50){
123     G4double y = 1.1;                             115     G4double y = 1.1; 
124     G4double x = fG4pow->Z13(A);                  116     G4double x = fG4pow->Z13(A);
125     R = y*(x - 1./x);                             117     R = y*(x - 1./x);
126   }                                               118   }
127   return R*CLHEP::fermi;                          119   return R*CLHEP::fermi;
128 }                                                 120 }
129                                                   121 
130 G4double G4NuclearRadii::RadiusHNGG(G4int A)      122 G4double G4NuclearRadii::RadiusHNGG(G4int A)
131 {                                                 123 {
132   G4double R = CLHEP::fermi;                      124   G4double R = CLHEP::fermi;
133   if(A > 20) {                                    125   if(A > 20) {
134     R *= 1.08*fG4pow->Z13(A)*(0.8 + 0.2*G4Exp(    126     R *= 1.08*fG4pow->Z13(A)*(0.8 + 0.2*G4Exp(-(G4double)(A - 20)/20.));
135   } else {                                        127   } else {
136     R *= 1.08*fG4pow->Z13(A)*(1.0 + 0.1*G4Exp(    128     R *= 1.08*fG4pow->Z13(A)*(1.0 + 0.1*G4Exp(-(G4double)(A - 20)/20.));
137   }                                               129   }
138   return R;                                       130   return R;
139 }                                                 131 }
140                                                   132 
141 G4double G4NuclearRadii::RadiusKNGG(G4int A)      133 G4double G4NuclearRadii::RadiusKNGG(G4int A)
142 {                                                 134 {
143   return 1.3*CLHEP::fermi*fG4pow->Z13(A);         135   return 1.3*CLHEP::fermi*fG4pow->Z13(A);
144 }                                                 136 }
145                                                   137 
146 G4double G4NuclearRadii::RadiusND(G4int A)        138 G4double G4NuclearRadii::RadiusND(G4int A)
147 {                                                 139 {
148   G4double R = CLHEP::fermi;                      140   G4double R = CLHEP::fermi;
149   if (1 == A) {                                << 141   if(1 == A) { return R*0.895; }
150     R *= 0.895;                                << 142 //  G4double x = R*fG4pow->Z13(A);
151   } else {                                     << 143 //  if(A <= 3.) { x *= 0.8; }
152     R *= fG4pow->Z13(A);                       << 144 //  else { x *= 1.7; }
153     if (A <= 3.) { R *= 0.8; }                 << 
154     else { R *= 1.7; }                         << 
155   }                                            << 
156   return R;                                       145   return R;
157 }                                                 146 }
158                                                   147 
159 G4double G4NuclearRadii::RadiusCB(G4int Z, G4i    148 G4double G4NuclearRadii::RadiusCB(G4int Z, G4int A)
160 {                                                 149 {
161   G4double R = ExplicitRadius(Z, A);              150   G4double R = ExplicitRadius(Z, A);
162   if(0.0 == R) {                                  151   if(0.0 == R) {
163     G4int z = std::min(Z, 92);                    152     G4int z = std::min(Z, 92);
164     R = r0[z]*fG4pow->Z13(A)*CLHEP::fermi;        153     R = r0[z]*fG4pow->Z13(A)*CLHEP::fermi;
165   }                                               154   }
166   return R;                                       155   return R;
167 }                                                 156 }
168                                                   157 
169 G4double G4NuclearRadii::ParticleRadius(const     158 G4double G4NuclearRadii::ParticleRadius(const G4ParticleDefinition* p)
170 {                                                 159 {
171   G4double R = CLHEP::fermi;                      160   G4double R = CLHEP::fermi;
172   G4int pdg = std::abs(p->GetPDGEncoding());      161   G4int pdg = std::abs(p->GetPDGEncoding());
173   if(pdg == 2112 || pdg == 2212)   { R *= 0.89    162   if(pdg == 2112 || pdg == 2212)   { R *= 0.895; }
174   else if(pdg == 211)  { R *= 0.663; }            163   else if(pdg == 211)  { R *= 0.663; }
175   else if(pdg == 321)  { R *= 0.340; }            164   else if(pdg == 321)  { R *= 0.340; }
176   else { R *= 0.5; }                              165   else { R *= 0.5; }
177   return R;                                       166   return R;
178 }                                                 167 }
179                                                   168 
180 G4double G4NuclearRadii::CoulombFactor(           169 G4double G4NuclearRadii::CoulombFactor(
181          const G4ParticleDefinition* thePartic    170          const G4ParticleDefinition* theParticle, 
182    const G4ParticleDefinition* nucleon,           171    const G4ParticleDefinition* nucleon, 
183    G4double ekin)                                 172    G4double ekin)
184 {                                                 173 {
185   G4double tR = 0.895*CLHEP::fermi;               174   G4double tR = 0.895*CLHEP::fermi;
186   G4double pR = ParticleRadius(theParticle);      175   G4double pR = ParticleRadius(theParticle);
187                                                   176 
188   G4double pZ = theParticle->GetPDGCharge()*fI    177   G4double pZ = theParticle->GetPDGCharge()*fInvep;
189   G4double tZ = nucleon->GetPDGCharge()*fInvep    178   G4double tZ = nucleon->GetPDGCharge()*fInvep;
190                                                   179 
191   G4double pM = theParticle->GetPDGMass();        180   G4double pM = theParticle->GetPDGMass(); 
192   G4double tM = nucleon->GetPDGMass();            181   G4double tM = nucleon->GetPDGMass();
193                                                   182 
194   G4double pElab = ekin + pM;                     183   G4double pElab = ekin + pM;
195   G4double totTcm  = std::sqrt(pM*pM + tM*tM +    184   G4double totTcm  = std::sqrt(pM*pM + tM*tM + 2.*pElab*tM) - pM -tM;
196                                                   185   
197   G4double bC = fAlpha*pZ*tZ/(pR + tR);           186   G4double bC = fAlpha*pZ*tZ/(pR + tR);
198   return (totTcm > bC) ? 1. - bC/totTcm : 0.0;    187   return (totTcm > bC) ? 1. - bC/totTcm : 0.0;
199 }                                                 188 }
200                                                   189 
201 G4double G4NuclearRadii::CoulombFactor(           190 G4double G4NuclearRadii::CoulombFactor(
202          G4int Z, G4int A,                        191          G4int Z, G4int A,
203    const G4ParticleDefinition* theParticle,       192    const G4ParticleDefinition* theParticle, 
204    G4double ekin)                                 193    G4double ekin)
205 {                                                 194 {
206   G4double tR = RadiusCB(Z, A);                   195   G4double tR = RadiusCB(Z, A);
207   G4double pR = ParticleRadius(theParticle);      196   G4double pR = ParticleRadius(theParticle);
208                                                   197 
209   G4double pZ = theParticle->GetPDGCharge()*fI    198   G4double pZ = theParticle->GetPDGCharge()*fInvep;
210                                                   199 
211   G4double pM = theParticle->GetPDGMass();        200   G4double pM = theParticle->GetPDGMass(); 
212   G4double tM = G4NucleiProperties::GetNuclear    201   G4double tM = G4NucleiProperties::GetNuclearMass(A, Z);
213                                                   202 
214   G4double pElab = ekin + pM;                     203   G4double pElab = ekin + pM;
215   G4double totTcm  = std::sqrt(pM*pM + tM*tM +    204   G4double totTcm  = std::sqrt(pM*pM + tM*tM + 2.*pElab*tM) - pM -tM;
216                                                   205   
217   G4double bC = fAlpha*pZ*Z/(pR + tR);            206   G4double bC = fAlpha*pZ*Z/(pR + tR);
218   return (totTcm > bC) ? 1. - bC/totTcm : 0.0;    207   return (totTcm > bC) ? 1. - bC/totTcm : 0.0;
219 }                                              << 
220                                                << 
221 G4double                                       << 
222 G4NuclearRadii::NeutronInelasticShape(G4int Z, << 
223 {                                              << 
224   G4double A = (Z < 100) ? aeff[Z] : aeff[100] << 
225   G4double elog = G4Log(ekin/CLHEP::GeV)/llog1 << 
226   G4double p3 = 0.6 + 13./A - 0.0005*A;        << 
227   G4double p4 = 7.2449 - 0.018242*A;           << 
228   G4double p5 = 1.36 + 1.8/A + 0.0005*A;       << 
229   G4double p6 = 1. + 200./A + 0.02*A;          << 
230   G4double p7 = 3.0 - (A - 70.)*(A - 200.)/110 << 
231                                                << 
232   G4double firstexp = G4Exp(-p4*(elog + p5));  << 
233   G4double secondexp = G4Exp(-p6*(elog + p7)); << 
234                                                << 
235   return (1. + p3*firstexp/(1. + firstexp))/(1 << 
236 }                                              << 
237                                                << 
238 G4double                                       << 
239 G4NuclearRadii::ProtonInelasticShape(G4int Z,  << 
240 {                                              << 
241   G4double A = (Z < 100) ? aeff[Z] : aeff[100] << 
242   G4double elog = G4Log(ekin/CLHEP::GeV)/llog1 << 
243   G4double ff1 = 5.6  - 0.016*A; // slope of t << 
244   G4double ff2 = 1.37 + 1.37/A;  // start of t << 
245   G4double ff3 = 0.8  + 18./A - 0.002*A;   //  << 
246   G4double res = (1.0 + ff3*(1.0 - (1.0/(1+G4E << 
247   ff1 = 8. - 8./A - 0.008*A; // slope of the r << 
248   ff2 = 2.34 - 5.4/A - 0.0028*A; // start of t << 
249   res /= (1.0 + G4Exp(-ff1*(elog + ff2)));     << 
250   return res;                                  << 
251 }                                                 208 }
252                                                   209 
253 const G4double G4NuclearRadii::r0[] = {           210 const G4double G4NuclearRadii::r0[] = {
254  1.2,                                             211  1.2,
255  1.3, 1.3, 1.3, 1.3,1.17,1.54,1.65,1.71, 1.7,1    212  1.3, 1.3, 1.3, 1.3,1.17,1.54,1.65,1.71, 1.7,1.75, // 1-10
256  1.7,1.57,1.53, 1.4, 1.3,1.30,1.44, 1.4, 1.4,     213  1.7,1.57,1.53, 1.4, 1.3,1.30,1.44, 1.4, 1.4, 1.4, //11-20
257  1.4, 1.4,1.46, 1.4, 1.4,1.46,1.55, 1.5,1.38,1    214  1.4, 1.4,1.46, 1.4, 1.4,1.46,1.55, 1.5,1.38,1.48, //21-30
258  1.4, 1.4, 1.4,1.46, 1.4, 1.4, 1.4, 1.4, 1.4,1    215  1.4, 1.4, 1.4,1.46, 1.4, 1.4, 1.4, 1.4, 1.4,1.45, //31-40
259  1.4, 1.4, 1.4, 1.4, 1.4, 1.4,1.45,1.48, 1.4,1    216  1.4, 1.4, 1.4, 1.4, 1.4, 1.4,1.45,1.48, 1.4,1.52, //41-50
260 1.46, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4,     217 1.46, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.4, 1.5, //51-60
261  1.4, 1.4, 1.4, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3,     218  1.4, 1.4, 1.4, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.4, //61-70
262  1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3,1.33,1    219  1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3,1.33,1.43, //71-80
263  1.3,1.32,1.34, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3,     220  1.3,1.32,1.34, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, 1.3, //81-90
264  1.3, 1.3};                                       221  1.3, 1.3};
265                                                   222