G4PAIySection Class Reference

#include <G4PAIySection.hh>

Public Member Functions
	G4PAIySection ()
	~G4PAIySection ()
void	Initialize (const G4Material *material, G4double maxEnergyTransfer, G4double betaGammaSq)
void	ComputeLowEnergyCof (const G4Material *material)
void	InitPAI ()
void	NormShift (G4double betaGammaSq)
void	SplainPAI (G4double betaGammaSq)
G4double	RutherfordIntegral (G4int intervalNumber, G4double limitLow, G4double limitHigh)
G4double	ImPartDielectricConst (G4int intervalNumber, G4double energy)
G4double	RePartDielectricConst (G4double energy)
G4double	DifPAIySection (G4int intervalNumber, G4double betaGammaSq)
G4double	PAIdNdxCerenkov (G4int intervalNumber, G4double betaGammaSq)
G4double	PAIdNdxPlasmon (G4int intervalNumber, G4double betaGammaSq)
void	IntegralPAIySection ()
void	IntegralCerenkov ()
void	IntegralPlasmon ()
G4double	SumOverInterval (G4int intervalNumber)
G4double	SumOverIntervaldEdx (G4int intervalNumber)
G4double	SumOverInterCerenkov (G4int intervalNumber)
G4double	SumOverInterPlasmon (G4int intervalNumber)
G4double	SumOverBorder (G4int intervalNumber, G4double energy)
G4double	SumOverBorderdEdx (G4int intervalNumber, G4double energy)
G4double	SumOverBordCerenkov (G4int intervalNumber, G4double energy)
G4double	SumOverBordPlasmon (G4int intervalNumber, G4double energy)
G4double	GetStepEnergyLoss (G4double step)
G4double	GetStepCerenkovLoss (G4double step)
G4double	GetStepPlasmonLoss (G4double step)
G4int	GetNumberOfGammas () const
G4int	GetSplineSize () const
G4int	GetIntervalNumber () const
G4double	GetEnergyInterval (G4int i)
G4double	GetDifPAIySection (G4int i)
G4double	GetPAIdNdxCrenkov (G4int i)
G4double	GetPAIdNdxPlasmon (G4int i)
G4double	GetMeanEnergyLoss () const
G4double	GetMeanCerenkovLoss () const
G4double	GetMeanPlasmonLoss () const
G4double	GetNormalizationCof () const
G4double	GetPAItable (G4int i, G4int j) const
G4double	GetLorentzFactor (G4int i) const
G4double	GetSplineEnergy (G4int i) const
G4double	GetIntegralPAIySection (G4int i) const
G4double	GetIntegralPAIdEdx (G4int i) const
G4double	GetIntegralCerenkov (G4int i) const
G4double	GetIntegralPlasmon (G4int i) const
void	SetVerbose (G4int v)

Detailed Description

Definition at line 52 of file G4PAIySection.hh.

Constructor & Destructor Documentation

G4PAIySection()

G4PAIySection::G4PAIySection

(

)

Definition at line 73 of file G4PAIySection.cc.

{
  fSandia = 0;
  fDensity = fElectronDensity = fNormalizationCof = fLowEnergyCof = 0.0;
  fIntervalNumber = fSplineNumber = 0;
  fVerbose = 0;
  for(G4int i=0; i<500; ++i) {
    fSplineEnergy[i] = 0.0;
    fRePartDielectricConst[i] = 0.0;
    fImPartDielectricConst[i] = 0.0;
    fIntegralTerm[i] = 0.0;
    fDifPAIySection[i] = 0.0;
    fdNdxCerenkov[i] = 0.0;
    fdNdxPlasmon[i] = 0.0;
    fIntegralPAIySection[i] = 0.0;
    fIntegralPAIdEdx[i] = 0.0;
    fIntegralCerenkov[i] = 0.0;
    fIntegralPlasmon[i] = 0.0;
    for(G4int j=0; j<112; ++j) { fPAItable[i][j] = 0.0; }
  }
}

~G4PAIySection()

G4PAIySection::~G4PAIySection

(

)

Definition at line 99 of file G4PAIySection.cc.

{}

Member Function Documentation

ComputeLowEnergyCof()

void G4PAIySection::ComputeLowEnergyCof

(

const G4Material *

material

)

Definition at line 215 of file G4PAIySection.cc.

{    
  G4int i, numberOfElements = material->GetNumberOfElements();
  G4double sumZ = 0., sumCof = 0.; 
 
  const G4double p0 =  1.20923e+00; 
  const G4double p1 =  3.53256e-01; 
  const G4double p2 = -1.45052e-03; 
  
  G4double* thisMaterialZ   = new G4double[numberOfElements];
  G4double* thisMaterialCof = new G4double[numberOfElements];
   
  for( i = 0; i < numberOfElements; i++ )
  {
    thisMaterialZ[i] = material->GetElement(i)->GetZ();
    sumZ += thisMaterialZ[i];
    thisMaterialCof[i] = p0+p1*thisMaterialZ[i]+p2*thisMaterialZ[i]*thisMaterialZ[i];   
  }
  for( i = 0; i < numberOfElements; i++ )
  {
    sumCof += thisMaterialCof[i]*thisMaterialZ[i]/sumZ;
  }
  fLowEnergyCof = sumCof;
  delete [] thisMaterialZ;
  delete [] thisMaterialCof;
  // G4cout<<"fLowEnergyCof = "<<fLowEnergyCof<<G4endl;
}

Referenced by Initialize().

DifPAIySection()

G4double G4PAIySection::DifPAIySection	(	G4int	intervalNumber,
		G4double	betaGammaSq
	)

Definition at line 558 of file G4PAIySection.cc.

{        
  G4double beta, be2,cof,x1,x2,x3,x4,x5,x6,x7,x8,result;
   //G4double beta, be4;
   //G4double be4;
   G4double betaBohr = fine_structure_const;
   // G4double betaBohr2 = fine_structure_const*fine_structure_const;
   // G4double betaBohr4 = betaBohr2*betaBohr2*4.0;
   be2 = betaGammaSq/(1 + betaGammaSq);
   //be4 = be2*be2;
   beta = sqrt(be2);
   cof = 1;
   x1 = log(2*electron_mass_c2/fSplineEnergy[i]);
 
   if( betaGammaSq < 0.01 ) x2 = log(be2);
   else
   {
     x2 = -log( (1/betaGammaSq - fRePartDielectricConst[i])*
            (1/betaGammaSq - fRePartDielectricConst[i]) + 
            fImPartDielectricConst[i]*fImPartDielectricConst[i] )/2;
   }
   if( fImPartDielectricConst[i] == 0.0 ||betaGammaSq < 0.01 )
   {
     x6=0;
   }
   else
   {
     x3 = -fRePartDielectricConst[i] + 1/betaGammaSq;
     x5 = -1 - fRePartDielectricConst[i] +
          be2*((1 +fRePartDielectricConst[i])*(1 + fRePartDielectricConst[i]) +
      fImPartDielectricConst[i]*fImPartDielectricConst[i]);
 
     x7 = atan2(fImPartDielectricConst[i],x3);
     x6 = x5 * x7;
   }
    // if(fImPartDielectricConst[i] == 0) x6 = 0;
   
   x4 = ((x1 + x2)*fImPartDielectricConst[i] + x6)/hbarc;
   //   if( x4 < 0.0 ) x4 = 0.0;
   x8 = (1 + fRePartDielectricConst[i])*(1 + fRePartDielectricConst[i]) + 
        fImPartDielectricConst[i]*fImPartDielectricConst[i];
 
   result = (x4 + cof*fIntegralTerm[i]/fSplineEnergy[i]/fSplineEnergy[i]);
   if(result < 1.0e-8) result = 1.0e-8;
   result *= fine_structure_const/be2/pi;
   // low energy correction
 
   G4double lowCof = fLowEnergyCof; // 6.0 ; // Ar ~ 4.; -> fLowCof as f(Z1,Z2)? 
 
   result *= (1 - exp(-beta/betaBohr/lowCof));
 
   //   result *= (1-exp(-beta/betaBohr))*(1-exp(-beta/betaBohr));
   //  result *= (1-exp(-be2/betaBohr2));
   // result *= (1-exp(-be4/betaBohr4));
   //   if(fDensity >= 0.1)
   if(x8 > 0.)
   { 
      result /= x8;
   }
   return result;
 
} // end of DifPAIySection 

Referenced by Initialize(), InitPAI(), NormShift(), and SplainPAI().

GetDifPAIySection()

G4double G4PAIySection::GetDifPAIySection ( G4int  i )

inline

Definition at line 123 of file G4PAIySection.hh.

{ return fDifPAIySection[i]; } 

GetEnergyInterval()

G4double G4PAIySection::GetEnergyInterval ( G4int  i )

inline

Definition at line 121 of file G4PAIySection.hh.

{ return fEnergyInterval[i]; } 

GetIntegralCerenkov()

G4double G4PAIySection::GetIntegralCerenkov ( G4int  i ) const

inline

Definition at line 230 of file G4PAIySection.hh.

{
  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralCerenkov"); }
  return fIntegralCerenkov[i];
}

GetIntegralPAIdEdx()

G4double G4PAIySection::GetIntegralPAIdEdx ( G4int  i ) const

inline

Definition at line 224 of file G4PAIySection.hh.

{
  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralPAIdEdx"); }
  return fIntegralPAIdEdx[i];
}

Referenced by G4PAIModel::BuildPAIonisationTable().

GetIntegralPAIySection()

G4double G4PAIySection::GetIntegralPAIySection ( G4int  i ) const

inline

Definition at line 218 of file G4PAIySection.hh.

{
  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralPAIySection"); }
  return fIntegralPAIySection[i];
}

Referenced by G4PAIModel::BuildPAIonisationTable().

GetIntegralPlasmon()

G4double G4PAIySection::GetIntegralPlasmon ( G4int  i ) const

inline

Definition at line 236 of file G4PAIySection.hh.

{
  if(i < 1 || i > fSplineNumber) { CallError(i, "GetIntegralPlasmon"); }
  return fIntegralPlasmon[i];
}

GetIntervalNumber()

G4int G4PAIySection::GetIntervalNumber ( ) const

inline

Definition at line 119 of file G4PAIySection.hh.

{ return fIntervalNumber; }

GetLorentzFactor()

G4double G4PAIySection::GetLorentzFactor ( G4int  i ) const

inline

Definition at line 207 of file G4PAIySection.hh.

{
   return fLorentzFactor[j];
}

GetMeanCerenkovLoss()

G4double G4PAIySection::GetMeanCerenkovLoss ( ) const

inline

Definition at line 128 of file G4PAIySection.hh.

{return fIntegralCerenkov[0]; }

GetMeanEnergyLoss()

G4double G4PAIySection::GetMeanEnergyLoss ( ) const

inline

Definition at line 127 of file G4PAIySection.hh.

{return fIntegralPAIySection[0]; }

Referenced by G4PAIModel::BuildPAIonisationTable().

GetMeanPlasmonLoss()

G4double G4PAIySection::GetMeanPlasmonLoss ( ) const

inline

Definition at line 129 of file G4PAIySection.hh.

{return fIntegralPlasmon[0]; }

GetNormalizationCof()

G4double G4PAIySection::GetNormalizationCof ( ) const

inline

Definition at line 131 of file G4PAIySection.hh.

{ return fNormalizationCof; }

GetNumberOfGammas()

G4int G4PAIySection::GetNumberOfGammas ( ) const

inline

Definition at line 115 of file G4PAIySection.hh.

{ return fNumberOfGammas; }

GetPAIdNdxCrenkov()

G4double G4PAIySection::GetPAIdNdxCrenkov ( G4int  i )

inline

Definition at line 124 of file G4PAIySection.hh.

{ return fdNdxCerenkov[i]; } 

GetPAIdNdxPlasmon()

G4double G4PAIySection::GetPAIdNdxPlasmon ( G4int  i )

inline

Definition at line 125 of file G4PAIySection.hh.

{ return fdNdxPlasmon[i]; } 

GetPAItable()

G4double G4PAIySection::GetPAItable ( G4int  i, G4int  j  ) const

inline

Definition at line 202 of file G4PAIySection.hh.

{
   return fPAItable[i][j];
}

GetSplineEnergy()

G4double G4PAIySection::GetSplineEnergy ( G4int  i ) const

inline

Definition at line 212 of file G4PAIySection.hh.

{
  if(i < 1 || i > fSplineNumber) { CallError(i, "GetSplineEnergy"); }
  return fSplineEnergy[i];
}

