97 std::vector<G4double*>::iterator
pos;
98 for (pos=
CST.begin(); pos<
CST.end(); pos++)
101 for (pos=
PAR.begin(); pos<
PAR.end(); pos++)
104 for (pos=
SST.begin(); pos<
SST.end(); pos++)
107 for (pos=
S1T.begin(); pos<
S1T.end(); pos++)
110 for (pos=
B1T.begin(); pos<
B1T.end(); pos++)
113 for (pos=
S2T.begin(); pos<
S2T.end(); pos++)
116 for (pos=
B2T.begin(); pos<
B2T.end(); pos++)
119 for (pos=
S3T.begin(); pos<
S3T.end(); pos++)
122 for (pos=
B3T.begin(); pos<
B3T.end(); pos++)
125 for (pos=
S4T.begin(); pos<
S4T.end(); pos++)
128 for (pos=
B4T.begin(); pos<
B4T.end(); pos++)
157 static G4ThreadLocal std::vector <G4int> *colN_G4MT_TLS_ = 0 ;
if (!colN_G4MT_TLS_) colN_G4MT_TLS_ =
new std::vector <G4int> ; std::vector <G4int> &colN = *colN_G4MT_TLS_;
158 static G4ThreadLocal std::vector <G4int> *colZ_G4MT_TLS_ = 0 ;
if (!colZ_G4MT_TLS_) colZ_G4MT_TLS_ =
new std::vector <G4int> ; std::vector <G4int> &colZ = *colZ_G4MT_TLS_;
159 static G4ThreadLocal std::vector <G4double> *colP_G4MT_TLS_ = 0 ;
if (!colP_G4MT_TLS_) colP_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &colP = *colP_G4MT_TLS_;
160 static G4ThreadLocal std::vector <G4double> *colTH_G4MT_TLS_ = 0 ;
if (!colTH_G4MT_TLS_) colTH_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &colTH = *colTH_G4MT_TLS_;
161 static G4ThreadLocal std::vector <G4double> *colCS_G4MT_TLS_ = 0 ;
if (!colCS_G4MT_TLS_) colCS_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &colCS = *colCS_G4MT_TLS_;
174 if(colN[i]==tgN && colZ[i]==tgZ)
193 if(lastCS<=0. && pEn>
lastTH)
214 colP.push_back(pMom);
233 static G4ThreadLocal std::vector <G4double> *PIN_G4MT_TLS_ = 0 ;
if (!PIN_G4MT_TLS_) PIN_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &PIN = *PIN_G4MT_TLS_;
300 G4int blast=
static_cast<int>(shift);
319 G4int blast=
static_cast<int>(shift);
360 const G4int n_npel=24;
361 const G4int n_ppel=32;
363 G4double np_el[n_npel]={12.,.05,.0001,5.,.35,6.75,.14,19.,.6,6.75,.14,13.,.14,.6,.00013,
364 75.,.001,7.2,4.32,.012,2.5,0.0,12.,.34};
367 G4double pp_el[n_ppel]={2.865,18.9,.6461,3.,9.,.425,.4276,.0022,5.,74.,3.,3.4,.2,.17,
368 .001,8.,.055,3.64,5.e-5,4000.,1500.,.46,1.2e6,3.5e6,5.e-5,1.e10,
369 8.5e8,1.e10,1.1,3.4e6,6.8e6,0.};
373 static const G4int N0=1;
374 static const G4double pZ0N1[7]={0., 0., 0., 0., 0., 0., 0.};
375 static const std::pair<G4int, const G4double*> Z0N1(1,pZ0N1);
376 static const std::pair<G4int, const G4double*> Z0[N0]={Z0N1};
378 static const G4int N1=3;
379 static const G4double pZ1N0[7]={0., 0., 0., 0., 0., 0., 0.};
380 static const std::pair<G4int, const G4double*> Z1N0(0,pZ1N0);
381 static const G4double pZ1N1[7]={6.E-5, 4., .055, 1.1E-8, .008, 1.2E-8, .019};
382 static const std::pair<G4int, const G4double*> Z1N1(1,pZ1N1);
383 static const G4double pZ1N2[7]={6.E-5, 2.2, .051, 1.E-8, .04, 9.E-8, .0075};
384 static const std::pair<G4int, const G4double*> Z1N2(2,pZ1N2);
385 static const std::pair<G4int, const G4double*> Z1[N1]={Z1N0, Z1N1, Z1N2};
387 static const G4int N2=2;
388 static const G4double pZ2N1[7]={6.E-5, 3., .06, 4.E-9, .03, 7.E-8, .015};
389 static const std::pair<G4int, const G4double*> Z2N1(1,pZ2N1);
390 static const G4double pZ2N2[7]={3.E-4, .23, 1., 1.5E-9, 2.E-02, 1.E-8, .003};
391 static const std::pair<G4int, const G4double*> Z2N2(2,pZ2N2);
392 static const std::pair<G4int, const G4double*> Z2[N2]={Z2N1, Z2N2};
394 static const G4int N3=2;
395 static const G4double pZ3N3[7]={3.1E-7, 1.7, 1.3E-4, 1.E-8, .02, 1.1E-7, .0023};
396 static const std::pair<G4int, const G4double*> Z3N1(3,pZ3N3);
397 static const G4double pZ3N4[7]={1.3E-6, 1.8, 7.6E-4, 9.E-9, .03, 1.E-7, .0029};
398 static const std::pair<G4int, const G4double*> Z3N2(4,pZ3N4);
399 static const std::pair<G4int, const G4double*> Z3[N3]={Z3N1, Z3N2};
401 static const G4int N4=2;
402 static const G4double pZ4N3[7]={2.E-4, 1.4, 2.7, 0., .02, 5.E-8, 0.};
403 static const std::pair<G4int, const G4double*> Z4N3(3,pZ4N3);
404 static const G4double pZ4N5[7]={1.E-6, 5.7, .0011, 3.E-9, .007, 2.E-8, .016};
405 static const std::pair<G4int, const G4double*> Z4N5(5,pZ4N5);
406 static const std::pair<G4int, const G4double*> Z4[N4]={Z4N3,Z4N5};
408 static const G4int N5=2;
409 static const G4double pZ5N5[7]={8.E-7, 5., 3.4E-4, 7.E-9, 1.E-02, 1.E-07, .0053};
410 static const std::pair<G4int, const G4double*> Z5N5(5,pZ5N5);
411 static const G4double pZ5N6[7]={4.8E-6, 6.6, .0035, 4.E-9, .003, 1.E-8, .012};
412 static const std::pair<G4int, const G4double*> Z5N6(6,pZ5N6);
413 static const std::pair<G4int, const G4double*> Z5[N5]={Z5N5, Z5N6};
415 static const G4int N6=2;
416 static const G4double pZ6N6[7]={4.9E-6, 6.6, .0035, 4.E-9, .002, 6.E-9, .011};
417 static const std::pair<G4int, const G4double*> Z6N6(6,pZ6N6);
418 static const G4double pZ6N7[7]={4.9E-6, 6.6, .0035, 4.E-9, .002, 6.E-9, .011};
419 static const std::pair<G4int, const G4double*> Z6N7(7,pZ6N7);
420 static const std::pair<G4int, const G4double*> Z6[N6]={Z6N6, Z6N7};
422 static const G4int N7=2;
423 static const G4double pZ7N7[7]={4.9E-6, 1.6, .03, .4E-9, .02, 6.E-8, .021};
424 static const std::pair<G4int, const G4double*> Z7N7(7,pZ7N7);
425 static const G4double pZ7N8[7]={2.5E-6, 5., .0021, 2.5E-9, .015, 5.E-8, .009};
426 static const std::pair<G4int, const G4double*> Z7N8(8,pZ7N8);
427 static const std::pair<G4int, const G4double*> Z7[N7]={Z7N7, Z7N8};
429 static const G4int N8=3;
430 static const G4double pZ8N8[7]={2.5E-6, 5.3, .0018, 3.E-9, .01, 1.5E-8, .0075};
431 static const std::pair<G4int, const G4double*> Z8N8(8,pZ8N8);
432 static const G4double pZ8N9[7]={1.4E-6, 2.1, .0025, 1.3E-9, .02, 1.3E-7, .0072};
433 static const std::pair<G4int, const G4double*> Z8N9(9,pZ8N9);
434 static const G4double pZ8N10[7]={1.4E-6, 2.1, .0025, 1.3E-9, .02, 1.3E-7, .0072};
435 static const std::pair<G4int, const G4double*> Z8N10(10,pZ8N10);
436 static const std::pair<G4int, const G4double*> Z8[N8]={Z8N8, Z8N9, Z8N10};
438 static const G4int N9=1;
439 static const G4double pZ9N10[7]={1.4E-6, 7.5, 6.7E-4, 4.E-9, 2.E-5, 7.E-12, .0066};
440 static const std::pair<G4int, const G4double*> Z9N10(10,pZ9N10);
441 static const std::pair<G4int, const G4double*> Z9[N9]={Z9N10};
443 static const G4int N10=3;
444 static const G4double pZ10N10[7]={1.4E-5, 7., .01, 0., 3.E-11, 7.E-24, .12};
445 static const std::pair<G4int, const G4double*> Z10N10(10,pZ10N10);
446 static const G4double pZ10N11[7]={1.4E-6, 7.5, 6.7E-4, 4.E-9, 2.E-5, 7.E-12, .0066};
447 static const std::pair<G4int, const G4double*> Z10N11(11,pZ10N11);
448 static const G4double pZ10N12[7]={1.4E-5, 7., .01, 0., 3.E-11, 7.E-24, .12};
449 static const std::pair<G4int, const G4double*> Z10N12(12,pZ10N12);
450 static const std::pair<G4int, const G4double*> Z10[N10]={Z10N10, Z10N11, Z10N12};
452 static const G4int N11=2;
453 static const G4double pZ11N11[7]={1.4E-5, 7., .01, 0., 3.E-11, 7.E-24, .12};
454 static const std::pair<G4int, const G4double*> Z11N11(11,pZ11N11);
455 static const G4double pZ11N12[7]={1.4E-6, 7.6, 6.E-4, 5.E-9, 7.E-9, 3.E-18, .0056};
456 static const std::pair<G4int, const G4double*> Z11N12(12,pZ11N12);
457 static const std::pair<G4int, const G4double*> Z11[N11]={Z11N11, Z11N12};
459 static const G4int N12=3;
460 static const G4double pZ12N12[7]={8.E-7, 3., .001, 1.8E-9, .0015, .2E-9, .006};
461 static const std::pair<G4int, const G4double*> Z12N12(12,pZ12N12);
462 static const G4double pZ12N13[7]={8.E-7, 7., 3.E-4, 6.E-9, .006, 4.E-8, .0042};
463 static const std::pair<G4int, const G4double*> Z12N13(13,pZ12N13);
464 static const G4double pZ12N14[7]={1.2E-6, 6.8, 5.E-4, 5.E-9, .007, 2.E-8, .0044};
465 static const std::pair<G4int, const G4double*> Z12N14(14,pZ12N14);
466 static const std::pair<G4int, const G4double*> Z12[N12]={Z12N12, Z12N13, Z12N14};
468 static const G4int N13=1;
469 static const G4double pZ13N14[7]={3.E-7, 5., 8.4E-5, 7.E-9, .008, 2.E-8, .0022};
470 static const std::pair<G4int, const G4double*> Z13N14(14,pZ13N14);
471 static const std::pair<G4int, const G4double*> Z13[N13]={Z13N14};
473 static const G4int N14=3;
474 static const G4double pZ14N14[7]={1.2E-6, 6., 4.E-4, 6.E-9, .012, 8.E-8, .0029};
475 static const std::pair<G4int, const G4double*> Z14N14(14,pZ14N14);
476 static const G4double pZ14N15[7]={2.4E-6, 4., .0016, 3.E-9, .018, 6.E-8, .0037};
477 static const std::pair<G4int, const G4double*> Z14N15(15,pZ14N15);
478 static const G4double pZ14N16[7]={6.E-7, 4., 3.7E-4, 3.E-9, .018, 6.E-8, .0036};
479 static const std::pair<G4int, const G4double*> Z14N16(16,pZ14N16);
480 static const std::pair<G4int, const G4double*> Z14[N14]={Z14N14, Z14N15, Z14N16};
482 static const G4int N15=1;
483 static const G4double pZ15N16[7]={6.E-7, 3., 8.2E-4, 1.4E-9, .03, 8.E-8, .0059};
484 static const std::pair<G4int, const G4double*> Z15N16(16,pZ15N16);
485 static const std::pair<G4int, const G4double*> Z15[N15]={Z15N16};
487 static const G4int N16=4;
488 static const G4double pZ16N16[7]={6.E-7, 3., 1.9E-4, 5.E-9, .03, 6.E-8, .0013};
489 static const std::pair<G4int, const G4double*> Z16N16(16,pZ16N16);
490 static const G4double pZ16N17[7]={2.4E-6, 3., .0023, 2.E-9, .03, 6.5E-8, .004};
491 static const std::pair<G4int, const G4double*> Z16N17(17,pZ16N17);
492 static const G4double pZ16N18[7]={2.4E-6, 1.6, .0031, 1.4E-9, .03, 4.E-08, .0028};
493 static const std::pair<G4int, const G4double*> Z16N18(18,pZ16N18);
494 static const G4double pZ16N20[7]={2.4E-6, 3.1, .0017, 2.5E-9, .03, 5.E-08, .0029};
495 static const std::pair<G4int, const G4double*> Z16N20(20,pZ16N20);
496 static const std::pair<G4int, const G4double*> Z16[N16]={Z16N16, Z16N17, Z16N18, Z16N20};
498 static const G4int N17=2;
499 static const G4double pZ17N18[7]={1.2E-7, .04, .062, 3.E-12, 3.E-02, 3.E-08, .027};
500 static const std::pair<G4int, const G4double*> Z17N18(18,pZ17N18);
501 static const G4double pZ17N20[7]={1.2E-7, 2., 6.8E-5, 2.7E-9, .03, 4.E-8, .0015};
502 static const std::pair<G4int, const G4double*> Z17N20(20,pZ17N20);
503 static const std::pair<G4int, const G4double*> Z17[N17]={Z17N18, Z17N20};
505 static const G4int N18=3;
506 static const G4double pZ18N18[7]={1.2E-7, .52, .017, 1.1E-11, .03, 3.E-8, .095};
507 static const std::pair<G4int, const G4double*> Z18N18(18,pZ18N18);
508 static const G4double pZ18N20[7]={1.2E-07, .09, .012, 1.8E-11, .03, 3.E-8, .011};
509 static const std::pair<G4int, const G4double*> Z18N20(20,pZ18N20);
510 static const G4double pZ18N22[7]={1.2E-7, .65, 1.2E-4, 1.5E-9, .03, 5.E-8, 8.E-4};
511 static const std::pair<G4int, const G4double*> Z18N22(22,pZ18N22);
512 static const std::pair<G4int, const G4double*> Z18[N18]={Z18N18, Z18N20, Z18N22};
514 static const G4int N19=3;
515 static const G4double pZ19N20[7]={1.2E-7, 1.3, 1.9E-4, .9E-9, .04, 5.5E-8, .0026};
516 static const std::pair<G4int, const G4double*> Z19N20(20,pZ19N20);
517 static const G4double pZ19N21[7]={1.6E-7, 1.2, 3.7E-4, .8E-9, .04, 6.5E-8, .0034};
518 static const std::pair<G4int, const G4double*> Z19N21(21,pZ19N21);
519 static const G4double pZ19N22[7]={6.E-8, 1.3, 1.2E-4, .9E-9, .04, 6.E-8, .0031};
520 static const std::pair<G4int, const G4double*> Z19N22(22,pZ19N22);
521 static const std::pair<G4int, const G4double*> Z19[N19]={Z19N20, Z19N21, Z19N22};
523 static const G4int N20=6;
524 static const G4double pZ20N20[7]={2.4E-7, 3.4, 2.1E-4, 1.5E-9, .035, 6.E-8, .0037};
525 static const std::pair<G4int, const G4double*> Z20N20(20,pZ20N20);
526 static const G4double pZ20N22[7]={6.E-8, 2.7, 2.7E-5, 3.E-9, .035, 6.E-8, .0014};
527 static const std::pair<G4int, const G4double*> Z20N22(22,pZ20N22);
528 static const G4double pZ20N23[7]={1.5E-8, 1.8, 3.4E-5, .6E-9, .04, 6.E-8, .0049};
529 static const std::pair<G4int, const G4double*> Z20N23(23,pZ20N23);
530 static const G4double pZ20N24[7]={3.E-6, 5., .002, 2.E-9, .03, 7.E-8, .0038};
531 static const std::pair<G4int, const G4double*> Z20N24(24,pZ20N24);
532 static const G4double pZ20N26[7]={1.7E-5, 18., .0027, 1.E-8, 2.E-7, 7.E-17, .0047};
533 static const std::pair<G4int, const G4double*> Z20N26(26,pZ20N26);
534 static const G4double pZ20N28[7]={7.6E-6, .4, .07, .13E-9, .05, 4.E-8, .0042};
535 static const std::pair<G4int, const G4double*> Z20N28(28,pZ20N28);
536 static const std::pair<G4int, const G4double*> Z20[N20]={Z20N20, Z20N22, Z20N23,
537 Z20N24, Z20N26, Z20N28};
539 static const G4int N21=1;
540 static const G4double pZ21N24[7]={3.6E-9, 1.5, 5.2E-5, .1E-9, .05, 1.E-7, .025};
541 static const std::pair<G4int, const G4double*> Z21N24(24,pZ21N24);
542 static const std::pair<G4int, const G4double*> Z21[N21]={Z21N24};
544 static const G4int N22=5;
545 static const G4double pZ22N24[7]={2.8E-8, 1.8, 5.6E-5, .6E-9, .05, 8.E-8, .0042};
546 static const std::pair<G4int, const G4double*> Z22N24(24,pZ22N24);
547 static const G4double pZ22N25[7]={3.1E-9, 1.6, 6.E-6, .8E-9, .04, 8.E-8, .0036};
548 static const std::pair<G4int, const G4double*> Z22N25(25,pZ22N25);
549 static const G4double pZ22N26[7]={3.E-9, 4., 3.2E-6, 1.4E-9, .05, 2.E-7, .0048};
550 static const std::pair<G4int, const G4double*> Z22N26(26,pZ22N26);
551 static const G4double pZ22N27[7]={1.E-8, 2., 3.4E-6, 4.5E-9, .05, 8.E-8, 7.7E-4};
552 static const std::pair<G4int, const G4double*> Z22N27(27,pZ22N27);
553 static const G4double pZ22N28[7]={4.E-7, 4., 3.7E-4, 1.E-09, .05, 1.E-7, .0041};
554 static const std::pair<G4int, const G4double*> Z22N28(28,pZ22N28);
555 static const std::pair<G4int, const G4double*> Z22[N22]={Z22N24, Z22N25, Z22N26,
558 static const G4int N23=2;
559 static const G4double pZ23N27[7]={.3E-9, 2., 8.E-7, .55E-9, .07, 1.7E-7, .0055};
560 static const std::pair<G4int, const G4double*> Z23N27(27,pZ23N27);
561 static const G4double pZ23N28[7]={.3E-9, 2., 8.E-7, .55E-9, .07, 1.7E-7, .0055};
562 static const std::pair<G4int, const G4double*> Z23N28(28,pZ23N28);
563 static const std::pair<G4int, const G4double*> Z23[N23]={Z23N27, Z23N28};
565 static const G4int N24=4;
566 static const G4double pZ24N26[7]={1.2E-9, 2.8, 1.E-6, 1.7E-9, .07, 1.7E-7, .0026};
567 static const std::pair<G4int, const G4double*> Z24N26(26,pZ24N26);
568 static const G4double pZ24N28[7]={4.4E-6, 11., .0012, 5.E-9, .04, 3.E-7, .0032};
569 static const std::pair<G4int, const G4double*> Z24N28(28,pZ24N28);
570 static const G4double pZ24N29[7]={1.8E-9, 2.4, 6.3E-6, .5E-9, .07, 2.E-7, .0085};
571 static const std::pair<G4int, const G4double*> Z24N29(29,pZ24N29);
572 static const G4double pZ24N30[7]={4.8E-8, 2.8, 4.4E-5, 1.4E-9, .07, 2.E-7, .0027};
573 static const std::pair<G4int, const G4double*> Z24N30(30,pZ24N30);
574 static const std::pair<G4int, const G4double*> Z24[N24]={Z24N26, Z24N28, Z24N29, Z24N30};
576 static const G4int N25=1;
577 static const G4double pZ25N30[7]={6.5E-11, 1.4, 1.E-7, .8E-9, .07, 1.7E-7, .0022};
578 static const std::pair<G4int, const G4double*> Z25N30(30,pZ25N30);
579 static const std::pair<G4int, const G4double*> Z25[N25]={Z25N30};
581 static const G4int N26=4;
582 static const G4double pZ26N28[7]={3.9E-8, 5., 1.7E-5, 3.E-9, .07, 3.E-7, .0023};
583 static const std::pair<G4int, const G4double*> Z26N28(28,pZ26N28);
584 static const G4double pZ26N30[7]={5.E-9, .4, 1.5E-4, 4.E-11, .1, 3.E-7, .012};
585 static const std::pair<G4int, const G4double*> Z26N30(30,pZ26N30);
586 static const G4double pZ26N31[7]={.5E-9, .5, 2.6E-6, .3E-9, .11, 5.E-7, .0027};
587 static const std::pair<G4int, const G4double*> Z26N31(31,pZ26N31);
588 static const G4double pZ26N32[7]={1.E-7, 3.1, 1.E-4, 1.3E-9, .11, 5.E-7, .0031};
589 static const std::pair<G4int, const G4double*> Z26N32(32,pZ26N32);
590 static const std::pair<G4int, const G4double*> Z26[N26]={Z26N28, Z26N30, Z26N31, Z26N32};
592 static const G4int N27=2;
593 static const G4double pZ27N31[7]={4.E-7, 3., .004, 0., .11, 4.5E-7, .07};
594 static const std::pair<G4int, const G4double*> Z27N31(31,pZ27N31);
595 static const G4double pZ27N32[7]={4.E-7, 5., 5.E-4, 1.2E-9, .13, 6.E-7, .006};
596 static const std::pair<G4int, const G4double*> Z27N32(32,pZ27N32);
597 static const std::pair<G4int, const G4double*> Z27[N27]={Z27N32, Z27N31};
599 static const G4int N28=6;
600 static const G4double pZ28N30[7]={1.E-7, 2.5, .001, .14E-9, .13, 6.E-7, .025};
601 static const std::pair<G4int, const G4double*> Z28N30(30,pZ28N30);
602 static const G4double pZ28N31[7]={1.E-7, 19., 1.2E-5, 1.E-8, 4.E-12, 3.E-22, .0024};
603 static const std::pair<G4int, const G4double*> Z28N31(31,pZ28N31);
604 static const G4double pZ28N32[7]={1.E-8, 2.5, 3.9E-6, 3.5E-9, .13, 6.E-7, .001};
605 static const std::pair<G4int, const G4double*> Z28N32(32,pZ28N32);
606 static const G4double pZ28N33[7]={5.E-9, 2.6, 1.5E-5, .42E-9, .13, 7.E-7, .008};
607 static const std::pair<G4int, const G4double*> Z28N33(33,pZ28N33);
608 static const G4double pZ28N34[7]={.24E-9, 2., 1.2E-6, .25E-9, .13, 6.E-7, .0094};
609 static const std::pair<G4int, const G4double*> Z28N34(34,pZ28N34);
610 static const G4double pZ28N36[7]={1.E-8, 3., 5.5E-8, 2.8E-7, .12, 6.E-7, 1.6E-5};
611 static const std::pair<G4int, const G4double*> Z28N36(36,pZ28N36);
612 static const std::pair<G4int, const G4double*> Z28[N28]={Z28N30, Z28N31, Z28N32, Z28N33,
615 static const G4int N29=2;
616 static const G4double pZ29N34[7]={1.1E-7, 3.5, 1.6E-4, .9E-9, .13, 7.E-7, .005};
617 static const std::pair<G4int, const G4double*> Z29N34(34,pZ29N34);
618 static const G4double pZ29N36[7]={1.1E-7, 3.5, 4.3E-4, .3E-9, .13, 7.E-7, .013};
619 static const std::pair<G4int, const G4double*> Z29N36(36,pZ29N36);
620 static const std::pair<G4int, const G4double*> Z29[N29]={Z29N34, Z29N36};
622 static const G4int N30=5;
623 static const G4double pZ30N34[7]={1.1E-7, 4., 1.2E-4, 1.2E-9, .17, 7.E-7, .004};
624 static const std::pair<G4int, const G4double*> Z30N34(34,pZ30N34);
625 static const G4double pZ30N36[7]={1.1E-7, 4., 1.2E-4, 1.2E-9, .17, 7.E-7, .004};
626 static const std::pair<G4int, const G4double*> Z30N36(36,pZ30N36);
627 static const G4double pZ30N37[7]={1.1E-7, 4., 1.2E-4, 1.2E-9, .17, 7.E-7, .004};
628 static const std::pair<G4int, const G4double*> Z30N37(37,pZ30N37);
629 static const G4double pZ30N38[7]={1.1E-7, 4., 1.2E-4, 1.2E-9, .17, 7.E-7, .004};
630 static const std::pair<G4int, const G4double*> Z30N38(38,pZ30N38);
631 static const G4double pZ30N40[7]={1.1E-7, 4., 1.2E-4, 1.2E-9, .17, 7.E-7, .004};
632 static const std::pair<G4int, const G4double*> Z30N40(40,pZ30N40);
633 static const std::pair<G4int, const G4double*> Z30[N30]={Z30N34, Z30N36, Z30N37,
636 static const G4int N31=2;
637 static const G4double pZ31N38[7]={5.E-8, 3.7, 1.1E-4, .55E-9, .17, 8.4E-7, .0076};
638 static const std::pair<G4int, const G4double*> Z31N38(38,pZ31N38);
639 static const G4double pZ31N40[7]={1.E-8, 3.1, 1.7E-5, .7E-9, .17, 9.E-7, .0048};
640 static const std::pair<G4int, const G4double*> Z31N40(40,pZ31N40);
641 static const std::pair<G4int, const G4double*> Z31[N31]={Z31N38, Z31N40};
643 static const G4int N32=5;
644 static const G4double pZ32N38[7]={5.E-5, 4., .17, .35E-9, .17, 9.E-7, .013};
645 static const std::pair<G4int, const G4double*> Z32N38(38,pZ32N38);
646 static const G4double pZ32N40[7]={5.E-7, 4.4, .001, .6E-9, .17, 9.E-7, .008};
647 static const std::pair<G4int, const G4double*> Z32N40(40,pZ32N40);
648 static const G4double pZ32N41[7]={5.E-9, 3., 8.E-6, .7E-9, .17, 1.E-6, .0043};
649 static const std::pair<G4int, const G4double*> Z32N41(41,pZ32N41);
650 static const G4double pZ32N42[7]={1.E-7, 4.2, 1.7E-4, .7E-9, .17, 1.E-6, .0065};
651 static const std::pair<G4int, const G4double*> Z32N42(42,pZ32N42);
652 static const G4double pZ32N44[7]={1.E-6, 4.6, .0018, .6E-9, .17, 1.E-6, .0073};
653 static const std::pair<G4int, const G4double*> Z32N44(44,pZ32N44);
654 static const std::pair<G4int, const G4double*> Z32[N32]={Z32N38, Z32N40, Z32N41,
657 static const G4int N33=2;
658 static const G4double pZ33N41[7]={1.E-8, 3.4, 1.5E-5, .72E-9, .17, 1.E-6, .0045};
659 static const std::pair<G4int, const G4double*> Z33N41(41,pZ33N41);
660 static const G4double pZ33N42[7]={1.E-8, 4.1, 1.3E-5, .75E-9, .2, 1.2E-6, .0048};
661 static const std::pair<G4int, const G4double*> Z33N42(42,pZ33N42);
662 static const std::pair<G4int, const G4double*> Z33[N33]={Z33N41, Z33N42};
664 static const G4int N34=7;
665 static const G4double pZ34N40[7]={6.E-8, 7.2, 6.E-5, 1.E-9, .32, 2.E-6, .0063};
666 static const std::pair<G4int, const G4double*> Z34N40(40,pZ34N40);
667 static const G4double pZ34N42[7]={4.E-5, 7.4, .1, .43E-9, .34, 2.1E-6, .016};
668 static const std::pair<G4int, const G4double*> Z34N42(42,pZ34N42);
669 static const G4double pZ34N43[7]={1.E-7, 6.2, 1.4E-4, .9E-9, .34, 2.1E-6, .0075};
670 static const std::pair<G4int, const G4double*> Z34N43(43,pZ34N43);
671 static const G4double pZ34N44[7]={1.E-7, 6.6, 1.3E-4, .9E-9, .34, 2.1E-6, .0075};
672 static const std::pair<G4int, const G4double*> Z34N44(44,pZ34N44);
673 static const G4double pZ34N45[7]={5.E-8, 6.6, 4.8E-5, 1.2E-9, .4, 2.6E-6, .0055};
674 static const std::pair<G4int, const G4double*> Z34N45(45,pZ34N45);
675 static const G4double pZ34N46[7]={2.E-7, 7.7, 1.3E-4, 1.7E-9, .34, 2.1E-6, .0043};
676 static const std::pair<G4int, const G4double*> Z34N46(46,pZ34N46);
677 static const G4double pZ34N48[7]={2.E-7, 8.3, 1.2E-4, 1.7E-9, .34, 2.1E-6, .0043};
678 static const std::pair<G4int, const G4double*> Z34N48(48,pZ34N48);
679 static const std::pair<G4int, const G4double*> Z34[N34]={Z34N40, Z34N42, Z34N43, Z34N44,
680 Z34N45, Z34N46, Z34N48};
682 static const G4int N35=2;
683 static const G4double pZ35N44[7]={5.E-8, 6., 2.8E-5, 2.E-9, .34, 2.1E-6, .0028};
684 static const std::pair<G4int, const G4double*> Z35N44(44,pZ35N44);
685 static const G4double pZ35N46[7]={4.E-8, 6.2, 3.7E-5, 1.1E-9, .34, 2.1E-6, .0049};
686 static const std::pair<G4int, const G4double*> Z35N46(46,pZ35N46);
687 static const std::pair<G4int, const G4double*> Z35[N35]={Z35N44, Z35N46};
689 static const G4int N36=7;
690 static const G4double pZ36N42[7]={1.6E-7, 6.8, 2.E-4, .8E-9, .35, 2.1E-6, .0076};
691 static const std::pair<G4int, const G4double*> Z36N42(42,pZ36N42);
692 static const G4double pZ36N44[7]={1.6E-7, 7.3, 1.6E-4, 1.E-9, .35, 2.1E-6, .0062};
693 static const std::pair<G4int, const G4double*> Z36N44(44,pZ36N44);
694 static const G4double pZ36N46[7]={1.6E-7, 7.3, 3.3E-4, .7E-9, .35, 2.1E-6, .013};
695 static const std::pair<G4int, const G4double*> Z36N46(46,pZ36N46);
696 static const G4double pZ36N47[7]={1.6E-6, 7.3, .003, .6E-9, .35, 2.1E-6, .011};
697 static const std::pair<G4int, const G4double*> Z36N47(47,pZ36N47);
698 static const G4double pZ36N48[7]={1.6E-7, 7.8, 7.6E-5, 2.E-9, .35, 2.1E-6, .0031};
699 static const std::pair<G4int, const G4double*> Z36N48(48,pZ36N48);
700 static const G4double pZ36N49[7]={6.E-7, 8., 4.8E-4, 1.4E-9, .27, 2.1E-6, .0053};
701 static const std::pair<G4int, const G4double*> Z36N49(49,pZ36N49);
702 static const G4double pZ36N50[7]={4.E-7, 8.1, 2.7E-4, 1.6E-9, .35, 2.1E-6, .0045};
703 static const std::pair<G4int, const G4double*> Z36N50(50,pZ36N50);
704 static const std::pair<G4int, const G4double*> Z36[N36]={Z36N42, Z36N44, Z36N46,
705 Z36N47, Z36N48, Z36N49, Z36N50};
707 static const G4int N37=3;
708 static const G4double pZ37N48[7]={1.6E-7, 7.2, 1.4E-4, 1.2E-9, .35, 2.1E-6, .0052};
709 static const std::pair<G4int, const G4double*> Z37N48(48,pZ37N48);
710 static const G4double pZ37N49[7]={8.E-8, 7.1, 4.7E-5, 1.6E-9, .27, 2.1E-6, .0034};
711 static const std::pair<G4int, const G4double*> Z37N49(49,pZ37N49);
712 static const G4double pZ37N50[7]={1.E-7, 8., 5.5E-5, 1.9E-9, .27, 1.5E-6, .0036};
713 static const std::pair<G4int, const G4double*> Z37N50(50,pZ37N50);
714 static const std::pair<G4int, const G4double*> Z37[N37]={Z37N48, Z37N49, Z37N50};
716 static const G4int N38=6;
717 static const G4double pZ38N46[7]={8.E-8, 7.3, 6.E-5, 1.3E-9, .27, 2.E-6, .0045};
718 static const std::pair<G4int, const G4double*> Z38N46(46,pZ38N46);
719 static const G4double pZ38N48[7]={8.E-8, 9.7, 2.3E-5, 3.5E-9, .4, 2.7E-6, .0023};
720 static const std::pair<G4int, const G4double*> Z38N48(48,pZ38N48);
721 static const G4double pZ38N49[7]={2.6E-7, 9.5, 1.9E-4, 1.5E-9, .4, 2.7E-6, .0057};
722 static const std::pair<G4int, const G4double*> Z38N49(49,pZ38N49);
723 static const G4double pZ38N50[7]={2.6E-7, 9.5, 2.E-4, 1.4E-9, .37, 3.2E-6, .0059};
724 static const std::pair<G4int, const G4double*> Z38N50(50,pZ38N50);
725 static const G4double pZ38N51[7]={1.3E-7, 9.9, 7.5E-5, 1.7E-9, .37, 3.2E-6, .0046};
726 static const std::pair<G4int, const G4double*> Z38N51(51,pZ38N51);
727 static const G4double pZ38N52[7]={2.6E-7, 9.6, 1.6E-4, 1.8E-9, .37, 2.7E06, .0047};
728 static const std::pair<G4int, const G4double*> Z38N52(52,pZ38N52);
729 static const std::pair<G4int, const G4double*> Z38[N38]={Z38N46, Z38N48, Z38N49, Z38N50,
732 static const G4int N39=3;
733 static const G4double pZ39N50[7]={2.6E-7, 9.9, 2.E-4, 1.1E-9, .37, 3.2E-6, .0062};
734 static const std::pair<G4int, const G4double*> Z39N50(50,pZ39N50);
735 static const G4double pZ39N51[7]={2.7E-5, 20., .013, 2.2E-9, .37, 1.9E-6, .0078};
736 static const std::pair<G4int, const G4double*> Z39N51(51,pZ39N51);
737 static const G4double pZ39N52[7]={2.E-7, 9.6, 1.2E-4, 1.7E-9, .37, 3.E-6, .0046};
738 static const std::pair<G4int, const G4double*> Z39N52(52,pZ39N52);
739 static const std::pair<G4int, const G4double*> Z39[N39]={Z39N50, Z39N51, Z39N52};
741 static const G4int N40=7;
742 static const G4double pZ40N50[7]={1.E-7, 9., 6.2E-5, 1.7E-9, .3, 3.E-6, .0044};
743 static const std::pair<G4int, const G4double*> Z40N50(50,pZ40N50);
744 static const G4double pZ40N51[7]={5.E-7, 9.8, 5.E-4, 1.E-9, .25, 2.1E-6, .0079};
745 static const std::pair<G4int, const G4double*> Z40N51(51,pZ40N51);
746 static const G4double pZ40N52[7]={3.E-7, 9.6, 2.2E-4, 1.2E-9, .25, 2.2E-6, .0056};
747 static const std::pair<G4int, const G4double*> Z40N52(52,pZ40N52);
748 static const G4double pZ40N53[7]={2.E-7, 9.6, 1.2E-4, 1.6E-9, .38, 2.9E-6, .0046};
749 static const std::pair<G4int, const G4double*> Z40N53(53,pZ40N53);
750 static const G4double pZ40N54[7]={2.E-7, 9.6, 1.8E-4, 1.1E-9, .25, 2.1E-6, .0067};
751 static const std::pair<G4int, const G4double*> Z40N54(54,pZ40N54);
752 static const G4double pZ40N55[7]={1.8E-7, 9.4, 1.1E-4, 1.6E-9, .33, 2.3E-6, .0045};
753 static const std::pair<G4int, const G4double*> Z40N55(55,pZ40N55);
754 static const G4double pZ40N56[7]={1.8E-7, 9.4, 1.1E-4, 1.6E-9, .2, 1.5E-6, .0045};
755 static const std::pair<G4int, const G4double*> Z40N56(56,pZ40N56);
756 static const std::pair<G4int, const G4double*> Z40[N40]={Z40N50, Z40N51, Z40N52, Z40N53,
757 Z40N54, Z40N55, Z40N56};
759 static const G4int N41=3;
760 static const G4double pZ41N52[7]={2.6E-7, 8.3, 2.E-4, 1.2E-9, .4, 4.E-6, .0051};
761 static const std::pair<G4int, const G4double*> Z41N52(52,pZ41N52);
762 static const G4double pZ41N53[7]={2.E-7, 8.3, 1.6E-4, 1.4E-9, .35, 2.5E-6, .0051};
763 static const std::pair<G4int, const G4double*> Z41N53(53,pZ41N53);
764 static const G4double pZ41N54[7]={1.5E-7, 8.6, 1.E-4, 1.5E-9, .35, 2.5E-6, .0045};
765 static const std::pair<G4int, const G4double*> Z41N54(54,pZ41N54);
766 static const std::pair<G4int, const G4double*> Z41[N41]={Z41N52, Z41N52, Z41N52};
768 static const G4int N42=8;
769 static const G4double pZ42N50[7]={2.E-7, 10., 1.1E-4, 1.8E-9, .3, 2.7E-6, .0044};
770 static const std::pair<G4int, const G4double*> Z42N50(50,pZ42N50);
771 static const G4double pZ42N52[7]={2.1E-7, 10., 1.2E-4, 1.7E-9, .3, 2.8E-6, .0046};
772 static const std::pair<G4int, const G4double*> Z42N52(52,pZ42N52);
773 static const G4double pZ42N53[7]={3.E-7, 10., 1.9E-4, 1.5E-9, .29, 3.E-6, .005};
774 static const std::pair<G4int, const G4double*> Z42N53(53,pZ42N53);
775 static const G4double pZ42N54[7]={1.5E-7, 10., 7.1E-5, 2.1E-9, .29, 2.9E-6, .0037};
776 static const std::pair<G4int, const G4double*> Z42N54(54,pZ42N54);
777 static const G4double pZ42N55[7]={1.9E-7, 9., 1.4E-4, 1.3E-9, .29, 2.8E-6, .0052};
778 static const std::pair<G4int, const G4double*> Z42N55(55,pZ42N55);
779 static const G4double pZ42N56[7]={1.9E-7, 9.9, 1.1E-4, 1.8E-9, .29, 2.4E-6, .0044};
780 static const std::pair<G4int, const G4double*> Z42N56(56,pZ42N56);
781 static const G4double pZ42N57[7]={1.4E-7, 8., 1.E-4, 1.4E-9, .34, 2.5E-6, .0044};
782 static const std::pair<G4int, const G4double*> Z42N57(57,pZ42N57);
783 static const G4double pZ42N58[7]={1.8E-7, 9.5, 1.E-4, 1.7E-9, .27, 2.2E-6, .0041};
784 static const std::pair<G4int, const G4double*> Z42N58(58,pZ42N58);
785 static const std::pair<G4int, const G4double*> Z42[N42]={Z42N50, Z42N52, Z42N53, Z42N54,
786 Z42N55, Z42N56, Z42N57, Z42N58};
788 static const G4int N43=1;
789 static const G4double pZ43N56[7]={1.E-7, 8., 7.2E-5, 1.4E-9, .24, 2.5E-6, .0044};
790 static const std::pair<G4int, const G4double*> Z43N56(56,pZ43N56);
791 static const std::pair<G4int, const G4double*> Z43[N43]={Z43N56};
793 static const G4int N44=9;
794 static const G4double pZ44N52[7]={1.9E-7, 10., 1.E-4, 2.1E-9, .4, 4.E-6, .004};
795 static const std::pair<G4int, const G4double*> Z44N52(52,pZ44N52);
796 static const G4double pZ44N54[7]={1.5E-7, 10., 7.7E-5, 2.1E-9, .29, 2.9E-6, .004};
797 static const std::pair<G4int, const G4double*> Z44N54(54,pZ44N54);
798 static const G4double pZ44N55[7]={1.8E-7, 10., 6.6E-5, 2.6E-9, .47, 4.6E-6, .0028};
799 static const std::pair<G4int, const G4double*> Z44N55(55,pZ44N55);
800 static const G4double pZ44N56[7]={1.8E-6, 10., .0017, 1.1E-9, .47, 6.E-6, .0073};
801 static const std::pair<G4int, const G4double*> Z44N56(56,pZ44N56);
802 static const G4double pZ44N57[7]={1.8E-7, 7.8, 1.3E-4, 1.4E-9, .42, 5.E-6, .0043};
803 static const std::pair<G4int, const G4double*> Z44N57(57,pZ44N57);
804 static const G4double pZ44N58[7]={1.7E-6, 9.8, .0015, 1.2E-9, .32, 4.E-6, .0065};
805 static const std::pair<G4int, const G4double*> Z44N58(58,pZ44N58);
806 static const G4double pZ44N59[7]={3.3E-7, 8.7, 1.9E-4, 1.6E-9, .32, 3.8E-6, .0038};
807 static const std::pair<G4int, const G4double*> Z44N59(59,pZ44N59);
808 static const G4double pZ44N60[7]={3.E-7, 8.7, 1.8E-4, 1.6E-9, .3, 3.2E-6, .004};
809 static const std::pair<G4int, const G4double*> Z44N60(60,pZ44N60);
810 static const G4double pZ44N61[7]={3.E-7, 8.8, 1.4E-4, 1.9E-9, .3, 3.2E-6, .003};
811 static const std::pair<G4int, const G4double*> Z44N61(61,pZ44N61);
812 static const std::pair<G4int, const G4double*> Z44[N44]={Z44N52, Z44N54, Z44N55, Z44N56,
813 Z44N57, Z44N58, Z44N59, Z44N60,
816 static const G4int N45=2;
817 static const G4double pZ45N58[7]={8.E-8, 8.7, 4.E-5, 1.8E-9, .29, 2.9E-6, .0033};
818 static const std::pair<G4int, const G4double*> Z45N58(58,pZ45N58);
819 static const G4double pZ45N60[7]={8.E-8, 8.7, .09, 1.3E-12, .29, 2.9E-6, 7.};
820 static const std::pair<G4int, const G4double*> Z45N60(60,pZ45N60);
821 static const std::pair<G4int, const G4double*> Z45[N45]={Z45N58, Z45N60};
823 static const G4int N46=7;
824 static const G4double pZ46N56[7]={2.E-7, 9.9, 1.2E-4, 1.5E-9, .35, 3.3E-6, .0045};
825 static const std::pair<G4int, const G4double*> Z46N56(56,pZ46N56);
826 static const G4double pZ46N58[7]={2.E-7, 9.9, 9.5E-5, 1.8E-9, .4, 4.E-6, .0036};
827 static const std::pair<G4int, const G4double*> Z46N58(58,pZ46N58);
828 static const G4double pZ46N59[7]={5.6E-7, 9., 4.6E-4, 1.2E-9, .5, 4.8E-6, .0056};
829 static const std::pair<G4int, const G4double*> Z46N59(59,pZ46N59);
830 static const G4double pZ46N60[7]={2.4E-7, 9.2, 1.2E-4, 1.8E-9, .47, 4.6E-6, .0035};
831 static const std::pair<G4int, const G4double*> Z46N60(60,pZ46N60);
832 static const G4double pZ46N61[7]={1.2E-7, 9.2, 4.4E-5, 2.8E-9, .5, 4.3E-6, .0025};
833 static const std::pair<G4int, const G4double*> Z46N61(61,pZ46N61);
834 static const G4double pZ46N62[7]={1.2E-7, 9.2, 3.2E-5, 3.4E-9, .48, 4.5E-6, .0018};
835 static const std::pair<G4int, const G4double*> Z46N62(62,pZ46N62);
836 static const G4double pZ46N64[7]={4.E-7, 9.1, 2.5E-4, 1.5E-9, .48, 4.7E-6, .0042};
837 static const std::pair<G4int, const G4double*> Z46N64(64,pZ46N64);
838 static const std::pair<G4int, const G4double*> Z46[N46]={Z46N56, Z46N58, Z46N59,
839 Z46N60, Z46N61, Z46N62, Z46N64};
841 static const G4int N47=4;
842 static const G4double pZ47N60[7]={1.4E-6, 9.7, .0011, 1.4E-9, .55, 5.E-6, .0056};
843 static const std::pair<G4int, const G4double*> Z47N60(60,pZ47N60);
844 static const G4double pZ47N62[7]={3.E-8, 8.7, 8.5E-6, 3.5E-9, .6, 5.6E-6, .0018};
845 static const std::pair<G4int, const G4double*> Z47N62(62,pZ47N62);
846 static const G4double pZ47N63[7]={3.E-6, 9.5, .002, 1.5E-9, .58, 5.E-6, .0047};
847 static const std::pair<G4int, const G4double*> Z47N63(63,pZ47N63);
848 static const G4double pZ47N64[7]={1.5E-7, 9., 9.E-5, 1.7E-9, .58, 5.6E-6, .0039};
849 static const std::pair<G4int, const G4double*> Z47N64(64,pZ47N64);
850 static const std::pair<G4int, const G4double*> Z47[N47]={Z47N60, Z47N62, Z47N63, Z47N64};
852 static const G4int N48=9;
853 static const G4double pZ48N58[7]={2.9E-7, 10., 1.3E-4, 1.9E-9, .4, 3.8E-6, .0034};
854 static const std::pair<G4int, const G4double*> Z48N58(58,pZ48N58);
855 static const G4double pZ48N60[7]={2.3E-7, 10., 8.2E-5, 2.5E-9, .5, 4.7E-6, .0026};
856 static const std::pair<G4int, const G4double*> Z48N60(60,pZ48N60);
857 static const G4double pZ48N62[7]={2.3E-7, 10., 9.9E-5, 2.5E-9, .5, 4.7E-6, .0031};
858 static const std::pair<G4int, const G4double*> Z48N62(62,pZ48N62);
859 static const G4double pZ48N63[7]={8.4E-7, 11., 4.3E-4, 1.8E-9, .5, 4.5E-6, .0042};
860 static const std::pair<G4int, const G4double*> Z48N63(63,pZ48N63);
861 static const G4double pZ48N64[7]={4.E-7, 11., 1.8E-4, 1.8E-9, .5, 4.6E-6, .0036};
862 static const std::pair<G4int, const G4double*> Z48N64(64,pZ48N64);
863 static const G4double pZ48N65[7]={1.6E-6, 12., .001, 0., .5, 4.6E-6, .013};
864 static const std::pair<G4int, const G4double*> Z48N65(65,pZ48N65);
865 static const G4double pZ48N66[7]={3.E-7, 11., 1.2E-4, 1.9E-9, .5, 4.6E-6, .0031};
866 static const std::pair<G4int, const G4double*> Z48N66(66,pZ48N66);
867 static const G4double pZ48N67[7]={3.8E-7, 11., 1.7E-4, 2.E-9, .6, 6.6E-6, .0035};
868 static const std::pair<G4int, const G4double*> Z48N67(67,pZ48N67);
869 static const G4double pZ48N68[7]={6.E-7, 11., 3.3E-4, 1.9E-9, .5, 4.6E-6, .0043};
870 static const std::pair<G4int, const G4double*> Z48N68(68,pZ48N68);
871 static const std::pair<G4int, const G4double*> Z48[N48]={Z48N58, Z48N60, Z48N62, Z48N63,
872 Z48N64, Z48N65, Z48N66, Z48N67,
875 static const G4int N49=2;
876 static const G4double pZ49N64[7]={2.7E-7, 12., 8.1E-5, 2.7E-9, .5, 5.E-6, .0026};
877 static const std::pair<G4int, const G4double*> Z49N64(64,pZ49N64);
878 static const G4double pZ49N66[7]={2.7E-7, 12., 5.5E-5, 4.E-9, .5, 5.E-6, .0018};
879 static const std::pair<G4int, const G4double*> Z49N66(66,pZ49N66);
880 static const std::pair<G4int, const G4double*> Z49[N49]={Z49N64, Z49N66};
882 static const G4int N50=14;
883 static const G4double pZ50N62[7]={4.E-7, 11., 1.6E-4, 2.2E-9, .5, 4.5E-6, .0032};
884 static const std::pair<G4int, const G4double*> Z50N62(62,pZ50N62);
885 static const G4double pZ50N63[7]={4.1E-7, 11., 1.6E-4, 2.4E-9, .54, 6.E-6, .0031};
886 static const std::pair<G4int, const G4double*> Z50N63(63,pZ50N63);
887 static const G4double pZ50N64[7]={5.E-7, 12., 1.9E-4, 2.2E-9, .5, 4.4E-6, .0032};
888 static const std::pair<G4int, const G4double*> Z50N64(64,pZ50N64);
889 static const G4double pZ50N65[7]={1.E-5, 12., .0077, 1.4E-9, .5, 5.E-6, .0066};
890 static const std::pair<G4int, const G4double*> Z50N65(65,pZ50N65);
891 static const G4double pZ50N66[7]={5.E-7, 12., 1.8E-4, 2.4E-9, .5, 5.E-6, .0031};
892 static const std::pair<G4int, const G4double*> Z50N66(66,pZ50N66);
893 static const G4double pZ50N67[7]={1.E-6, 12., 4.4E-4, 1.8E-9, .5, 5.E-6, .0037};
894 static const std::pair<G4int, const G4double*> Z50N67(67,pZ50N67);
895 static const G4double pZ50N68[7]={5.E-7, 12., 2.E-4, 2.4E-9, .5, 5.E-6, .0033};
896 static const std::pair<G4int, const G4double*> Z50N68(68,pZ50N68);
897 static const G4double pZ50N69[7]={6.E-7, 12., 2.5E-4, 2.E-9, .5, 5.E-6, .0035};
898 static const std::pair<G4int, const G4double*> Z50N69(69,pZ50N69);
899 static const G4double pZ50N70[7]={1.E-6, 12., 4.7E-4, 2.E-9, .5, 5.E-6, .0039};
900 static const std::pair<G4int, const G4double*> Z50N70(70,pZ50N70);
901 static const G4double pZ50N72[7]={1.E-6, 12., 3.7E-4, 2.2E-9, .5, 5.E-6, .0031};
902 static const std::pair<G4int, const G4double*> Z50N72(72,pZ50N72);
903 static const G4double pZ50N73[7]={5.E-7, 12., 1.7E-4, 2.8E-9, .5, 5.E-6, .0028};
904 static const std::pair<G4int, const G4double*> Z50N73(73,pZ50N73);
905 static const G4double pZ50N74[7]={5.E-7, 12., 2.E-4, 2.E-9, .5, 5.E-6, .0033};
906 static const std::pair<G4int, const G4double*> Z50N74(74,pZ50N74);
907 static const G4double pZ50N75[7]={5.E-7, 12., 1.9E-4, 2.8E-9, .5, 5.E-6, .003};
908 static const std::pair<G4int, const G4double*> Z50N75(75,pZ50N75);
909 static const G4double pZ50N76[7]={5.E-7, 12., 1.7E-4, 2.8E-9, .5, 5.E-6, .0028};
910 static const std::pair<G4int, const G4double*> Z50N76(76,pZ50N76);
911 static const std::pair<G4int, const G4double*> Z50[N50]={Z50N62, Z50N63, Z50N64, Z50N65,
912 Z50N66, Z50N67, Z50N68, Z50N69,
913 Z50N70, Z50N72, Z50N73, Z50N74,
916 static const G4int N51=5;
917 static const G4double pZ51N70[7]={6.E-7, 12., 2.E-4, 2.8E-9, .5, 5.E-6, .0028};
918 static const std::pair<G4int, const G4double*> Z51N70(70,pZ51N70);
919 static const G4double pZ51N72[7]={6.E-7, 12., 1.9E-4, 3.E-9, .5, 5.E-6, .0025};
920 static const std::pair<G4int, const G4double*> Z51N72(72,pZ51N72);
921 static const G4double pZ51N73[7]={1.1E-6, 12., 3.5E-4, 2.9E-9, .5, 5.E-6, .0026};
922 static const std::pair<G4int, const G4double*> Z51N73(73,pZ51N73);
923 static const G4double pZ51N74[7]={5.5E-7, 12., 1.9E-4, 2.9E-9, .5, 5.E-6, .0027};
924 static const std::pair<G4int, const G4double*> Z51N74(74,pZ51N74);
925 static const G4double pZ51N75[7]={6.E-7, 12., 2.E-4, 2.9E-9, .5, 5.E-6, .0027};
926 static const std::pair<G4int, const G4double*> Z51N75(75,pZ51N75);
927 static const std::pair<G4int, const G4double*> Z51[N51]={Z51N70, Z51N72, Z51N73, Z51N74,
930 static const G4int N52=11;
931 static const G4double pZ52N68[7]={2.7E-7, 12., 8.4E-5, 3.2E-9, 1., 8.E-6, .0026};
932 static const std::pair<G4int, const G4double*> Z52N68(68,pZ52N68);
933 static const G4double pZ52N70[7]={2.7E-7, 12., 3.8E-5, 6.E-9, 1., 8.E-6, .0012};
934 static const std::pair<G4int, const G4double*> Z52N70(70,pZ52N70);
935 static const G4double pZ52N71[7]={2.7E-8, 12., 1.8E-6, 2.E-8, 1., 8.E-6, 4.8E-4};
936 static const std::pair<G4int, const G4double*> Z52N71(71,pZ52N71);
937 static const G4double pZ52N72[7]={2.6E-6, 14., .0014, 2.E-9, 1., 9.E-6, .005};
938 static const std::pair<G4int, const G4double*> Z52N72(72,pZ52N72);
939 static const G4double pZ52N73[7]={1.E-6, 14., 2.4E-4, 3.9E-9, 1., 9.E-6, .0022};
940 static const std::pair<G4int, const G4double*> Z52N73(73,pZ52N73);
941 static const G4double pZ52N74[7]={8.E-7, 14., 2.3E-4, 3.6E-9, 1.4, 1.3E-5, .0028};
942 static const std::pair<G4int, const G4double*> Z52N74(74,pZ52N74);
943 static const G4double pZ52N75[7]={8.E-7, 14., 2.1E-4, 3.6E-9, 1.4, 1.3E-5, .0025};
944 static const std::pair<G4int, const G4double*> Z52N75(75,pZ52N75);
945 static const G4double pZ52N76[7]={8.E-7, 14., 2.5E-4, 3.E-9, 1.4, 1.3E-5, .003};
946 static const std::pair<G4int, const G4double*> Z52N76(76,pZ52N76);
947 static const G4double pZ52N77[7]={5.E-7, 15., 1.2E-4, 4.3E-9, 1.4, 1.4E-5, .0023};
948 static const std::pair<G4int, const G4double*> Z52N77(77,pZ52N77);
949 static const G4double pZ52N78[7]={8.E-7, 14., 2.7E-4, 2.7E-9, 1.4, 1.3E-5, .0031};
950 static const std::pair<G4int, const G4double*> Z52N78(78,pZ52N78);
951 static const G4double pZ52N80[7]={4.7E-7, 14., 1.8E-4, 2.2E-9, .83, 1.E-5, .0036};
952 static const std::pair<G4int, const G4double*> Z52N80(80,pZ52N80);
953 static const std::pair<G4int, const G4double*> Z52[N52]={Z52N68, Z52N70, Z52N71, Z52N72,
954 Z52N73, Z52N74, Z52N75, Z52N76,
955 Z52N77, Z52N78, Z52N80};
957 static const G4int N53=5;
958 static const G4double pZ53N74[7]={9.4E-7, 14., 2.5E-4, 3.E-9, .7, 7.3E-6, .0025};
959 static const std::pair<G4int, const G4double*> Z53N74(74,pZ53N74);
960 static const G4double pZ53N76[7]={2.1E-5, 14., .015, 1.1E-9, 1.1, 1.E-5, .007};
961 static const std::pair<G4int, const G4double*> Z53N76(76,pZ53N76);
962 static const G4double pZ53N77[7]={1.1E-6, 14., 2.4E-4, 3.3E-9, .9, 1.E-5, .0021};
963 static const std::pair<G4int, const G4double*> Z53N77(77,pZ53N77);
964 static const G4double pZ53N78[7]={5.5E-7, 14., 1.5E-4, 3.7E-9, 1.2, 1.2E-5, .0024};
965 static const std::pair<G4int, const G4double*> Z53N78(78,pZ53N78);
966 static const G4double pZ53N82[7]={3.2E-6, 14., .0017, 1.8E-9, .8, 8.E-6, .0024};
967 static const std::pair<G4int, const G4double*> Z53N82(82,pZ53N82);
968 static const std::pair<G4int, const G4double*> Z53[N53]={Z53N74, Z53N76, Z53N77, Z53N78,
971 static const G4int N54=12;
972 static const G4double pZ54N69[7]={3.E-6, 14., 8.E-4, 3.7E-9, .9, 1.1E-5, .15};
973 static const std::pair<G4int, const G4double*> Z54N69(69,pZ54N69);
974 static const G4double pZ54N70[7]={1.5E-7, 14., 1.4E-6, 9.E-8, .7, 8.E-6, 9.5E-5};
975 static const std::pair<G4int, const G4double*> Z54N70(70,pZ54N70);
976 static const G4double pZ54N72[7]={1.5E-6, 14., 5.6E-4, 3.E-9, 1.2, 1.1E-5, .0036};
977 static const std::pair<G4int, const G4double*> Z54N72(72,pZ54N72);
978 static const G4double pZ54N74[7]={1.8E-6, 14., 8.8E-4, 2.E-9, 1.3, 1.2E-5, .0047};
979 static const std::pair<G4int, const G4double*> Z54N74(74,pZ54N74);
980 static const G4double pZ54N75[7]={1.E-6, 14., 2.6E-4, 3.7E-9, 1.5, 1.4E-5, .0024};
981 static const std::pair<G4int, const G4double*> Z54N75(75,pZ54N75);
982 static const G4double pZ54N76[7]={1.8E-6, 14., 8.E-4, 2.E-9, 1.2, 1.4E-5, .0042};
983 static const std::pair<G4int, const G4double*> Z54N76(76,pZ54N76);
984 static const G4double pZ54N77[7]={2.3E-7, 14., 1.9E-5, 9.E-9, 1.2, 1.4E-5, 7.7E-4};
985 static const std::pair<G4int, const G4double*> Z54N77(77,pZ54N77);
986 static const G4double pZ54N78[7]={6.E-7, 14., 1.6E-4, 3.E-9, 1.2, 1.4E-5, .0025};
987 static const std::pair<G4int, const G4double*> Z54N78(78,pZ54N78);
988 static const G4double pZ54N79[7]={6.E-7, 14., 1.6E-4, 3.3E-9, 1.6, 1.5E-5, .0024};
989 static const std::pair<G4int, const G4double*> Z54N79(79,pZ54N79);
990 static const G4double pZ54N80[7]={6.6E-7, 14., 2.1E-4, 2.5E-9, 1.2, 1.4E-5, .003};
991 static const std::pair<G4int, const G4double*> Z54N80(80,pZ54N80);
992 static const G4double pZ54N81[7]={.03, 40., 2.1, 2.5E-9, 1.E-16, 6.E-36, 140.};
993 static const std::pair<G4int, const G4double*> Z54N81(81,pZ54N81);
994 static const G4double pZ54N82[7]={3.1E-6, 14., .0019, 1.6E-9, 1., 1.3E-5, .0054};
995 static const std::pair<G4int, const G4double*> Z54N82(82,pZ54N82);
996 static const std::pair<G4int, const G4double*> Z54[N54]={Z54N69, Z54N70, Z54N72, Z54N74,
997 Z54N75, Z54N76, Z54N77, Z54N78,
998 Z54N79, Z54N80, Z54N81, Z54N82};
1000 static const G4int N55=5;
1001 static const G4double pZ55N78[7]={1.4E-6, 14., 4.E-4, 3.E-9, 1.2, 1.4E-5, .0026};
1002 static const std::pair<G4int, const G4double*> Z55N78(78,pZ55N78);
1003 static const G4double pZ55N79[7]={.028, 14., 44., .5E-9, 1.2, 1.3E-5, .015};
1004 static const std::pair<G4int, const G4double*> Z55N79(79,pZ55N79);
1005 static const G4double pZ55N80[7]={2.E-6, 14., 9.5E-4, 2.E-9, 1.2, 1.4E-5, .0042};
1006 static const std::pair<G4int, const G4double*> Z55N80(80,pZ55N80);
1007 static const G4double pZ55N81[7]={2.5E-7, 14., 6.5E-5, 3.8E-9, 1.2, 1.4E-5, .0023};
1008 static const std::pair<G4int, const G4double*> Z55N81(81,pZ55N81);
1009 static const G4double pZ55N82[7]={2.5E-7, 14., 6.5E-5, 3.8E-9, 1.2, 1.4E-5, .0023};
1010 static const std::pair<G4int, const G4double*> Z55N82(82,pZ55N82);
1011 static const std::pair<G4int, const G4double*> Z55[N55]={Z55N78, Z55N79, Z55N80, Z55N81,
1014 static const G4int N56=9;
1015 static const G4double pZ56N74[7]={4.E-7, 14., 2.8E-5, 1.2E-8, 1.2, 1.5E-5, 6.6E-4};
1016 static const std::pair<G4int, const G4double*> Z56N74(74,pZ56N74);
1017 static const G4double pZ56N76[7]={4.E-6, 14., .0022, 1.4E-9, 1.3, 1.6E-5, .0053};
1018 static const std::pair<G4int, const G4double*> Z56N76(76,pZ56N76);
1019 static const G4double pZ56N77[7]={2.E-7, 14., 3.7E-5, 5.E-9, 1.1, 1.5E-5, .0016};
1020 static const std::pair<G4int, const G4double*> Z56N77(77,pZ56N77);
1021 static const G4double pZ56N78[7]={1.6E-6, 14., 6.E-4, 2.E-9, 1.3, 1.6E-5, .0033};
1022 static const std::pair<G4int, const G4double*> Z56N78(78,pZ56N78);
1023 static const G4double pZ56N79[7]={5.E-7, 17., 8.E-5, 4.5E-9, 1.3, 1.6E-5, .0018};
1024 static const std::pair<G4int, const G4double*> Z56N79(79,pZ56N79);
1025 static const G4double pZ56N80[7]={2.E-6, 20., 3.E-4, 6.E-9, 1.3, 1.8E-5, .0019};
1026 static const std::pair<G4int, const G4double*> Z56N80(80,pZ56N80);
1027 static const G4double pZ56N81[7]={5.8E-6, 20., .0018, 3.E-9, 1.3, 1.7E-5, .0041};
1028 static const std::pair<G4int, const G4double*> Z56N81(81,pZ56N81);
1029 static const G4double pZ56N82[7]={2.7E-6, 20., 5.5E-4, 4.E-9, 1.4, 2.E-5, .0027};
1030 static const std::pair<G4int, const G4double*> Z56N82(82,pZ56N82);
1031 static const G4double pZ56N84[7]={1.1E-6, 21., 1.E-4, 9.E-9, 2., 2.7E-5, .0012};
1032 static const std::pair<G4int, const G4double*> Z56N84(84,pZ56N84);
1033 static const std::pair<G4int, const G4double*> Z56[N56]={Z56N74, Z56N76, Z56N77, Z56N78,
1034 Z56N79, Z56N80, Z56N81, Z56N82,
1037 static const G4int N57=3;
1038 static const G4double pZ57N81[7]={2.7E-6, 20., .0017, 1.1E-9, 1.4, 2.E-5, .0083};
1039 static const std::pair<G4int, const G4double*> Z57N81(81,pZ57N81);
1040 static const G4double pZ57N82[7]={5.4E-6, 20., .0027, 1.5E-9, 1., 1.5E-5, .0065};
1041 static const std::pair<G4int, const G4double*> Z57N82(82,pZ57N82);
1042 static const G4double pZ57N83[7]={2.7E-6, 20., 2.6E-4, 6.E-9, 1.4, 2.E-5, .0012};
1043 static const std::pair<G4int, const G4double*> Z57N83(83,pZ57N83);
1044 static const std::pair<G4int, const G4double*> Z57[N57]={Z57N81, Z57N82, Z57N83};
1046 static const G4int N58=8;
1047 static const G4double pZ58N78[7]={1.8E-6, 20., 3.7E-4, 4.E-9, 1.4, 2.E-5, .0028};
1048 static const std::pair<G4int, const G4double*> Z58N78(78,pZ58N78);
1049 static const G4double pZ58N80[7]={1.8E-6, 18., 2.6E-4, 6.E-9, 1.3, 2.1E-5, .0017};
1050 static const std::pair<G4int, const G4double*> Z58N80(80,pZ58N80);
1051 static const G4double pZ58N81[7]={.0018, 18., 2.9, .6E-9, 1.3, 2.E-5, .02};
1052 static const std::pair<G4int, const G4double*> Z58N81(81,pZ58N81);
1053 static const G4double pZ58N82[7]={1.8E-6, 18., 3.7E-4, 5.E-9, 1.3, 2.E-5, .0024};
1054 static const std::pair<G4int, const G4double*> Z58N82(82,pZ58N82);
1055 static const G4double pZ58N83[7]={7.2E-6, 20., .0025, 2.3E-9, 1.1, 1.7E-5, .0045};
1056 static const std::pair<G4int, const G4double*> Z58N83(83,pZ58N83);
1057 static const G4double pZ58N84[7]={1.2E-6, 18., 2.E-4, 6.E-9, 1.5, 1.9E-5, .0018};
1058 static const std::pair<G4int, const G4double*> Z58N84(84,pZ58N84);
1059 static const G4double pZ58N85[7]={1.2E-6, 16., 6.E-4, 1.7E-9, 1.4, 2.E-5, .0053};
1060 static const std::pair<G4int, const G4double*> Z58N85(85,pZ58N85);
1061 static const G4double pZ58N86[7]={6.E-7, 18., 1.E-4, 6.E-9, 1.5, 1.9E-5, .0018};
1062 static const std::pair<G4int, const G4double*> Z58N86(86,pZ58N86);
1063 static const std::pair<G4int, const G4double*> Z58[N58]={Z58N78, Z58N80, Z58N81, Z58N82,
1064 Z58N83, Z58N84, Z58N85, Z58N86};
1066 static const G4int N59=3;
1067 static const G4double pZ59N82[7]={9.5E-7, 16., 1.6E-4, 4.E-9, 1.4, 2.E-5, .0017};
1068 static const std::pair<G4int, const G4double*> Z59N82(82,pZ59N82);
1069 static const G4double pZ59N83[7]={9.5E-7, 16., 1.9E-4, 4.E-9, 1.4, 1.9E-5, .0021};
1070 static const std::pair<G4int, const G4double*> Z59N83(83,pZ59N83);
1071 static const G4double pZ59N84[7]={9.5E-6, 16., .019, .4E-9, 2., 2.4E-5, .021};
1072 static const std::pair<G4int, const G4double*> Z59N84(84,pZ59N84);
1073 static const std::pair<G4int, const G4double*> Z59[N59]={Z59N82, Z59N83, Z59N84};
1075 static const G4int N60=8;
1076 static const G4double pZ60N82[7]={9.6E-6, 21., .0036, 2.3E-9, 1.4, 2.E-5, .0052};
1077 static const std::pair<G4int, const G4double*> Z60N82(82,pZ60N82);
1078 static const G4double pZ60N83[7]={9.6E-4, 20., 4., .25E-9, 1.4, 2.E-5, .052};
1079 static const std::pair<G4int, const G4double*> Z60N83(83,pZ60N83);
1080 static const G4double pZ60N84[7]={4.8E-7, 21., 3.3E-5, 1.E-08, 1.3, 2.2E-5, 9.E-4};
1081 static const std::pair<G4int, const G4double*> Z60N84(84,pZ60N84);
1082 static const G4double pZ60N85[7]={.0048, 20., 4.5, .9E-9, 1.3, 2.E-5, .012};
1083 static const std::pair<G4int, const G4double*> Z60N85(85,pZ60N85);
1084 static const G4double pZ60N86[7]={1.2E-6, 16., 7.7E-4, 1.5E-9, 1.3, 1.8E-5, .0066};
1085 static const std::pair<G4int, const G4double*> Z60N86(86,pZ60N86);
1086 static const G4double pZ60N87[7]={.0012, 15., 8.4, .1E-9, 1.3, 1.5E-5, .071};
1087 static const std::pair<G4int, const G4double*> Z60N87(87,pZ60N87);
1088 static const G4double pZ60N88[7]={1.5E-7, 16., 4.1E-5, 2.5E-9, 1.3, 1.6E-5, .0027};
1089 static const std::pair<G4int, const G4double*> Z60N88(88,pZ60N88);
1090 static const G4double pZ60N90[7]={1.5E-7, 16., 4.3E-5, 2.5E-9, 1.3, 1.6E-5, .0029};
1091 static const std::pair<G4int, const G4double*> Z60N90(90,pZ60N90);
1092 static const std::pair<G4int, const G4double*> Z60[N60]={Z60N82, Z60N83, Z60N84, Z60N85,
1093 Z60N86, Z60N87, Z60N88, Z60N90};
1095 static const G4int N61=3;
1096 static const G4double pZ61N86[7]={6.E-7, 16., 8.E-4, .6E-9, 3., 2.8E-5, .014};
1097 static const std::pair<G4int, const G4double*> Z61N86(86,pZ61N86);
1098 static const G4double pZ61N87[7]={6.2E-8, 16., 1.2E-5, 4.E-9, 2.2, 2.5E-5, .0019};
1099 static const std::pair<G4int, const G4double*> Z61N87(87,pZ61N87);
1100 static const G4double pZ61N88[7]={3.2E-8, 16., 6.4E-6, 4.E-9, 2.2, 2.5E-5, .002};
1101 static const std::pair<G4int, const G4double*> Z61N88(88,pZ61N88);
1102 static const std::pair<G4int, const G4double*> Z61[N61]={Z61N86, Z61N87, Z61N88};
1104 static const G4int N62=9;
1105 static const G4double pZ62N82[7]={1.2E-7, 16., 2.1E-5, 5.E-9, 1.4, 2.E-5, .0017};
1106 static const std::pair<G4int, const G4double*> Z62N82(82,pZ62N82);
1107 static const G4double pZ62N85[7]={1.2E-7, 16., 5.3E-5, 1.5E-9, 1.3, 1.7E-5, .0045};
1108 static const std::pair<G4int, const G4double*> Z62N85(85,pZ62N85);
1109 static const G4double pZ62N86[7]={6.E-8, 16., 1.7E-5, 3.E-9, 1.3, 1.7E-5, .0028};
1110 static const std::pair<G4int, const G4double*> Z62N86(86,pZ62N86);
1111 static const G4double pZ62N87[7]={6.E-8, 15., 5.2E-4, .11E-9, 1.3, 1.5E-5, .074};
1112 static const std::pair<G4int, const G4double*> Z62N87(87,pZ62N87);
1113 static const G4double pZ62N88[7]={6.E-7, 16., 8.6E-4, .7E-9, 1.3, 1.7E-5, .015};
1114 static const std::pair<G4int, const G4double*> Z62N88(88,pZ62N88);
1115 static const G4double pZ62N89[7]={6.E-7, 16., 5.E-4, 0., 1.3, 1.6E-5, .053};
1116 static const std::pair<G4int, const G4double*> Z62N89(89,pZ62N89);
1117 static const G4double pZ62N90[7]={6.E-8, 15., 1.3E-5, 4.5E-9, 1.3, 1.7E-5, .0019};
1118 static const std::pair<G4int, const G4double*> Z62N90(90,pZ62N90);
1119 static const G4double pZ62N91[7]={6.E-8, 15., 1.5E-5, 2.E-9, 1.3, 1.6E-5, .0024};
1120 static const std::pair<G4int, const G4double*> Z62N91(91,pZ62N91);
1121 static const G4double pZ62N92[7]={1.2E-7, 15., 8.6E-5, 1.2E-9, 1.3, 1.6E-5, .007};
1122 static const std::pair<G4int, const G4double*> Z62N92(92,pZ62N92);
1123 static const std::pair<G4int, const G4double*> Z62[N62]={Z62N82, Z62N85, Z62N86, Z62N87,
1124 Z62N88, Z62N89, Z62N90, Z62N91,
1127 static const G4int N63=7;
1128 static const G4double pZ63N88[7]={6.E-8, 15., 2.8E-5, 2.E-9, 1.3, 1.5E-5, .0046};
1129 static const std::pair<G4int, const G4double*> Z63N88(88,pZ63N88);
1130 static const G4double pZ63N89[7]={6.E-7, 15., .0011, .5E-9, 2.4, 2.4E-5, .017};
1131 static const std::pair<G4int, const G4double*> Z63N89(89,pZ63N89);
1132 static const G4double pZ63N90[7]={3.E-7, 15., 1.8E-4, 1.1E-9, 1., 1.2E-5, .0054};
1133 static const std::pair<G4int, const G4double*> Z63N90(90,pZ63N90);
1134 static const G4double pZ63N91[7]={4.1E-7, 15., 1.4E-4, 1.9E-9, 1., 1.4E-5, .0032};
1135 static const std::pair<G4int, const G4double*> Z63N91(91,pZ63N91);
1136 static const G4double pZ63N92[7]={5.E-8, 15., 2.4E-5, 2.8E-9, 1., 1.3E-5, .0037};
1137 static const std::pair<G4int, const G4double*> Z63N92(92,pZ63N92);
1138 static const G4double pZ63N93[7]={4.1E-8, 17., 1.6E-5, 2.E-9, 3.3, 3.4E-5, .004};
1139 static const std::pair<G4int, const G4double*> Z63N93(93,pZ63N93);
1140 static const G4double pZ63N94[7]={4.2E-8, 17., 1.6E-5, 2.E-9, 1.2, 1.6E-5, .004};
1141 static const std::pair<G4int, const G4double*> Z63N94(94,pZ63N94);
1142 static const std::pair<G4int, const G4double*> Z63[N63]={Z63N88, Z63N89, Z63N90, Z63N91,
1143 Z63N92, Z63N93, Z63N94};
1145 static const G4int N64=8;
1146 static const G4double pZ64N88[7]={4.2E-8, 14., 2.E-4, 0., 1.2, 1.3E-5, .19};
1147 static const std::pair<G4int, const G4double*> Z64N88(88,pZ64N88);
1148 static const G4double pZ64N89[7]={2.E-6, 14., .0016, 1.4E-9, 1.6, 1.6E-5, .0057};
1149 static const std::pair<G4int, const G4double*> Z64N89(89,pZ64N89);
1150 static const G4double pZ64N90[7]={1.7E-7, 12., 8.4E-5, 2.E-9, 1.8, 2.2E-5, .0035};
1151 static const std::pair<G4int, const G4double*> Z64N90(90,pZ64N90);
1152 static const G4double pZ64N91[7]={1.7E-7, 13., 5.E-4, .3E-9, 1.7, 1.9E-5, .026};
1153 static const std::pair<G4int, const G4double*> Z64N91(91,pZ64N91);
1154 static const G4double pZ64N92[7]={1.7E-7, 13., 7.E-5, 2.5E-9, 1.8, 2.2E-5, .003};
1155 static const std::pair<G4int, const G4double*> Z64N92(92,pZ64N92);
1156 static const G4double pZ64N93[7]={1.7E-6, 12., .002, 0., 1.7, 1.8E-5, .47};
1157 static const std::pair<G4int, const G4double*> Z64N93(93,pZ64N93);
1158 static const G4double pZ64N94[7]={3.4E-7, 13., 1.5E-4, 2.E-9, 1.8, 2.3E-5, .0034};
1159 static const std::pair<G4int, const G4double*> Z64N94(94,pZ64N94);
1160 static const G4double pZ64N96[7]={2.6E-6, 13., .0019, 1.2E-9, 1., 1.2E-5, .0056};
1161 static const std::pair<G4int, const G4double*> Z64N96(96,pZ64N96);
1162 static const std::pair<G4int, const G4double*> Z64[N64]={Z64N88, Z64N89, Z64N90, Z64N91,
1163 Z64N92, Z64N93, Z64N94, Z64N96};
1165 static const G4int N65=2;
1166 static const G4double pZ65N94[7]={9.E-7, 16., 3.9E-4, 1.7E-9, 2., 2.2E-5, .0042};
1167 static const std::pair<G4int, const G4double*> Z65N94(94,pZ65N94);
1168 static const G4double pZ65N95[7]={4.5E-7, 16., 1.1E-4, 3.E-9, 1.7, 2.2E-5, .0024};
1169 static const std::pair<G4int, const G4double*> Z65N95(95,pZ65N95);
1170 static const std::pair<G4int, const G4double*> Z65[N65]={Z65N94, Z65N94};
1172 static const G4int N66=7;
1173 static const G4double pZ66N90[7]={1.2E-7, 13., 4.E-5, 3.E-9, 1.2, 1.4E-5, .0025};
1174 static const std::pair<G4int, const G4double*> Z66N90(90,pZ66N90);
1175 static const G4double pZ66N92[7]={1.2E-7, 13., 6.7E-5, 2.E-9, 1., 1.1E-5, .004};
1176 static const std::pair<G4int, const G4double*> Z66N92(92,pZ66N92);
1177 static const G4double pZ66N94[7]={1.2E-7, 13., 5.3E-5, 1.6E-9, 1., 1.1E-5, .0034};
1178 static const std::pair<G4int, const G4double*> Z66N94(94,pZ66N94);
1179 static const G4double pZ66N95[7]={1.2E-6, 13., .0017, .7E-9, 1.3, 1.3E-5, .011};
1180 static const std::pair<G4int, const G4double*> Z66N95(95,pZ66N95);
1181 static const G4double pZ66N96[7]={1.2E-7, 13., 8.E-6, 1.5E-7, 1., 1.1E-5, 1.E-4};
1182 static const std::pair<G4int, const G4double*> Z66N96(96,pZ66N96);
1183 static const G4double pZ66N97[7]={1.5E-7, 13., 4.E-5, 4.E-9, 1.3, 1.3E-5, .002};
1184 static const std::pair<G4int, const G4double*> Z66N97(97,pZ66N97);
1185 static const G4double pZ66N98[7]={3.E-7, 13., .001, 4.E-9, 1.3, 1.3E-5, .23};
1186 static const std::pair<G4int, const G4double*> Z66N98(98,pZ66N98);
1187 static const std::pair<G4int, const G4double*> Z66[N66]={Z66N90, Z66N92, Z66N94, Z66N95,
1188 Z66N96, Z66N97, Z66N98};
1190 static const G4int N67=2;
1191 static const G4double pZ67N98[7]={3.E-7, 13., 2.2E-4, 1.5E-9, 1., 1.E-5, .0054};
1192 static const std::pair<G4int, const G4double*> Z67N98(98,pZ67N98);
1193 static const G4double pZ67N99[7]={7.5E-8, 13., 2.6E-5, 4.5E-9, 1.5, 1.5E-5, .0021};
1194 static const std::pair<G4int, const G4double*> Z67N99(99,pZ67N99);
1195 static const std::pair<G4int, const G4double*> Z67[N67]={Z67N98, Z67N99};
1197 static const G4int N68=6;
1198 static const G4double pZ68N94[7]={1.2E-7, 13., 7.8E-5, 1.6E-9, .9, 9.E-6, .005};
1199 static const std::pair<G4int, const G4double*> Z68N94(94,pZ68N94);
1200 static const G4double pZ68N96[7]={1.2E-7, 13., 8.5E-5, 1.2E-9, .9, 8.E-6, .0055};
1201 static const std::pair<G4int, const G4double*> Z68N96(96,pZ68N96);
1202 static const G4double pZ68N98[7]={1.E-6, 13., .0011, .8E-9, .9, 8.E-6, .0087};
1203 static const std::pair<G4int, const G4double*> Z68N98(98,pZ68N98);
1204 static const G4double pZ68N99[7]={1.2E-7, 13., 2.5E-5, 4.5E-9, .9, 9.E-6, .0015};
1205 static const std::pair<G4int, const G4double*> Z68N99(99,pZ68N99);
1206 static const G4double pZ68N100[7]={2.E-6, 13., .0015, 1.1E-9, .9, 8.E-6, .0058};
1207 static const std::pair<G4int, const G4double*> Z68N100(100,pZ68N100);
1208 static const G4double pZ68N102[7]={2.E-6, 13., .0018, 1.E-9, .9, 8.E-6, .007};
1209 static const std::pair<G4int, const G4double*> Z68N102(102,pZ68N102);
1210 static const std::pair<G4int, const G4double*> Z68[N68]={Z68N94, Z68N96, Z68N98,
1211 Z68N99, Z68N100, Z68N102};
1213 static const G4int N69=1;
1214 static const G4double pZ69N100[7]={1.2E-7, 13., 2.5E-5, 4.5E-9, .9, 9.E-6, .0015};
1215 static const std::pair<G4int, const G4double*> Z69N100(100,pZ69N100);
1216 static const std::pair<G4int, const G4double*> Z69[N69]={Z69N100};
1219 static const G4int N70=7;
1220 static const G4double pZ70N98[7]={1.E-6, 13., .0011, .8E-9, .9, 8.E-6, .0087};
1221 static const std::pair<G4int, const G4double*> Z70N98(98,pZ70N98);
1222 static const G4double pZ70N100[7]={2.E-6, 13., .0015, 1.1E-9, .9, 8.E-6, .0058};
1223 static const std::pair<G4int, const G4double*> Z70N100(100,pZ70N100);
1224 static const G4double pZ70N101[7]={1.2E-7, 13., 2.5E-5, 4.5E-9, .9, 9.E-6, .0015};
1225 static const std::pair<G4int, const G4double*> Z70N101(101,pZ70N101);
1226 static const G4double pZ70N102[7]={2.E-6, 13., .0018, 1.E-9, .9, 8.E-6, .007};
1227 static const std::pair<G4int, const G4double*> Z70N102(102,pZ70N102);
1228 static const G4double pZ70N103[7]={5.E-7, 18., .01, 1.7E-6, 1.2, 1.4E-5, 1.4E-5};
1229 static const std::pair<G4int, const G4double*> Z70N103(103,pZ70N103);
1230 static const G4double pZ70N104[7]={5.E-7, 18., 1.9E-4, 2.5E-9, 1.2, 1.4E-5, .004};
1231 static const std::pair<G4int, const G4double*> Z70N104(104,pZ70N104);
1232 static const G4double pZ70N106[7]={5.E-7, 18., 1.3E-4, 2.E-9, 1.2, 1.4E-5, .0027};
1233 static const std::pair<G4int, const G4double*> Z70N106(106,pZ70N106);
1234 static const std::pair<G4int, const G4double*> Z70[N70]={Z70N98, Z70N100, Z70N101,
1235 Z70N102, Z70N103, Z70N104,
1238 static const G4int N71=2;
1239 static const G4double pZ71N104[7]={5.E-7, 18., 1.8E-4, 2.E-9, .9, 9.E-6, .0036};
1240 static const std::pair<G4int, const G4double*> Z71N104(104,pZ71N104);
1241 static const G4double pZ71N105[7]={2.5E-7, 18., 9.E-5, 1.E-8, .9, 9.E-6, .0016};
1242 static const std::pair<G4int, const G4double*> Z71N105(105,pZ71N105);
1243 static const std::pair<G4int, const G4double*> Z71[N71]={Z71N104, Z71N105};
1245 static const G4int N72=6;
1246 static const G4double pZ72N102[7]={1.E-6, 18., 8.8E-4, .8E-9, 1., 1.1E-5, .0092};
1247 static const std::pair<G4int, const G4double*> Z72N102(102,pZ72N102);
1248 static const G4double pZ72N104[7]={5.E-7, 18., 1.9E-4, 2.5E-9, 1.2, 1.4E-5, .004};
1249 static const std::pair<G4int, const G4double*> Z72N104(104,pZ72N104);
1250 static const G4double pZ72N105[7]={5.E-7, 18., .01, 1.7E-6, 1.2, 1.4E-5, 1.4E-5};
1251 static const std::pair<G4int, const G4double*> Z72N105(105,pZ72N105);
1252 static const G4double pZ72N106[7]={5.E-7, 18., 1.3E-4, 2.E-9, 1.2, 1.4E-5, .0027};
1253 static const std::pair<G4int, const G4double*> Z72N106(106,pZ72N106);
1254 static const G4double pZ72N107[7]={2.5E-7, 18., 1.E-4, 2.E-9, 1.2, 1.5E-5, .0041};
1255 static const std::pair<G4int, const G4double*> Z72N107(107,pZ72N107);
1256 static const G4double pZ72N108[7]={1.E-6, 18., .0012, .6E-9, 1.2, 1.5E-5, .012};
1257 static const std::pair<G4int, const G4double*> Z72N108(108,pZ72N108);
1258 static const std::pair<G4int, const G4double*> Z72[N72]={Z72N102, Z72N104, Z72N105,
1259 Z72N106, Z72N107, Z72N108};
1261 static const G4int N73=2;
1262 static const G4double pZ73N108[7]={5.E-7, 18., 1.7E-4, 2.E-9, 1.2, 1.4E-5, .0035};
1263 static const std::pair<G4int, const G4double*> Z73N108(108,pZ73N108);
1264 static const G4double pZ73N109[7]={1.E-6, 14., .002, .3E-9, 1.3, 1.5E-5, .016};
1265 static const std::pair<G4int, const G4double*> Z73N109(109,pZ73N109);
1266 static const std::pair<G4int, const G4double*> Z73[N73]={Z73N108, Z73N108};
1268 static const G4int N74=5;
1269 static const G4double pZ74N106[7]={5.E-7, 18., 1.3E-4, 2.E-9, 1.2, 1.4E-5, .0027};
1270 static const std::pair<G4int, const G4double*> Z74N106(106,pZ74N106);
1271 static const G4double pZ74N108[7]={4.E-6, 14., .0034, .9E-9, 1.3, 1.4E-5, .0067};
1272 static const std::pair<G4int, const G4double*> Z74N108(108,pZ74N108);
1273 static const G4double pZ74N109[7]={1.2E-7, 14., 3.E-5, 3.E-9, 1.3, 1.4E-5, .0019};
1274 static const std::pair<G4int, const G4double*> Z74N109(109,pZ74N109);
1275 static const G4double pZ74N110[7]={1.2E-7, 14., 3.6E-5, 2.E-9, 1.3, 1.3E-5, .0024};
1276 static const std::pair<G4int, const G4double*> Z74N110(110,pZ74N110);
1277 static const G4double pZ74N112[7]={1.2E-7, 14., 3.E-5, 1.3E-7, 1.3, 1.3E-5, 1.4E-4};
1278 static const std::pair<G4int, const G4double*> Z74N112(112,pZ74N112);
1279 static const std::pair<G4int, const G4double*> Z74[N74]={Z74N106, Z74N108, Z74N109,
1282 static const G4int N75=2;
1283 static const G4double pZ75N110[7]={1.2E-7, 14., 8.E-5, 1.2E-9, 1.3, 1.5E-5, .005};
1284 static const std::pair<G4int, const G4double*> Z75N110(110,pZ75N110);
1285 static const G4double pZ75N112[7]={1.2E-7, 14., 8.8E-5, 1.1E-9, 1.3, 1.5E-5, .0055};
1286 static const std::pair<G4int, const G4double*> Z75N112(112,pZ75N112);
1287 static const std::pair<G4int, const G4double*> Z75[N75]={Z75N110, Z75N112};
1290 static const G4int N76=7;
1291 static const G4double pZ76N108[7]={4.E-6, 14., .0034, .9E-9, 1.3, 1.4E-5, .0067};
1292 static const std::pair<G4int, const G4double*> Z76N108(108,pZ76N108);
1293 static const G4double pZ76N110[7]={1.2E-7, 14., 3.6E-5, 2.E-9, 1.3, 1.3E-5, .0024};
1294 static const std::pair<G4int, const G4double*> Z76N110(110,pZ76N110);
1295 static const G4double pZ76N111[7]={1.2E-7, 14., 8.8E-5, 1.1E-9, 1.3, 1.5E-5, .0055};
1296 static const std::pair<G4int, const G4double*> Z76N111(111,pZ76N111);
1297 static const G4double pZ76N112[7]={1.2E-7, 14., 3.E-5, 1.3E-7, 1.3, 1.3E-5, 1.4E-4};
1298 static const std::pair<G4int, const G4double*> Z76N112(112,pZ76N112);
1299 static const G4double pZ76N113[7]={1.2E-7, 14., 8.8E-5, 1.1E-9, 1.3, 1.5E-5, .0055};
1300 static const std::pair<G4int, const G4double*> Z76N113(113,pZ76N113);
1301 static const G4double pZ76N114[7]={1.2E-7, 14., 3.E-5, 1.3E-7, 1.3, 1.3E-5, 1.4E-4};
1302 static const std::pair<G4int, const G4double*> Z76N114(114,pZ76N114);
1303 static const G4double pZ76N116[7]={1.2E-7, 14., 3.E-5, 1.3E-7, 1.3, 1.3E-5, 1.4E-4};
1304 static const std::pair<G4int, const G4double*> Z76N116(116,pZ76N116);
1305 static const std::pair<G4int, const G4double*> Z76[N76]={Z76N108, Z76N110, Z76N111,
1306 Z76N112, Z76N113, Z76N114,
1309 static const G4int N77=2;
1310 static const G4double pZ77N114[7]={4.8E-7, 14., 5.2E-4, .7E-9, 1.5, 1.7E-5, .0082};
1311 static const std::pair<G4int, const G4double*> Z77N114(114,pZ77N114);
1312 static const G4double pZ77N116[7]={4.8E-7, 14., 4.5E-4, .8E-9, 1.8, 2.3E-5, .0073};
1313 static const std::pair<G4int, const G4double*> Z77N116(116,pZ77N116);
1314 static const std::pair<G4int, const G4double*> Z77[N77]={Z77N114, Z77N116};
1317 static const G4int N78=6;
1318 static const G4double pZ78N112[7]={6.E-8, 19., 3.3E-4, .1E-9, 1.6, 1.8E-5, .06};
1319 static const std::pair<G4int, const G4double*> Z78N112(112,pZ78N112);
1320 static const G4double pZ78N114[7]={6.E-8, 19., 3.3E-4, .1E-9, 1.6, 1.8E-5, .06};
1321 static const std::pair<G4int, const G4double*> Z78N114(114,pZ78N114);
1322 static const G4double pZ78N116[7]={6.E-8, 19., 3.3E-4, .1E-9, 1.6, 1.8E-5, .06};
1323 static const std::pair<G4int, const G4double*> Z78N116(116,pZ78N116);
1324 static const G4double pZ78N117[7]={9.6E-7, 20., .001, .2E-9, 1.6, 2.E-5, .037};
1325 static const std::pair<G4int, const G4double*> Z78N117(117,pZ78N117);
1326 static const G4double pZ78N118[7]={2.4E-7, 20., 1.6E-4, 1.3E-9, 1.6, 1.8E-5, .007};
1327 static const std::pair<G4int, const G4double*> Z78N118(118,pZ78N118);
1328 static const G4double pZ78N120[7]={2.E-6, 19., .0015, .9E-9, 1.6, 1.8E-5, .0078};
1329 static const std::pair<G4int, const G4double*> Z78N120(120,pZ78N120);
1330 static const std::pair<G4int, const G4double*> Z78[N78]={Z78N112, Z78N114, Z78N116,
1331 Z78N117, Z78N118, Z78N120};
1333 static const G4int N79=1;
1334 static const G4double pZ79N118[7]={2.4E-7, 19., 1.E-4, 1.4E-9, 1.3, 1.7E-5, .0042};
1335 static const std::pair<G4int, const G4double*> Z79N118(118,pZ79N118);
1336 static const std::pair<G4int, const G4double*> Z79[N79]={Z79N118};
1338 static const G4int N80=7;
1339 static const G4double pZ80N116[7]={6.E-8, 19., 3.3E-4, .1E-9, 1.6, 1.8E-5, .06};
1340 static const std::pair<G4int, const G4double*> Z80N116(116,pZ80N116);
1341 static const G4double pZ80N118[7]={2.4E-7, 20., 1.6E-4, 1.3E-9, 1.6, 1.8E-5, .007};
1342 static const std::pair<G4int, const G4double*> Z80N118(118,pZ80N118);
1343 static const G4double pZ80N119[7]={9.6E-7, 20., .001, .2E-9, 1.6, 2.E-5, .037};
1344 static const std::pair<G4int, const G4double*> Z80N119(119,pZ80N119);
1345 static const G4double pZ80N120[7]={2.E-6, 19., .0015, .9E-9, 1.6, 1.8E-5, .0078};
1346 static const std::pair<G4int, const G4double*> Z80N120(120,pZ80N120);
1347 static const G4double pZ80N121[7]={1.E-6, 20., 7.E-4, 1.E-9, 1.6, 1.8E-5, .0076};
1348 static const std::pair<G4int, const G4double*> Z80N121(121,pZ80N121);
1349 static const G4double pZ80N122[7]={2.E-6, 18., .0016, .8E-9, 1.6, 1.8E-5, .0078};
1350 static const std::pair<G4int, const G4double*> Z80N122(122,pZ80N122);
1351 static const G4double pZ80N124[7]={2.0E-6, 18., .0032, .4E-9, 1.6, 1.8E-5, .016};
1352 static const std::pair<G4int, const G4double*> Z80N124(124,pZ80N124);
1353 static const std::pair<G4int, const G4double*> Z80[N80]={Z80N116, Z80N118, Z80N119,
1354 Z80N120, Z80N121, Z80N122,
1357 static const G4int N81=2;
1358 static const G4double pZ81N122[7]={2.4E-7, 19., 1.E-4, 1.4E-9, 1.3, 1.7E-5, .0042};
1359 static const std::pair<G4int, const G4double*> Z81N122(122,pZ81N122);
1360 static const G4double pZ81N124[7]={};
1361 static const std::pair<G4int, const G4double*> Z81N124(124,pZ81N124);
1362 static const std::pair<G4int, const G4double*> Z81[N81]={Z81N122, Z81N124};
1364 static const G4int N82=4;
1365 static const G4double pZ82N122[7]={4.E-6, 20., .0022, 1.E-9, 1.6, 1.8E-5, .0058};
1366 static const std::pair<G4int, const G4double*> Z82N122(122,pZ82N122);
1367 static const G4double pZ82N124[7]={4.E-6, 20., .0022, 1.E-9, 1.6, 1.8E-5, .0058};
1368 static const std::pair<G4int, const G4double*> Z82N124(124,pZ82N124);
1369 static const G4double pZ82N125[7]={2.E-6, 20., .0011, 1.2E-9, 1.6, 1.8E-5, .0056};
1370 static const std::pair<G4int, const G4double*> Z82N125(125,pZ82N125);
1371 static const G4double pZ82N126[7]={4.E-6, 20., .0023, 1.2E-9, 1.6, 1.8E-5, .0058};
1372 static const std::pair<G4int, const G4double*> Z82N126(126,pZ82N126);
1373 static const std::pair<G4int, const G4double*> Z82[N82]={Z82N122, Z82N124, Z82N125,
1376 static const G4int N83=1;
1377 static const G4double pZ83N126[7]={8.E-7, 23., 3.3E-4, 1.8E-9, 1.6, 1.8E-5, .005};
1378 static const std::pair<G4int, const G4double*> Z83N126(126,pZ83N126);
1379 static const std::pair<G4int, const G4double*> Z83[N83]={Z83N126};
1381 static const G4int N84=1;
1382 static const G4double pZ84N125[7]={2.E-6, 20., .0011, 1.2E-9, 1.6, 1.8E-5, .0056};
1383 static const std::pair<G4int, const G4double*> Z84N125(125,pZ84N125);
1384 static const std::pair<G4int, const G4double*> Z84[N84]={Z84N125};
1386 static const G4int N85=1;
1387 static const G4double pZ85N125[7]={2.E-6, 20., .0011, 1.2E-9, 1.6, 1.8E-5, .0056};
1388 static const std::pair<G4int, const G4double*> Z85N125(125,pZ85N125);
1389 static const std::pair<G4int, const G4double*> Z85[N85]={Z85N125};
1391 static const G4int N86=1;
1392 static const G4double pZ86N136[7]={1.E-7, 23., 5.5E-5, 1.2E-9, 1.6, 1.8E-5, .0062};
1393 static const std::pair<G4int, const G4double*> Z86N136(136,pZ86N136);
1394 static const std::pair<G4int, const G4double*> Z86[N86]={Z86N136};
1396 static const G4int N87=1;
1397 static const G4double pZ87N136[7]={2.E-7, 23., 1.1E-4, 1.2E-9, 1.6, 1.8E-5, .0062};
1398 static const std::pair<G4int, const G4double*> Z87N136(136,pZ87N136);
1399 static const std::pair<G4int, const G4double*> Z87[N87]={Z87N136};
1401 static const G4int N88=4;
1402 static const G4double pZ88N135[7]={1.E-7, 23., 5.5E-5, 1.2E-9, 1.6, 1.8E-5, .0062};
1403 static const std::pair<G4int, const G4double*> Z88N135(135,pZ88N135);
1404 static const G4double pZ88N136[7]={1.E-7, 23., 5.5E-5, 1.2E-9, 1.6, 1.8E-5, .0062};
1405 static const std::pair<G4int, const G4double*> Z88N136(136,pZ88N136);
1406 static const G4double pZ88N137[7]={1.E-7, 23., 5.5E-5, 1.2E-9, 1.6, 1.8E-5, .0062};
1407 static const std::pair<G4int, const G4double*> Z88N137(137,pZ88N137);
1408 static const G4double pZ88N138[7]={4.E-7, 23., 1.7E-4, 1.5E-9, 1.6, 1.8E-5, .005};
1409 static const std::pair<G4int, const G4double*> Z88N138(138,pZ88N138);
1410 static const std::pair<G4int, const G4double*> Z88[N88]={Z88N135, Z88N136, Z88N137,
1413 static const G4int N89=3;
1414 static const G4double pZ89N136[7]={2.E-7, 23., 1.1E-4, 1.2E-9, 1.6, 1.8E-5, .0062};
1415 static const std::pair<G4int, const G4double*> Z89N136(136,pZ89N136);
1416 static const G4double pZ89N137[7]={4.E-7, 23., 2.2E-4, 1.2E-9, 1.6, 1.8E-5, .0062};
1417 static const std::pair<G4int, const G4double*> Z89N137(137,pZ89N137);
1418 static const G4double pZ89N138[7]={1.E-7, 23., 5.5E-5, 1.2E-9, 1.6, 1.8E-5, .0062};
1419 static const std::pair<G4int, const G4double*> Z89N138(138,pZ89N138);
1420 static const std::pair<G4int, const G4double*> Z89[N89]={Z89N136, Z89N137, Z89N138};
1422 static const G4int N90=7;
1423 static const G4double pZ90N137[7]={4.E-7, 23., 2.2E-4, 1.2E-9, 1.6, 1.8E-5, .0062};
1424 static const std::pair<G4int, const G4double*> Z90N137(137,pZ90N137);
1425 static const G4double pZ90N138[7]={1.E-6, 23., .0016, .4E-9, 3., 3.E-5, .019};
1426 static const std::pair<G4int, const G4double*> Z90N138(138,pZ90N138);
1427 static const G4double pZ90N139[7]={2.5E-7, 23., 1.1E-4, 1.4E-9, 2.4, 2.7E-5, .0049};
1428 static const std::pair<G4int, const G4double*> Z90N139(139,pZ90N139);
1429 static const G4double pZ90N140[7]={1.2E-7, 23., 3.E-5, 2.E-9, 3., 3.E-5, .003};
1430 static const std::pair<G4int, const G4double*> Z90N140(140,pZ90N140);
1431 static const G4double pZ90N142[7]={4.E-6, 23., .0023, 1.1E-9, 1.8, 2.3E-5, .0064};
1432 static const std::pair<G4int, const G4double*> Z90N142(142,pZ90N142);
1433 static const G4double pZ90N143[7]={9.4E-7, 23., 5.4E-4, 1.1E-9, 3., 3.E-5, .0066};
1434 static const std::pair<G4int, const G4double*> Z90N143(143,pZ90N143);
1435 static const G4double pZ90N144[7]={2.5E-7, 23., 1.4E-4, 1.1E-9, 3., 3.E-5, .0066};
1436 static const std::pair<G4int, const G4double*> Z90N144(144,pZ90N144);
1437 static const std::pair<G4int, const G4double*> Z90[N90]={Z90N137, Z90N138, Z90N139,
1438 Z90N140, Z90N142, Z90N143,
1441 static const G4int N91=3;
1442 static const G4double pZ91N140[7]={1.E-5, 23., .0052, 1.6E-9, 1.8, 2.3E-5, .0057};
1443 static const std::pair<G4int, const G4double*> Z91N140(140,pZ91N140);
1444 static const G4double pZ91N141[7]={8.E-6, 23., .006, 0., 3.5, 3.5E-5, .021};
1445 static const std::pair<G4int, const G4double*> Z91N141(141,pZ91N141);
1446 static const G4double pZ91N142[7]={8.E-6, 23., .0042, 1.E-9, 2., 2.5E-5, .006};
1447 static const std::pair<G4int, const G4double*> Z91N142(142,pZ91N142);
1448 static const std::pair<G4int, const G4double*> Z91[N91]={Z91N140, Z91N141, Z91N142};
1450 static const G4int N92=10;
1451 static const G4double pZ92N140[7]={1.4E-6, 20., 8.E-4, 1.5E-9, 2.5, 2.8E-5, .0055};
1452 static const std::pair<G4int, const G4double*> Z92N140(140,pZ92N140);
1453 static const G4double pZ92N141[7]={5.6E-6, 20., .0033, 1.E-9, 2.5, 2.8E-5, .006};
1454 static const std::pair<G4int, const G4double*> Z92N141(141,pZ92N141);
1455 static const G4double pZ92N142[7]={5.6E-6, 20., .0034, 0., 2.5, 2.8E-5, .0072};
1456 static const std::pair<G4int, const G4double*> Z92N142(142,pZ92N142);
1457 static const G4double pZ92N143[7]={5.6E-6, 20., .0032, 0., 2., 2.3E-5, .006};
1458 static const std::pair<G4int, const G4double*> Z92N143(143,pZ92N143);
1459 static const G4double pZ92N144[7]={3.6E-7, 20., 1.6E-4, 1.3E-9, 2.2, 2.7E-5, .0043};
1460 static const std::pair<G4int, const G4double*> Z92N144(144,pZ92N144);
1461 static const G4double pZ92N145[7]={3.6E-6, 20., .003, 0., 2.2, 2.7E-5, .045};
1462 static const std::pair<G4int, const G4double*> Z92N145(145,pZ92N145);
1463 static const G4double pZ92N146[7]={3.6E-7, 20., 1.6E-4, 1.3E-9, 2.2, 2.7E-5, .0043};
1464 static const std::pair<G4int, const G4double*> Z92N146(146,pZ92N146);
1465 static const G4double pZ92N147[7]={3.6E-6, 20., .0014, 1.3E-9, 2.2, 2.7E-5, 12.};
1466 static const std::pair<G4int, const G4double*> Z92N147(147,pZ92N147);
1467 static const G4double pZ92N148[7]={3.4E-7, 20., 1.3E-4, 1.3E-9, 2.2, 2.8E-5, .0036};
1468 static const std::pair<G4int, const G4double*> Z92N148(148,pZ92N148);
1469 static const G4double pZ92N149[7]={3.3E-7, 20., 1.5E-4, 1.2E-9, 3., 3.4E-5, .0044};
1470 static const std::pair<G4int, const G4double*> Z92N149(149,pZ92N149);
1471 static const std::pair<G4int, const G4double*> Z92[N92]={Z92N140, Z92N141, Z92N142,
1472 Z92N143, Z92N144, Z92N145,
1473 Z92N146, Z92N147, Z92N148,
1476 static const G4int N93=5;
1477 static const G4double pZ93N142[7]={3.4E-6, 20., .002, 1.3E-9, 3., 3.3E-5, .0056};
1478 static const std::pair<G4int, const G4double*> Z93N142(142,pZ93N142);
1479 static const G4double pZ93N143[7]={3.4E-6, 20., .002, 1.6E-9, 3.5, 3.6E-5, .005};
1480 static const std::pair<G4int, const G4double*> Z93N143(143,pZ93N143);
1481 static const G4double pZ93N144[7]={6.8E-6, 18., .0052, .8E-9, 2.4, 3.E-5, .0072};
1482 static const std::pair<G4int, const G4double*> Z93N144(144,pZ93N144);
1483 static const G4double pZ93N145[7]={3.4E-6, 20., .002, 1.E-9, 3.5, 3.6E-5, .006};
1484 static const std::pair<G4int, const G4double*> Z93N145(145,pZ93N145);
1485 static const G4double pZ93N146[7]={3.4E-6, 20., .002, 1.5E-9, 3.5, 3.6E-5, .0053};
1486 static const std::pair<G4int, const G4double*> Z93N146(146,pZ93N146);
1487 static const std::pair<G4int, const G4double*> Z93[N93]={Z93N142, Z93N143, Z93N144,
1490 static const G4int N94=10;
1491 static const G4double pZ94N142[7]={6.8E-7, 16., 4.5E-4, 1.7E-9, 2.6, 3.E-5, .0047};
1492 static const std::pair<G4int, const G4double*> Z94N142(142,pZ94N142);
1493 static const G4double pZ94N143[7]={6.8E-6, 18., .0044, .9E-9, 3.3, 3.5E-5, .0058};
1494 static const std::pair<G4int, const G4double*> Z94N143(143,pZ94N143);
1495 static const G4double pZ94N144[7]={6.8E-7, 16., 6.E-4, 0., 2.7, 2.6E-5, .0082};
1496 static const std::pair<G4int, const G4double*> Z94N144(144,pZ94N144);
1497 static const G4double pZ94N145[7]={2.6E-6, 16., .0017, 1.8E-9, 1.8, 2.E-5, .004};
1498 static const std::pair<G4int, const G4double*> Z94N145(145,pZ94N145);
1499 static const G4double pZ94N146[7]={2.5E-7, 20., 9.E-5, 3.6E-8, 3.4, 3.8E-5, 5.4E-4};
1500 static const std::pair<G4int, const G4double*> Z94N146(146,pZ94N146);
1501 static const G4double pZ94N147[7]={1.4E-5, 16., .01, .8E-9, 2.7, 2.6E-5, .0055};
1502 static const std::pair<G4int, const G4double*> Z94N147(147,pZ94N147);
1503 static const G4double pZ94N148[7]={3.4E-7, 20., 1.3E-4, 1.2E-9, 3.2, 3.E-5, .0036};
1504 static const std::pair<G4int, const G4double*> Z94N148(148,pZ94N148);
1505 static const G4double pZ94N149[7]={5.2E-6, 20., .0035, .4E-9, 2.3, 3.E-5, .0095};
1506 static const std::pair<G4int, const G4double*> Z94N149(149,pZ94N149);
1507 static const G4double pZ94N150[7]={3.3E-7, 20., 1.6E-4, 1.2E-9, 3., 3.E-5, .0046};
1508 static const std::pair<G4int, const G4double*> Z94N150(150,pZ94N150);
1509 static const G4double pZ94N152[7]={2.5E-6, 16., .0018, 1.2E-9, 3., 3.1E-5, .0052};
1510 static const std::pair<G4int, const G4double*> Z94N152(152,pZ94N152);
1511 static const std::pair<G4int, const G4double*> Z94[N94]={Z94N142, Z94N143, Z94N144,
1512 Z94N145, Z94N146, Z94N147,
1513 Z94N148, Z94N149, Z94N150,
1516 static const G4int N95=4;
1517 static const G4double pZ95N156[7]={2.5E-6, 18., .0016, .9E-9, 2., 2.3E-5, .0058};
1518 static const std::pair<G4int, const G4double*> Z95N156(156,pZ95N156);
1519 static const G4double pZ95N157[7]={5.E-6, 18., .003, 2.7E-9, 2., 2.3E-5, .0039};
1520 static const std::pair<G4int, const G4double*> Z95N157(157,pZ95N157);
1521 static const G4double pZ95N158[7]={5.E-6, 19., .0033, 2.6E-9, 2., 2.3E-5, .0044};
1522 static const std::pair<G4int, const G4double*> Z95N158(158,pZ95N158);
1523 static const G4double pZ95N159[7]={5.E-5, 20., .029, 1.1E-9, 2., 2.3E-5, .0057};
1524 static const std::pair<G4int, const G4double*> Z95N159(159,pZ95N159);
1525 static const std::pair<G4int, const G4double*> Z95[N95]={Z95N156, Z95N157, Z95N158,
1528 static const G4int N96=10;
1529 static const G4double pZ96N145[7]={5.E-5, 22., .027, 1.1E-9, 2.2, 2.2E-5, .006};
1530 static const std::pair<G4int, const G4double*> Z96N145(145,pZ96N145);
1531 static const G4double pZ96N146[7]={5.E-5, 24., .027, 2.E-9, 2.2, 2.2E-5, .0055};
1532 static const std::pair<G4int, const G4double*> Z96N146(146,pZ96N146);
1533 static const G4double pZ96N147[7]={5.E-5, 22., .025, 2.5E-9, 2.2, 2.4E-5, .0044};
1534 static const std::pair<G4int, const G4double*> Z96N147(147,pZ96N147);
1535 static const G4double pZ96N148[7]={5.E-5, 23., .028, 1.9E-9, 2.2, 3.E-5, .0055};
1536 static const std::pair<G4int, const G4double*> Z96N148(148,pZ96N148);
1537 static const G4double pZ96N149[7]={5.E-5, 23., .025, 1.6E-9, 3., 3.5E-5, .0054};
1538 static const std::pair<G4int, const G4double*> Z96N149(149,pZ96N149);
1539 static const G4double pZ96N150[7]={5.E-5, 24., .026, 2.E-9, 3., 3.6E-5, .0045};
1540 static const std::pair<G4int, const G4double*> Z96N150(150,pZ96N150);
1541 static const G4double pZ96N151[7]={5.E-5, 24., .022, 2.4E-9, 3., 3.6E-5, .0039};
1542 static const std::pair<G4int, const G4double*> Z96N151(151,pZ96N151);
1543 static const G4double pZ96N152[7]={6.5E-7, 25., 2.E-4, 3.4E-9, 3., 3.6E-5, .003};
1544 static const std::pair<G4int, const G4double*> Z96N152(152,pZ96N152);
1545 static const G4double pZ96N153[7]={1.6E-6, 21., 7.E-4, 1.4E-9, 3., 3.6E-5, .0045};
1546 static const std::pair<G4int, const G4double*> Z96N153(153,pZ96N153);
1547 static const G4double pZ96N154[7]={1.3E-5, 16., .016, 0., 3., 3.6E-5, .017};
1548 static const std::pair<G4int, const G4double*> Z96N154(154,pZ96N154);
1549 static const std::pair<G4int, const G4double*> Z96[N96]={Z96N145, Z96N146, Z96N147,
1550 Z96N148, Z96N149, Z96N150,
1551 Z96N151, Z96N152, Z96N153,
1554 static const G4int N97=2;
1555 static const G4double pZ97N152[7]={6.5E-7, 22., 3.5E-4, 2.7E-9, 3., 4.E-5, .004};
1556 static const std::pair<G4int, const G4double*> Z97N152(152,pZ97N152);
1557 static const G4double pZ97N153[7]={6.5E-6, 22., .0036, 1.E-9, 2.7, 4.E-5, .006};
1558 static const std::pair<G4int, const G4double*> Z97N153(153,pZ97N153);
1559 static const std::pair<G4int, const G4double*> Z97[N97]={Z97N152, Z97N153};
1561 static const G4int N98=6;
1562 static const G4double pZ98N151[7]={6.5E-6, 22., .0035, .9E-9, 3., 4.E-5, .0068};
1563 static const std::pair<G4int, const G4double*> Z98N151(151,pZ98N151);
1564 static const G4double pZ98N152[7]={1.3E-6, 22., 7.E-4, 2.E-9, 2.7, 4.E-5, .0045};
1565 static const std::pair<G4int, const G4double*> Z98N152(152,pZ98N152);
1566 static const G4double pZ98N153[7]={2.6E-6, 22., .0014, 2.1E-9, 2.7, 4.E-5, .0044};
1567 static const std::pair<G4int, const G4double*> Z98N153(153,pZ98N153);
1568 static const G4double pZ98N154[7]={2.6E-6, 22., .0014, 1.3E-9, 2.7, 4.E-5, .0054};
1569 static const std::pair<G4int, const G4double*> Z98N154(154,pZ98N154);
1570 static const G4double pZ98N155[7]={2.6E-5, 22., .03, 0., 2.7, 4.E-5, .03};
1571 static const std::pair<G4int, const G4double*> Z98N155(155,pZ98N155);
1572 static const G4double pZ98N156[7]={5.2E-7, 22., 2.6E-4, 1.3E-9, 2.7, 4.E-5, .005};
1573 static const std::pair<G4int, const G4double*> Z98N156(156,pZ98N156);
1574 static const std::pair<G4int, const G4double*> Z98[N98]={Z98N151, Z98N152, Z98N153,
1575 Z98N154, Z98N155, Z98N156};
1577 static const G4int NZ=99;
1578 static const std::pair<G4int, const G4double*>* Pars[NZ]={Z0,Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z8,Z9,
1579 Z10,Z11,Z12,Z13,Z14,Z15,Z16,Z17,Z18,Z19,Z20,Z21,Z22,Z23,Z24,Z25,Z26,Z27,Z28,Z29,Z30,
1580 Z31,Z32,Z33,Z34,Z35,Z36,Z37,Z38,Z39,Z40,Z41,Z42,Z43,Z44,Z45,Z46,Z47,Z48,Z49,Z50,Z51,
1581 Z52,Z53,Z54,Z55,Z56,Z57,Z58,Z59,Z60,Z61,Z62,Z63,Z64,Z65,Z66,Z67,Z68,Z69,Z70,Z71,Z72,
1582 Z73,Z74,Z75,Z76,Z77,Z78,Z79,Z80,Z81,Z82,Z83,Z84,Z85,Z86,Z87,Z88,Z89,Z90,Z91,Z92,Z93,
1583 Z94,Z95,Z96,Z97,Z98};
1584 static const G4int NIso[NZ]={N0,N1,N2,N3,N4,N5,N6,N7,N8,N9,N10,N11,N12,N13,N14,N15,N16,
1585 N17,N18,N19,N20,N21,N22,N23,N24,N25,N26,N27,N28,N29,N30,N31,N32,N33,N34,N35,N36,N37,
1586 N38,N39,N40,N41,N42,N43,N44,N45,N46,N47,N48,N49,N50,N51,N52,N53,N54,N55,N56,N57,N58,
1587 N59,N60,N61,N62,N63,N64,N65,N66,N67,N68,N69,N70,N71,N72,N73,N74,N75,N76,N77,N78,N79,
1588 N80,N81,N82,N83,N84,N85,N86,N87,N88,N89,N90,N91,N92,N93,N94,N95,N96,N97,N98};
1602 if ( tgZ == 1 && tgN == 0 )
1604 for (
G4int ip=0; ip<n_npel; ip++)
lastPAR[ip]=np_el[ip];
1606 else if ( tgZ == 0 && tgN == 1 )
1608 for (
G4int ip=0; ip<n_ppel; ip++)
lastPAR[ip]=pp_el[ip];
1633 lastPAR[ 1]=4.8*std::exp(ala*1.14)/(1.+3.6/
a3);
1634 lastPAR[ 2]=1./(1.+.004*
a4)+2.E-6*a3/(1.+1.3E-6*a3);
1637 lastPAR[ 6]=5.5E-6*std::exp(ala*1.3);
1642 if(nn)
for (
G4int in=0; in<
nn; in++)
1644 std::pair<G4int, const G4double*> curIs=Pars[tgZ][in];
1645 G4int cn=curIs.first;
1662 G4cout<<
"-Warning-G4ChipsNeutronElasticXS::CalcCS: Z="<<tgZ<<
", N="<<tgN
1663 <<
" isotope is not implemented in CHIPS"<<
G4endl;
1679 lastPAR[17]=.7/(1.+4.e-12*a16);
1680 lastPAR[18]=2.5/a8/(a4+1.e-16*a32);
1688 lastPAR[25]=(42.+2.7e-11*a16)/(1.+.14*a);
1692 lastPAR[28]=.0024*a8/(1.+2.6e-4*a7);
1693 lastPAR[29]=3.5e-36*a32*a8/(1.+5.e-15*a32/
a);
1694 lastPAR[30]=1.e5/(a8+2.5e12/a16);
1695 lastPAR[31]=8.e7/(a12+1.e-27*a28*a28);
1704 lastPAR[38]=1.5e-4/(1.+5.e-6*a12);
1712 lastPAR[44]=1.e8*std::exp(.32*asa);
1713 lastPAR[45]=20.*std::exp(.45*asa);
1715 lastPAR[47]=2.5e5*std::exp(.085*a3);
1729 lastPAR[15]=4.5*std::pow(a,1.15);
1730 lastPAR[16]=.06*std::pow(a,.6);
1731 lastPAR[17]=.6*a/(1.+2.e15/a16);
1732 lastPAR[18]=.17/(a+9.e5/a3+1.5e33/a32);
1733 lastPAR[19]=(.001+7.e-11*
a5)/(1.+4.4e-11*a5);
1734 lastPAR[20]=(p1a10*p1a10+2.e-29)/(1.+2.e-22*a12);
1736 lastPAR[21]=400./a12+2.e-22*a9;
1737 lastPAR[22]=1.e-32*a12/(1.+5.e22/a14);
1738 lastPAR[23]=1000./a2+9.5*sa*ssa;
1739 lastPAR[24]=4.e-6*a*asa+1.e11/a16;
1740 lastPAR[25]=(120./a+.002*
a2)/(1.+2.e14/a16);
1747 lastPAR[30]=.0011*asa/(1.+3.e34/a32/
a4);
1748 lastPAR[31]=1.e-5*a2+2.e14/a16;
1749 lastPAR[32]=1.2e-11*a2/(1.+1.5e19/a12);
1750 lastPAR[33]=.016*asa/(1.+5.e16/a16);
1752 lastPAR[34]=.002*a4/(1.+7.e7/std::pow(a-6.83,14));
1753 lastPAR[35]=2.e6/a6+7.2/std::pow(a,.11);
1754 lastPAR[36]=11.*a3/(1.+7.e23/a16/a8);
1757 lastPAR[38]=(.1+4.4e-5*
a2)/(1.+5.e5/a4);
1758 lastPAR[39]=3.5e-4*a2/(1.+1.e8/a8);
1761 lastPAR[42]=1.e-9/a+s4a16*s4a16;
1773 lastPAR[52]=3.e32/(a32+1.e32);
1781 lastPAR[57]=1.e15+2.e27/a4/(1.+2.e-18*a16);
1811 for(
G4int ip=ini; ip<=fin; ip++)
1830 else G4cout<<
"*Warning*G4ChipsNeutronElasticXS::GetPTables: PDG="<<PDG
1831 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
" > fin="<<fin<<
", LP="<<LP
1832 <<
" > ILP="<<ILP<<
" nothing is done!"<<
G4endl;
1834 else G4cout<<
"*Warning*G4ChipsNeutronElasticXS::GetPTables: PDG="<<PDG<<
", Z="
1835 <<tgZ<<
", N="<<tgN<<
", i="<<ini<<
">= max="<<
nPoints<<
", LP="<<LP
1836 <<
" > ILP="<<ILP<<
", lPMax="<<
lPMax<<
" nothing is done!"<<
G4endl;
1845 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
1846 <<
", while it is defined only for PDG=2112 (n)" <<
G4endl;
1847 G4Exception(
"G4ChipsNeutronElasticXS::GetPTables()",
"HAD_CHPS_0000",
1860 if(PDG!=2112)
G4cout<<
"*Warning*G4ChipsNeutronElasticXS::GetExT:PDG="<<PDG<<
G4endl;
1864 if(tgZ==1 && tgN==0)
1880 q2=-std::log(1.-ran)/
theB1;
1900 if(a>6.5)E3*=tm2*tm2*tm2;
1915 q2=-std::log(1.-ran)/
theB1;
1916 if(std::fabs(tss)>1.e-7) q2=(std::sqrt(
theB1*(
theB1+(tss+tss)*q2))-
theB1)/tss;
1922 q2=-std::log(1.-ran)/
theB2;
1924 if(a<6.5) q2=std::pow(q2,third);
1925 else q2=std::pow(q2,fifth);
1931 q2=-std::log(1.-ran)/
theB3;
1933 if(a>6.5) q2=std::pow(q2,sevth);
1939 q2=-std::log(1.-ran)/
theB4;
1940 if(a<6.5) q2=lastTM-q2;
1944 if(!(q2>=-1.||q2<=1.))
G4cout<<
"*NAN*G4QNeutronElCroSect::GetExchangeT: -t="<<q2<<
G4endl;
1958 if(
lastLP<-4.3)
return 0.;
1965 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
1966 <<
", while it is defined only for PDG=2112 (n) " <<
G4endl;
1967 G4Exception(
"G4ChipsNeutronElasticXS::GetSlope()",
"HAD_CHPS_0000",
1986 if(PDG!=2112)
G4cout<<
"*Warning*G4ChipsNeutronElasticXS::GetTaV:PDG="<<PDG<<
G4endl;
1989 G4cout<<
"*Warning*G4QNElasticCrS::GetTabValue: (1-92) No isotopes for Z="<<tgZ<<
G4endl;
2004 if ( tgZ == 1 && tgN == 0 )
2010 theS1=(lastPAR[9]+lastPAR[10]*dl1*dl1+lastPAR[11]/p)/(1.+lastPAR[12]/p4)
2011 +lastPAR[13]/(p4+lastPAR[14]);
2012 theB1=(lastPAR[17]+lastPAR[18]/(p4*p4+lastPAR[19]*p3))/(1.+lastPAR[20]/p4);
2013 theS2=(lastPAR[15]+lastPAR[16]/p4/p)/p3;
2014 theB2=lastPAR[22]/(p*sp+lastPAR[23]);
2020 return lastPAR[0]/(p2s+lastPAR[1]*p+lastPAR[2]/ssp)+lastPAR[4]/p
2021 +(lastPAR[5]+lastPAR[6]*dl1*dl1+lastPAR[7]/p)/(1.+lastPAR[8]/p4);
2048 lastPAR[46]/(1.+lastPAR[47]*std::pow(p,lastPAR[48])));
2049 theB4=lastPAR[49]*pa/p2/(1.+pa*lastPAR[50]);
2084 static const G4double mNeut2= mNeut*mNeut;
2087 if(tgZ==0 && tgN==1)
2089 G4double tMid=std::sqrt(pP2+mNeut2)*mNeut-mNeut2;
2097 G4double mds=dmt*std::sqrt(pP2+mNeut2)+mNeut2+mt*mt;
2098 return dmt*dmt*pP2/mds;
2106 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N =" << tgN
2107 <<
", while it is defined only for n projectiles & Z_target>0" <<
G4endl;
2108 G4Exception(
"G4ChipsNeutronElasticXS::GetQ2max()",
"HAD_CHPS_0000",
G4double GetTabValues(G4double lp, G4int pPDG, G4int tgZ, G4int tgN)
std::ostringstream G4ExceptionDescription
std::vector< G4double * > PAR
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
std::vector< G4double * > SST
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
virtual G4bool IsIsoApplicable(const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
std::vector< G4double * > B3T
G4ChipsNeutronElasticXS()
G4ParticleDefinition * GetDefinition() const
std::vector< G4double * > B2T
std::vector< G4double * > B1T
G4double GetTotalMomentum() const
std::vector< G4double * > S3T
G4IonTable * GetIonTable() const
G4GLOB_DLL std::ostream G4cout
G4double GetPTables(G4double lpP, G4double lPm, G4int PDG, G4int tZ, G4int tN)
std::vector< G4double * > B4T
static G4Proton * Proton()
static G4Neutron * Neutron()
G4_DECLARE_XS_FACTORY(G4ChipsNeutronElasticXS)
static const G4double A[nN]
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
std::vector< G4double * > S2T
std::vector< G4double * > S4T
G4double GetQ2max(G4int pPDG, G4int tgZ, G4int tgN, G4double pP)
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
std::vector< G4double * > CST
static const double gigaelectronvolt
static const double millibarn
~G4ChipsNeutronElasticXS()
G4double GetSlope(G4int tZ, G4int tN, G4int pPDG)
G4double CalculateCrossSection(G4bool CS, G4int F, G4int I, G4int pPDG, G4int Z, G4int N, G4double pP)
static const G4double pos
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
std::vector< G4double * > S1T