Context Navigation

source: trunk/source/materials/src/G4MaterialStoppingICRU73.cc @ 1240

Last change on this file since 1240 was 1196, checked in by garnier, 15 years ago
update CVS release candidate geant4.9.3.01
File size: 161.6 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	// $Id: G4MaterialStoppingICRU73.cc,v 1.10 2009/11/09 16:51:07 vnivanch Exp $
28	// GEANT4 tag $Name: materials-V09-02-18 $
29
30	//---------------------------------------------------------------------------
31	//
32	// GEANT4 Class file
33	//
34	// Description: Data on stopping power
35	//
36	// Author: A.Ivantchenko 07.08.2008
37	//
38	// in the framework of the ESA Technology Research Programme
39	// (ESA contract 21435/08/NL/AT)
40	//
41	// Modifications:
42	// 29.04.2009 A.Ivantchenko added revised data for G4WATER, G4WATER_VAPOR,
43	// G4_NYLON-6/6 provided by Prof.P.Sigmund Univ. Southern Denmark
44	// 01.11.2009 A. Lechner: Extended energy range of G4_WATER tables up to
45	// 1 GeV/u
46	// 03.11.2009 A. Lechner: Changed material name from G4_NYLON-6/6 to
47	// G4_NYLON-6-6. Added new methods BuildPhysicsVector according
48	// to interface changes in base class G4VIonDEDXTable.
49	//
50	//----------------------------------------------------------------------------
51	//
52	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53
54	#include "G4MaterialStoppingICRU73.hh"
55	#include "G4LPhysicsFreeVector.hh"
56
57	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
58
59	G4MaterialStoppingICRU73::G4MaterialStoppingICRU73(G4bool splineFlag)
60	{
61	spline = splineFlag;
62
63	minIonAtomicNmb = 3;
64	maxIonAtomicNmb = 18;
65
66	G4String matName[31] = {"G4_A-150_TISSUE","G4_ADIPOSE_TISSUE_ICRP","G4_AIR","G4_ALUMINUM_OXIDE","G4_BONE_COMPACT_ICRU","G4_BONE_CORTICAL_ICRP","G4_C-552","G4_CALCIUM_FLUORIDE","G4_CARBON_DIOXIDE","G4_Pyrex_Glass","G4_KAPTON","G4_LITHIUM_FLUORIDE","G4_LITHIUM_TETRABORATE","G4_METHANE","G4_MUSCLE_STRIATED_ICRU","G4_NYLON-6-6","G4_PHOTO_EMULSION","G4_PLASTIC_SC_VINYLTOLUENE","G4_POLYCARBONATE","G4_POLYETHYLENE","G4_MYLAR","G4_LUCITE","G4_POLYSTYRENE","G4_TEFLON","G4_PROPANE","G4_SILICON_DIOXIDE","G4_SODIUM_IODIDE","G4_TISSUE-METHANE","G4_TISSUE-PROPANE","G4_WATER","G4_WATER_VAPOR"};
67
68	for( size_t i = 0; i < 31; i++ ) {
69	namesMat.push_back( matName[i] );
70	}
71	}
72
73	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
74
75	G4MaterialStoppingICRU73::~G4MaterialStoppingICRU73()
76	{
77	size_t nmb = dedx.size();
78
79	for(size_t i = 0; i < nmb; i++) delete dedx[i];
80
81	namesMat.clear();
82	dedxKeys.clear();
83	dedx.clear();
84	}
85
86	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
87	G4bool G4MaterialStoppingICRU73::BuildPhysicsVector(G4int ionZ,
88	G4int matZ)
89	{
90	return IsApplicable( ionZ, matZ );
91	}
92
93	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
94
95	G4bool G4MaterialStoppingICRU73::BuildPhysicsVector(G4int ionZ,
96	const G4String& matName)
97	{
98	return IsApplicable( ionZ, matName );
99	}
100
101	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
102
103	G4bool G4MaterialStoppingICRU73::IsApplicable(G4int, G4int)
104	{
105	return false;
106	}
107
108	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
109
110	G4bool G4MaterialStoppingICRU73::IsApplicable(G4int ionZ,
111	const G4String& matName)
112	{
113	G4bool isApplicable = false;
114
115	if( ionZ >= minIonAtomicNmb && ionZ <= maxIonAtomicNmb ) {
116
117	size_t nmbMat = namesMat.size();
118
119	for( size_t i = 0; i < nmbMat; i++ ) {
120
121	if( namesMat[i] == matName ) { isApplicable = true; break; }
122	}
123	}
124
125	return isApplicable;
126	}
127
128	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
129
130	G4PhysicsVector* G4MaterialStoppingICRU73::CreatePhysicsVector(
131	G4double* energy,
132	G4double* stoppower,
133	G4double factor)
134	{
135	G4LPhysicsFreeVector* pv = new G4LPhysicsFreeVector(31,energy[0],energy[30]);
136	pv->SetSpline(spline);
137
138	for( G4int i = 0; i < 31; i++ ) {
139	pv->PutValues(i,energy[i],stoppower[i]*factor);
140	}
141
142	return pv;
143	}
144
145	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
146
147	G4PhysicsVector* G4MaterialStoppingICRU73::CreatePhysicsVectorFullRange(
148	G4double* energy,
149	G4double* stoppower,
150	G4double factor)
151	{
152	G4LPhysicsFreeVector* pv = new G4LPhysicsFreeVector(53,energy[0],energy[52]);
153	pv->SetSpline(spline);
154
155	for( G4int i = 0; i < 53; i++ ) {
156	pv->PutValues(i,energy[i],stoppower[i]*factor);
157	}
158
159	return pv;
160	}
161
162	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
163
164	G4double G4MaterialStoppingICRU73::GetDEDX(G4double, G4int, G4int) {
165
166	G4double dedxValue = -1.0;
167
168	return dedxValue;
169	}
170
171	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
172
173	G4double G4MaterialStoppingICRU73::GetDEDX(G4double kinEnergyPerNucleon,
174	G4int ionZ,
175	const G4String& matName)
176	{
177	G4double dedxValue = -1.0;
178
179	G4PhysicsVector* physicsVector = GetPhysicsVector(ionZ, matName);
180
181	if(physicsVector != 0) {
182
183	size_t nmbBins = physicsVector -> GetVectorLength();
184	G4double lowerEnergyEdge = physicsVector -> GetLowEdgeEnergy(0);
185	G4double upperEnergyEdge = physicsVector -> GetLowEdgeEnergy(nmbBins-1);
186
187	if(kinEnergyPerNucleon <= upperEnergyEdge) {
188	G4bool b;
189
190	if(kinEnergyPerNucleon < lowerEnergyEdge) {
191	dedxValue = physicsVector -> GetValue(lowerEnergyEdge, b) *
192	std::sqrt(kinEnergyPerNucleon / lowerEnergyEdge);
193	}
194	else {
195	dedxValue = physicsVector -> GetValue(kinEnergyPerNucleon, b);
196	}
197	}
198	}
199
200	return dedxValue;
201	}
202
203	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
204
205	G4PhysicsVector* G4MaterialStoppingICRU73::GetPhysicsVector(
206	G4int ionZ,
207	const G4String& matName)
208	{
209	G4PhysicsVector* physicsVector = 0;
210
211	DEDXKey key = std::make_pair(ionZ, matName);
212
213	// Checking if dE/dx vector was already created
214	size_t nmbKeys = dedxKeys.size();
215	for( size_t i = 0; i < nmbKeys; i++ ) {
216
217	if( dedxKeys[i] == key ) { physicsVector = dedx[i]; break; }
218	}
219	if(physicsVector != 0) { return physicsVector; }
220
221	// Creating dE/dx vector if ion-material pair covered
222	G4double factor = MeV * cm2 / (0.001 * g);
223
224	G4double E[53] = {.025,.03,.04,.05,.06,.07,0.08,0.09,0.1,0.15,0.2,0.25,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.5,2,2.5,3,4,5,6,7,8,9,10,15,20,25,30,40,50,60,70,80,90,100,150,200,250,300,400,500,600,700,800,900,1000};
225
226	for( G4int i = 0; i < 53; i++ ) E[i] *= MeV;
227
228	if( matName == namesMat[0] ) {
229	if( ionZ == 3 ) {
230	G4double A_3[31]={2.748,2.992,3.386,3.674,3.877,4.016,4.108,4.166,4.2,4.186,4.059,3.906,3.751,3.458,3.2,2.974,2.777,2.603,2.449,2.313,1.807,1.485,1.263,1.1,0.8801,0.7372,0.6368,0.5623,0.5046,0.4586,0.4209};
231	physicsVector = CreatePhysicsVector(E,A_3,factor);
232	}
233	else if( ionZ == 4 ) {
234	G4double A_4[31]={3.391,3.718,4.275,4.722,5.07,5.333,5.527,5.668,5.769,5.944,5.889,5.762,5.61,5.291,4.986,4.706,4.45,4.219,4.008,3.817,3.075,2.571,2.21,1.938,1.56,1.31,1.132,1,0.8972,0.8151,0.748};
235	physicsVector = CreatePhysicsVector(E,A_4,factor);
236	}
237	else if( ionZ == 5 ) {
238	G4double A_5[31]={3.899,4.3,5.003,5.593,6.081,6.474,6.784,7.025,7.21,7.643,7.704,7.633,7.509,7.202,6.881,6.57,6.277,6.004,5.75,5.515,4.562,3.881,3.373,2.982,2.423,2.043,1.769,1.564,1.403,1.275,1.169};
239	physicsVector = CreatePhysicsVector(E,A_5,factor);
240	}
241	else if( ionZ == 6 ) {
242	G4double A_6[31]={4.3,4.767,5.598,6.316,6.933,7.453,7.883,8.233,8.514,9.262,9.476,9.487,9.409,9.143,8.827,8.504,8.189,7.888,7.603,7.333,6.205,5.361,4.712,4.2,3.448,2.925,2.542,2.251,2.022,1.837,1.686};
243	physicsVector = CreatePhysicsVector(E,A_6,factor);
244	}
245	else if( ionZ == 7 ) {
246	G4double A_7[31]={4.67,5.188,6.12,6.941,7.664,8.295,8.836,9.293,9.673,10.77,11.17,11.29,11.28,11.08,10.79,10.48,10.15,9.843,9.537,9.243,7.973,6.983,6.2,5.569,4.621,3.946,3.444,3.057,2.751,2.502,2.297};
247	physicsVector = CreatePhysicsVector(E,A_7,factor);
248	}
249	else if( ionZ == 8 ) {
250	G4double A_8[31]={4.999,5.564,6.589,7.5,8.319,9.049,9.693,10.25,10.73,12.2,12.82,13.07,13.14,13.04,12.79,12.5,12.18,11.86,11.54,11.23,9.852,8.734,7.828,7.082,5.938,5.105,4.475,3.984,3.591,3.271,3.005};
251	physicsVector = CreatePhysicsVector(E,A_8,factor);
252	}
253	else if( ionZ == 9 ) {
254	G4double A_9[31]={5.251,5.856,6.959,7.948,8.846,9.66,10.39,11.04,11.61,13.46,14.31,14.69,14.85,14.85,14.65,14.38,14.07,13.75,13.43,13.11,11.65,10.44,9.434,8.594,7.282,6.309,5.562,4.973,4.497,4.106,3.779};
255	physicsVector = CreatePhysicsVector(E,A_9,factor);
256	}
257	else if( ionZ == 10 ) {
258	G4double A_10[31]={5.481,6.119,7.29,8.347,9.312,10.19,11,11.74,12.39,14.62,15.71,16.25,16.5,16.61,16.48,16.24,15.95,15.64,15.32,15,13.47,12.17,11.08,10.15,8.689,7.581,6.721,6.034,5.475,5.012,4.623};
259	physicsVector = CreatePhysicsVector(E,A_10,factor);
260	}
261	else if( ionZ == 11 ) {
262	G4double A_11[31]={5.684,6.355,7.597,8.728,9.768,10.73,11.62,12.44,13.19,15.86,17.28,18.04,18.46,18.78,18.78,18.63,18.41,18.13,17.84,17.54,16.01,14.62,13.42,12.38,10.69,9.381,8.343,7.504,6.815,6.239,5.753};
263	physicsVector = CreatePhysicsVector(E,A_11,factor);
264	}
265	else if( ionZ == 12 ) {
266	G4double A_12[31]={5.906,6.584,7.845,9.005,10.07,11.07,12,12.87,13.66,16.63,18.27,19.18,19.69,20.11,20.15,20.01,19.79,19.51,19.21,18.9,17.31,15.88,14.63,13.55,11.78,10.41,9.323,8.435,7.7,7.083,6.558};
267	physicsVector = CreatePhysicsVector(E,A_12,factor);
268	}
269	else if( ionZ == 13 ) {
270	G4double A_13[31]={6.121,6.817,8.116,9.314,10.42,11.46,12.44,13.35,14.2,17.49,19.41,20.51,21.16,21.76,21.9,21.84,21.65,21.41,21.12,20.82,19.23,17.75,16.44,15.29,13.4,11.91,10.71,9.726,8.906,8.212,7.619};
271	physicsVector = CreatePhysicsVector(E,A_13,factor);
272	}
273	else if( ionZ == 14 ) {
274	G4double A_14[31]={6.335,7.047,8.378,9.61,10.75,11.83,12.84,13.8,14.69,18.27,20.45,21.74,22.53,23.31,23.57,23.57,23.44,23.23,22.96,22.67,21.09,19.58,18.22,17.01,15,13.4,12.1,11.03,10.13,9.368,8.711};
275	physicsVector = CreatePhysicsVector(E,A_14,factor);
276	}
277	else if( ionZ == 15 ) {
278	G4double A_15[31]={6.589,7.318,8.682,9.948,11.13,12.24,13.28,14.27,15.21,19.05,21.5,23,23.94,24.93,25.31,25.41,25.34,25.17,24.94,24.68,23.14,21.6,20.19,18.93,16.8,15.08,13.67,12.5,11.5,10.66,9.93};
279	physicsVector = CreatePhysicsVector(E,A_15,factor);
280	}
281	else if( ionZ == 16 ) {
282	G4double A_16[31]={6.793,7.544,8.947,10.25,11.47,12.61,13.69,14.72,15.7,19.79,22.5,24.22,25.31,26.52,27.05,27.24,27.24,27.12,26.93,26.69,25.22,23.67,22.23,20.92,18.67,16.84,15.32,14.04,12.96,12.03,11.22};
283	physicsVector = CreatePhysicsVector(E,A_16,factor);
284	}
285	else if( ionZ == 17 ) {
286	G4double A_17[31]={7.055,7.821,9.248,10.57,11.82,12.99,14.1,15.15,16.16,20.46,23.42,25.34,26.59,28.01,28.67,28.95,29.02,28.95,28.8,28.59,27.18,25.63,24.15,22.8,20.46,18.53,16.92,15.55,14.39,13.38,12.51};
287	physicsVector = CreatePhysicsVector(E,A_17,factor);
288	}
289	else if( ionZ == 18 ) {
290	G4double A_18[31]={7.196,7.978,9.43,10.77,12.04,13.24,14.37,15.44,16.47,20.95,24.14,26.26,27.67,29.29,30.09,30.46,30.59,30.57,30.45,30.27,28.91,27.36,25.87,24.48,22.07,20.06,18.38,16.94,15.72,14.65,13.72};
291	physicsVector = CreatePhysicsVector(E,A_18,factor);
292	}
293	}
294
295	else if( matName == namesMat[1] ) {
296	if( ionZ == 3 ) {
297	G4double AA_3[31]={2.852,3.097,3.49,3.778,3.983,4.125,4.221,4.283,4.321,4.313,4.182,4.022,3.859,3.553,3.284,3.05,2.845,2.665,2.507,2.366,1.846,1.516,1.289,1.123,0.898,0.7521,0.6496,0.5735,0.5147,0.4677,0.4292};
298	physicsVector = CreatePhysicsVector(E,AA_3,factor);
299	}
300	else if( ionZ == 4 ) {
301	G4double AA_4[31]={3.53,3.861,4.421,4.867,5.214,5.479,5.677,5.823,5.929,6.119,6.064,5.93,5.769,5.434,5.115,4.822,4.557,4.317,4.1,3.902,3.139,2.623,2.254,1.977,1.591,1.336,1.155,1.019,0.915,0.8313,0.7627};
302	physicsVector = CreatePhysicsVector(E,AA_4,factor);
303	}
304	else if( ionZ == 5 ) {
305	G4double AA_5[31]={4.062,4.473,5.185,5.776,6.262,6.655,6.968,7.213,7.405,7.862,7.928,7.854,7.721,7.396,7.057,6.731,6.426,6.142,5.879,5.635,4.655,3.957,3.438,3.04,2.469,2.082,1.804,1.594,1.431,1.3,1.192};
306	physicsVector = CreatePhysicsVector(E,AA_5,factor);
307	}
308	else if( ionZ == 6 ) {
309	G4double AA_6[31]={4.476,4.958,5.807,6.532,7.149,7.668,8.099,8.452,8.739,9.519,9.746,9.757,9.672,9.386,9.051,8.71,8.381,8.067,7.77,7.491,6.329,5.464,4.801,4.279,3.512,2.98,2.59,2.294,2.061,1.873,1.718};
310	physicsVector = CreatePhysicsVector(E,AA_6,factor);
311	}
312	else if( ionZ == 7 ) {
313	G4double AA_7[31]={4.857,5.394,6.352,7.185,7.913,8.544,9.085,9.543,9.927,11.06,11.48,11.61,11.59,11.38,11.07,10.73,10.39,10.06,9.744,9.439,8.128,7.114,6.314,5.672,4.706,4.019,3.508,3.115,2.803,2.55,2.341};
314	physicsVector = CreatePhysicsVector(E,AA_7,factor);
315	}
316	else if( ionZ == 8 ) {
317	G4double AA_8[31]={5.192,5.778,6.836,7.768,8.596,9.33,9.974,10.53,11.01,12.52,13.17,13.43,13.5,13.38,13.11,12.79,12.46,12.12,11.79,11.47,10.04,8.895,7.969,7.209,6.044,5.197,4.557,4.057,3.658,3.332,3.062};
318	physicsVector = CreatePhysicsVector(E,AA_8,factor);
319	}
320	else if( ionZ == 9 ) {
321	G4double AA_9[31]={5.447,6.073,7.214,8.232,9.145,9.967,10.7,11.35,11.92,13.81,14.69,15.09,15.25,15.24,15.02,14.73,14.4,14.06,13.73,13.39,11.88,10.63,9.606,8.75,7.413,6.423,5.664,5.064,4.581,4.183,3.85};
322	physicsVector = CreatePhysicsVector(E,AA_9,factor);
323	}
324	else if( ionZ == 10 ) {
325	G4double AA_10[31]={5.677,6.338,7.55,8.641,9.63,10.52,11.34,12.07,12.73,14.99,16.12,16.68,16.95,17.05,16.9,16.64,16.33,16,15.66,15.32,13.74,12.4,11.28,10.34,8.848,7.721,6.845,6.146,5.578,5.107,4.711};
326	physicsVector = CreatePhysicsVector(E,AA_10,factor);
327	}
328	else if( ionZ == 11 ) {
329	G4double AA_11[31]={5.885,6.579,7.866,9.036,10.1,11.08,11.99,12.81,13.56,16.27,17.73,18.53,18.96,19.28,19.27,19.1,18.85,18.55,18.24,17.92,16.33,14.9,13.67,12.6,10.87,9.542,8.487,7.634,6.934,6.349,5.854};
330	physicsVector = CreatePhysicsVector(E,AA_11,factor);
331	}
332	else if( ionZ == 12 ) {
333	G4double AA_12[31]={6.112,6.813,8.118,9.318,10.42,11.44,12.39,13.26,14.06,17.05,18.74,19.68,20.21,20.64,20.67,20.51,20.27,19.97,19.65,19.31,17.66,16.18,14.9,13.8,12,10.6,9.493,8.59,7.842,7.215,6.681};
334	physicsVector = CreatePhysicsVector(E,AA_12,factor);
335	}
336	else if( ionZ == 13 ) {
337	G4double AA_13[31]={6.332,7.05,8.393,9.633,10.78,11.85,12.84,13.77,14.63,17.94,19.9,21.05,21.72,22.33,22.47,22.38,22.18,21.91,21.61,21.29,19.63,18.1,16.75,15.58,13.64,12.12,10.9,9.904,9.069,8.364,7.761};
338	physicsVector = CreatePhysicsVector(E,AA_13,factor);
339	}
340	else if( ionZ == 14 ) {
341	G4double AA_14[31]={6.551,7.286,8.66,9.935,11.12,12.23,13.26,14.23,15.14,18.74,20.97,22.31,23.13,23.93,24.18,24.17,24.02,23.78,23.5,23.19,21.53,19.96,18.57,17.33,15.28,13.65,12.33,11.23,10.32,9.542,8.873};
342	physicsVector = CreatePhysicsVector(E,AA_14,factor);
343	}
344	else if( ionZ == 15 ) {
345	G4double AA_15[31]={6.812,7.564,8.97,10.27,11.5,12.64,13.72,14.73,15.68,19.54,22.04,23.59,24.57,25.59,25.97,26.05,25.96,25.77,25.52,25.23,23.62,22.02,20.58,19.29,17.11,15.35,13.92,12.72,11.71,10.85,10.11};
346	physicsVector = CreatePhysicsVector(E,AA_15,factor);
347	}
348	else if( ionZ == 16 ) {
349	G4double AA_16[31]={7.022,7.795,9.24,10.58,11.84,13.02,14.14,15.19,16.19,20.31,23.06,24.83,25.98,27.22,27.75,27.92,27.9,27.77,27.55,27.3,25.73,24.13,22.64,21.3,19.01,17.14,15.59,14.29,13.19,12.25,11.43};
350	physicsVector = CreatePhysicsVector(E,AA_16,factor);
351	}
352	else if( ionZ == 17 ) {
353	G4double AA_17[31]={7.291,8.078,9.546,10.91,12.2,13.41,14.55,15.64,16.67,21,24,25.98,27.28,28.74,29.41,29.68,29.73,29.64,29.47,29.24,27.74,26.13,24.61,23.22,20.83,18.86,17.22,15.83,14.65,13.62,12.73};
354	physicsVector = CreatePhysicsVector(E,AA_17,factor);
355	}
356	else if( ionZ == 18 ) {
357	G4double AA_18[31]={7.435,8.239,9.732,11.11,12.42,13.66,14.83,15.94,16.99,21.52,24.74,26.92,28.38,30.07,30.87,31.24,31.35,31.31,31.17,30.97,29.52,27.91,26.37,24.95,22.48,20.43,18.71,17.25,16,14.92,13.97};
358	physicsVector = CreatePhysicsVector(E,AA_18,factor);
359	}
360	}
361
362	else if( matName == namesMat[2] ) {
363	if( ionZ == 3 ) {
364	G4double B_3[31]={2.062,2.238,2.532,2.764,2.943,3.077,3.176,3.247,3.296,3.354,3.282,3.176,3.061,2.839,2.639,2.461,2.305,2.167,2.044,1.933,1.523,1.257,1.072,0.9372,0.752,0.6313,0.5463,0.4831,0.4341,0.3949,0.3627};
365	physicsVector = CreatePhysicsVector(E,B_3,factor);
366	}
367	else if( ionZ == 4 ) {
368	G4double B_4[31]={2.56,2.797,3.205,3.546,3.829,4.06,4.245,4.39,4.502,4.751,4.76,4.686,4.581,4.346,4.113,3.896,3.696,3.513,3.346,3.193,2.592,2.179,1.879,1.653,1.335,1.123,0.9729,0.8601,0.7724,0.7024,0.645};
369	physicsVector = CreatePhysicsVector(E,B_4,factor);
370	}
371	else if( ionZ == 5 ) {
372	G4double B_5[31]={2.959,3.253,3.766,4.204,4.582,4.905,5.178,5.404,5.588,6.086,6.216,6.203,6.13,5.915,5.676,5.439,5.212,4.998,4.799,4.612,3.848,3.292,2.872,2.547,2.076,1.755,1.522,1.347,1.21,1.1,1.009};
373	physicsVector = CreatePhysicsVector(E,B_5,factor);
374	}
375	else if( ionZ == 6 ) {
376	G4double B_6[31]={3.28,3.628,4.24,4.769,5.234,5.642,5.998,6.305,6.566,7.35,7.635,7.707,7.683,7.514,7.286,7.045,6.805,6.572,6.35,6.138,5.239,4.553,4.019,3.593,2.962,2.518,2.192,1.943,1.746,1.588,1.458};
377	physicsVector = CreatePhysicsVector(E,B_6,factor);
378	}
379	else if( ionZ == 7 ) {
380	G4double B_7[31]={3.554,3.946,4.645,5.254,5.796,6.28,6.712,7.094,7.428,8.517,8.989,9.172,9.216,9.119,8.923,8.692,8.451,8.21,7.974,7.746,6.741,5.941,5.298,4.775,3.979,3.407,2.978,2.646,2.382,2.168,1.99};
381	physicsVector = CreatePhysicsVector(E,B_7,factor);
382	}
383	else if( ionZ == 8 ) {
384	G4double B_8[31]={3.796,4.228,5.005,5.69,6.302,6.856,7.357,7.808,8.212,9.609,10.28,10.59,10.72,10.72,10.57,10.36,10.13,9.893,9.653,9.415,8.332,7.434,6.695,6.079,5.121,4.415,3.877,3.456,3.117,2.84,2.61};
385	physicsVector = CreatePhysicsVector(E,B_8,factor);
386	}
387	else if( ionZ == 9 ) {
388	G4double B_9[31]={4.002,4.467,5.312,6.064,6.739,7.354,7.915,8.427,8.892,10.58,11.47,11.91,12.13,12.23,12.13,11.95,11.73,11.5,11.25,11.01,9.875,8.901,8.082,7.389,6.291,5.466,4.828,4.321,3.91,3.572,3.288};
389	physicsVector = CreatePhysicsVector(E,B_9,factor);
390	}
391	else if( ionZ == 10 ) {
392	G4double B_10[31]={4.19,4.683,5.587,6.399,7.133,7.803,8.418,8.986,9.507,11.47,12.58,13.18,13.49,13.71,13.67,13.52,13.32,13.1,12.86,12.61,11.43,10.39,9.506,8.744,7.517,6.578,5.842,5.252,4.768,4.367,4.029};
393	physicsVector = CreatePhysicsVector(E,B_10,factor);
394	}
395	else if( ionZ == 11 ) {
396	G4double B_11[31]={4.337,4.854,5.813,6.684,7.477,8.206,8.878,9.502,10.08,12.34,13.71,14.5,14.94,15.33,15.4,15.33,15.18,14.98,14.77,14.54,13.38,12.3,11.34,10.51,9.131,8.049,7.183,6.477,5.894,5.404,4.988};
397	physicsVector = CreatePhysicsVector(E,B_11,factor);
398	}
399	else if( ionZ == 12 ) {
400	G4double B_12[31]={4.499,5.029,6.02,6.929,7.763,8.531,9.243,9.905,10.52,13,14.57,15.51,16.07,16.58,16.7,16.66,16.52,16.33,16.11,15.87,14.66,13.53,12.52,11.64,10.18,9.031,8.103,7.343,6.711,6.178,5.723};
401	physicsVector = CreatePhysicsVector(E,B_12,factor);
402	}
403	else if( ionZ == 13 ) {
404	G4double B_13[31]={4.647,5.193,6.215,7.16,8.035,8.846,9.599,10.3,10.96,13.65,15.45,16.57,17.25,17.93,18.16,18.17,18.08,17.91,17.71,17.49,16.28,15.12,14.07,13.14,11.57,10.32,9.31,8.47,7.765,7.166,6.653};
405	physicsVector = CreatePhysicsVector(E,B_13,factor);
406	}
407	else if( ionZ == 14 ) {
408	G4double B_14[31]={4.794,5.352,6.401,7.378,8.289,9.137,9.928,10.66,11.36,14.25,16.25,17.55,18.37,19.22,19.55,19.64,19.58,19.45,19.27,19.06,17.87,16.68,15.6,14.62,12.96,11.62,10.52,9.613,8.84,8.18,7.611};
409	physicsVector = CreatePhysicsVector(E,B_14,factor);
410	}
411	else if( ionZ == 15 ) {
412	G4double B_15[31]={4.946,5.515,6.588,7.595,8.539,9.422,10.24,11.02,11.75,14.83,17.04,18.51,19.48,20.53,20.99,21.16,21.17,21.08,20.93,20.75,19.6,18.41,17.29,16.28,14.53,13.09,11.9,10.9,10.05,9.324,8.69};
413	physicsVector = CreatePhysicsVector(E,B_15,factor);
414	}
415	else if( ionZ == 16 ) {
416	G4double B_16[31]={5.082,5.665,6.766,7.802,8.779,9.697,10.55,11.36,12.13,15.38,17.79,19.45,20.56,21.82,22.41,22.67,22.74,22.7,22.59,22.43,21.36,20.17,19.03,17.98,16.15,14.62,13.34,12.26,11.33,10.53,9.839};
417	physicsVector = CreatePhysicsVector(E,B_16,factor);
418	}
419	else if( ionZ == 17 ) {
420	G4double B_17[31]={5.242,5.836,6.957,8.015,9.017,9.963,10.85,11.69,12.48,15.89,18.48,20.31,21.57,23.04,23.77,24.11,24.25,24.25,24.18,24.05,23.04,21.86,20.7,19.62,17.71,16.11,14.75,13.6,12.6,11.73,10.98};
421	physicsVector = CreatePhysicsVector(E,B_17,factor);
422	}
423	else if( ionZ == 18 ) {
424	G4double B_18[31]={5.351,5.954,7.092,8.167,9.189,10.15,11.07,11.93,12.75,16.29,19.03,21.03,22.44,24.11,24.96,25.39,25.59,25.64,25.59,25.49,24.53,23.35,22.18,21.08,19.11,17.45,16.03,14.81,13.76,12.85,12.04};
425	physicsVector = CreatePhysicsVector(E,B_18,factor);
426	}
427	}
428
429	else if( matName == namesMat[3] ) {
430	if( ionZ == 3 ) {
431	G4double C_3[31]={1.355,1.466,1.653,1.806,1.931,2.032,2.115,2.181,2.234,2.369,2.391,2.366,2.32,2.204,2.082,1.966,1.859,1.76,1.671,1.59,1.276,1.066,0.9168,0.8057,0.6519,0.5505,0.4785,0.4246,0.3826,0.3489,0.3212};
432	physicsVector = CreatePhysicsVector(E,C_3,factor);
433	}
434	else if( ionZ == 4 ) {
435	G4double C_4[31]={1.73,1.882,2.144,2.364,2.553,2.714,2.851,2.967,3.065,3.355,3.462,3.485,3.466,3.368,3.24,3.105,2.973,2.847,2.728,2.617,2.166,1.843,1.603,1.419,1.156,0.979,0.8516,0.7556,0.6807,0.6205,0.5712};
436	physicsVector = CreatePhysicsVector(E,C_4,factor);
437	}
438	else if( ionZ == 5 ) {
439	G4double C_5[31]={2.04,2.23,2.562,2.845,3.093,3.311,3.503,3.671,3.816,4.292,4.511,4.603,4.627,4.576,4.464,4.328,4.185,4.043,3.905,3.773,3.207,2.777,2.444,2.182,1.796,1.528,1.332,1.183,1.066,0.9714,0.8937};
440	physicsVector = CreatePhysicsVector(E,C_5,factor);
441	}
442	else if( ionZ == 6 ) {
443	G4double C_6[31]={2.295,2.523,2.922,3.267,3.572,3.844,4.087,4.306,4.5,5.177,5.527,5.704,5.786,5.803,5.722,5.599,5.458,5.31,5.161,5.015,4.36,3.835,3.415,3.074,2.559,2.191,1.917,1.706,1.538,1.403,1.29};
444	physicsVector = CreatePhysicsVector(E,C_6,factor);
445	}
446	else if( ionZ == 7 ) {
447	G4double C_7[31]={2.515,2.773,3.232,3.633,3.99,4.312,4.604,4.869,5.109,5.992,6.491,6.769,6.923,7.029,6.998,6.901,6.771,6.626,6.475,6.322,5.603,4.998,4.497,4.081,3.436,2.963,2.603,2.323,2.098,1.915,1.762};
448	physicsVector = CreatePhysicsVector(E,C_7,factor);
449	}
450	else if( ionZ == 8 ) {
451	G4double C_8[31]={2.707,2.993,3.506,3.959,4.366,4.735,5.073,5.383,5.667,6.753,7.409,7.8,8.035,8.247,8.278,8.218,8.11,7.976,7.829,7.675,6.915,6.244,5.673,5.188,4.416,3.836,3.388,3.033,2.745,2.509,2.311};
452	physicsVector = CreatePhysicsVector(E,C_8,factor);
453	}
454	else if( ionZ == 9 ) {
455	G4double C_9[31]={2.876,3.185,3.746,4.245,4.696,5.108,5.487,5.838,6.162,7.438,8.249,8.755,9.075,9.399,9.497,9.478,9.395,9.276,9.136,8.985,8.201,7.48,6.852,6.309,5.429,4.753,4.222,3.795,3.447,3.157,2.914};
456	physicsVector = CreatePhysicsVector(E,C_9,factor);
457	}
458	else if( ionZ == 10 ) {
459	G4double C_10[31]={3.029,3.359,3.961,4.502,4.994,5.446,5.863,6.251,6.613,8.069,9.037,9.665,10.07,10.52,10.69,10.72,10.67,10.57,10.44,10.3,9.503,8.742,8.065,7.471,6.491,5.724,5.113,4.617,4.207,3.864,3.574};
460	physicsVector = CreatePhysicsVector(E,C_10,factor);
461	}
462	else if( ionZ == 11 ) {
463	G4double C_11[31]={3.151,3.498,4.137,4.716,5.247,5.736,6.191,6.615,7.012,8.647,9.778,10.54,11.05,11.63,11.89,11.98,11.97,11.9,11.78,11.65,10.86,10.07,9.354,8.713,7.637,6.781,6.089,5.52,5.046,4.647,4.306};
464	physicsVector = CreatePhysicsVector(E,C_11,factor);
465	}
466	else if( ionZ == 12 ) {
467	G4double C_12[31]={3.282,3.642,4.308,4.919,5.482,6.005,6.491,6.947,7.375,9.168,10.45,11.34,11.95,12.68,13.04,13.19,13.22,13.17,13.08,12.96,12.18,11.37,10.62,9.939,8.782,7.848,7.083,6.449,5.916,5.463,5.073};
468	physicsVector = CreatePhysicsVector(E,C_12,factor);
469	}
470	else if( ionZ == 13 ) {
471	G4double C_13[31]={3.4,3.772,4.463,5.1,5.691,6.244,6.762,7.248,7.705,9.646,11.08,12.1,12.81,13.7,14.16,14.38,14.46,14.45,14.38,14.28,13.53,12.71,11.93,11.21,9.982,8.972,8.136,7.436,6.844,6.336,5.898};
472	physicsVector = CreatePhysicsVector(E,C_13,factor);
473	}
474	else if( ionZ == 14 ) {
475	G4double C_14[31]={3.513,3.896,4.608,5.268,5.884,6.463,7.008,7.521,8.006,10.08,11.65,12.8,13.62,14.66,15.23,15.53,15.66,15.69,15.65,15.57,14.85,14.02,13.22,12.48,11.18,10.1,9.201,8.441,7.792,7.234,6.748};
476	physicsVector = CreatePhysicsVector(E,C_14,factor);
477	}
478	else if( ionZ == 15 ) {
479	G4double C_15[31]={3.621,4.013,4.743,5.423,6.063,6.666,7.237,7.776,8.287,10.49,12.2,13.48,14.41,15.63,16.32,16.71,16.91,16.99,16.98,16.93,16.29,15.47,14.65,13.88,12.52,11.37,10.4,9.571,8.86,8.243,7.705};
480	physicsVector = CreatePhysicsVector(E,C_15,factor);
481	}
482	else if( ionZ == 16 ) {
483	G4double C_16[31]={3.72,4.122,4.874,5.575,6.238,6.865,7.46,8.024,8.561,10.89,12.73,14.14,15.19,16.58,17.39,17.87,18.15,18.28,18.32,18.29,17.73,16.93,16.11,15.32,13.9,12.69,11.65,10.76,9.984,9.311,8.719};
484	physicsVector = CreatePhysicsVector(E,C_16,factor);
485	}
486	else if( ionZ == 17 ) {
487	G4double C_17[31]={3.832,4.242,5.008,5.726,6.406,7.054,7.669,8.255,8.814,11.26,13.22,14.75,15.91,17.47,18.41,18.99,19.33,19.52,19.6,19.61,19.13,18.35,17.52,16.72,15.25,13.98,12.88,11.93,11.1,10.38,9.736};
488	physicsVector = CreatePhysicsVector(E,C_17,factor);
489	}
490	else if( ionZ == 18 ) {
491	G4double C_18[31]={3.921,4.338,5.118,5.849,6.544,7.208,7.842,8.448,9.026,11.58,13.65,15.29,16.55,18.29,19.35,20.02,20.43,20.67,20.8,20.84,20.43,19.67,18.84,18.02,16.51,15.18,14.03,13.04,12.16,11.39,10.71};
492	physicsVector = CreatePhysicsVector(E,C_18,factor);
493	}
494	}
495
496	else if( matName == namesMat[4] ) {
497	if( ionZ == 3 ) {
498	G4double D_3[31]={2.204,2.403,2.723,2.962,3.137,3.265,3.355,3.418,3.459,3.491,3.406,3.291,3.171,2.939,2.731,2.548,2.387,2.244,2.117,2.004,1.581,1.307,1.117,0.9771,0.7854,0.6601,0.5718,0.5059,0.4547,0.4138,0.3803};
499	physicsVector = CreatePhysicsVector(E,D_3,factor);
500	}
501	else if( ionZ == 4 ) {
502	G4double D_4[31]={2.712,2.981,3.438,3.805,4.096,4.324,4.5,4.636,4.738,4.953,4.942,4.856,4.742,4.494,4.252,4.026,3.82,3.631,3.459,3.301,2.684,2.259,1.951,1.719,1.391,1.172,1.016,0.8994,0.8084,0.7355,0.6757};
503	physicsVector = CreatePhysicsVector(E,D_4,factor);
504	}
505	else if( ionZ == 5 ) {
506	G4double D_5[31]={3.108,3.438,4.018,4.502,4.904,5.236,5.506,5.725,5.9,6.355,6.459,6.43,6.345,6.115,5.863,5.616,5.381,5.16,4.954,4.762,3.975,3.404,2.974,2.64,2.158,1.827,1.588,1.406,1.264,1.15,1.057};
507	physicsVector = CreatePhysicsVector(E,D_5,factor);
508	}
509	else if( ionZ == 6 ) {
510	G4double D_6[31]={3.419,3.802,4.489,5.081,5.589,6.02,6.385,6.691,6.946,7.683,7.935,7.988,7.949,7.761,7.52,7.267,7.017,6.776,6.546,6.328,5.402,4.698,4.151,3.715,3.069,2.615,2.28,2.024,1.822,1.658,1.524};
511	physicsVector = CreatePhysicsVector(E,D_6,factor);
512	}
513	else if( ionZ == 7 ) {
514	G4double D_7[31]={3.699,4.122,4.892,5.573,6.171,6.693,7.146,7.536,7.87,8.913,9.345,9.503,9.529,9.41,9.198,8.955,8.703,8.452,8.208,7.972,6.937,6.116,5.46,4.925,4.113,3.528,3.09,2.75,2.479,2.259,2.077};
515	physicsVector = CreatePhysicsVector(E,D_7,factor);
516	}
517	else if( ionZ == 8 ) {
518	G4double D_8[31]={3.948,4.407,5.25,6.009,6.688,7.293,7.829,8.3,8.711,10.07,10.7,10.99,11.09,11.06,10.9,10.67,10.43,10.18,9.928,9.682,8.563,7.642,6.886,6.257,5.281,4.562,4.014,3.583,3.237,2.954,2.718};
519	physicsVector = CreatePhysicsVector(E,D_8,factor);
520	}
521	else if( ionZ == 9 ) {
522	G4double D_9[31]={4.148,4.637,5.539,6.363,7.111,7.788,8.396,8.94,9.423,11.09,11.94,12.35,12.54,12.61,12.5,12.3,12.07,11.82,11.57,11.32,10.15,9.147,8.308,7.602,6.483,5.644,4.994,4.477,4.058,3.711,3.421};
523	physicsVector = CreatePhysicsVector(E,D_9,factor);
524	}
525	else if( ionZ == 10 ) {
526	G4double D_10[31]={4.332,4.847,5.801,6.677,7.485,8.224,8.897,9.508,10.06,12.03,13.09,13.65,13.94,14.13,14.07,13.91,13.7,13.46,13.22,12.96,11.75,10.68,9.773,8.995,7.745,6.79,6.041,5.439,4.946,4.535,4.189};
527	physicsVector = CreatePhysicsVector(E,D_10,factor);
528	}
529	else if( ionZ == 11 ) {
530	G4double D_11[31]={4.494,5.033,6.04,6.97,7.838,8.642,9.383,10.06,10.69,13,14.34,15.1,15.53,15.9,15.95,15.87,15.71,15.5,15.28,15.04,13.83,12.71,11.72,10.86,9.432,8.319,7.429,6.705,6.106,5.604,5.177};
531	physicsVector = CreatePhysicsVector(E,D_11,factor);
532	}
533	else if( ionZ == 12 ) {
534	G4double D_12[31]={4.672,5.218,6.244,7.198,8.094,8.932,9.712,10.43,11.1,13.66,15.2,16.1,16.62,17.1,17.21,17.15,17,16.8,16.58,16.34,15.09,13.92,12.89,11.99,10.5,9.319,8.373,7.597,6.952,6.407,5.942};
535	physicsVector = CreatePhysicsVector(E,D_12,factor);
536	}
537	else if( ionZ == 13 ) {
538	G4double D_13[31]={4.842,5.403,6.458,7.44,8.367,9.243,10.06,10.83,11.55,14.36,16.13,17.2,17.86,18.51,18.72,18.72,18.62,18.44,18.24,18.01,16.76,15.57,14.49,13.53,11.93,10.66,9.617,8.759,8.039,7.428,6.903};
539	physicsVector = CreatePhysicsVector(E,D_13,factor);
540	}
541	else if( ionZ == 14 ) {
542	G4double D_14[31]={5.012,5.585,6.664,7.673,8.626,9.531,10.39,11.19,11.95,14.99,16.98,18.23,19.02,19.84,20.16,20.23,20.17,20.03,19.84,19.63,18.4,17.18,16.06,15.06,13.37,12,10.88,9.941,9.152,8.477,7.895};
543	physicsVector = CreatePhysicsVector(E,D_14,factor);
544	}
545	else if( ionZ == 15 ) {
546	G4double D_15[31]={5.2,5.786,6.889,7.923,8.9,9.832,10.72,11.56,12.35,15.61,17.82,19.26,20.19,21.2,21.64,21.8,21.79,21.7,21.54,21.35,20.17,18.94,17.8,16.75,14.96,13.5,12.28,11.26,10.39,9.648,9.002};
547	physicsVector = CreatePhysicsVector(E,D_15,factor);
548	}
549	else if( ionZ == 16 ) {
550	G4double D_16[31]={5.36,5.962,7.094,8.155,9.159,10.12,11.03,11.91,12.73,16.2,18.62,20.24,21.32,22.54,23.11,23.36,23.42,23.37,23.25,23.09,21.97,20.74,19.57,18.49,16.61,15.06,13.75,12.65,11.7,10.89,10.18};
551	physicsVector = CreatePhysicsVector(E,D_16,factor);
552	}
553	else if( ionZ == 17 ) {
554	G4double D_17[31]={5.554,6.169,7.319,8.4,9.423,10.4,11.34,12.24,13.09,16.74,19.36,21.16,22.38,23.8,24.5,24.83,24.96,24.96,24.88,24.74,23.69,22.47,21.28,20.17,18.22,16.58,15.19,14.02,13,12.12,11.35};
555	physicsVector = CreatePhysicsVector(E,D_17,factor);
556	}
557	else if( ionZ == 18 ) {
558	G4double D_18[31]={5.672,6.299,7.469,8.566,9.605,10.6,11.55,12.47,13.35,17.15,19.95,21.92,23.28,24.91,25.73,26.16,26.34,26.39,26.34,26.22,25.23,24.02,22.82,21.68,19.67,17.97,16.52,15.28,14.21,13.28,12.46};
559	physicsVector = CreatePhysicsVector(E,D_18,factor);
560	}
561	}
562
563	else if( matName == namesMat[5] ) {
564	if( ionZ == 3 ) {
565	G4double E_3[31]={2.041,2.227,2.528,2.754,2.923,3.048,3.138,3.201,3.244,3.289,3.217,3.113,3.003,2.789,2.596,2.425,2.275,2.141,2.021,1.915,1.515,1.256,1.074,0.9407,0.7572,0.6371,0.5522,0.4889,0.4397,0.4003,0.3679};
566	physicsVector = CreatePhysicsVector(E,E_3,factor);
567	}
568	else if( ionZ == 4 ) {
569	G4double E_4[31]={2.507,2.76,3.189,3.535,3.811,4.031,4.203,4.337,4.44,4.665,4.667,4.593,4.491,4.264,4.041,3.831,3.639,3.463,3.301,3.153,2.572,2.169,1.876,1.654,1.341,1.132,0.9817,0.8693,0.7817,0.7114,0.6538};
570	physicsVector = CreatePhysicsVector(E,E_4,factor);
571	}
572	else if( ionZ == 5 ) {
573	G4double E_5[31]={2.869,3.179,3.723,4.178,4.559,4.875,5.135,5.349,5.521,5.982,6.098,6.081,6.009,5.801,5.571,5.343,5.125,4.92,4.727,4.547,3.808,3.268,2.859,2.541,2.08,1.764,1.534,1.359,1.223,1.113,1.023};
574	physicsVector = CreatePhysicsVector(E,E_5,factor);
575	}
576	else if( ionZ == 6 ) {
577	G4double E_6[31]={3.155,3.512,4.157,4.714,5.192,5.601,5.95,6.245,6.492,7.226,7.489,7.554,7.528,7.364,7.145,6.914,6.683,6.46,6.246,6.042,5.174,4.51,3.991,3.576,2.959,2.524,2.203,1.957,1.762,1.605,1.475};
578	physicsVector = CreatePhysicsVector(E,E_6,factor);
579	}
580	else if( ionZ == 7 ) {
581	G4double E_7[31]={3.408,3.802,4.524,5.165,5.729,6.222,6.653,7.026,7.348,8.376,8.816,8.984,9.023,8.928,8.74,8.519,8.288,8.057,7.831,7.612,6.644,5.871,5.249,4.741,3.966,3.407,2.986,2.659,2.399,2.186,2.011};
582	physicsVector = CreatePhysicsVector(E,E_7,factor);
583	}
584	else if( ionZ == 8 ) {
585	G4double E_8[31]={3.634,4.062,4.85,5.564,6.204,6.775,7.283,7.731,8.126,9.455,10.09,10.38,10.5,10.5,10.35,10.15,9.932,9.702,9.471,9.243,8.2,7.334,6.619,6.023,5.092,4.405,3.879,3.465,3.133,2.859,2.632};
586	physicsVector = CreatePhysicsVector(E,E_8,factor);
587	}
588	else if( ionZ == 9 ) {
589	G4double E_9[31]={3.82,4.275,5.117,5.891,6.596,7.235,7.81,8.326,8.787,10.41,11.25,11.67,11.88,11.97,11.88,11.71,11.5,11.27,11.04,10.81,9.719,8.781,7.99,7.319,6.254,5.451,4.828,4.331,3.928,3.594,3.314};
590	physicsVector = CreatePhysicsVector(E,E_9,factor);
591	}
592	else if( ionZ == 10 ) {
593	G4double E_10[31]={3.991,4.47,5.359,6.182,6.943,7.64,8.277,8.855,9.378,11.29,12.34,12.91,13.21,13.41,13.38,13.24,13.06,12.84,12.62,12.38,11.26,10.26,9.4,8.663,7.473,6.56,5.842,5.263,4.789,4.393,4.059};
594	physicsVector = CreatePhysicsVector(E,E_10,factor);
595	}
596	else if( ionZ == 11 ) {
597	G4double E_11[31]={4.141,4.641,5.578,6.448,7.264,8.022,8.721,9.364,9.953,12.18,13.49,14.24,14.68,15.06,15.14,15.07,14.93,14.75,14.55,14.33,13.21,12.17,11.24,10.42,9.075,8.016,7.167,6.475,5.901,5.419,5.009};
598	physicsVector = CreatePhysicsVector(E,E_11,factor);
599	}
600	else if( ionZ == 12 ) {
601	G4double E_12[31]={4.304,4.812,5.769,6.662,7.505,8.297,9.034,9.718,10.35,12.8,14.31,15.21,15.74,16.24,16.37,16.33,16.2,16.03,15.83,15.61,14.45,13.37,12.4,11.54,10.13,9.002,8.096,7.352,6.731,6.207,5.758};
602	physicsVector = CreatePhysicsVector(E,E_12,factor);
603	}
604	else if( ionZ == 13 ) {
605	G4double E_13[31]={4.461,4.982,5.965,6.884,7.756,8.583,9.36,10.09,10.76,13.45,15.18,16.25,16.91,17.57,17.8,17.83,17.74,17.59,17.41,17.2,16.06,14.94,13.93,13.03,11.51,10.29,9.299,8.476,7.784,7.196,6.69};
606	physicsVector = CreatePhysicsVector(E,E_13,factor);
607	}
608	else if( ionZ == 14 ) {
609	G4double E_14[31]={4.617,5.149,6.155,7.098,7.993,8.848,9.658,10.42,11.14,14.04,15.98,17.22,18.01,18.84,19.17,19.26,19.23,19.11,18.94,18.75,17.63,16.49,15.45,14.5,12.89,11.59,10.52,9.621,8.862,8.214,7.653};
610	physicsVector = CreatePhysicsVector(E,E_14,factor);
611	}
612	else if( ionZ == 15 ) {
613	G4double E_15[31]={4.786,5.33,6.357,7.322,8.239,9.119,9.959,10.75,11.51,14.61,16.76,18.17,19.1,20.13,20.58,20.76,20.77,20.7,20.57,20.4,19.33,18.19,17.11,16.13,14.43,13.04,11.88,10.9,10.07,9.349,8.727};
614	physicsVector = CreatePhysicsVector(E,E_15,factor);
615	}
616	else if( ionZ == 16 ) {
617	G4double E_16[31]={4.931,5.49,6.543,7.534,8.475,9.379,10.25,11.07,11.86,15.16,17.5,19.1,20.17,21.39,21.98,22.24,22.32,22.29,22.2,22.05,21.04,19.91,18.82,17.8,16.03,14.54,13.3,12.24,11.34,10.55,9.865};
618	physicsVector = CreatePhysicsVector(E,E_16,factor);
619	}
620	else if( ionZ == 17 ) {
621	G4double E_17[31]={5.106,5.676,6.746,7.755,8.714,9.635,10.52,11.37,12.19,15.65,18.18,19.95,21.16,22.59,23.3,23.65,23.79,23.81,23.75,23.63,22.69,21.57,20.46,19.42,17.57,16.02,14.7,13.57,12.6,11.75,11.01};
622	physicsVector = CreatePhysicsVector(E,E_17,factor);
623	}
624	else if( ionZ == 18 ) {
625	G4double E_18[31]={5.216,5.797,6.886,7.91,8.884,9.818,10.72,11.59,12.43,16.04,18.74,20.67,22.02,23.64,24.47,24.91,25.12,25.18,25.15,25.06,24.18,23.06,21.94,20.88,18.98,17.36,15.98,14.8,13.77,12.87,12.09};
626	physicsVector = CreatePhysicsVector(E,E_18,factor);
627	}
628	}
629
630	else if( matName == namesMat[6] ) {
631	if( ionZ == 3 ) {
632	G4double F_3[31]={1.926,2.103,2.397,2.623,2.794,2.92,3.013,3.079,3.126,3.192,3.136,3.044,2.941,2.737,2.551,2.384,2.236,2.105,1.987,1.882,1.486,1.229,1.05,0.9178,0.7373,0.6196,0.5365,0.4747,0.4267,0.3883,0.3569};
633	physicsVector = CreatePhysicsVector(E,F_3,factor);
634	}
635	else if( ionZ == 4 ) {
636	G4double F_4[31]={2.371,2.606,3.017,3.357,3.634,3.855,4.029,4.165,4.271,4.519,4.545,4.488,4.397,4.186,3.972,3.77,3.582,3.409,3.251,3.105,2.528,2.129,1.838,1.618,1.309,1.102,0.955,0.8448,0.7592,0.6906,0.6344};
637	physicsVector = CreatePhysicsVector(E,F_4,factor);
638	}
639	else if( ionZ == 5 ) {
640	G4double F_5[31]={2.721,3.009,3.521,3.961,4.337,4.653,4.915,5.129,5.304,5.787,5.931,5.937,5.881,5.695,5.478,5.26,5.048,4.848,4.659,4.482,3.751,3.214,2.808,2.492,2.034,1.721,1.494,1.323,1.189,1.081,0.9929};
641	physicsVector = CreatePhysicsVector(E,F_5,factor);
642	}
643	else if( ionZ == 6 ) {
644	G4double F_6[31]={2.999,3.334,3.936,4.466,4.933,5.339,5.687,5.983,6.233,6.986,7.28,7.372,7.366,7.23,7.029,6.809,6.587,6.37,6.161,5.961,5.103,4.443,3.926,3.513,2.899,2.468,2.15,1.907,1.715,1.561,1.433};
645	physicsVector = CreatePhysicsVector(E,F_6,factor);
646	}
647	else if( ionZ == 7 ) {
648	G4double F_7[31]={3.243,3.614,4.29,4.894,5.437,5.922,6.35,6.724,7.048,8.093,8.565,8.766,8.83,8.769,8.602,8.395,8.174,7.951,7.731,7.516,6.56,5.791,5.171,4.665,3.893,3.336,2.919,2.595,2.338,2.129,1.956};
649	physicsVector = CreatePhysicsVector(E,F_7,factor);
650	}
651	else if( ionZ == 8 ) {
652	G4double F_8[31]={3.463,3.867,4.607,5.277,5.887,6.443,6.944,7.391,7.787,9.131,9.799,10.13,10.27,10.31,10.19,10.01,9.799,9.579,9.355,9.133,8.105,7.245,6.531,5.936,5.007,4.322,3.798,3.388,3.058,2.788,2.563};
653	physicsVector = CreatePhysicsVector(E,F_8,factor);
654	}
655	else if( ionZ == 9 ) {
656	G4double F_9[31]={3.645,4.075,4.87,5.595,6.263,6.877,7.441,7.953,8.413,10.05,10.91,11.37,11.61,11.75,11.68,11.53,11.34,11.12,10.9,10.68,9.598,8.667,7.878,7.21,6.147,5.347,4.727,4.234,3.834,3.504,3.228};
657	physicsVector = CreatePhysicsVector(E,F_9,factor);
658	}
659	else if( ionZ == 10 ) {
660	G4double F_10[31]={3.812,4.265,5.107,5.881,6.597,7.263,7.881,8.45,8.97,10.88,11.96,12.56,12.9,13.15,13.15,13.04,12.86,12.66,12.44,12.22,11.11,10.11,9.261,8.527,7.341,6.432,5.718,5.144,4.674,4.283,3.953};
661	physicsVector = CreatePhysicsVector(E,F_10,factor);
662	}
663	else if( ionZ == 11 ) {
664	G4double F_11[31]={3.948,4.421,5.31,6.132,6.899,7.618,8.292,8.921,9.504,11.73,13.06,13.85,14.32,14.76,14.87,14.83,14.71,14.54,14.35,14.14,13.05,12.02,11.1,10.29,8.949,7.893,7.048,6.358,5.787,5.308,4.901};
665	physicsVector = CreatePhysicsVector(E,F_11,factor);
666	}
667	else if( ionZ == 12 ) {
668	G4double F_12[31]={4.1,4.581,5.489,6.337,7.131,7.879,8.584,9.248,9.869,12.32,13.84,14.77,15.34,15.89,16.06,16.05,15.94,15.78,15.58,15.37,14.24,13.17,12.2,11.36,9.946,8.829,7.929,7.191,6.576,6.057,5.613};
669	physicsVector = CreatePhysicsVector(E,F_12,factor);
670	}
671	else if( ionZ == 13 ) {
672	G4double F_13[31]={4.243,4.737,5.669,6.544,7.369,8.148,8.886,9.587,10.25,12.93,14.66,15.76,16.46,17.18,17.45,17.51,17.44,17.31,17.13,16.93,15.81,14.71,13.71,12.81,11.31,10.09,9.109,8.292,7.606,7.024,6.524};
673	physicsVector = CreatePhysicsVector(E,F_13,factor);
674	}
675	else if( ionZ == 14 ) {
676	G4double F_14[31]={4.384,4.889,5.844,6.742,7.591,8.397,9.164,9.895,10.59,13.47,15.41,16.68,17.5,18.4,18.78,18.9,18.88,18.78,18.63,18.45,17.35,16.22,15.19,14.25,12.66,11.36,10.3,9.409,8.658,8.016,7.462};
677	physicsVector = CreatePhysicsVector(E,F_14,factor);
678	}
679	else if( ionZ == 15 ) {
680	G4double F_15[31]={4.538,5.054,6.03,6.951,7.825,8.656,9.45,10.21,10.93,14.01,16.15,17.59,18.56,19.65,20.16,20.36,20.41,20.35,20.23,20.07,19.03,17.9,16.84,15.86,14.18,12.79,11.64,10.67,9.844,9.133,8.517};
681	physicsVector = CreatePhysicsVector(E,F_15,factor);
682	}
683	else if( ionZ == 16 ) {
684	G4double F_16[31]={4.67,5.2,6.202,7.148,8.047,8.904,9.723,10.51,11.26,14.52,16.86,18.47,19.58,20.88,21.52,21.82,21.93,21.92,21.84,21.71,20.73,19.61,18.53,17.52,15.76,14.29,13.05,12,11.1,10.32,9.639};
685	physicsVector = CreatePhysicsVector(E,F_16,factor);
686	}
687	else if( ionZ == 17 ) {
688	G4double F_17[31]={4.832,5.372,6.391,7.353,8.271,9.148,9.989,10.8,11.57,14.98,17.51,19.29,20.53,22.03,22.8,23.19,23.37,23.41,23.36,23.26,22.35,21.24,20.14,19.11,17.28,15.73,14.42,13.3,12.33,11.49,10.75};
689	physicsVector = CreatePhysicsVector(E,F_17,factor);
690	}
691	else if( ionZ == 18 ) {
692	G4double F_18[31]={4.931,5.481,6.515,7.491,8.423,9.316,10.17,11,11.79,15.33,18.02,19.96,21.33,23.03,23.93,24.41,24.64,24.73,24.71,24.64,23.79,22.68,21.58,20.52,18.64,17.03,15.66,14.49,13.47,12.58,11.8};
693	physicsVector = CreatePhysicsVector(E,F_18,factor);
694	}
695	}
696
697	else if( matName == namesMat[7] ) {
698	if( ionZ == 3 ) {
699	G4double G_3[31]={1.514,1.638,1.842,2.004,2.133,2.236,2.317,2.379,2.427,2.52,2.5,2.443,2.372,2.226,2.088,1.963,1.851,1.75,1.659,1.577,1.265,1.057,0.9105,0.801,0.6489,0.5484,0.4768,0.4231,0.3813,0.3477,0.32};
700	physicsVector = CreatePhysicsVector(E,G_3,factor);
701	}
702	else if( ionZ == 4 ) {
703	G4double G_4[31]={1.913,2.085,2.375,2.615,2.815,2.982,3.122,3.237,3.331,3.582,3.636,3.612,3.556,3.409,3.253,3.102,2.96,2.829,2.707,2.594,2.143,1.824,1.587,1.406,1.148,0.9733,0.8473,0.7521,0.6777,0.6179,0.5687};
704	physicsVector = CreatePhysicsVector(E,G_4,factor);
705	}
706	else if( ionZ == 5 ) {
707	G4double G_5[31]={2.238,2.453,2.825,3.138,3.406,3.637,3.836,4.007,4.152,4.594,4.753,4.786,4.762,4.642,4.488,4.327,4.169,4.018,3.874,3.738,3.169,2.743,2.415,2.158,1.779,1.516,1.323,1.176,1.06,0.9664,0.8893};
708	physicsVector = CreatePhysicsVector(E,G_5,factor);
709	}
710	else if( ionZ == 6 ) {
711	G4double G_6[31]={2.506,2.761,3.209,3.592,3.925,4.217,4.474,4.7,4.898,5.549,5.838,5.948,5.97,5.897,5.761,5.603,5.44,5.278,5.12,4.968,4.304,3.783,3.369,3.034,2.529,2.169,1.9,1.692,1.527,1.394,1.283};
712	physicsVector = CreatePhysicsVector(E,G_6,factor);
713	}
714	else if( ionZ == 7 ) {
715	G4double G_7[31]={2.732,3.021,3.534,3.981,4.375,4.724,5.036,5.315,5.563,6.432,6.87,7.074,7.158,7.154,7.051,6.909,6.75,6.586,6.422,6.261,5.528,4.925,4.431,4.022,3.39,2.927,2.576,2.301,2.08,1.9,1.75};
716	physicsVector = CreatePhysicsVector(E,G_7,factor);
717	}
718	else if( ionZ == 8 ) {
719	G4double G_8[31]={2.934,3.251,3.824,4.329,4.779,5.183,5.547,5.876,6.173,7.259,7.855,8.168,8.324,8.406,8.348,8.231,8.086,7.927,7.764,7.599,6.818,6.147,5.582,5.105,4.349,3.783,3.345,2.998,2.717,2.485,2.291};
720	physicsVector = CreatePhysicsVector(E,G_8,factor);
721	}
722	else if( ionZ == 9 ) {
723	G4double G_9[31]={3.114,3.456,4.078,4.635,5.135,5.589,6.001,6.376,6.718,8.011,8.77,9.196,9.431,9.608,9.602,9.514,9.386,9.236,9.076,8.911,8.097,7.374,6.75,6.214,5.349,4.687,4.168,3.751,3.41,3.127,2.888};
724	physicsVector = CreatePhysicsVector(E,G_9,factor);
725	}
726	else if( ionZ == 10 ) {
727	G4double G_10[31]={3.28,3.643,4.309,4.911,5.457,5.956,6.412,6.831,7.216,8.707,9.63,10.18,10.5,10.78,10.84,10.79,10.68,10.54,10.39,10.23,9.397,8.63,7.956,7.367,6.401,5.649,5.05,4.565,4.163,3.827,3.542};
728	physicsVector = CreatePhysicsVector(E,G_10,factor);
729	}
730	else if( ionZ == 11 ) {
731	G4double G_11[31]={3.422,3.803,4.507,5.148,5.736,6.277,6.775,7.236,7.661,9.347,10.44,11.12,11.54,11.95,12.08,12.07,11.99,11.87,11.73,11.58,10.75,9.946,9.227,8.591,7.53,6.689,6.011,5.454,4.99,4.598,4.264};
732	physicsVector = CreatePhysicsVector(E,G_11,factor);
733	}
734	else if( ionZ == 12 ) {
735	G4double G_12[31]={3.568,3.963,4.697,5.371,5.994,6.572,7.107,7.604,8.067,9.931,11.19,12,12.52,13.06,13.27,13.31,13.27,13.17,13.04,12.9,12.07,11.24,10.49,9.811,8.666,7.745,6.994,6.372,5.849,5.405,5.023};
736	physicsVector = CreatePhysicsVector(E,G_12,factor);
737	}
738	else if( ionZ == 13 ) {
739	G4double G_13[31]={3.704,4.112,4.873,5.576,6.23,6.84,7.41,7.942,8.439,10.47,11.89,12.83,13.46,14.15,14.44,14.54,14.54,14.47,14.36,14.23,13.42,12.57,11.79,11.08,9.853,8.856,8.034,7.347,6.765,6.268,5.837};
740	physicsVector = CreatePhysicsVector(E,G_13,factor);
741	}
742	else if( ionZ == 14 ) {
743	G4double G_14[31]={3.835,4.255,5.039,5.767,6.447,7.086,7.686,8.25,8.779,10.97,12.54,13.62,14.35,15.18,15.57,15.73,15.77,15.73,15.64,15.53,14.75,13.89,13.08,12.34,11.05,9.98,9.092,8.344,7.707,7.158,6.681};
744	physicsVector = CreatePhysicsVector(E,G_14,factor);
745	}
746	else if( ionZ == 15 ) {
747	G4double G_15[31]={3.959,4.391,5.195,5.945,6.65,7.315,7.943,8.535,9.094,11.44,13.15,14.37,15.21,16.21,16.71,16.95,17.04,17.04,16.98,16.89,16.16,15.32,14.49,13.72,12.36,11.22,10.27,9.452,8.753,8.148,7.62};
748	physicsVector = CreatePhysicsVector(E,G_15,factor);
749	}
750	else if( ionZ == 16 ) {
751	G4double G_16[31]={4.077,4.521,5.348,6.12,6.848,7.538,8.193,8.812,9.398,11.88,13.74,15.09,16.05,17.23,17.83,18.15,18.3,18.34,18.32,18.25,17.59,16.76,15.92,15.13,13.71,12.51,11.49,10.61,9.853,9.192,8.613};
752	physicsVector = CreatePhysicsVector(E,G_16,factor);
753	}
754	else if( ionZ == 17 ) {
755	G4double G_17[31]={4.203,4.656,5.5,6.289,7.036,7.747,8.424,9.068,9.679,12.3,14.29,15.77,16.84,18.19,18.91,19.31,19.52,19.6,19.61,19.57,18.98,18.17,17.33,16.52,15.05,13.79,12.71,11.77,10.96,10.25,9.617};
756	physicsVector = CreatePhysicsVector(E,G_17,factor);
757	}
758	else if( ionZ == 18 ) {
759	G4double G_18[31]={4.313,4.776,5.638,6.443,7.208,7.938,8.635,9.299,9.933,12.67,14.79,16.39,17.58,19.09,19.93,20.41,20.67,20.81,20.85,20.83,20.31,19.51,18.66,17.83,16.31,15,13.86,12.88,12.02,11.26,10.59};
760	physicsVector = CreatePhysicsVector(E,G_18,factor);
761	}
762	}
763
764	else if( matName == namesMat[8] ) {
765	if( ionZ == 3 ) {
766	G4double H_3[31]={1.909,2.083,2.374,2.602,2.778,2.91,3.008,3.079,3.128,3.196,3.137,3.042,2.938,2.733,2.547,2.38,2.232,2.101,1.984,1.878,1.484,1.227,1.048,0.9166,0.7364,0.6188,0.5359,0.4741,0.4262,0.3879,0.3565};
767	physicsVector = CreatePhysicsVector(E,H_3,factor);
768	}
769	else if( ionZ == 4 ) {
770	G4double H_4[31]={2.35,2.582,2.986,3.325,3.605,3.832,4.014,4.157,4.267,4.525,4.547,4.487,4.395,4.182,3.967,3.765,3.577,3.404,3.246,3.1,2.525,2.127,1.837,1.617,1.308,1.101,0.9543,0.8441,0.7585,0.69,0.6338};
771	physicsVector = CreatePhysicsVector(E,H_4,factor);
772	}
773	else if( ionZ == 5 ) {
774	G4double H_5[31]={2.696,2.981,3.484,3.919,4.295,4.615,4.885,5.108,5.29,5.791,5.935,5.937,5.878,5.689,5.472,5.253,5.042,4.842,4.653,4.477,3.747,3.212,2.807,2.492,2.034,1.721,1.494,1.322,1.188,1.081,0.9923};
775	physicsVector = CreatePhysicsVector(E,H_5,factor);
776	}
777	else if( ionZ == 6 ) {
778	G4double H_6[31]={2.972,3.304,3.899,4.422,4.885,5.292,5.645,5.949,6.207,6.988,7.286,7.374,7.365,7.224,7.022,6.802,6.58,6.364,6.155,5.956,5.101,4.443,3.927,3.515,2.901,2.469,2.151,1.907,1.715,1.56,1.433};
779	physicsVector = CreatePhysicsVector(E,H_6,factor);
780	}
781	else if( ionZ == 7 ) {
782	G4double H_7[31]={3.208,3.579,4.25,4.848,5.386,5.869,6.3,6.679,7.011,8.091,8.573,8.771,8.831,8.765,8.596,8.389,8.169,7.947,7.727,7.514,6.562,5.795,5.177,4.671,3.898,3.34,2.922,2.598,2.34,2.13,1.957};
783	physicsVector = CreatePhysicsVector(E,H_7,factor);
784	}
785	else if( ionZ == 8 ) {
786	G4double H_8[31]={3.419,3.823,4.562,5.228,5.833,6.384,6.885,7.336,7.737,9.121,9.806,10.13,10.27,10.3,10.18,10,9.794,9.574,9.351,9.13,8.108,7.251,6.54,5.945,5.016,4.33,3.805,3.393,3.062,2.791,2.566};
787	physicsVector = CreatePhysicsVector(E,H_8,factor);
788	}
789	else if( ionZ == 9 ) {
790	G4double H_9[31]={3.605,4.037,4.835,5.56,6.223,6.833,7.395,7.907,8.372,10.05,10.94,11.4,11.63,11.76,11.69,11.54,11.34,11.13,10.9,10.68,9.605,8.676,7.889,7.221,6.158,5.357,4.735,4.241,3.84,3.509,3.231};
791	physicsVector = CreatePhysicsVector(E,H_9,factor);
792	}
793	else if( ionZ == 10 ) {
794	G4double H_10[31]={3.767,4.221,5.068,5.843,6.556,7.218,7.831,8.399,8.92,10.87,11.98,12.59,12.92,13.17,13.16,13.05,12.87,12.67,12.45,12.23,11.12,10.14,9.283,8.549,7.362,6.451,5.733,5.157,4.685,4.292,3.961};
795	physicsVector = CreatePhysicsVector(E,H_10,factor);
796	}
797	else if( ionZ == 11 ) {
798	G4double H_11[31]={3.897,4.371,5.263,6.088,6.854,7.567,8.235,8.858,9.438,11.69,13.05,13.84,14.31,14.73,14.83,14.79,14.66,14.5,14.3,14.1,13.01,11.99,11.08,10.28,8.943,7.892,7.049,6.36,5.79,5.311,4.904};
799	physicsVector = CreatePhysicsVector(E,H_11,factor);
800	}
801	else if( ionZ == 12 ) {
802	G4double H_12[31]={4.043,4.528,5.444,6.3,7.098,7.845,8.547,9.207,9.826,12.29,13.86,14.8,15.37,15.92,16.08,16.06,15.95,15.79,15.59,15.38,14.26,13.19,12.23,11.38,9.973,8.854,7.952,7.211,6.593,6.072,5.627};
803	physicsVector = CreatePhysicsVector(E,H_12,factor);
804	}
805	else if( ionZ == 13 ) {
806	G4double H_13[31]={4.179,4.676,5.618,6.503,7.334,8.116,8.855,9.552,10.21,12.9,14.68,15.8,16.5,17.21,17.48,17.52,17.45,17.32,17.15,16.95,15.83,14.74,13.74,12.84,11.34,10.12,9.136,8.316,7.628,7.044,6.541};
807	physicsVector = CreatePhysicsVector(E,H_13,factor);
808	}
809	else if( ionZ == 14 ) {
810	G4double H_14[31]={4.312,4.82,5.785,6.695,7.555,8.368,9.137,9.867,10.56,13.45,15.43,16.72,17.56,18.45,18.81,18.93,18.9,18.8,18.65,18.46,17.37,16.25,15.22,14.29,12.69,11.4,10.33,9.439,8.685,8.04,7.483};
811	physicsVector = CreatePhysicsVector(E,H_14,factor);
812	}
813	else if( ionZ == 15 ) {
814	G4double H_15[31]={4.45,4.968,5.953,6.887,7.774,8.615,9.415,10.18,10.9,13.98,16.16,17.64,18.61,19.7,20.19,20.39,20.43,20.37,20.25,20.09,19.05,17.94,16.88,15.91,14.23,12.84,11.68,10.71,9.878,9.164,8.545};
815	physicsVector = CreatePhysicsVector(E,H_15,factor);
816	}
817	else if( ionZ == 16 ) {
818	G4double H_16[31]={4.574,5.105,6.114,7.072,7.987,8.857,9.685,10.47,11.23,14.48,16.86,18.51,19.64,20.93,21.55,21.84,21.95,21.94,21.86,21.73,20.76,19.65,18.57,17.57,15.81,14.34,13.1,12.04,11.14,10.36,9.674};
819	physicsVector = CreatePhysicsVector(E,H_16,factor);
820	}
821	else if( ionZ == 17 ) {
822	G4double H_17[31]={4.722,5.263,6.289,7.266,8.201,9.094,9.946,10.76,11.54,14.94,17.5,19.33,20.59,22.09,22.85,23.23,23.4,23.44,23.39,23.29,22.39,21.29,20.2,19.17,17.34,15.79,14.48,13.35,12.38,11.54,10.8};
823	physicsVector = CreatePhysicsVector(E,H_17,factor);
824	}
825	else if( ionZ == 18 ) {
826	G4double H_18[31]={4.82,5.37,6.411,7.401,8.352,9.262,10.13,10.97,11.77,15.3,18.02,20,21.4,23.1,23.99,24.46,24.68,24.77,24.75,24.67,23.83,22.74,21.64,20.59,18.71,17.1,15.73,14.55,13.52,12.63,11.84};
827	physicsVector = CreatePhysicsVector(E,H_18,factor);
828	}
829	}
830
831	else if( matName == namesMat[9] ) {
832	if( ionZ == 3 ) {
833	G4double I_3[31]={1.707,1.866,2.117,2.301,2.438,2.539,2.615,2.67,2.71,2.774,2.74,2.674,2.597,2.437,2.284,2.145,2.02,1.907,1.805,1.714,1.366,1.136,0.9741,0.8543,0.6893,0.5809,0.5042,0.4469,0.4023,0.3666,0.3372};
834	physicsVector = CreatePhysicsVector(E,I_3,factor);
835	}
836	else if( ionZ == 4 ) {
837	G4double I_4[31]={2.09,2.31,2.678,2.965,3.192,3.371,3.514,3.626,3.715,3.934,3.972,3.941,3.881,3.724,3.555,3.389,3.232,3.085,2.949,2.823,2.319,1.964,1.703,1.504,1.222,1.033,0.8972,0.7952,0.7156,0.6518,0.5994};
838	physicsVector = CreatePhysicsVector(E,I_4,factor);
839	}
840	else if( ionZ == 5 ) {
841	G4double I_5[31]={2.389,2.659,3.13,3.515,3.829,4.088,4.302,4.479,4.625,5.041,5.183,5.211,5.185,5.061,4.899,4.725,4.551,4.383,4.223,4.071,3.434,2.96,2.597,2.313,1.898,1.612,1.403,1.245,1.12,1.02,0.9381};
842	physicsVector = CreatePhysicsVector(E,I_5,factor);
843	}
844	else if( ionZ == 6 ) {
845	G4double I_6[31]={2.629,2.94,3.5,3.976,4.376,4.714,5,5.243,5.449,6.09,6.361,6.466,6.489,6.42,6.281,6.113,5.936,5.757,5.581,5.412,4.669,4.088,3.629,3.259,2.704,2.311,2.019,1.794,1.617,1.473,1.354};
846	physicsVector = CreatePhysicsVector(E,I_6,factor);
847	}
848	else if( ionZ == 7 ) {
849	G4double I_7[31]={2.835,3.177,3.806,4.359,4.837,5.249,5.604,5.911,6.178,7.059,7.483,7.684,7.771,7.78,7.681,7.533,7.363,7.184,7.002,6.822,6.001,5.327,4.779,4.327,3.63,3.124,2.741,2.443,2.204,2.01,1.849};
850	physicsVector = CreatePhysicsVector(E,I_7,factor);
851	}
852	else if( ionZ == 8 ) {
853	G4double I_8[31]={3.025,3.392,4.077,4.696,5.245,5.726,6.147,6.517,6.843,7.97,8.559,8.87,9.032,9.134,9.089,8.973,8.82,8.648,8.467,8.284,7.408,6.658,6.029,5.5,4.666,4.044,3.566,3.188,2.883,2.633,2.424};
854	physicsVector = CreatePhysicsVector(E,I_8,factor);
855	}
856	else if( ionZ == 9 ) {
857	G4double I_9[31]={3.194,3.582,4.312,4.983,5.591,6.134,6.615,7.043,7.425,8.789,9.548,9.974,10.22,10.42,10.43,10.35,10.22,10.06,9.882,9.698,8.787,7.977,7.283,6.689,5.736,5.01,4.443,3.989,3.619,3.313,3.055};
858	physicsVector = CreatePhysicsVector(E,I_9,factor);
859	}
860	else if( ionZ == 10 ) {
861	G4double I_10[31]={3.352,3.759,4.526,5.24,5.899,6.496,7.032,7.515,7.948,9.545,10.48,11.03,11.36,11.67,11.75,11.71,11.61,11.46,11.3,11.12,10.18,9.324,8.573,7.922,6.858,6.034,5.38,4.852,4.417,4.054,3.746};
862	physicsVector = CreatePhysicsVector(E,I_10,factor);
863	}
864	else if( ionZ == 11 ) {
865	G4double I_11[31]={3.483,3.906,4.707,5.459,6.163,6.813,7.405,7.942,8.43,10.28,11.41,12.12,12.56,13.03,13.21,13.24,13.18,13.07,12.93,12.77,11.86,10.98,10.18,9.467,8.279,7.335,6.573,5.948,5.428,4.989,4.615};
866	physicsVector = CreatePhysicsVector(E,I_11,factor);
867	}
868	else if( ionZ == 12 ) {
869	G4double I_12[31]={3.626,4.061,4.882,5.657,6.391,7.076,7.708,8.287,8.816,10.86,12.16,12.98,13.52,14.1,14.35,14.42,14.38,14.29,14.15,13.99,13.06,12.14,11.3,10.55,9.285,8.277,7.457,6.779,6.212,5.732,5.319};
870	physicsVector = CreatePhysicsVector(E,I_12,factor);
871	}
872	else if( ionZ == 13 ) {
873	G4double I_13[31]={3.759,4.208,5.051,5.846,6.605,7.323,7.994,8.614,9.185,11.43,12.91,13.87,14.51,15.25,15.59,15.73,15.74,15.68,15.56,15.42,14.51,13.56,12.69,11.9,10.55,9.462,8.565,7.817,7.186,6.647,6.183};
874	physicsVector = CreatePhysicsVector(E,I_13,factor);
875	}
876	else if( ionZ == 14 ) {
877	G4double I_14[31]={3.889,4.35,5.213,6.027,6.805,7.55,8.253,8.908,9.516,11.94,13.59,14.7,15.45,16.34,16.78,16.99,17.05,17.02,16.94,16.81,15.93,14.97,14.07,13.24,11.82,10.66,9.687,8.874,8.183,7.589,7.075};
878	physicsVector = CreatePhysicsVector(E,I_14,factor);
879	}
880	else if( ionZ == 15 ) {
881	G4double I_15[31]={4.017,4.488,5.371,6.202,6.998,7.765,8.496,9.184,9.827,12.43,14.25,15.51,16.38,17.44,18,18.29,18.41,18.43,18.38,18.28,17.47,16.51,15.59,14.73,13.24,11.99,10.95,10.06,9.303,8.647,8.076};
882	physicsVector = CreatePhysicsVector(E,I_15,factor);
883	}
884	else if( ionZ == 16 ) {
885	G4double I_16[31]={4.134,4.62,5.527,6.377,7.191,7.978,8.733,9.45,10.13,12.9,14.89,16.29,17.28,18.51,19.2,19.57,19.77,19.84,19.83,19.76,19.02,18.08,17.14,16.26,14.7,13.38,12.26,11.31,10.48,9.766,9.137};
886	physicsVector = CreatePhysicsVector(E,I_16,factor);
887	}
888	else if( ionZ == 17 ) {
889	G4double I_17[31]={4.266,4.762,5.689,6.554,7.382,8.184,8.958,9.698,10.4,13.33,15.47,17.01,18.12,19.54,20.34,20.81,21.07,21.19,21.21,21.18,20.52,19.59,18.65,17.75,16.13,14.75,13.56,12.54,11.66,10.88,10.2};
890	physicsVector = CreatePhysicsVector(E,I_17,factor);
891	}
892	else if( ionZ == 18 ) {
893	G4double I_18[31]={4.361,4.869,5.813,6.692,7.531,8.343,9.13,9.889,10.61,13.68,15.96,17.63,18.86,20.44,21.36,21.91,22.23,22.4,22.47,22.46,21.86,20.95,20,19.08,17.42,15.98,14.75,13.68,12.74,11.92,11.2};
894	physicsVector = CreatePhysicsVector(E,I_18,factor);
895	}
896	}
897
898	else if( matName == namesMat[10] ) {
899	if( ionZ == 3 ) {
900	G4double J_3[31]={2.149,2.349,2.681,2.933,3.118,3.25,3.342,3.405,3.445,3.476,3.393,3.281,3.163,2.935,2.729,2.547,2.386,2.243,2.116,2.002,1.576,1.3,1.109,0.9689,0.777,0.6521,0.5641,0.4987,0.448,0.4075,0.3743};
901	physicsVector = CreatePhysicsVector(E,J_3,factor);
902	}
903	else if( ionZ == 4 ) {
904	G4double J_4[31]={2.637,2.902,3.364,3.746,4.054,4.294,4.477,4.616,4.718,4.931,4.923,4.841,4.733,4.492,4.254,4.031,3.825,3.637,3.464,3.306,2.684,2.255,1.944,1.709,1.379,1.16,1.004,0.8878,0.7972,0.7248,0.6655};
905	physicsVector = CreatePhysicsVector(E,J_4,factor);
906	}
907	else if( ionZ == 5 ) {
908	G4double J_5[31]={3.024,3.346,3.92,4.415,4.836,5.185,5.468,5.695,5.874,6.329,6.435,6.411,6.333,6.113,5.87,5.627,5.395,5.175,4.969,4.777,3.984,3.406,2.971,2.633,2.145,1.812,1.571,1.39,1.248,1.134,1.041};
909	physicsVector = CreatePhysicsVector(E,J_5,factor);
910	}
911	else if( ionZ == 6 ) {
912	G4double J_6[31]={3.334,3.707,4.38,4.974,5.497,5.949,6.332,6.651,6.914,7.658,7.912,7.969,7.939,7.764,7.533,7.287,7.041,6.803,6.574,6.356,5.424,4.711,4.156,3.713,3.058,2.599,2.261,2.004,1.801,1.638,1.503};
913	physicsVector = CreatePhysicsVector(E,J_6,factor);
914	}
915	else if( ionZ == 7 ) {
916	G4double J_7[31]={3.611,4.025,4.777,5.451,6.058,6.599,7.073,7.482,7.829,8.893,9.326,9.489,9.524,9.421,9.222,8.987,8.741,8.495,8.252,8.017,6.976,6.144,5.476,4.932,4.106,3.513,3.07,2.727,2.455,2.234,2.051};
917	physicsVector = CreatePhysicsVector(E,J_7,factor);
918	}
919	else if( ionZ == 8 ) {
920	G4double J_8[31]={3.86,4.31,5.135,5.88,6.56,7.179,7.735,8.228,8.658,10.05,10.68,10.97,11.09,11.08,10.92,10.71,10.48,10.23,9.99,9.746,8.623,7.689,6.919,6.278,5.283,4.552,3.996,3.56,3.211,2.925,2.688};
921	physicsVector = CreatePhysicsVector(E,J_8,factor);
922	}
923	else if( ionZ == 9 ) {
924	G4double J_9[31]={4.066,4.547,5.433,6.238,6.978,7.662,8.287,8.853,9.357,11.07,11.92,12.34,12.54,12.62,12.52,12.34,12.12,11.88,11.63,11.38,10.2,9.196,8.344,7.624,6.485,5.631,4.971,4.448,4.025,3.676,3.384};
925	physicsVector = CreatePhysicsVector(E,J_9,factor);
926	}
927	else if( ionZ == 10 ) {
928	G4double J_10[31]={4.255,4.762,5.702,6.56,7.353,8.091,8.776,9.406,9.978,12.01,13.08,13.64,13.94,14.14,14.1,13.95,13.75,13.52,13.27,13.02,11.8,10.72,9.805,9.014,7.741,6.771,6.011,5.402,4.904,4.491,4.143};
929	physicsVector = CreatePhysicsVector(E,J_10,factor);
930	}
931	else if( ionZ == 11 ) {
932	G4double J_11[31]={4.408,4.939,5.931,6.845,7.693,8.488,9.234,9.93,10.57,12.96,14.31,15.06,15.5,15.88,15.95,15.88,15.74,15.55,15.33,15.1,13.9,12.77,11.78,10.9,9.461,8.331,7.427,6.693,6.086,5.578,5.147};
933	physicsVector = CreatePhysicsVector(E,J_11,factor);
934	}
935	else if( ionZ == 12 ) {
936	G4double J_12[31]={4.579,5.119,6.133,7.075,7.954,8.779,9.557,10.29,10.97,13.62,15.17,16.07,16.6,17.09,17.22,17.17,17.03,16.84,16.62,16.39,15.14,13.96,12.92,12,10.49,9.297,8.338,7.554,6.902,6.353,5.884};
937	physicsVector = CreatePhysicsVector(E,J_12,factor);
938	}
939	else if( ionZ == 13 ) {
940	G4double J_13[31]={4.741,5.295,6.339,7.313,8.224,9.083,9.897,10.66,11.39,14.3,16.1,17.17,17.83,18.49,18.71,18.73,18.63,18.47,18.27,18.05,16.81,15.6,14.52,13.55,11.92,10.63,9.58,8.712,7.985,7.368,6.839};
941	physicsVector = CreatePhysicsVector(E,J_13,factor);
942	}
943	else if( ionZ == 14 ) {
944	G4double J_14[31]={4.901,5.468,6.538,7.539,8.48,9.369,10.21,11.01,11.77,14.91,16.94,18.19,18.98,19.81,20.14,20.22,20.17,20.04,19.87,19.66,18.44,17.21,16.08,15.07,13.35,11.96,10.83,9.885,9.087,8.408,7.822};
945	physicsVector = CreatePhysicsVector(E,J_14,factor);
946	}
947	else if( ionZ == 15 ) {
948	G4double J_15[31]={5.077,5.658,6.753,7.782,8.752,9.669,10.54,11.37,12.16,15.52,17.77,19.21,20.14,21.17,21.62,21.79,21.81,21.72,21.58,21.39,20.23,18.99,17.83,16.78,14.96,13.47,12.24,11.21,10.33,9.58,8.927};
949	physicsVector = CreatePhysicsVector(E,J_15,factor);
950	}
951	else if( ionZ == 16 ) {
952	G4double J_16[31]={5.227,5.824,6.95,8.008,9.008,9.954,10.85,11.71,12.54,16.09,18.57,20.2,21.28,22.51,23.09,23.35,23.43,23.4,23.29,23.14,22.04,20.81,19.63,18.54,16.63,15.05,13.72,12.6,11.65,10.82,10.1};
953	physicsVector = CreatePhysicsVector(E,J_16,factor);
954	}
955	else if( ionZ == 17 ) {
956	G4double J_17[31]={5.413,6.021,7.168,8.248,9.27,10.24,11.16,12.04,12.89,16.61,19.3,21.12,22.34,23.77,24.48,24.83,24.98,24.99,24.92,24.79,23.76,22.54,21.34,20.22,18.23,16.57,15.16,13.97,12.94,12.05,11.27};
957	physicsVector = CreatePhysicsVector(E,J_17,factor);
958	}
959	else if( ionZ == 18 ) {
960	G4double J_18[31]={5.524,6.145,7.31,8.406,9.446,10.43,11.37,12.27,13.14,16.99,19.88,21.86,23.23,24.85,25.69,26.13,26.33,26.39,26.35,26.25,25.28,24.07,22.85,21.7,19.66,17.93,16.47,15.21,14.13,13.18,12.35};
961	physicsVector = CreatePhysicsVector(E,J_18,factor);
962	}
963	}
964
965	else if( matName == namesMat[11] ) {
966	if( ionZ == 3 ) {
967	G4double K_3[31]={1.622,1.757,1.984,2.165,2.312,2.429,2.522,2.595,2.652,2.774,2.762,2.702,2.624,2.457,2.297,2.152,2.022,1.905,1.8,1.705,1.349,1.116,0.9538,0.8343,0.6707,0.5639,0.4886,0.4325,0.3889,0.354,0.3255};
968	physicsVector = CreatePhysicsVector(E,K_3,factor);
969	}
970	else if( ionZ == 4 ) {
971	G4double K_4[31]={2.038,2.224,2.541,2.806,3.029,3.218,3.376,3.507,3.616,3.918,3.997,3.981,3.922,3.757,3.577,3.402,3.237,3.084,2.943,2.813,2.294,1.933,1.67,1.471,1.19,1.003,0.8695,0.7696,0.6918,0.6295,0.5785};
972	physicsVector = CreatePhysicsVector(E,K_4,factor);
973	}
974	else if( ionZ == 5 ) {
975	G4double K_5[31]={2.368,2.601,3.005,3.348,3.645,3.903,4.127,4.32,4.486,5.003,5.207,5.26,5.241,5.109,4.932,4.746,4.561,4.384,4.217,4.059,3.403,2.918,2.551,2.264,1.85,1.566,1.36,1.205,1.083,0.9854,0.9052};
976	physicsVector = CreatePhysicsVector(E,K_5,factor);
977	}
978	else if( ionZ == 6 ) {
979	G4double K_6[31]={2.629,2.905,3.391,3.81,4.178,4.503,4.791,5.046,5.27,6.024,6.377,6.522,6.559,6.483,6.326,6.142,5.95,5.76,5.575,5.397,4.628,4.032,3.565,3.192,2.636,2.245,1.957,1.736,1.562,1.422,1.306};
980	physicsVector = CreatePhysicsVector(E,K_6,factor);
981	}
982	else if( ionZ == 7 ) {
983	G4double K_7[31]={2.856,3.166,3.721,4.207,4.639,5.025,5.373,5.686,5.966,6.963,7.487,7.744,7.854,7.86,7.74,7.571,7.382,7.187,6.993,6.803,5.947,5.254,4.695,4.237,3.538,3.034,2.655,2.362,2.129,1.939,1.782};
984	physicsVector = CreatePhysicsVector(E,K_7,factor);
985	}
986	else if( ionZ == 8 ) {
987	G4double K_8[31]={3.056,3.396,4.013,4.561,5.053,5.498,5.902,6.269,6.603,7.841,8.547,8.93,9.126,9.234,9.166,9.024,8.849,8.658,8.462,8.265,7.346,6.572,5.928,5.39,4.55,3.929,3.455,3.083,2.784,2.539,2.335};
988	physicsVector = CreatePhysicsVector(E,K_8,factor);
989	}
990	else if( ionZ == 9 ) {
991	G4double K_9[31]={3.218,3.582,4.251,4.853,5.398,5.894,6.349,6.765,7.148,8.61,9.495,10.01,10.3,10.51,10.5,10.39,10.23,10.05,9.855,9.656,8.694,7.857,7.147,6.544,5.584,4.861,4.299,3.853,3.491,3.191,2.941};
992	physicsVector = CreatePhysicsVector(E,K_9,factor);
993	}
994	else if( ionZ == 10 ) {
995	G4double K_10[31]={3.364,3.749,4.461,5.11,5.702,6.246,6.747,7.21,7.637,9.314,10.38,11.03,11.42,11.76,11.81,11.74,11.6,11.43,11.24,11.04,10.05,9.166,8.398,7.737,6.667,5.846,5.199,4.68,4.254,3.9,3.602};
996	physicsVector = CreatePhysicsVector(E,K_10,factor);
997	}
998	else if( ionZ == 11 ) {
999	G4double K_11[31]={3.476,3.878,4.632,5.327,5.97,6.565,7.116,7.629,8.107,10.03,11.31,12.14,12.66,13.18,13.35,13.35,13.27,13.13,12.97,12.79,11.83,10.9,10.08,9.351,8.14,7.185,6.419,5.793,5.275,4.84,4.47};
1000	physicsVector = CreatePhysicsVector(E,K_11,factor);
1001	}
1002	else if( ionZ == 12 ) {
1003	G4double K_12[31]={3.604,4.015,4.786,5.505,6.174,6.798,7.38,7.924,8.434,10.52,11.96,12.91,13.53,14.18,14.41,14.44,14.37,14.24,14.08,13.9,12.89,11.93,11.07,10.31,9.035,8.027,7.213,6.544,5.987,5.517,5.115};
1004	physicsVector = CreatePhysicsVector(E,K_12,factor);
1005	}
1006	else if( ionZ == 13 ) {
1007	G4double K_13[31]={3.723,4.144,4.937,5.679,6.378,7.034,7.649,8.228,8.772,11.03,12.64,13.75,14.49,15.31,15.65,15.75,15.72,15.62,15.48,15.31,14.32,13.33,12.43,11.63,10.27,9.178,8.286,7.547,6.926,6.398,5.945};
1008	physicsVector = CreatePhysicsVector(E,K_13,factor);
1009	}
1010	else if( ionZ == 14 ) {
1011	G4double K_14[31]={3.84,4.27,5.081,5.844,6.566,7.248,7.893,8.502,9.077,11.5,13.27,14.52,15.39,16.38,16.83,17,17.02,16.95,16.83,16.67,15.7,14.7,13.77,12.93,11.5,10.33,9.367,8.563,7.883,7.301,6.798};
1012	physicsVector = CreatePhysicsVector(E,K_14,factor);
1013	}
1014	else if( ionZ == 15 ) {
1015	G4double K_15[31]={3.971,4.409,5.238,6.02,6.764,7.472,8.144,8.782,9.387,11.96,13.89,15.29,16.29,17.48,18.05,18.31,18.39,18.36,18.27,18.14,17.22,16.22,15.27,14.39,12.88,11.63,10.59,9.708,8.96,8.316,7.758};
1016	physicsVector = CreatePhysicsVector(E,K_15,factor);
1017	}
1018	else if( ionZ == 16 ) {
1019	G4double K_16[31]={4.08,4.53,5.379,6.182,6.95,7.683,8.382,9.047,9.68,12.41,14.49,16.04,17.16,18.55,19.26,19.61,19.76,19.78,19.73,19.62,18.77,17.77,16.8,15.9,14.31,12.98,11.87,10.91,10.1,9.394,8.779};
1020	physicsVector = CreatePhysicsVector(E,K_16,factor);
1021	}
1022	else if( ionZ == 17 ) {
1023	G4double K_17[31]={4.22,4.677,5.538,6.353,7.136,7.889,8.609,9.298,9.956,12.81,15.04,16.72,17.97,19.55,20.4,20.84,21.05,21.12,21.1,21.02,20.23,19.24,18.26,17.34,15.69,14.29,13.11,12.1,11.22,10.46,9.792};
1024	physicsVector = CreatePhysicsVector(E,K_17,factor);
1025	}
1026	else if( ionZ == 18 ) {
1027	G4double K_18[31]={4.299,4.762,5.633,6.458,7.252,8.018,8.755,9.463,10.14,13.11,15.46,17.27,18.64,20.41,21.38,21.91,22.18,22.3,22.31,22.26,21.52,20.55,19.55,18.61,16.92,15.47,14.24,13.18,12.25,11.45,10.74};
1028	physicsVector = CreatePhysicsVector(E,K_18,factor);
1029	}
1030	}
1031
1032	else if( matName == namesMat[12] ) {
1033	if( ionZ == 3 ) {
1034	G4double L_3[31]={1.968,2.15,2.444,2.663,2.824,2.943,3.029,3.091,3.133,3.185,3.123,3.027,2.922,2.717,2.53,2.363,2.215,2.083,1.966,1.861,1.466,1.211,1.033,0.9028,0.7245,0.6083,0.5265,0.4656,0.4184,0.3807,0.3498};
1035	physicsVector = CreatePhysicsVector(E,L_3,factor);
1036	}
1037	else if( ionZ == 4 ) {
1038	G4double L_4[31]={2.403,2.652,3.075,3.412,3.679,3.89,4.055,4.184,4.283,4.511,4.526,4.464,4.371,4.157,3.942,3.739,3.55,3.377,3.218,3.072,2.497,2.1,1.811,1.593,1.286,1.083,0.9375,0.8289,0.7446,0.6771,0.6219};
1039	physicsVector = CreatePhysicsVector(E,L_4,factor);
1040	}
1041	else if( ionZ == 5 ) {
1042	G4double L_5[31]={2.739,3.045,3.58,4.026,4.396,4.701,4.952,5.157,5.324,5.779,5.908,5.906,5.846,5.656,5.438,5.218,5.006,4.804,4.615,4.437,3.706,3.171,2.768,2.454,2.001,1.691,1.467,1.298,1.166,1.06,0.9734};
1043	physicsVector = CreatePhysicsVector(E,L_5,factor);
1044	}
1045	else if( ionZ == 6 ) {
1046	G4double L_6[31]={3.003,3.358,3.992,4.538,5.005,5.403,5.74,6.024,6.262,6.98,7.255,7.335,7.324,7.181,6.978,6.756,6.533,6.314,6.104,5.904,5.045,4.386,3.872,3.462,2.853,2.426,2.112,1.872,1.683,1.531,1.405};
1047	physicsVector = CreatePhysicsVector(E,L_6,factor);
1048	}
1049	else if( ionZ == 7 ) {
1050	G4double L_7[31]={3.237,3.628,4.341,4.97,5.522,6.003,6.421,6.781,7.092,8.093,8.541,8.726,8.781,8.712,8.541,8.332,8.109,7.884,7.662,7.447,6.488,5.72,5.102,4.599,3.832,3.28,2.868,2.548,2.294,2.088,1.918};
1051	physicsVector = CreatePhysicsVector(E,L_7,factor);
1052	}
1053	else if( ionZ == 8 ) {
1054	G4double L_8[31]={3.448,3.869,4.647,5.349,5.977,6.536,7.03,7.464,7.845,9.135,9.774,10.08,10.22,10.24,10.12,9.934,9.723,9.5,9.274,9.051,8.019,7.157,6.446,5.853,4.931,4.251,3.733,3.327,3.002,2.735,2.514};
1055	physicsVector = CreatePhysicsVector(E,L_8,factor);
1056	}
1057	else if( ionZ == 9 ) {
1058	G4double L_9[31]={3.63,4.073,4.905,5.665,6.357,6.982,7.543,8.043,8.489,10.06,10.89,11.33,11.55,11.67,11.6,11.45,11.25,11.03,10.81,10.58,9.495,8.562,7.775,7.109,6.053,5.26,4.646,4.158,3.764,3.438,3.166};
1059	physicsVector = CreatePhysicsVector(E,L_9,factor);
1060	}
1061	else if( ionZ == 10 ) {
1062	G4double L_10[31]={3.799,4.263,5.137,5.946,6.692,7.375,7.996,8.558,9.064,10.91,11.94,12.52,12.83,13.07,13.06,12.94,12.76,12.55,12.33,12.1,10.98,9.992,9.139,8.407,7.228,6.326,5.619,5.051,4.587,4.202,3.877};
1063	physicsVector = CreatePhysicsVector(E,L_10,factor);
1064	}
1065	else if( ionZ == 11 ) {
1066	G4double L_11[31]={3.933,4.415,5.328,6.184,6.983,7.724,8.407,9.032,9.603,11.75,13.03,13.78,14.23,14.64,14.74,14.69,14.57,14.4,14.2,13.99,12.89,11.86,10.95,10.15,8.813,7.768,6.931,6.25,5.686,5.213,4.812};
1067	physicsVector = CreatePhysicsVector(E,L_11,factor);
1068	}
1069	else if( ionZ == 12 ) {
1070	G4double L_12[31]={4.086,4.576,5.507,6.386,7.213,7.988,8.709,9.376,9.991,12.35,13.82,14.72,15.26,15.79,15.94,15.92,15.81,15.64,15.44,15.23,14.08,13,12.04,11.2,9.795,8.687,7.795,7.065,6.457,5.945,5.508};
1071	physicsVector = CreatePhysicsVector(E,L_12,factor);
1072	}
1073	else if( ionZ == 13 ) {
1074	G4double L_13[31]={4.23,4.732,5.684,6.586,7.441,8.25,9.011,9.721,10.38,12.97,14.65,15.71,16.38,17.07,17.32,17.36,17.29,17.15,16.97,16.77,15.64,14.53,13.53,12.64,11.14,9.934,8.956,8.148,7.471,6.896,6.402};
1075	physicsVector = CreatePhysicsVector(E,L_13,factor);
1076	}
1077	else if( ionZ == 14 ) {
1078	G4double L_14[31]={4.374,4.887,5.858,6.779,7.658,8.495,9.287,10.03,10.73,13.53,15.41,16.63,17.42,18.28,18.64,18.75,18.72,18.61,18.46,18.27,17.15,16.03,14.99,14.06,12.47,11.18,10.13,9.246,8.504,7.87,7.323};
1079	physicsVector = CreatePhysicsVector(E,L_14,factor);
1080	}
1081	else if( ionZ == 15 ) {
1082	G4double L_15[31]={4.527,5.052,6.042,6.982,7.882,8.744,9.566,10.34,11.08,14.08,16.15,17.55,18.48,19.53,20.01,20.2,20.23,20.17,20.04,19.88,18.82,17.69,16.62,15.65,13.97,12.59,11.45,10.49,9.67,8.968,8.359};
1083	physicsVector = CreatePhysicsVector(E,L_15,factor);
1084	}
1085	else if( ionZ == 16 ) {
1086	G4double L_16[31]={4.659,5.201,6.216,7.177,8.1,8.985,9.833,10.64,11.41,14.61,16.87,18.43,19.5,20.75,21.36,21.64,21.74,21.72,21.64,21.5,20.5,19.38,18.29,17.29,15.53,14.06,12.84,11.79,10.9,10.13,9.462};
1087	physicsVector = CreatePhysicsVector(E,L_16,factor);
1088	}
1089	else if( ionZ == 17 ) {
1090	G4double L_17[31]={4.821,5.374,6.409,7.386,8.323,9.227,10.1,10.93,11.72,15.09,17.52,19.25,20.46,21.9,22.64,23.01,23.17,23.2,23.15,23.04,22.11,20.99,19.89,18.86,17.03,15.49,14.19,13.07,12.12,11.28,10.56};
1091	physicsVector = CreatePhysicsVector(E,L_17,factor);
1092	}
1093	else if( ionZ == 18 ) {
1094	G4double L_18[31]={4.923,5.487,6.539,7.527,8.475,9.391,10.27,11.12,11.94,15.45,18.04,19.92,21.26,22.89,23.75,24.21,24.43,24.51,24.49,24.4,23.52,22.41,21.3,20.24,18.36,16.77,15.4,14.24,13.23,12.35,11.58};
1095	physicsVector = CreatePhysicsVector(E,L_18,factor);
1096	}
1097	}
1098
1099	else if( matName == namesMat[13] ) {
1100	if( ionZ == 3 ) {
1101	G4double M_3[31]={3.507,3.775,4.211,4.539,4.78,4.952,5.072,5.152,5.203,5.217,5.066,4.872,4.671,4.291,3.955,3.664,3.41,3.188,2.992,2.819,2.187,1.79,1.517,1.32,1.052,0.8799,0.759,0.6694,0.6001,0.5449,0.4998};
1102	physicsVector = CreatePhysicsVector(E,M_3,factor);
1103	}
1104	else if( ionZ == 4 ) {
1105	G4double M_4[31]={4.48,4.839,5.448,5.94,6.331,6.635,6.866,7.04,7.169,7.41,7.349,7.185,6.985,6.562,6.159,5.791,5.459,5.16,4.891,4.647,3.713,3.091,2.648,2.319,1.862,1.561,1.348,1.189,1.066,0.9678,0.8876};
1106	physicsVector = CreatePhysicsVector(E,M_4,factor);
1107	}
1108	else if( ionZ == 5 ) {
1109	G4double M_5[31]={5.276,5.728,6.503,7.151,7.691,8.133,8.489,8.772,8.995,9.539,9.619,9.522,9.352,8.935,8.501,8.087,7.7,7.342,7.013,6.71,5.503,4.657,4.034,3.559,2.884,2.428,2.102,1.856,1.665,1.512,1.386};
1110	physicsVector = CreatePhysicsVector(E,M_5,factor);
1111	}
1112	else if( ionZ == 6 ) {
1113	G4double M_6[31]={5.909,6.449,7.385,8.181,8.863,9.444,9.93,10.33,10.66,11.57,11.83,11.83,11.71,11.34,10.9,10.46,10.04,9.64,9.265,8.915,7.474,6.421,5.624,5,4.094,3.468,3.012,2.666,2.394,2.175,1.995};
1114	physicsVector = CreatePhysicsVector(E,M_6,factor);
1115	}
1116	else if( ionZ == 7 ) {
1117	G4double M_7[31]={6.512,7.125,8.193,9.111,9.914,10.61,11.22,11.74,12.17,13.47,13.96,14.09,14.05,13.74,13.33,12.89,12.44,12.02,11.61,11.23,9.59,8.349,7.384,6.615,5.473,4.667,4.07,3.612,3.249,2.956,2.713};
1118	physicsVector = CreatePhysicsVector(E,M_7,factor);
1119	}
1120	else if( ionZ == 8 ) {
1121	G4double M_8[31]={7.02,7.704,8.902,9.939,10.85,11.67,12.39,13.02,13.56,15.28,16.02,16.3,16.36,16.16,15.79,15.37,14.92,14.48,14.05,13.64,11.84,10.43,9.31,8.399,7.017,6.023,5.275,4.695,4.232,3.855,3.542};
1122	physicsVector = CreatePhysicsVector(E,M_8,factor);
1123	}
1124	else if( ionZ == 9 ) {
1125	G4double M_9[31]={7.39,8.136,9.451,10.59,11.61,12.53,13.35,14.08,14.72,16.86,17.86,18.3,18.47,18.39,18.07,17.67,17.23,16.79,16.35,15.92,14.01,12.47,11.22,10.19,8.599,7.436,6.549,5.853,5.292,4.832,4.448};
1126	physicsVector = CreatePhysicsVector(E,M_9,factor);
1127	}
1128	else if( ionZ == 10 ) {
1129	G4double M_10[31]={7.713,8.512,9.933,11.17,12.28,13.29,14.2,15.03,15.76,18.32,19.61,20.23,20.51,20.56,20.32,19.95,19.53,19.09,18.64,18.2,16.2,14.54,13.18,12.04,10.26,8.933,7.909,7.098,6.439,5.894,5.436};
1130	physicsVector = CreatePhysicsVector(E,M_10,factor);
1131	}
1132	else if( ionZ == 11 ) {
1133	G4double M_11[31]={7.985,8.827,10.34,11.67,12.87,13.95,14.95,15.87,16.7,19.7,21.3,22.14,22.57,22.8,22.65,22.34,21.95,21.52,21.08,20.64,18.56,16.8,15.32,14.07,12.09,10.59,9.417,8.48,7.713,7.076,6.537};
1134	physicsVector = CreatePhysicsVector(E,M_11,factor);
1135	}
1136	else if( ionZ == 12 ) {
1137	G4double M_12[31]={8.334,9.2,10.77,12.16,13.42,14.58,15.64,16.62,17.53,20.92,22.84,23.9,24.48,24.91,24.87,24.62,24.26,23.86,23.42,22.98,20.86,19,17.43,16.08,13.93,12.27,10.97,9.911,9.042,8.315,7.698};
1138	physicsVector = CreatePhysicsVector(E,M_12,factor);
1139	}
1140	else if( ionZ == 13 ) {
1141	G4double M_13[31]={8.64,9.535,11.16,12.62,13.95,15.16,16.28,17.33,18.3,22.05,24.28,25.57,26.32,26.96,27.05,26.87,26.57,26.19,25.78,25.35,23.19,21.27,19.61,18.17,15.84,14.04,12.6,11.43,10.46,9.637,8.938};
1142	physicsVector = CreatePhysicsVector(E,M_13,factor);
1143	}
1144	else if( ionZ == 14 ) {
1145	G4double M_14[31]={8.94,9.859,11.54,13.06,14.44,15.7,16.88,17.97,19,23.07,25.6,27.12,28.03,28.89,29.11,29.02,28.76,28.42,28.03,27.61,25.45,23.47,21.73,20.22,17.75,15.8,14.24,12.96,11.89,10.99,10.21};
1146	physicsVector = CreatePhysicsVector(E,M_14,factor);
1147	}
1148	else if( ionZ == 15 ) {
1149	G4double M_15[31]={9.314,10.26,11.99,13.56,15,16.32,17.54,18.69,19.76,24.14,26.97,28.71,29.8,30.9,31.27,31.27,31.09,30.8,30.44,30.05,27.91,25.88,24.08,22.49,19.85,17.76,16.07,14.66,13.49,12.49,11.63};
1150	physicsVector = CreatePhysicsVector(E,M_15,factor);
1151	}
1152	else if( ionZ == 16 ) {
1153	G4double M_16[31]={9.598,10.57,12.37,14,15.49,16.87,18.15,19.34,20.47,25.14,28.27,30.26,31.54,32.89,33.42,33.54,33.43,33.19,32.88,32.51,30.42,28.36,26.49,24.84,22.05,19.81,17.98,16.46,15.17,14.07,13.12};
1154	physicsVector = CreatePhysicsVector(E,M_16,factor);
1155	}
1156	else if( ionZ == 17 ) {
1157	G4double M_17[31]={9.983,10.97,12.8,14.46,16,17.42,18.74,19.98,21.15,26.06,29.46,31.68,33.13,34.73,35.42,35.64,35.61,35.43,35.15,34.81,32.77,30.69,28.78,27.06,24.15,21.79,19.84,18.21,16.83,15.64,14.61};
1158	physicsVector = CreatePhysicsVector(E,M_17,factor);
1159	}
1160	else if( ionZ == 18 ) {
1161	G4double M_18[31]={10.18,11.19,13.05,14.75,16.32,17.78,19.14,20.42,21.63,26.77,30.43,32.88,34.52,36.37,37.21,37.55,37.59,37.46,37.22,36.91,34.93,32.84,30.89,29.12,26.11,23.64,21.6,19.88,18.42,17.15,16.05};
1162	physicsVector = CreatePhysicsVector(E,M_18,factor);
1163	}
1164	}
1165
1166	else if( matName == namesMat[14] ) {
1167	if( ionZ == 3 ) {
1168	G4double N_3[31]={2.654,2.875,3.232,3.499,3.697,3.84,3.942,4.013,4.059,4.086,3.975,3.828,3.676,3.388,3.134,2.912,2.719,2.548,2.398,2.264,1.77,1.455,1.238,1.08,0.8645,0.7247,0.6264,0.5534,0.4968,0.4516,0.4146};
1169	physicsVector = CreatePhysicsVector(E,N_3,factor);
1170	}
1171	else if( ionZ == 4 ) {
1172	G4double N_4[31]={3.29,3.594,4.103,4.51,4.833,5.088,5.286,5.438,5.552,5.786,5.758,5.641,5.493,5.179,4.878,4.602,4.351,4.124,3.918,3.731,3.007,2.517,2.165,1.901,1.532,1.287,1.114,0.984,0.8833,0.8027,0.7368};
1173	physicsVector = CreatePhysicsVector(E,N_4,factor);
1174	}
1175	else if( ionZ == 5 ) {
1176	G4double N_5[31]={3.785,4.166,4.817,5.355,5.801,6.169,6.471,6.715,6.912,7.417,7.519,7.465,7.347,7.046,6.727,6.42,6.131,5.863,5.614,5.384,4.456,3.794,3.301,2.922,2.377,2.007,1.74,1.539,1.381,1.255,1.152};
1177	physicsVector = CreatePhysicsVector(E,N_5,factor);
1178	}
1179	else if( ionZ == 6 ) {
1180	G4double N_6[31]={4.166,4.615,5.399,6.062,6.626,7.106,7.512,7.852,8.137,8.959,9.23,9.267,9.2,8.94,8.626,8.306,7.995,7.699,7.419,7.155,6.057,5.238,4.609,4.112,3.382,2.873,2.499,2.214,1.99,1.81,1.661};
1181	physicsVector = CreatePhysicsVector(E,N_6,factor);
1182	}
1183	else if( ionZ == 7 ) {
1184	G4double N_7[31]={4.507,5.009,5.902,6.671,7.339,7.919,8.421,8.854,9.225,10.39,10.86,11.02,11.02,10.84,10.55,10.23,9.915,9.603,9.302,9.015,7.779,6.82,6.063,5.452,4.532,3.876,3.386,3.008,2.708,2.465,2.264};
1185	physicsVector = CreatePhysicsVector(E,N_7,factor);
1186	}
1187	else if( ionZ == 8 ) {
1188	G4double N_8[31]={4.804,5.352,6.342,7.21,7.973,8.648,9.242,9.765,10.22,11.74,12.44,12.73,12.83,12.74,12.5,12.2,11.89,11.57,11.26,10.95,9.606,8.524,7.649,6.928,5.82,5.012,4.399,3.92,3.536,3.223,2.962};
1189	physicsVector = CreatePhysicsVector(E,N_8,factor);
1190	}
1191	else if( ionZ == 9 ) {
1192	G4double N_9[31]={5.032,5.618,6.689,7.641,8.489,9.247,9.924,10.53,11.06,12.92,13.86,14.3,14.49,14.51,14.33,14.05,13.75,13.43,13.11,12.8,11.38,10.2,9.227,8.415,7.144,6.199,5.471,4.896,4.431,4.048,3.727};
1193	physicsVector = CreatePhysicsVector(E,N_9,factor);
1194	}
1195	else if( ionZ == 10 ) {
1196	G4double N_10[31]={5.239,5.856,6.994,8.019,8.942,9.776,10.53,11.21,11.82,14.01,15.19,15.8,16.1,16.25,16.12,15.89,15.6,15.29,14.97,14.65,13.17,11.91,10.85,9.956,8.532,7.456,6.616,5.946,5.399,4.945,4.562};
1197	physicsVector = CreatePhysicsVector(E,N_10,factor);
1198	}
1199	else if( ionZ == 11 ) {
1200	G4double N_11[31]={5.425,6.072,7.277,8.376,9.379,10.29,11.13,11.89,12.58,15.17,16.66,17.51,17.97,18.33,18.33,18.18,17.95,17.67,17.38,17.07,15.58,14.24,13.08,12.07,10.44,9.172,8.168,7.356,6.687,6.127,5.653};
1201	physicsVector = CreatePhysicsVector(E,N_11,factor);
1202	}
1203	else if( ionZ == 12 ) {
1204	G4double N_12[31]={5.628,6.282,7.507,8.639,9.68,10.64,11.52,12.33,13.07,15.91,17.62,18.62,19.19,19.67,19.72,19.59,19.36,19.09,18.79,18.47,16.92,15.53,14.32,13.28,11.57,10.24,9.174,8.308,7.589,6.985,6.471};
1205	physicsVector = CreatePhysicsVector(E,N_12,factor);
1206	}
1207	else if( ionZ == 13 ) {
1208	G4double N_13[31]={5.824,6.493,7.751,8.919,10,11.01,11.94,12.81,13.61,16.73,18.7,19.89,20.61,21.28,21.45,21.38,21.2,20.95,20.67,20.36,18.8,17.37,16.1,14.99,13.15,11.71,10.54,9.579,8.777,8.098,7.517};
1209	physicsVector = CreatePhysicsVector(E,N_13,factor);
1210	}
1211	else if( ionZ == 14 ) {
1212	G4double N_14[31]={6.02,6.703,7.987,9.187,10.31,11.36,12.34,13.25,14.09,17.47,19.68,21.07,21.94,22.81,23.09,23.1,22.96,22.75,22.48,22.19,20.64,19.16,17.85,16.68,14.73,13.18,11.92,10.87,9.99,9.241,8.596};
1213	physicsVector = CreatePhysicsVector(E,N_14,factor);
1214	}
1215	else if( ionZ == 15 ) {
1216	G4double N_15[31]={6.246,6.943,8.254,9.482,10.64,11.72,12.74,13.7,14.59,18.21,20.66,22.25,23.28,24.37,24.79,24.89,24.82,24.64,24.41,24.15,22.63,21.13,19.77,18.56,16.5,14.83,13.46,12.31,11.35,10.52,9.8};
1217	physicsVector = CreatePhysicsVector(E,N_15,factor);
1218	}
1219	else if( ionZ == 16 ) {
1220	G4double N_16[31]={6.433,7.148,8.491,9.751,10.94,12.07,13.12,14.12,15.05,18.91,21.59,23.4,24.59,25.91,26.47,26.67,26.67,26.55,26.35,26.11,24.65,23.15,21.75,20.49,18.32,16.55,15.07,13.83,12.78,11.87,11.08};
1221	physicsVector = CreatePhysicsVector(E,N_16,factor);
1222	}
1223	else if( ionZ == 17 ) {
1224	G4double N_17[31]={6.668,7.395,8.757,10.04,11.25,12.4,13.49,14.52,15.49,19.55,22.45,24.45,25.81,27.35,28.06,28.36,28.42,28.35,28.19,27.98,26.57,25.07,23.64,22.34,20.08,18.21,16.65,15.32,14.19,13.21,12.35};
1225	physicsVector = CreatePhysicsVector(E,N_17,factor);
1226	}
1227	else if( ionZ == 18 ) {
1228	G4double N_18[31]={6.799,7.541,8.924,10.22,11.45,12.62,13.74,14.8,15.8,20.04,23.14,25.33,26.84,28.61,29.46,29.86,29.99,29.96,29.84,29.65,28.3,26.79,25.35,24.01,21.68,19.74,18.1,16.71,15.51,14.47,13.56};
1229	physicsVector = CreatePhysicsVector(E,N_18,factor);
1230	}
1231	}
1232
1233	else if( matName == namesMat[15] ) {
1234	if( ionZ == 3 ) {
1235	G4double O_3[31]={2.784,3.026,3.417,3.703,3.906,4.047,4.141,4.202,4.238,4.231,4.103,3.948,3.79,3.492,3.23,3.001,2.801,2.625,2.469,2.331,1.821,1.496,1.272,1.108,0.8863,0.7423,0.6413,0.5662,0.5082,0.4618,0.4239};
1236	physicsVector = CreatePhysicsVector(E,O_3,factor);
1237	}
1238	else if( ionZ == 4 ) {
1239	G4double O_4[31]={3.441,3.767,4.322,4.765,5.11,5.374,5.57,5.714,5.818,6.004,5.951,5.822,5.666,5.341,5.031,4.746,4.487,4.252,4.039,3.845,3.096,2.589,2.225,1.952,1.571,1.319,1.14,1.007,0.9035,0.8208,0.7532};
1240	physicsVector = CreatePhysicsVector(E,O_4,factor);
1241	}
1242	else if( ionZ == 5 ) {
1243	G4double O_5[31]={3.957,4.361,5.063,5.649,6.133,6.524,6.834,7.078,7.266,7.715,7.782,7.711,7.583,7.27,6.941,6.625,6.327,6.05,5.792,5.554,4.592,3.905,3.395,3.001,2.438,2.057,1.782,1.575,1.413,1.284,1.178};
1244	physicsVector = CreatePhysicsVector(E,O_5,factor);
1245	}
1246	else if( ionZ == 6 ) {
1247	G4double O_6[31]={4.362,4.834,5.668,6.385,6.997,7.514,7.942,8.292,8.576,9.344,9.568,9.581,9.501,9.227,8.903,8.573,8.253,7.947,7.657,7.384,6.245,5.394,4.741,4.226,3.47,2.944,2.559,2.266,2.036,1.85,1.698};
1248	physicsVector = CreatePhysicsVector(E,O_6,factor);
1249	}
1250	else if( ionZ == 7 ) {
1251	G4double O_7[31]={4.74,5.265,6.205,7.026,7.747,8.374,8.912,9.367,9.748,10.87,11.29,11.41,11.4,11.19,10.89,10.57,10.24,9.915,9.604,9.307,8.022,7.024,6.237,5.602,4.65,3.971,3.466,3.077,2.769,2.519,2.313};
1252	physicsVector = CreatePhysicsVector(E,O_7,factor);
1253	}
1254	else if( ionZ == 8 ) {
1255	G4double O_8[31]={5.063,5.635,6.671,7.587,8.405,9.133,9.773,10.33,10.81,12.3,12.94,13.2,13.27,13.16,12.91,12.6,12.27,11.95,11.62,11.31,9.911,8.784,7.871,7.122,5.972,5.136,4.503,4.009,3.615,3.293,3.025};
1256	physicsVector = CreatePhysicsVector(E,O_8,factor);
1257	}
1258	else if( ionZ == 9 ) {
1259	G4double O_9[31]={7.121,7.892,9.331,10.67,11.93,13.12,14.24,15.31,16.32,20.61,23.59,25.53,26.82,28.26,28.93,29.21,29.27,29.19,29.03,28.81,27.36,25.79,24.3,22.94,20.59,18.64,17.02,15.65,14.48,13.47,12.59};
1260	physicsVector = CreatePhysicsVector(E,O_9,factor);
1261	}
1262	else if( ionZ == 10 ) {
1263	G4double O_10[31]={5.543,6.188,7.373,8.442,9.413,10.3,11.11,11.84,12.49,14.73,15.84,16.4,16.66,16.77,16.63,16.38,16.08,15.76,15.44,15.11,13.56,12.25,11.15,10.22,8.743,7.63,6.764,6.074,5.512,5.046,4.655};
1264	physicsVector = CreatePhysicsVector(E,O_10,factor);
1265	}
1266	else if( ionZ == 11 ) {
1267	G4double O_11[31]={5.732,6.403,7.645,8.776,9.811,10.76,11.64,12.45,13.18,15.8,17.2,17.93,18.32,18.59,18.53,18.33,18.06,17.76,17.44,17.11,15.53,14.14,12.95,11.93,10.3,9.044,8.057,7.261,6.609,6.064,5.603};
1268	physicsVector = CreatePhysicsVector(E,O_11,factor);
1269	}
1270	else if( ionZ == 12 ) {
1271	G4double O_12[31]={5.969,6.654,7.93,9.104,10.19,11.19,12.12,12.99,13.78,16.75,18.42,19.34,19.87,20.29,20.33,20.19,19.95,19.67,19.36,19.04,17.42,15.98,14.72,13.63,11.86,10.48,9.379,8.487,7.749,7.129,6.601};
1272	physicsVector = CreatePhysicsVector(E,O_12,factor);
1273	}
1274	else if( ionZ == 13 ) {
1275	G4double O_13[31]={6.184,6.887,8.201,9.414,10.54,11.59,12.57,13.48,14.33,17.62,19.56,20.69,21.35,21.96,22.1,22.03,21.84,21.58,21.29,20.98,19.36,17.87,16.55,15.39,13.48,11.98,10.78,9.788,8.963,8.266,7.67};
1276	physicsVector = CreatePhysicsVector(E,O_13,factor);
1277	}
1278	else if( ionZ == 14 ) {
1279	G4double O_14[31]={6.399,7.118,8.463,9.71,10.87,11.96,12.98,13.94,14.83,18.4,20.6,21.93,22.73,23.53,23.79,23.79,23.64,23.42,23.15,22.85,21.24,19.71,18.33,17.12,15.1,13.49,12.18,11.1,10.2,9.43,8.769};
1280	physicsVector = CreatePhysicsVector(E,O_14,factor);
1281	}
1282	else if( ionZ == 15 ) {
1283	G4double O_15[31]={6.653,7.389,8.766,10.05,11.24,12.37,13.42,14.42,15.36,19.19,21.66,23.19,24.15,25.16,25.55,25.63,25.56,25.38,25.14,24.87,23.3,21.74,20.32,19.05,16.91,15.18,13.76,12.58,11.58,10.73,9.995};
1284	physicsVector = CreatePhysicsVector(E,O_15,factor);
1285	}
1286	else if( ionZ == 16 ) {
1287	G4double O_16[31]={6.858,7.616,9.031,10.35,11.58,12.74,13.84,14.87,15.85,19.93,22.66,24.41,25.53,26.76,27.29,27.48,27.47,27.34,27.15,26.9,25.39,23.82,22.36,21.05,18.79,16.94,15.42,14.14,13.05,12.11,11.3};
1288	physicsVector = CreatePhysicsVector(E,O_16,factor);
1289	}
1290	else if( ionZ == 17 ) {
1291	G4double O_17[31]={7.121,7.892,9.331,10.67,11.93,13.12,14.24,15.31,16.32,20.61,23.59,25.53,26.82,28.26,28.93,29.21,29.27,29.19,29.03,28.81,27.36,25.79,24.3,22.94,20.59,18.64,17.02,15.65,14.48,13.47,12.59};
1292	physicsVector = CreatePhysicsVector(E,O_17,factor);
1293	}
1294	else if( ionZ == 18 ) {
1295	G4double O_18[31]={7.278,8.066,9.53,10.89,12.17,13.38,14.53,15.62,16.66,21.14,24.34,26.49,27.93,29.6,30.41,30.78,30.91,30.88,30.76,30.57,29.17,27.6,26.09,24.69,22.25,20.22,18.52,17.08,15.84,14.77,13.83};
1296	physicsVector = CreatePhysicsVector(E,O_18,factor);
1297	}
1298	}
1299
1300	else if( matName == namesMat[16] ) {
1301	if( ionZ == 3 ) {
1302	G4double P_3[31]={1.023,1.127,1.298,1.427,1.522,1.591,1.641,1.677,1.704,1.755,1.743,1.711,1.671,1.585,1.5,1.42,1.346,1.279,1.218,1.162,0.9451,0.7987,0.6937,0.6149,0.5043,0.4302,0.3768,0.3364,0.3046,0.2788,0.2575};
1303	physicsVector = CreatePhysicsVector(E,P_3,factor);
1304	}
1305	else if( ionZ == 4 ) {
1306	G4double P_4[31]={1.245,1.386,1.628,1.825,1.983,2.106,2.202,2.276,2.335,2.487,2.524,2.517,2.489,2.41,2.319,2.227,2.138,2.053,1.972,1.897,1.589,1.367,1.201,1.073,0.8873,0.7599,0.667,0.596,0.54,0.4945,0.4569};
1307	physicsVector = CreatePhysicsVector(E,P_4,factor);
1308	}
1309	else if( ionZ == 5 ) {
1310	G4double P_5[31]={1.418,1.587,1.888,2.146,2.363,2.543,2.689,2.808,2.904,3.184,3.29,3.323,3.32,3.266,3.184,3.091,2.996,2.9,2.808,2.719,2.336,2.043,1.816,1.635,1.366,1.177,1.036,0.9275,0.8412,0.7709,0.7124};
1311	physicsVector = CreatePhysicsVector(E,P_5,factor);
1312	}
1313	else if( ionZ == 6 ) {
1314	G4double P_6[31]={1.557,1.747,2.097,2.406,2.676,2.91,3.107,3.273,3.411,3.837,4.029,4.114,4.146,4.132,4.07,3.986,3.892,3.794,3.695,3.598,3.158,2.803,2.517,2.285,1.93,1.674,1.48,1.328,1.207,1.107,1.024};
1315	physicsVector = CreatePhysicsVector(E,P_6,factor);
1316	}
1317	else if( ionZ == 7 ) {
1318	G4double P_7[31]={1.678,1.883,2.268,2.618,2.933,3.213,3.458,3.669,3.85,4.436,4.726,4.876,4.951,4.991,4.96,4.894,4.809,4.715,4.616,4.515,4.036,3.63,3.292,3.011,2.572,2.246,1.995,1.796,1.635,1.502,1.391};
1319	physicsVector = CreatePhysicsVector(E,P_7,factor);
1320	}
1321	else if( ionZ == 8 ) {
1322	G4double P_8[31]={1.791,2.007,2.419,2.803,3.155,3.476,3.764,4.019,4.242,4.994,5.393,5.614,5.74,5.845,5.854,5.813,5.744,5.658,5.563,5.463,4.961,4.512,4.129,3.802,3.281,2.886,2.576,2.328,2.125,1.955,1.812};
1323	physicsVector = CreatePhysicsVector(E,P_8,factor);
1324	}
1325	else if( ionZ == 9 ) {
1326	G4double P_9[31]={1.889,2.115,2.548,2.958,3.34,3.694,4.02,4.313,4.575,5.493,6.004,6.304,6.485,6.663,6.716,6.704,6.654,6.58,6.493,6.397,5.885,5.404,4.984,4.62,4.026,3.567,3.201,2.905,2.659,2.453,2.278};
1327	physicsVector = CreatePhysicsVector(E,P_9,factor);
1328	}
1329	else if( ionZ == 10 ) {
1330	G4double P_10[31]={1.981,2.217,2.667,3.097,3.504,3.885,4.241,4.569,4.867,5.948,6.577,6.96,7.202,7.46,7.564,7.585,7.559,7.501,7.424,7.335,6.824,6.321,5.868,5.47,4.811,4.29,3.871,3.527,3.239,2.996,2.788};
1331	physicsVector = CreatePhysicsVector(E,P_10,factor);
1332	}
1333	else if( ionZ == 11 ) {
1334	G4double P_11[31]={2.062,2.309,2.776,3.225,3.654,4.061,4.446,4.806,5.138,6.389,7.15,7.631,7.948,8.31,8.483,8.553,8.561,8.529,8.471,8.395,7.906,7.387,6.906,6.472,5.739,5.149,4.665,4.264,3.925,3.637,3.388};
1335	physicsVector = CreatePhysicsVector(E,P_11,factor);
1336	}
1337	else if( ionZ == 12 ) {
1338	G4double P_12[31]={2.147,2.4,2.877,3.336,3.778,4.201,4.604,4.986,5.344,6.735,7.608,8.174,8.553,9,9.226,9.334,9.368,9.354,9.308,9.241,8.762,8.231,7.73,7.275,6.498,5.865,5.342,4.905,4.532,4.213,3.937};
1339	physicsVector = CreatePhysicsVector(E,P_12,factor);
1340	}
1341	else if( ionZ == 13 ) {
1342	G4double P_13[31]={2.227,2.489,2.98,3.448,3.903,4.341,4.761,5.163,5.544,7.074,8.069,8.729,9.181,9.73,10.02,10.18,10.25,10.26,10.24,10.19,9.739,9.205,8.686,8.209,7.381,6.698,6.128,5.645,5.233,4.877,4.566};
1343	physicsVector = CreatePhysicsVector(E,P_13,factor);
1344	}
1345	else if( ionZ == 14 ) {
1346	G4double P_14[31]={2.303,2.574,3.079,3.557,4.021,4.47,4.905,5.322,5.722,7.378,8.491,9.246,9.773,10.43,10.79,11,11.11,11.15,11.15,11.11,10.71,10.17,9.643,9.147,8.275,7.546,6.931,6.407,5.957,5.565,5.222};
1347	physicsVector = CreatePhysicsVector(E,P_14,factor);
1348	}
1349	else if( ionZ == 15 ) {
1350	G4double P_15[31]={2.379,2.659,3.178,3.665,4.137,4.597,5.043,5.474,5.89,7.66,8.89,9.743,10.35,11.12,11.56,11.83,11.98,12.06,12.08,12.07,11.72,11.2,10.66,10.15,9.236,8.462,7.802,7.235,6.743,6.314,5.935};
1351	physicsVector = CreatePhysicsVector(E,P_15,factor);
1352	}
1353	else if( ionZ == 16 ) {
1354	G4double P_16[31]={1.791,2.007,2.419,2.803,3.155,3.476,3.764,4.019,4.242,4.994,5.393,5.614,5.74,5.845,5.854,5.813,5.744,5.658,5.563,5.463,4.961,4.512,4.129,3.802,3.281,2.886,2.576,2.328,2.125,1.955,1.812};
1355	physicsVector = CreatePhysicsVector(E,P_16,factor);
1356	}
1357	else if( ionZ == 17 ) {
1358	G4double P_17[31]={2.525,2.82,3.37,3.879,4.368,4.844,5.308,5.761,6.202,8.163,9.615,10.66,11.42,12.43,13.04,13.43,13.68,13.83,13.92,13.95,13.74,13.26,12.72,12.19,11.21,10.35,9.614,8.967,8.4,7.898,7.452};
1359	physicsVector = CreatePhysicsVector(E,P_17,factor);
1360	}
1361	else if( ionZ == 18 ) {
1362	G4double P_18[31]={2.583,2.886,3.449,3.97,4.466,4.948,5.419,5.879,6.328,8.364,9.915,11.05,11.89,13.02,13.71,14.16,14.46,14.65,14.76,14.82,14.67,14.22,13.68,13.15,12.14,11.26,10.48,9.805,9.206,8.674,8.2};
1363	physicsVector = CreatePhysicsVector(E,P_18,factor);
1364	}
1365	}
1366
1367	else if( matName == namesMat[17] ) {
1368	if( ionZ == 3 ) {
1369	G4double Q_3[31]={2.597,2.831,3.212,3.493,3.691,3.828,3.919,3.977,4.012,4.006,3.89,3.747,3.601,3.326,3.081,2.867,2.679,2.513,2.366,2.235,1.75,1.439,1.224,1.068,0.8543,0.7158,0.6185,0.5463,0.4904,0.4457,0.4092};
1370	physicsVector = CreatePhysicsVector(E,Q_3,factor);
1371	}
1372	else if( ionZ == 4 ) {
1373	G4double Q_4[31]={3.198,3.511,4.048,4.481,4.821,5.078,5.269,5.409,5.509,5.688,5.643,5.528,5.387,5.088,4.801,4.536,4.294,4.073,3.872,3.689,2.977,2.493,2.144,1.881,1.515,1.273,1.1,0.9718,0.8721,0.7924,0.7272};
1374	physicsVector = CreatePhysicsVector(E,Q_4,factor);
1375	}
1376	else if( ionZ == 5 ) {
1377	G4double Q_5[31]={3.674,4.056,4.73,5.3,5.774,6.157,6.461,6.698,6.881,7.313,7.382,7.323,7.21,6.926,6.626,6.333,6.056,5.797,5.555,5.331,4.419,3.763,3.273,2.896,2.354,1.985,1.72,1.52,1.365,1.24,1.137};
1378	physicsVector = CreatePhysicsVector(E,Q_5,factor);
1379	}
1380	else if( ionZ == 6 ) {
1381	G4double Q_6[31]={4.051,4.495,5.289,5.98,6.577,7.083,7.502,7.845,8.121,8.862,9.081,9.103,9.036,8.794,8.5,8.198,7.901,7.617,7.346,7.09,6.012,5.2,4.574,4.079,3.351,2.844,2.472,2.189,1.967,1.787,1.64};
1382	physicsVector = CreatePhysicsVector(E,Q_6,factor);
1383	}
1384	else if( ionZ == 7 ) {
1385	G4double Q_7[31]={4.399,4.891,5.779,6.565,7.264,7.876,8.403,8.849,9.221,10.31,10.71,10.84,10.84,10.67,10.4,10.11,9.804,9.506,9.216,8.938,7.725,6.775,6.021,5.411,4.493,3.838,3.351,2.975,2.677,2.435,2.235};
1386	physicsVector = CreatePhysicsVector(E,Q_7,factor);
1387	}
1388	else if( ionZ == 8 ) {
1389	G4double Q_8[31]={4.708,5.245,6.219,7.091,7.878,8.585,9.21,9.755,10.22,11.68,12.29,12.55,12.63,12.55,12.33,12.05,11.76,11.46,11.16,10.87,9.548,8.476,7.603,6.883,5.775,4.967,4.355,3.878,3.496,3.184,2.926};
1390	physicsVector = CreatePhysicsVector(E,Q_8,factor);
1391	}
1392	else if( ionZ == 9 ) {
1393	G4double Q_9[31]={4.95,5.523,6.572,7.516,8.376,9.161,9.871,10.5,11.06,12.88,13.72,14.1,14.27,14.29,14.12,13.87,13.59,13.29,12.99,12.69,11.3,10.13,9.163,8.353,7.083,6.139,5.414,4.841,4.378,3.998,3.68};
1394	physicsVector = CreatePhysicsVector(E,Q_9,factor);
1395	}
1396	else if( ionZ == 10 ) {
1397	G4double Q_10[31]={5.171,5.776,6.888,7.895,8.819,9.67,10.45,11.16,11.8,13.98,15.06,15.6,15.86,15.99,15.88,15.67,15.4,15.11,14.81,14.51,13.06,11.82,10.76,9.872,8.452,7.378,6.542,5.875,5.331,4.881,4.502};
1398	physicsVector = CreatePhysicsVector(E,Q_10,factor);
1399	}
1400	else if( ionZ == 11 ) {
1401	G4double Q_11[31]={5.362,5.997,7.176,8.252,9.245,10.17,11.03,11.82,12.54,15.16,16.55,17.3,17.72,18.06,18.08,17.95,17.74,17.49,17.22,16.93,15.49,14.17,13.02,12.02,10.39,9.122,8.117,7.303,6.634,6.075,5.602};
1402	physicsVector = CreatePhysicsVector(E,Q_11,factor);
1403	}
1404	else if( ionZ == 12 ) {
1405	G4double Q_12[31]={5.572,6.215,7.413,8.517,9.541,10.5,11.39,12.23,13,15.9,17.51,18.4,18.91,19.35,19.41,19.3,19.1,18.85,18.57,18.27,16.77,15.4,14.21,13.17,11.46,10.13,9.076,8.213,7.499,6.899,6.388};
1406	physicsVector = CreatePhysicsVector(E,Q_12,factor);
1407	}
1408	else if( ionZ == 13 ) {
1409	G4double Q_13[31]={5.774,6.434,7.668,8.809,9.871,10.87,11.8,12.68,13.51,16.71,18.59,19.68,20.32,20.93,21.1,21.05,20.9,20.67,20.41,20.13,18.63,17.22,15.96,14.86,13.03,11.59,10.43,9.471,8.674,7.999,7.422};
1410	physicsVector = CreatePhysicsVector(E,Q_13,factor);
1411	}
1412	else if( ionZ == 14 ) {
1413	G4double Q_14[31]={5.976,6.651,7.915,9.089,10.18,11.21,12.18,13.1,13.97,17.45,19.59,20.86,21.64,22.43,22.7,22.73,22.62,22.43,22.19,21.92,20.43,18.99,17.69,16.53,14.59,13.05,11.79,10.74,9.87,9.126,8.486};
1414	physicsVector = CreatePhysicsVector(E,Q_14,factor);
1415	}
1416	else if( ionZ == 15 ) {
1417	G4double Q_15[31]={6.211,6.903,8.199,9.405,10.53,11.59,12.6,13.55,14.45,18.19,20.58,22.06,22.99,23.98,24.38,24.49,24.45,24.3,24.1,23.86,22.42,20.95,19.61,18.4,16.34,14.68,13.32,12.18,11.21,10.39,9.676};
1418	physicsVector = CreatePhysicsVector(E,Q_15,factor);
1419	}
1420	else if( ionZ == 16 ) {
1421	G4double Q_16[31]={6.402,7.115,8.448,9.69,10.85,11.95,12.98,13.97,14.91,18.88,21.54,23.22,24.31,25.51,26.05,26.26,26.28,26.19,26.02,25.8,24.43,22.96,21.58,20.33,18.17,16.39,14.92,13.69,12.63,11.73,10.94};
1422	physicsVector = CreatePhysicsVector(E,Q_16,factor);
1423	}
1424	else if( ionZ == 17 ) {
1425	G4double Q_17[31]={6.647,7.373,8.73,9.996,11.18,12.3,13.37,14.38,15.35,19.51,22.41,24.29,25.53,26.94,27.61,27.91,28,27.96,27.83,27.64,26.33,24.86,23.46,22.16,19.91,18.04,16.48,15.16,14.03,13.05,12.2};
1426	physicsVector = CreatePhysicsVector(E,Q_17,factor);
1427	}
1428	else if( ionZ == 18 ) {
1429	G4double Q_18[31]={6.781,7.522,8.903,10.19,11.4,12.54,13.62,14.65,15.64,19.98,23.1,25.17,26.56,28.17,28.98,29.37,29.52,29.52,29.42,29.26,28.01,26.54,25.12,23.79,21.47,19.53,17.9,16.51,15.32,14.28,13.38};
1430	physicsVector = CreatePhysicsVector(E,Q_18,factor);
1431	}
1432	}
1433
1434	else if( matName == namesMat[18] ) {
1435	if( ionZ == 3 ) {
1436	G4double R_3[31]={2.358,2.573,2.927,3.191,3.383,3.517,3.609,3.671,3.709,3.726,3.628,3.501,3.37,3.119,2.895,2.698,2.524,2.37,2.234,2.112,1.658,1.366,1.163,1.015,0.8137,0.6824,0.59,0.5214,0.4682,0.4258,0.391};
1437	physicsVector = CreatePhysicsVector(E,R_3,factor);
1438	}
1439	else if( ionZ == 4 ) {
1440	G4double R_4[31]={2.899,3.186,3.682,4.086,4.408,4.656,4.843,4.983,5.085,5.286,5.263,5.166,5.041,4.773,4.512,4.269,4.046,3.842,3.656,3.486,2.822,2.367,2.038,1.79,1.443,1.213,1.05,0.9278,0.8329,0.7571,0.695};
1441	physicsVector = CreatePhysicsVector(E,R_4,factor);
1442	}
1443	else if( ionZ == 5 ) {
1444	G4double R_5[31]={3.326,3.677,4.296,4.825,5.269,5.633,5.926,6.158,6.34,6.789,6.88,6.841,6.747,6.497,6.226,5.959,5.706,5.467,5.245,5.037,4.188,3.574,3.113,2.757,2.243,1.894,1.642,1.452,1.303,1.184,1.087};
1445	physicsVector = CreatePhysicsVector(E,R_5,factor);
1446	}
1447	else if( ionZ == 6 ) {
1448	G4double R_6[31]={3.666,4.073,4.802,5.44,5.996,6.471,6.871,7.201,7.471,8.219,8.46,8.503,8.456,8.249,7.988,7.715,7.445,7.185,6.936,6.7,5.699,4.94,4.352,3.886,3.196,2.715,2.361,2.092,1.88,1.709,1.568};
1449	physicsVector = CreatePhysicsVector(E,R_6,factor);
1450	}
1451	else if( ionZ == 7 ) {
1452	G4double R_7[31]={3.974,4.425,5.241,5.967,6.615,7.188,7.685,8.11,8.47,9.55,9.975,10.12,10.14,10,9.776,9.513,9.24,8.969,8.704,8.448,7.327,6.439,5.732,5.158,4.289,3.667,3.203,2.844,2.56,2.329,2.139};
1453	physicsVector = CreatePhysicsVector(E,R_7,factor);
1454	}
1455	else if( ionZ == 8 ) {
1456	G4double R_8[31]={4.249,4.74,5.635,6.439,7.168,7.827,8.414,8.93,9.377,10.8,11.43,11.71,11.81,11.77,11.58,11.34,11.08,10.8,10.53,10.26,9.055,8.057,7.239,6.562,5.515,4.748,4.166,3.711,3.346,3.048,2.801};
1457	physicsVector = CreatePhysicsVector(E,R_8,factor);
1458	}
1459	else if( ionZ == 9 ) {
1460	G4double R_9[31]={4.469,4.993,5.956,6.826,7.623,8.352,9.016,9.611,10.13,11.9,12.75,13.16,13.35,13.4,13.27,13.06,12.81,12.54,12.26,11.99,10.71,9.634,8.728,7.967,6.767,5.871,5.181,4.635,4.193,3.829,3.525};
1461	physicsVector = CreatePhysicsVector(E,R_9,factor);
1462	}
1463	else if( ionZ == 10 ) {
1464	G4double R_10[31]={4.672,5.224,6.244,7.173,8.027,8.818,9.547,10.21,10.81,12.92,14,14.55,14.84,15.01,14.93,14.75,14.52,14.26,13.99,13.72,12.39,11.23,10.25,9.418,8.077,7.058,6.263,5.627,5.108,4.677,4.314};
1465	physicsVector = CreatePhysicsVector(E,R_10,factor);
1466	}
1467	else if( ionZ == 11 ) {
1468	G4double R_11[31]={4.841,5.42,6.5,7.49,8.406,9.261,10.06,10.8,11.47,13.97,15.35,16.11,16.54,16.9,16.95,16.85,16.68,16.46,16.21,15.96,14.64,13.43,12.36,11.42,9.896,8.702,7.752,6.981,6.345,5.813,5.362};
1469	physicsVector = CreatePhysicsVector(E,R_11,factor);
1470	}
1471	else if( ionZ == 12 ) {
1472	G4double R_12[31]={5.029,5.617,6.716,7.735,8.682,9.568,10.4,11.18,11.9,14.66,16.25,17.16,17.68,18.15,18.25,18.17,18,17.78,17.53,17.27,15.9,14.64,13.52,12.55,10.95,9.694,8.688,7.867,7.186,6.613,6.125};
1473	physicsVector = CreatePhysicsVector(E,R_12,factor);
1474	}
1475	else if( ionZ == 13 ) {
1476	G4double R_13[31]={5.209,5.812,6.943,7.995,8.978,9.901,10.77,11.59,12.36,15.41,17.25,18.34,18.99,19.64,19.83,19.82,19.7,19.51,19.28,19.02,17.66,16.36,15.19,14.16,12.45,11.08,9.982,9.073,8.313,7.669,7.118};
1477	physicsVector = CreatePhysicsVector(E,R_13,factor);
1478	}
1479	else if( ionZ == 14 ) {
1480	G4double R_14[31]={5.387,6.004,7.163,8.245,9.259,10.21,11.11,11.97,12.78,16.07,18.16,19.43,20.22,21.04,21.34,21.4,21.32,21.16,20.96,20.72,19.37,18.04,16.83,15.76,13.93,12.47,11.28,10.29,9.46,8.75,8.139};
1481	physicsVector = CreatePhysicsVector(E,R_14,factor);
1482	}
1483	else if( ionZ == 15 ) {
1484	G4double R_15[31]={5.591,6.222,7.409,8.52,9.564,10.54,11.48,12.37,13.21,16.74,19.07,20.53,21.47,22.49,22.92,23.06,23.05,22.94,22.76,22.55,21.25,19.91,18.67,17.54,15.61,14.04,12.75,11.66,10.75,9.965,9.285};
1485	physicsVector = CreatePhysicsVector(E,R_15,factor);
1486	}
1487	else if( ionZ == 16 ) {
1488	G4double R_16[31]={5.759,6.409,7.629,8.771,9.848,10.86,11.82,12.74,13.62,17.37,19.93,21.6,22.69,23.91,24.48,24.72,24.77,24.71,24.57,24.39,23.16,21.82,20.55,19.38,17.35,15.68,14.29,13.12,12.11,11.25,10.5};
1489	physicsVector = CreatePhysicsVector(E,R_16,factor);
1490	}
1491	else if( ionZ == 17 ) {
1492	G4double R_17[31]={5.972,6.633,7.875,9.04,10.14,11.18,12.17,13.11,14.02,17.94,20.73,22.58,23.82,25.25,25.95,26.28,26.4,26.38,26.28,26.13,24.97,23.63,22.33,21.13,19.02,17.26,15.79,14.53,13.45,12.52,11.71};
1493	physicsVector = CreatePhysicsVector(E,R_17,factor);
1494	}
1495	else if( ionZ == 18 ) {
1496	G4double R_18[31]={6.093,6.768,8.031,9.213,10.33,11.39,12.39,13.36,14.28,18.36,21.35,23.39,24.78,26.41,27.24,27.65,27.83,27.86,27.8,27.67,26.56,25.23,23.92,22.68,20.52,18.69,17.14,15.83,14.69,13.71,12.84};
1497	physicsVector = CreatePhysicsVector(E,R_18,factor);
1498	}
1499	}
1500
1501	else if( matName == namesMat[19] ) {
1502	if( ionZ == 3 ) {
1503	G4double S_3[31]={3.038,3.3,3.719,4.021,4.232,4.375,4.469,4.527,4.56,4.531,4.384,4.211,4.037,3.712,3.428,3.181,2.966,2.777,2.61,2.462,1.918,1.574,1.337,1.164,0.9296,0.778,0.6717,0.5928,0.5318,0.4832,0.4433};
1504	physicsVector = CreatePhysicsVector(E,S_3,factor);
1505	}
1506	else if( ionZ == 4 ) {
1507	G4double S_4[31]={3.761,4.114,4.711,5.183,5.547,5.82,6.02,6.164,6.267,6.432,6.358,6.209,6.036,5.678,5.34,5.03,4.751,4.498,4.269,4.062,3.261,2.722,2.337,2.048,1.647,1.382,1.194,1.054,0.9453,0.8586,0.7877};
1508	physicsVector = CreatePhysicsVector(E,S_4,factor);
1509	}
1510	else if( ionZ == 5 ) {
1511	G4double S_5[31]={4.329,4.766,5.525,6.153,6.667,7.077,7.399,7.648,7.839,8.273,8.318,8.226,8.079,7.729,7.368,7.023,6.7,6.4,6.123,5.866,4.837,4.106,3.565,3.149,2.555,2.154,1.865,1.647,1.478,1.342,1.231};
1512	physicsVector = CreatePhysicsVector(E,S_5,factor);
1513	}
1514	else if( ionZ == 6 ) {
1515	G4double S_6[31]={4.774,5.286,6.188,6.959,7.613,8.16,8.609,8.972,9.264,10.03,10.23,10.22,10.12,9.809,9.45,9.088,8.738,8.406,8.093,7.799,6.576,5.669,4.976,4.431,3.634,3.081,2.677,2.369,2.128,1.934,1.774};
1516	physicsVector = CreatePhysicsVector(E,S_6,factor);
1517	}
1518	else if( ionZ == 7 ) {
1519	G4double S_7[31]={5.194,5.763,6.779,7.664,8.436,9.103,9.671,10.15,10.54,11.67,12.07,12.17,12.14,11.9,11.56,11.2,10.84,10.49,10.15,9.827,8.445,7.38,6.543,5.872,4.867,4.154,3.624,3.216,2.893,2.632,2.416};
1520	physicsVector = CreatePhysicsVector(E,S_7,factor);
1521	}
1522	else if( ionZ == 8 ) {
1523	G4double S_8[31]={5.553,6.174,7.294,8.281,9.159,9.936,10.61,11.2,11.7,13.23,13.85,14.09,14.14,13.99,13.7,13.35,13,12.64,12.28,11.94,10.43,9.228,8.257,7.463,6.249,5.369,4.705,4.188,3.775,3.438,3.158};
1524	physicsVector = CreatePhysicsVector(E,S_8,factor);
1525	}
1526	else if( ionZ == 9 ) {
1527	G4double S_9[31]={5.826,6.49,7.698,8.774,9.741,10.61,11.39,12.07,12.67,14.59,15.45,15.83,15.97,15.92,15.68,15.36,15.01,14.65,14.3,13.94,12.34,11.03,9.95,9.055,7.662,6.633,5.845,5.225,4.725,4.313,3.97};
1528	physicsVector = CreatePhysicsVector(E,S_9,factor);
1529	}
1530	else if( ionZ == 10 ) {
1531	G4double S_10[31]={6.073,6.774,8.058,9.21,10.25,11.21,12.07,12.84,13.53,15.85,16.96,17.5,17.75,17.82,17.64,17.35,17.02,16.66,16.3,15.94,14.27,12.87,11.69,10.7,9.143,7.971,7.062,6.339,5.751,5.264,4.855};
1532	physicsVector = CreatePhysicsVector(E,S_10,factor);
1533	}
1534	else if( ionZ == 11 ) {
1535	G4double S_11[31]={6.28,7.01,8.355,9.574,10.69,11.71,12.65,13.51,14.28,17.02,18.42,19.15,19.52,19.75,19.66,19.42,19.11,18.77,18.42,18.06,16.34,14.85,13.58,12.5,10.77,9.448,8.411,7.577,6.894,6.324,5.843};
1536	physicsVector = CreatePhysicsVector(E,S_11,factor);
1537	}
1538	else if( ionZ == 12 ) {
1539	G4double S_12[31]={6.544,7.288,8.67,9.936,11.1,12.18,13.17,14.1,14.94,18.04,19.74,20.67,21.17,21.57,21.57,21.39,21.12,20.8,20.45,20.1,18.35,16.78,15.44,14.28,12.4,10.95,9.793,8.857,8.083,7.434,6.883};
1540	physicsVector = CreatePhysicsVector(E,S_12,factor);
1541	}
1542	else if( ionZ == 13 ) {
1543	G4double S_13[31]={6.783,7.547,8.969,10.28,11.49,12.61,13.66,14.64,15.54,18.99,20.98,22.11,22.76,23.34,23.45,23.34,23.11,22.82,22.5,22.15,20.39,18.77,17.36,16.13,14.11,12.52,11.25,10.21,9.349,8.619,7.996};
1544	physicsVector = CreatePhysicsVector(E,S_13,factor);
1545	}
1546	else if( ionZ == 14 ) {
1547	G4double S_14[31]={7.021,7.803,9.26,10.61,11.86,13.02,14.11,15.13,16.09,19.85,22.11,23.44,24.24,25.01,25.24,25.2,25.02,24.76,24.46,24.13,22.37,20.71,19.24,17.94,15.8,14.1,12.72,11.59,10.64,9.832,9.14};
1548	physicsVector = CreatePhysicsVector(E,S_14,factor);
1549	}
1550	else if( ionZ == 15 ) {
1551	G4double S_15[31]={7.309,8.11,9.602,10.98,12.27,13.47,14.6,15.67,16.67,20.71,23.26,24.8,25.76,26.75,27.11,27.16,27.05,26.84,26.57,26.26,24.53,22.85,21.32,19.97,17.68,15.86,14.36,13.12,12.08,11.18,10.41};
1552	physicsVector = CreatePhysicsVector(E,S_15,factor);
1553	}
1554	else if( ionZ == 16 ) {
1555	G4double S_16[31]={7.538,8.363,9.898,11.32,12.65,13.89,15.06,16.17,17.21,21.53,24.36,26.12,27.25,28.46,28.96,29.12,29.08,28.92,28.68,28.4,26.74,25.03,23.47,22.05,19.65,17.7,16.09,14.74,13.6,12.62,11.77};
1556	physicsVector = CreatePhysicsVector(E,S_16,factor);
1557	}
1558	else if( ionZ == 17 ) {
1559	G4double S_17[31]={7.835,8.674,10.23,11.68,13.04,14.31,15.51,16.65,17.73,22.27,25.36,27.34,28.63,30.06,30.7,30.95,30.98,30.87,30.67,30.42,28.81,27.1,25.5,24.04,21.53,19.47,17.76,16.32,15.09,14.03,13.11};
1560	physicsVector = CreatePhysicsVector(E,S_17,factor);
1561	}
1562	else if( ionZ == 18 ) {
1563	G4double S_18[31]={8.007,8.865,10.45,11.92,13.3,14.6,15.82,16.99,18.1,22.84,26.18,28.37,29.83,31.48,32.27,32.61,32.71,32.65,32.49,32.26,30.71,28.99,27.36,25.86,23.27,21.12,19.33,17.81,16.51,15.38,14.4};
1564	physicsVector = CreatePhysicsVector(E,S_18,factor);
1565	}
1566	}
1567
1568	else if( matName == namesMat[20] ) {
1569	if( ionZ == 3 ) {
1570	G4double T_3[31]={2.242,2.448,2.786,3.041,3.229,3.363,3.457,3.52,3.562,3.591,3.503,3.385,3.26,3.021,2.806,2.617,2.45,2.302,2.17,2.052,1.613,1.331,1.134,0.9906,0.7941,0.6663,0.5764,0.5095,0.4576,0.4162,0.3823};
1571	physicsVector = CreatePhysicsVector(E,T_3,factor);
1572	}
1573	else if( ionZ == 4 ) {
1574	G4double T_4[31]={2.756,3.03,3.503,3.891,4.203,4.446,4.632,4.772,4.877,5.093,5.081,4.993,4.877,4.622,4.373,4.14,3.926,3.73,3.552,3.388,2.746,2.306,1.987,1.746,1.409,1.185,1.025,0.9067,0.8141,0.7401,0.6796};
1575	physicsVector = CreatePhysicsVector(E,T_4,factor);
1576	}
1577	else if( ionZ == 5 ) {
1578	G4double T_5[31]={3.161,3.495,4.086,4.592,5.019,5.372,5.66,5.889,6.072,6.535,6.639,6.611,6.525,6.291,6.033,5.779,5.536,5.308,5.094,4.894,4.076,3.482,3.035,2.689,2.19,1.85,1.604,1.419,1.274,1.158,1.063};
1579	physicsVector = CreatePhysicsVector(E,T_5,factor);
1580	}
1581	else if( ionZ == 6 ) {
1582	G4double T_6[31]={3.482,3.87,4.566,5.176,5.709,6.168,6.556,6.88,7.147,7.905,8.161,8.215,8.177,7.987,7.741,7.482,7.224,6.975,6.737,6.51,5.547,4.814,4.244,3.791,3.121,2.653,2.308,2.045,1.838,1.671,1.534};
1583	physicsVector = CreatePhysicsVector(E,T_6,factor);
1584	}
1585	else if( ionZ == 7 ) {
1586	G4double T_7[31]={3.769,4.2,4.981,5.675,6.296,6.847,7.328,7.742,8.095,9.177,9.617,9.778,9.808,9.69,9.475,9.225,8.966,8.707,8.454,8.209,7.131,6.275,5.59,5.034,4.19,3.584,3.132,2.782,2.504,2.279,2.093};
1587	physicsVector = CreatePhysicsVector(E,T_7,factor);
1588	}
1589	else if( ionZ == 8 ) {
1590	G4double T_8[31]={4.026,4.495,5.352,6.122,6.821,7.453,8.019,8.519,8.955,10.37,11.02,11.3,11.41,11.39,11.22,11,10.75,10.49,10.23,9.978,8.813,7.851,7.061,6.405,5.388,4.642,4.074,3.63,3.274,2.983,2.742};
1591	physicsVector = CreatePhysicsVector(E,T_8,factor);
1592	}
1593	else if( ionZ == 9 ) {
1594	G4double T_9[31]={4.235,4.735,5.656,6.491,7.254,7.955,8.593,9.168,9.679,11.42,12.28,12.71,12.9,12.98,12.86,12.67,12.43,12.18,11.91,11.65,10.43,9.391,8.515,7.778,6.613,5.742,5.069,4.536,4.104,3.749,3.451};
1595	physicsVector = CreatePhysicsVector(E,T_9,factor);
1596	}
1597	else if( ionZ == 10 ) {
1598	G4double T_10[31]={4.427,4.954,5.93,6.821,7.641,8.4,9.1,9.742,10.32,12.39,13.48,14.05,14.34,14.53,14.48,14.32,14.1,13.85,13.6,13.33,12.07,10.95,10,9.197,7.895,6.904,6.129,5.508,5.001,4.58,4.225};
1599	physicsVector = CreatePhysicsVector(E,T_10,factor);
1600	}
1601	else if( ionZ == 11 ) {
1602	G4double T_11[31]={4.585,5.138,6.169,7.118,7.997,8.817,9.583,10.29,10.94,13.38,14.75,15.53,15.97,16.35,16.41,16.32,16.16,15.95,15.72,15.48,14.22,13.06,12.03,11.13,9.654,8.497,7.574,6.824,6.205,5.686,5.247};
1603	physicsVector = CreatePhysicsVector(E,T_11,factor);
1604	}
1605	else if( ionZ == 12 ) {
1606	G4double T_12[31]={4.762,5.323,6.375,7.352,8.263,9.115,9.915,10.66,11.36,14.05,15.64,16.55,17.09,17.58,17.69,17.63,17.48,17.27,17.04,16.79,15.48,14.27,13.19,12.25,10.7,9.481,8.502,7.701,7.036,6.476,5.999};
1607	physicsVector = CreatePhysicsVector(E,T_12,factor);
1608	}
1609	else if( ionZ == 13 ) {
1610	G4double T_13[31]={4.93,5.505,6.587,7.597,8.542,9.43,10.26,11.06,11.8,14.76,16.59,17.68,18.35,19.01,19.23,19.23,19.12,18.94,18.73,18.49,17.19,15.94,14.82,13.83,12.16,10.84,9.767,8.881,8.14,7.511,6.972};
1611	physicsVector = CreatePhysicsVector(E,T_13,factor);
1612	}
1613	else if( ionZ == 14 ) {
1614	G4double T_14[31]={5.097,5.685,6.793,7.83,8.806,9.725,10.59,11.42,12.2,15.39,17.45,18.73,19.54,20.37,20.69,20.77,20.7,20.56,20.36,20.14,18.86,17.59,16.42,15.38,13.62,12.2,11.04,10.07,9.264,8.571,7.973};
1615	physicsVector = CreatePhysicsVector(E,T_14,factor);
1616	}
1617	else if( ionZ == 15 ) {
1618	G4double T_15[31]={5.283,5.884,7.018,8.083,9.087,10.03,10.93,11.79,12.61,16.02,18.32,19.79,20.74,21.77,22.22,22.38,22.38,22.28,22.12,21.92,20.69,19.41,18.21,17.12,15.26,13.74,12.48,11.42,10.53,9.763,9.098};
1619	physicsVector = CreatePhysicsVector(E,T_15,factor);
1620	}
1621	else if( ionZ == 16 ) {
1622	G4double T_16[31]={5.44,6.058,7.223,8.317,9.352,10.33,11.26,12.15,12.99,16.62,19.14,20.8,21.9,23.15,23.73,23.98,24.05,24,23.88,23.71,22.54,21.26,20.04,18.91,16.96,15.34,13.99,12.84,11.87,11.02,10.29};
1623	physicsVector = CreatePhysicsVector(E,T_16,factor);
1624	}
1625	else if( ionZ == 17 ) {
1626	G4double T_17[31]={5.636,6.266,7.45,8.566,9.623,10.62,11.58,12.49,13.36,17.16,19.89,21.74,22.99,24.44,25.16,25.5,25.63,25.63,25.54,25.4,24.3,23.03,21.79,20.63,18.59,16.89,15.45,14.23,13.18,12.27,11.48};
1627	physicsVector = CreatePhysicsVector(E,T_17,factor);
1628	}
1629	else if( ionZ == 18 ) {
1630	G4double T_18[31]={5.751,6.393,7.597,8.728,9.803,10.82,11.79,12.73,13.62,17.56,20.49,22.52,23.91,25.56,26.4,26.83,27.02,27.07,27.02,26.9,25.86,24.59,23.33,22.15,20.05,18.28,16.78,15.5,14.4,13.43,12.59};
1631	physicsVector = CreatePhysicsVector(E,T_18,factor);
1632	}
1633	}
1634
1635	else if( matName == namesMat[21] ) {
1636	if( ionZ == 3 ) {
1637	G4double U_3[31]={2.534,2.757,3.121,3.391,3.587,3.725,3.82,3.884,3.924,3.937,3.829,3.69,3.547,3.276,3.035,2.824,2.639,2.476,2.331,2.202,1.725,1.42,1.208,1.054,0.8438,0.7073,0.6114,0.5401,0.485,0.4409,0.4048};
1638	physicsVector = CreatePhysicsVector(E,U_3,factor);
1639	}
1640	else if( ionZ == 4 ) {
1641	G4double U_4[31]={3.127,3.427,3.941,4.355,4.683,4.936,5.128,5.272,5.378,5.583,5.551,5.442,5.304,5.011,4.728,4.466,4.228,4.011,3.814,3.634,2.934,2.458,2.114,1.857,1.496,1.257,1.088,0.9608,0.8624,0.7838,0.7194};
1642	physicsVector = CreatePhysicsVector(E,U_4,factor);
1643	}
1644	else if( ionZ == 5 ) {
1645	G4double U_5[31]={3.592,3.963,4.61,5.155,5.609,5.981,6.28,6.518,6.705,7.168,7.256,7.206,7.097,6.819,6.523,6.234,5.961,5.706,5.468,5.248,4.352,3.709,3.228,2.857,2.324,1.961,1.7,1.503,1.349,1.226,1.125};
1646	physicsVector = CreatePhysicsVector(E,U_5,factor);
1647	}
1648	else if( ionZ == 6 ) {
1649	G4double U_6[31]={3.958,4.39,5.159,5.822,6.393,6.88,7.288,7.625,7.902,8.673,8.918,8.952,8.892,8.656,8.367,8.068,7.776,7.496,7.23,6.978,5.92,5.124,4.51,4.024,3.308,2.809,2.443,2.164,1.945,1.768,1.623};
1650	physicsVector = CreatePhysicsVector(E,U_6,factor);
1651	}
1652	else if( ionZ == 7 ) {
1653	G4double U_7[31]={4.294,4.776,5.641,6.401,7.071,7.659,8.167,8.602,8.969,10.08,10.51,10.66,10.67,10.5,10.24,9.947,9.649,9.355,9.07,8.796,7.607,6.675,5.936,5.338,4.436,3.792,3.312,2.941,2.647,2.409,2.212};
1654	physicsVector = CreatePhysicsVector(E,U_7,factor);
1655	}
1656	else if( ionZ == 8 ) {
1657	G4double U_8[31]={4.583,5.108,6.06,6.907,7.666,8.345,8.947,9.474,9.931,11.4,12.05,12.32,12.41,12.34,12.13,11.86,11.57,11.27,10.98,10.69,9.398,8.348,7.493,6.788,5.701,4.907,4.305,3.835,3.458,3.151,2.895};
1658	physicsVector = CreatePhysicsVector(E,U_8,factor);
1659	}
1660	else if( ionZ == 9 ) {
1661	G4double U_9[31]={4.813,5.373,6.399,7.321,8.155,8.912,9.594,10.2,10.74,12.56,13.43,13.85,14.03,14.06,13.9,13.65,13.37,13.08,12.78,12.49,11.12,9.983,9.035,8.241,6.995,6.067,5.353,4.789,4.333,3.957,3.643};
1662	physicsVector = CreatePhysicsVector(E,U_9,factor);
1663	}
1664	else if( ionZ == 10 ) {
1665	G4double U_10[31]={5.022,5.613,6.702,7.688,8.588,9.413,10.17,10.85,11.47,13.63,14.74,15.31,15.59,15.74,15.64,15.43,15.17,14.88,14.59,14.29,12.87,11.65,10.62,9.745,8.351,7.295,6.472,5.814,5.278,4.833,4.458};
1666	physicsVector = CreatePhysicsVector(E,U_10,factor);
1667	}
1668	else if( ionZ == 11 ) {
1669	G4double U_11[31]={5.194,5.808,6.949,7.992,8.951,9.836,10.65,11.41,12.1,14.61,15.99,16.74,17.14,17.44,17.42,17.26,17.03,16.76,16.47,16.18,14.73,13.44,12.33,11.38,9.835,8.647,7.709,6.952,6.329,5.808,5.368};
1670	physicsVector = CreatePhysicsVector(E,U_11,factor);
1671	}
1672	else if( ionZ == 12 ) {
1673	G4double U_12[31]={5.404,6.031,7.205,8.289,9.292,10.22,11.09,11.9,12.65,15.48,17.11,18.05,18.58,19.05,19.12,19.01,18.81,18.56,18.28,17.99,16.52,15.18,14.01,12.99,11.32,10.02,8.974,8.125,7.421,6.829,6.325};
1674	physicsVector = CreatePhysicsVector(E,U_12,factor);
1675	}
1676	else if( ionZ == 13 ) {
1677	G4double U_13[31]={5.596,6.24,7.447,8.566,9.609,10.58,11.5,12.35,13.15,16.27,18.17,19.29,19.97,20.61,20.78,20.74,20.59,20.37,20.11,19.83,18.36,16.97,15.75,14.67,12.87,11.46,10.31,9.371,8.585,7.92,7.35};
1678	physicsVector = CreatePhysicsVector(E,U_13,factor);
1679	}
1680	else if( ionZ == 14 ) {
1681	G4double U_14[31]={5.788,6.446,7.681,8.832,9.908,10.92,11.87,12.76,13.6,16.98,19.13,20.44,21.26,22.08,22.37,22.4,22.29,22.1,21.87,21.6,20.14,18.72,17.45,16.32,14.42,12.9,11.66,10.63,9.77,9.036,8.404};
1682	physicsVector = CreatePhysicsVector(E,U_14,factor);
1683	}
1684	else if( ionZ == 15 ) {
1685	G4double U_15[31]={6.009,6.682,7.945,9.126,10.23,11.28,12.26,13.19,14.07,17.7,20.09,21.6,22.57,23.6,24.02,24.14,24.09,23.95,23.75,23.51,22.09,20.66,19.34,18.16,16.14,14.51,13.17,12.05,11.1,10.29,9.584};
1686	physicsVector = CreatePhysicsVector(E,U_15,factor);
1687	}
1688	else if( ionZ == 16 ) {
1689	G4double U_16[31]={6.191,6.883,8.18,9.394,10.54,11.61,12.63,13.6,14.52,18.37,21.01,22.72,23.85,25.1,25.65,25.87,25.89,25.8,25.64,25.42,24.07,22.63,21.28,20.06,17.94,16.2,14.76,13.54,12.51,11.61,10.84};
1690	physicsVector = CreatePhysicsVector(E,U_16,factor);
1691	}
1692	else if( ionZ == 17 ) {
1693	G4double U_17[31]={6.421,7.126,8.445,9.68,10.85,11.95,13,13.99,14.94,18.98,21.85,23.76,25.04,26.5,27.19,27.5,27.59,27.55,27.42,27.24,25.94,24.51,23.14,21.87,19.66,17.83,16.3,15,13.89,12.92,12.08};
1694	physicsVector = CreatePhysicsVector(E,U_17,factor);
1695	}
1696	else if( ionZ == 18 ) {
1697	G4double U_18[31]={6.564,7.283,8.624,9.878,11.06,12.19,13.26,14.28,15.25,19.47,22.54,24.64,26.08,27.75,28.59,28.99,29.15,29.15,29.06,28.9,27.66,26.23,24.83,23.53,21.25,19.35,17.74,16.37,15.19,14.17,13.27};
1698	physicsVector = CreatePhysicsVector(E,U_18,factor);
1699	}
1700	}
1701
1702	else if( matName == namesMat[22] ) {
1703	if( ionZ == 3 ) {
1704	G4double V_3[31]={2.541,2.772,3.148,3.426,3.623,3.759,3.849,3.907,3.942,3.94,3.827,3.688,3.546,3.276,3.037,2.827,2.642,2.479,2.335,2.206,1.728,1.422,1.21,1.055,0.8447,0.7079,0.6118,0.5404,0.4851,0.441,0.4049};
1705	physicsVector = CreatePhysicsVector(E,V_3,factor);
1706	}
1707	else if( ionZ == 4 ) {
1708	G4double V_4[31]={3.127,3.434,3.963,4.392,4.729,4.984,5.174,5.313,5.413,5.594,5.553,5.441,5.304,5.014,4.733,4.473,4.236,4.019,3.822,3.642,2.941,2.463,2.119,1.86,1.498,1.259,1.088,0.9614,0.8628,0.784,0.7195};
1709	physicsVector = CreatePhysicsVector(E,V_4,factor);
1710	}
1711	else if( ionZ == 5 ) {
1712	G4double V_5[31]={3.591,3.966,4.629,5.192,5.661,6.041,6.342,6.577,6.759,7.191,7.263,7.208,7.1,6.825,6.532,6.246,5.974,5.72,5.483,5.263,4.366,3.719,3.236,2.863,2.328,1.964,1.702,1.504,1.35,1.227,1.126};
1713	physicsVector = CreatePhysicsVector(E,V_5,factor);
1714	}
1715	else if( ionZ == 6 ) {
1716	G4double V_6[31]={3.96,4.395,5.175,5.856,6.445,6.946,7.362,7.702,7.977,8.714,8.935,8.961,8.899,8.665,8.38,8.085,7.795,7.517,7.251,7,5.94,5.14,4.523,4.035,3.315,2.814,2.446,2.166,1.946,1.769,1.623};
1717	physicsVector = CreatePhysicsVector(E,V_6,factor);
1718	}
1719	else if( ionZ == 7 ) {
1720	G4double V_7[31]={4.304,4.787,5.659,6.433,7.122,7.727,8.249,8.691,9.061,10.14,10.55,10.68,10.68,10.51,10.26,9.969,9.674,9.383,9.099,8.826,7.635,6.699,5.955,5.353,4.446,3.799,3.316,2.944,2.649,2.41,2.213};
1721	physicsVector = CreatePhysicsVector(E,V_7,factor);
1722	}
1723	else if( ionZ == 8 ) {
1724	G4double V_8[31]={4.602,5.127,6.083,6.94,7.716,8.414,9.032,9.572,10.04,11.48,12.09,12.35,12.44,12.36,12.15,11.89,11.6,11.31,11.01,10.73,9.436,8.381,7.52,6.81,5.715,4.916,4.311,3.839,3.461,3.152,2.896};
1725	physicsVector = CreatePhysicsVector(E,V_8,factor);
1726	}
1727	else if( ionZ == 9 ) {
1728	G4double V_9[31]={4.84,5.401,6.429,7.356,8.203,8.977,9.679,10.3,10.85,12.66,13.5,13.88,14.05,14.08,13.92,13.68,13.41,13.12,12.82,12.53,11.16,10.02,9.063,8.264,7.01,6.077,5.359,4.792,4.334,3.958,3.643};
1729	physicsVector = CreatePhysicsVector(E,V_9,factor);
1730	}
1731	else if( ionZ == 10 ) {
1732	G4double V_10[31]={5.057,5.649,6.739,7.728,8.637,9.476,10.25,10.95,11.58,13.75,14.82,15.35,15.62,15.76,15.66,15.46,15.2,14.92,14.62,14.33,12.91,11.68,10.65,9.768,8.365,7.304,6.477,5.817,5.278,4.832,4.457};
1733	physicsVector = CreatePhysicsVector(E,V_10,factor);
1734	}
1735	else if( ionZ == 11 ) {
1736	G4double V_11[31]={5.234,5.85,6.992,8.036,8.999,9.896,10.73,11.5,12.21,14.75,16.08,16.8,17.18,17.46,17.44,17.29,17.06,16.8,16.51,16.22,14.77,13.48,12.36,11.4,9.853,8.658,7.716,6.955,6.33,5.809,5.367};
1737	physicsVector = CreatePhysicsVector(E,V_11,factor);
1738	}
1739	else if( ionZ == 12 ) {
1740	G4double V_12[31]={5.449,6.079,7.253,8.337,9.343,10.28,11.16,11.99,12.75,15.62,17.22,18.12,18.62,19.06,19.13,19.03,18.84,18.59,18.32,18.04,16.57,15.22,14.05,13.02,11.34,10.03,8.983,8.13,7.423,6.83,6.324};
1741	physicsVector = CreatePhysicsVector(E,V_12,factor);
1742	}
1743	else if( ionZ == 13 ) {
1744	G4double V_13[31]={5.647,6.293,7.503,8.623,9.666,10.64,11.57,12.43,13.25,16.42,18.29,19.37,20.01,20.63,20.8,20.76,20.61,20.4,20.15,19.87,18.41,17.02,15.79,14.7,12.9,11.47,10.32,9.377,8.588,7.921,7.35};
1745	physicsVector = CreatePhysicsVector(E,V_13,factor);
1746	}
1747	else if( ionZ == 14 ) {
1748	G4double V_14[31]={5.843,6.505,7.745,8.897,9.973,10.98,11.94,12.84,13.7,17.14,19.26,20.53,21.31,22.1,22.38,22.42,22.31,22.13,21.9,21.64,20.19,18.77,17.49,16.35,14.44,12.91,11.67,10.64,9.772,9.036,8.403};
1749	physicsVector = CreatePhysicsVector(E,V_14,factor);
1750	}
1751	else if( ionZ == 15 ) {
1752	G4double V_15[31]={6.072,6.75,8.021,9.205,10.31,11.35,12.34,13.28,14.17,17.87,20.25,21.71,22.64,23.63,24.04,24.16,24.12,23.98,23.79,23.55,22.15,20.71,19.39,18.2,16.17,14.53,13.18,12.06,11.1,10.29,9.583};
1753	physicsVector = CreatePhysicsVector(E,V_15,factor);
1754	}
1755	else if( ionZ == 16 ) {
1756	G4double V_16[31]={6.259,6.957,8.265,9.483,10.63,11.7,12.72,13.69,14.62,18.55,21.18,22.85,23.93,25.13,25.68,25.89,25.92,25.84,25.68,25.47,24.14,22.7,21.35,20.11,17.98,16.23,14.78,13.55,12.51,11.62,10.84};
1757	physicsVector = CreatePhysicsVector(E,V_16,factor);
1758	}
1759	else if( ionZ == 17 ) {
1760	G4double V_17[31]={6.496,7.208,8.54,9.782,10.95,12.05,13.09,14.09,15.04,19.17,22.04,23.91,25.14,26.54,27.22,27.53,27.62,27.59,27.47,27.29,26.01,24.58,23.2,21.92,19.7,17.86,16.32,15.01,13.89,12.93,12.08};
1761	physicsVector = CreatePhysicsVector(E,V_17,factor);
1762	}
1763	else if( ionZ == 18 ) {
1764	G4double V_18[31]={6.642,7.369,8.724,9.986,11.17,12.3,13.36,14.38,15.35,19.64,22.74,24.8,26.18,27.8,28.61,29,29.17,29.18,29.09,28.94,27.73,26.29,24.89,23.58,21.29,19.37,17.75,16.37,15.19,14.16,13.27};
1765	physicsVector = CreatePhysicsVector(E,V_18,factor);
1766	}
1767	}
1768
1769	else if( matName == namesMat[23] ) {
1770	if( ionZ == 3 ) {
1771	G4double W_3[31]={1.603,1.75,1.999,2.196,2.352,2.474,2.568,2.64,2.694,2.799,2.774,2.705,2.623,2.454,2.295,2.151,2.023,1.908,1.804,1.711,1.358,1.127,0.965,0.8452,0.6807,0.573,0.4969,0.4401,0.396,0.3606,0.3316};
1772	physicsVector = CreatePhysicsVector(E,W_3,factor);
1773	}
1774	else if( ionZ == 4 ) {
1775	G4double W_4[31]={1.976,2.173,2.518,2.808,3.052,3.255,3.421,3.556,3.665,3.952,4.014,3.986,3.921,3.751,3.572,3.4,3.238,3.088,2.949,2.821,2.31,1.952,1.69,1.49,1.209,1.02,0.8848,0.7835,0.7046,0.6415,0.5896};
1776	physicsVector = CreatePhysicsVector(E,W_4,factor);
1777	}
1778	else if( ionZ == 5 ) {
1779	G4double W_5[31]={2.267,2.509,2.939,3.312,3.636,3.916,4.156,4.36,4.531,5.044,5.229,5.268,5.24,5.101,4.924,4.74,4.56,4.387,4.224,4.07,3.425,2.946,2.581,2.295,1.879,1.593,1.385,1.227,1.104,1.005,0.923};
1780	physicsVector = CreatePhysicsVector(E,W_5,factor);
1781	}
1782	else if( ionZ == 6 ) {
1783	G4double W_6[31]={2.498,2.779,3.287,3.735,4.133,4.486,4.798,5.069,5.305,6.07,6.406,6.535,6.56,6.474,6.317,6.136,5.949,5.764,5.584,5.412,4.66,4.073,3.609,3.237,2.679,2.285,1.994,1.77,1.594,1.451,1.333};
1784	physicsVector = CreatePhysicsVector(E,W_6,factor);
1785	}
1786	else if( ionZ == 7 ) {
1787	G4double W_7[31]={2.696,3.008,3.579,4.091,4.554,4.973,5.349,5.685,5.983,7.012,7.524,7.764,7.86,7.853,7.732,7.566,7.384,7.196,7.008,6.824,5.991,5.311,4.757,4.3,3.6,3.092,2.709,2.411,2.174,1.981,1.821};
1788	physicsVector = CreatePhysicsVector(E,W_7,factor);
1789	}
1790	else if( ionZ == 8 ) {
1791	G4double W_8[31]={2.869,3.207,3.834,4.403,4.924,5.401,5.838,6.234,6.591,7.882,8.584,8.951,9.133,9.225,9.155,9.017,8.848,8.665,8.477,8.289,7.399,6.641,6.006,5.472,4.632,4.007,3.527,3.149,2.845,2.595,2.387};
1792	physicsVector = CreatePhysicsVector(E,W_8,factor);
1793	}
1794	else if( ionZ == 9 ) {
1795	G4double W_9[31]={3.02,3.38,4.053,4.67,5.241,5.769,6.256,6.705,7.116,8.654,9.546,10.05,10.32,10.52,10.5,10.4,10.25,10.07,9.886,9.696,8.768,7.95,7.249,6.65,5.69,4.961,4.393,3.939,3.57,3.265,3.009};
1796	physicsVector = CreatePhysicsVector(E,W_9,factor);
1797	}
1798	else if( ionZ == 10 ) {
1799	G4double W_10[31]={3.16,3.539,4.25,4.909,5.523,6.095,6.629,7.125,7.584,9.359,10.44,11.09,11.46,11.78,11.83,11.76,11.63,11.47,11.29,11.1,10.15,9.283,8.527,7.87,6.8,5.972,5.317,4.789,4.355,3.993,3.688};
1800	physicsVector = CreatePhysicsVector(E,W_10,factor);
1801	}
1802	else if( ionZ == 11 ) {
1803	G4double W_11[31]={3.269,3.662,4.406,5.102,5.754,6.367,6.942,7.482,7.987,9.996,11.28,12.08,12.57,13.04,13.17,13.15,13.05,12.91,12.74,12.56,11.6,10.7,9.889,9.18,8.003,7.076,6.334,5.727,5.225,4.803,4.444};
1804	physicsVector = CreatePhysicsVector(E,W_11,factor);
1805	}
1806	else if( ionZ == 12 ) {
1807	G4double W_12[31]={3.392,3.795,4.561,5.284,5.967,6.612,7.222,7.798,8.341,10.56,12.04,12.99,13.6,14.22,14.44,14.47,14.41,14.29,14.14,13.96,13,12.07,11.22,10.47,9.203,8.191,7.369,6.692,6.126,5.647,5.237};
1808	physicsVector = CreatePhysicsVector(E,W_12,factor);
1809	}
1810	else if( ionZ == 13 ) {
1811	G4double W_13[31]={3.506,3.919,4.705,5.45,6.159,6.833,7.473,8.081,8.659,11.06,12.73,13.85,14.58,15.37,15.69,15.79,15.77,15.68,15.54,15.38,14.44,13.49,12.61,11.81,10.46,9.366,8.467,7.719,7.088,6.551,6.089};
1812	physicsVector = CreatePhysicsVector(E,W_13,factor);
1813	}
1814	else if( ionZ == 14 ) {
1815	G4double W_14[31]={3.618,4.04,4.843,5.607,6.337,7.034,7.701,8.337,8.945,11.52,13.37,14.63,15.49,16.46,16.89,17.05,17.08,17.02,16.91,16.76,15.84,14.88,13.97,13.14,11.71,10.54,9.573,8.76,8.07,7.478,6.966};
1816	physicsVector = CreatePhysicsVector(E,W_14,factor);
1817	}
1818	else if( ionZ == 15 ) {
1819	G4double W_15[31]={3.732,4.163,4.982,5.763,6.513,7.233,7.924,8.587,9.221,11.96,13.98,15.41,16.4,17.56,18.11,18.37,18.45,18.43,18.36,18.24,17.38,16.41,15.48,14.63,13.12,11.87,10.82,9.936,9.178,8.524,7.954};
1820	physicsVector = CreatePhysicsVector(E,W_15,factor);
1821	}
1822	else if( ionZ == 16 ) {
1823	G4double W_16[31]={3.835,4.276,5.116,5.916,6.686,7.428,8.142,8.829,9.488,12.38,14.56,16.15,17.28,18.64,19.33,19.67,19.82,19.85,19.81,19.72,18.92,17.97,17.03,16.15,14.58,13.26,12.13,11.17,10.35,9.631,9.005};
1824	physicsVector = CreatePhysicsVector(E,W_16,factor);
1825	}
1826	else if( ionZ == 17 ) {
1827	G4double W_17[31]={3.958,4.407,5.258,6.072,6.857,7.616,8.349,9.056,9.738,12.76,15.1,16.83,18.09,19.65,20.48,20.91,21.12,21.2,21.2,21.13,20.41,19.47,18.52,17.62,15.99,14.6,13.41,12.39,11.5,10.73,10.05};
1828	physicsVector = CreatePhysicsVector(E,W_17,factor);
1829	}
1830	else if( ionZ == 18 ) {
1831	G4double W_18[31]={4.048,4.504,5.369,6.194,6.991,7.764,8.512,9.236,9.936,13.07,15.55,17.43,18.81,20.57,21.52,22.05,22.33,22.45,22.48,22.44,21.78,20.85,19.89,18.97,17.29,15.85,14.6,13.53,12.59,11.77,11.05};
1832	physicsVector = CreatePhysicsVector(E,W_18,factor);
1833	}
1834	}
1835
1836	else if( matName == namesMat[24] ) {
1837	if( ionZ == 3 ) {
1838	G4double X_3[31]={3.332,3.612,4.056,4.373,4.592,4.739,4.835,4.893,4.924,4.88,4.713,4.519,4.326,3.97,3.659,3.39,3.157,2.952,2.773,2.613,2.031,1.663,1.411,1.228,0.9798,0.8194,0.7071,0.6238,0.5594,0.508,0.466};
1839	physicsVector = CreatePhysicsVector(E,X_3,factor);
1840	}
1841	else if( ionZ == 4 ) {
1842	G4double X_4[31]={4.135,4.515,5.152,5.649,6.03,6.313,6.519,6.667,6.771,6.928,6.834,6.663,6.468,6.07,5.698,5.36,5.055,4.781,4.533,4.309,3.45,2.875,2.465,2.159,1.735,1.455,1.256,1.109,0.9941,0.9026,0.8279};
1843	physicsVector = CreatePhysicsVector(E,X_4,factor);
1844	}
1845	else if( ionZ == 5 ) {
1846	G4double X_5[31]={4.766,5.239,6.053,6.721,7.261,7.69,8.024,8.281,8.476,8.912,8.94,8.827,8.657,8.263,7.862,7.483,7.129,6.802,6.501,6.223,5.115,4.335,3.759,3.318,2.69,2.266,1.961,1.732,1.554,1.411,1.294};
1847	physicsVector = CreatePhysicsVector(E,X_5,factor);
1848	}
1849	else if( ionZ == 6 ) {
1850	G4double X_6[31]={5.255,5.811,6.786,7.61,8.302,8.877,9.345,9.723,10.02,10.8,10.99,10.97,10.84,10.48,10.08,9.68,9.295,8.932,8.59,8.271,6.952,5.981,5.244,4.666,3.822,3.239,2.813,2.489,2.236,2.031,1.863};
1851	physicsVector = CreatePhysicsVector(E,X_6,factor);
1852	}
1853	else if( ionZ == 7 ) {
1854	G4double X_7[31]={5.72,6.34,7.441,8.391,9.212,9.916,10.51,11.01,11.42,12.58,12.97,13.06,13,12.71,12.33,11.93,11.53,11.14,10.77,10.42,8.924,7.783,6.891,6.178,5.115,4.363,3.805,3.376,3.037,2.763,2.536};
1855	physicsVector = CreatePhysicsVector(E,X_7,factor);
1856	}
1857	else if( ionZ == 8 ) {
1858	G4double X_8[31]={6.114,6.792,8.008,9.073,10.01,10.83,11.55,12.16,12.68,14.26,14.89,15.11,15.14,14.95,14.61,14.22,13.82,13.42,13.03,12.66,11.02,9.728,8.693,7.849,6.564,5.636,4.937,4.394,3.96,3.606,3.313};
1859	physicsVector = CreatePhysicsVector(E,X_8,factor);
1860	}
1861	else if( ionZ == 9 ) {
1862	G4double X_9[31]={6.408,7.134,8.448,9.611,10.65,11.58,12.4,13.11,13.74,15.73,16.61,16.98,17.1,17.01,16.72,16.36,15.96,15.56,15.17,14.78,13.04,11.63,10.47,9.522,8.047,6.961,6.132,5.48,4.955,4.523,4.163};
1863	physicsVector = CreatePhysicsVector(E,X_9,factor);
1864	}
1865	else if( ionZ == 10 ) {
1866	G4double X_10[31]={6.674,7.439,8.836,10.09,11.21,12.23,13.14,13.96,14.68,17.09,18.23,18.77,19,19.03,18.81,18.47,18.09,17.7,17.3,16.9,15.08,13.56,12.31,11.26,9.601,8.365,7.408,6.648,6.03,5.519,5.09};
1867	physicsVector = CreatePhysicsVector(E,X_10,factor);
1868	}
1869	else if( ionZ == 11 ) {
1870	G4double X_11[31]={6.899,7.695,9.16,10.48,11.69,12.78,13.78,14.69,15.51,18.36,19.81,20.54,20.91,21.1,20.96,20.68,20.33,19.94,19.55,19.15,17.27,15.66,14.3,13.15,11.31,9.914,8.822,7.945,7.227,6.629,6.124};
1871	physicsVector = CreatePhysicsVector(E,X_11,factor);
1872	}
1873	else if( ionZ == 12 ) {
1874	G4double X_12[31]={7.191,8.003,9.507,10.88,12.14,13.3,14.36,15.34,16.23,19.48,21.23,22.17,22.68,23.05,23.01,22.78,22.46,22.1,21.71,21.32,19.39,17.71,16.26,15.03,13.03,11.49,10.27,9.286,8.473,7.791,7.213};
1875	physicsVector = CreatePhysicsVector(E,X_12,factor);
1876	}
1877	else if( ionZ == 13 ) {
1878	G4double X_13[31]={7.454,8.287,9.835,11.26,12.57,13.78,14.9,15.94,16.9,20.51,22.57,23.72,24.38,24.95,25.02,24.86,24.59,24.26,23.89,23.5,21.56,19.81,18.29,16.97,14.82,13.14,11.8,10.71,9.799,9.032,8.377};
1879	physicsVector = CreatePhysicsVector(E,X_13,factor);
1880	}
1881	else if( ionZ == 14 ) {
1882	G4double X_14[31]={7.717,8.569,10.15,11.61,12.97,14.22,15.39,16.49,17.5,21.44,23.79,25.16,25.97,26.73,26.92,26.84,26.62,26.32,25.97,25.6,23.65,21.85,20.27,18.88,16.6,14.8,13.34,12.15,11.15,10.3,9.576};
1883	physicsVector = CreatePhysicsVector(E,X_14,factor);
1884	}
1885	else if( ionZ == 15 ) {
1886	G4double X_15[31]={8.04,8.913,10.54,12.04,13.43,14.73,15.94,17.08,18.14,22.39,25.04,26.63,27.61,28.59,28.92,28.93,28.78,28.52,28.21,27.86,25.94,24.11,22.46,21.01,18.58,16.64,15.06,13.75,12.65,11.71,10.91};
1887	physicsVector = CreatePhysicsVector(E,X_15,factor);
1888	}
1889	else if( ionZ == 16 ) {
1890	G4double X_16[31]={8.294,9.193,10.86,12.4,13.84,15.18,16.44,17.63,18.74,23.29,26.23,28.06,29.21,30.42,30.9,31.02,30.94,30.74,30.46,30.13,28.27,26.41,24.72,23.2,20.63,18.56,16.86,15.44,14.24,13.21,12.32};
1891	physicsVector = CreatePhysicsVector(E,X_16,factor);
1892	}
1893	else if( ionZ == 17 ) {
1894	G4double X_17[31]={8.625,9.54,11.24,12.8,14.27,15.64,16.94,18.16,19.31,24.1,27.32,29.37,30.69,32.13,32.76,32.97,32.96,32.81,32.57,32.27,30.46,28.59,26.85,25.28,22.6,20.42,18.61,17.09,15.8,14.69,13.72};
1895	physicsVector = CreatePhysicsVector(E,X_17,factor);
1896	}
1897	else if( ionZ == 18 ) {
1898	G4double X_18[31]={8.813,9.748,11.47,13.07,14.56,15.96,17.28,18.53,19.71,24.73,28.21,30.48,31.98,33.65,34.43,34.74,34.8,34.7,34.49,34.23,32.47,30.59,28.82,27.21,24.43,22.15,20.26,18.65,17.29,16.1,15.07};
1899	physicsVector = CreatePhysicsVector(E,X_18,factor);
1900	}
1901	}
1902
1903	else if( matName == namesMat[25] ) {
1904	if( ionZ == 3 ) {
1905	G4double Y_3[31]={1.481,1.6,1.8,1.961,2.091,2.197,2.28,2.347,2.399,2.52,2.526,2.488,2.432,2.302,2.17,2.046,1.932,1.829,1.735,1.649,1.321,1.102,0.9471,0.8318,0.6725,0.5676,0.4931,0.4374,0.394,0.3592,0.3306};
1906	physicsVector = CreatePhysicsVector(E,Y_3,factor);
1907	}
1908	else if( ionZ == 4 ) {
1909	G4double Y_4[31]={1.894,2.058,2.337,2.571,2.77,2.938,3.08,3.198,3.296,3.576,3.664,3.67,3.638,3.521,3.379,3.233,3.092,2.958,2.833,2.716,2.242,1.905,1.656,1.464,1.193,1.009,0.8775,0.7783,0.7009,0.6388,0.5877};
1910	physicsVector = CreatePhysicsVector(E,Y_4,factor);
1911	}
1912	else if( ionZ == 5 ) {
1913	G4double Y_5[31]={2.237,2.443,2.798,3.099,3.361,3.59,3.79,3.963,4.112,4.582,4.781,4.853,4.861,4.786,4.657,4.508,4.355,4.203,4.057,3.917,3.321,2.871,2.525,2.252,1.852,1.575,1.372,1.218,1.097,1,0.9198};
1914	physicsVector = CreatePhysicsVector(E,Y_5,factor);
1915	}
1916	else if( ionZ == 6 ) {
1917	G4double Y_6[31]={2.522,2.767,3.196,3.564,3.887,4.174,4.429,4.656,4.856,5.536,5.867,6.022,6.085,6.072,5.972,5.835,5.681,5.521,5.362,5.207,4.515,3.965,3.527,3.172,2.638,2.258,1.974,1.756,1.584,1.444,1.328};
1918	physicsVector = CreatePhysicsVector(E,Y_6,factor);
1919	}
1920	else if( ionZ == 7 ) {
1921	G4double Y_7[31]={2.767,3.047,3.541,3.97,4.349,4.69,4.997,5.274,5.523,6.419,6.901,7.156,7.287,7.359,7.305,7.192,7.049,6.892,6.729,6.565,5.804,5.168,4.645,4.212,3.542,3.052,2.681,2.391,2.16,1.97,1.813};
1922	physicsVector = CreatePhysicsVector(E,Y_7,factor);
1923	}
1924	else if( ionZ == 8 ) {
1925	G4double Y_8[31]={2.982,3.293,3.847,4.333,4.766,5.158,5.514,5.839,6.135,7.245,7.888,8.255,8.465,8.639,8.644,8.566,8.443,8.297,8.137,7.972,7.164,6.457,5.86,5.354,4.552,3.952,3.488,3.121,2.825,2.581,2.377};
1926	physicsVector = CreatePhysicsVector(E,Y_8,factor);
1927	}
1928	else if( ionZ == 9 ) {
1929	G4double Y_9[31]={3.17,3.508,4.116,4.653,5.135,5.573,5.973,6.341,6.68,7.991,8.796,9.279,9.572,9.852,9.921,9.882,9.783,9.65,9.498,9.335,8.499,7.738,7.079,6.512,5.596,4.896,4.347,3.906,3.547,3.248,2.997};
1930	physicsVector = CreatePhysicsVector(E,Y_9,factor);
1931	}
1932	else if( ionZ == 10 ) {
1933	G4double Y_10[31]={3.34,3.702,4.357,4.942,5.469,5.95,6.392,6.801,7.18,8.682,9.65,10.26,10.64,11.04,11.18,11.18,11.11,11,10.86,10.7,9.851,9.046,8.335,7.713,6.693,5.897,5.265,4.752,4.329,3.975,3.676};
1934	physicsVector = CreatePhysicsVector(E,Y_10,factor);
1935	}
1936	else if( ionZ == 11 ) {
1937	G4double Y_11[31]={3.478,3.859,4.556,5.183,5.753,6.275,6.757,7.205,7.623,9.316,10.45,11.19,11.68,12.21,12.44,12.5,12.47,12.38,12.25,12.11,11.26,10.42,9.668,8.996,7.875,6.986,6.27,5.682,5.193,4.78,4.428};
1938	physicsVector = CreatePhysicsVector(E,Y_11,factor);
1939	}
1940	else if( ionZ == 12 ) {
1941	G4double Y_12[31]={3.623,4.019,4.749,5.412,6.018,6.578,7.096,7.578,8.029,9.891,11.19,12.06,12.65,13.32,13.64,13.76,13.77,13.71,13.6,13.47,12.63,11.77,10.98,10.27,9.058,8.087,7.294,6.638,6.087,5.62,5.218};
1942	physicsVector = CreatePhysicsVector(E,Y_12,factor);
1943	}
1944	else if( ionZ == 13 ) {
1945	G4double Y_13[31]={3.754,4.165,4.923,5.616,6.256,6.849,7.401,7.917,8.4,10.42,11.88,12.88,13.58,14.4,14.82,15.01,15.07,15.04,14.96,14.85,14.04,13.16,12.33,11.58,10.3,9.246,8.38,7.655,7.042,6.518,6.066};
1946	physicsVector = CreatePhysicsVector(E,Y_13,factor);
1947	}
1948	else if( ionZ == 14 ) {
1949	G4double Y_14[31]={3.88,4.303,5.086,5.806,6.474,7.097,7.68,8.225,8.738,10.9,12.51,13.65,14.45,15.43,15.95,16.21,16.32,16.34,16.28,16.18,15.41,14.52,13.67,12.89,11.53,10.41,9.478,8.69,8.019,7.442,6.941};
1950	physicsVector = CreatePhysicsVector(E,Y_14,factor);
1951	}
1952	else if( ionZ == 15 ) {
1953	G4double Y_15[31]={4,4.433,5.238,5.982,6.678,7.329,7.94,8.515,9.057,11.36,13.11,14.39,15.31,16.46,17.1,17.45,17.63,17.69,17.67,17.6,16.9,16.02,15.15,14.34,12.92,11.72,10.71,9.853,9.117,8.48,7.924};
1954	physicsVector = CreatePhysicsVector(E,Y_15,factor);
1955	}
1956	else if( ionZ == 16 ) {
1957	G4double Y_16[31]={4.111,4.556,5.384,6.154,6.876,7.554,8.193,8.796,9.366,11.81,13.7,15.11,16.14,17.47,18.23,18.68,18.92,19.03,19.06,19.02,18.4,17.54,16.66,15.83,14.34,13.08,12,11.07,10.27,9.578,8.966};
1958	physicsVector = CreatePhysicsVector(E,Y_16,factor);
1959	}
1960	else if( ionZ == 17 ) {
1961	G4double Y_17[31]={4.234,4.688,5.534,6.322,7.066,7.768,8.431,9.058,9.652,12.22,14.24,15.78,16.92,18.43,19.31,19.85,20.16,20.33,20.39,20.39,19.85,19.01,18.13,17.28,15.74,14.41,13.27,12.28,11.43,10.67,10.01};
1962	physicsVector = CreatePhysicsVector(E,Y_17,factor);
1963	}
1964	else if( ionZ == 18 ) {
1965	G4double Y_18[31]={4.333,4.796,5.658,6.463,7.224,7.946,8.63,9.28,9.896,12.58,14.71,16.37,17.63,19.31,20.31,20.94,21.32,21.54,21.64,21.67,21.21,20.38,19.49,18.63,17.04,15.65,14.46,13.42,12.51,11.72,11.01};
1966	physicsVector = CreatePhysicsVector(E,Y_18,factor);
1967	}
1968	}
1969
1970	else if( matName == namesMat[26] ) {
1971	if( ionZ == 3 ) {
1972	G4double Z_3[31]={1.572,1.764,2.123,2.452,2.743,2.998,3.22,3.417,3.593,4.244,4.656,4.924,5.094,5.266,5.319,5.315,5.281,5.229,5.167,5.099,4.73,4.375,4.056,3.775,3.31,2.946,2.656,2.419,2.222,2.057,1.915};
1973	physicsVector = CreatePhysicsVector(E,Z_3,factor);
1974	}
1975	else if( ionZ == 4 ) {
1976	G4double Z_4[31]={1.032,1.139,1.311,1.439,1.539,1.618,1.683,1.737,1.781,1.917,1.966,1.971,1.955,1.896,1.828,1.759,1.693,1.631,1.571,1.516,1.284,1.113,0.983,0.8813,0.7334,0.6312,0.5563,0.4989,0.4534,0.4163,0.3856};
1977	physicsVector = CreatePhysicsVector(E,Z_4,factor);
1978	}
1979	else if( ionZ == 5 ) {
1980	G4double Z_5[31]={1.181,1.315,1.542,1.72,1.861,1.977,2.073,2.154,2.224,2.452,2.56,2.603,2.61,2.574,2.512,2.443,2.372,2.302,2.234,2.169,1.883,1.66,1.483,1.34,1.127,0.9759,0.863,0.7755,0.7057,0.6485,0.6009};
1981	physicsVector = CreatePhysicsVector(E,Z_5,factor);
1982	}
1983	else if( ionZ == 6 ) {
1984	G4double Z_6[31]={1.3,1.455,1.729,1.954,2.137,2.29,2.419,2.529,2.625,2.955,3.132,3.222,3.26,3.26,3.214,3.152,3.083,3.011,2.939,2.868,2.542,2.273,2.053,1.871,1.591,1.386,1.231,1.11,1.011,0.9306,0.8628};
1985	physicsVector = CreatePhysicsVector(E,Z_6,factor);
1986	}
1987	else if( ionZ == 7 ) {
1988	G4double Z_7[31]={1.399,1.569,1.88,2.146,2.369,2.557,2.718,2.858,2.981,3.417,3.669,3.814,3.892,3.94,3.921,3.873,3.811,3.742,3.67,3.597,3.246,2.941,2.682,2.463,2.117,1.858,1.658,1.499,1.369,1.262,1.171};
1989	physicsVector = CreatePhysicsVector(E,Z_7,factor);
1990	}
1991	else if( ionZ == 8 ) {
1992	G4double Z_8[31]={1.49,1.672,2.01,2.311,2.571,2.793,2.986,3.155,3.305,3.849,4.181,4.385,4.508,4.614,4.628,4.601,4.551,4.49,4.421,4.35,3.985,3.65,3.358,3.105,2.697,2.384,2.138,1.939,1.776,1.64,1.524};
1993	physicsVector = CreatePhysicsVector(E,Z_8,factor);
1994	}
1995	else if( ionZ == 9 ) {
1996	G4double Z_9[31]={1.572,1.764,2.123,2.452,2.743,2.998,3.22,3.417,3.593,4.244,4.656,4.924,5.094,5.266,5.319,5.315,5.281,5.229,5.167,5.099,4.73,4.375,4.056,3.775,3.31,2.946,2.656,2.419,2.222,2.057,1.915};
1997	physicsVector = CreatePhysicsVector(E,Z_9,factor);
1998	}
1999	else if( ionZ == 10 ) {
2000	G4double Z_10[31]={1.652,1.85,2.226,2.576,2.895,3.178,3.429,3.652,3.853,4.609,5.103,5.435,5.657,5.901,5.999,6.024,6.009,5.97,5.917,5.855,5.49,5.121,4.779,4.473,3.957,3.545,3.212,2.937,2.707,2.511,2.343};
2001	physicsVector = CreatePhysicsVector(E,Z_10,factor);
2002	}
2003	else if( ionZ == 11 ) {
2004	G4double Z_11[31]={1.725,1.93,2.319,2.685,3.026,3.335,3.611,3.86,4.085,4.946,5.522,5.921,6.199,6.522,6.671,6.731,6.739,6.717,6.675,6.621,6.271,5.893,5.534,5.205,4.642,4.185,3.81,3.497,3.232,3.006,2.81};
2005	physicsVector = CreatePhysicsVector(E,Z_11,factor);
2006	}
2007	else if( ionZ == 12 ) {
2008	G4double Z_12[31]={1.8,2.011,2.41,2.788,3.146,3.477,3.777,4.049,4.297,5.258,5.915,6.381,6.713,7.118,7.322,7.42,7.454,7.451,7.423,7.379,7.047,6.665,6.292,5.945,5.34,4.842,4.428,4.079,3.782,3.527,3.304};
2009	physicsVector = CreatePhysicsVector(E,Z_12,factor);
2010	}
2011	else if( ionZ == 13 ) {
2012	G4double Z_13[31]={1.872,2.089,2.498,2.885,3.256,3.604,3.925,4.219,4.487,5.546,6.283,6.817,7.206,7.697,7.96,8.099,8.163,8.181,8.17,8.138,7.835,7.455,7.072,6.71,6.067,5.53,5.078,4.695,4.366,4.08,3.83};
2013	physicsVector = CreatePhysicsVector(E,Z_13,factor);
2014	}
2015	else if( ionZ == 14 ) {
2016	G4double Z_14[31]={1.942,2.166,2.583,2.978,3.358,3.72,4.059,4.371,4.66,5.81,6.625,7.226,7.673,8.253,8.577,8.76,8.856,8.897,8.904,8.885,8.619,8.245,7.856,7.481,6.805,6.232,5.746,5.329,4.97,4.656,4.38};
2017	physicsVector = CreatePhysicsVector(E,Z_14,factor);
2018	}
2019	else if( ionZ == 15 ) {
2020	G4double Z_15[31]={2.007,2.237,2.663,3.065,3.452,3.825,4.178,4.508,4.814,6.052,6.945,7.613,8.119,8.792,9.184,9.417,9.55,9.62,9.648,9.648,9.429,9.07,8.681,8.297,7.592,6.984,6.463,6.013,5.622,5.278,4.975};
2021	physicsVector = CreatePhysicsVector(E,Z_15,factor);
2022	}
2023	else if( ionZ == 16 ) {
2024	G4double Z_16[31]={2.071,2.309,2.745,3.153,3.546,3.928,4.293,4.638,4.961,6.283,7.25,7.983,8.548,9.316,9.78,10.06,10.24,10.34,10.39,10.41,10.24,9.904,9.518,9.129,8.4,7.761,7.207,6.725,6.303,5.93,5.599};
2025	physicsVector = CreatePhysicsVector(E,Z_16,factor);
2026	}
2027	else if( ionZ == 17 ) {
2028	G4double Z_17[31]={2.134,2.38,2.826,3.24,3.638,4.025,4.401,4.758,5.095,6.494,7.533,8.331,8.953,9.818,10.35,10.69,10.91,11.04,11.11,11.15,11.05,10.73,10.35,9.958,9.209,8.542,7.959,7.447,6.996,6.597,6.24};
2029	physicsVector = CreatePhysicsVector(E,Z_17,factor);
2030	}
2031	else if( ionZ == 18 ) {
2032	G4double Z_18[31]={2.199,2.452,2.911,3.332,3.735,4.128,4.511,4.879,5.229,6.701,7.809,8.67,9.348,10.31,10.92,11.31,11.56,11.73,11.83,11.89,11.83,11.54,11.17,10.77,10,9.31,8.699,8.16,7.683,7.258,6.877};
2033	physicsVector = CreatePhysicsVector(E,Z_18,factor);
2034	}
2035	}
2036
2037	else if( matName == namesMat[27] ) {
2038	if( ionZ == 3 ) {
2039	G4double XX_3[31]={2.663,2.892,3.261,3.535,3.734,3.875,3.973,4.038,4.08,4.092,3.977,3.829,3.678,3.391,3.138,2.917,2.724,2.553,2.403,2.269,1.774,1.459,1.241,1.082,0.8659,0.7257,0.6271,0.5539,0.4973,0.452,0.415};
2040	physicsVector = CreatePhysicsVector(E,XX_3,factor);
2041	}
2042	else if( ionZ == 4 ) {
2043	G4double XX_4[31]={3.295,3.605,4.129,4.549,4.881,5.138,5.334,5.481,5.59,5.801,5.764,5.645,5.498,5.186,4.887,4.611,4.361,4.135,3.929,3.741,3.016,2.524,2.17,1.905,1.535,1.289,1.115,0.9851,0.8841,0.8034,0.7374};
2044	physicsVector = CreatePhysicsVector(E,XX_4,factor);
2045	}
2046	else if( ionZ == 5 ) {
2047	G4double XX_5[31]={3.789,4.174,4.84,5.395,5.855,6.232,6.535,6.778,6.969,7.445,7.532,7.474,7.356,7.057,6.741,6.436,6.148,5.88,5.632,5.402,4.472,3.807,3.311,2.93,2.382,2.011,1.743,1.541,1.383,1.257,1.153};
2048	physicsVector = CreatePhysicsVector(E,XX_5,factor);
2049	}
2050	else if( ionZ == 6 ) {
2051	G4double XX_6[31]={4.173,4.625,5.421,6.101,6.682,7.175,7.589,7.932,8.214,9.005,9.254,9.283,9.213,8.956,8.645,8.328,8.019,7.723,7.444,7.181,6.08,5.257,4.624,4.125,3.39,2.878,2.503,2.217,1.993,1.812,1.663};
2052	physicsVector = CreatePhysicsVector(E,XX_6,factor);
2053	}
2054	else if( ionZ == 7 ) {
2055	G4double XX_7[31]={4.521,5.025,5.925,6.708,7.394,7.99,8.506,8.946,9.32,10.46,10.9,11.04,11.04,10.86,10.57,10.26,9.946,9.635,9.336,9.049,7.81,6.846,6.084,5.47,4.544,3.884,3.392,3.012,2.711,2.467,2.266};
2056	physicsVector = CreatePhysicsVector(E,XX_7,factor);
2057	}
2058	else if( ionZ == 8 ) {
2059	G4double XX_8[31]={4.828,5.377,6.371,7.249,8.029,8.721,9.332,9.866,10.33,11.83,12.49,12.77,12.86,12.77,12.53,12.24,11.92,11.61,11.3,11,9.647,8.56,7.678,6.953,5.837,5.023,4.407,3.926,3.54,3.226,2.964};
2060	physicsVector = CreatePhysicsVector(E,XX_8,factor);
2061	}
2062	else if( ionZ == 9 ) {
2063	G4double XX_9[31]={5.063,5.651,6.724,7.682,8.543,9.318,10.01,10.63,11.18,13.03,13.93,14.35,14.53,14.54,14.36,14.09,13.79,13.47,13.16,12.85,11.42,10.24,9.259,8.442,7.162,6.211,5.48,4.902,4.435,4.051,3.729};
2064	physicsVector = CreatePhysicsVector(E,XX_9,factor);
2065	}
2066	else if( ionZ == 10 ) {
2067	G4double XX_10[31]={5.278,5.897,7.037,8.065,8.998,9.846,10.62,11.31,11.94,14.13,15.28,15.86,16.14,16.28,16.16,15.92,15.64,15.33,15.02,14.7,13.21,11.95,10.88,9.983,8.55,7.468,6.625,5.951,5.402,4.947,4.564};
2068	physicsVector = CreatePhysicsVector(E,XX_10,factor);
2069	}
2070	else if( ionZ == 11 ) {
2071	G4double XX_11[31]={5.468,6.118,7.326,8.428,9.437,10.36,11.22,12,12.71,15.31,16.77,17.58,18.03,18.38,18.38,18.23,18.01,17.74,17.45,17.15,15.65,14.31,13.14,12.13,10.48,9.207,8.196,7.378,6.704,6.141,5.665};
2072	physicsVector = CreatePhysicsVector(E,XX_11,factor);
2073	}
2074	else if( ionZ == 12 ) {
2075	G4double XX_12[31]={5.677,6.333,7.561,8.693,9.738,10.7,11.6,12.42,13.18,16.06,17.74,18.7,19.24,19.7,19.76,19.63,19.41,19.13,18.84,18.53,16.98,15.58,14.37,13.32,11.6,10.26,9.187,8.317,7.596,6.99,6.474};
2076	physicsVector = CreatePhysicsVector(E,XX_12,factor);
2077	}
2078	else if( ionZ == 13 ) {
2079	G4double XX_13[31]={5.878,6.551,7.813,8.982,10.07,11.08,12.02,12.9,13.72,16.89,18.83,19.98,20.67,21.32,21.48,21.42,21.24,21,20.72,20.42,18.86,17.42,16.15,15.03,13.18,11.73,10.55,9.59,8.785,8.104,7.521};
2080	physicsVector = CreatePhysicsVector(E,XX_13,factor);
2081	}
2082	else if( ionZ == 14 ) {
2083	G4double XX_14[31]={6.078,6.766,8.057,9.258,10.38,11.43,12.42,13.34,14.2,17.63,19.83,21.18,22.01,22.84,23.12,23.13,23,22.79,22.53,22.24,20.7,19.22,17.9,16.72,14.77,13.2,11.93,10.88,9.998,9.247,8.6};
2084	physicsVector = CreatePhysicsVector(E,XX_14,factor);
2085	}
2086	else if( ionZ == 15 ) {
2087	G4double XX_15[31]={6.312,7.015,8.334,9.566,10.72,11.81,12.83,13.79,14.7,18.38,20.82,22.38,23.37,24.41,24.82,24.93,24.86,24.69,24.47,24.21,22.7,21.2,19.83,18.61,16.53,14.86,13.48,12.33,11.36,10.52,9.805};
2088	physicsVector = CreatePhysicsVector(E,XX_15,factor);
2089	}
2090	else if( ionZ == 16 ) {
2091	G4double XX_16[31]={6.503,7.226,8.58,9.845,11.04,12.16,13.22,14.22,15.16,19.09,21.78,23.54,24.69,25.96,26.52,26.72,26.72,26.6,26.41,26.18,24.73,23.22,21.82,20.55,18.37,16.58,15.1,13.85,12.79,11.88,11.08};
2092	physicsVector = CreatePhysicsVector(E,XX_16,factor);
2093	}
2094	else if( ionZ == 17 ) {
2095	G4double XX_17[31]={6.746,7.481,8.856,10.14,11.36,12.51,13.6,14.63,15.61,19.74,22.65,24.61,25.92,27.41,28.11,28.4,28.47,28.4,28.25,28.05,26.66,25.15,23.72,22.41,20.13,18.25,16.67,15.34,14.2,13.22,12.36};
2096	physicsVector = CreatePhysicsVector(E,XX_17,factor);
2097	}
2098	else if( ionZ == 18 ) {
2099	G4double XX_18[31]={6.88,7.63,9.027,10.33,11.57,12.74,13.85,14.91,15.91,20.22,23.34,25.5,26.96,28.67,29.51,29.89,30.03,30.01,29.89,29.71,28.38,26.87,25.42,24.07,21.73,19.77,18.12,16.72,15.52,14.47,13.56};
2100	physicsVector = CreatePhysicsVector(E,XX_18,factor);
2101	}
2102	}
2103
2104	else if( matName == namesMat[28] ) {
2105	if( ionZ == 3 ) {
2106	G4double YY_3[31]={2.7,2.933,3.311,3.589,3.789,3.93,4.026,4.09,4.129,4.134,4.014,3.864,3.711,3.421,3.166,2.942,2.747,2.575,2.423,2.288,1.788,1.47,1.25,1.09,0.8721,0.7307,0.6314,0.5576,0.5005,0.4549,0.4176};
2107	physicsVector = CreatePhysicsVector(E,YY_3,factor);
2108	}
2109	else if( ionZ == 4 ) {
2110	G4double YY_4[31]={3.338,3.654,4.189,4.618,4.955,5.214,5.409,5.555,5.662,5.863,5.82,5.698,5.548,5.232,4.93,4.652,4.4,4.17,3.962,3.773,3.04,2.544,2.187,1.919,1.546,1.298,1.123,0.9916,0.8899,0.8086,0.7421};
2111	physicsVector = CreatePhysicsVector(E,YY_4,factor);
2112	}
2113	else if( ionZ == 5 ) {
2114	G4double YY_5[31]={3.839,4.23,4.909,5.475,5.944,6.325,6.631,6.874,7.064,7.528,7.607,7.545,7.424,7.121,6.802,6.493,6.203,5.932,5.681,5.449,4.509,3.837,3.337,2.952,2.399,2.025,1.754,1.551,1.392,1.265,1.16};
2115	physicsVector = CreatePhysicsVector(E,YY_5,factor);
2116	}
2117	else if( ionZ == 6 ) {
2118	G4double YY_6[31]={4.229,4.687,5.496,6.189,6.782,7.283,7.702,8.048,8.33,9.11,9.349,9.373,9.3,9.037,8.723,8.403,8.09,7.792,7.51,7.244,6.131,5.3,4.66,4.156,3.414,2.898,2.52,2.232,2.006,1.823,1.673};
2119	physicsVector = CreatePhysicsVector(E,YY_6,factor);
2120	}
2121	else if( ionZ == 7 ) {
2122	G4double YY_7[31]={4.586,5.096,6.008,6.805,7.503,8.11,8.634,9.08,9.455,10.58,11.02,11.15,11.15,10.96,10.67,10.36,10.04,9.722,9.419,9.129,7.876,6.902,6.131,5.511,4.576,3.911,3.415,3.032,2.729,2.483,2.28};
2123	physicsVector = CreatePhysicsVector(E,YY_7,factor);
2124	}
2125	else if( ionZ == 8 ) {
2126	G4double YY_8[31]={4.9,5.456,6.463,7.354,8.146,8.851,9.473,10.01,10.48,11.98,12.63,12.9,12.98,12.89,12.64,12.35,12.03,11.71,11.4,11.09,9.73,8.63,7.739,7.006,5.879,5.058,4.436,3.951,3.563,3.246,2.983};
2127	physicsVector = CreatePhysicsVector(E,YY_8,factor);
2128	}
2129	else if( ionZ == 9 ) {
2130	G4double YY_9[31]={5.141,5.736,6.822,7.793,8.667,9.455,10.16,10.79,11.35,13.2,14.08,14.5,14.67,14.68,14.49,14.22,13.91,13.59,13.27,12.96,11.52,10.32,9.33,8.504,7.212,6.253,5.515,4.933,4.463,4.075,3.752};
2131	physicsVector = CreatePhysicsVector(E,YY_9,factor);
2132	}
2133	else if( ionZ == 10 ) {
2134	G4double YY_10[31]={5.361,5.988,7.141,8.182,9.127,9.988,10.77,11.48,12.11,14.32,15.45,16.02,16.3,16.43,16.3,16.06,15.77,15.46,15.14,14.82,13.32,12.04,10.97,10.06,8.609,7.517,6.667,5.988,5.435,4.976,4.591};
2135	physicsVector = CreatePhysicsVector(E,YY_10,factor);
2136	}
2137	else if( ionZ == 11 ) {
2138	G4double YY_11[31]={5.555,6.214,7.436,8.551,9.573,10.51,11.38,12.17,12.89,15.53,16.98,17.77,18.21,18.55,18.56,18.4,18.17,17.9,17.61,17.31,15.8,14.44,13.25,12.23,10.57,9.277,8.256,7.43,6.75,6.182,5.702};
2139	physicsVector = CreatePhysicsVector(E,YY_11,factor);
2140	}
2141	else if( ionZ == 12 ) {
2142	G4double YY_12[31]={5.769,6.434,7.677,8.821,9.878,10.86,11.76,12.6,13.37,16.28,17.95,18.9,19.44,19.88,19.93,19.8,19.57,19.3,19,18.68,17.12,15.7,14.48,13.41,11.68,10.32,9.246,8.369,7.642,7.032,6.512};
2143	physicsVector = CreatePhysicsVector(E,YY_12,factor);
2144	}
2145	else if( ionZ == 13 ) {
2146	G4double YY_13[31]={5.974,6.657,7.934,9.117,10.22,11.24,12.19,13.08,13.91,17.12,19.06,20.2,20.88,21.51,21.67,21.6,21.42,21.17,20.89,20.59,19.02,17.56,16.27,15.14,13.28,11.81,10.62,9.65,8.838,8.152,7.565};
2147	physicsVector = CreatePhysicsVector(E,YY_13,factor);
2148	}
2149	else if( ionZ == 14 ) {
2150	G4double YY_14[31]={6.18,6.877,8.185,9.4,10.54,11.6,12.59,13.52,14.4,17.88,20.07,21.41,22.23,23.05,23.32,23.33,23.19,22.98,22.72,22.43,20.86,19.37,18.03,16.85,14.87,13.29,12.01,10.95,10.06,9.301,8.65};
2151	physicsVector = CreatePhysicsVector(E,YY_14,factor);
2152	}
2153	else if( ionZ == 15 ) {
2154	G4double YY_15[31]={6.421,7.134,8.471,9.718,10.89,11.98,13.02,13.99,14.9,18.64,21.09,22.63,23.61,24.64,25.04,25.14,25.07,24.9,24.67,24.41,22.88,21.37,19.99,18.75,16.65,14.95,13.56,12.4,11.42,10.59,9.861};
2155	physicsVector = CreatePhysicsVector(E,YY_15,factor);
2156	}
2157	else if( ionZ == 16 ) {
2158	G4double YY_16[31]={6.617,7.35,8.723,10,11.21,12.34,13.41,14.42,15.38,19.36,22.06,23.81,24.95,26.21,26.75,26.95,26.95,26.83,26.64,26.4,24.93,23.41,21.99,20.7,18.5,16.69,15.19,13.94,12.87,11.95,11.15};
2159	physicsVector = CreatePhysicsVector(E,YY_16,factor);
2160	}
2161	else if( ionZ == 17 ) {
2162	G4double YY_17[31]={6.866,7.613,9.008,10.31,11.54,12.7,13.8,14.84,15.83,20.02,22.95,24.9,26.2,27.67,28.36,28.65,28.71,28.64,28.49,28.28,26.87,25.35,23.9,22.57,20.27,18.37,16.78,15.44,14.29,13.29,12.43};
2163	physicsVector = CreatePhysicsVector(E,YY_17,factor);
2164	}
2165	else if( ionZ == 18 ) {
2166	G4double YY_18[31]={7.003,7.765,9.183,10.51,11.75,12.94,14.06,15.12,16.14,20.5,23.65,25.8,27.25,28.94,29.76,30.15,30.28,30.26,30.14,29.95,28.6,27.08,25.61,24.25,21.87,19.9,18.23,16.82,15.61,14.56,13.64};
2167	physicsVector = CreatePhysicsVector(E,YY_18,factor);
2168	}
2169	}
2170
2171	else if( matName == namesMat[29] ) {
2172	if( ionZ == 3) {
2173	G4double ZZ_3[53]={2.3193,2.5198,2.8539,3.1164,3.3203,3.4756,3.5914,3.6755,3.7347,3.8125,3.7349,3.6134,3.4818,3.2258,2.9949,2.7909,2.611,2.4517,2.3103,2.1841,1.7151,1.4139,1.2053,1.0525,0.84417,0.70862,0.61317,0.54214,0.48708,0.44305,0.40697,0.29312,0.23208,0.19364,0.16706,0.13252,0.11092,0.09608,0.08522,0.076915,0.07035,0.065026,0.048615,0.040137,0.034964,0.03149,0.027148,0.024579,0.022911,0.021761,0.020937,0.020327,0.019862};
2174	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_3,factor);
2175	}
2176	else if( ionZ == 4) {
2177	G4double ZZ_4[53]={2.872,3.1439,3.6102,3.9967,4.3169,4.5788,4.7897,4.9568,5.0872,5.387,5.4028,5.3185,5.1971,4.9243,4.6549,4.4036,4.1732,3.9629,3.771,3.5957,2.9117,2.4439,2.1062,1.8518,1.4956,1.2587,1.0901,0.96393,0.86589,0.78742,0.72313,0.52053,0.41214,0.34394,0.2968,0.2355,0.19717,0.17081,0.15152,0.13676,0.1251,0.11564,0.086471,0.071399,0.062202,0.056023,0.048303,0.043734,0.040767,0.038723,0.037256,0.036172,0.035345};
2178	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_4,factor);
2179	}
2180	else if( ionZ == 5) {
2181	G4double ZZ_5[53]={3.2922,3.6315,4.2226,4.7242,5.1543,5.5219,5.8323,6.0911,6.3043,6.8888,7.0451,7.0309,6.9445,6.6925,6.4129,6.1372,5.8747,5.6283,5.3983,5.1841,4.3121,3.6826,3.2109,2.8459,2.3203,1.9619,1.7028,1.5072,1.3543,1.2315,1.1307,0.81305,0.64344,0.53693,0.46337,0.36777,0.30797,0.26684,0.23674,0.21371,0.1955,0.18072,0.13517,0.11163,0.097256,0.087601,0.075535,0.068395,0.063757,0.060561,0.058268,0.056574,0.05528};
2182	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_5,factor);
2183	}
2184	else if( ionZ == 6) {
2185	G4double ZZ_6[53]={3.6037,4.0035,4.7125,5.3248,5.8601,6.3287,6.7363,7.0875,7.3872,8.2986,8.635,8.7189,8.6879,8.485,8.2162,7.9331,7.6534,7.3839,7.1272,6.884,5.8573,5.0814,4.4808,4.0044,3.3005,2.808,2.4458,2.1691,1.9511,1.7751,1.6302,1.1714,0.9263,0.77269,0.66676,0.52925,0.44328,0.38415,0.34086,0.30773,0.28154,0.26028,0.19473,0.16083,0.14014,0.12624,0.10886,0.098575,0.091894,0.08729,0.083986,0.081545,0.079682};
2186	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_6,factor);
2187	}
2188	else if( ionZ == 7) {
2189	G4double ZZ_7[53]={3.8821,4.3278,5.1314,5.8371,6.4632,7.021,7.5168,7.9545,8.3378,9.5943,10.145,10.356,10.402,10.278,10.041,9.7677,9.4845,9.2035,8.9301,8.6668,7.5173,6.6126,5.8919,5.308,4.4226,3.7883,3.3138,2.9467,2.655,2.4179,2.2217,1.5965,1.2614,1.0516,0.90715,0.71995,0.60305,0.52268,0.46384,0.41882,0.3832,0.35431,0.26515,0.21903,0.19087,0.17194,0.14829,0.13429,0.12519,0.11892,0.11442,0.1111,0.10856};
2190	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_7,factor);
2191	}
2192	else if( ionZ == 8) {
2193	G4double ZZ_8[53]={4.1215,4.6063,5.4945,6.2868,6.9978,7.6391,8.2175,8.7372,9.201,10.808,11.596,11.955,12.096,12.077,11.89,11.639,11.364,11.081,10.799,10.523,9.2787,8.2615,7.4307,6.7431,5.6787,4.8981,4.3045,3.8393,3.4663,3.161,2.9069,2.0903,1.6501,1.3745,1.1851,0.94004,0.78733,0.68244,0.60568,0.54694,0.50048,0.46278,0.34643,0.28622,0.24945,0.22474,0.19384,0.17555,0.16367,0.15547,0.14959,0.14525,0.14194};
2194	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_8,factor);
2195	}
2196	else if( ionZ == 9) {
2197	G4double ZZ_9[53]={4.2951,4.8107,5.7678,6.6342,7.4196,8.1343,8.7858,9.3786,9.9152,11.857,12.89,13.408,13.652,13.749,13.62,13.398,13.136,12.857,12.573,12.291,10.982,9.88,8.9601,8.1871,6.9687,6.0574,5.3535,4.7954,4.3434,3.9704,3.658,2.642,2.0878,1.7393,1.4994,1.1892,0.99601,0.86336,0.76632,0.69207,0.63334,0.58568,0.43857,0.36242,0.3159,0.28462,0.24552,0.22237,0.20733,0.19696,0.18951,0.18401,0.17981};
2198	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_9,factor);
2199	}
2200	else if( ionZ == 10) {
2201	G4double ZZ_10[53]={4.4513,4.991,6.004,6.9338,7.7852,8.5662,9.284,9.9435,10.547,12.813,14.099,14.791,15.151,15.382,15.321,15.136,14.895,14.626,14.345,14.061,12.705,11.53,10.532,9.6823,8.3208,7.2846,6.4735,5.8234,5.2919,4.8502,4.4779,3.252,2.5745,2.146,1.8503,1.4675,1.2291,1.0654,0.94575,0.85419,0.78178,0.72301,0.54158,0.44763,0.39022,0.35162,0.30335,0.27477,0.25619,0.24339,0.23419,0.2274,0.22222};
2202	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_10,factor);
2203	}
2204	else if( ionZ == 11) {
2205	G4double ZZ_11[53]={4.5914,5.1553,6.2244,7.2203,8.1435,8.9982,9.7906,10.526,11.206,13.848,15.453,16.384,16.916,17.369,17.442,17.344,17.16,16.93,16.675,16.407,15.045,13.799,12.706,11.753,10.187,8.9672,7.9956,7.2072,6.5562,6.0112,5.5493,4.0154,3.1635,2.6262,2.2573,1.7832,1.4904,1.2907,1.1451,1.0339,0.94615,0.87498,0.65548,0.54186,0.47243,0.42574,0.36734,0.33275,0.31027,0.29478,0.28364,0.27542,0.26915};
2206	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_11,factor);
2207	}
2208	else if( ionZ == 12) {
2209	G4double ZZ_12[53]={4.7537,5.3178,6.3962,7.4137,8.3663,9.2554,10.085,10.859,11.581,14.455,16.279,17.373,18.018,18.598,18.727,18.654,18.479,18.25,17.99,17.716,16.313,15.026,13.895,12.907,11.277,9.9981,8.9727,8.1344,7.4376,6.8507,6.35,4.6625,3.7049,3.0916,2.666,2.1135,1.7695,1.5336,1.3613,1.2296,1.1255,1.041,0.7802,0.64511,0.56252,0.50698,0.43749,0.39634,0.36958,0.35114,0.33789,0.3281,0.32063};
2210	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_12,factor);
2211	}
2212	else if( ionZ == 13) {
2213	G4double ZZ_13[53]={4.911,5.4856,6.5852,7.6302,8.6193,9.551,10.426,11.248,12.018,15.152,17.23,18.531,19.33,20.11,20.354,20.352,20.224,20.024,19.785,19.521,18.12,16.795,15.611,14.565,12.819,11.431,10.305,9.3772,8.6003,7.9413,7.3761,5.45,4.342,3.6272,3.129,2.4808,2.0767,1.7997,1.5975,1.4429,1.3207,1.2216,0.91581,0.75739,0.66053,0.59537,0.51383,0.46554,0.43413,0.41248,0.39693,0.38544,0.37667};
2214	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_13,factor);
2215	}
2216	else if( ionZ == 14) {
2217	G4double ZZ_14[53]={5.0697,5.6529,6.77,7.8376,8.8564,9.8229,10.736,11.597,12.409,15.773,18.09,19.594,20.549,21.535,21.9,21.975,21.896,21.73,21.513,21.265,19.878,18.525,17.298,16.202,14.354,12.867,11.651,10.64,9.7876,9.0607,8.4341,6.2756,5.0169,4.1981,3.6247,2.8758,2.4078,2.0868,1.8523,1.6731,1.5315,1.4167,1.0623,0.8787,0.76644,0.69091,0.59637,0.54036,0.50394,0.47883,0.46078,0.44746,0.43728};
2218	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_14,factor);
2219	}
2220	else if( ionZ == 15) {
2221	G4double ZZ_15[53]={5.2616,5.8538,6.9867,8.074,9.1192,10.117,11.065,11.964,12.815,16.396,18.945,20.657,21.779,22.993,23.502,23.672,23.659,23.54,23.356,23.132,21.792,20.425,19.162,18.02,16.066,14.472,13.155,12.051,11.114,10.311,9.6149,7.1917,5.7614,4.8249,4.1667,3.305,2.7661,2.3966,2.1269,1.9215,1.7603,1.6263,1.2197,1.0091,0.88027,0.79361,0.68512,0.62083,0.57902,0.55019,0.52947,0.51417,0.50248};
2222	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_15,factor);
2223	}
2224	else if( ionZ == 16) {
2225	G4double ZZ_16[53]={5.4129,6.0193,7.1761,8.2871,9.36,10.39,11.373,12.308,13.196,16.986,19.762,21.683,22.976,24.431,25.091,25.363,25.421,25.354,25.208,25.013,23.734,22.366,21.074,19.891,17.84,16.146,14.731,13.536,12.516,11.636,10.869,8.1727,6.5617,5.4997,4.7505,3.7671,3.1514,2.7293,2.4215,2.1866,2.0012,1.8509,1.3881,1.1485,1.002,0.90348,0.78008,0.70695,0.65938,0.62657,0.60299,0.58558,0.57229};
2226	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_16,factor);
2227	}
2228	else if( ionZ == 17) {
2229	G4double ZZ_17[53]={5.6171,6.2307,7.3984,8.5209,9.6097,10.661,11.669,12.632,13.551,17.518,20.497,22.615,24.076,25.769,26.58,26.953,27.082,27.066,26.958,26.791,25.579,24.214,22.901,21.685,19.553,17.771,16.271,14.996,13.9,12.949,12.118,9.1688,7.385,6.2003,5.3604,4.2535,3.5588,3.0821,2.7343,2.4689,2.2595,2.0898,1.5673,1.297,1.1318,1.0205,0.88128,0.79874,0.74504,0.708,0.68138,0.66172,0.64671};
2230	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_17,factor);
2231	}
2232	else if( ionZ == 18) {
2233	G4double ZZ_18[53]={5.7158,6.3394,7.5204,8.6525,9.7528,10.82,11.849,12.836,13.78,17.904,21.07,23.375,24.999,26.928,27.889,28.361,28.559,28.592,28.521,28.381,27.228,25.869,24.541,23.299,21.103,19.255,17.686,16.346,15.19,14.183,13.3,10.137,8.2021,6.9062,5.9819,4.757,3.9841,3.4523,3.0636,2.7668,2.5324,2.3425,1.7573,1.4546,1.2694,1.1448,0.98872,0.8962,0.83601,0.79448,0.76464,0.74259,0.72575};
2234	physicsVector = CreatePhysicsVectorFullRange(E,ZZ_18,factor);
2235	}
2236	}
2237
2238	else if( matName == namesMat[30] ) {
2239	if( ionZ == 3) {
2240	G4double Z3_3[31]={2.5756,2.7949,3.1551,3.4304,3.6374,3.7896,3.899,3.9751,4.0258,4.0654,3.9592,3.8156,3.6658,3.3814,3.1293,2.9089,2.716,2.5462,2.3961,2.2626,1.7699,1.456,1.2393,1.0811,0.86577,0.72598,0.6277,0.55463,0.49803,0.45281,0.41577};
2241	physicsVector = CreatePhysicsVector(E,Z3_3,factor);
2242	}
2243	else if( ionZ == 4) {
2244	G4double Z3_4[31]={3.1808,3.4804,3.989,4.4026,4.7367,5.0026,5.2109,5.3716,5.4936,5.7485,5.7287,5.6167,5.4725,5.1634,4.8658,4.5919,4.343,4.1174,3.9126,3.7263,3.0052,2.5166,2.1655,1.9019,1.5337,1.2895,1.1158,0.98607,0.88532,0.80474,0.73875};
2245	physicsVector = CreatePhysicsVector(E,Z3_4,factor);
2246	}
2247	else if( ionZ == 5) {
2248	G4double Z3_5[31]={3.6406,4.016,4.6644,5.207,5.6635,6.0445,6.3589,6.6152,6.8219,7.3586,7.4733,7.4265,7.3134,7.0185,6.705,6.4016,6.1159,5.8498,5.603,5.3742,4.4516,3.7924,3.3012,2.9226,2.379,2.0095,1.7428,1.5416,1.3845,1.2585,1.1551};
2249	physicsVector = CreatePhysicsVector(E,Z3_5,factor);
2250	}
2251	else if( ionZ == 6) {
2252	G4double Z3_6[31]={3.9822,4.4258,5.206,5.8723,6.4462,6.9391,7.359,7.7136,8.0105,8.8754,9.165,9.212,9.1507,8.8994,8.592,8.2768,7.9698,7.6768,7.3997,7.1387,6.0479,5.2336,4.6072,4.1123,3.3838,2.8758,2.5029,2.2184,1.9945,1.8139,1.6652};
2253	physicsVector = CreatePhysicsVector(E,Z3_6,factor);
2254	}
2255	else if( ionZ == 7) {
2256	G4double Z3_7[31]={4.2846,4.7807,5.6683,6.4396,7.1157,7.7088,8.2263,8.6746,9.0603,10.276,10.776,10.946,10.958,10.781,10.502,10.192,9.8785,9.5706,9.2737,8.9896,7.7641,6.8117,6.0585,5.4511,4.534,3.8794,3.3908,3.0134,2.7137,2.4704,2.2692};
2257	physicsVector = CreatePhysicsVector(E,Z3_7,factor);
2258	}
2259	else if( ionZ == 8) {
2260	G4double Z3_8[31]={4.5438,5.0849,6.0686,6.9369,7.7081,8.3949,9.0047,9.5432,10.016,11.593,12.327,12.642,12.746,12.669,12.436,12.146,11.836,11.524,11.216,10.917,9.5846,8.511,7.6411,6.9248,5.8213,5.0153,4.4038,3.9255,3.5425,3.2292,2.9686};
2261	physicsVector = CreatePhysicsVector(E,Z3_8,factor);
2262	}
2263	else if( ionZ == 9) {
2264	G4double Z3_9[31]={4.7328,5.3092,6.3724,7.3248,8.1798,8.9497,9.6421,10.262,10.815,12.742,13.72,14.19,14.394,14.429,14.25,13.986,13.688,13.377,13.065,12.757,11.35,10.183,9.2173,8.4105,7.1453,6.2032,5.4776,4.9034,4.4389,4.0561,3.7356};
2265	physicsVector = CreatePhysicsVector(E,Z3_9,factor);
2266	}
2267	else if( ionZ == 10) {
2268	G4double Z3_10[31]={4.9049,5.5088,6.6366,7.6612,8.5905,9.4349,10.202,10.898,11.525,13.795,15.027,15.667,15.984,16.149,16.035,15.807,15.527,15.223,14.912,14.6,13.136,11.888,10.839,9.9499,8.5338,7.4615,6.6246,5.9551,5.4088,4.9551,4.573};
2269	physicsVector = CreatePhysicsVector(E,Z3_10,factor);
2270	}
2271	else if( ionZ == 11) {
2272	G4double Z3_11[31]={5.0608,5.6915,6.8824,7.9813,8.9898,9.9152,10.765,11.544,12.255,14.923,16.477,17.354,17.84,18.223,18.241,18.096,17.872,17.605,17.317,17.019,15.541,14.216,13.064,12.067,10.439,9.1772,8.1756,7.3646,6.6961,6.137,5.6636};
2273	physicsVector = CreatePhysicsVector(E,Z3_11,factor);
2274	}
2275	else if( ionZ == 12) {
2276	G4double Z3_12[31]={5.2373,5.8701,7.0762,8.204,9.249,10.215,11.11,11.935,12.696,15.62,17.397,18.435,19.032,19.538,19.609,19.488,19.27,19.003,18.709,18.403,16.875,15.501,14.306,13.269,11.57,10.243,9.1835,8.3192,7.6021,6.9989,6.4848};
2277	physicsVector = CreatePhysicsVector(E,Z3_12,factor);
2278	}
2279	else if( ionZ == 13) {
2280	G4double Z3_13[31]={5.4086,6.0533,7.2853,8.4469,9.535,10.55,11.496,12.376,13.191,16.403,18.44,19.685,20.434,21.135,21.319,21.267,21.094,20.855,20.579,20.283,18.75,17.331,16.078,14.978,13.154,11.713,10.548,9.5911,8.791,8.1134,7.5329};
2281	physicsVector = CreatePhysicsVector(E,Z3_13,factor);
2282	}
2283	else if( ionZ == 14) {
2284	G4double Z3_14[31]={5.5814,6.2362,7.4891,8.6787,9.8027,10.858,11.847,12.773,13.637,17.106,19.39,20.838,21.741,22.645,22.947,22.97,22.845,22.638,22.383,22.1,20.575,19.122,17.82,16.666,14.733,13.187,11.928,10.884,10.006,9.2581,8.6142};
2285	physicsVector = CreatePhysicsVector(E,Z3_14,factor);
2286	}
2287	else if( ionZ == 15) {
2288	G4double Z3_15[31]={5.7894,6.4545,7.7257,8.9393,10.095,11.187,12.216,13.184,14.093,17.806,20.334,21.99,23.059,24.187,24.63,24.746,24.686,24.524,24.301,24.042,22.557,21.085,19.742,18.537,16.491,14.832,13.468,12.327,11.362,10.535,9.8194};
2289	physicsVector = CreatePhysicsVector(E,Z3_15,factor);
2290	}
2291	else if( ionZ == 16) {
2292	G4double Z3_16[31]={5.958,6.6384,7.9356,9.1763,10.364,11.493,12.561,13.57,14.521,18.469,21.236,23.104,24.343,25.708,26.3,26.516,26.526,26.414,26.229,25.997,24.568,23.089,21.712,20.462,18.312,16.547,15.081,13.846,12.794,11.888,11.099};
2293	physicsVector = CreatePhysicsVector(E,Z3_16,factor);
2294	}
2295	else if( ionZ == 17) {
2296	G4double Z3_17[31]={6.1773,6.8664,8.1767,9.4315,10.639,11.794,12.892,13.933,14.919,19.07,22.054,24.123,25.528,27.128,27.869,28.184,28.264,28.202,28.053,27.849,26.481,25,23.598,22.311,20.072,18.215,16.659,15.34,14.209,13.23,12.375};
2297	physicsVector = CreatePhysicsVector(E,Z3_17,factor);
2298	}
2299	else if( ionZ == 18) {
2300	G4double Z3_18[31]={6.2903,6.991,8.3169,9.5836,10.806,11.981,13.104,14.172,15.188,19.521,22.708,24.968,26.536,28.37,29.258,29.67,29.817,29.804,29.691,29.512,28.199,26.719,25.297,23.98,21.671,19.742,18.113,16.726,15.532,14.494,13.585};
2301	physicsVector = CreatePhysicsVector(E,Z3_18,factor);
2302	}
2303	}
2304
2305	if(physicsVector != 0) {
2306	dedxKeys.push_back( key );
2307	dedx.push_back( physicsVector );
2308	}
2309
2310	return physicsVector;
2311	}
2312
2313	//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
2314
2315	G4PhysicsVector* G4MaterialStoppingICRU73::GetPhysicsVector(G4int, G4int)
2316	{
2317	G4PhysicsVector* physicsVector = 0;
2318	return physicsVector;
2319	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats: