#include <G4UPiNuclearCrossSection.hh>

Inheritance diagram for G4UPiNuclearCrossSection:

Public Member Functions
	G4UPiNuclearCrossSection ()

virtual	~G4UPiNuclearCrossSection ()

virtual G4bool	IsElementApplicable (const G4DynamicParticle aParticle, G4int Z, const G4Material )

G4double	GetElasticCrossSection (const G4DynamicParticle *aParticle, G4int Z, G4int A)

G4double	GetInelasticCrossSection (const G4DynamicParticle *aParticle, G4int Z, G4int A)

void	BuildPhysicsTable (const G4ParticleDefinition &)

void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) const

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

virtual G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element elm=0, const G4Material *mat=0)

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat=0)

virtual G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=0, const G4Element elm=0, const G4Material mat=0)

virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

const G4String &	GetName () const

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 47 of file G4UPiNuclearCrossSection.hh.

Constructor & Destructor Documentation

G4UPiNuclearCrossSection::G4UPiNuclearCrossSection ( )

Definition at line 47 of file G4UPiNuclearCrossSection.cc.

References G4NistManager::GetAtomicMassAmu(), G4NistManager::Instance(), python.hepunit::MeV, G4PionMinus::PionMinus(), and G4PionPlus::PionPlus().

  : G4VCrossSectionDataSet("G4UPiNuclearCrossSection")
 {
   isInitialized = false;
   piPlusElastic = piPlusInelastic = piMinusElastic = piMinusInelastic = 0;
   piPlus  = G4PionPlus::PionPlus();
   piMinus = G4PionMinus::PionMinus();
 
   NZ = 16;
   aPower  = 0.75;
   elow    = 20.0*MeV;
   elowest = MeV;
   G4NistManager* nist = G4NistManager::Instance();
   for(G4int i=1; i<93; ++i) {
     APower[i] = std::pow(nist->GetAtomicMassAmu(i),aPower);
   }
 }

G4UPiNuclearCrossSection::~G4UPiNuclearCrossSection ( )

virtual

Definition at line 65 of file G4UPiNuclearCrossSection.cc.

 {
   delete piPlusElastic;
   delete piPlusInelastic;
   delete piMinusElastic;
   delete piMinusInelastic;
 }

Member Function Documentation

void G4UPiNuclearCrossSection::BuildPhysicsTable ( const G4ParticleDefinition & p )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 192 of file G4UPiNuclearCrossSection.cc.

References G4NistManager::GetAtomicMassAmu(), G4NistManager::Instance(), and n.

