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().

                                                               : theR(aR)
{
  G4ResonanceID theID;
  G4int pdg=theR->GetPDGEncoding();
  data=0;
       if(theID.IsDelta1232(pdg)) data = delta;
  else if(theID.IsDelta1600(pdg)) data = delta1600;
  else if(theID.IsDelta1620(pdg)) data = delta1620;
  else if(theID.IsDelta1700(pdg)) data = delta1700;
  else if(theID.IsDelta1900(pdg)) data = delta1900;
  else if(theID.IsDelta1905(pdg)) data = delta1905;
  else if(theID.IsDelta1910(pdg)) data = delta1910;
  else if(theID.IsDelta1920(pdg)) data = delta1920;
  else if(theID.IsDelta1930(pdg)) data = delta1930;
  else if(theID.IsDelta1950(pdg)) data = delta1950;
  
  else if(theID.IsN1440(pdg)) data = N1440;
  else if(theID.IsN1520(pdg)) data = N1520;
  else if(theID.IsN1535(pdg)) data = N1535;
  else if(theID.IsN1650(pdg)) data = N1650;
  else if(theID.IsN1675(pdg)) data = N1675;
  else if(theID.IsN1680(pdg)) data = N1680;
  else if(theID.IsN1700(pdg)) data = N1700;
  else if(theID.IsN1710(pdg)) data = N1710;
  else if(theID.IsN1720(pdg)) data = N1720;
  else if(theID.IsN1900(pdg)) data = N1900;
  else if(theID.IsN1990(pdg)) data = N1990;
  else if(theID.IsN2090(pdg)) data = N2090;
  else if(theID.IsN2190(pdg)) data = N2190;
  else if(theID.IsN2220(pdg)) data = N2220;
  else if(theID.IsN2250(pdg)) data = N2250;
}

Member Function Documentation

G4double G4DetailedBalancePhaseSpaceIntegral::GetPhaseSpaceIntegral

(

G4double

sqrts

)

Definition at line 146 of file G4DetailedBalancePhaseSpaceIntegral.cc.

References python.hepunit::GeV.

Referenced by G4VXResonance::DetailedBalance().

{
  G4double result = 0;
  G4int it=0;
  for(G4int ie=0; ie<120; ie++)
  {
    if(sqrts[ie]*GeV>sqs) break;
    it = ie;
  }
  G4double x1=sqrts[it]*GeV;
  G4double x2=sqrts[it+1]*GeV;
  G4double y1=data[it];
  G4double y2=data[it+1];
  result = y1+(sqs-x1)*(y2-y1)/(x2-x1);
  return result;
}

---

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

• G4DetailedBalancePhaseSpaceIntegral.hh
• G4DetailedBalancePhaseSpaceIntegral.cc