Referenced by G4PAIModel::BuildPAIonisationTable().

GetSplineSize()

G4int G4PAIySection::GetSplineSize ( ) const

inline

Definition at line 117 of file G4PAIySection.hh.

{ return fSplineNumber; }

Referenced by G4PAIModel::BuildPAIonisationTable().

GetStepCerenkovLoss()

G4double G4PAIySection::GetStepCerenkovLoss

(

G4double

step

)

Definition at line 1238 of file G4PAIySection.cc.

{  
  G4int iTransfer ;
  G4long numOfCollisions;
  G4double loss = 0.0;
  G4double meanNumber, position;
 
  // G4cout<<" G4PAIySection::GetStepCreLosnkovs "<<G4endl;
 
 
 
  meanNumber = fIntegralCerenkov[1]*step;
  numOfCollisions = G4Poisson(meanNumber);
 
  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;
 
  while(numOfCollisions)
  {
    position = fIntegralCerenkov[1]*G4UniformRand();
 
    for( iTransfer=1; iTransfer<=fSplineNumber; iTransfer++ )
    {
        if( position >= fIntegralCerenkov[iTransfer] ) break;
    }
    loss += fSplineEnergy[iTransfer] ;
    numOfCollisions--;
  }
  // G4cout<<"PAI Cerenkov loss = "<<loss/keV<<" keV"<<G4endl; 
 
  return loss;
}

GetStepEnergyLoss()

G4double G4PAIySection::GetStepEnergyLoss

(

G4double

step

)

Definition at line 1202 of file G4PAIySection.cc.

{  
  G4int iTransfer ;
  G4long numOfCollisions;
  G4double loss = 0.0;
  G4double meanNumber, position;
 
  // G4cout<<" G4PAIySection::GetStepEnergyLoss "<<G4endl;
 
 
 
  meanNumber = fIntegralPAIySection[1]*step;
  numOfCollisions = G4Poisson(meanNumber);
 
  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;
 
  while(numOfCollisions)
  {
    position = fIntegralPAIySection[1]*G4UniformRand();
 
    for( iTransfer=1; iTransfer<=fSplineNumber; iTransfer++ )
    {
        if( position >= fIntegralPAIySection[iTransfer] ) break;
    }
    loss += fSplineEnergy[iTransfer] ;
    numOfCollisions--;
  }
  // G4cout<<"PAI energy loss = "<<loss/keV<<" keV"<<G4endl; 
 
  return loss;
}

GetStepPlasmonLoss()

G4double G4PAIySection::GetStepPlasmonLoss

(

G4double

step

)

Definition at line 1274 of file G4PAIySection.cc.

{  
  G4int iTransfer ;
  G4long numOfCollisions;
  G4double loss = 0.0;
  G4double meanNumber, position;
 
  // G4cout<<" G4PAIySection::GetStepCreLosnkovs "<<G4endl;
 
 
 
  meanNumber = fIntegralPlasmon[1]*step;
  numOfCollisions = G4Poisson(meanNumber);
 
  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;
 
  while(numOfCollisions)
  {
    position = fIntegralPlasmon[1]*G4UniformRand();
 
    for( iTransfer=1; iTransfer<=fSplineNumber; iTransfer++ )
    {
        if( position >= fIntegralPlasmon[iTransfer] ) break;
    }
    loss += fSplineEnergy[iTransfer] ;
    numOfCollisions--;
  }
  // G4cout<<"PAI Plasmon loss = "<<loss/keV<<" keV"<<G4endl; 
 
  return loss;
}

ImPartDielectricConst()

G4double G4PAIySection::ImPartDielectricConst	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 481 of file G4PAIySection.cc.

{
   G4double energy2,energy3,energy4,result;
 
   energy2 = energy1*energy1;
   energy3 = energy2*energy1;
   energy4 = energy3*energy1;
   
   result = fA1[k]/energy1+fA2[k]/energy2+fA3[k]/energy3+fA4[k]/energy4;  
   result *=hbarc/energy1;
   
   return result;
 
}  // end of ImPartDielectricConst 

Referenced by NormShift(), and SplainPAI().

Initialize()

void G4PAIySection::Initialize	(	const G4Material *	material,
		G4double	maxEnergyTransfer,
		G4double	betaGammaSq
	)

Definition at line 106 of file G4PAIySection.cc.

{
  G4int i, j;
  G4double energy;   
  // fVerbose = 1;   
  fDensity         = material->GetDensity();
  fElectronDensity = material->GetElectronDensity();
  //G4int numberOfElements = material->GetNumberOfElements();
 
  fSandia = material->GetSandiaTable();
 
  fIntervalNumber = fSandia->GetMaxInterval();
 
  fIntervalNumber--;
 
  for( i = 1; i <= fIntervalNumber; i++ )
  {
    energy = fSandia->GetSandiaMatTablePAI(i-1,0); //vmg 20.11.10
 
    if( energy >= maxEnergyTransfer || i > fIntervalNumber )
    {
      fEnergyInterval[i] = maxEnergyTransfer;
      fIntervalNumber = i;
      break;
    }
    fEnergyInterval[i] = energy;
    fA1[i]             = fSandia->GetSandiaMatTablePAI(i-1,1);
    fA2[i]             = fSandia->GetSandiaMatTablePAI(i-1,2);
    fA3[i]             = fSandia->GetSandiaMatTablePAI(i-1,3);
    fA4[i]             = fSandia->GetSandiaMatTablePAI(i-1,4);
 
    if( fVerbose > 0 && std::fabs( betaGammaSq - 8. ) < 0.4 )
    {
    G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<" keV \t"<<fA1[i]<<"\t"<<fA2[i] <<"\t"<<fA3[i] <<"\t"<<fA4[i]<<G4endl;
    }
  } 
 
  
  if( fEnergyInterval[fIntervalNumber] != maxEnergyTransfer )
  {
    fIntervalNumber++;
    fEnergyInterval[fIntervalNumber] = maxEnergyTransfer;
    fA1[fIntervalNumber] = fA1[fIntervalNumber-1];
    fA2[fIntervalNumber] = fA2[fIntervalNumber-1];
    fA3[fIntervalNumber] = fA3[fIntervalNumber-1];
    fA4[fIntervalNumber] = fA4[fIntervalNumber-1];
  }
 
   // Now checking, if two borders are too close together
  for( i = 1; i < fIntervalNumber; i++ )
  {
    // G4cout<<fEnergyInterval[i]<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
    //   <<fA3[i]<<"\t"<<fA4[i]<<G4endl;
    if(fEnergyInterval[i+1]-fEnergyInterval[i] <
           1.5*fDelta*(fEnergyInterval[i+1]+fEnergyInterval[i]))
    {
      for( j = i; j < fIntervalNumber; j++ )
      {
    fEnergyInterval[j] = fEnergyInterval[j+1];
    fA1[j] = fA1[j+1];
    fA2[j] = fA2[j+1];
    fA3[j] = fA3[j+1];
    fA4[j] = fA4[j+1];
      }
      fIntervalNumber--;
      i--;
    }
  }
  if( fVerbose > 0 && std::fabs( betaGammaSq - 8. ) < 0.4 )
  {
    G4cout<<"Sandia cofs in G4PAIySection::Initialize(), mat = "<<material->GetName()<<G4endl;
    G4cout<<"for bg2 = "<<betaGammaSq<<", Tmax = "<< maxEnergyTransfer/keV<<" keV"<<G4endl;
    G4cout<<"energy \t"<<"a1 \t"<<"a2 \t"<<"a3 \t"<<"a4 \t"<<G4endl;
 
      for( j = 1; j < fIntervalNumber; j++ )
      {
    G4cout<<j<<"\t"<<fEnergyInterval[j]/keV<<" keV \t"<<fA1[j]<<"\t"<<fA2[j] <<"\t"<<fA3[j] <<"\t"<<fA4[j]<<G4endl;
      }
  }
 
      // Preparation of fSplineEnergy array corresponding to min ionisation, G~4
      
  G4double   betaGammaSqRef = 
     fLorentzFactor[fRefGammaNumber]*fLorentzFactor[fRefGammaNumber] - 1;
      
  ComputeLowEnergyCof(material);            
  NormShift(betaGammaSqRef);             
  SplainPAI(betaGammaSqRef);
      
   // Preparation of integral PAI cross section for input betaGammaSq
   
  for( i = 1; i <= fSplineNumber; i++ )
  {
    fDifPAIySection[i] = DifPAIySection(i,betaGammaSq);
    fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
    fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
  }
  IntegralPAIySection();
  IntegralCerenkov();
  IntegralPlasmon();
}

Referenced by G4PAIModel::BuildPAIonisationTable().

InitPAI()

void G4PAIySection::InitPAI

(

)

Definition at line 251 of file G4PAIySection.cc.

{    
   G4int i;
   G4double betaGammaSq = fLorentzFactor[fRefGammaNumber]*
                          fLorentzFactor[fRefGammaNumber] - 1;
 
   // Preparation of integral PAI cross section for reference gamma
   
   NormShift(betaGammaSq);             
   SplainPAI(betaGammaSq);
 
   IntegralPAIySection();
   IntegralCerenkov();
   IntegralPlasmon();
 
   for(i = 0; i<=fSplineNumber; i++)
   {
      fPAItable[i][fRefGammaNumber] = fIntegralPAIySection[i];
      if(i != 0) 
      {
     fPAItable[i][0] = fSplineEnergy[i];
      }
   }
   fPAItable[0][0] = fSplineNumber;
   
   for(G4int j = 1; j < 112; j++)       // for other gammas
   {
      if( j == fRefGammaNumber ) continue;
      
      betaGammaSq = fLorentzFactor[j]*fLorentzFactor[j] - 1;
      
      for(i = 1; i <= fSplineNumber; i++)
      {
         fDifPAIySection[i] = DifPAIySection(i,betaGammaSq);
         fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
         fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
      }
      IntegralPAIySection();
      IntegralCerenkov();
      IntegralPlasmon();
      
      for(i = 0; i <= fSplineNumber; i++)
      {
         fPAItable[i][j] = fIntegralPAIySection[i];
      }
   } 
 
}  

IntegralCerenkov()

void G4PAIySection::IntegralCerenkov

(

)

Definition at line 762 of file G4PAIySection.cc.

{
  G4int i, k;
   fIntegralCerenkov[fSplineNumber] = 0;
   fIntegralCerenkov[0] = 0;
   k = fIntervalNumber -1;
 
   for( i = fSplineNumber-1; i >= 1; i-- )
   {
      if(fSplineEnergy[i] >= fEnergyInterval[k])
      {
        fIntegralCerenkov[i] = fIntegralCerenkov[i+1] + SumOverInterCerenkov(i);
    // G4cout<<"int: i = "<<i<<"; sumC = "<<fIntegralCerenkov[i]<<G4endl;
      }
      else
      {
        fIntegralCerenkov[i] = fIntegralCerenkov[i+1] + 
                               SumOverBordCerenkov(i+1,fEnergyInterval[k]);
    k--;
    // G4cout<<"bord: i = "<<i<<"; sumC = "<<fIntegralCerenkov[i]<<G4endl;
      }
   }
 
}   // end of IntegralCerenkov 

Referenced by Initialize(), and InitPAI().

IntegralPAIySection()

void G4PAIySection::IntegralPAIySection

(

)

Definition at line 731 of file G4PAIySection.cc.

{
  fIntegralPAIySection[fSplineNumber] = 0;
  fIntegralPAIdEdx[fSplineNumber]     = 0;
  fIntegralPAIySection[0]             = 0;
  G4int k = fIntervalNumber -1;
 
  for(G4int i = fSplineNumber-1; i >= 1; i--)
  {
    if(fSplineEnergy[i] >= fEnergyInterval[k])
    {
      fIntegralPAIySection[i] = fIntegralPAIySection[i+1] + SumOverInterval(i);
      fIntegralPAIdEdx[i] = fIntegralPAIdEdx[i+1] + SumOverIntervaldEdx(i);
    }
    else
    {
      fIntegralPAIySection[i] = fIntegralPAIySection[i+1] + 
                               SumOverBorder(i+1,fEnergyInterval[k]);
      fIntegralPAIdEdx[i] = fIntegralPAIdEdx[i+1] + 
                               SumOverBorderdEdx(i+1,fEnergyInterval[k]);
      k--;
    }
  }
}   // end of IntegralPAIySection 

Referenced by Initialize(), and InitPAI().

IntegralPlasmon()

void G4PAIySection::IntegralPlasmon

(

)

Definition at line 793 of file G4PAIySection.cc.

{
   fIntegralPlasmon[fSplineNumber] = 0;
   fIntegralPlasmon[0] = 0;
   G4int k = fIntervalNumber -1;
   for(G4int i=fSplineNumber-1;i>=1;i--)
   {
      if(fSplineEnergy[i] >= fEnergyInterval[k])
      {
        fIntegralPlasmon[i] = fIntegralPlasmon[i+1] + SumOverInterPlasmon(i);
      }
      else
      {
        fIntegralPlasmon[i] = fIntegralPlasmon[i+1] + 
                               SumOverBordPlasmon(i+1,fEnergyInterval[k]);
    k--;
      }
   }
 
}   // end of IntegralPlasmon

Referenced by Initialize(), and InitPAI().

NormShift()

void G4PAIySection::NormShift

(

G4double

betaGammaSq

)

Definition at line 305 of file G4PAIySection.cc.

{
  G4int i, j;
 
  for( i = 1; i <= fIntervalNumber-1; i++ )
  {
    for( j = 1; j <= 2; j++ )
    {
      fSplineNumber = (i-1)*2 + j;
 
      if( j == 1 ) fSplineEnergy[fSplineNumber] = fEnergyInterval[i  ]*(1+fDelta);
      else         fSplineEnergy[fSplineNumber] = fEnergyInterval[i+1]*(1-fDelta); 
      //    G4cout<<"cn = "<<fSplineNumber<<"; "<<"energy = "
      //  <<fSplineEnergy[fSplineNumber]<<G4endl;
    }
  }
  fIntegralTerm[1]=RutherfordIntegral(1,fEnergyInterval[1],fSplineEnergy[1]);
 
  j = 1;
 
  for(i=2;i<=fSplineNumber;i++)
  {
    if(fSplineEnergy[i]<fEnergyInterval[j+1])
    {
         fIntegralTerm[i] = fIntegralTerm[i-1] + 
                        RutherfordIntegral(j,fSplineEnergy[i-1],
                                                 fSplineEnergy[i]   );
    }
    else
    {
       G4double x = RutherfordIntegral(j,fSplineEnergy[i-1],
                                           fEnergyInterval[j+1]   );
         j++;
         fIntegralTerm[i] = fIntegralTerm[i-1] + x + 
                        RutherfordIntegral(j,fEnergyInterval[j],
                                                 fSplineEnergy[i]    );
    }
    // G4cout<<i<<"\t"<<fSplineEnergy[i]<<"\t"<<fIntegralTerm[i]<<"\n"<<G4endl;
  } 
  fNormalizationCof = 2*pi*pi*hbarc*hbarc*fine_structure_const/electron_mass_c2;
  fNormalizationCof *= fElectronDensity/fIntegralTerm[fSplineNumber];
 
  // G4cout<<"fNormalizationCof = "<<fNormalizationCof<<G4endl;
 
      // Calculation of PAI differrential cross-section (1/(keV*cm))
      // in the energy points near borders of energy intervals
 
   for(G4int k=1;k<=fIntervalNumber-1;k++)
   {
      for(j=1;j<=2;j++)
      {
         i = (k-1)*2 + j;
         fImPartDielectricConst[i] = fNormalizationCof*
                                 ImPartDielectricConst(k,fSplineEnergy[i]);
         fRePartDielectricConst[i] = fNormalizationCof*
                                 RePartDielectricConst(fSplineEnergy[i]);
         fIntegralTerm[i] *= fNormalizationCof;
 
         fDifPAIySection[i] = DifPAIySection(i,betaGammaSq);
         fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
         fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
      }
   }
 
}  // end of NormShift 

Referenced by Initialize(), and InitPAI().

PAIdNdxCerenkov()

G4double G4PAIySection::PAIdNdxCerenkov	(	G4int	intervalNumber,
		G4double	betaGammaSq
	)

Definition at line 626 of file G4PAIySection.cc.

{        
   G4double logarithm, x3, x5, argument, modul2, dNdxC; 
   G4double be2, be4, betaBohr2,betaBohr4,cofBetaBohr;
 
   //G4double cof         = 1.0;
   cofBetaBohr = 4.0;
   betaBohr2   = fine_structure_const*fine_structure_const;
   betaBohr4   = betaBohr2*betaBohr2*cofBetaBohr;
 
   be2 = betaGammaSq/(1 + betaGammaSq);
   be4 = be2*be2;
 
   if( betaGammaSq < 0.01 ) logarithm = log(1.0+betaGammaSq); // 0.0;
   else
   {
     logarithm  = -log( (1/betaGammaSq - fRePartDielectricConst[i])*
                    (1/betaGammaSq - fRePartDielectricConst[i]) + 
                    fImPartDielectricConst[i]*fImPartDielectricConst[i] )*0.5;
     logarithm += log(1+1.0/betaGammaSq);
   }
 
   if( fImPartDielectricConst[i] == 0.0 || betaGammaSq < 0.01 )
   {
     argument = 0.0;
   }
   else
   {
     x3 = -fRePartDielectricConst[i] + 1.0/betaGammaSq;
     x5 = -1.0 - fRePartDielectricConst[i] +
          be2*((1.0 +fRePartDielectricConst[i])*(1.0 + fRePartDielectricConst[i]) +
      fImPartDielectricConst[i]*fImPartDielectricConst[i]);
     if( x3 == 0.0 ) argument = 0.5*pi;
     else            argument = atan2(fImPartDielectricConst[i],x3);
     argument *= x5 ;
   }   
   dNdxC = ( logarithm*fImPartDielectricConst[i] + argument )/hbarc;
  
   if(dNdxC < 1.0e-8) dNdxC = 1.0e-8;
 
   dNdxC *= fine_structure_const/be2/pi;
 
   dNdxC *= (1-exp(-be4/betaBohr4));
 
   //   if(fDensity >= 0.1)
   // { 
   modul2 = (1.0 + fRePartDielectricConst[i])*(1.0 + fRePartDielectricConst[i]) + 
                    fImPartDielectricConst[i]*fImPartDielectricConst[i];
   if(modul2 > 0.)
     {
       dNdxC /= modul2;
     }
   return dNdxC;
 
} // end of PAIdNdxCerenkov 

Referenced by Initialize(), InitPAI(), NormShift(), and SplainPAI().

PAIdNdxPlasmon()

G4double G4PAIySection::PAIdNdxPlasmon	(	G4int	intervalNumber,
		G4double	betaGammaSq
	)

Definition at line 688 of file G4PAIySection.cc.

{        
   G4double cof, resonance, modul2, dNdxP;
   G4double be2, be4, betaBohr2, betaBohr4, cofBetaBohr;
 
   cof = 1;
   cofBetaBohr = 4.0;
   betaBohr2   = fine_structure_const*fine_structure_const;
   betaBohr4   = betaBohr2*betaBohr2*cofBetaBohr;
 
   be2 = betaGammaSq/(1 + betaGammaSq);
   be4 = be2*be2;
 
   resonance = log(2*electron_mass_c2*be2/fSplineEnergy[i]);  
   resonance *= fImPartDielectricConst[i]/hbarc;
 
 
   dNdxP = ( resonance + cof*fIntegralTerm[i]/fSplineEnergy[i]/fSplineEnergy[i] );
 
   if( dNdxP < 1.0e-8 ) dNdxP = 1.0e-8;
 
   dNdxP *= fine_structure_const/be2/pi;
   dNdxP *= (1-exp(-be4/betaBohr4));
 
//   if( fDensity >= 0.1 )
//   { 
   modul2 = (1 + fRePartDielectricConst[i])*(1 + fRePartDielectricConst[i]) + 
     fImPartDielectricConst[i]*fImPartDielectricConst[i];
   if(modul2 > 0.)
     { 
       dNdxP /= modul2;
     }
   return dNdxP;
 
} // end of PAIdNdxPlasmon 

Referenced by Initialize(), InitPAI(), NormShift(), and SplainPAI().

RePartDielectricConst()

G4double G4PAIySection::RePartDielectricConst

(

G4double

energy

)

Definition at line 504 of file G4PAIySection.cc.

{       
   G4double x0, x02, x03, x04, x05, x1, x2, xx1 ,xx2 , xx12,
            c1, c2, c3, cof1, cof2, xln1, xln2, xln3, result;
 
   x0 = enb;
   result = 0;
   
   for(G4int i=1;i<=fIntervalNumber-1;i++)
   {
      x1 = fEnergyInterval[i];
      x2 = fEnergyInterval[i+1];
      xx1 = x1 - x0;
      xx2 = x2 - x0;
      xx12 = xx2/xx1;
      
      if(xx12<0)
      {
     xx12 = -xx12;
      }
      xln1 = log(x2/x1);
      xln2 = log(xx12);
      xln3 = log((x2 + x0)/(x1 + x0));
      x02 = x0*x0;
      x03 = x02*x0;
      x04 = x03*x0;
      x05 = x04*x0;
      c1  = (x2 - x1)/x1/x2;
      c2  = (x2 - x1)*(x2 +x1)/x1/x1/x2/x2;
      c3  = (x2 -x1)*(x1*x1 + x1*x2 + x2*x2)/x1/x1/x1/x2/x2/x2;
 
      result -= (fA1[i]/x02 + fA3[i]/x04)*xln1;
      result -= (fA2[i]/x02 + fA4[i]/x04)*c1;
      result -= fA3[i]*c2/2/x02;
      result -= fA4[i]*c3/3/x02;
 
      cof1 = fA1[i]/x02 + fA3[i]/x04;
      cof2 = fA2[i]/x03 + fA4[i]/x05;
 
      result += 0.5*(cof1 +cof2)*xln2;
      result += 0.5*(cof1 - cof2)*xln3;
   } 
   result *= 2*hbarc/pi;
   
   return result;
 
}   // end of RePartDielectricConst 

Referenced by NormShift(), and SplainPAI().

RutherfordIntegral()

G4double G4PAIySection::RutherfordIntegral	(	G4int	intervalNumber,
		G4double	limitLow,
		G4double	limitHigh
	)

Definition at line 460 of file G4PAIySection.cc.

{
   G4double  c1, c2, c3;
   // G4cout<<"RI: x1 = "<<x1<<"; "<<"x2 = "<<x2<<G4endl;   
   c1 = (x2 - x1)/x1/x2;
   c2 = (x2 - x1)*(x2 + x1)/x1/x1/x2/x2;
   c3 = (x2 - x1)*(x1*x1 + x1*x2 + x2*x2)/x1/x1/x1/x2/x2/x2;
   // G4cout<<" RI: c1 = "<<c1<<"; "<<"c2 = "<<c2<<"; "<<"c3 = "<<c3<<G4endl;   
   
   return  fA1[k]*log(x2/x1) + fA2[k]*c1 + fA3[k]*c2/2 + fA4[k]*c3/3;
 
}   // end of RutherfordIntegral 

Referenced by NormShift(), and SplainPAI().

SetVerbose()

void G4PAIySection::SetVerbose ( G4int  v )

inline

Definition at line 144 of file G4PAIySection.hh.

{ fVerbose = v; };

Referenced by G4PAIModel::Initialise().

SplainPAI()

void G4PAIySection::SplainPAI

(

G4double

betaGammaSq

)

Definition at line 377 of file G4PAIySection.cc.

{
   G4int k = 1;
   G4int i = 1;
 
   while ( (i < fSplineNumber) && (fSplineNumber < fMaxSplineSize-1) )
   {
      if(fSplineEnergy[i+1] > fEnergyInterval[k+1])
      {
          k++;   // Here next energy point is in next energy interval
      i++;
          continue;
      }
                   // Shifting of arrayes for inserting the geometrical 
               // average of 'i' and 'i+1' energy points to 'i+1' place
      fSplineNumber++;
 
      for(G4int j = fSplineNumber; j >= i+2; j-- )
      {
         fSplineEnergy[j]          = fSplineEnergy[j-1];
         fImPartDielectricConst[j] = fImPartDielectricConst[j-1];
     fRePartDielectricConst[j] = fRePartDielectricConst[j-1];
     fIntegralTerm[j]          = fIntegralTerm[j-1];
 
     fDifPAIySection[j] = fDifPAIySection[j-1];
         fdNdxCerenkov[j]   = fdNdxCerenkov[j-1];
         fdNdxPlasmon[j]    = fdNdxPlasmon[j-1];
      }
      G4double x1  = fSplineEnergy[i];
      G4double x2  = fSplineEnergy[i+1];
      G4double yy1 = fDifPAIySection[i];
      G4double y2  = fDifPAIySection[i+1];
 
      G4double en1 = sqrt(x1*x2);
      fSplineEnergy[i+1] = en1;
 
         // Calculation of logarithmic linear approximation
         // in this (enr) energy point, which number is 'i+1' now
 
      G4double a = log10(y2/yy1)/log10(x2/x1);
      G4double b = log10(yy1) - a*log10(x1);
      G4double y = a*log10(en1) + b;
      y = pow(10.,y);
 
         // Calculation of the PAI dif. cross-section at this point
 
      fImPartDielectricConst[i+1] = fNormalizationCof*
                                ImPartDielectricConst(k,fSplineEnergy[i+1]);
      fRePartDielectricConst[i+1] = fNormalizationCof*
                                RePartDielectricConst(fSplineEnergy[i+1]);
      fIntegralTerm[i+1] = fIntegralTerm[i] + fNormalizationCof*
                       RutherfordIntegral(k,fSplineEnergy[i],
                                                fSplineEnergy[i+1]);
 
      fDifPAIySection[i+1] = DifPAIySection(i+1,betaGammaSq);
      fdNdxCerenkov[i+1]   = PAIdNdxCerenkov(i+1,betaGammaSq);
      fdNdxPlasmon[i+1]    = PAIdNdxPlasmon(i+1,betaGammaSq);
 
          // Condition for next division of this segment or to pass
          // to higher energies
 
      G4double x = 2*(fDifPAIySection[i+1] - y)/(fDifPAIySection[i+1] + y);
 
      if( x < 0 ) 
      {
     x = -x;
      }
      if( x > fError && fSplineNumber < fMaxSplineSize-1 )
      {
     continue;  // next division
      }
      i += 2;  // pass to next segment
 
   }   // close 'while'
 
}  // end of SplainPAI 

Referenced by Initialize(), and InitPAI().

SumOverBordCerenkov()

G4double G4PAIySection::SumOverBordCerenkov	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 1078 of file G4PAIySection.cc.

{               
   G4double x0,x1,y0,yy1,a,e0,c,d,result;
 
   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fdNdxCerenkov[i];
   yy1 = fdNdxCerenkov[i+1];
 
   //  G4cout<<G4endl;
   //G4cout<<"SumBordC, i = "<<i<<"; en0 = "<<en0<<"; x0 ="<<x0<<"; x1 = "<<x1
   //     <<"; y0 = "<<y0<<"; yy1 = "<<yy1<<G4endl;
   c = x1/x0;
   d = e0/x0;
   a = log10(yy1/y0)/log10(c);
 
   G4double b = 0.0;
   if(a < 20.) b = y0/pow(x0,a);
   
   a += 1.0;
   if( a == 0 ) result = b*log(x0/e0);
   else         result = y0*(x0 - e0*pow(d,a-1))/a;   
   a += 1.0;
 
   if( a == 0 ) fIntegralCerenkov[0] += b*log(x0/e0);
   else         fIntegralCerenkov[0] += y0*(x0*x0 - e0*e0*pow(d,a-2))/a;
 
   //G4cout<<"a = "<<a<<"; b = "<<b<<"; result = "<<result<<G4endl;
   
   x0  = fSplineEnergy[i - 1];
   x1  = fSplineEnergy[i - 2];
   y0  = fdNdxCerenkov[i - 1];
   yy1 = fdNdxCerenkov[i - 2];
 
   //G4cout<<"x0 ="<<x0<<"; x1 = "<<x1
   //    <<"; y0 = "<<y0<<"; yy1 = "<<yy1<<G4endl;
 
   c = x1/x0;
   d = e0/x0;
   a  = log10(yy1/y0)/log10(x1/x0);
  
   //   G4cout << "a= " << a << G4endl;
   if(a < 20.) b = y0/pow(x0,a);
 
   if(a > 20.0) b = 0.0;
   else         b = y0/pow(x0,a);  // pow(10.,b0);
 
   //G4cout << "b= " << b << G4endl;
 
   a += 1.0;
   if( a == 0 ) result += b*log(e0/x0);
   else         result += y0*(e0*pow(d,a-1) - x0 )/a;
   a += 1.0;
   //G4cout << "result= " << result << G4endl;
 
   if( a == 0 )   fIntegralCerenkov[0] += b*log(e0/x0);
   else           fIntegralCerenkov[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
 
   //G4cout<<"a = "<<a<<"; b = "<<b<<"; result = "<<result<<G4endl;    
 
   return result;
 
} 

Referenced by IntegralCerenkov().

SumOverBorder()

G4double G4PAIySection::SumOverBorder	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 952 of file G4PAIySection.cc.

{               
  G4double x0,x1,y0,yy1,a,/*c,*/d,e0,result;
 
   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fDifPAIySection[i];
   yy1 = fDifPAIySection[i+1];
 
   //c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(x1/x0);
 
   G4double b = 0.0;
   if(a < 20.) b = y0/pow(x0,a);
   
   a += 1;
   if(a == 0)
   {
      result = b*log(x0/e0);
   }
   else
   {
      result = y0*(x0 - e0*pow(d,a-1))/a;
   }
   a++;
   if(a == 0)
   {
      fIntegralPAIySection[0] += b*log(x0/e0);
   }
   else 
   {
      fIntegralPAIySection[0] += y0*(x0*x0 - e0*e0*pow(d,a-2))/a;
   }
   x0 = fSplineEnergy[i - 1];
   x1 = fSplineEnergy[i - 2];
   y0 = fDifPAIySection[i - 1];
   yy1 = fDifPAIySection[i - 2];
 
   //c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(x1/x0);
   //  b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);
   a += 1;
   if(a == 0)
   {
      result += b*log(e0/x0);
   }
   else
   {
      result += y0*(e0*pow(d,a-1) - x0)/a;
   }
   a++;
   if(a == 0) 
   {
      fIntegralPAIySection[0] += b*log(e0/x0);
   }
   else
   {
      fIntegralPAIySection[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
   }
   return result;
 
} 

Referenced by IntegralPAIySection().

SumOverBorderdEdx()

G4double G4PAIySection::SumOverBorderdEdx	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 1022 of file G4PAIySection.cc.

{               
  G4double x0,x1,y0,yy1,a,/*c,*/d,e0,result;
 
   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fDifPAIySection[i];
   yy1 = fDifPAIySection[i+1];
 
   //c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(x1/x0);
   
   G4double b = 0.0;
   if(a < 20.) b = y0/pow(x0,a);
   
   a += 2;
   if(a == 0)
   {
      result = b*log(x0/e0);
   }
   else 
   {
      result = y0*(x0*x0 - e0*e0*pow(d,a-2))/a;
   }
   x0 = fSplineEnergy[i - 1];
   x1 = fSplineEnergy[i - 2];
   y0 = fDifPAIySection[i - 1];
   yy1 = fDifPAIySection[i - 2];
 
   //c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(x1/x0);
 
   if(a < 20.) b = y0/pow(x0,a);
 
   a += 2;
   if(a == 0) 
   {
      result += b*log(e0/x0);
   }
   else
   {
      result += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
   }
   return result;
 
} 

Referenced by IntegralPAIySection().

SumOverBordPlasmon()

G4double G4PAIySection::SumOverBordPlasmon	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 1149 of file G4PAIySection.cc.

{               
   G4double x0,x1,y0,yy1,a,c,d,e0,result;
 
   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fdNdxPlasmon[i];
   yy1 = fdNdxPlasmon[i+1];
 
   c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(c);
 
   G4double b = 0.0;
   if(a < 20.) b = y0/pow(x0,a);
   
   a += 1.0;
   if( a == 0 ) result = b*log(x0/e0);
   else         result = y0*(x0 - e0*pow(d,a-1))/a;   
   a += 1.0;
 
   if( a == 0 ) fIntegralPlasmon[0] += b*log(x0/e0);
   else         fIntegralPlasmon[0] += y0*(x0*x0 - e0*e0*pow(d,a-2))/a;
   
   x0 = fSplineEnergy[i - 1];
   x1 = fSplineEnergy[i - 2];
   y0 = fdNdxPlasmon[i - 1];
   yy1 = fdNdxPlasmon[i - 2];
 
   c = x1/x0;
   d = e0/x0;
   a = log10(yy1/y0)/log10(c);
 
   if(a < 20.) b = y0/pow(x0,a);
 
   a += 1.0;
   if( a == 0 ) result += b*log(e0/x0);
   else         result += y0*(e0*pow(d,a-1) - x0)/a;
   a += 1.0;
 
   if( a == 0 )   fIntegralPlasmon[0] += b*log(e0/x0);
   else           fIntegralPlasmon[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
   
   return result;
 
} 

Referenced by IntegralPlasmon().

SumOverInterCerenkov()

G4double G4PAIySection::SumOverInterCerenkov

(

G4int

intervalNumber

)

Definition at line 887 of file G4PAIySection.cc.

{         
   G4double x0,x1,y0,yy1,a,c,result;
 
   x0  = fSplineEnergy[i];
   x1  = fSplineEnergy[i+1];
   y0  = fdNdxCerenkov[i];
   yy1 = fdNdxCerenkov[i+1];
   // G4cout<<"SumC, i = "<<i<<"; x0 ="<<x0<<"; x1 = "<<x1
   //   <<"; y0 = "<<y0<<"; yy1 = "<<yy1<<G4endl;
 
   c = x1/x0;
   a = log10(yy1/y0)/log10(c);
   G4double b = 0.0;
   if(a < 20.) b = y0/pow(x0,a);
 
   a += 1.0;
   if(a == 0) result = b*log(c);
   else       result = y0*(x1*pow(c,a-1) - x0)/a;   
   a += 1.0;
 
   if( a == 0 ) fIntegralCerenkov[0] += b*log(x1/x0);
   else         fIntegralCerenkov[0] += y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   //  G4cout<<"a = "<<a<<"; b = "<<b<<"; result = "<<result<<G4endl;   
   return result;
 
} //  end of SumOverInterCerenkov

Referenced by IntegralCerenkov().

SumOverInterPlasmon()

G4double G4PAIySection::SumOverInterPlasmon

(

G4int

intervalNumber

)

Definition at line 921 of file G4PAIySection.cc.

{         
  G4double x0,x1,y0,yy1,a,c,result;
 
   x0  = fSplineEnergy[i];
   x1  = fSplineEnergy[i+1];
   y0  = fdNdxPlasmon[i];
   yy1 = fdNdxPlasmon[i+1];
   c = x1/x0;
   a = log10(yy1/y0)/log10(c);
 
   G4double b = 0.0;
   if(a < 20.) b = y0/pow(x0,a);
 
   a += 1.0;
   if(a == 0) result = b*log(x1/x0);
   else       result = y0*(x1*pow(c,a-1) - x0)/a;   
   a += 1.0;
 
   if( a == 0 ) fIntegralPlasmon[0] += b*log(x1/x0);
   else         fIntegralPlasmon[0] += y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   
   return result;
 
} //  end of SumOverInterPlasmon

Referenced by IntegralPlasmon().

SumOverInterval()

G4double G4PAIySection::SumOverInterval

(

G4int

intervalNumber

)

Definition at line 820 of file G4PAIySection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;
 
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fDifPAIySection[i];
   yy1 = fDifPAIySection[i+1];
   c = x1/x0;
   a = log10(yy1/y0)/log10(c);
   // b = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);
   a += 1;
   if(a == 0) 
   {
      result = b*log(x1/x0);
   }
   else
   {
      result = y0*(x1*pow(c,a-1) - x0)/a;
   }
   a++;
   if(a == 0) 
   {
      fIntegralPAIySection[0] += b*log(x1/x0);
   }
   else
   {
      fIntegralPAIySection[0] += y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   }
   return result;
 
} //  end of SumOverInterval

Referenced by IntegralPAIySection().

SumOverIntervaldEdx()

G4double G4PAIySection::SumOverIntervaldEdx

(

G4int

intervalNumber

)

Definition at line 856 of file G4PAIySection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;
 
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fDifPAIySection[i];
   yy1 = fDifPAIySection[i+1];
   c = x1/x0;
   a = log10(yy1/y0)/log10(c);
   // b = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);
   a += 2;
   if(a == 0) 
   {
     result = b*log(x1/x0);
   }
   else
   {
     result = y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   }
   return result;
 
} //  end of SumOverInterval

Referenced by IntegralPAIySection().

---

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

• G4PAIySection.hh
• G4PAIySection.cc