Referenced by G4BGGPionElasticXS::BuildPhysicsTable(), and G4BGGPionInelasticXS::BuildPhysicsTable().

 {
   if(isInitialized) { return; }
   if(&p != piPlus && &p != piMinus) { 
     throw G4HadronicException(__FILE__, __LINE__,"Is applicable only for pions");
     return;
   }
   isInitialized = true;
 
   const G4int n = 16;
   const G4int iz[n] = {2,4,6,7,8,11,13,20,26,29,42,48,50,74,82,92};
   NZ = n;
   theZ.reserve(n);
   theA.reserve(n);
 
   G4NistManager* nist = G4NistManager::Instance();
   G4int i;
   for(i=0; i<n; ++i) {
     theZ.push_back(iz[i]);
     theA.push_back(nist->GetAtomicMassAmu(iz[i]));
   }
 
   piPlusElastic    = new G4PhysicsTable();
   piPlusInelastic  = new G4PhysicsTable();
   piMinusElastic   = new G4PhysicsTable();
   piMinusInelastic = new G4PhysicsTable();
 
   const G4double e1[38] = {
     0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.13, 0.14, 0.15, 0.16, 
     0.17, 0.18, 0.19, 0.2,  0.22,0.24, 0.26, 0.28, 0.3,  0.35, 
     0.4,  0.45, 0.5,  0.55, 0.6, 0.7,  0.8,  0.9,  1,    2, 
     3,    5,    10,   20,   50,  100,  500, 1000};
   const G4double e2[39] = {
     0.02, 0.04, 0.06, 0.08, 0.1,  0.11, 0.12, 0.13, 0.14, 0.15, 
     0.16, 0.17, 0.18, 0.2,  0.22, 0.24, 0.26, 0.28, 0.3,  0.35, 
     0.4,  0.45, 0.5,  0.55, 0.575,0.6,  0.7,  0.8,  0.9,  1, 
     2,    3,    5,   10,   20,    50,   100,  500,  1000};
   const G4double e3[31] = {
     0.02, 0.04, 0.06, 0.08, 0.1,  0.12, 0.14, 0.16, 0.18, 0.2, 
     0.22, 0.25, 0.3,  0.35, 0.4,  0.45, 0.5,  0.6,  0.7,  0.8, 
     0.9,  1,    2,    3,    5,    10,    20,   50,   100, 500, 1000};
   const G4double e4[32] = {
     0.02, 0.04, 0.06, 0.08, 0.1,  0.12, 0.14, 0.16, 0.18,  0.2, 
     0.22, 0.25, 0.3,  0.35, 0.4,  0.45, 0.5, 0.55,  0.6,   0.7,  
     0.8,  0.9,    1,    2,    3,    5,   10,   20,   50,  100,  500,  1000};
   const G4double e5[34] = {
    0.02,  0.04, 0.05, 0.06, 0.07, 0.08, 0.09,  0.1, 0.12, 0.14, 
    0.16, 0.18,  0.2,  0.22, 0.25, 0.3,  0.35,  0.4, 0.45,  0.5,  
    0.6,  0.7,   0.8,  0.9,    1,    2,    3,   5,   10,    20, 50, 100, 500, 1000};
   const G4double e6[35] = {
    0.02,  0.04, 0.05, 0.06, 0.07, 0.08, 0.09,  0.1, 0.12, 0.14, 
    0.16, 0.18,  0.2,  0.22, 0.25,  0.3, 0.35,  0.4, 0.45,  0.5, 
    0.55,  0.6,  0.7,  0.8,  0.9,    1,    2,    3,    5,   10, 20, 50, 100, 500, 1000};
 
   const G4double he_t[38] = {
    40,  70, 108,   152, 208,   276, 300, 320,   329, 333, 
   332, 328, 322,   310, 288,   260, 240, 216,   196, 144, 
   125, 112, 108.5, 109, 110.5, 117, 123, 128.5, 135, 110, 
    96,  87,  85,  83.5, 83.5, 83.5, 83.5, 83.5};
   const G4double he_in[38] = {
    18,  38,  62,    98, 136,   176, 190, 200,   209, 212, 
   212, 208, 204,   196, 176,   164, 150, 134,   124, 97.5, 
    90,  85,82.5,  83.5, 86.5,   93, 97.5,100,   102, 83, 
    77,  75,  74,  72.5, 72.5, 72.5, 72.5, 72.5};
   const G4double be_m_t[38] = {
   150, 210, 294,   396, 520,   600, 623, 635,   642, 640, 
   630, 615, 600,   576, 540,   504, 470, 435,   400, 340, 
   294, 258, 236,   230, 233,   244, 257, 270,   276, 250, 
   230, 215, 205,   194, 188,   186, 186, 186};
   const G4double be_m_in[38] = {
    90, 126, 177,   240, 320,   380, 400, 410,   414, 410, 
   400, 387, 371,   360, 333,   312, 285, 260,   237, 216, 
   198, 187, 182,   180, 182,   187, 193, 203,   207, 179, 
   172, 165, 159,   155, 144,   144, 144, 144};
   const G4double be_p_t[38] = {
    96, 150, 222,   320, 430,   514, 545, 565,   574, 574, 
   564, 552, 535,   522, 490,   462, 432, 398,   367, 314, 
    276, 248, 232,  230, 233,   244, 257, 270,   276, 250, 
   230, 215, 205,   194, 188,   186, 186, 186};
   const G4double be_p_in[38] = {
    60,  95, 142,   194, 262,   319, 345, 361,   364, 364, 
   354, 350, 330,   319, 298,   280, 258, 237,   216, 200, 
   189, 183, 182,   180, 182,   187, 193, 203,   207, 179, 
   172, 165, 159,   155, 144,   144, 144, 144};
 
   const G4double c_m_t[39] = {
   204, 260, 366,   517, 630,   673, 694, 704, 710, 711, 
   706, 694, 676,   648, 616,   584, 548, 518, 489, 426, 
   376, 342, 323,   310, 312,   313, 319, 333, 342, 348, 
   310, 290, 268,   250, 245,   237, 234, 234, 234};
   const G4double c_m_in[39] = {
   128, 160, 224,   315, 388,   416, 430, 438, 444, 445, 
   440, 432, 416,   400, 380,   354, 320, 304, 288, 264, 
   246, 240, 233,   232, 233,   234, 238, 246, 252, 256, 
   220, 210, 198,   187, 183,   176, 174, 174, 174};
   const G4double c_p_t[39] = {
   140, 192, 294,   428, 594,   642, 662, 687, 685, 688, 
   684, 672, 656,   630, 598,   567, 533, 504, 474, 416, 
   369, 336, 319,   310, 312,   313, 319, 333, 342, 348, 
   310, 290, 268,   250, 245,   237, 234, 234, 234};
   const G4double c_p_in[39] = {
    94, 132, 184,   260, 370,   398, 408, 420, 426, 428, 
   424, 416, 400,   386, 366,   340, 308, 294, 280, 257, 
   241, 236, 231,   232, 233,   234, 238, 246, 252, 256, 
   220, 210, 198,   187, 183,   176, 174, 174, 174};
   const G4double n_m_t[39] = {
   246, 308, 424,   590, 729,   776, 800, 821, 822, 817, 
   800, 778, 768,   728, 690,   654, 615, 584, 556, 480, 
   430, 393, 373,   367, 368,   370, 375, 388, 390, 397, 
   364, 337, 310,   291, 275,   268, 268, 268, 268};
   const G4double n_m_in[39] = {
   155, 188, 256,   360, 456,   492, 512, 526, 526, 520, 
   504, 491, 475,   450, 425,   396, 376, 360, 340, 300, 
   282, 270, 265,   265, 266,   268, 273, 280, 288, 288, 
   256, 237, 226,   218, 208,   202, 202, 202, 202};
   const G4double n_p_t[39] = {
   150, 212, 328,   500, 680,   735, 762, 781, 782, 779, 
   770, 748, 740,   706, 672,   633, 600, 569, 541, 467, 
   419, 385, 368,   364, 366,   368, 375, 388, 390, 397, 
   364, 337, 310,   291, 275,   268, 268, 268, 268};
   const G4double n_p_in[39] = {
    90, 140, 208,   300, 426,   467, 490, 504, 504, 500, 
   484, 474, 460,   437, 413,   381, 365, 350, 330, 292, 
   276, 267, 263,   264, 265,   267, 273, 280, 288, 288, 
   256, 237, 226,   218, 208,   202, 202, 202, 202};
 
   const G4double o_m_t[31] = {
   280, 360, 500,   685, 812,   861, 870, 865, 835, 800, 
   755, 700, 600,   537, 493,   468, 441, 436, 443, 449, 
   460, 463, 432,   385, 350,   325, 312, 307, 303, 303, 303};
   const G4double o_m_in[31] = {
   190, 207, 300,   420, 500,   540, 550, 542, 520, 490, 
   460, 423, 360,   339, 321,   314, 312, 314, 319, 324, 
   328, 330, 300,   275, 250,   240, 229, 225, 222, 222, 222};
   const G4double o_p_t[31] = {
   170, 240, 390,   570, 740,   818, 830, 822, 800, 765, 
   725, 675, 585,   525, 483,   458, 444, 447, 453, 449,
   460, 463, 432,   385, 350,   325, 312, 307, 303, 303, 303};
   const G4double o_p_in[31] = {
   100, 145, 240,   340, 470,   518, 530, 522, 505, 477, 
   448, 412, 350,   330, 316,   310, 308, 311, 317, 324,
   328, 330, 300,   275, 250,   240, 229, 225, 222, 222, 222};
   const G4double na_m_t[31] = {
   450, 545, 705,   910, 1020, 1075, 1087, 1080, 1042, 987, 
   943, 885, 790,   700,  650,  610,  585,  575,  585, 595, 
   600, 610, 556,   524,  494,  458,  445,  429,  427, 427, 427};
   const G4double na_m_in[31] = {
   275, 315, 413,   545,  620,  660,  670,  662,  630, 593, 
   570, 520, 465,   420,  410,  395,  390,  400,  410, 418, 
   420, 422, 372,   348,  330,  320,  310,  294,  292, 292, 292};
   const G4double na_p_t[31] = {
   210, 320, 530,   795,  960, 1035, 1050, 1040, 1007, 957, 
   918, 865, 773,   685,  636,  598,  575,  565,  578, 590, 
   598, 610, 556,   524,  494,  458,  445,  429,  427, 427, 427};
   const G4double na_p_in[31] = {
   115, 210, 340,   495,  585,  630,  645,  637,  605, 572, 
   550, 505, 455,   410,  401,  388,  383,  393,  405, 414, 
   418, 422, 372,   348,  330,  320,  310,  294,  292, 292, 292};
   const G4double al_m_t[31] = {
   532, 637, 832,  1057, 1207, 1230, 1210, 1174, 1133, 1095, 
  1038, 970, 890,   807,  750,  710,  675,  665,  670,  673, 
   678, 682, 618,   574,  546,  520,  507,  495,  488,  488, 488};
   const G4double al_m_in[31] = {
   300, 360, 495,   665,  750,  765,  750,  730,  700,  660, 
   615, 570, 520,   490,  470,  450,  448,  450,  450,  452, 
   456, 460, 408,   392,  376,  356,  347,  338,  332,  332, 332};
   const G4double al_p_t[31] = {
   225, 350, 616,   945, 1122, 1175, 1157, 1128, 1088, 1045, 
   988, 935, 870,   787,  730,  690,  660,  652,  660,  668, 
   678, 682, 618,   574,  546,  520,  507,  495,  488,  488, 488};
   const G4double al_p_in[31] = {
   120, 238, 390,   610,  712,  735,  720,  703,  655,  635, 
   590, 550, 505,   475,  455,  438,  440,  445,  445,  450, 
   456, 460, 408,   392,  376,  356,  347,  338,  332,  332, 332};
 
   const G4double ca_m_t[31] = {
   800, 980, 1240, 1460, 1570, 1600, 1580, 1535, 1475, 1425, 
  1375,1295, 1200, 1083, 1000,  948,  915,  895,  900,  908, 
   915, 922, 856,   795,  740,  705,  682,  660,  660,  660, 660};
   const G4double ca_m_in[31] = {
   470, 550, 620,   860,  955,  960,  920,  860,  820,  780, 
   740, 665, 637,   615,  600,  590,  580,  580,  600,  608,  
   610, 615, 550,   525,  510,  488,  470,  450,  450,  450,  450};
   const G4double ca_p_t[31] = {
   275, 445, 790,  1195, 1440, 1485, 1475, 1435, 1385, 1335, 
  1295,1245,1160,  1050,  970,  923,  895,  877,  887,  897, 
   904, 913, 855,   795,  740,  705,  682,  660,  660,  660, 660};
   const G4double ca_p_in[31] = {
   160, 315, 500,   745,  870,  905,  900,  860,  810,  770, 
   740, 710, 640,   617,  595,  585,  575,  575,  590,  600, 
   602, 608, 550,   525,  510,  488,  470,  450,  450,  450,  450};
   // ca data may have typo
 
   const G4double fe_m_t[32] = {
  1175, 1363, 1670, 1950, 2050, 2040, 1975, 1886, 1834, 1773, 
  1720, 1635, 1474, 1380, 1269, 1225, 1182, 1162, 1159, 1162, 
  1178, 1190, 1197, 1102, 1135,  975,  945,  925,  905,  905,  
   905,  905};
   const G4double fe_m_in[32] = {
   625,  725,  910, 1180, 1275, 1250, 1200, 1150, 1100, 1040,  
   995,  925,  825,  810,  780,  760,  745,  740,  740,  740,  
   750,  760,  765,  690,  660,  635,  615,  600,  585,  585,  
   585,  585};
   const G4double fe_p_t[32] =  {
   330,  575, 1010, 1500, 1837, 1875, 1820, 1751, 1691, 1636, 
  1690, 1450, 1396, 1305, 1219, 1190, 1148, 1138, 1134, 1144, 
   1163, 1175, 1183, 1198, 1135, 975,  945,  925,  905,  905,  
   905,  905};
   const G4double fe_p_in[32] = {
   210,  410,  707, 1010, 1125, 1150, 1100, 1070, 1010,  960,  
   920,  776,  780,  760,  750,  740,  720,  725,  725,  730,  
   740,  750,  755,  690,  660,  635,  615,  600,  585,  585,  
   585,  585};
   const G4double cu_m_t[32] = {
  1400, 1600, 1875, 2088, 2200, 2220, 2175, 2125, 2075, 2012, 
  1950, 1855, 1670, 1530, 1430, 1370, 1315, 1315, 1315, 1330, 
  1345, 1360, 1365, 1250, 1185, 1128, 1070, 1035, 1010, 1010, 
  1010, 1010};
   const G4double cu_m_in[32] = {
   725,  840, 1020, 1200, 1295, 1300, 1267, 1240, 1213, 1175, 
  1125, 1042,  950,  900,  860,  840,  830,  832,  835,  840,  
   850,  860,  865,  785,  735,  705,  680,  650,  630,  630,  
   630,  630};
   const G4double cu_p_t[32] = {
   355,  605, 1120, 1630, 1940, 2010, 2010, 1980, 1925, 1895, 
  1830, 1730, 1585, 1490, 1400, 1340, 1290, 1290, 1290, 1310, 
  1330, 1345, 1350, 1240, 1185, 1128, 1070, 1035, 1010, 1010, 
  1010, 1010};
   const G4double cu_p_in[32] = {
   230,  425,  780,  1025, 1155, 1190, 1190, 1180, 1125, 1100, 
  1050, 1000,  900,   870,  835,  815,  810,  812,  815,  825,  
   840,  850,  855,   780,  735,  705,  680,  650,  630,  630,  
   630,  630};
 
   const G4double mo_m_t[34] = {
   2430, 2610, 2710, 2790, 2880, 2940, 2965, 2970, 2970, 2920, 
   2840, 2720, 2570, 2500, 2365, 2200, 2050, 1926, 1825, 1768, 
   1749, 1750, 1778, 1789, 1808, 1690, 1645, 1530, 1492, 1450, 
   1425, 1425, 1425, 1425};
   const G4double mo_m_in[34] = {
    925, 1125, 1250, 1375, 1500, 1600, 1680, 1750, 1770, 1730, 
   1660, 1580, 1500, 1450, 1330, 1250, 1190, 1140, 1100, 1075, 
   1075, 1070, 1088, 1095, 1110, 1035, 1005,  940,  917,  880, 
    860,  860,  860, 860};
   const G4double mo_p_t[34] = {
    410,  730, 1110, 1530, 1920, 2200, 2385, 2520, 2600, 2630, 
   2575, 2470, 2320, 2285, 2185, 2053, 1945, 1852, 1776, 1719, 
   1710, 1716, 1746, 1759, 1778, 1675, 1645, 1530, 1492, 1450, 
   1425, 1425, 1425, 1425};
   const G4double mo_p_in[34] = {
    270,  540,  825,  975, 1140, 1285, 1400, 1480, 1555, 1580, 
   1525, 1470, 1360, 1340, 1255, 1160, 1120, 1085, 1060, 1045, 
   1045, 1045, 1065, 1075, 1090, 1025, 1005,  940,  917,  880, 
    860,  860,  860, 860};
   const G4double cd_m_t[34] = {
   3060, 3125, 3170, 3220, 3255, 3280, 3290, 3260, 3270, 3200, 
   3120, 3080, 3090, 2920, 2810, 2640, 2362, 2230, 2115, 2050, 
   2020, 2025, 2040, 2070, 2100, 1900, 1795, 1740, 1675, 1645, 
   1625, 1620, 1620, 1620};
   const G4double cd_m_in[34]= {
   1025, 1275, 1440, 1625, 1740, 1800, 1880, 1920, 1980, 1920, 
   1850, 1810, 1720, 1650, 1560, 1450, 1330, 1290, 1245, 1210, 
   1200, 1200, 1205, 1205, 1230, 1130, 1085, 1060, 1000,  985, 
    975,  970,  970, 970};
   const G4double cd_p_t[34] = {
    455,  780, 1170, 1700, 2120, 2400, 2600, 2720, 2820, 2840, 
   2800, 2760, 2720, 2640, 2560, 2450, 2252, 2130, 2035, 1985, 
   1970, 1975, 2005, 2035, 2070, 1880, 1795, 1740, 1675, 1645, 
   1625, 1620, 1620, 1620};
   const G4double cd_p_in[34] = {
    310,  580,  880, 1060, 1270, 1400, 1530, 1610, 1660, 1680, 
   1640, 1600, 1560, 1500, 1430, 1330, 1280, 1230, 1200, 1180, 
   1170, 1175, 1180, 1180, 1210, 1120, 1085, 1060, 1000,  985, 
    975,  970,  970, 970};
 
   const G4double sn_m_t[35] = {
   3000, 3180, 3250, 3300, 3300, 3410, 3470, 3450, 3410, 3350, 
   3280, 3200, 3120, 3050, 2900, 2630, 2500, 2325, 2190, 2100, 
   2060, 2055, 2055, 2055, 2067, 2085, 2000, 1900, 1835, 1770, 
   1720, 1700, 1695, 1695, 1695};
   const G4double sn_m_in[35] = {
   1050, 1350, 1520, 1650, 1800, 1980, 2070, 2120, 2090, 2050, 
   1980, 1920, 1830, 1770, 1670, 1500, 1435, 1350, 1300, 1230, 
   1220, 1235, 1235, 1235, 1237, 1240, 1160, 1120, 1090, 1065, 
   1040, 1020, 1015, 1015, 1015};
   const G4double sn_p_t[35] =  { 
    465,  800, 1200, 1760, 2170, 2480, 2730, 2885, 2970, 2980, 
   2970, 2890, 2840, 2790, 2620, 2450, 2335, 2205, 2080, 2020, 
   2010, 1990, 1990, 2015, 2030, 2045, 1980, 1890, 1835, 1770, 
   1720, 1700, 1695, 1695, 1695};
   const G4double sn_p_in[35] = { 
    315,  590,  880, 1220, 1460, 1580, 1700, 1770, 1810, 1810, 
   1800, 1730, 1680, 1630, 1530, 1400, 1335, 1270, 1210, 1180, 
   1190, 1190, 1190, 1205, 1210, 1210, 1150, 1115, 1090, 1065, 
   1040, 1020, 1015, 1015, 1015};
   const G4double w_m_t[35] = {
   5200, 5115, 5025, 4975, 4900, 4850, 4780, 4725, 4600, 4490, 
   4355, 4255, 4125, 4040, 3830, 3580, 3330, 3110, 2955, 2860, 
   2852, 2845, 2885, 2900, 2915, 2940, 2800, 2660, 2570, 2490, 
   2460, 2425, 2420, 2420, 2420};
   const G4double w_m_in[35] = {
   1450, 1850, 2100, 2350, 2550, 2700, 2825, 2900, 2850, 2750, 
   2630, 2525, 2400, 2300, 2200, 2070, 1880, 1770, 1715, 1680, 
   1680, 1680, 1685, 1690, 1700, 1720, 1635, 1560, 1530, 1460, 
   1440, 1410, 1410, 1410, 1410};
   const G4double w_p_t[35] = { 
    480,  900, 1500, 2350, 3020, 3420, 3650, 3775, 3875, 3830, 
   3750, 3700, 3630, 3550, 3550, 3290, 3070, 2890, 2840, 2730, 
   2725, 2720, 2770, 2805, 2828, 2865, 2770, 2640, 2570, 2490,
   2460, 2425, 2420, 2420, 2420};
   const G4double w_p_in[35] = { 
    325,  680,  990, 1500, 1850, 2150, 2250, 2300, 2350, 2330, 
   2280, 2230, 2200, 2120, 2130, 1900, 1780, 1670, 1635, 1600, 
   1602, 1605, 1610, 1615, 1630, 1660, 1620, 1550, 1530, 1460,
   1440, 1410, 1410, 1410, 1410};
 
   const G4double pb_m_t[35] = {
   5890, 5700, 5610, 5580, 5550, 5480, 5400, 5300, 5100, 4930, 
   4750, 4600, 4400, 4280, 4170, 3915, 3650, 3470, 3260, 3150, 
   3120, 3070, 3085, 3100, 3120, 3160, 3070, 2930, 2820, 2750, 
   2710, 2655, 2640, 2640, 2640};
   const G4double pb_m_in[35] = {
   1575, 2025, 2300, 2575, 2850, 3000, 3115, 3180, 3080, 2940, 
   2800, 2670, 2550, 2450, 2370, 2220, 2110, 2000, 1920, 1880, 
   1850, 1800, 1805, 1810, 1820, 1840, 1800, 1720, 1640, 1620, 
   1570, 1530, 1530, 1530, 1530};
   const G4double pb_p_t[35] = { 
    515,  940, 1500, 2400, 3270, 3750, 4050, 4140, 4260, 4200, 
   4080, 3990, 3990, 3810, 3730, 3520, 3370, 3186, 3110, 3010, 
   2990, 2985, 3005, 3020, 3040, 3080, 3020, 2905, 2790, 2750,
   2710, 2655, 2640, 2640, 2640};
   const G4double pb_p_in[35] = { 
    348,  707, 1040, 1650, 2100, 2400, 2580, 2640, 2650, 2520, 
   2410, 2300, 2250, 2190, 2130, 2000, 1930, 1870, 1830, 1790, 
   1770, 1765, 1775, 1780, 1790, 1800, 1775, 1710, 1620, 1620,
   1570, 1530, 1530, 1530, 1530};
   const G4double u_m_t[35] = {
   7080, 6830, 6650, 6530, 6400, 6280, 6100, 5840, 5660, 5520, 
   5330, 5160, 4990, 4810, 4630, 4323, 4130, 3870, 3700, 3550, 
   3490, 3465, 3467, 3475, 3495, 3515, 3440, 3360, 3150, 3040, 
   2985, 2955, 2940, 2940, 2940};
   const G4double u_m_in[35] = {
   1740, 2220, 2500, 2820, 3080, 3300, 3420, 3500, 3420, 3330, 
   3200, 3060, 2940, 2850, 2710, 2470, 2380, 2250, 2160, 2080, 
   2040, 2045, 2047, 2050, 2055, 2060, 2010, 1980, 1830, 1780, 
   1735, 1710, 1700, 1700, 1700};
   const G4double u_p_t[35] = { 
    485,  960, 1580, 2700, 3550, 4050, 4320, 4420, 4620, 4660, 
   4580, 4470, 4350, 4295, 4187, 3938, 3755, 3573, 3450, 3342, 
   3310, 3295, 3310, 3330, 3375, 3405, 3350, 3338, 3135, 3040,
   2985, 2955, 2940, 2940, 2940};
   const G4double u_p_in[35] = { 
    334,  720, 1020, 1560, 2100, 2300, 2550, 2700, 2880, 2880, 
   2760, 2660, 2550, 2510, 2430, 2270, 2130, 2060, 2000, 1970, 
   1950, 1950, 1960, 1960, 1970, 1980, 1950, 1978, 1830, 1780,
   1735, 1710, 1700, 1700, 1700};
 
   AddDataSet("pi-",he_t,   he_in,  e1, 38);
   AddDataSet("pi+",he_t,   he_in,  e1, 38);
   AddDataSet("pi-",be_m_t, be_m_in, e1, 38);
   AddDataSet("pi+",be_p_t, be_p_in, e1, 38);
   AddDataSet("pi-",c_m_t,  c_m_in,  e2, 39);
   AddDataSet("pi+",c_p_t,  c_p_in,  e2, 39);
   AddDataSet("pi-",n_m_t,  n_m_in,  e2, 39);
   AddDataSet("pi+",n_p_t,  n_p_in,  e2, 39);
   AddDataSet("pi-",o_m_t,  o_m_in,  e3, 31);
   AddDataSet("pi+",o_p_t,  o_p_in,  e3, 31);
   AddDataSet("pi-",na_m_t, na_m_in, e3, 31);
   AddDataSet("pi+",na_p_t, na_p_in, e3, 31);
   AddDataSet("pi-",al_m_t, al_m_in, e3, 31);
   AddDataSet("pi+",al_p_t, al_p_in, e3, 31);
   AddDataSet("pi-",ca_m_t, ca_m_in, e3, 31);
   AddDataSet("pi+",ca_p_t, ca_p_in, e3, 31);
   AddDataSet("pi-",fe_m_t, fe_m_in, e4, 32);
   AddDataSet("pi+",fe_p_t, fe_p_in, e4, 32);
   AddDataSet("pi-",cu_m_t, cu_m_in, e4, 32);
   AddDataSet("pi+",cu_p_t, cu_p_in, e4, 32);
   AddDataSet("pi-",mo_m_t, mo_m_in, e5, 34);
   AddDataSet("pi+",mo_p_t, mo_p_in, e5, 34);
   AddDataSet("pi-",cd_m_t, cd_m_in, e5, 34);
   AddDataSet("pi+",cd_p_t, cd_p_in, e5, 34);
   AddDataSet("pi-",sn_m_t, sn_m_in, e6, 35);
   AddDataSet("pi+",sn_p_t, sn_p_in, e6, 35);
   AddDataSet("pi-",w_m_t,  w_m_in,  e6, 35);
   AddDataSet("pi+",w_p_t,  w_p_in,  e6, 35);
   AddDataSet("pi-",pb_m_t, pb_m_in, e6, 35);
   AddDataSet("pi+",pb_p_t, pb_p_in, e6, 35);
   AddDataSet("pi-",u_m_t,  u_m_in,  e6, 35);
   AddDataSet("pi+",u_p_t,  u_p_in,  e6, 35);
 }

