G4PAIxSection Class Reference

#include <G4PAIxSection.hh>

Collaboration diagram for G4PAIxSection:

Public Member Functions
	G4PAIxSection ()
	G4PAIxSection (G4MaterialCutsCouple *matCC)
	G4PAIxSection (G4int materialIndex, G4double maxEnergyTransfer)
	G4PAIxSection (G4int materialIndex, G4double maxEnergyTransfer, G4double betaGammaSq, G4double **photoAbsCof, G4int intNumber)
	G4PAIxSection (G4int materialIndex, G4double maxEnergyTransfer, G4double betaGammaSq)
	~G4PAIxSection ()
void	Initialize (const G4Material *material, G4double maxEnergyTransfer, G4double betaGammaSq, G4SandiaTable *)
void	ComputeLowEnergyCof (const G4Material *material)
void	ComputeLowEnergyCof ()
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	GetPhotonRange (G4double energy)
G4double	GetElectronRange (G4double energy)
G4double	RePartDielectricConst (G4double energy)
G4double	DifPAIxSection (G4int intervalNumber, G4double betaGammaSq)
G4double	PAIdNdxCerenkov (G4int intervalNumber, G4double betaGammaSq)
G4double	PAIdNdxMM (G4int intervalNumber, G4double betaGammaSq)
G4double	PAIdNdxPlasmon (G4int intervalNumber, G4double betaGammaSq)
G4double	PAIdNdxResonance (G4int intervalNumber, G4double betaGammaSq)
void	IntegralPAIxSection ()
void	IntegralCerenkov ()
void	IntegralMM ()
void	IntegralPlasmon ()
void	IntegralResonance ()
G4double	SumOverInterval (G4int intervalNumber)
G4double	SumOverIntervaldEdx (G4int intervalNumber)
G4double	SumOverInterCerenkov (G4int intervalNumber)
G4double	SumOverInterMM (G4int intervalNumber)
G4double	SumOverInterPlasmon (G4int intervalNumber)
G4double	SumOverInterResonance (G4int intervalNumber)
G4double	SumOverBorder (G4int intervalNumber, G4double energy)
G4double	SumOverBorderdEdx (G4int intervalNumber, G4double energy)
G4double	SumOverBordCerenkov (G4int intervalNumber, G4double energy)
G4double	SumOverBordMM (G4int intervalNumber, G4double energy)
G4double	SumOverBordPlasmon (G4int intervalNumber, G4double energy)
G4double	SumOverBordResonance (G4int intervalNumber, G4double energy)
G4double	GetStepEnergyLoss (G4double step)
G4double	GetStepCerenkovLoss (G4double step)
G4double	GetStepMMLoss (G4double step)
G4double	GetStepPlasmonLoss (G4double step)
G4double	GetStepResonanceLoss (G4double step)
G4double	GetEnergyTransfer ()
G4double	GetCerenkovEnergyTransfer ()
G4double	GetMMEnergyTransfer ()
G4double	GetPlasmonEnergyTransfer ()
G4double	GetResonanceEnergyTransfer ()
G4double	GetRutherfordEnergyTransfer ()
G4int	GetNumberOfGammas () const
G4int	GetSplineSize () const
G4int	GetIntervalNumber () const
G4double	GetEnergyInterval (G4int i)
G4double	GetDifPAIxSection (G4int i)
G4double	GetPAIdNdxCerenkov (G4int i)
G4double	GetPAIdNdxMM (G4int i)
G4double	GetPAIdNdxPlasmon (G4int i)
G4double	GetPAIdNdxResonance (G4int i)
G4double	GetMeanEnergyLoss () const
G4double	GetMeanCerenkovLoss () const
G4double	GetMeanMMLoss () const
G4double	GetMeanPlasmonLoss () const
G4double	GetMeanResonanceLoss () const
G4double	GetNormalizationCof () const
G4double	GetLowEnergyCof () const
void	SetVerbose (G4int v)
G4double	GetPAItable (G4int i, G4int j) const
G4double	GetLorentzFactor (G4int i) const
G4double	GetSplineEnergy (G4int i) const
G4double	GetIntegralPAIxSection (G4int i) const
G4double	GetIntegralPAIdEdx (G4int i) const
G4double	GetIntegralCerenkov (G4int i) const
G4double	GetIntegralMM (G4int i) const
G4double	GetIntegralPlasmon (G4int i) const
G4double	GetIntegralResonance (G4int i) const

Private Member Functions
void	CallError (G4int i, const G4String &methodName) const
G4PAIxSection &	operator= (const G4PAIxSection &right)
	G4PAIxSection (const G4PAIxSection &)

Private Attributes
G4int	fIntervalNumber
G4double	fNormalizationCof
G4int	fMaterialIndex
G4double	fDensity
G4double	fElectronDensity
G4double	fLowEnergyCof
G4int	fSplineNumber
G4int	fVerbose
G4OrderedTable *	fMatSandiaMatrix
G4SandiaTable *	fSandia
G4DataVector	fEnergyInterval
G4DataVector	fA1
G4DataVector	fA2
G4DataVector	fA3
G4DataVector	fA4
G4DataVector	fSplineEnergy
G4DataVector	fRePartDielectricConst
G4DataVector	fImPartDielectricConst
G4DataVector	fIntegralTerm
G4DataVector	fDifPAIxSection
G4DataVector	fdNdxCerenkov
G4DataVector	fdNdxPlasmon
G4DataVector	fdNdxMM
G4DataVector	fdNdxResonance
G4DataVector	fIntegralPAIxSection
G4DataVector	fIntegralPAIdEdx
G4DataVector	fIntegralCerenkov
G4DataVector	fIntegralPlasmon
G4DataVector	fIntegralMM
G4DataVector	fIntegralResonance
G4double	fPAItable [500][112]

Static Private Attributes
static const G4double	fDelta = 0.005
static const G4double	fError = 0.005
static G4int	fNumberOfGammas = 111
static const G4double	fLorentzFactor [112]
static const G4int	fRefGammaNumber = 29
static const G4int	fMaxSplineSize = 1000

Detailed Description

Definition at line 68 of file G4PAIxSection.hh.

Constructor & Destructor Documentation

G4PAIxSection() [1/6]

G4PAIxSection::G4PAIxSection

(

)

Definition at line 92 of file G4PAIxSection.cc.

{
  fSandia = 0;
  fMatSandiaMatrix = 0;
  fDensity = fElectronDensity = fNormalizationCof = fLowEnergyCof = 0.0;
  fIntervalNumber = fSplineNumber = 0;
  fVerbose = 0;
    
  fSplineEnergy          = G4DataVector(fMaxSplineSize,0.0);
  fRePartDielectricConst = G4DataVector(fMaxSplineSize,0.0);
  fImPartDielectricConst = G4DataVector(fMaxSplineSize,0.0);
  fIntegralTerm          = G4DataVector(fMaxSplineSize,0.0);
  fDifPAIxSection        = G4DataVector(fMaxSplineSize,0.0);
  fdNdxCerenkov          = G4DataVector(fMaxSplineSize,0.0);
  fdNdxPlasmon           = G4DataVector(fMaxSplineSize,0.0);
  fdNdxMM                = G4DataVector(fMaxSplineSize,0.0);
  fdNdxResonance         = G4DataVector(fMaxSplineSize,0.0);
  fIntegralPAIxSection   = G4DataVector(fMaxSplineSize,0.0);
  fIntegralPAIdEdx       = G4DataVector(fMaxSplineSize,0.0);
  fIntegralCerenkov      = G4DataVector(fMaxSplineSize,0.0);
  fIntegralPlasmon       = G4DataVector(fMaxSplineSize,0.0);
  fIntegralMM            = G4DataVector(fMaxSplineSize,0.0);
  fIntegralResonance     = G4DataVector(fMaxSplineSize,0.0);

  fMaterialIndex = 0;   

  for( G4int i = 0; i < 500; ++i ) 
  {
    for( G4int j = 0; j < 112; ++j )  fPAItable[i][j] = 0.0; 
  }
}

Here is the caller graph for this function:

G4PAIxSection() [2/6]

G4PAIxSection::G4PAIxSection

(

G4MaterialCutsCouple *

matCC

)

Definition at line 129 of file G4PAIxSection.cc.

{
  fDensity       = matCC->GetMaterial()->GetDensity();
  G4int matIndex = matCC->GetMaterial()->GetIndex();
  fMaterialIndex = matIndex;   

  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
  fSandia = (*theMaterialTable)[matIndex]->GetSandiaTable();

  fVerbose = 0;

  G4int i, j; 
  fMatSandiaMatrix = new G4OrderedTable();
 
  for (i = 0; i < fSandia->GetMaxInterval()-1; i++)
  {
     fMatSandiaMatrix->push_back(new G4DataVector(5,0.));
  }                         
  for (i = 0; i < fSandia->GetMaxInterval()-1; i++)
  {
    (*(*fMatSandiaMatrix)[i])[0] = fSandia->GetSandiaMatTable(i,0);

    for(j = 1; j < 5; j++)
    {
      (*(*fMatSandiaMatrix)[i])[j] = fSandia->GetSandiaMatTable(i,j)*fDensity;
    }     
  }
  ComputeLowEnergyCof();                               
  //  fEnergyInterval = fA1 = fA2 = fA3 = fA4 = 0;
}

Here is the call graph for this function:

G4PAIxSection() [3/6]

G4PAIxSection::G4PAIxSection	(	G4int	materialIndex,
		G4double	maxEnergyTransfer
	)

Definition at line 162 of file G4PAIxSection.cc.

{
  fSandia = 0;
  fMatSandiaMatrix = 0;
  fVerbose = 0;
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
  G4int i, j;   

  fMaterialIndex   = materialIndex;   
  fDensity                = (*theMaterialTable)[materialIndex]->GetDensity();
  fElectronDensity        = (*theMaterialTable)[materialIndex]->
                             GetElectronDensity();
  fIntervalNumber         = (*theMaterialTable)[materialIndex]->
                             GetSandiaTable()->GetMatNbOfIntervals();
  fIntervalNumber--;      
  // G4cout<<fDensity<<"\t"<<fElectronDensity<<"\t"<<fIntervalNumber<<G4endl;

  fEnergyInterval = G4DataVector(fIntervalNumber+2,0.0);
  fA1             = G4DataVector(fIntervalNumber+2,0.0);
  fA2             = G4DataVector(fIntervalNumber+2,0.0);
  fA3             = G4DataVector(fIntervalNumber+2,0.0);
  fA4             = G4DataVector(fIntervalNumber+2,0.0);

  for(i = 1; i <= fIntervalNumber; i++ )
    {
      if(((*theMaterialTable)[materialIndex]->
    GetSandiaTable()->GetSandiaCofForMaterial(i-1,0) >= maxEnergyTransfer) ||
              i > fIntervalNumber               )
        {
          fEnergyInterval[i] = maxEnergyTransfer;
          fIntervalNumber = i;
          break;
        }
         fEnergyInterval[i] = (*theMaterialTable)[materialIndex]->
                              GetSandiaTable()->GetSandiaCofForMaterial(i-1,0);
         fA1[i]             = (*theMaterialTable)[materialIndex]->
                              GetSandiaTable()->GetSandiaCofForMaterial(i-1,1);
         fA2[i]             = (*theMaterialTable)[materialIndex]->
                              GetSandiaTable()->GetSandiaCofForMaterial(i-1,2);
         fA3[i]             = (*theMaterialTable)[materialIndex]->
                              GetSandiaTable()->GetSandiaCofForMaterial(i-1,3);
         fA4[i]             = (*theMaterialTable)[materialIndex]->
                              GetSandiaTable()->GetSandiaCofForMaterial(i-1,4);
         // G4cout<<i<<"\t"<<fEnergyInterval[i]<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
         //                               <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }   
  if(fEnergyInterval[fIntervalNumber] != maxEnergyTransfer)
    {
         fIntervalNumber++;
         fEnergyInterval[fIntervalNumber] = maxEnergyTransfer;
    }

  // Now checking, if two borders are too close together

  for(i=1;i<fIntervalNumber;i++)
    {
        if(fEnergyInterval[i+1]-fEnergyInterval[i] >
           1.5*fDelta*(fEnergyInterval[i+1]+fEnergyInterval[i]))
        {
          continue;
        }
        else
        {
          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--;
        }
    }


      /* *********************************

      fSplineEnergy          = new G4double[fMaxSplineSize];   
      fRePartDielectricConst = new G4double[fMaxSplineSize];   
      fImPartDielectricConst = new G4double[fMaxSplineSize];   
      fIntegralTerm          = new G4double[fMaxSplineSize];   
      fDifPAIxSection        = new G4double[fMaxSplineSize];   
      fIntegralPAIxSection   = new G4double[fMaxSplineSize];   
      
      for(i=0;i<fMaxSplineSize;i++)
      {
         fSplineEnergy[i]          = 0.0;   
         fRePartDielectricConst[i] = 0.0;   
         fImPartDielectricConst[i] = 0.0;   
         fIntegralTerm[i]          = 0.0;   
         fDifPAIxSection[i]        = 0.0;   
         fIntegralPAIxSection[i]   = 0.0;   
      }
      **************************************************  */   
      ComputeLowEnergyCof();      
      InitPAI();  // create arrays allocated above
      /*     
      delete[] fEnergyInterval;
      delete[] fA1;
      delete[] fA2;
      delete[] fA3;
      delete[] fA4; 
      */   
}

Here is the call graph for this function:

G4PAIxSection() [4/6]

G4PAIxSection::G4PAIxSection	(	G4int	materialIndex,
		G4double	maxEnergyTransfer,
		G4double	betaGammaSq,
		G4double **	photoAbsCof,
		G4int	intNumber
	)

Definition at line 274 of file G4PAIxSection.cc.

{
  fSandia = 0;
  fDensity = fElectronDensity = fNormalizationCof = fLowEnergyCof = 0.0;
  fIntervalNumber = fSplineNumber = 0;
  fVerbose = 0;
    
  fSplineEnergy          = G4DataVector(500,0.0);
  fRePartDielectricConst = G4DataVector(500,0.0);
  fImPartDielectricConst = G4DataVector(500,0.0);
  fIntegralTerm          = G4DataVector(500,0.0);
  fDifPAIxSection        = G4DataVector(500,0.0);
  fdNdxCerenkov          = G4DataVector(500,0.0);
  fdNdxPlasmon           = G4DataVector(500,0.0);
  fdNdxMM                = G4DataVector(500,0.0);
  fdNdxResonance         = G4DataVector(500,0.0);
  fIntegralPAIxSection   = G4DataVector(500,0.0);
  fIntegralPAIdEdx       = G4DataVector(500,0.0);
  fIntegralCerenkov      = G4DataVector(500,0.0);
  fIntegralPlasmon       = G4DataVector(500,0.0);
  fIntegralMM            = G4DataVector(500,0.0);
  fIntegralResonance     = G4DataVector(500,0.0);

  for( G4int i = 0; i < 500; ++i ) 
  {
    for( G4int j = 0; j < 112; ++j )  fPAItable[i][j] = 0.0; 
  }

  fSandia = 0;
  fMatSandiaMatrix = 0;
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
  G4int i, j; 
  
  fMaterialIndex   = materialIndex;      
  fDensity         = (*theMaterialTable)[materialIndex]->GetDensity();
  fElectronDensity = (*theMaterialTable)[materialIndex]->GetElectronDensity();

  fIntervalNumber         = intNumber;
  fIntervalNumber--;
  //   G4cout<<fDensity<<"\t"<<fElectronDensity<<"\t"<<fIntervalNumber<<G4endl;
  
  fEnergyInterval = G4DataVector(fIntervalNumber+2,0.0);
  fA1             = G4DataVector(fIntervalNumber+2,0.0);
  fA2             = G4DataVector(fIntervalNumber+2,0.0);
  fA3             = G4DataVector(fIntervalNumber+2,0.0);
  fA4             = G4DataVector(fIntervalNumber+2,0.0);


  /*
      fEnergyInterval = new G4double[fIntervalNumber+2];
      fA1             = new G4double[fIntervalNumber+2];
      fA2             = new G4double[fIntervalNumber+2];
      fA3             = new G4double[fIntervalNumber+2];
      fA4             = new G4double[fIntervalNumber+2];
  */
  for( i = 1; i <= fIntervalNumber; i++ )
    {
         if( ( photoAbsCof[i-1][0] >= maxEnergyTransfer ) ||
             i > fIntervalNumber )
         {
            fEnergyInterval[i] = maxEnergyTransfer;
            fIntervalNumber = i;
            break;
         }
         fEnergyInterval[i] = photoAbsCof[i-1][0];
         fA1[i]             = photoAbsCof[i-1][1];
         fA2[i]             = photoAbsCof[i-1][2];
         fA3[i]             = photoAbsCof[i-1][3];
         fA4[i]             = photoAbsCof[i-1][4];
         // G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
         //      <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }
      // G4cout<<"i last = "<<i<<"; "<<"fIntervalNumber = "<<fIntervalNumber<<G4endl; 
  
  if(fEnergyInterval[fIntervalNumber] != maxEnergyTransfer)
    {
         fIntervalNumber++;
         fEnergyInterval[fIntervalNumber] = maxEnergyTransfer;
    }
      // G4cout<<"after check of max energy transfer"<<G4endl;

  for( i = 1; i <= fIntervalNumber; i++ )
    {
        // G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
        //   <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }
      // Now checking, if two borders are too close together

  for( i = 1; i < fIntervalNumber; i++ )
    {
        if(fEnergyInterval[i+1]-fEnergyInterval[i] >
           1.5*fDelta*(fEnergyInterval[i+1]+fEnergyInterval[i]))
        {
          continue;
        }
        else
        {
          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--;
        }
    }
  // G4cout<<"after check of close borders"<<G4endl;

  for( i = 1; i <= fIntervalNumber; i++ )
    {
        // G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
        //  <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }

  // Preparation of fSplineEnergy array corresponding to min ionisation, G~4

  ComputeLowEnergyCof();            
  G4double   betaGammaSqRef = 
    fLorentzFactor[fRefGammaNumber]*fLorentzFactor[fRefGammaNumber] - 1;

  NormShift(betaGammaSqRef);             
  SplainPAI(betaGammaSqRef);
      
  // Preparation of integral PAI cross section for input betaGammaSq
   
  for(i = 1; i <= fSplineNumber; i++)
    {
         fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
         fdNdxMM[i]   = PAIdNdxMM(i,betaGammaSq);
         fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
         fdNdxResonance[i]  = PAIdNdxResonance(i,betaGammaSq);
         fDifPAIxSection[i] = DifPAIxSection(i,betaGammaSq);

         // G4cout<<i<<"; dNdxC = "<<fdNdxCerenkov[i]<<"; dNdxP = "<<fdNdxPlasmon[i]
         //    <<"; dNdxPAI = "<<fDifPAIxSection[i]<<G4endl;
    }
  IntegralCerenkov();
  IntegralMM();
  IntegralPlasmon();
  IntegralResonance();
  IntegralPAIxSection();
  /*      
      delete[] fEnergyInterval;
      delete[] fA1;
      delete[] fA2;
      delete[] fA3;
      delete[] fA4;
  */    
}

Here is the call graph for this function:

G4PAIxSection() [5/6]

G4PAIxSection::G4PAIxSection	(	G4int	materialIndex,
		G4double	maxEnergyTransfer,
		G4double	betaGammaSq
	)

Definition at line 435 of file G4PAIxSection.cc.

{
  fSandia = 0;
  fMatSandiaMatrix = 0;
  fVerbose = 0;
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();

  G4int i, j, numberOfElements;   

  fMaterialIndex   = materialIndex;   
  fDensity         = (*theMaterialTable)[materialIndex]->GetDensity();
  fElectronDensity = (*theMaterialTable)[materialIndex]->GetElectronDensity();
  numberOfElements = (*theMaterialTable)[materialIndex]->GetNumberOfElements();

  G4int* thisMaterialZ = new G4int[numberOfElements];
   
  for( i = 0; i < numberOfElements; i++ )
   {
         thisMaterialZ[i] = (G4int)(*theMaterialTable)[materialIndex]->
                                      GetElement(i)->GetZ();
   }
  // fSandia = new G4SandiaTable(materialIndex);
  fSandia = (*theMaterialTable)[materialIndex]->GetSandiaTable();
  G4SandiaTable     thisMaterialSandiaTable(materialIndex);
  fIntervalNumber = thisMaterialSandiaTable.SandiaIntervals(thisMaterialZ,
                                                            numberOfElements);
  fIntervalNumber = thisMaterialSandiaTable.SandiaMixing
                           ( thisMaterialZ ,
                      (*theMaterialTable)[materialIndex]->GetFractionVector() ,
                             numberOfElements,fIntervalNumber);

  fIntervalNumber--;

  fEnergyInterval = G4DataVector(fIntervalNumber+2,0.0);
  fA1             = G4DataVector(fIntervalNumber+2,0.0);
  fA2             = G4DataVector(fIntervalNumber+2,0.0);
  fA3             = G4DataVector(fIntervalNumber+2,0.0);
  fA4             = G4DataVector(fIntervalNumber+2,0.0);

  /*
      fEnergyInterval = new G4double[fIntervalNumber+2];
      fA1             = new G4double[fIntervalNumber+2];
      fA2             = new G4double[fIntervalNumber+2];
      fA3             = new G4double[fIntervalNumber+2];
      fA4             = new G4double[fIntervalNumber+2];
  */
  for( i = 1; i <= fIntervalNumber; i++ )
    {
  if((thisMaterialSandiaTable.GetPhotoAbsorpCof(i,0) >= maxEnergyTransfer) ||
          i > fIntervalNumber)
         {
            fEnergyInterval[i] = maxEnergyTransfer;
            fIntervalNumber = i;
            break;
         }
   fEnergyInterval[i] = thisMaterialSandiaTable.GetPhotoAbsorpCof(i,0);
   fA1[i]             = thisMaterialSandiaTable.GetPhotoAbsorpCof(i,1)*fDensity;
   fA2[i]             = thisMaterialSandiaTable.GetPhotoAbsorpCof(i,2)*fDensity;
   fA3[i]             = thisMaterialSandiaTable.GetPhotoAbsorpCof(i,3)*fDensity;
   fA4[i]             = thisMaterialSandiaTable.GetPhotoAbsorpCof(i,4)*fDensity;

    }   
  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];
    }
  for(i=1;i<=fIntervalNumber;i++)
    {
        // G4cout<<fEnergyInterval[i]<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
        //   <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }
  // Now checking, if two borders are too close together

  for( i = 1; i < fIntervalNumber; i++ )
    {
        if(fEnergyInterval[i+1]-fEnergyInterval[i] >
           1.5*fDelta*(fEnergyInterval[i+1]+fEnergyInterval[i]))
        {
          continue;
        }
        else
        {
          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--;
        }
    }

      /* *********************************
      fSplineEnergy          = new G4double[fMaxSplineSize];   
      fRePartDielectricConst = new G4double[fMaxSplineSize];   
      fImPartDielectricConst = new G4double[fMaxSplineSize];   
      fIntegralTerm          = new G4double[fMaxSplineSize];   
      fDifPAIxSection        = new G4double[fMaxSplineSize];   
      fIntegralPAIxSection   = new G4double[fMaxSplineSize];   
      
      for(i=0;i<fMaxSplineSize;i++)
      {
         fSplineEnergy[i]          = 0.0;   
         fRePartDielectricConst[i] = 0.0;   
         fImPartDielectricConst[i] = 0.0;   
         fIntegralTerm[i]          = 0.0;   
         fDifPAIxSection[i]        = 0.0;   
         fIntegralPAIxSection[i]   = 0.0;   
      }
      */ 

      // Preparation of fSplineEnergy array corresponding to min ionisation, G~4

  ComputeLowEnergyCof();      
  G4double   betaGammaSqRef = 
    fLorentzFactor[fRefGammaNumber]*fLorentzFactor[fRefGammaNumber] - 1;

  NormShift(betaGammaSqRef);             
  SplainPAI(betaGammaSqRef);
      
  // Preparation of integral PAI cross section for input betaGammaSq
   
  for(i = 1; i <= fSplineNumber; i++)
    {
         fDifPAIxSection[i] = DifPAIxSection(i,betaGammaSq);
         fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
         fdNdxMM[i]   = PAIdNdxMM(i,betaGammaSq);
         fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
         fdNdxResonance[i]  = PAIdNdxResonance(i,betaGammaSq);
    }
  IntegralPAIxSection();
  IntegralCerenkov();
  IntegralMM();
  IntegralPlasmon();
  IntegralResonance();    
}

Here is the call graph for this function:

~G4PAIxSection()

G4PAIxSection::~G4PAIxSection

(

)

Definition at line 586 of file G4PAIxSection.cc.

{
   /* ************************
   delete[] fSplineEnergy         ;   
   delete[] fRePartDielectricConst;   
   delete[] fImPartDielectricConst;   
   delete[] fIntegralTerm         ;   
   delete[] fDifPAIxSection       ;   
   delete[] fIntegralPAIxSection  ;
   */ 
  delete fMatSandiaMatrix;
}

G4PAIxSection() [6/6]

G4PAIxSection::G4PAIxSection ( const G4PAIxSection &  )

private

Member Function Documentation

CallError()

void G4PAIxSection::CallError ( G4int  i, const G4String &  methodName  ) const

private

Definition at line 2525 of file G4PAIxSection.cc.

{
  G4String head = "G4PAIxSection::" + methodName + "()";
  G4ExceptionDescription ed;
  ed << "Wrong index " << i << " fSplineNumber= " << fSplineNumber;
  G4Exception(head,"pai001",FatalException,ed);
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeLowEnergyCof() [1/2]

void G4PAIxSection::ComputeLowEnergyCof

(

const G4Material *

material

)

Definition at line 746 of file G4PAIxSection.cc.

{    
  G4int i, numberOfElements = material->GetNumberOfElements();
  G4double sumZ = 0., sumCof = 0.; 

  static const G4double p0 =  1.20923e+00; 
  static const G4double p1 =  3.53256e-01; 
  static 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;
}

Here is the call graph for this function:

ComputeLowEnergyCof() [2/2]

void G4PAIxSection::ComputeLowEnergyCof

(

)

Definition at line 779 of file G4PAIxSection.cc.

{    
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
  G4int i, numberOfElements = (*theMaterialTable)[fMaterialIndex]->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] = (*theMaterialTable)[fMaterialIndex]->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;
  // G4cout<<"fLowEnergyCof = "<<fLowEnergyCof<<G4endl;
  delete [] thisMaterialZ;
  delete [] thisMaterialCof;
}

Here is the call graph for this function:

DifPAIxSection()

G4double G4PAIxSection::DifPAIxSection	(	G4int	intervalNumber,
		G4double	betaGammaSq
	)

Definition at line 1225 of file G4PAIxSection.cc.

{        
   G4double cof,x1,x2,x3,x4,x5,x6,x7,x8,result;

   G4double betaBohr  = fine_structure_const;
   // G4double betaBohr2 = fine_structure_const*fine_structure_const;
   // G4double betaBohr3 = betaBohr*betaBohr2; // *4.0;

   G4double be2  = betaGammaSq/(1 + betaGammaSq);
   G4double beta = sqrt(be2);
   // G4double be3 = beta*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(-be2/betaBohr2/lowCof));

   // result *= (1 - exp(-be3/betaBohr3/lowCof)); // ~ be for be<<betaBohr

   // result *= (1 - exp(-be4/betaBohr4/lowCof));

   if(fDensity >= 0.1)
   { 
      result /= x8;
   }
   return result;

} // end of DifPAIxSection 

GetCerenkovEnergyTransfer()

G4double G4PAIxSection::GetCerenkovEnergyTransfer

(

)

Definition at line 2335 of file G4PAIxSection.cc.

{  
  G4int iTransfer ;

  G4double energyTransfer, position;

  position = fIntegralCerenkov[1]*G4UniformRand();

  for( iTransfer = 1; iTransfer <= fSplineNumber; iTransfer++ )
  {
        if( position >= fIntegralCerenkov[iTransfer] ) break;
  }
  if(iTransfer > fSplineNumber) iTransfer--;
 
  energyTransfer = fSplineEnergy[iTransfer];

  if(iTransfer > 1)
  {
    energyTransfer -= (fSplineEnergy[iTransfer]-fSplineEnergy[iTransfer-1])*G4UniformRand();
  }
  return energyTransfer;
}

GetDifPAIxSection()

G4double G4PAIxSection::GetDifPAIxSection ( G4int  i )

inline

Definition at line 193 of file G4PAIxSection.hh.

{ return fDifPAIxSection[i]; } 

GetElectronRange()

G4double G4PAIxSection::GetElectronRange

(

G4double

energy

)

Definition at line 1129 of file G4PAIxSection.cc.

{
  G4double range;
  /*
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();

  G4double Z = (*theMaterialTable)[fMaterialIndex]->GetIonisation()->GetZeffective();
  G4double A = (*theMaterialTable)[fMaterialIndex]->GetA();

  energy /= keV; // energy in keV in parametrised formula

  if (energy < 10.)
  {
    range = 3.872e-3*A/Z;
    range *= pow(energy,1.492);
  }
  else
  {
    range = 6.97e-3*pow(energy,1.6);
  }
  */
  // Blum&Rolandi Particle Detection with Drift Chambers, p. 7

  G4double cofA = 5.37e-4*g/cm2/keV;
  G4double cofB = 0.9815;
  G4double cofC = 3.123e-3/keV;
  // energy /= keV;

  range = cofA*energy*( 1 - cofB/(1 + cofC*energy) ); 

  // range *= g/cm2;
  range /= fDensity;

  return range;
}

Here is the call graph for this function:

GetEnergyInterval()

G4double G4PAIxSection::GetEnergyInterval ( G4int  i )

inline

Definition at line 191 of file G4PAIxSection.hh.

{ return fEnergyInterval[i]; } 

GetEnergyTransfer()

G4double G4PAIxSection::GetEnergyTransfer

(

)

Definition at line 2254 of file G4PAIxSection.cc.

{  
  G4int iTransfer ;

  G4double energyTransfer, position;

  position = fIntegralPAIxSection[1]*G4UniformRand();

  for( iTransfer = 1; iTransfer <= fSplineNumber; iTransfer++ )
  {
        if( position >= fIntegralPAIxSection[iTransfer] ) break;
  }
  if(iTransfer > fSplineNumber) iTransfer--;
 
  energyTransfer = fSplineEnergy[iTransfer];

  if(iTransfer > 1)
  {
    energyTransfer -= (fSplineEnergy[iTransfer]-fSplineEnergy[iTransfer-1])*G4UniformRand();
  }
  return energyTransfer;
}

GetIntegralCerenkov()

G4double G4PAIxSection::GetIntegralCerenkov ( G4int  i ) const

inline

Definition at line 324 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetIntegralMM()

G4double G4PAIxSection::GetIntegralMM ( G4int  i ) const

inline

Definition at line 330 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetIntegralPAIdEdx()

G4double G4PAIxSection::GetIntegralPAIdEdx ( G4int  i ) const

inline

Definition at line 318 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetIntegralPAIxSection()

G4double G4PAIxSection::GetIntegralPAIxSection ( G4int  i ) const

inline

Definition at line 312 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetIntegralPlasmon()

G4double G4PAIxSection::GetIntegralPlasmon ( G4int  i ) const

inline

Definition at line 336 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetIntegralResonance()

G4double G4PAIxSection::GetIntegralResonance ( G4int  i ) const

inline

Definition at line 342 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetIntervalNumber()

G4int G4PAIxSection::GetIntervalNumber ( ) const

inline

Definition at line 189 of file G4PAIxSection.hh.

{ return fIntervalNumber; }

GetLorentzFactor()

G4double G4PAIxSection::GetLorentzFactor ( G4int  i ) const

inline

Definition at line 301 of file G4PAIxSection.hh.

{
   return fLorentzFactor[j];
}

Here is the caller graph for this function:

GetLowEnergyCof()

G4double G4PAIxSection::GetLowEnergyCof ( ) const

inline

Definition at line 207 of file G4PAIxSection.hh.

{ return fLowEnergyCof; }

GetMeanCerenkovLoss()

G4double G4PAIxSection::GetMeanCerenkovLoss ( ) const

inline

Definition at line 200 of file G4PAIxSection.hh.

{return fIntegralCerenkov[0]; }

GetMeanEnergyLoss()

G4double G4PAIxSection::GetMeanEnergyLoss ( ) const

inline

Definition at line 199 of file G4PAIxSection.hh.

{return fIntegralPAIxSection[0]; }

GetMeanMMLoss()

G4double G4PAIxSection::GetMeanMMLoss ( ) const

inline

Definition at line 201 of file G4PAIxSection.hh.

{return fIntegralMM[0]; }

GetMeanPlasmonLoss()

G4double G4PAIxSection::GetMeanPlasmonLoss ( ) const

inline

Definition at line 202 of file G4PAIxSection.hh.

{return fIntegralPlasmon[0]; }

GetMeanResonanceLoss()

G4double G4PAIxSection::GetMeanResonanceLoss ( ) const

inline

Definition at line 203 of file G4PAIxSection.hh.

{return fIntegralResonance[0]; }

GetMMEnergyTransfer()

G4double G4PAIxSection::GetMMEnergyTransfer

(

)

Definition at line 2362 of file G4PAIxSection.cc.

{  
  G4int iTransfer ;

  G4double energyTransfer, position;

  position = fIntegralMM[1]*G4UniformRand();

  for( iTransfer = 1; iTransfer <= fSplineNumber; iTransfer++ )
  {
        if( position >= fIntegralMM[iTransfer] ) break;
  }
  if(iTransfer > fSplineNumber) iTransfer--;
 
  energyTransfer = fSplineEnergy[iTransfer];

  if(iTransfer > 1)
  {
    energyTransfer -= (fSplineEnergy[iTransfer]-fSplineEnergy[iTransfer-1])*G4UniformRand();
  }
  return energyTransfer;
}

GetNormalizationCof()

G4double G4PAIxSection::GetNormalizationCof ( ) const

inline

Definition at line 205 of file G4PAIxSection.hh.

{ return fNormalizationCof; }

GetNumberOfGammas()

G4int G4PAIxSection::GetNumberOfGammas ( ) const

inline

Definition at line 185 of file G4PAIxSection.hh.

{ return fNumberOfGammas; }

GetPAIdNdxCerenkov()

G4double G4PAIxSection::GetPAIdNdxCerenkov ( G4int  i )

inline

Definition at line 194 of file G4PAIxSection.hh.

{ return fdNdxCerenkov[i]; } 

GetPAIdNdxMM()

G4double G4PAIxSection::GetPAIdNdxMM ( G4int  i )

inline

Definition at line 195 of file G4PAIxSection.hh.

{ return fdNdxMM[i]; } 

GetPAIdNdxPlasmon()

G4double G4PAIxSection::GetPAIdNdxPlasmon ( G4int  i )

inline

Definition at line 196 of file G4PAIxSection.hh.

{ return fdNdxPlasmon[i]; } 

GetPAIdNdxResonance()

G4double G4PAIxSection::GetPAIdNdxResonance ( G4int  i )

inline

Definition at line 197 of file G4PAIxSection.hh.

{ return fdNdxResonance[i]; } 

GetPAItable()

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

inline

Definition at line 296 of file G4PAIxSection.hh.

{
   return fPAItable[i][j];
}

Here is the caller graph for this function:

GetPhotonRange()

G4double G4PAIxSection::GetPhotonRange

(

G4double

energy

)

Definition at line 1096 of file G4PAIxSection.cc.

{
  G4int i;
  G4double energy2, energy3, energy4, result, lambda;

  energy2 = energy1*energy1;
  energy3 = energy2*energy1;
  energy4 = energy3*energy1;

  // G4double* SandiaCof = fSandia->GetSandiaCofForMaterialPAI(energy1);
  // result = SandiaCof[0]/energy1+SandiaCof[1]/energy2+SandiaCof[2]/energy3+SandiaCof[3]/energy4;
  // result *= fDensity;

  for( i = 1; i <= fIntervalNumber; i++ )
  {
     if( energy1 < fEnergyInterval[i]) break;
  }
  i--;
  if(i == 0) i = 1;

  result = fA1[i]/energy1+fA2[i]/energy2+fA3[i]/energy3+fA4[i]/energy4;  

  if( result > DBL_MIN ) lambda = 1./result;
  else                   lambda = DBL_MAX;
   
  return lambda;
}  

GetPlasmonEnergyTransfer()

G4double G4PAIxSection::GetPlasmonEnergyTransfer

(

)

Definition at line 2416 of file G4PAIxSection.cc.

{  
  G4int iTransfer ;

  G4double energyTransfer, position;

  position = fIntegralPlasmon[1]*G4UniformRand();

  for( iTransfer = 1; iTransfer <= fSplineNumber; iTransfer++ )
  {
        if( position >= fIntegralPlasmon[iTransfer] ) break;
  }
  if(iTransfer > fSplineNumber) iTransfer--;
 
  energyTransfer = fSplineEnergy[iTransfer];

  if(iTransfer > 1)
  {
    energyTransfer -= (fSplineEnergy[iTransfer]-fSplineEnergy[iTransfer-1])*G4UniformRand();
  }
  return energyTransfer;
}

GetResonanceEnergyTransfer()

G4double G4PAIxSection::GetResonanceEnergyTransfer

(

)

Definition at line 2471 of file G4PAIxSection.cc.

{  
  G4int iTransfer ;

  G4double energyTransfer, position;

  position = fIntegralResonance[1]*G4UniformRand();

  for( iTransfer = 1; iTransfer <= fSplineNumber; iTransfer++ )
  {
        if( position >= fIntegralResonance[iTransfer] ) break;
  }
  if(iTransfer > fSplineNumber) iTransfer--;
 
  energyTransfer = fSplineEnergy[iTransfer];

  if(iTransfer > 1)
  {
    energyTransfer -= (fSplineEnergy[iTransfer]-fSplineEnergy[iTransfer-1])*G4UniformRand();
  }
  return energyTransfer;
}

GetRutherfordEnergyTransfer()

G4double G4PAIxSection::GetRutherfordEnergyTransfer

(

)

Definition at line 2499 of file G4PAIxSection.cc.

{  
  G4int iTransfer ;

  G4double energyTransfer, position;

  position = (fIntegralPlasmon[1]-fIntegralResonance[1])*G4UniformRand();

  for( iTransfer = 1; iTransfer <= fSplineNumber; iTransfer++ )
  {
        if( position >= (fIntegralPlasmon[iTransfer]-fIntegralResonance[iTransfer]) ) break;
  }
  if(iTransfer > fSplineNumber) iTransfer--;
 
  energyTransfer = fSplineEnergy[iTransfer];

  if(iTransfer > 1)
  {
    energyTransfer -= (fSplineEnergy[iTransfer]-fSplineEnergy[iTransfer-1])*G4UniformRand();
  }
  return energyTransfer;
}

GetSplineEnergy()

G4double G4PAIxSection::GetSplineEnergy ( G4int  i ) const

inline

Definition at line 306 of file G4PAIxSection.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetSplineSize()

G4int G4PAIxSection::GetSplineSize ( ) const

inline

Definition at line 187 of file G4PAIxSection.hh.

{ return fSplineNumber; }

GetStepCerenkovLoss()

G4double G4PAIxSection::GetStepCerenkovLoss

(

G4double

step

)

Definition at line 2281 of file G4PAIxSection.cc.

{  
  G4long numOfCollisions;
  G4double meanNumber, loss = 0.0;

  // G4cout<<" G4PAIxSection::GetStepCerenkovLoss "<<G4endl;

  meanNumber = fIntegralCerenkov[1]*step;
  numOfCollisions = G4Poisson(meanNumber);

  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;

  while(numOfCollisions)
  {
    loss += GetCerenkovEnergyTransfer();
    numOfCollisions--;
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  }
  // G4cout<<"PAI Cerenkov loss = "<<loss/keV<<" keV"<<G4endl; 

  return loss;
}

Here is the call graph for this function:

GetStepEnergyLoss()

G4double G4PAIxSection::GetStepEnergyLoss

(

G4double

step

)

Definition at line 2227 of file G4PAIxSection.cc.

{  
  G4long numOfCollisions;
  G4double meanNumber, loss = 0.0;

  // G4cout<<" G4PAIxSection::GetStepEnergyLoss "<<G4endl;

  meanNumber = fIntegralPAIxSection[1]*step;
  numOfCollisions = G4Poisson(meanNumber);

  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;

  while(numOfCollisions)
  {
    loss += GetEnergyTransfer();
    numOfCollisions--;
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  }
  // G4cout<<"PAI energy loss = "<<loss/keV<<" keV"<<G4endl; 

  return loss;
}

Here is the call graph for this function:

GetStepMMLoss()

G4double G4PAIxSection::GetStepMMLoss

(

G4double

step

)

Definition at line 2308 of file G4PAIxSection.cc.

{  
  G4long numOfCollisions;
  G4double meanNumber, loss = 0.0;

  // G4cout<<" G4PAIxSection::GetStepMMLoss "<<G4endl;

  meanNumber = fIntegralMM[1]*step;
  numOfCollisions = G4Poisson(meanNumber);

  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;

  while(numOfCollisions)
  {
    loss += GetMMEnergyTransfer();
    numOfCollisions--;
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  }
  // G4cout<<"PAI MM-Cerenkov loss = "<<loss/keV<<" keV"<<G4endl; 

  return loss;
}

Here is the call graph for this function:

GetStepPlasmonLoss()

G4double G4PAIxSection::GetStepPlasmonLoss

(

G4double

step

)

Definition at line 2389 of file G4PAIxSection.cc.

{  
  G4long numOfCollisions;
  G4double  meanNumber, loss = 0.0;

  // G4cout<<" G4PAIxSection::GetStepPlasmonLoss "<<G4endl;

  meanNumber = fIntegralPlasmon[1]*step;
  numOfCollisions = G4Poisson(meanNumber);

  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;

  while(numOfCollisions)
  {
    loss += GetPlasmonEnergyTransfer();
    numOfCollisions--;
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  }
  // G4cout<<"PAI Plasmon loss = "<<loss/keV<<" keV"<<G4endl; 

  return loss;
}

Here is the call graph for this function:

GetStepResonanceLoss()

G4double G4PAIxSection::GetStepResonanceLoss

(

G4double

step

)

Definition at line 2443 of file G4PAIxSection.cc.

{  
  G4long numOfCollisions;
  G4double meanNumber, loss = 0.0;

  // G4cout<<" G4PAIxSection::GetStepCreLosnkovs "<<G4endl;

  meanNumber = fIntegralResonance[1]*step;
  numOfCollisions = G4Poisson(meanNumber);

  //   G4cout<<"numOfCollisions = "<<numOfCollisions<<G4endl;

  while(numOfCollisions)
  {
    loss += GetResonanceEnergyTransfer();
    numOfCollisions--;
    // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  }
  // G4cout<<"PAI resonance loss = "<<loss/keV<<" keV"<<G4endl; 

  return loss;
}

Here is the call graph for this function:

ImPartDielectricConst()

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

Definition at line 1076 of file G4PAIxSection.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 

Initialize()

void G4PAIxSection::Initialize	(	const G4Material *	material,
		G4double	maxEnergyTransfer,
		G4double	betaGammaSq,
		G4SandiaTable *	sandia
	)

Definition at line 606 of file G4PAIxSection.cc.

{
  if(fVerbose > 0)
  {
    G4cout<<G4endl;
    G4cout<<"G4PAIxSection::Initialize(...,G4SandiaTable* sandia)"<<G4endl;
    G4cout<<G4endl;
  }
  G4int i, j;

  fSandia          = sandia;
  fIntervalNumber  = sandia->GetMaxInterval();
  fDensity         = material->GetDensity();
  fElectronDensity = material->GetElectronDensity();

  // fIntervalNumber--;

  if( fVerbose > 0 )
  {
    G4cout<<"fDensity = "<<fDensity<<"\t"<<fElectronDensity<<"\t fIntervalNumber = "<<fIntervalNumber<<G4endl;
  }  
  fEnergyInterval = G4DataVector(fIntervalNumber+2,0.0);
  fA1             = G4DataVector(fIntervalNumber+2,0.0);
  fA2             = G4DataVector(fIntervalNumber+2,0.0);
  fA3             = G4DataVector(fIntervalNumber+2,0.0);
  fA4             = G4DataVector(fIntervalNumber+2,0.0);

  for( i = 1; i <= fIntervalNumber; i++ ) 
  {
    if ( sandia->GetSandiaMatTablePAI(i-1,0) < 1.*eV && sandia->GetLowerI1() == false) 
    { 
      fIntervalNumber--;
      continue;
    }
    if( ( sandia->GetSandiaMatTablePAI(i-1,0) >= maxEnergyTransfer ) || i >= fIntervalNumber ) 
    {
      fEnergyInterval[i] = maxEnergyTransfer;
      fIntervalNumber = i;
      break;
    }
    fEnergyInterval[i] = sandia->GetSandiaMatTablePAI(i-1,0);
    fA1[i]             = sandia->GetSandiaMatTablePAI(i-1,1);
    fA2[i]             = sandia->GetSandiaMatTablePAI(i-1,2);
    fA3[i]             = sandia->GetSandiaMatTablePAI(i-1,3);
    fA4[i]             = sandia->GetSandiaMatTablePAI(i-1,4);

      if( fVerbose > 0 ) 
      {
        G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
             <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
      }
  }
  if( fVerbose > 0 ) G4cout<<"last i = "<<i<<"; "<<"fIntervalNumber = "<<fIntervalNumber<<G4endl;   

  if( fEnergyInterval[fIntervalNumber] != maxEnergyTransfer )
  {
      fIntervalNumber++;
      fEnergyInterval[fIntervalNumber] = maxEnergyTransfer;
  }
  if( fVerbose > 0 )
  {  
    for( i = 1; i <= fIntervalNumber; i++ )
    {
      G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
        <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }
  }  
  if( fVerbose > 0 )    G4cout<<"Now checking, if two borders are too close together"<<G4endl;

  for( i = 1; i < fIntervalNumber; i++ )
  {
    if( fEnergyInterval[i+1]-fEnergyInterval[i] >
         1.5*fDelta*(fEnergyInterval[i+1]+fEnergyInterval[i]) && fEnergyInterval[i] > 0.) continue;
    else
    {
      if( fVerbose > 0 )  G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fEnergyInterval[i+1]/keV;

      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 )
  {
    for( i = 1; i <= fIntervalNumber; i++ )
    {
      G4cout<<i<<"\t"<<fEnergyInterval[i]/keV<<"\t"<<fA1[i]<<"\t"<<fA2[i]<<"\t"
        <<fA3[i]<<"\t"<<fA4[i]<<"\t"<<G4endl;
    }
  }
  // Preparation of fSplineEnergy array corresponding to min ionisation, G~4

  ComputeLowEnergyCof(material);
            
  G4double   betaGammaSqRef = 
    fLorentzFactor[fRefGammaNumber]*fLorentzFactor[fRefGammaNumber] - 1;

  NormShift(betaGammaSqRef);             
  SplainPAI(betaGammaSqRef);
      
  // Preparation of integral PAI cross section for input betaGammaSq
   
  for( i = 1; i <= fSplineNumber; i++ )
  {
     fDifPAIxSection[i] = DifPAIxSection(i,betaGammaSq);


     fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
     fdNdxMM[i]   = PAIdNdxMM(i,betaGammaSq);
     fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
     fdNdxResonance[i]  = PAIdNdxResonance(i,betaGammaSq);
  }
  IntegralPAIxSection();   
  IntegralCerenkov();
  IntegralMM();
  IntegralPlasmon();
  IntegralResonance();
   
  for( i = 1; i <= fSplineNumber; i++ )
  {
    if(fVerbose>0) G4cout<<i<<"; w = "<<fSplineEnergy[i]/keV<<" keV; dN/dx_>w = "<<fIntegralPAIxSection[i]<<" 1/mm"<<G4endl;
  }
}

Here is the call graph for this function:

InitPAI()

void G4PAIxSection::InitPAI

(

)

Definition at line 813 of file G4PAIxSection.cc.

{    
   G4int i;
   G4double betaGammaSq = fLorentzFactor[fRefGammaNumber]*
                          fLorentzFactor[fRefGammaNumber] - 1;

   // Preparation of integral PAI cross section for reference gamma
   
   NormShift(betaGammaSq);             
   SplainPAI(betaGammaSq);

   IntegralPAIxSection();
   IntegralCerenkov();
   IntegralMM();
   IntegralPlasmon();
   IntegralResonance();

   for(i = 0; i<= fSplineNumber; i++)
   {
      fPAItable[i][fRefGammaNumber] = fIntegralPAIxSection[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++)
      {
         fDifPAIxSection[i] = DifPAIxSection(i,betaGammaSq);
         fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
         fdNdxMM[i]   = PAIdNdxMM(i,betaGammaSq);
         fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
         fdNdxResonance[i]  = PAIdNdxResonance(i,betaGammaSq);
      }
      IntegralPAIxSection();
      IntegralCerenkov();
      IntegralMM();
      IntegralPlasmon();
      IntegralResonance();
      
      for(i = 0; i <= fSplineNumber; i++)
      {
         fPAItable[i][j] = fIntegralPAIxSection[i];
      }
   } 

}  

IntegralCerenkov()

void G4PAIxSection::IntegralCerenkov

(

)

Definition at line 1522 of file G4PAIxSection.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 

IntegralMM()

void G4PAIxSection::IntegralMM

(

)

Definition at line 1553 of file G4PAIxSection.cc.

{
  G4int i, k;
   fIntegralMM[fSplineNumber] = 0;
   fIntegralMM[0] = 0;
   k = fIntervalNumber -1;

   for( i = fSplineNumber-1; i >= 1; i-- )
   {
      if(fSplineEnergy[i] >= fEnergyInterval[k])
      {
        fIntegralMM[i] = fIntegralMM[i+1] + SumOverInterMM(i);
        // G4cout<<"int: i = "<<i<<"; sumC = "<<fIntegralMM[i]<<G4endl;
      }
      else
      {
        fIntegralMM[i] = fIntegralMM[i+1] + 
                                   SumOverBordMM(i+1,fEnergyInterval[k]);
        k--;
        // G4cout<<"bord: i = "<<i<<"; sumC = "<<fIntegralMM[i]<<G4endl;
      }
   }

}   // end of IntegralMM 

IntegralPAIxSection()

void G4PAIxSection::IntegralPAIxSection

(

)

Definition at line 1490 of file G4PAIxSection.cc.

{
  fIntegralPAIxSection[fSplineNumber] = 0;
  fIntegralPAIdEdx[fSplineNumber]     = 0;
  fIntegralPAIxSection[0]             = 0;
  G4int i, k = fIntervalNumber -1;

  for( i = fSplineNumber-1; i >= 1; i--)
  {
    if(fSplineEnergy[i] >= fEnergyInterval[k])
    {
      fIntegralPAIxSection[i] = fIntegralPAIxSection[i+1] + SumOverInterval(i);
      fIntegralPAIdEdx[i] = fIntegralPAIdEdx[i+1] + SumOverIntervaldEdx(i);
    }
    else
    {
      fIntegralPAIxSection[i] = fIntegralPAIxSection[i+1] + 
                                   SumOverBorder(i+1,fEnergyInterval[k]);
      fIntegralPAIdEdx[i] = fIntegralPAIdEdx[i+1] + 
                                   SumOverBorderdEdx(i+1,fEnergyInterval[k]);
      k--;
    }
    if(fVerbose>0) G4cout<<"i = "<<i<<"; k = "<<k<<"; intPAIxsc[i] = "<<fIntegralPAIxSection[i]<<G4endl;
  }
}   // end of IntegralPAIxSection 

IntegralPlasmon()

void G4PAIxSection::IntegralPlasmon

(

)

Definition at line 1584 of file G4PAIxSection.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

IntegralResonance()

void G4PAIxSection::IntegralResonance

(

)

Definition at line 1611 of file G4PAIxSection.cc.

{
   fIntegralResonance[fSplineNumber] = 0;
   fIntegralResonance[0] = 0;
   G4int k = fIntervalNumber -1;
   for(G4int i=fSplineNumber-1;i>=1;i--)
   {
      if(fSplineEnergy[i] >= fEnergyInterval[k])
      {
        fIntegralResonance[i] = fIntegralResonance[i+1] + SumOverInterResonance(i);
      }
      else
      {
        fIntegralResonance[i] = fIntegralResonance[i+1] + 
                                   SumOverBordResonance(i+1,fEnergyInterval[k]);
        k--;
      }
   }

}   // end of IntegralResonance

NormShift()

void G4PAIxSection::NormShift

(

G4double

betaGammaSq

)

Definition at line 873 of file G4PAIxSection.cc.

{
  G4int i, j;

  if(fVerbose>0) G4cout<<"      G4PAIxSection::NormShift call "<<G4endl;


  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); 
      if(fVerbose>0) G4cout<<"cn = "<<fSplineNumber<<"; "<<"w = "<<fSplineEnergy[fSplineNumber]/keV<<" keV"<<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]    );
    }
   if(fVerbose>0)  G4cout<<i<<"  Shift: w = "<<fSplineEnergy[i]/keV<<" keV \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;

         fDifPAIxSection[i] = DifPAIxSection(i,betaGammaSq);
         fdNdxCerenkov[i]   = PAIdNdxCerenkov(i,betaGammaSq);
         fdNdxMM[i]   = PAIdNdxMM(i,betaGammaSq);
         fdNdxPlasmon[i]    = PAIdNdxPlasmon(i,betaGammaSq);
         fdNdxResonance[i]    = PAIdNdxResonance(i,betaGammaSq);
   if(fVerbose>0)  G4cout<<i<<"  Shift: w = "<<fSplineEnergy[i]/keV<<" keV, xsc = "<<fDifPAIxSection[i]<<"\n"<<G4endl;
      }
   }

}  // end of NormShift 

operator=()

G4PAIxSection& G4PAIxSection::operator= ( const G4PAIxSection &  right )

private

Here is the caller graph for this function:

PAIdNdxCerenkov()

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

Definition at line 1302 of file G4PAIxSection.cc.

{        
   G4double logarithm, x3, x5, argument, modul2, dNdxC; 
   G4double be2, betaBohr2, cofBetaBohr;

   cofBetaBohr = 4.0;
   betaBohr2   = fine_structure_const*fine_structure_const;
   G4double betaBohr4   = betaBohr2*betaBohr2*cofBetaBohr;

   be2 = betaGammaSq/(1 + betaGammaSq);
   G4double 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];
      dNdxC /= modul2;
   }
   return dNdxC;

} // end of PAIdNdxCerenkov 

PAIdNdxMM()

G4double G4PAIxSection::PAIdNdxMM	(	G4int	intervalNumber,
		G4double	betaGammaSq
	)

Definition at line 1360 of file G4PAIxSection.cc.

{        
   G4double logarithm, x3, x5, argument, dNdxC; 
   G4double be2, be4, betaBohr2,betaBohr4,cofBetaBohr;

   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 = be2*( 1.0 + fRePartDielectricConst[i] ) - 1.0;
     if( x3 == 0.0 ) argument = 0.5*pi;
     else            argument = atan2(fImPartDielectricConst[i],x3);
     argument *= x5 ;
   }   
   dNdxC = ( logarithm*fImPartDielectricConst[i]*be2 + argument )/hbarc;
  
   if(dNdxC < 1.0e-8) dNdxC = 1.0e-8;

   dNdxC *= fine_structure_const/be2/pi;

   dNdxC *= (1-exp(-be4/betaBohr4));
   return dNdxC;

} // end of PAIdNdxMM 

PAIdNdxPlasmon()

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

Definition at line 1410 of file G4PAIxSection.cc.

{        
   G4double resonance, modul2, dNdxP, cof = 1.;
   G4double be2, betaBohr;
  
   betaBohr   = fine_structure_const;
   be2 = betaGammaSq/(1 + betaGammaSq);

   G4double beta = sqrt(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(-beta/betaBohr/fLowEnergyCof));

   // dNdxP *= (1-exp(-be4/betaBohr4));

   if( fDensity >= 0.1 )
   { 
     modul2 = (1 + fRePartDielectricConst[i])*(1 + fRePartDielectricConst[i]) + 
        fImPartDielectricConst[i]*fImPartDielectricConst[i];
     dNdxP /= modul2;
   }
   return dNdxP;

} // end of PAIdNdxPlasmon 

PAIdNdxResonance()

G4double G4PAIxSection::PAIdNdxResonance	(	G4int	intervalNumber,
		G4double	betaGammaSq
	)

Definition at line 1450 of file G4PAIxSection.cc.

{        
   G4double resonance, modul2, dNdxP;
   G4double be2, be4, betaBohr2, betaBohr4, cofBetaBohr;

   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;

   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];
     dNdxP /= modul2;
   }
   return dNdxP;

} // end of PAIdNdxResonance 

RePartDielectricConst()

G4double G4PAIxSection::RePartDielectricConst

(

G4double

energy

)

Definition at line 1171 of file G4PAIxSection.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 

RutherfordIntegral()

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

Definition at line 1055 of file G4PAIxSection.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 

SetVerbose()

void G4PAIxSection::SetVerbose ( G4int  v )

inline

Definition at line 209 of file G4PAIxSection.hh.

{fVerbose=v;};

Here is the call graph for this function:

SplainPAI()

void G4PAIxSection::SplainPAI

(

G4double

betaGammaSq

)

Definition at line 950 of file G4PAIxSection.cc.

{
  G4int j, k = 1, i = 1;

  if(fVerbose>0) G4cout<<"                   G4PAIxSection::SplainPAI call "<<G4endl;

  while ( (i < fSplineNumber) && (fSplineNumber < fMaxSplineSize-1) )
  {
     // if( std::abs(fSplineEnergy[i+1]-fEnergyInterval[k+1]) > (fSplineEnergy[i+1]+fEnergyInterval[k+1])*5.e-7 )
     if( fSplineEnergy[i+1] > fEnergyInterval[k+1] )
     {
          k++;   // Here next energy point is in next energy interval
          i++;
          if(fVerbose>0) G4cout<<"                     in if: i = "<<i<<"; k = "<<k<<G4endl;
          continue;
     }
     if(fVerbose>0) G4cout<<"       out if: i = "<<i<<"; k = "<<k<<G4endl;

                        // Shifting of arrayes for inserting the geometrical 
                       // average of 'i' and 'i+1' energy points to 'i+1' place
     fSplineNumber++;

     for( 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];

         fDifPAIxSection[j] = fDifPAIxSection[j-1];
         fdNdxCerenkov[j]   = fdNdxCerenkov[j-1];
         fdNdxMM[j]         = fdNdxMM[j-1];
         fdNdxPlasmon[j]    = fdNdxPlasmon[j-1];
         fdNdxResonance[j]  = fdNdxResonance[j-1];
     }
      G4double x1  = fSplineEnergy[i];
      G4double x2  = fSplineEnergy[i+1];
      G4double yy1 = fDifPAIxSection[i];
      G4double y2  = fDifPAIxSection[i+1];

      if(fVerbose>0) G4cout<<"Spline: x1 = "<<x1<<"; x2 = "<<x2<<", yy1 = "<<yy1<<"; y2 = "<<y2<<G4endl;


      G4double en1 = sqrt(x1*x2);
      // G4double    en1 = 0.5*(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]);

      fDifPAIxSection[i+1] = DifPAIxSection(i+1,betaGammaSq);
      fdNdxCerenkov[i+1]   = PAIdNdxCerenkov(i+1,betaGammaSq);
      fdNdxMM[i+1]         = PAIdNdxMM(i+1,betaGammaSq);
      fdNdxPlasmon[i+1]    = PAIdNdxPlasmon(i+1,betaGammaSq);
      fdNdxResonance[i+1]  = PAIdNdxResonance(i+1,betaGammaSq);

                  // Condition for next division of this segment or to pass

    if(fVerbose>0) G4cout<<"Spline, a = "<<a<<"; b = "<<b<<"; new xsc = "<<y<<"; compxsc = "<<fDifPAIxSection[i+1]<<G4endl;

                  // to higher energies

      G4double x = 2*(fDifPAIxSection[i+1] - y)/(fDifPAIxSection[i+1] + y);

      G4double delta = 2.*(fSplineEnergy[i+1]-fSplineEnergy[i])/(fSplineEnergy[i+1]+fSplineEnergy[i]);

      if( x < 0 ) 
      {
         x = -x;
      }
      if( x > fError && fSplineNumber < fMaxSplineSize-1 && delta > 2.*fDelta )
      {
         continue;  // next division
      }
      i += 2;  // pass to next segment

      // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
  }   // close 'while'

}  // end of SplainPAI 

SumOverBordCerenkov()

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

Definition at line 1990 of file G4PAIxSection.cc.

{               
   G4double x0,x1,y0,yy1,a,b,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);
   if(a > 10.0) return 0.;  
   // b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a); // pow(10.,b0);   
   
   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<<"; b0 = "<<b0<<"; 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);
   // b0 = log10(y0) - a*log10(x0);
   b  =  y0/pow(x0,a);  // pow(10.,b0);

   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 )   fIntegralCerenkov[0] += b*log(e0/x0);
   else           fIntegralCerenkov[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;

   // G4cout<<"a = "<<a<<"; b0 = "<<b0<<"; b = "
   // <<b<<"; result = "<<result<<G4endl;    

   return result;

} 

SumOverBorder()

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

Definition at line 1862 of file G4PAIxSection.cc.

{               
  G4double x0,x1,y0,yy1,a,b,/*c,*/d,e0,result;

   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fDifPAIxSection[i];
   yy1 = fDifPAIxSection[i+1];

   //c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(x1/x0);
   if(a > 10.0) return 0.;  

   if(fVerbose>0) G4cout<<"SumOverBorder, a = "<<a<<G4endl;

   // b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);  // pow(10.,b);
   
   a += 1.;
   if( std::fabs(a) < 1.e-6 )
   {
      result = b*log(x0/e0);
   }
   else
   {
      result = y0*(x0 - e0*pow(d,a-1))/a;
   }
   a += 1.;
   if( std::fabs(a) < 1.e-6 )
   {
      fIntegralPAIxSection[0] += b*log(x0/e0);
   }
   else 
   {
      fIntegralPAIxSection[0] += y0*(x0*x0 - e0*e0*pow(d,a-2))/a;
   }
   x0 = fSplineEnergy[i - 1];
   x1 = fSplineEnergy[i - 2];
   y0 = fDifPAIxSection[i - 1];
   yy1 = fDifPAIxSection[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( std::fabs(a) < 1.e-6 )
   {
      result += b*log(e0/x0);
   }
   else
   {
      result += y0*(e0*pow(d,a-1) - x0)/a;
   }
   a += 1.;
   if( std::fabs(a) < 1.e-6 ) 
   {
      fIntegralPAIxSection[0] += b*log(e0/x0);
   }
   else
   {
      fIntegralPAIxSection[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
   }
   return result;

} 

SumOverBorderdEdx()

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

Definition at line 1935 of file G4PAIxSection.cc.

{               
  G4double x0,x1,y0,yy1,a,b,/*c,*/d,e0,result;

   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fDifPAIxSection[i];
   yy1 = fDifPAIxSection[i+1];

   //c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(x1/x0);
   if(a > 10.0) return 0.;  
   // b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);  // pow(10.,b);
   
   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 = fDifPAIxSection[i - 1];
   yy1 = fDifPAIxSection[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 += 2;
   if(a == 0) 
   {
      result += b*log(e0/x0);
   }
   else
   {
      result += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
   }
   return result;

} 

SumOverBordMM()

G4double G4PAIxSection::SumOverBordMM	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 2055 of file G4PAIxSection.cc.

{               
   G4double x0,x1,y0,yy1,a,b,e0,c,d,result;

   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fdNdxMM[i];
   yy1 = fdNdxMM[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);
   if(a > 10.0) return 0.;  
   // b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a); // pow(10.,b0);   
   
   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 ) fIntegralMM[0] += b*log(x0/e0);
   else         fIntegralMM[0] += y0*(x0*x0 - e0*e0*pow(d,a-2))/a;

// G4cout<<"a = "<<a<<"; b0 = "<<b0<<"; b = "<<b<<"; result = "<<result<<G4endl;
   
   x0  = fSplineEnergy[i - 1];
   x1  = fSplineEnergy[i - 2];
   y0  = fdNdxMM[i - 1];
   yy1 = fdNdxMM[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);
   // b0 = log10(y0) - a*log10(x0);
   b  =  y0/pow(x0,a);  // pow(10.,b0);

   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 )   fIntegralMM[0] += b*log(e0/x0);
   else           fIntegralMM[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;

   // G4cout<<"a = "<<a<<"; b0 = "<<b0<<"; b = "
   // <<b<<"; result = "<<result<<G4endl;    

   return result;

} 

SumOverBordPlasmon()

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

Definition at line 2120 of file G4PAIxSection.cc.

{               
   G4double x0,x1,y0,yy1,a,b,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);
   if(a > 10.0) return 0.;  
   //  b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a); //pow(10.,b);
   
   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);
   // b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);// pow(10.,b0);

   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;

} 

SumOverBordResonance()

G4double G4PAIxSection::SumOverBordResonance	(	G4int	intervalNumber,
		G4double	energy
	)

Definition at line 2174 of file G4PAIxSection.cc.

{               
   G4double x0,x1,y0,yy1,a,b,c,d,e0,result;

   e0 = en0;
   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   y0 = fdNdxResonance[i];
   yy1 = fdNdxResonance[i+1];

   c = x1/x0;
   d = e0/x0;   
   a = log10(yy1/y0)/log10(c);
   if(a > 10.0) return 0.;  
   //  b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a); //pow(10.,b);
   
   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 ) fIntegralResonance[0] += b*log(x0/e0);
   else         fIntegralResonance[0] += y0*(x0*x0 - e0*e0*pow(d,a-2))/a;
   
   x0 = fSplineEnergy[i - 1];
   x1 = fSplineEnergy[i - 2];
   y0 = fdNdxResonance[i - 1];
   yy1 = fdNdxResonance[i - 2];

   c = x1/x0;
   d = e0/x0;
   a = log10(yy1/y0)/log10(c);
   // b0 = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);// pow(10.,b0);

   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 )   fIntegralResonance[0] += b*log(e0/x0);
   else           fIntegralResonance[0] += y0*(e0*e0*pow(d,a-2) - x0*x0)/a;
   
   return result;

} 

SumOverInterCerenkov()

G4double G4PAIxSection::SumOverInterCerenkov

(

G4int

intervalNumber

)

Definition at line 1719 of file G4PAIxSection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;

   x0  = fSplineEnergy[i];
   x1  = fSplineEnergy[i+1];

   if(x1+x0 <= 0.0 || std::abs( 2.*(x1-x0)/(x1+x0) ) < 1.e-6) return 0.;

   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);
   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

SumOverInterMM()

G4double G4PAIxSection::SumOverInterMM

(

G4int

intervalNumber

)

Definition at line 1755 of file G4PAIxSection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;

   x0  = fSplineEnergy[i];
   x1  = fSplineEnergy[i+1];

   if(x1+x0 <= 0.0 || std::abs( 2.*(x1-x0)/(x1+x0) ) < 1.e-6) return 0.;

   y0  = fdNdxMM[i];
   yy1 = fdNdxMM[i+1];
   //G4cout<<"SumC, i = "<<i<<"; x0 ="<<x0<<"; x1 = "<<x1
   //   <<"; y0 = "<<y0<<"; yy1 = "<<yy1<<G4endl;

   c = x1/x0;
   //G4cout<<" c = "<<c<< " yy1/y0= " << yy1/y0 <<G4endl;   
   a = log10(yy1/y0)/log10(c);
   if(a > 10.0) return 0.;  
   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 ) fIntegralMM[0] += b*log(c);
   else         fIntegralMM[0] += y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   //G4cout<<"a = "<<a<<"; b = "<<b<<"; result = "<<result<<G4endl;   
   return result;

} //  end of SumOverInterMM

SumOverInterPlasmon()

G4double G4PAIxSection::SumOverInterPlasmon

(

G4int

intervalNumber

)

Definition at line 1793 of file G4PAIxSection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;

   x0  = fSplineEnergy[i];
   x1  = fSplineEnergy[i+1];

   if(x1+x0 <= 0.0 || std::abs( 2.*(x1-x0)/(x1+x0) ) < 1.e-6) return 0.;

   y0  = fdNdxPlasmon[i];
   yy1 = fdNdxPlasmon[i+1];
   c =x1/x0;
   a = log10(yy1/y0)/log10(c);
   if(a > 10.0) return 0.;  
   // b = log10(y0) - a*log10(x0);
   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

SumOverInterResonance()

G4double G4PAIxSection::SumOverInterResonance

(

G4int

intervalNumber

)

Definition at line 1828 of file G4PAIxSection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;

   x0  = fSplineEnergy[i];
   x1  = fSplineEnergy[i+1];

   if(x1+x0 <= 0.0 || std::abs( 2.*(x1-x0)/(x1+x0) ) < 1.e-6) return 0.;

   y0  = fdNdxResonance[i];
   yy1 = fdNdxResonance[i+1];
   c =x1/x0;
   a = log10(yy1/y0)/log10(c);
   if(a > 10.0) return 0.;  
   // b = log10(y0) - a*log10(x0);
   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 ) fIntegralResonance[0] += b*log(x1/x0);
   else         fIntegralResonance[0] += y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   
   return result;

} //  end of SumOverInterResonance

SumOverInterval()

G4double G4PAIxSection::SumOverInterval

(

G4int

intervalNumber

)

Definition at line 1638 of file G4PAIxSection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;

   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];
   if(fVerbose>0) G4cout<<"SumOverInterval i= " << i << " x0 = "<<x0<<"; x1 = "<<x1<<G4endl;

   if( x1+x0 <= 0.0 || std::abs( 2.*(x1-x0)/(x1+x0) ) < 1.e-6) return 0.;

   y0 = fDifPAIxSection[i];
   yy1 = fDifPAIxSection[i+1];

   if(fVerbose>0) G4cout<<"x0 = "<<x0<<"; x1 = "<<x1<<", y0 = "<<y0<<"; yy1 = "<<yy1<<G4endl;

   c = x1/x0;
   a = log10(yy1/y0)/log10(c);

   if(fVerbose>0) G4cout<<"SumOverInterval, a = "<<a<<"; c = "<<c<<G4endl;

   // b = log10(y0) - a*log10(x0);
   b = y0/pow(x0,a);
   a += 1.;
   if( std::fabs(a) < 1.e-6 ) 
   {
      result = b*log(x1/x0);
   }
   else
   {
      result = y0*(x1*pow(c,a-1) - x0)/a;
   }
   a += 1.;
   if( std::fabs(a) < 1.e-6 ) 
   {
      fIntegralPAIxSection[0] += b*log(x1/x0);
   }
   else
   {
      fIntegralPAIxSection[0] += y0*(x1*x1*pow(c,a-2) - x0*x0)/a;
   }
   if(fVerbose>0) G4cout<<"SumOverInterval, result = "<<result<<G4endl;
   return result;

} //  end of SumOverInterval

SumOverIntervaldEdx()

G4double G4PAIxSection::SumOverIntervaldEdx

(

G4int

intervalNumber

)

Definition at line 1685 of file G4PAIxSection.cc.

{         
   G4double x0,x1,y0,yy1,a,b,c,result;

   x0 = fSplineEnergy[i];
   x1 = fSplineEnergy[i+1];

   if(x1+x0 <= 0.0 || std::abs( 2.*(x1-x0)/(x1+x0) ) < 1.e-6) return 0.;

   y0 = fDifPAIxSection[i];
   yy1 = fDifPAIxSection[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

Member Data Documentation

fA1

G4DataVector G4PAIxSection::fA1

private

Definition at line 261 of file G4PAIxSection.hh.

fA2

G4DataVector G4PAIxSection::fA2

private

Definition at line 262 of file G4PAIxSection.hh.

fA3

G4DataVector G4PAIxSection::fA3

private

Definition at line 263 of file G4PAIxSection.hh.

fA4

G4DataVector G4PAIxSection::fA4

private

Definition at line 264 of file G4PAIxSection.hh.

fDelta

const G4double G4PAIxSection::fDelta = 0.005

staticprivate

Definition at line 233 of file G4PAIxSection.hh.

fDensity

G4double G4PAIxSection::fDensity

private

Definition at line 248 of file G4PAIxSection.hh.

fDifPAIxSection

G4DataVector G4PAIxSection::fDifPAIxSection

private

Definition at line 274 of file G4PAIxSection.hh.

fdNdxCerenkov

G4DataVector G4PAIxSection::fdNdxCerenkov

private

Definition at line 275 of file G4PAIxSection.hh.

fdNdxMM

G4DataVector G4PAIxSection::fdNdxMM

private

Definition at line 277 of file G4PAIxSection.hh.

fdNdxPlasmon

G4DataVector G4PAIxSection::fdNdxPlasmon

private

Definition at line 276 of file G4PAIxSection.hh.

fdNdxResonance

G4DataVector G4PAIxSection::fdNdxResonance

private

Definition at line 278 of file G4PAIxSection.hh.

fElectronDensity

G4double G4PAIxSection::fElectronDensity

private

Definition at line 249 of file G4PAIxSection.hh.

fEnergyInterval

G4DataVector G4PAIxSection::fEnergyInterval

private

Definition at line 260 of file G4PAIxSection.hh.

fError

const G4double G4PAIxSection::fError = 0.005

staticprivate

Definition at line 234 of file G4PAIxSection.hh.

fImPartDielectricConst

G4DataVector G4PAIxSection::fImPartDielectricConst

private

Definition at line 272 of file G4PAIxSection.hh.

fIntegralCerenkov

G4DataVector G4PAIxSection::fIntegralCerenkov

private

Definition at line 282 of file G4PAIxSection.hh.

fIntegralMM

G4DataVector G4PAIxSection::fIntegralMM

private

Definition at line 284 of file G4PAIxSection.hh.

fIntegralPAIdEdx

G4DataVector G4PAIxSection::fIntegralPAIdEdx

private

Definition at line 281 of file G4PAIxSection.hh.

fIntegralPAIxSection

G4DataVector G4PAIxSection::fIntegralPAIxSection

private

Definition at line 280 of file G4PAIxSection.hh.

fIntegralPlasmon

G4DataVector G4PAIxSection::fIntegralPlasmon

private

Definition at line 283 of file G4PAIxSection.hh.

fIntegralResonance

G4DataVector G4PAIxSection::fIntegralResonance

private

Definition at line 285 of file G4PAIxSection.hh.

fIntegralTerm

G4DataVector G4PAIxSection::fIntegralTerm

private

Definition at line 273 of file G4PAIxSection.hh.

fIntervalNumber

G4int G4PAIxSection::fIntervalNumber

private

Definition at line 242 of file G4PAIxSection.hh.

fLorentzFactor

const G4double G4PAIxSection::fLorentzFactor

staticprivate

Initial value:

=     
{
0.0,
1.094989e+00, 1.107813e+00, 1.122369e+00, 1.138890e+00, 1.157642e+00,
1.178925e+00, 1.203082e+00, 1.230500e+00, 1.261620e+00, 1.296942e+00, 
1.337032e+00, 1.382535e+00, 1.434181e+00, 1.492800e+00, 1.559334e+00,
1.634850e+00, 1.720562e+00, 1.817845e+00, 1.928263e+00, 2.053589e+00, 
2.195835e+00, 2.357285e+00, 2.540533e+00, 2.748522e+00, 2.984591e+00,
3.252533e+00, 3.556649e+00, 3.901824e+00, 4.293602e+00, 4.738274e+00, 
5.242981e+00, 5.815829e+00, 6.466019e+00, 7.203990e+00, 8.041596e+00,
8.992288e+00, 1.007133e+01, 1.129606e+01, 1.268614e+01, 1.426390e+01, 
1.605467e+01, 1.808721e+01, 2.039417e+01, 2.301259e+01, 2.598453e+01,
2.935771e+01, 3.318630e+01, 3.753180e+01, 4.246399e+01, 4.806208e+01, 
5.441597e+01, 6.162770e+01, 6.981310e+01, 7.910361e+01, 8.964844e+01,
1.016169e+02, 1.152013e+02, 1.306197e+02, 1.481198e+02, 1.679826e+02, 
1.905270e+02, 2.161152e+02, 2.451581e+02, 2.781221e+02, 3.155365e+02,
3.580024e+02, 4.062016e+02, 4.609081e+02, 5.230007e+02, 5.934765e+02, 
6.734672e+02, 7.642575e+02, 8.673056e+02, 9.842662e+02, 1.117018e+03,
1.267692e+03, 1.438709e+03, 1.632816e+03, 1.853128e+03, 2.103186e+03, 
2.387004e+03, 2.709140e+03, 3.074768e+03, 3.489760e+03, 3.960780e+03,
4.495394e+03, 5.102185e+03, 5.790900e+03, 6.572600e+03, 7.459837e+03, 
8.466860e+03, 9.609843e+03, 1.090714e+04, 1.237959e+04, 1.405083e+04,
1.594771e+04, 1.810069e+04, 2.054434e+04, 2.331792e+04, 2.646595e+04, 
3.003901e+04, 3.409446e+04, 3.869745e+04, 4.392189e+04, 4.985168e+04,
5.658206e+04, 6.422112e+04, 7.289153e+04, 8.273254e+04, 9.390219e+04, 
1.065799e+05
}

Definition at line 237 of file G4PAIxSection.hh.

fLowEnergyCof

G4double G4PAIxSection::fLowEnergyCof

private

Definition at line 250 of file G4PAIxSection.hh.

fMaterialIndex

G4int G4PAIxSection::fMaterialIndex

private

Definition at line 247 of file G4PAIxSection.hh.

fMatSandiaMatrix

G4OrderedTable* G4PAIxSection::fMatSandiaMatrix

private

Definition at line 256 of file G4PAIxSection.hh.

fMaxSplineSize

const G4int G4PAIxSection::fMaxSplineSize = 1000

staticprivate

Definition at line 267 of file G4PAIxSection.hh.

fNormalizationCof

G4double G4PAIxSection::fNormalizationCof

private

Definition at line 243 of file G4PAIxSection.hh.

fNumberOfGammas

G4int G4PAIxSection::fNumberOfGammas = 111

staticprivate

Definition at line 236 of file G4PAIxSection.hh.

fPAItable

G4double G4PAIxSection::fPAItable[500][112]

private

Definition at line 289 of file G4PAIxSection.hh.

fRefGammaNumber

const G4int G4PAIxSection::fRefGammaNumber = 29

staticprivate

Definition at line 240 of file G4PAIxSection.hh.

fRePartDielectricConst

G4DataVector G4PAIxSection::fRePartDielectricConst

private

Definition at line 271 of file G4PAIxSection.hh.

fSandia

G4SandiaTable* G4PAIxSection::fSandia

private

Definition at line 258 of file G4PAIxSection.hh.

fSplineEnergy

G4DataVector G4PAIxSection::fSplineEnergy

private

Definition at line 270 of file G4PAIxSection.hh.

fSplineNumber

G4int G4PAIxSection::fSplineNumber

private

Definition at line 251 of file G4PAIxSection.hh.

fVerbose

G4int G4PAIxSection::fVerbose

private

Definition at line 252 of file G4PAIxSection.hh.

---

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

• G4PAIxSection.hh
• G4PAIxSection.cc