#include <G4GSPWACorrections.hh>

Public Member Functions
	G4GSPWACorrections (G4bool iselectron=true)

	~G4GSPWACorrections ()

void	Initialise ()

void	GetPWACorrectionFactors (G4double logekin, G4double beta2, G4int matindx, G4double &corToScr, G4double &corToQ1, G4double &corToG2PerG1)

Detailed Description

Definition at line 61 of file G4GSPWACorrections.hh.

Constructor & Destructor Documentation

◆ G4GSPWACorrections()

G4GSPWACorrections::G4GSPWACorrections ( G4bool iselectron = true )

Definition at line 65 of file G4GSPWACorrections.cc.

                                                        : fIsElectron(iselectron) {
  // init grids related data member values
  fMaxEkin        = CLHEP::electron_mass_c2*(1./std::sqrt(1.-gMaxBeta2)-1.);
  fLogMinEkin     = G4Log(gMinEkin);
  fInvLogDelEkin  = (gNumEkin-gNumBeta2)/G4Log(gMidEkin/gMinEkin);
  G4double pt2    = gMidEkin*(gMidEkin+2.0*CLHEP::electron_mass_c2);
  fMinBeta2       = pt2/(pt2+CLHEP::electron_mass_c2*CLHEP::electron_mass_c2);
  fInvDelBeta2    = (gNumBeta2-1.)/(gMaxBeta2-fMinBeta2);
}

◆ ~G4GSPWACorrections()

G4GSPWACorrections::~G4GSPWACorrections ( )

Definition at line 76 of file G4GSPWACorrections.cc.

                                        {
  ClearDataPerElement();
  ClearDataPerMaterial();
}

Member Function Documentation

◆ GetPWACorrectionFactors()

void G4GSPWACorrections::GetPWACorrectionFactors	(	G4double	logekin,
		G4double	beta2,
		G4int	matindx,
		G4double &	corToScr,
		G4double &	corToQ1,
		G4double &	corToG2PerG1 )

Definition at line 82 of file G4GSPWACorrections.cc.

                                                                                                                 {
  G4int    ekinIndxLow = 0;
  G4double remRfaction = 0.;
  if (beta2>=gMaxBeta2) {
    ekinIndxLow = gNumEkin - 1;
    // remRfaction = -1.
  } else if (beta2>=fMinBeta2) {  // linear interpolation on \beta^2
    remRfaction   = (beta2 - fMinBeta2) * fInvDelBeta2;
    ekinIndxLow   = (G4int)remRfaction;
    remRfaction  -= ekinIndxLow;
    ekinIndxLow  += (gNumEkin - gNumBeta2);
  } else if (logekin>=fLogMinEkin) {
    remRfaction   = (logekin - fLogMinEkin) * fInvLogDelEkin;
    ekinIndxLow   = (G4int)remRfaction;
    remRfaction  -= ekinIndxLow;
  } // the defaults otherwise i.e. use the lowest energy values when ekin is smaller than the minum ekin
  //
  DataPerMaterial *data = fDataPerMaterial[matindx];
  corToScr      = data->fCorScreening[ekinIndxLow];
  corToQ1       = data->fCorFirstMoment[ekinIndxLow];
  corToG2PerG1  = data->fCorSecondMoment[ekinIndxLow];
  if (remRfaction>0.) {
    corToScr      += remRfaction*(data->fCorScreening[ekinIndxLow+1]    - data->fCorScreening[ekinIndxLow]);
    corToQ1       += remRfaction*(data->fCorFirstMoment[ekinIndxLow+1]  - data->fCorFirstMoment[ekinIndxLow]);
    corToG2PerG1  += remRfaction*(data->fCorSecondMoment[ekinIndxLow+1] - data->fCorSecondMoment[ekinIndxLow]);
  }
}

Referenced by G4GoudsmitSaundersonMscModel::CrossSectionPerVolume(), and G4GoudsmitSaundersonMscModel::GetTransportMeanFreePath().

◆ Initialise()

void G4GSPWACorrections::Initialise ( )

Definition at line 112 of file G4GSPWACorrections.cc.

                                     {
  // load PWA correction data for each elements that belongs to materials that are used in the detector
  InitDataPerElement();
  // clear  PWA correction data per material
  ClearDataPerMaterial();
  // initialise PWA correction data for the materials that are used in the detector
  InitDataPerMaterials();
}

Referenced by G4GoudsmitSaundersonMscModel::Initialise().

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