void G4UPiNuclearCrossSection::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 583 of file G4UPiNuclearCrossSection.cc.

 {
   outFile << "G4UPiNuclearCrossSection calculates the total, elastic and\n"
           << "inelastic cross sections for pion scattering from nuclei\n"
           << "heavier than hydrogen.  It is based on the Barashenkov\n"
           << "parameterization and is valid for all incident energies.\n";
 }

void G4UPiNuclearCrossSection::DumpPhysicsTable ( const G4ParticleDefinition & p )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 177 of file G4UPiNuclearCrossSection.cc.

References G4cout, and G4endl.

 {
   if(&p == piPlus) {
     G4cout << "### G4UPiNuclearCrossSection Elastic data for pi+" << G4endl;
     G4cout << *piPlusElastic << G4endl;
     G4cout << "### G4UPiNuclearCrossSection Inelastic data for pi+" << G4endl;
     G4cout << *piPlusInelastic << G4endl;
   } else if(&p == piMinus) {
     G4cout << "### G4UPiNuclearCrossSection Elastic data for pi-" << G4endl;
     G4cout << *piMinusElastic << G4endl;
     G4cout << "### G4UPiNuclearCrossSection Inelastic data for pi-" << G4endl;
     G4cout << *piMinusInelastic << G4endl;
   }
 }

G4double G4UPiNuclearCrossSection::GetElasticCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		G4int	A
	)

Definition at line 81 of file G4UPiNuclearCrossSection.cc.

References G4DynamicParticle::GetDefinition(), and G4DynamicParticle::GetKineticEnergy().

Referenced by G4BGGPionElasticXS::BuildPhysicsTable(), and G4BGGPionElasticXS::GetElementCrossSection().

 {
   G4double cross = 0.0;
   G4PhysicsTable* table = 0;
   const G4ParticleDefinition* part = dp->GetDefinition();
   if(part == piPlus) { table = piPlusElastic; }
   else if(part == piMinus) { table = piMinusElastic; }
   if(table) {
     cross = Interpolate(Z, A, dp->GetKineticEnergy(),table);
   }
   return cross;
 }

G4double G4UPiNuclearCrossSection::GetInelasticCrossSection	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		G4int	A
	)

Definition at line 96 of file G4UPiNuclearCrossSection.cc.

References G4DynamicParticle::GetDefinition(), and G4DynamicParticle::GetKineticEnergy().

Referenced by G4BGGPionInelasticXS::BuildPhysicsTable(), and G4BGGPionInelasticXS::GetElementCrossSection().

 {
   G4double cross = 0.0;
   G4double fact  = 1.0;
   G4double ekin  = dp->GetKineticEnergy();
   G4PhysicsTable* table = 0;
   const G4ParticleDefinition* part = dp->GetDefinition();
 
   // Coulomb barrier
   if(part == piPlus) {
     if(ekin > elowest) {
       table = piPlusInelastic;
       if(ekin < elow) {
         fact = std::sqrt((ekin-elowest)/(elow-elowest));
         ekin = elow;
       }
     }
   } else if(part == piMinus) {
     table = piMinusInelastic;
     if(ekin < elow) { ekin = elow; }
   }
   if(table) {
     cross = fact*Interpolate(Z, A, ekin, table);
   }
   return cross;
 }

G4bool G4UPiNuclearCrossSection::IsElementApplicable	(	const G4DynamicParticle *	aParticle,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 74 of file G4UPiNuclearCrossSection.cc.

 {
   return (1 < Z);
 }

The documentation for this class was generated from the following files:

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation