64 LEN =
new std::vector<G4double*>;
65 HEN =
new std::vector<G4double*>;
71 for(
G4int i=0; i<lens; ++i)
delete[] (*
LEN)[i];
74 for(
G4int i=0; i<hens; ++i)
delete[] (*
HEN)[i];
126 if(lastCS<=0. && pMom>
lastTH)
145 colP.push_back(pMom);
174 static const G4double THmiG=THmin*.001;
178 static const G4double Pmin=THmin+(nL-1)*dP;
181 static const G4double milP=std::log(Pmin);
182 static const G4double malP=std::log(Pmax);
183 static const G4double dlP=(malP-milP)/(nH-1);
184 static const G4double milPG=std::log(.001*Pmin);
196 if(sync<=I)
G4cerr<<
"*!*G4ChipsNetronNuclCS::CalcCrossSect:Sync="<<sync<<
"<="<<I<<
G4endl;
222 G4cerr<<
"***G4ChipsNetronNuclearCS::CalcCrossSect: Sync="<<sync<<
"#"<<I<<
", Z=" <<targZ
223 <<
", N="<<targN<<
", F="<<F<<
G4endl;
231 if (Momentum<
lastTH)
return 0.;
232 else if (Momentum<Pmin)
236 else if (Momentum<Pmax)
246 if(sigma<0.)
return 0.;
254 static const G4int N0=1;
255 static const G4double pZ0N1[4]={1., 0., 0., 1.};
256 static const std::pair<G4int, const G4double*> Z0N1(1,pZ0N1);
257 static const std::pair<G4int, const G4double*> Z0[N0]={Z0N1};
259 static const G4int N1=2;
260 static const G4double pZ1N1[4]={6.E-8, 0., 0., 1.};
261 static const std::pair<G4int, const G4double*> Z1N1(1,pZ1N1);
262 static const G4double pZ1N2[4]={9.E-8, 0., 0., 1.};
263 static const std::pair<G4int, const G4double*> Z1N2(2,pZ1N2);
264 static const std::pair<G4int, const G4double*> Z1[N1]={Z1N1, Z1N2};
266 static const G4int N2=2;
267 static const G4double pZ2N1[4]={1.E-13, 9000., 1.E-4, 2.E-4};
268 static const std::pair<G4int, const G4double*> Z2N1(1,pZ2N1);
269 static const G4double pZ2N2[4]={7.E-4, 0., 0., 1.};
270 static const std::pair<G4int, const G4double*> Z2N2(2,pZ2N2);
271 static const std::pair<G4int, const G4double*> Z2[N2]={Z2N1, Z2N2};
273 static const G4int N3=2;
274 static const G4double pZ3N3[4]={1.E-9, 3200., .021, 1.5E-5};
275 static const std::pair<G4int, const G4double*> Z3N1(3,pZ3N3);
276 static const G4double pZ3N4[4]={3.E-9, 200., .051, 2.5E-4};
277 static const std::pair<G4int, const G4double*> Z3N2(4,pZ3N4);
278 static const std::pair<G4int, const G4double*> Z3[N3]={Z3N1, Z3N2};
280 static const G4int N4=1;
281 static const G4double pZ4N5[4]={9.E-9, 400., .088, 4.E-4};
282 static const std::pair<G4int, const G4double*> Z4N5(5,pZ4N5);
283 static const std::pair<G4int, const G4double*> Z4[N4]={Z4N5};
285 static const G4int N5=2;
286 static const G4double pZ5N5[4]={2.E-10, 2700., .009, 4.E-4};
287 static const std::pair<G4int, const G4double*> Z5N5(5,pZ5N5);
288 static const G4double pZ5N6[4]={2.E-8, 110., .030, 1.E-4};
289 static const std::pair<G4int, const G4double*> Z5N6(6,pZ5N6);
290 static const std::pair<G4int, const G4double*> Z5[N5]={Z5N5, Z5N6};
292 static const G4int N6=2;
293 static const G4double pZ6N6[4]={1.5E-7, 300., .129, 5.E-4};
294 static const std::pair<G4int, const G4double*> Z6N6(6,pZ6N6);
295 static const G4double pZ6N7[4]={1.5E-7, 300., .129, 5.E-4};
296 static const std::pair<G4int, const G4double*> Z6N7(7,pZ6N7);
297 static const std::pair<G4int, const G4double*> Z6[N6]={Z6N6, Z6N7};
299 static const G4int N7=2;
300 static const G4double pZ7N7[4]={5.E-8, 500., .085, 2.E-4};
301 static const std::pair<G4int, const G4double*> Z7N7(7,pZ7N7);
302 static const G4double pZ7N8[4]={5.E-8, 140., .15, 9.E-4};
303 static const std::pair<G4int, const G4double*> Z7N8(8,pZ7N8);
304 static const std::pair<G4int, const G4double*> Z7[N7]={Z7N7, Z7N8};
306 static const G4int N8=3;
307 static const G4double pZ8N8[4]={7.E-8, 0., .021, 1.5E-5};
308 static const std::pair<G4int, const G4double*> Z8N8(8,pZ8N8);
309 static const G4double pZ8N9[4]={2.E-8, 170., .062, 1.E-3};
310 static const std::pair<G4int, const G4double*> Z8N9(9,pZ8N9);
311 static const G4double pZ8N10[4]={1.E-9, 0., .051, 2.5E-4};
312 static const std::pair<G4int, const G4double*> Z8N10(10,pZ8N10);
313 static const std::pair<G4int, const G4double*> Z8[N8]={Z8N8, Z8N9, Z8N10};
315 static const G4int N9=1;
316 static const G4double pZ9N10[4]={1.E-11, 3000., .026, 3.E-5};
317 static const std::pair<G4int, const G4double*> Z9N10(10,pZ9N10);
318 static const std::pair<G4int, const G4double*> Z9[N9]={Z9N10};
320 static const G4int N10=3;
321 static const G4double pZ10N10[4]={4.E-8, 0., .021, 1.5E-5};
322 static const std::pair<G4int, const G4double*> Z10N10(10,pZ10N10);
323 static const G4double pZ10N11[4]={4.E-8, 0., .021, 1.5E-5};
324 static const std::pair<G4int, const G4double*> Z10N11(11,pZ10N11);
325 static const G4double pZ10N12[4]={4.E-8, 0., .051, 2.5E-4};
326 static const std::pair<G4int, const G4double*> Z10N12(12,pZ10N12);
327 static const std::pair<G4int, const G4double*> Z10[N10]={Z10N10, Z10N11, Z10N12};
329 static const G4int N11=1;
330 static const G4double pZ11N12[4]={8.E-10, 500., .05, 3.E-4};
331 static const std::pair<G4int, const G4double*> Z11N12(12,pZ11N12);
332 static const std::pair<G4int, const G4double*> Z11[N11]={Z11N12};
334 static const G4int N12=3;
335 static const G4double pZ12N12[4]={2.E-9, 350., .065, 3.E-4};
336 static const std::pair<G4int, const G4double*> Z12N12(12,pZ12N12);
337 static const G4double pZ12N13[4]={2.E-9, 350., .068, 2.E-4};
338 static const std::pair<G4int, const G4double*> Z12N13(13,pZ12N13);
339 static const G4double pZ12N14[4]={2.E-9, 0., .051, 2.5E-4};
340 static const std::pair<G4int, const G4double*> Z12N14(14,pZ12N14);
341 static const std::pair<G4int, const G4double*> Z12[N12]={Z12N12, Z12N13, Z12N14};
343 static const G4int N13=1;
344 static const G4double pZ13N14[4]={9.E-9, 500., .075, 4.E-4};
345 static const std::pair<G4int, const G4double*> Z13N14(14,pZ13N14);
346 static const std::pair<G4int, const G4double*> Z13[N13]={Z13N14};
348 static const G4int N14=3;
349 static const G4double pZ14N14[4]={4.E-9, 200., .076, 1.E-4};
350 static const std::pair<G4int, const G4double*> Z14N14(14,pZ14N14);
351 static const G4double pZ14N15[4]={6.E-9, 500., .073, 4.E-4};
352 static const std::pair<G4int, const G4double*> Z14N15(15,pZ14N15);
353 static const G4double pZ14N16[4]={4.E-9, 200., .076, 1.E-4};
354 static const std::pair<G4int, const G4double*> Z14N16(16,pZ14N16);
355 static const std::pair<G4int, const G4double*> Z14[N14]={Z14N14, Z14N15, Z14N16};
357 static const G4int N15=1;
358 static const G4double pZ15N16[4]={6.E-9, 550., .077, 2.E-4};
359 static const std::pair<G4int, const G4double*> Z15N16(16,pZ15N16);
360 static const std::pair<G4int, const G4double*> Z15[N15]={Z15N16};
362 static const G4int N16=4;
363 static const G4double pZ16N16[4]={1.5E-8, 500., .087, 5.E-4};
364 static const std::pair<G4int, const G4double*> Z16N16(16,pZ16N16);
365 static const G4double pZ16N17[4]={1.E-8, 300., .07, 4.E-3};
366 static const std::pair<G4int, const G4double*> Z16N17(17,pZ16N17);
367 static const G4double pZ16N18[4]={2.E-8, 300., .094, 3.E-4};
368 static const std::pair<G4int, const G4double*> Z16N18(18,pZ16N18);
369 static const G4double pZ16N20[4]={2.E-8, 200., .11, 3.E-4};
370 static const std::pair<G4int, const G4double*> Z16N20(20,pZ16N20);
371 static const std::pair<G4int, const G4double*> Z16[N16]={Z16N16, Z16N17, Z16N18, Z16N20};
373 static const G4int N17=2;
374 static const G4double pZ17N18[4]={3.E-9, 300., .072, 4.E-4};
375 static const std::pair<G4int, const G4double*> Z17N18(18,pZ17N18);
376 static const G4double pZ17N20[4]={5.E-9, 0., .051, 2.5E-4};
377 static const std::pair<G4int, const G4double*> Z17N20(20,pZ17N20);
378 static const std::pair<G4int, const G4double*> Z17[N17]={Z17N18, Z17N20};
380 static const G4int N18=3;
381 static const G4double pZ18N18[4]={2.5E-9, 300., .074, 2.E-4};
382 static const std::pair<G4int, const G4double*> Z18N18(18,pZ18N18);
383 static const G4double pZ18N20[4]={2.E-8, 400., .084, 4.E-4};
384 static const std::pair<G4int, const G4double*> Z18N20(20,pZ18N20);
385 static const G4double pZ18N22[4]={1.E-9, 100., .065, 2.E-4};
386 static const std::pair<G4int, const G4double*> Z18N22(22,pZ18N22);
387 static const std::pair<G4int, const G4double*> Z18[N18]={Z18N18, Z18N20, Z18N22};
389 static const G4int N19=3;
390 static const G4double pZ19N20[4]={3.E-9, 4., .02, 2.E-4};
391 static const std::pair<G4int, const G4double*> Z19N20(20,pZ19N20);
392 static const G4double pZ19N21[4]={3.E-9, 500., .062, 7.E-4};
393 static const std::pair<G4int, const G4double*> Z19N21(21,pZ19N21);
394 static const G4double pZ19N22[4]={3.E-9, 400., .073, 3.E-4};
395 static const std::pair<G4int, const G4double*> Z19N22(22,pZ19N22);
396 static const std::pair<G4int, const G4double*> Z19[N19]={Z19N20, Z19N21, Z19N22};
398 static const G4int N20=6;
399 static const G4double pZ20N20[4]={3.E-9, 0., .021, 1.5E-5};
400 static const std::pair<G4int, const G4double*> Z20N20(20,pZ20N20);
401 static const G4double pZ20N22[4]={2.E-9, 400., .072, 4.E-4};
402 static const std::pair<G4int, const G4double*> Z20N22(22,pZ20N22);
403 static const G4double pZ20N23[4]={.3E-9, 280., .042, 2.5E-4};
404 static const std::pair<G4int, const G4double*> Z20N23(23,pZ20N23);
405 static const G4double pZ20N24[4]={1.E-9, 300., .062, 2.E-4};
406 static const std::pair<G4int, const G4double*> Z20N24(24,pZ20N24);
407 static const G4double pZ20N26[4]={1.5E-8, 400., .064, 2.E-4};
408 static const std::pair<G4int, const G4double*> Z20N26(26,pZ20N26);
409 static const G4double pZ20N28[4]={7.E-9, 0., .051, 2.5E-4};
410 static const std::pair<G4int, const G4double*> Z20N28(28,pZ20N28);
411 static const std::pair<G4int, const G4double*> Z20[N20]={Z20N20, Z20N22, Z20N23,
412 Z20N24, Z20N26, Z20N28};
414 static const G4int N21=1;
415 static const G4double pZ21N24[4]={5.E-9, 1000., .068, 6.E-4};
416 static const std::pair<G4int, const G4double*> Z21N24(24,pZ21N24);
417 static const std::pair<G4int, const G4double*> Z21[N21]={Z21N24};
419 static const G4int N22=5;
420 static const G4double pZ22N24[4]={4.E-9, 900., .065, 6.E-4};
421 static const std::pair<G4int, const G4double*> Z22N24(24,pZ22N24);
422 static const G4double pZ22N25[4]={4.E-9, 1000., .065, 1.E-3};
423 static const std::pair<G4int, const G4double*> Z22N25(25,pZ22N25);
424 static const G4double pZ22N26[4]={4.E-9, 900., .066, 4.E-4};
425 static const std::pair<G4int, const G4double*> Z22N26(26,pZ22N26);
426 static const G4double pZ22N27[4]={4.E-9, 800., .021, 3.E-4};
427 static const std::pair<G4int, const G4double*> Z22N27(27,pZ22N27);
428 static const G4double pZ22N28[4]={4.E-9, 550., .067, 2.E-4};
429 static const std::pair<G4int, const G4double*> Z22N28(28,pZ22N28);
430 static const std::pair<G4int, const G4double*> Z22[N22]={Z22N24, Z22N25, Z22N26,
433 static const G4int N23=2;
434 static const G4double pZ23N27[4]={4.E-9, 700., .065, 1.E-3};
435 static const std::pair<G4int, const G4double*> Z23N27(27,pZ23N27);
436 static const G4double pZ23N28[4]={4.E-9, 700., .065, 1.E-3};
437 static const std::pair<G4int, const G4double*> Z23N28(28,pZ23N28);
438 static const std::pair<G4int, const G4double*> Z23[N23]={Z23N27, Z23N28};
440 static const G4int N24=4;
441 static const G4double pZ24N26[4]={1.E-9, 750., .056, 2.E-4};
442 static const std::pair<G4int, const G4double*> Z24N26(26,pZ24N26);
443 static const G4double pZ24N28[4]={1.E-9, 350., .061, 1.E-4};
444 static const std::pair<G4int, const G4double*> Z24N28(28,pZ24N28);
445 static const G4double pZ24N29[4]={.4E-9, 650., .056, 1.5E-4};
446 static const std::pair<G4int, const G4double*> Z24N29(29,pZ24N29);
447 static const G4double pZ24N30[4]={1.E-9, 700., .054, 3.E-4};
448 static const std::pair<G4int, const G4double*> Z24N30(30,pZ24N30);
449 static const std::pair<G4int, const G4double*> Z24[N24]={Z24N26, Z24N28, Z24N29, Z24N30};
451 static const G4int N25=1;
452 static const G4double pZ25N30[4]={.3E-9, 650., .042, 3.5E-4};
453 static const std::pair<G4int, const G4double*> Z25N30(30,pZ25N30);
454 static const std::pair<G4int, const G4double*> Z25[N25]={Z25N30};
456 static const G4int N26=4;
457 static const G4double pZ26N28[4]={.9E-9, 200., .062, 1.E-4};
458 static const std::pair<G4int, const G4double*> Z26N28(28,pZ26N28);
459 static const G4double pZ26N30[4]={.9E-9, 1500., .055, 5.E-5};
460 static const std::pair<G4int, const G4double*> Z26N30(30,pZ26N30);
461 static const G4double pZ26N31[4]={.9E-9, 1100., .048, 9.E-4};
462 static const std::pair<G4int, const G4double*> Z26N31(31,pZ26N31);
463 static const G4double pZ26N32[4]={.9E-9, 500., .055, 2.E-4};
464 static const std::pair<G4int, const G4double*> Z26N32(32,pZ26N32);
465 static const std::pair<G4int, const G4double*> Z26[N26]={Z26N28, Z26N30, Z26N31, Z26N32};
467 static const G4int N27=1;
468 static const G4double pZ27N32[4]={.2E-9, 21., .008, 3.E-6};
469 static const std::pair<G4int, const G4double*> Z27N32(32,pZ27N32);
470 static const std::pair<G4int, const G4double*> Z27[N27]={Z27N32};
472 static const G4int N28=5;
473 static const G4double pZ28N30[4]={.3E-9, 0., .021, 1.5E-5};
474 static const std::pair<G4int, const G4double*> Z28N30(30,pZ28N30);
475 static const G4double pZ28N32[4]={.3E-9, 0., .021, 1.5E-5};
476 static const std::pair<G4int, const G4double*> Z28N32(32,pZ28N32);
477 static const G4double pZ28N33[4]={.3E-9, 0., .021, 1.5E-5};
478 static const std::pair<G4int, const G4double*> Z28N33(33,pZ28N33);
479 static const G4double pZ28N34[4]={.3E-9, 700., .0065, 2.E-6};
480 static const std::pair<G4int, const G4double*> Z28N34(34,pZ28N34);
481 static const G4double pZ28N36[4]={.3E-9, 75., .0107, 4.E-6};
482 static const std::pair<G4int, const G4double*> Z28N36(36,pZ28N36);
483 static const std::pair<G4int, const G4double*> Z28[N28]={Z28N30, Z28N32, Z28N33,
486 static const G4int N29=2;
487 static const G4double pZ29N34[4]={.1E-9, 35., .005, 6.E-4};
488 static const std::pair<G4int, const G4double*> Z29N34(34,pZ29N34);
489 static const G4double pZ29N36[4]={.2E-9, 23., .01, 2.E-4};
490 static const std::pair<G4int, const G4double*> Z29N36(36,pZ29N36);
491 static const std::pair<G4int, const G4double*> Z29[N29]={Z29N34, Z29N36};
493 static const G4int N30=5;
494 static const G4double pZ30N34[4]={.2E-9, 140., .02, 8.E-6};
495 static const std::pair<G4int, const G4double*> Z30N34(34,pZ30N34);
496 static const G4double pZ30N36[4]={.2E-9, 140., .02, 8.E-6};
497 static const std::pair<G4int, const G4double*> Z30N36(36,pZ30N36);
498 static const G4double pZ30N37[4]={.2E-9, 140., .02, 8.E-6};
499 static const std::pair<G4int, const G4double*> Z30N37(37,pZ30N37);
500 static const G4double pZ30N38[4]={.2E-9, 140., .02, 8.E-6};
501 static const std::pair<G4int, const G4double*> Z30N38(38,pZ30N38);
502 static const G4double pZ30N40[4]={.2E-9, 140., .02, 8.E-6};
503 static const std::pair<G4int, const G4double*> Z30N40(40,pZ30N40);
504 static const std::pair<G4int, const G4double*> Z30[N30]={Z30N34, Z30N36, Z30N37,
507 static const G4int N31=2;
508 static const G4double pZ31N38[4]={.3E-9, 450., .050, 3.E-4};
509 static const std::pair<G4int, const G4double*> Z31N38(38,pZ31N38);
510 static const G4double pZ31N40[4]={.3E-9, 600., .048, 2.E-4};
511 static const std::pair<G4int, const G4double*> Z31N40(40,pZ31N40);
512 static const std::pair<G4int, const G4double*> Z31[N31]={Z31N38, Z31N40};
514 static const G4int N32=5;
515 static const G4double pZ32N38[4]={.2E-9, 200., .05, 2.E-4};
516 static const std::pair<G4int, const G4double*> Z32N38(38,pZ32N38);
517 static const G4double pZ32N40[4]={.2E-9, 600., .05, 2.E-4};
518 static const std::pair<G4int, const G4double*> Z32N40(40,pZ32N40);
519 static const G4double pZ32N41[4]={1.5E-11, 600., .028, 3.E-4};
520 static const std::pair<G4int, const G4double*> Z32N41(41,pZ32N41);
521 static const G4double pZ32N42[4]={9.E-11, 400., .048, 3.E-4};
522 static const std::pair<G4int, const G4double*> Z32N42(42,pZ32N42);
523 static const G4double pZ32N44[4]={9.E-11, 400., .043, 3.E-4};
524 static const std::pair<G4int, const G4double*> Z32N44(44,pZ32N44);
525 static const std::pair<G4int, const G4double*> Z32[N32]={Z32N38, Z32N40, Z32N41,
528 static const G4int N33=1;
529 static const G4double pZ33N42[4]={1.E-11, 1000., .032, 1.E-4};
530 static const std::pair<G4int, const G4double*> Z33N42(42,pZ33N42);
531 static const std::pair<G4int, const G4double*> Z33[N33]={Z33N42};
533 static const G4int N34=6;
534 static const G4double pZ34N40[4]={3.E-12, 500., .01, 2.5E-5};
535 static const std::pair<G4int, const G4double*> Z34N40(40,pZ34N40);
536 static const G4double pZ34N42[4]={3.E-12, 500., .01, 2.5E-5};
537 static const std::pair<G4int, const G4double*> Z34N42(42,pZ34N42);
538 static const G4double pZ34N43[4]={3.E-12, 500., .01, 2.5E-5};
539 static const std::pair<G4int, const G4double*> Z34N43(43,pZ34N43);
540 static const G4double pZ34N44[4]={3.E-12, 500., .01, 2.5E-5};
541 static const std::pair<G4int, const G4double*> Z34N44(44,pZ34N44);
542 static const G4double pZ34N46[4]={3.E-12, 500., .01, 2.5E-5};
543 static const std::pair<G4int, const G4double*> Z34N46(46,pZ34N46);
544 static const G4double pZ34N48[4]={3.E-12, 500., .01, 2.5E-4};
545 static const std::pair<G4int, const G4double*> Z34N48(48,pZ34N48);
546 static const std::pair<G4int, const G4double*> Z34[N34]={Z34N40, Z34N42, Z34N43,
547 Z34N44, Z34N46, Z34N48};
549 static const G4int N35=2;
550 static const G4double pZ35N44[4]={3.E-12, 500., .01, 2.5E-5};
551 static const std::pair<G4int, const G4double*> Z35N44(44,pZ35N44);
552 static const G4double pZ35N46[4]={3.E-12, 500., .01, 2.5E-4};
553 static const std::pair<G4int, const G4double*> Z35N46(46,pZ35N46);
554 static const std::pair<G4int, const G4double*> Z35[N35]={Z35N44, Z35N46};
556 static const G4int N36=6;
557 static const G4double pZ36N42[4]={3.E-12, 500., .01, 2.5E-5};
558 static const std::pair<G4int, const G4double*> Z36N42(42,pZ36N42);
559 static const G4double pZ36N44[4]={3.E-12, 500., .01, 2.5E-5};
560 static const std::pair<G4int, const G4double*> Z36N44(44,pZ36N44);
561 static const G4double pZ36N46[4]={3.E-12, 500., .01, 2.5E-5};
562 static const std::pair<G4int, const G4double*> Z36N46(46,pZ36N46);
563 static const G4double pZ36N47[4]={3.E-12, 500., .01, 2.5E-5};
564 static const std::pair<G4int, const G4double*> Z36N47(47,pZ36N47);
565 static const G4double pZ36N48[4]={3.E-12, 500., .01, 2.5E-5};
566 static const std::pair<G4int, const G4double*> Z36N48(48,pZ36N48);
567 static const G4double pZ36N50[4]={3.E-12, 500., .01, 2.5E-4};
568 static const std::pair<G4int, const G4double*> Z36N50(50,pZ36N50);
569 static const std::pair<G4int, const G4double*> Z36[N36]={Z36N42, Z36N44, Z36N46,
570 Z36N47, Z36N48, Z36N50};
572 static const G4int N37=2;
573 static const G4double pZ37N48[4]={3.E-12, 500., .01, 2.5E-5};
574 static const std::pair<G4int, const G4double*> Z37N48(48,pZ37N48);
575 static const G4double pZ37N50[4]={3.E-12, 500., .01, 2.5E-4};
576 static const std::pair<G4int, const G4double*> Z37N50(50,pZ37N50);
577 static const std::pair<G4int, const G4double*> Z37[N37]={Z37N48, Z37N50};
579 static const G4int N38=4;
580 static const G4double pZ38N46[4]={3.E-12, 500., .01, 2.5E-5};
581 static const std::pair<G4int, const G4double*> Z38N46(46,pZ38N46);
582 static const G4double pZ38N48[4]={3.E-12, 500., .01, 2.5E-5};
583 static const std::pair<G4int, const G4double*> Z38N48(48,pZ38N48);
584 static const G4double pZ38N49[4]={3.E-12, 500., .01, 2.5E-5};
585 static const std::pair<G4int, const G4double*> Z38N49(49,pZ38N49);
586 static const G4double pZ38N50[4]={3.E-12, 500., .01, 2.5E-4};
587 static const std::pair<G4int, const G4double*> Z38N50(50,pZ38N50);
588 static const std::pair<G4int, const G4double*> Z38[N38]={Z38N46, Z38N48, Z38N49, Z38N50};
590 static const G4int N39=1;
591 static const G4double pZ39N50[4]={3.E-12, 500., .01, 2.5E-4};
592 static const std::pair<G4int, const G4double*> Z39N50(50,pZ39N50);
593 static const std::pair<G4int, const G4double*> Z39[N39]={Z39N50};
595 static const G4int N40=5;
596 static const G4double pZ40N50[4]={3.E-12, 500., .01, 2.5E-5};
597 static const std::pair<G4int, const G4double*> Z40N50(50,pZ40N50);
598 static const G4double pZ40N51[4]={3.E-12, 500., .01, 2.5E-5};
599 static const std::pair<G4int, const G4double*> Z40N51(51,pZ40N51);
600 static const G4double pZ40N52[4]={3.E-12, 500., .01, 2.5E-5};
601 static const std::pair<G4int, const G4double*> Z40N52(52,pZ40N52);
602 static const G4double pZ40N54[4]={3.E-12, 500., .01, 2.5E-5};
603 static const std::pair<G4int, const G4double*> Z40N54(54,pZ40N54);
604 static const G4double pZ40N56[4]={3.E-12, 500., .01, 2.5E-4};
605 static const std::pair<G4int, const G4double*> Z40N56(56,pZ40N56);
606 static const std::pair<G4int, const G4double*> Z40[N40]={Z40N50, Z40N51, Z40N52,
609 static const G4int N41=1;
610 static const G4double pZ41N52[4]={3.E-12, 500., .01, 2.5E-4};
611 static const std::pair<G4int, const G4double*> Z41N52(52,pZ41N52);
612 static const std::pair<G4int, const G4double*> Z41[N41]={Z41N52};
614 static const G4int N42=7;
615 static const G4double pZ42N50[4]={3.E-12, 500., .01, 2.5E-5};
616 static const std::pair<G4int, const G4double*> Z42N50(50,pZ42N50);
617 static const G4double pZ42N52[4]={3.E-12, 500., .01, 2.5E-5};
618 static const std::pair<G4int, const G4double*> Z42N52(52,pZ42N52);
619 static const G4double pZ42N53[4]={3.E-12, 500., .01, 2.5E-5};
620 static const std::pair<G4int, const G4double*> Z42N53(53,pZ42N53);
621 static const G4double pZ42N54[4]={3.E-12, 500., .01, 2.5E-5};
622 static const std::pair<G4int, const G4double*> Z42N54(54,pZ42N54);
623 static const G4double pZ42N55[4]={3.E-12, 500., .01, 2.5E-5};
624 static const std::pair<G4int, const G4double*> Z42N55(55,pZ42N55);
625 static const G4double pZ42N56[4]={3.E-12, 500., .01, 2.5E-5};
626 static const std::pair<G4int, const G4double*> Z42N56(56,pZ42N56);
627 static const G4double pZ42N58[4]={3.E-12, 500., .01, 2.5E-4};
628 static const std::pair<G4int, const G4double*> Z42N58(58,pZ42N58);
629 static const std::pair<G4int, const G4double*> Z42[N42]={Z42N50, Z42N52, Z42N53, Z42N54,
630 Z42N55, Z42N56, Z42N58};
632 static const G4int N43=1;
633 static const G4double pZ43N0[4]={3.E-12, 500., .01, 2.5E-4};
634 static const std::pair<G4int, const G4double*> Z43N0(0,pZ43N0);
635 static const std::pair<G4int, const G4double*> Z43[N43]={Z43N0};
637 static const G4int N44=7;
638 static const G4double pZ44N52[4]={3.E-12, 500., .01, 2.5E-5};
639 static const std::pair<G4int, const G4double*> Z44N52(52,pZ44N52);
640 static const G4double pZ44N54[4]={3.E-12, 500., .01, 2.5E-5};
641 static const std::pair<G4int, const G4double*> Z44N54(54,pZ44N54);
642 static const G4double pZ44N55[4]={3.E-12, 500., .01, 2.5E-5};
643 static const std::pair<G4int, const G4double*> Z44N55(55,pZ44N55);
644 static const G4double pZ44N56[4]={3.E-12, 500., .01, 2.5E-5};
645 static const std::pair<G4int, const G4double*> Z44N56(56,pZ44N56);
646 static const G4double pZ44N57[4]={3.E-12, 500., .01, 2.5E-5};
647 static const std::pair<G4int, const G4double*> Z44N57(57,pZ44N57);
648 static const G4double pZ44N58[4]={3.E-12, 500., .01, 2.5E-5};
649 static const std::pair<G4int, const G4double*> Z44N58(58,pZ44N58);
650 static const G4double pZ44N60[4]={3.E-12, 500., .01, 2.5E-4};
651 static const std::pair<G4int, const G4double*> Z44N60(60,pZ44N60);
652 static const std::pair<G4int, const G4double*> Z44[N44]={Z44N52, Z44N54, Z44N55, Z44N56,
653 Z44N57, Z44N58, Z44N60};
655 static const G4int N45=1;
656 static const G4double pZ45N58[4]={3.E-12, 500., .01, 2.5E-4};
657 static const std::pair<G4int, const G4double*> Z45N58(58,pZ45N58);
658 static const std::pair<G4int, const G4double*> Z45[N45]={Z45N58};
660 static const G4int N46=6;
661 static const G4double pZ46N56[4]={3.E-12, 500., .01, 2.5E-5};
662 static const std::pair<G4int, const G4double*> Z46N56(56,pZ46N56);
663 static const G4double pZ46N58[4]={3.E-12, 500., .01, 2.5E-5};
664 static const std::pair<G4int, const G4double*> Z46N58(58,pZ46N58);
665 static const G4double pZ46N59[4]={3.E-12, 500., .01, 2.5E-5};
666 static const std::pair<G4int, const G4double*> Z46N59(59,pZ46N59);
667 static const G4double pZ46N60[4]={3.E-12, 500., .01, 2.5E-5};
668 static const std::pair<G4int, const G4double*> Z46N60(60,pZ46N60);
669 static const G4double pZ46N62[4]={3.E-12, 500., .01, 2.5E-5};
670 static const std::pair<G4int, const G4double*> Z46N62(62,pZ46N62);
671 static const G4double pZ46N64[4]={3.E-12, 500., .01, 2.5E-4};
672 static const std::pair<G4int, const G4double*> Z46N64(64,pZ46N64);
673 static const std::pair<G4int, const G4double*> Z46[N46]={Z46N56, Z46N58, Z46N59,
674 Z46N60, Z46N62, Z46N64};
676 static const G4int N47=2;
677 static const G4double pZ47N60[4]={3.E-12, 500., .01, 2.7E-5};
678 static const std::pair<G4int, const G4double*> Z47N60(60,pZ47N60);
679 static const G4double pZ47N62[4]={3.E-12, 480., .01, 2.5E-4};
680 static const std::pair<G4int, const G4double*> Z47N62(62,pZ47N62);
681 static const std::pair<G4int, const G4double*> Z47[N47]={Z47N60, Z47N62};
683 static const G4int N48=8;
684 static const G4double pZ48N58[4]={3.E-12, 500., .01, 2.5E-5};
685 static const std::pair<G4int, const G4double*> Z48N58(58,pZ48N58);
686 static const G4double pZ48N60[4]={3.E-12, 500., .01, 2.5E-5};
687 static const std::pair<G4int, const G4double*> Z48N60(60,pZ48N60);
688 static const G4double pZ48N62[4]={3.E-12, 500., .01, 2.5E-5};
689 static const std::pair<G4int, const G4double*> Z48N62(62,pZ48N62);
690 static const G4double pZ48N63[4]={3.E-12, 500., .01, 2.5E-5};
691 static const std::pair<G4int, const G4double*> Z48N63(63,pZ48N63);
692 static const G4double pZ48N64[4]={3.E-12, 500., .01, 2.5E-5};
693 static const std::pair<G4int, const G4double*> Z48N64(64,pZ48N64);
694 static const G4double pZ48N65[4]={3.E-12, 500., .01, 2.5E-5};
695 static const std::pair<G4int, const G4double*> Z48N65(65,pZ48N65);
696 static const G4double pZ48N66[4]={3.E-12, 500., .01, 2.5E-5};
697 static const std::pair<G4int, const G4double*> Z48N66(66,pZ48N66);
698 static const G4double pZ48N68[4]={3.E-12, 500., .01, 2.5E-4};
699 static const std::pair<G4int, const G4double*> Z48N68(68,pZ48N68);
700 static const std::pair<G4int, const G4double*> Z48[N48]={Z48N58, Z48N60, Z48N62, Z48N63,
701 Z48N64, Z48N65, Z48N66, Z48N68};
703 static const G4int N49=2;
704 static const G4double pZ49N64[4]={3.E-12, 500., .01, 2.5E-5};
705 static const std::pair<G4int, const G4double*> Z49N64(64,pZ49N64);
706 static const G4double pZ49N66[4]={3.E-12, 500., .01, 2.5E-4};
707 static const std::pair<G4int, const G4double*> Z49N66(66,pZ49N66);
708 static const std::pair<G4int, const G4double*> Z49[N49]={Z49N64, Z49N66};
710 static const G4int N50=10;
711 static const G4double pZ50N62[4]={3.E-12, 500., .01, 2.5E-5};
712 static const std::pair<G4int, const G4double*> Z50N62(62,pZ50N62);
713 static const G4double pZ50N64[4]={3.E-12, 500., .01, 2.5E-5};
714 static const std::pair<G4int, const G4double*> Z50N64(64,pZ50N64);
715 static const G4double pZ50N65[4]={3.E-12, 500., .01, 2.5E-5};
716 static const std::pair<G4int, const G4double*> Z50N65(65,pZ50N65);
717 static const G4double pZ50N66[4]={3.E-12, 500., .01, 2.5E-5};
718 static const std::pair<G4int, const G4double*> Z50N66(66,pZ50N66);
719 static const G4double pZ50N67[4]={3.E-12, 500., .01, 2.5E-5};
720 static const std::pair<G4int, const G4double*> Z50N67(67,pZ50N67);
721 static const G4double pZ50N68[4]={3.E-12, 500., .01, 2.5E-5};
722 static const std::pair<G4int, const G4double*> Z50N68(68,pZ50N68);
723 static const G4double pZ50N69[4]={3.E-12, 500., .01, 2.5E-5};
724 static const std::pair<G4int, const G4double*> Z50N69(69,pZ50N69);
725 static const G4double pZ50N70[4]={3.E-12, 500., .01, 2.5E-5};
726 static const std::pair<G4int, const G4double*> Z50N70(70,pZ50N70);
727 static const G4double pZ50N72[4]={3.E-12, 500., .01, 2.5E-5};
728 static const std::pair<G4int, const G4double*> Z50N72(72,pZ50N72);
729 static const G4double pZ50N74[4]={3.E-12, 500., .01, 2.5E-4};
730 static const std::pair<G4int, const G4double*> Z50N74(74,pZ50N74);
731 static const std::pair<G4int, const G4double*> Z50[N50]={Z50N62, Z50N64, Z50N65, Z50N66,
732 Z50N67, Z50N68, Z50N69, Z50N70,
735 static const G4int N51=2;
736 static const G4double pZ51N70[4]={3.E-12, 500., .01, 2.5E-5};
737 static const std::pair<G4int, const G4double*> Z51N70(70,pZ51N70);
738 static const G4double pZ51N72[4]={3.E-12, 500., .01, 2.5E-4};
739 static const std::pair<G4int, const G4double*> Z51N72(72,pZ51N72);
740 static const std::pair<G4int, const G4double*> Z51[N51]={Z51N70, Z51N72};
742 static const G4int N52=8;
743 static const G4double pZ52N68[4]={3.E-12, 500., .01, 2.5E-5};
744 static const std::pair<G4int, const G4double*> Z52N68(68,pZ52N68);
745 static const G4double pZ52N70[4]={3.E-12, 500., .01, 2.5E-5};
746 static const std::pair<G4int, const G4double*> Z52N70(70,pZ52N70);
747 static const G4double pZ52N71[4]={3.E-12, 500., .01, 2.5E-5};
748 static const std::pair<G4int, const G4double*> Z52N71(71,pZ52N71);
749 static const G4double pZ52N72[4]={3.E-12, 500., .01, 2.5E-5};
750 static const std::pair<G4int, const G4double*> Z52N72(72,pZ52N72);
751 static const G4double pZ52N73[4]={3.E-12, 500., .01, 2.5E-5};
752 static const std::pair<G4int, const G4double*> Z52N73(73,pZ52N73);
753 static const G4double pZ52N74[4]={3.E-12, 500., .01, 2.5E-5};
754 static const std::pair<G4int, const G4double*> Z52N74(74,pZ52N74);
755 static const G4double pZ52N76[4]={3.E-12, 500., .01, 2.5E-5};
756 static const std::pair<G4int, const G4double*> Z52N76(76,pZ52N76);
757 static const G4double pZ52N78[4]={3.E-12, 500., .01, 2.5E-4};
758 static const std::pair<G4int, const G4double*> Z52N78(78,pZ52N78);
759 static const std::pair<G4int, const G4double*> Z52[N52]={Z52N68, Z52N70, Z52N71, Z52N72,
760 Z52N73, Z52N74, Z52N76, Z52N78};
762 static const G4int N53=1;
763 static const G4double pZ53N74[4]={3.E-12, 500., .01, 2.5E-4};
764 static const std::pair<G4int, const G4double*> Z53N74(74,pZ53N74);
765 static const std::pair<G4int, const G4double*> Z53[N53]={Z53N74};
767 static const G4int N54=9;
768 static const G4double pZ54N70[4]={3.E-12, 500., .01, 2.5E-5};
769 static const std::pair<G4int, const G4double*> Z54N70(70,pZ54N70);
770 static const G4double pZ54N72[4]={3.E-12, 500., .01, 2.5E-5};
771 static const std::pair<G4int, const G4double*> Z54N72(72,pZ54N72);
772 static const G4double pZ54N74[4]={3.E-12, 500., .01, 2.5E-5};
773 static const std::pair<G4int, const G4double*> Z54N74(74,pZ54N74);
774 static const G4double pZ54N75[4]={3.E-12, 500., .01, 2.5E-5};
775 static const std::pair<G4int, const G4double*> Z54N75(75,pZ54N75);
776 static const G4double pZ54N76[4]={3.E-12, 500., .01, 2.5E-5};
777 static const std::pair<G4int, const G4double*> Z54N76(76,pZ54N76);
778 static const G4double pZ54N77[4]={3.E-12, 500., .01, 2.5E-5};
779 static const std::pair<G4int, const G4double*> Z54N77(77,pZ54N77);
780 static const G4double pZ54N78[4]={3.E-12, 500., .01, 2.5E-5};
781 static const std::pair<G4int, const G4double*> Z54N78(78,pZ54N78);
782 static const G4double pZ54N80[4]={3.E-12, 500., .01, 2.5E-5};
783 static const std::pair<G4int, const G4double*> Z54N80(80,pZ54N80);
784 static const G4double pZ54N82[4]={3.E-12, 500., .01, 2.5E-4};
785 static const std::pair<G4int, const G4double*> Z54N82(82,pZ54N82);
786 static const std::pair<G4int, const G4double*> Z54[N54]={Z54N70, Z54N72, Z54N74,
787 Z54N75, Z54N76, Z54N77,
788 Z54N78, Z54N80, Z54N82};
790 static const G4int N55=1;
791 static const G4double pZ55N78[4]={3.E-12, 500., .01, 2.5E-4};
792 static const std::pair<G4int, const G4double*> Z55N78(78,pZ55N78);
793 static const std::pair<G4int, const G4double*> Z55[N55]={Z55N78};
795 static const G4int N56=7;
796 static const G4double pZ56N74[4]={3.E-12, 500., .01, 2.5E-5};
797 static const std::pair<G4int, const G4double*> Z56N74(74,pZ56N74);
798 static const G4double pZ56N76[4]={3.E-12, 500., .01, 2.5E-5};
799 static const std::pair<G4int, const G4double*> Z56N76(76,pZ56N76);
800 static const G4double pZ56N78[4]={3.E-12, 500., .01, 2.5E-5};
801 static const std::pair<G4int, const G4double*> Z56N78(78,pZ56N78);
802 static const G4double pZ56N79[4]={3.E-12, 500., .01, 2.5E-5};
803 static const std::pair<G4int, const G4double*> Z56N79(79,pZ56N79);
804 static const G4double pZ56N80[4]={3.E-12, 500., .01, 2.5E-5};
805 static const std::pair<G4int, const G4double*> Z56N80(80,pZ56N80);
806 static const G4double pZ56N81[4]={3.E-12, 500., .01, 2.5E-5};
807 static const std::pair<G4int, const G4double*> Z56N81(81,pZ56N81);
808 static const G4double pZ56N82[4]={3.E-12, 500., .01, 2.5E-4};
809 static const std::pair<G4int, const G4double*> Z56N82(82,pZ56N82);
810 static const std::pair<G4int, const G4double*> Z56[N56]={Z56N74, Z56N76, Z56N78, Z56N79,
811 Z56N80, Z56N81, Z56N82};
813 static const G4int N57=2;
814 static const G4double pZ57N81[4]={3.E-12, 500., .01, 2.5E-5};
815 static const std::pair<G4int, const G4double*> Z57N81(81,pZ57N81);
816 static const G4double pZ57N82[4]={3.E-12, 500., .01, 2.5E-4};
817 static const std::pair<G4int, const G4double*> Z57N82(82,pZ57N82);
818 static const std::pair<G4int, const G4double*> Z57[N57]={Z57N81, Z57N82};
820 static const G4int N58=4;
821 static const G4double pZ58N78[4]={3.E-12, 500., .01, 2.5E-5};
822 static const std::pair<G4int, const G4double*> Z58N78(78,pZ58N78);
823 static const G4double pZ58N80[4]={3.E-12, 500., .01, 2.5E-5};
824 static const std::pair<G4int, const G4double*> Z58N80(80,pZ58N80);
825 static const G4double pZ58N82[4]={3.E-12, 500., .01, 2.5E-5};
826 static const std::pair<G4int, const G4double*> Z58N82(82,pZ58N82);
827 static const G4double pZ58N84[4]={3.E-12, 500., .01, 2.5E-4};
828 static const std::pair<G4int, const G4double*> Z58N84(84,pZ58N84);
829 static const std::pair<G4int, const G4double*> Z58[N58]={Z58N78, Z58N80, Z58N82, Z58N84};
831 static const G4int N59=1;
832 static const G4double pZ59N82[4]={3.E-12, 500., .01, 2.5E-4};
833 static const std::pair<G4int, const G4double*> Z59N82(82,pZ59N82);
834 static const std::pair<G4int, const G4double*> Z59[N59]={Z59N82};
836 static const G4int N60=7;
837 static const G4double pZ60N82[4]={3.E-12, 500., .01, 2.5E-5};
838 static const std::pair<G4int, const G4double*> Z60N82(82,pZ60N82);
839 static const G4double pZ60N83[4]={3.E-12, 500., .01, 2.5E-5};
840 static const std::pair<G4int, const G4double*> Z60N83(83,pZ60N83);
841 static const G4double pZ60N84[4]={3.E-12, 500., .01, 2.5E-5};
842 static const std::pair<G4int, const G4double*> Z60N84(84,pZ60N84);
843 static const G4double pZ60N85[4]={3.E-12, 500., .01, 2.5E-5};
844 static const std::pair<G4int, const G4double*> Z60N85(85,pZ60N85);
845 static const G4double pZ60N86[4]={3.E-12, 500., .01, 2.5E-5};
846 static const std::pair<G4int, const G4double*> Z60N86(86,pZ60N86);
847 static const G4double pZ60N88[4]={3.E-12, 500., .01, 2.5E-5};
848 static const std::pair<G4int, const G4double*> Z60N88(88,pZ60N88);
849 static const G4double pZ60N90[4]={3.E-12, 500., .01, 2.5E-4};
850 static const std::pair<G4int, const G4double*> Z60N90(90,pZ60N90);
851 static const std::pair<G4int, const G4double*> Z60[N60]={Z60N82, Z60N83, Z60N84, Z60N85,
852 Z60N86, Z60N88, Z60N90};
854 static const G4int N61=1;
855 static const G4double pZ61N0[4]={3.E-12, 500., .01, 2.5E-4};
856 static const std::pair<G4int, const G4double*> Z61N0(0,pZ61N0);
857 static const std::pair<G4int, const G4double*> Z61[N61]={Z61N0};
859 static const G4int N62=7;
860 static const G4double pZ62N82[4]={3.E-12, 500., .01, 2.5E-5};
861 static const std::pair<G4int, const G4double*> Z62N82(82,pZ62N82);
862 static const G4double pZ62N85[4]={3.E-12, 500., .01, 2.5E-5};
863 static const std::pair<G4int, const G4double*> Z62N85(85,pZ62N85);
864 static const G4double pZ62N86[4]={3.E-12, 500., .01, 2.5E-5};
865 static const std::pair<G4int, const G4double*> Z62N86(86,pZ62N86);
866 static const G4double pZ62N87[4]={3.E-12, 500., .01, 2.5E-5};
867 static const std::pair<G4int, const G4double*> Z62N87(87,pZ62N87);
868 static const G4double pZ62N88[4]={3.E-12, 500., .01, 2.5E-5};
869 static const std::pair<G4int, const G4double*> Z62N88(88,pZ62N88);
870 static const G4double pZ62N90[4]={3.E-12, 500., .01, 2.5E-5};
871 static const std::pair<G4int, const G4double*> Z62N90(90,pZ62N90);
872 static const G4double pZ62N92[4]={3.E-12, 500., .01, 2.5E-4};
873 static const std::pair<G4int, const G4double*> Z62N92(92,pZ62N92);
874 static const std::pair<G4int, const G4double*> Z62[N62]={Z62N82, Z62N85, Z62N86, Z62N87,
875 Z62N88, Z62N90, Z62N92};
877 static const G4int N63=2;
878 static const G4double pZ63N88[4]={3.E-12, 500., .01, 2.5E-5};
879 static const std::pair<G4int, const G4double*> Z63N88(88,pZ63N88);
880 static const G4double pZ63N90[4]={3.E-12, 500., .01, 2.5E-4};
881 static const std::pair<G4int, const G4double*> Z63N90(90,pZ63N90);
882 static const std::pair<G4int, const G4double*> Z63[N63]={Z63N88, Z63N90};
884 static const G4int N64=7;
885 static const G4double pZ64N88[4]={3.E-12, 500., .01, 2.5E-5};
886 static const std::pair<G4int, const G4double*> Z64N88(88,pZ64N88);
887 static const G4double pZ64N90[4]={3.E-12, 500., .01, 2.5E-5};
888 static const std::pair<G4int, const G4double*> Z64N90(90,pZ64N90);
889 static const G4double pZ64N91[4]={3.E-12, 500., .01, 2.5E-5};
890 static const std::pair<G4int, const G4double*> Z64N91(91,pZ64N91);
891 static const G4double pZ64N92[4]={3.E-12, 500., .01, 2.5E-5};
892 static const std::pair<G4int, const G4double*> Z64N92(92,pZ64N92);
893 static const G4double pZ64N93[4]={3.E-12, 500., .01, 2.5E-5};
894 static const std::pair<G4int, const G4double*> Z64N93(93,pZ64N93);
895 static const G4double pZ64N94[4]={3.E-12, 500., .01, 2.5E-5};
896 static const std::pair<G4int, const G4double*> Z64N94(94,pZ64N94);
897 static const G4double pZ64N96[4]={3.E-12, 500., .01, 2.5E-4};
898 static const std::pair<G4int, const G4double*> Z64N96(96,pZ64N96);
899 static const std::pair<G4int, const G4double*> Z64[N64]={Z64N88, Z64N90, Z64N91, Z64N92,
900 Z64N93, Z64N94, Z64N96};
902 static const G4int N65=1;
903 static const G4double pZ65N94[4]={3.E-12, 500., .01, 2.5E-4};
904 static const std::pair<G4int, const G4double*> Z65N94(82,pZ65N94);
905 static const std::pair<G4int, const G4double*> Z65[N65]={Z65N94};
907 static const G4int N66=7;
908 static const G4double pZ66N90[4]={3.E-12, 500., .01, 2.5E-5};
909 static const std::pair<G4int, const G4double*> Z66N90(90,pZ66N90);
910 static const G4double pZ66N92[4]={3.E-12, 500., .01, 2.5E-5};
911 static const std::pair<G4int, const G4double*> Z66N92(92,pZ66N92);
912 static const G4double pZ66N94[4]={3.E-12, 500., .01, 2.5E-5};
913 static const std::pair<G4int, const G4double*> Z66N94(94,pZ66N94);
914 static const G4double pZ66N95[4]={3.E-12, 500., .01, 2.5E-5};
915 static const std::pair<G4int, const G4double*> Z66N95(95,pZ66N95);
916 static const G4double pZ66N96[4]={3.E-12, 500., .01, 2.5E-5};
917 static const std::pair<G4int, const G4double*> Z66N96(96,pZ66N96);
918 static const G4double pZ66N97[4]={3.E-12, 500., .01, 2.5E-5};
919 static const std::pair<G4int, const G4double*> Z66N97(97,pZ66N97);
920 static const G4double pZ66N98[4]={3.E-12, 500., .01, 2.5E-4};
921 static const std::pair<G4int, const G4double*> Z66N98(98,pZ66N98);
922 static const std::pair<G4int, const G4double*> Z66[N66]={Z66N90, Z66N92, Z66N94, Z66N95,
923 Z66N96, Z66N97, Z66N98};
925 static const G4int N67=1;
926 static const G4double pZ67N98[4]={3.E-12, 500., .01, 2.5E-4};
927 static const std::pair<G4int, const G4double*> Z67N98(98,pZ67N98);
928 static const std::pair<G4int, const G4double*> Z67[N67]={Z67N98};
930 static const G4int N68=6;
931 static const G4double pZ68N94[4]={3.E-12, 500., .01, 2.5E-5};
932 static const std::pair<G4int, const G4double*> Z68N94(94,pZ68N94);
933 static const G4double pZ68N96[4]={3.E-12, 500., .01, 2.5E-5};
934 static const std::pair<G4int, const G4double*> Z68N96(96,pZ68N96);
935 static const G4double pZ68N98[4]={3.E-12, 500., .01, 2.5E-5};
936 static const std::pair<G4int, const G4double*> Z68N98(98,pZ68N98);
937 static const G4double pZ68N99[4]={3.E-12, 500., .01, 2.5E-5};
938 static const std::pair<G4int, const G4double*> Z68N99(99,pZ68N99);
939 static const G4double pZ68N100[4]={3.E-12, 500., .01, 2.5E-5};
940 static const std::pair<G4int, const G4double*> Z68N100(100,pZ68N100);
941 static const G4double pZ68N102[4]={3.E-12, 500., .01, 2.5E-4};
942 static const std::pair<G4int, const G4double*> Z68N102(102,pZ68N102);
943 static const std::pair<G4int, const G4double*> Z68[N68]={Z68N94, Z68N96, Z68N98,
944 Z68N99, Z68N100, Z68N102};
946 static const G4int N69=1;
947 static const G4double pZ69N100[4]={3.E-12, 500., .01, 2.5E-4};
948 static const std::pair<G4int, const G4double*> Z69N100(100,pZ69N100);
949 static const std::pair<G4int, const G4double*> Z69[N69]={Z69N100};
951 static const G4int N70=7;
952 static const G4double pZ70N98[4]={3.E-12, 500., .01, 2.5E-5};
953 static const std::pair<G4int, const G4double*> Z70N98(98,pZ70N98);
954 static const G4double pZ70N100[4]={3.E-12, 500., .01, 2.5E-5};
955 static const std::pair<G4int, const G4double*> Z70N100(100,pZ70N100);
956 static const G4double pZ70N101[4]={3.E-12, 500., .01, 2.5E-5};
957 static const std::pair<G4int, const G4double*> Z70N101(101,pZ70N101);
958 static const G4double pZ70N102[4]={3.E-12, 500., .01, 2.5E-5};
959 static const std::pair<G4int, const G4double*> Z70N102(102,pZ70N102);
960 static const G4double pZ70N103[4]={3.E-12, 500., .01, 2.5E-5};
961 static const std::pair<G4int, const G4double*> Z70N103(103,pZ70N103);
962 static const G4double pZ70N104[4]={3.E-12, 500., .01, 2.5E-5};
963 static const std::pair<G4int, const G4double*> Z70N104(104,pZ70N104);
964 static const G4double pZ70N106[4]={3.E-12, 500., .01, 2.5E-4};
965 static const std::pair<G4int, const G4double*> Z70N106(106,pZ70N106);
966 static const std::pair<G4int, const G4double*> Z70[N70]={Z70N98, Z70N100, Z70N101,
967 Z70N102, Z70N103, Z70N104,
970 static const G4int N71=2;
971 static const G4double pZ71N104[4]={3.E-12, 500., .01, 2.5E-5};
972 static const std::pair<G4int, const G4double*> Z71N104(104,pZ71N104);
973 static const G4double pZ71N105[4]={3.E-12, 500., .01, 2.5E-4};
974 static const std::pair<G4int, const G4double*> Z71N105(105,pZ71N105);
975 static const std::pair<G4int, const G4double*> Z71[N71]={Z71N104, Z71N105};
977 static const G4int N72=6;
978 static const G4double pZ72N102[4]={3.E-12, 500., .01, 2.5E-5};
979 static const std::pair<G4int, const G4double*> Z72N102(102,pZ72N102);
980 static const G4double pZ72N104[4]={3.E-12, 500., .01, 2.5E-5};
981 static const std::pair<G4int, const G4double*> Z72N104(104,pZ72N104);
982 static const G4double pZ72N105[4]={3.E-12, 500., .01, 2.5E-5};
983 static const std::pair<G4int, const G4double*> Z72N105(105,pZ72N105);
984 static const G4double pZ72N106[4]={3.E-12, 500., .01, 2.5E-5};
985 static const std::pair<G4int, const G4double*> Z72N106(106,pZ72N106);
986 static const G4double pZ72N107[4]={3.E-12, 500., .01, 2.5E-5};
987 static const std::pair<G4int, const G4double*> Z72N107(107,pZ72N107);
988 static const G4double pZ72N108[4]={3.E-12, 500., .01, 2.5E-4};
989 static const std::pair<G4int, const G4double*> Z72N108(108,pZ72N108);
990 static const std::pair<G4int, const G4double*> Z72[N72]={Z72N102, Z72N104, Z72N105,
991 Z72N106, Z72N107, Z72N108};
993 static const G4int N73=1;
994 static const G4double pZ73N108[4]={4.E-12, 1100., .027, 1.E-3};
995 static const std::pair<G4int, const G4double*> Z73N108(108,pZ73N108);
996 static const std::pair<G4int, const G4double*> Z73[N73]={Z73N108};
998 static const G4int N74=5;
999 static const G4double pZ74N106[4]={7.E-12, 1000., .03, 2.E-4};
1000 static const std::pair<G4int, const G4double*> Z74N106(106,pZ74N106);
1001 static const G4double pZ74N108[4]={7.E-12, 1300., .03, 1.5E-4};
1002 static const std::pair<G4int, const G4double*> Z74N108(108,pZ74N108);
1003 static const G4double pZ74N109[4]={2.E-12, 1700., .023, 2.E-4};
1004 static const std::pair<G4int, const G4double*> Z74N109(109,pZ74N109);
1005 static const G4double pZ74N110[4]={7.E-12, 1100., .03, 1.5E-4};
1006 static const std::pair<G4int, const G4double*> Z74N110(110,pZ74N110);
1007 static const G4double pZ74N112[4]={7.E-12, 1100., .03, 1.5E-4};
1008 static const std::pair<G4int, const G4double*> Z74N112(112,pZ74N112);
1009 static const std::pair<G4int, const G4double*> Z74[N74]={Z74N106, Z74N108, Z74N109,
1012 static const G4int N75=2;
1013 static const G4double pZ75N110[4]={5.E-12, 1000., .025, 3.E-4};
1014 static const std::pair<G4int, const G4double*> Z75N110(110,pZ75N110);
1015 static const G4double pZ75N112[4]={5.E-12, 1000., .025, 3.E-4};
1016 static const std::pair<G4int, const G4double*> Z75N112(112,pZ75N112);
1017 static const std::pair<G4int, const G4double*> Z75[N75]={Z75N110, Z75N112};
1019 static const G4int N76=7;
1020 static const G4double pZ76N108[4]={3.E-12, 500., .01, 2.5E-5};
1021 static const std::pair<G4int, const G4double*> Z76N108(108,pZ76N108);
1022 static const G4double pZ76N110[4]={3.E-12, 500., .01, 2.5E-5};
1023 static const std::pair<G4int, const G4double*> Z76N110(110,pZ76N110);
1024 static const G4double pZ76N111[4]={3.E-12, 500., .01, 2.5E-5};
1025 static const std::pair<G4int, const G4double*> Z76N111(111,pZ76N111);
1026 static const G4double pZ76N112[4]={3.E-12, 500., .01, 2.5E-5};
1027 static const std::pair<G4int, const G4double*> Z76N112(112,pZ76N112);
1028 static const G4double pZ76N113[4]={3.E-12, 500., .01, 2.5E-5};
1029 static const std::pair<G4int, const G4double*> Z76N113(113,pZ76N113);
1030 static const G4double pZ76N114[4]={3.E-12, 500., .01, 2.5E-5};
1031 static const std::pair<G4int, const G4double*> Z76N114(114,pZ76N114);
1032 static const G4double pZ76N116[4]={3.E-12, 500., .01, 2.5E-4};
1033 static const std::pair<G4int, const G4double*> Z76N116(116,pZ76N116);
1034 static const std::pair<G4int, const G4double*> Z76[N76]={Z76N108, Z76N110, Z76N111,
1035 Z76N112, Z76N113, Z76N114,
1038 static const G4int N77=2;
1039 static const G4double pZ77N114[4]={4.E-12, 1700., .028, 2.E-4};
1040 static const std::pair<G4int, const G4double*> Z77N114(114,pZ77N114);
1041 static const G4double pZ77N116[4]={5.E-12, 1500., .028, 2.E-4};
1042 static const std::pair<G4int, const G4double*> Z77N116(116,pZ77N116);
1043 static const std::pair<G4int, const G4double*> Z77[N77]={Z77N114, Z77N116};
1045 static const G4int N78=6;
1046 static const G4double pZ78N112[4]={3.E-12, 500., .01, 2.5E-5};
1047 static const std::pair<G4int, const G4double*> Z78N112(112,pZ78N112);
1048 static const G4double pZ78N114[4]={3.E-12, 500., .01, 2.5E-5};
1049 static const std::pair<G4int, const G4double*> Z78N114(114,pZ78N114);
1050 static const G4double pZ78N116[4]={3.E-12, 500., .01, 2.5E-5};
1051 static const std::pair<G4int, const G4double*> Z78N116(116,pZ78N116);
1052 static const G4double pZ78N117[4]={3.E-12, 500., .01, 2.5E-5};
1053 static const std::pair<G4int, const G4double*> Z78N117(117,pZ78N117);
1054 static const G4double pZ78N118[4]={3.E-12, 500., .01, 2.5E-5};
1055 static const std::pair<G4int, const G4double*> Z78N118(118,pZ78N118);
1056 static const G4double pZ78N120[4]={3.E-12, 500., .01, 2.5E-4};
1057 static const std::pair<G4int, const G4double*> Z78N120(120,pZ78N120);
1058 static const std::pair<G4int, const G4double*> Z78[N78]={Z78N112, Z78N114, Z78N116,
1059 Z78N117, Z78N118, Z78N120};
1061 static const G4int N79=1;
1062 static const G4double pZ79N118[4]={.2E-9, 1600., .043, 5.E-4};
1063 static const std::pair<G4int, const G4double*> Z79N118(118,pZ79N118);
1064 static const std::pair<G4int, const G4double*> Z79[N79]={Z79N118};
1066 static const G4int N80=7;
1067 static const G4double pZ80N116[4]={6.E-8, 2500., .085, 2.E-3};
1068 static const std::pair<G4int, const G4double*> Z80N116(116,pZ80N116);
1069 static const G4double pZ80N118[4]={6.E-8, 2500., .083, 1.7E-3};
1070 static const std::pair<G4int, const G4double*> Z80N118(118,pZ80N118);
1071 static const G4double pZ80N119[4]={6.E-8, 2600., .073, 2.5E-3};
1072 static const std::pair<G4int, const G4double*> Z80N119(119,pZ80N119);
1073 static const G4double pZ80N120[4]={6.E-8, 2500., .084, 1.7E-3};
1074 static const std::pair<G4int, const G4double*> Z80N120(120,pZ80N120);
1075 static const G4double pZ80N121[4]={1.5E-7, 2600., .078, 4.E-3};
1076 static const std::pair<G4int, const G4double*> Z80N121(121,pZ80N121);
1077 static const G4double pZ80N122[4]={6.E-8, 2500., .083, 1.6E-3};
1078 static const std::pair<G4int, const G4double*> Z80N122(122,pZ80N122);
1079 static const G4double pZ80N124[4]={6.E-8, 2500., .083, 1.5E-3};
1080 static const std::pair<G4int, const G4double*> Z80N124(124,pZ80N124);
1081 static const std::pair<G4int, const G4double*> Z80[N80]={Z80N116, Z80N118, Z80N119,
1082 Z80N120, Z80N121, Z80N122,
1085 static const G4int N81=2;
1086 static const G4double pZ81N122[4]={3.E-12, 500., .01, 2.5E-5};
1087 static const std::pair<G4int, const G4double*> Z81N122(122,pZ81N122);
1088 static const G4double pZ81N124[4]={3.E-12, 500., .01, 2.5E-4};
1089 static const std::pair<G4int, const G4double*> Z81N124(124,pZ81N124);
1090 static const std::pair<G4int, const G4double*> Z81[N81]={Z81N122, Z81N124};
1092 static const G4int N82=4;
1093 static const G4double pZ82N122[4]={.2E-9, 40., .002, 6.E-4};
1094 static const std::pair<G4int, const G4double*> Z82N122(122,pZ82N122);
1095 static const G4double pZ82N124[4]={6.E-9, 1700., .076, 7.E-4};
1096 static const std::pair<G4int, const G4double*> Z82N124(124,pZ82N124);
1097 static const G4double pZ82N125[4]={.2E-9, 770., .057, 4.5E-4};
1098 static const std::pair<G4int, const G4double*> Z82N125(125,pZ82N125);
1099 static const G4double pZ82N126[4]={4.E-9, 0., .051, 2.E-4};
1100 static const std::pair<G4int, const G4double*> Z82N126(126,pZ82N126);
1101 static const std::pair<G4int, const G4double*> Z82[N82]={Z82N122, Z82N124, Z82N125,
1104 static const G4int N83=1;
1105 static const G4double pZ83N126[4]={1.5E-9, 150., .052, 5.E-5};
1106 static const std::pair<G4int, const G4double*> Z83N126(126,pZ83N126);
1107 static const std::pair<G4int, const G4double*> Z83[N83]={Z83N126};
1109 static const G4int N84=1;
1110 static const G4double pZ84N0[4]={3.E-12, 500., .01, 2.5E-4};
1111 static const std::pair<G4int, const G4double*> Z84N0(0,pZ84N0);
1112 static const std::pair<G4int, const G4double*> Z84[N84]={Z84N0};
1114 static const G4int N85=1;
1115 static const G4double pZ85N0[4]={3.E-12, 500., .01, 2.5E-4};
1116 static const std::pair<G4int, const G4double*> Z85N0(0,pZ85N0);
1117 static const std::pair<G4int, const G4double*> Z85[N85]={Z85N0};
1119 static const G4int N86=1;
1120 static const G4double pZ86N0[4]={3.E-12, 500., .01, 2.5E-4};
1121 static const std::pair<G4int, const G4double*> Z86N0(0,pZ86N0);
1122 static const std::pair<G4int, const G4double*> Z86[N86]={Z86N0};
1124 static const G4int N87=1;
1125 static const G4double pZ87N0[4]={3.E-12, 500., .01, 2.5E-4};
1126 static const std::pair<G4int, const G4double*> Z87N0(0,pZ87N0);
1127 static const std::pair<G4int, const G4double*> Z87[N87]={Z87N0};
1129 static const G4int N88=1;
1130 static const G4double pZ88N138[4]={3.E-9, 2200., .057, 1.2E-3};
1131 static const std::pair<G4int, const G4double*> Z88N138(138,pZ88N138);
1132 static const std::pair<G4int, const G4double*> Z88[N88]={Z88N138};
1134 static const G4int N89=1;
1135 static const G4double pZ89N0[4]={3.E-12, 500., .01, 2.5E-4};
1136 static const std::pair<G4int, const G4double*> Z89N0(0,pZ89N0);
1137 static const std::pair<G4int, const G4double*> Z89[N89]={Z89N0};
1139 static const G4int N90=1;
1140 static const G4double pZ90N142[4]={1.E-11, 1200., .028, 3.E-4};
1141 static const std::pair<G4int, const G4double*> Z90N142(142,pZ90N142);
1142 static const std::pair<G4int, const G4double*> Z90[N90]={Z90N142};
1144 static const G4int N91=1;
1145 static const G4double pZ91N0[4]={3.E-12, 500., .01, 2.5E-4};
1146 static const std::pair<G4int, const G4double*> Z91N0(0,pZ91N0);
1147 static const std::pair<G4int, const G4double*> Z91[N91]={Z91N0};
1149 static const G4int N92=2;
1150 static const G4double pZ92N143[4]={2.E-11, 2700., .026, 6.E-4};
1151 static const std::pair<G4int, const G4double*> Z92N143(143,pZ92N143);
1152 static const G4double pZ92N146[4]={1.E-11, 1700., .029, 2.5E-4};
1153 static const std::pair<G4int, const G4double*> Z92N146(146,pZ92N146);
1154 static const std::pair<G4int, const G4double*> Z92[N92]={Z92N143, Z92N146};
1156 static const G4int N93=1;
1157 static const G4double pZ93N144[4]={4.E-8, 3700., .066, 3.5E-3};
1158 static const std::pair<G4int, const G4double*> Z93N144(144,pZ93N144);
1159 static const std::pair<G4int, const G4double*> Z93[N93]={Z93N144};
1161 static const G4int N94=3;
1162 static const G4double pZ94N145[4]={8.E-11, 2900., .029, 1.3E-3};
1163 static const std::pair<G4int, const G4double*> Z94N145(145,pZ94N145);
1164 static const G4double pZ94N148[4]={9.E-12, 1400., .025, 3.E-4};
1165 static const std::pair<G4int, const G4double*> Z94N148(148,pZ94N148);
1166 static const G4double pZ94N150[4]={4.E-12, 1500., .023, 1.2E-4};
1167 static const std::pair<G4int, const G4double*> Z94N150(150,pZ94N150);
1168 static const std::pair<G4int, const G4double*> Z94[N94]={Z94N145, Z94N148, Z94N150};
1170 static const G4int N95=1;
1171 static const G4double pZ95N0[4]={3.E-12, 500., .01, 2.5E-4};
1172 static const std::pair<G4int, const G4double*> Z95N0(0,pZ95N0);
1173 static const std::pair<G4int, const G4double*> Z95[N95]={Z95N0};
1175 static const G4int N96=1;
1176 static const G4double pZ96N151[4]={1.5E-8, 3700., .055, 2.E-3};
1177 static const std::pair<G4int, const G4double*> Z96N151(151,pZ96N151);
1178 static const std::pair<G4int, const G4double*> Z96[N96]={Z96N151};
1180 static const G4int NZ=97;
1181 static const std::pair<G4int, const G4double*>* Pars[NZ]={Z0,Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z8,Z9,
1182 Z10,Z11,Z12,Z13,Z14,Z15,Z16,Z17,Z18,Z19,Z20,Z21,Z22,Z23,Z24,Z25,Z26,Z27,Z28,Z29,Z30,
1183 Z31,Z32,Z33,Z34,Z35,Z36,Z37,Z38,Z39,Z40,Z41,Z42,Z43,Z44,Z45,Z46,Z47,Z48,Z49,Z50,Z51,
1184 Z52,Z53,Z54,Z55,Z56,Z57,Z58,Z59,Z60,Z61,Z62,Z63,Z64,Z65,Z66,Z67,Z68,Z69,Z70,Z71,Z72,
1185 Z73,Z74,Z75,Z76,Z77,Z78,Z79,Z80,Z81,Z82,Z83,Z84,Z85,Z86,Z87,Z88,Z89,Z90,Z91,Z92,Z93,
1187 static const G4int NIso[NZ]={N0,N1,N2,N3,N4,N5,N6,N7,N8,N9,N10,N11,N12,N13,N14,N15,N16,
1188 N17,N18,N19,N20,N21,N22,N23,N24,N25,N26,N27,N28,N29,N30,N31,N32,N33,N34,N35,N36,N37,
1189 N38,N39,N40,N41,N42,N43,N44,N45,N46,N47,N48,N49,N50,N51,N52,N53,N54,N55,N56,N57,N58,
1190 N59,N60,N61,N62,N63,N64,N65,N66,N67,N68,N69,N70,N71,N72,N73,N74,N75,N76,N77,N78,N79,
1191 N80,N81,N82,N83,N84,N85,N86,N87,N88,N89,N90,N91,N92,N93,N94,N95,N96};
1198 else if(tZ<97 && tN<152)
1206 if(nn)
for (
G4int in=0; in<
nn; in++)
1208 std::pair<G4int, const G4double*> curIs=Pars[tZ][in];
1209 if(curIs.first == tN)
1220 if(nfound)
G4cout<<
"-Warning-G4ChipsNeutronInelasticXS::CSLin: Z="<<tZ<<
", N="
1221 <<tN<<
" isotope is not implemented in CHIPS"<<
G4endl;
1226 sigma+=pex*std::exp(-dp*dp/wid);
1231 G4cerr<<
"-Warning-G4ChipsNeutronNuclearCroSect::CSLin:*Bad A* Z="<<tZ<<
", N="<<tN<<
G4endl;
1234 if(sigma<0.)
return 0.;
1256 El=1./(0.00012+p2*(0.051+0.1*
p2));
1273 El=LE+(0.0557*lp2+6.72+30./P)/(1.+0.49*rp2/P);
1274 To=LE+(0.3 *lp2+38.2)/(1.+0.54*rp2*rp2);
1287 else if(tZ<97 && tN<152)
1306 G4double c=(170.+3600./sa2)/(1.+65./sa2);
1309 G4double r=.21+.62*dl2/(1.+.5*dl2);
1310 G4double gg=42.*(std::exp(al*0.8)+4.E-8*
a4)/(1.+28./a)/(1.+5.E-5*
a2);
1311 G4double e=5.*((a6+.021*a8)/(1.+.0013*a7)+.001*
a3)/(1.+.0007*a2);
1317 sigma=(c+d*d)/(1+r/p4)+(gg+e*std::exp(-ss*P))/(1+h/p4/p4);
1321 G4cerr<<
"-Warning-G4ChipsNeutronNuclearCroSect::CSForm:*Bad A* Z="<<tZ<<
", N="<<tN<<
G4endl;
1324 if(sigma<0.)
return 0.;
1332 G4cerr<<
"***G4ChipsNeutronInelasticXS::EquLinearFit: DX="<<DX<<
", N="<<N<<
G4endl;
1338 G4int jj=
static_cast<int>(d);
1340 else if(jj>N2) jj=N2;
std::vector< G4double > colCS
~G4ChipsNeutronInelasticXS()
std::vector< G4int > colZ
std::vector< G4double * > * HEN
G4_DECLARE_XS_FACTORY(G4ChipsNeutronInelasticXS)
G4double GetTotalMomentum() const
G4GLOB_DLL std::ostream G4cout
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
std::vector< G4int > colN
virtual G4bool IsIsoApplicable(const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
static const G4double A[nN]
G4double EquLinearFit(G4double X, G4int N, G4double X0, G4double DX, G4double *Y)
G4ChipsNeutronInelasticXS()
static const double millibarn
std::vector< G4double > colP
std::vector< G4double > colTH
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
G4double CrossSectionFormula(G4int targZ, G4int targN, G4double P, G4double lP)
G4double CrossSectionLog(G4int targZ, G4int targN, G4double lP)
static const G4double pos
static const G4double THmin
G4double CrossSectionLin(G4int targZ, G4int targN, G4double P)
G4GLOB_DLL std::ostream G4cerr
G4double CalculateCrossSection(G4int F, G4int I, G4int PDG, G4int Z, G4int N, G4double Momentum)
std::vector< G4double * > * LEN