source: trunk/source/processes/electromagnetic/lowenergy/src/G4QAOLowEnergyLoss.cc @ 1337

Last change on this file since 1337 was 819, checked in by garnier, 16 years ago

import all except CVS
File size: 24.4 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer                                           *
4// *                                                                  *
5// * The  Geant4 software  is  copyright of the Copyright Holders  of *
6// * the Geant4 Collaboration.  It is provided  under  the terms  and *
7// * conditions of the Geant4 Software License,  included in the file *
8// * LICENSE and available at  http://cern.ch/geant4/license .  These *
9// * include a list of copyright holders.                             *
10// *                                                                  *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work  make  any representation or  warranty, express or implied, *
14// * regarding  this  software system or assume any liability for its *
15// * use.  Please see the license in the file  LICENSE  and URL above *
16// * for the full disclaimer and the limitation of liability.         *
17// *                                                                  *
18// * This  code  implementation is the result of  the  scientific and *
19// * technical work of the GEANT4 collaboration.                      *
20// * By using,  copying,  modifying or  distributing the software (or *
21// * any work based  on the software)  you  agree  to acknowledge its *
22// * use  in  resulting  scientific  publications,  and indicate your *
23// * acceptance of all terms of the Geant4 Software license.          *
24// ********************************************************************
25//
26//
27// -------------------------------------------------------------
28//      GEANT 4 class implementation file
29//
30//      History: New Implementation
31//   
32//      ---------- G4QAOLowEnergyLoss physics process -------
33//                  by Stephane Chauvie, 5 May 2000
34// Modified:
35//
36// 24/05/2000 MGP  Modified to remove compilation warnings on Linux and DEC
37//                 Introduced sizes of L0, L1, L2 arrays
38// 23/05/2000 MGP  Made compliant to design
39// 02/08/2000 V.Ivanchenko Clean up according new design
40// 16/09/2000 S. Chauvie  Oscillator for all materials
41// 03/10/2000 V.Ivanchenko CodeWizard clean up
42// 05/11/2000 V.Ivanchenko "Aluminum" - correct name, end of cycle
43//            over shells, and two bugs from previous edition
44// 10/05/2001 V.Ivanchenko Clean up againist Linux compilation with -Wall
45// 13/05/2001 S. Chauvie corrected bugs
46// 01/06/2001 V.Ivanchenko replace names by Z, change the validity range
47//                         from 50 keV to 5 KeV and change sign of the
48//                         Barkas term
49// 4/06/2001 S. Chauvie  Corrected small bugs
50//
51// ************************************************************
52// It is the Quantal Harmonic Oscillator Model for energy loss
53// of slow antiproton
54// ************************************************************
55// --------------------------------------------------------------
56
57//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
58
59#include "G4QAOLowEnergyLoss.hh"
60#include "G4DynamicParticle.hh"
61#include "G4Material.hh"
62#include "G4ParticleDefinition.hh"
63#include "G4AntiProton.hh"
64
65//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
66
67G4QAOLowEnergyLoss::G4QAOLowEnergyLoss(const G4String& name)
68  : G4VLowEnergyModel(name)
69{
70  numberOfMaterials = 6;
71  sizeL0 = 67;
72  sizeL1 = 22;
73  sizeL2 = 14;
74}
75
76
77G4QAOLowEnergyLoss::~G4QAOLowEnergyLoss()
78{;}
79
80
81G4double G4QAOLowEnergyLoss::HighEnergyLimit(const G4ParticleDefinition* ,
82                                             const G4Material* ) const
83{
84  return 2.0*MeV ;
85}
86
87
88G4double G4QAOLowEnergyLoss::LowEnergyLimit(const G4ParticleDefinition* ,
89                                            const G4Material* ) const
90{
91  //  return 50.0*keV ;
92  return 5.0*keV ;
93}
94
95
96G4double G4QAOLowEnergyLoss::HighEnergyLimit(const G4ParticleDefinition* ) const
97{
98  return 2.0*MeV ;
99}
100
101G4double G4QAOLowEnergyLoss::LowEnergyLimit(const G4ParticleDefinition* ) const
102{
103  //  return 50.0*keV ;
104  return 5.0*keV ;
105}
106
107
108G4bool G4QAOLowEnergyLoss::IsInCharge(const G4DynamicParticle* particle,
109                                      const G4Material* material) const
110{
111  G4bool isInCharge = false;
112
113  G4bool hasMaterial = false;
114
115  if (material->GetNumberOfElements() == 1) hasMaterial = true;
116 
117  if ((particle->GetDefinition()) == (G4AntiProton::AntiProtonDefinition())
118               && hasMaterial) isInCharge = true;
119 
120  return isInCharge;
121
122}
123
124G4bool G4QAOLowEnergyLoss::IsInCharge(const G4ParticleDefinition* aParticle,
125                                      const G4Material* material) const
126{
127 
128  G4bool isInCharge = false;
129 
130  G4bool hasMaterial = false;
131 
132  if (material->GetNumberOfElements() == 1) hasMaterial = true;
133 
134
135  if (aParticle == (G4AntiProton::AntiProtonDefinition())
136                && hasMaterial) isInCharge = true;
137 
138  return isInCharge;
139 
140}
141
142
143G4double G4QAOLowEnergyLoss::TheValue(const G4DynamicParticle* particle,
144                                      const G4Material* material) 
145{
146  G4double zParticle = (G4int)(particle->GetCharge())/eplus;
147
148  G4double energy = particle->GetKineticEnergy() ;
149  G4double eloss  = EnergyLoss(material,energy,zParticle) ;
150
151  return eloss ;
152}
153
154
155G4double G4QAOLowEnergyLoss::TheValue(const G4ParticleDefinition* aParticle,
156                                      const G4Material* material,
157                                      G4double kineticEnergy) 
158{
159  G4double zParticle = (aParticle->GetPDGCharge())/eplus;
160
161  G4double eloss  = EnergyLoss(material,kineticEnergy,zParticle) ;
162
163  return eloss ;
164}
165
166
167
168G4double G4QAOLowEnergyLoss::EnergyLoss(const G4Material* material,
169                                        G4double kineticEnergy,
170                                        G4double zParticle) const 
171{
172  G4int nbOfShell = GetNumberOfShell(material);
173  if(nbOfShell < 1) nbOfShell = 1;
174 
175  //G4int iz = (G4int)(material->GetZ());
176  //if(iz < 1) iz = 1;
177  //else if(iz > 100) iz = 100;
178  //G4int nbOfShell = fNumberOfShells[iz];
179 
180  /*
181  if(material->GetName()=="Graphite"){
182  G4cout << " E(MeV)= " << kineticEnergy/MeV << " n= "
183         << nbOfShell
184         << " for " << material->GetZ() << G4endl;
185  }
186  */
187  G4double dedx=0;
188
189  G4double v = c_light * std::sqrt( 2.0 * kineticEnergy / proton_mass_c2 );
190  G4double coeff = twopi * proton_mass_c2 * 
191                  (material-> GetTotNbOfElectPerVolume()) / 
192                   electron_mass_c2 ;
193  G4double fBetheVelocity = fine_structure_const * c_light / v;
194  coeff *= fine_structure_const * fine_structure_const * hbarc_squared / 
195           kineticEnergy ;
196
197  //G4double beta = std::sqrt( 2.0 * kineticEnergy / proton_mass_c2 );
198  //G4double fBetheVelocity = std::sqrt( 2.0 * 25.0 * keV / proton_mass_c2 )/beta;
199  //G4double coeff= twopi_mc2_rcl2*(material->GetElectronDensity())/(beta*beta);
200
201 
202  G4double l0Term = 0, l1Term = 0, l2Term = 0;
203 
204  for (G4int nos = 0 ; nos < nbOfShell ; nos++){
205   
206    G4double l0 = 0, l1 = 0, l2 = 0;
207    G4double NormalizedEnergy = ( 2.0 * electron_mass_c2 * v * v  ) / 
208                          ( c_squared * GetShellEnergy(material,nos) );
209    //G4double NormalizedEnergy = ( 2.0 * electron_mass_c2 * beta * beta  ) /
210    //                              GetShellEnergy(material,nos);
211    G4double shStrength = GetShellStrength(material,nos);
212
213    l0 = GetL0(NormalizedEnergy);
214    l0Term += shStrength  * l0; 
215   
216    l1 = GetL1(NormalizedEnergy);
217    l1Term += shStrength * l1; 
218   
219    l2 = GetL2(NormalizedEnergy);
220    l2Term += shStrength * l2; 
221
222    /*
223    if(material->GetName()=="Graphite"){
224    G4cout << nos << ". "
225           << " E(MeV)= " << kineticEnergy/MeV
226           << " normE= "  << NormalizedEnergy
227           << " sh en= "  << GetShellEnergy(material,nos)
228           << " str= "  << shStrength
229           << " v0/v= " << fBetheVelocity
230           << " l0= " << l0Term
231           << " l1= " << l1Term
232           << " l2= " << l2Term
233           << G4endl;
234           }
235     */   
236  }
237       
238  dedx = coeff * zParticle * zParticle * (l0Term
239       + zParticle * fBetheVelocity * l1Term
240       + zParticle * zParticle * fBetheVelocity * fBetheVelocity * l2Term);
241
242  //  G4cout << " E(MeV)= " << kineticEnergy/MeV
243  //         << " dedx(Mev/mm)= " << dedx*mm/MeV << G4endl;
244         
245  return dedx ; 
246                           
247}
248
249
250G4int G4QAOLowEnergyLoss::GetNumberOfShell(const G4Material* material) const
251{
252  // Set return value from table
253  G4int Z = (G4int)(material->GetZ());
254  G4int nShell = 0;
255
256  // Set return value if in material available from Aahrus
257  for(G4int i=0; i<numberOfMaterials; i++) {
258
259    if(materialAvailable[i] == Z){
260        nShell = nbofShellForMaterial[i];
261        break;
262    }
263    else nShell = fNumberOfShells[Z];
264  }
265 
266   return nShell;
267}
268
269
270
271G4double G4QAOLowEnergyLoss::GetShellEnergy(const G4Material* material,
272                                            G4int nbOfTheShell) const
273{
274    //
275    G4double shellEnergy = alShellEnergy[0];
276
277    if(material->GetZ() == 13) shellEnergy =  alShellEnergy[nbOfTheShell];
278    else if(material->GetZ() == 14)shellEnergy = siShellEnergy[nbOfTheShell];
279    else if(material->GetZ() == 29)shellEnergy = cuShellEnergy[nbOfTheShell]; 
280    else if(material->GetZ() == 73)shellEnergy =  taShellEnergy[nbOfTheShell];
281    else if(material->GetZ() == 79)shellEnergy =  auShellEnergy[nbOfTheShell];
282    else if(material->GetZ() == 78)shellEnergy =  ptShellEnergy[nbOfTheShell];
283    else if  (material->GetNumberOfElements() == 1)
284      shellEnergy = GetOscillatorEnergy(material, nbOfTheShell);
285    else G4cout << "WARNING - G4QAOLowEnergyLoss::GetShellEnergy - "
286                << "The model is not available for "
287                << material->GetName() 
288                << G4endl;
289
290  return  shellEnergy;
291}
292
293
294G4double G4QAOLowEnergyLoss::GetOscillatorEnergy(const G4Material* material,
295                                                 G4int nbOfTheShell) const
296{ 
297 
298  const G4Element* element = material->GetElement(0);
299 
300  G4int Z = (G4int)(element->GetZ());
301   
302G4double squaredPlasmonEnergy = 28.816 * 28.816  * 1e-6 
303                                * material->GetDensity()/g/cm3
304                                * (Z/element->GetN()) ;
305//G4double squaredPlasmonEnergy = 28.16 * 28.16 * 1e-6 
306//                                * (material->GetDensity()) * (cm3/g)
307//                              *  Z / (element->GetN()) ;
308 
309  G4double plasmonTerm = 0.66667 * GetOccupationNumber(Z,nbOfTheShell) 
310                       * squaredPlasmonEnergy / (Z*Z) ; 
311 
312  G4double ionTerm = std::exp(0.5) * (element->GetAtomicShell(nbOfTheShell)) ;
313
314  ionTerm = ionTerm*ionTerm ;
315   
316  G4double oscShellEnergy = std::sqrt( ionTerm + plasmonTerm );
317 
318/*  if(material->GetName()=="Graphite"){
319    G4cout << "\t" << Z
320           << "\t" << nbOfTheShell
321           << "\t" << squaredPlasmonEnergy
322           << "\t" << plasmonTerm
323           << "\t" << ionTerm
324           << "\t" << GetOccupationNumber(Z,nbOfTheShell)
325           << "\t" << element->GetAtomicShell(nbOfTheShell)
326           << "\t" << oscShellEnergy << G4endl;}
327*/
328  return  oscShellEnergy;
329}
330
331
332G4double G4QAOLowEnergyLoss::GetShellStrength(const G4Material* material,
333                                              G4int nbOfTheShell) const
334{
335  G4double shellStrength = alShellStrength[0];
336 
337  if(material->GetZ() == 13) shellStrength = alShellStrength[nbOfTheShell];
338  else if(material->GetZ() == 14)shellStrength =siShellStrength[nbOfTheShell];
339  else if(material->GetZ() == 29)shellStrength =cuShellStrength[nbOfTheShell];
340  else if(material->GetZ() == 73)shellStrength =taShellStrength[nbOfTheShell];
341  else if(material->GetZ() == 79)shellStrength =auShellStrength[nbOfTheShell];
342  else if(material->GetZ() == 78)shellStrength =ptShellStrength[nbOfTheShell];
343  else if  (material->GetNumberOfElements() == 1){
344    G4int Z = (G4int)(material->GetZ());
345    shellStrength = GetOccupationNumber(Z,nbOfTheShell) / Z ;}
346  else G4cout << "WARNING - G4QAOLowEnergyLoss::GetShellEnergy - "
347              << "The model is not available for "
348              << material->GetName() 
349              << G4endl;
350
351  return shellStrength;
352}
353
354G4double G4QAOLowEnergyLoss::GetOccupationNumber(G4int Z, G4int ShellNb) const
355{
356
357  G4int indice = ShellNb ;
358  for (G4int z = 1 ; z < Z ; z++) {indice += fNumberOfShells[z];}
359
360  return nbOfElectronPerSubShell[indice+1];
361}
362
363
364G4double G4QAOLowEnergyLoss::GetL0(G4double normEnergy) const 
365{
366  G4int n;
367 
368  for(n = 0; n < sizeL0; n++) {
369    if( normEnergy < L0[n][0] ) break;
370  }
371  if(0 == n) n = 1 ;
372  if(n >= sizeL0) n = sizeL0 - 1 ;
373
374  G4double l0    = L0[n][1];
375  G4double l0p   = L0[n-1][1];
376  G4double bethe = l0p + (l0 - l0p) * ( normEnergy - L0[n-1][0]) / 
377                  (L0[n][0] - L0[n-1][0]);
378  return bethe ;
379 
380}
381
382G4double G4QAOLowEnergyLoss::GetL1(G4double normEnergy) const
383{
384  G4int n;
385
386  for(n = 0; n < sizeL1; n++) {
387    if( normEnergy < L1[n][0] ) break;
388  }
389  if(0 == n) n = 1 ;
390  if(n >= sizeL1) n = sizeL1 - 1 ;
391
392  G4double l1    = L1[n][1];
393  G4double l1p   = L1[n-1][1];
394  G4double barkas= l1p + (l1 - l1p) * ( normEnergy - L1[n-1][0]) / 
395                  (L1[n][0] - L1[n-1][0]);
396 
397  return barkas;
398 
399}
400
401
402G4double G4QAOLowEnergyLoss::GetL2(G4double normEnergy) const
403{
404  G4int n;
405  for(n = 0; n < sizeL2; n++) {
406    if( normEnergy < L2[n][0] ) break;
407  }
408  if(0 == n) n = 1 ;
409  if(n >= sizeL2) n = sizeL2 - 1 ;
410
411  G4double l2    = L2[n][1];
412  G4double l2p   = L2[n-1][1];
413  G4double bloch = l2p + (l2 - l2p) * ( normEnergy - L2[n-1][0]) / 
414                  (L2[n][0] - L2[n-1][0]);
415 
416  return bloch;
417}
418
419
420const G4int G4QAOLowEnergyLoss::materialAvailable[6] = {13,14,29,73,79,78};
421  /*
422  "Aluminum",
423  "Silicon",
424  "Copper",
425  "Tantalum",
426  "Gold",
427  "Platinum"};
428  */
429const G4int G4QAOLowEnergyLoss::nbofShellForMaterial[6] = {3,3,4,6,6,6 };
430             
431G4double G4QAOLowEnergyLoss::alShellEnergy[3]  ={ 2795e-6, 202e-6,  16.9e-6};
432G4double G4QAOLowEnergyLoss::alShellStrength[3]={ 0.1349, 0.6387, 0.2264};
433G4double G4QAOLowEnergyLoss::siShellEnergy[3]  ={ 3179e-6, 249e-6, 20.3e-6 };
434G4double G4QAOLowEnergyLoss::siShellStrength[3]={ 0.1222, 0.5972, 0.2806};
435G4double G4QAOLowEnergyLoss::cuShellEnergy[4]  ={ 16931e-6, 1930e-6, 199e-6, 39.6e-6};
436G4double G4QAOLowEnergyLoss::cuShellStrength[4]={ 0.0505, 0.2561, 0.4913, 0.2021};
437G4double G4QAOLowEnergyLoss::taShellEnergy[6]  ={ 88926e-6, 18012e-6, 3210e-6, 575e-6, 108.7e-6, 30.8e-6};
438G4double G4QAOLowEnergyLoss::taShellStrength[6]={ 0.0126, 0.0896, 0.2599, 0.3413, 0.2057, 0.0908};
439G4double G4QAOLowEnergyLoss::auShellEnergy[6]={ 96235e-6, 25918e-6, 4116e-6, 599e-6, 87.3e-6, 36.9e-6};
440G4double G4QAOLowEnergyLoss::auShellStrength[6]={ 0.0139, 0.0803, 0.2473, 0.423, 0.1124, 0.1231};
441G4double G4QAOLowEnergyLoss::ptShellEnergy[6]={ 95017e-6, 25590e-6, 4063e-6, 576e-6, 81.9e-6, 31.4e-6};
442G4double G4QAOLowEnergyLoss::ptShellStrength[6]={ 0.0129, 0.0745, 0.2295, 0.4627, 0.1324, 0.0879};
443
444
445const G4double G4QAOLowEnergyLoss::L0[67][2] =
446{
447  {0.00,        0.000001},
448  {0.10,        0.000001},
449  {0.12,        0.00001},
450  {0.14,        0.00005},
451  {0.16,        0.00014},
452  {0.18,        0.00030},
453  {0.20,        0.00057},
454  {0.25,        0.00189},
455  {0.30,        0.00429},
456  {0.35,        0.00784},
457  {0.40,        0.01248},
458  {0.45,        0.01811},
459  {0.50,        0.02462},
460  {0.60,        0.03980},
461  {0.70,        0.05731},
462  {0.80,        0.07662},
463  {0.90,        0.09733},
464  {1.00,        0.11916},
465  {1.20,        0.16532},
466  {1.40,        0.21376},
467  {1.60,        0.26362},
468  {1.80,        0.31428},
469  {2.00,        0.36532},
470  {2.50,        0.49272},
471  {3.00,        0.61765},
472  {3.50,        0.73863},
473  {4.00,        0.85496},
474  {4.50,        0.96634},
475  {5.00,        1.07272},
476  {6.00,        1.27086},
477  {7.00,        1.45075},
478  {8.00,        1.61412},
479  {9.00,        1.76277},
480  {10.00,       1.89836},
481  {12.00,       2.13625},
482  {14.00,       2.33787},
483  {16.00,       2.51093},
484  {18.00,       2.66134},
485  {20.00,       2.79358},
486  {25.00, 3.06539},
487  {30.00, 3.27902},
488  {35.00, 3.45430},
489  {40.00, 3.60281},
490  {45.00, 3.73167},
491  {50.00, 3.84555},
492  {60.00, 4.04011},
493  {70.00, 4.20264},
494  {80.00, 4.34229},
495  {90.00,  4.46474},
496  {100.00, 4.57378},
497  {120.00, 4.76155},
498  {140.00, 4.91953},
499  {160.00, 5.05590},
500  {180.00, 5.17588},
501  {200.00, 5.28299},
502  {250.00, 5.50925},
503  {300.00, 5.69364},
504  {350.00, 5.84926},
505  {400.00, 5.98388},
506  {450.00, 6.10252},
507  {500.00, 6.20856},
508  {600.00, 6.39189},
509  {700.00, 6.54677},
510  {800.00, 6.68084},
511  {900.00, 6.79905},
512  {1000.00, 6.90474}
513};
514
515
516const G4double G4QAOLowEnergyLoss::L1[22][2] =
517{
518  {0.00,       -0.000001},
519  {0.10,       -0.00001},
520  {0.20,       -0.00049},
521  {0.30,       -0.00084},
522  {0.40,        0.00085},
523  {0.50,        0.00519},
524  {0.60,        0.01198},
525  {0.70,        0.02074},
526  {0.80,        0.03133},
527  {0.90,        0.04369},
528  {1.00,        0.06035},
529  {2.00,        0.24023},
530  {3.00,        0.44284},
531  {4.00,        0.62012},
532  {5.00,        0.77031},
533  {6.00,        0.90390},
534  {7.00,        1.02705},
535  {8.00,        1.10867},
536  {9.00,        1.17546},
537  {10.00,       1.21599},
538  {15.00,       1.24349},
539  {20.00,       1.16752}
540};
541
542
543const G4double G4QAOLowEnergyLoss::L2[14][2] =
544{
545  {0.00,        0.000001},
546  {0.10,        0.00001},
547  {0.20,        0.00000},
548  {0.40,       -0.00120},
549  {0.60,       -0.00036},
550  {0.80,        0.00372},
551  {1.00,        0.01298},
552  {2.00,        0.08296},
553  {4.00,        0.21953},
554  {6.00,        0.23903},
555  {8.00,        0.20893},
556  {10.00,       0.10879},
557  {20.00,      -0.88409},       
558  {40.00,      -1.13902}
559};
560
561
562const G4int G4QAOLowEnergyLoss::nbOfElectronPerSubShell[1540] = 
563{
564  0, // consistency with G4AtomicShells
565  1,//------ H
566  2,//------ He
567  2,  1,//------ Li
568  2,  2,//------ Be
569  2,  2,  1,//------ B
570  2,  2,  2,//------ C
571  2,  2,  2,  1,//------ N
572  2,  2,  2,  2,//------ O
573  2,  2,  5,//------ F
574  2,  2,  2,  4,//------ Ne
575  2,  2,  2,  4,  1,//------ Na
576  2,  2,  2,  4,  2,//------ Mg
577  2,  2,  2,  4,  2,  1,//------ Al
578  2,  2,  2,  4,  2,  2,//------ Si
579  2,  2,  2,  4,  2,  3,//------ P
580  2,  2,  2,  4,  2,  4,//------
581  2,  2,  2,  4,  2,  5,//------
582  2,  2,  2,  4,  2,  2,  4,//------
583  2,  2,  2,  4,  2,  2,  4,  1,//------
584  2,  2,  2,  4,  2,  2,  4,  2,//------
585  2,  2,  2,  4,  2,  2,  4,  1,  2,//------
586  2,  2,  2,  4,  2,  2,  4,  2,  2,//------
587  2,  2,  2,  4,  2,  2,  4,  3,  2,//------
588  2,  2,  2,  4,  2,  2,  4,  4,  2,//------
589  2,  2,  2,  4,  2,  2,  4,  5,  2,//------
590  2,  2,  2,  4,  2,  2,  4,  6,  2,//------
591  2,  2,  2,  4,  2,  2,  4,  7,  2,//------
592  2,  2,  2,  4,  2,  2,  4,  4,  4,  2,//------
593  2,  2,  2,  4,  2,  2,  4,  4,  5,  2,//------
594  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,//------
595  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  1,//------
596  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,//------
597  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  3,//------
598  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  4,//------
599  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  5,//------
600  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,//------
601  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  1,//------
602  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  2,//------
603  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  1,  2,//------
604  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  2,  2,//------
605  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  3,  2,//------
606  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  2,//------
607  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  5,  2,//------
608  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  6,  2,//------
609  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  7,  2,//------
610  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  4,  2,//------
611  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  5,  2,//------
612  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,//------
613  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  1,//------
614  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  2,//------
615  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  3,//------
616  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  4,//------
617  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  5,//------
618  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  2,  4,//------
619  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  2,  4,  1,//------
620  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  2,  4,  2,//------
621  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  2,  4,  1,  2,//------
622  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  2,  2,  2,  4,  2,//------
623  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  3,  2,  2,  4,  2,//------
624  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  4,  2,  2,  4,  2,//------
625  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  5,  2,  2,  4,  2,//------
626  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  2,  2,  4,  2,//------
627  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  7,  2,  2,  4,  2,//------
628  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  7,  2,  2,  4,  1,  2,//------
629  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  9,  2,  2,  4,  2,//------
630  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6, 10,  2,  2,  4,  2,//------
631  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6, 11,  2,  2,  4,  2,//------
632  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6, 12,  2,  2,  4,  2,//------
633  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6, 13,  2,  2,  4,  2,//------
634  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  2,//------
635  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  1,  2,//------
636  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  2,  2,//------
637  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  3,  2,//------
638  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  2,//------
639  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  5,  2,//------
640  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  6,  2,//------
641  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  7,  2,//------
642  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  9,  1,//------
643  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  1,//------
644  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,//------
645  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  1,//------
646  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,//------
647  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  3,//------
648  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  4,//------
649  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  3,//------
650  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  4,//------
651  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  4,  1,//------
652  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  4,  2,//------
653  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  4,  1,  2,//------
654  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  4,  2,  2,//------
655  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  2,  2,  2,  4,  1,  2,//------
656  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  3,  2,  2,  4,  1,  2,//------
657  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  4,  2,  2,  4,  1,  2,//------
658  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  6,  2,  2,  4,  2,//------
659  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  7,  2,  2,  4,  2,//------
660  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  7,  2,  2,  4,  1,  2,//------
661  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6,  8,  2,  2,  4,  1,  2,//------
662  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6, 10,  2,  2,  4,  2,//------
663  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6, 11,  2,  2,  4,  2,//------
664  2,  2,  2,  4,  2,  2,  4,  4,  6,  2,  2,  4,  4,  6,  6,  8,  2,  2,  4,  4,  6, 12,  2,  2,  4,  2 //-----
665};
666
667const G4int G4QAOLowEnergyLoss::fNumberOfShells[101] =
668{
669 0 ,  // nonexisting zero element
670
671 1 ,  1 ,  2 ,  2 ,  3 ,    3 ,  4 ,  4 ,  3 ,  4 ,  //  1 - 10
672
673 5 ,  5 ,  6 ,  6 ,  6 ,    6 ,  6 ,  7 ,  8 ,  8 ,  // 11 - 20
674
675 9 ,  9 ,  9 ,  9 ,  9 ,    9 ,  9 , 10 , 10 , 10 ,  // 21 - 30
676
67711 , 11 , 11 , 11 , 11 ,   12 , 13 , 13 , 14 , 14 ,  // 31 - 40
678
67914 , 14 , 14 , 14 , 14 ,   15 , 15 , 15 , 16 , 16 ,  // 41 - 50
680
681// ----------------------------------------------------------
682
68316 , 16 , 16 , 17 , 18 ,   18 , 19 , 19 , 19 , 19 ,  // 51 - 60
684
68519 , 19 , 19 , 20 , 19 ,   19 , 19 , 19 , 19 , 20 ,  // 61 - 70
686
68721 , 21 , 21 , 21 , 21 ,   21 , 21 , 21 , 22 , 22 ,  // 71 - 80
688
68923 , 23 , 23 , 23 , 24 ,   24 , 25 , 25 , 26 , 26 ,  // 81 - 90
690
69127 , 27 , 27 , 26 , 26 ,   27 , 27 , 26 , 26 , 26    // 91 - 100
692
693};
694
695
Note: See TracBrowser for help on using the repository browser.