G4DetailedBalancePhaseSpaceIntegral Class Reference

#include <G4DetailedBalancePhaseSpaceIntegral.hh>


Public Member Functions

 G4DetailedBalancePhaseSpaceIntegral (G4ParticleDefinition *aR)
G4double GetPhaseSpaceIntegral (G4double sqrts)


Detailed Description

Definition at line 32 of file G4DetailedBalancePhaseSpaceIntegral.hh.


Constructor & Destructor Documentation

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 }


Member Function Documentation

G4double G4DetailedBalancePhaseSpaceIntegral::GetPhaseSpaceIntegral ( G4double  sqrts  ) 

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 }


The documentation for this class was generated from the following files:
Generated on Mon May 27 17:51:45 2013 for Geant4 by  doxygen 1.4.7