#include <G4DetailedBalancePhaseSpaceIntegral.hh>
Public Member Functions | |
G4DetailedBalancePhaseSpaceIntegral (G4ParticleDefinition *aR) | |
G4double | GetPhaseSpaceIntegral (G4double sqrts) |
Definition at line 32 of file G4DetailedBalancePhaseSpaceIntegral.hh.
G4DetailedBalancePhaseSpaceIntegral::G4DetailedBalancePhaseSpaceIntegral | ( | G4ParticleDefinition * | aR | ) |
Definition at line 112 of file G4DetailedBalancePhaseSpaceIntegral.cc.
References G4ParticleDefinition::GetPDGEncoding(), G4ResonanceID::IsDelta1232(), G4ResonanceID::IsDelta1600(), G4ResonanceID::IsDelta1620(), G4ResonanceID::IsDelta1700(), G4ResonanceID::IsDelta1900(), G4ResonanceID::IsDelta1905(), G4ResonanceID::IsDelta1910(), G4ResonanceID::IsDelta1920(), G4ResonanceID::IsDelta1930(), G4ResonanceID::IsDelta1950(), G4ResonanceID::IsN1440(), G4ResonanceID::IsN1520(), G4ResonanceID::IsN1535(), G4ResonanceID::IsN1650(), G4ResonanceID::IsN1675(), G4ResonanceID::IsN1680(), G4ResonanceID::IsN1700(), G4ResonanceID::IsN1710(), G4ResonanceID::IsN1720(), G4ResonanceID::IsN1900(), G4ResonanceID::IsN1990(), G4ResonanceID::IsN2090(), G4ResonanceID::IsN2190(), G4ResonanceID::IsN2220(), and G4ResonanceID::IsN2250().
00112 : theR(aR) 00113 { 00114 G4ResonanceID theID; 00115 G4int pdg=theR->GetPDGEncoding(); 00116 data=0; 00117 if(theID.IsDelta1232(pdg)) data = delta; 00118 else if(theID.IsDelta1600(pdg)) data = delta1600; 00119 else if(theID.IsDelta1620(pdg)) data = delta1620; 00120 else if(theID.IsDelta1700(pdg)) data = delta1700; 00121 else if(theID.IsDelta1900(pdg)) data = delta1900; 00122 else if(theID.IsDelta1905(pdg)) data = delta1905; 00123 else if(theID.IsDelta1910(pdg)) data = delta1910; 00124 else if(theID.IsDelta1920(pdg)) data = delta1920; 00125 else if(theID.IsDelta1930(pdg)) data = delta1930; 00126 else if(theID.IsDelta1950(pdg)) data = delta1950; 00127 00128 else if(theID.IsN1440(pdg)) data = N1440; 00129 else if(theID.IsN1520(pdg)) data = N1520; 00130 else if(theID.IsN1535(pdg)) data = N1535; 00131 else if(theID.IsN1650(pdg)) data = N1650; 00132 else if(theID.IsN1675(pdg)) data = N1675; 00133 else if(theID.IsN1680(pdg)) data = N1680; 00134 else if(theID.IsN1700(pdg)) data = N1700; 00135 else if(theID.IsN1710(pdg)) data = N1710; 00136 else if(theID.IsN1720(pdg)) data = N1720; 00137 else if(theID.IsN1900(pdg)) data = N1900; 00138 else if(theID.IsN1990(pdg)) data = N1990; 00139 else if(theID.IsN2090(pdg)) data = N2090; 00140 else if(theID.IsN2190(pdg)) data = N2190; 00141 else if(theID.IsN2220(pdg)) data = N2220; 00142 else if(theID.IsN2250(pdg)) data = N2250; 00143 }
Definition at line 146 of file G4DetailedBalancePhaseSpaceIntegral.cc.
00147 { 00148 G4double result = 0; 00149 G4int it=0; 00150 for(G4int ie=0; ie<120; ie++) 00151 { 00152 if(sqrts[ie]*GeV>sqs) break; 00153 it = ie; 00154 } 00155 G4double x1=sqrts[it]*GeV; 00156 G4double x2=sqrts[it+1]*GeV; 00157 G4double y1=data[it]; 00158 G4double y2=data[it+1]; 00159 result = y1+(sqs-x1)*(y2-y1)/(x2-x1); 00160 return result; 00161 }