G4StrawTubeXTRadiator Class Reference

#include <G4StrawTubeXTRadiator.hh>

Inheritance diagram for G4StrawTubeXTRadiator:

Public Member Functions
	G4StrawTubeXTRadiator (G4LogicalVolume *anEnvelope, G4Material *, G4Material *, G4double, G4double, G4Material *, G4bool unishut=false, const G4String &processName="StrawTubeXTRadiator")
	~G4StrawTubeXTRadiator ()
G4double	GetMediumFormationZone (G4double, G4double, G4double)
void	ComputeMediumPhotoAbsCof ()
G4double	GetMediumLinearPhotoAbs (G4double)
G4complex	GetMediumComplexFZ (G4double, G4double, G4double)
G4double	GetStackFactor (G4double energy, G4double gamma, G4double varAngle) override
void	ProcessDescription (std::ostream &) const override
void	DumpInfo () const override
Public Member Functions inherited from G4VXTRenergyLoss
	G4VXTRenergyLoss (G4LogicalVolume *anEnvelope, G4Material *, G4Material *, G4double, G4double, G4int, const G4String &processName="XTRenergyLoss", G4ProcessType type=fElectromagnetic)
virtual	~G4VXTRenergyLoss ()
	G4VXTRenergyLoss (G4VXTRenergyLoss &)=delete
G4VXTRenergyLoss &	operator= (const G4VXTRenergyLoss &right)=delete
virtual G4bool	IsApplicable (const G4ParticleDefinition &) override
virtual G4VParticleChange *	PostStepDoIt (const G4Track &aTrack, const G4Step &aStep) override
virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition) override
virtual void	BuildPhysicsTable (const G4ParticleDefinition &) override
void	BuildEnergyTable ()
void	BuildAngleForEnergyBank ()
void	BuildTable ()
void	BuildAngleTable ()
void	BuildGlobalAngleTable ()
G4complex	OneInterfaceXTRdEdx (G4double energy, G4double gamma, G4double varAngle)
G4double	SpectralAngleXTRdEdx (G4double varAngle)
virtual G4double	SpectralXTRdEdx (G4double energy)
G4double	AngleSpectralXTRdEdx (G4double energy)
G4double	AngleXTRdEdx (G4double varAngle)
G4double	OneBoundaryXTRNdensity (G4double energy, G4double gamma, G4double varAngle) const
G4double	XTRNSpectralAngleDensity (G4double varAngle)
G4double	XTRNSpectralDensity (G4double energy)
G4double	XTRNAngleSpectralDensity (G4double energy)
G4double	XTRNAngleDensity (G4double varAngle)
void	GetNumberOfPhotons ()
G4double	GetPlateFormationZone (G4double, G4double, G4double)
G4complex	GetPlateComplexFZ (G4double, G4double, G4double)
void	ComputePlatePhotoAbsCof ()
G4double	GetPlateLinearPhotoAbs (G4double)
void	GetPlateZmuProduct ()
G4double	GetPlateZmuProduct (G4double, G4double, G4double)
G4double	GetGasFormationZone (G4double, G4double, G4double)
G4complex	GetGasComplexFZ (G4double, G4double, G4double)
void	ComputeGasPhotoAbsCof ()
G4double	GetGasLinearPhotoAbs (G4double)
void	GetGasZmuProduct ()
G4double	GetGasZmuProduct (G4double, G4double, G4double)
G4double	GetPlateCompton (G4double)
G4double	GetGasCompton (G4double)
G4double	GetComptonPerAtom (G4double, G4double)
G4double	GetXTRrandomEnergy (G4double scaledTkin, G4int iTkin)
G4double	GetXTRenergy (G4int iPlace, G4double position, G4int iTransfer)
G4double	GetRandomAngle (G4double energyXTR, G4int iTkin)
G4double	GetAngleXTR (G4int iTR, G4double position, G4int iAngle)
void	SetGamma (G4double gamma)
G4double	GetGamma ()
void	SetEnergy (G4double energy)
G4double	GetEnergy ()
void	SetVarAngle (G4double varAngle)
G4double	GetVarAngle ()
void	SetCompton (G4bool pC)
G4bool	GetCompton ()
G4int	GetKrange ()
void	SetKrange (G4int kk)
void	SetAlphaGas (G4double ag)
G4double	GetAlphaGas ()
void	SetAlphaPlate (G4double ap)
G4double	GetAlphaPlate ()
void	SetTheMinEnergyTR (G4double minetr)
G4double	GetTheMinEnergyTR ()
void	SetTheMaxEnergyTR (G4double maxetr)
G4double	GetTheMaxEnergyTR ()
void	SetMinEnergyTR (G4double minetr)
G4double	GetMinEnergyTR ()
void	SetMaxEnergyTR (G4double maxetr)
G4double	GetMaxEnergyTR ()
void	SetTheMinAngle (G4double minang)
G4double	GetTheMinAngle ()
void	SetTheMaxAngle (G4double maxang)
G4double	GetTheMaxAngle ()
void	SetMinThetaTR (G4double minatr)
G4double	GetMinThetaTR ()
void	SetMaxThetaTR (G4double maxatr)
G4double	GetMaxThetaTR ()
void	SetFastAngle (G4bool fatr)
G4bool	GetFastAngle ()
void	SetAngleRadDistr (G4bool fatr)
G4bool	GetAngleRadDistr ()
G4PhysicsLogVector *	GetProtonVector ()
G4int	GetTotBin ()
G4PhysicsFreeVector *	GetAngleVector (G4double energy, G4int n)
Public Member Functions inherited from G4VDiscreteProcess
	G4VDiscreteProcess (const G4String &aName, G4ProcessType aType=fNotDefined)
	G4VDiscreteProcess (G4VDiscreteProcess &)
virtual	~G4VDiscreteProcess ()
G4VDiscreteProcess &	operator= (const G4VDiscreteProcess &)=delete
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)
virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)
virtual G4double	GetCrossSection (const G4double, const G4MaterialCutsCouple *)
virtual G4double	MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *)
Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
	G4VProcess (const G4VProcess &right)
virtual	~G4VProcess ()
G4VProcess &	operator= (const G4VProcess &)=delete
G4bool	operator== (const G4VProcess &right) const
G4bool	operator!= (const G4VProcess &right) const
G4double	GetCurrentInteractionLength () const
void	SetPILfactor (G4double value)
G4double	GetPILfactor () const
G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual void	PreparePhysicsTable (const G4ParticleDefinition &)
virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
const G4String &	GetProcessName () const
G4ProcessType	GetProcessType () const
void	SetProcessType (G4ProcessType)
G4int	GetProcessSubType () const
void	SetProcessSubType (G4int)
virtual const G4VProcess *	GetCreatorProcess () const
virtual void	StartTracking (G4Track *)
virtual void	EndTracking ()
virtual void	SetProcessManager (const G4ProcessManager *)
virtual const G4ProcessManager *	GetProcessManager ()
virtual void	ResetNumberOfInteractionLengthLeft ()
G4double	GetNumberOfInteractionLengthLeft () const
G4double	GetTotalNumberOfInteractionLengthTraversed () const
G4bool	isAtRestDoItIsEnabled () const
G4bool	isAlongStepDoItIsEnabled () const
G4bool	isPostStepDoItIsEnabled () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
virtual void	SetMasterProcess (G4VProcess *masterP)
const G4VProcess *	GetMasterProcess () const
virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Protected Attributes
G4SandiaTable *	fMediumPhotoAbsCof
G4double	fSigma3
G4int	fMatIndex3
Protected Attributes inherited from G4VXTRenergyLoss
G4double	fTheMinEnergyTR
G4double	fTheMaxEnergyTR
G4double	fTheMinAngle
G4double	fTheMaxAngle
G4int	fTotBin
G4int	fBinTR
G4int	fKrange
G4ParticleDefinition *	fPtrGamma
G4double *	fGammaCutInKineticEnergy
G4LogicalVolume *	fEnvelope
G4PhysicsTable *	fAngleDistrTable
G4PhysicsTable *	fEnergyDistrTable
G4PhysicsTable *	fAngleForEnergyTable
G4PhysicsLogVector *	fProtonEnergyVector
G4PhysicsLogVector *	fXTREnergyVector
G4SandiaTable *	fPlatePhotoAbsCof
G4SandiaTable *	fGasPhotoAbsCof
G4ParticleChange	fParticleChange
std::vector< G4PhysicsTable * >	fAngleBank
G4double	fGammaTkinCut
G4double	fMinEnergyTR
G4double	fMaxEnergyTR
G4double	fMinThetaTR
G4double	fMaxThetaTR
G4double	fTotalDist
G4double	fPlateThick
G4double	fGasThick
G4double	fAlphaPlate
G4double	fAlphaGas
G4double	fGamma
G4double	fEnergy
G4double	fVarAngle
G4double	fLambda
G4double	fSigma1
G4double	fSigma2
G4int	fMatIndex1
G4int	fMatIndex2
G4int	fPlateNumber
G4bool	fExitFlux
G4bool	fFastAngle
G4bool	fAngleRadDistr
G4bool	fCompton
G4int	secID = -1
Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager = nullptr
G4VParticleChange *	pParticleChange = nullptr
G4ParticleChange	aParticleChange
G4double	theNumberOfInteractionLengthLeft = -1.0
G4double	currentInteractionLength = -1.0
G4double	theInitialNumberOfInteractionLength = -1.0
G4String	theProcessName
G4String	thePhysicsTableFileName
G4ProcessType	theProcessType = fNotDefined
G4int	theProcessSubType = -1
G4double	thePILfactor = 1.0
G4int	verboseLevel = 0
G4bool	enableAtRestDoIt = true
G4bool	enableAlongStepDoIt = true
G4bool	enablePostStepDoIt = true

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)
Protected Member Functions inherited from G4VDiscreteProcess
Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)
void	ClearNumberOfInteractionLengthLeft ()
Static Protected Attributes inherited from G4VXTRenergyLoss
static constexpr G4double	fMinProtonTkin = 100. * CLHEP::GeV
static constexpr G4double	fMaxProtonTkin = 100. * CLHEP::TeV
static constexpr G4double	fPlasmaCof
static constexpr G4double	fCofTR = CLHEP::fine_structure_const / CLHEP::pi

Detailed Description

Definition at line 50 of file G4StrawTubeXTRadiator.hh.

Constructor & Destructor Documentation

G4StrawTubeXTRadiator()

G4StrawTubeXTRadiator::G4StrawTubeXTRadiator ( G4LogicalVolume * anEnvelope, G4Material * foilMat, G4Material * gasMat, G4double a, G4double b, G4Material * mediumMat, G4bool unishut = false, const G4String & processName = "StrawTubeXTRadiator" )

explicit

Definition at line 35 of file G4StrawTubeXTRadiator.cc.

  : G4VXTRenergyLoss(anEnvelope, foilMat, gasMat, a, b, 1, processName)
{
  if(verboseLevel > 0)
    G4cout << "Straw tube X-ray TR  radiator EM process is called" << G4endl;
 
  if(unishut)
  {
    fAlphaPlate = 1. / 3.;
    fAlphaGas   = 12.4;
    if(verboseLevel > 0)
      G4cout << "straw uniform shooting: "
             << "fAlphaPlate = " << fAlphaPlate
             << " ; fAlphaGas = " << fAlphaGas << G4endl;
  }
  else
  {
    fAlphaPlate = 0.5;
    fAlphaGas   = 5.;
    if(verboseLevel > 0)
      G4cout << "straw isotropical shooting: "
             << "fAlphaPlate = " << fAlphaPlate
             << " ; fAlphaGas = " << fAlphaGas << G4endl;
  }
 
  // index of medium material
  fMatIndex3 = (G4int)mediumMat->GetIndex();
  if(verboseLevel > 0)
    G4cout << "medium material = " << mediumMat->GetName() << G4endl;
 
  // plasma energy squared for plate material
  fSigma3 = fPlasmaCof * mediumMat->GetElectronDensity();
  if(verboseLevel > 0)
    G4cout << "medium plasma energy = " << std::sqrt(fSigma3) / eV << " eV"
           << G4endl;
 
  // Compute cofs for preparation of linear photo absorption in external medium
  ComputeMediumPhotoAbsCof();
}

~G4StrawTubeXTRadiator()

G4StrawTubeXTRadiator::~G4StrawTubeXTRadiator ( )

default

Member Function Documentation

ComputeMediumPhotoAbsCof()

void G4StrawTubeXTRadiator::ComputeMediumPhotoAbsCof

(

)

Definition at line 162 of file G4StrawTubeXTRadiator.cc.

{
  const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
  const G4Material* mat                   = (*theMaterialTable)[fMatIndex3];
  fMediumPhotoAbsCof                      = mat->GetSandiaTable();
}

Referenced by G4StrawTubeXTRadiator().

DumpInfo()

void G4StrawTubeXTRadiator::DumpInfo ( ) const

inlineoverridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 69 of file G4StrawTubeXTRadiator.hh.

{ ProcessDescription(G4cout); };

GetMediumComplexFZ()

G4complex G4StrawTubeXTRadiator::GetMediumComplexFZ	(	G4double	omega,
		G4double	gamma,
		G4double	varAngle )

Definition at line 142 of file G4StrawTubeXTRadiator.cc.

{
  G4double cof, length, delta, real_v, image_v;
 
  length = 0.5 * GetMediumFormationZone(omega, gamma, varAngle);
  delta  = length * GetMediumLinearPhotoAbs(omega);
  cof    = 1.0 / (1.0 + delta * delta);
 
  real_v  = length * cof;
  image_v = real_v * delta;
 
  G4complex zone(real_v, image_v);
  return zone;
}

Referenced by GetStackFactor().

GetMediumFormationZone()

G4double G4StrawTubeXTRadiator::GetMediumFormationZone	(	G4double	omega,
		G4double	gamma,
		G4double	varAngle )

Definition at line 130 of file G4StrawTubeXTRadiator.cc.

{
  G4double cof, lambda;
  lambda = 1.0 / gamma / gamma + varAngle + fSigma3 / omega / omega;
  cof    = 2.0 * hbarc / omega / lambda;
  return cof;
}

Referenced by GetMediumComplexFZ().

GetMediumLinearPhotoAbs()

G4double G4StrawTubeXTRadiator::GetMediumLinearPhotoAbs

(

G4double

omega

)

Definition at line 172 of file G4StrawTubeXTRadiator.cc.

{
  G4double omega2, omega3, omega4;
 
  omega2 = omega * omega;
  omega3 = omega2 * omega;
  omega4 = omega2 * omega2;
 
  const G4double* SandiaCof =
    fMediumPhotoAbsCof->GetSandiaCofForMaterial(omega);
 
  G4double cross = SandiaCof[0] / omega + SandiaCof[1] / omega2 +
                   SandiaCof[2] / omega3 + SandiaCof[3] / omega4;
  return cross;
}

Referenced by GetMediumComplexFZ().

GetStackFactor()

G4double G4StrawTubeXTRadiator::GetStackFactor ( G4double energy, G4double gamma, G4double varAngle )

overridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 95 of file G4StrawTubeXTRadiator.cc.

{
  G4double result, L2, L3, M2, M3;
 
  L2 = GetPlateFormationZone(energy, gamma, varAngle);
  L3 = GetGasFormationZone(energy, gamma, varAngle);
 
  M2 = GetPlateLinearPhotoAbs(energy);
  M3 = GetGasLinearPhotoAbs(energy);
 
  G4complex C2(1.0 + 0.5 * fPlateThick * M2 / fAlphaPlate,
               fPlateThick / L2 / fAlphaPlate);
  G4complex C3(1.0 + 0.5 * fGasThick * M3 / fAlphaGas,
               fGasThick / L3 / fAlphaGas);
 
  G4complex H2 = std::pow(C2, -fAlphaPlate);
  G4complex H3 = std::pow(C3, -fAlphaGas);
  G4complex H  = H2 * H3;
 
  G4complex Z1 = GetMediumComplexFZ(energy, gamma, varAngle);
  G4complex Z2 = GetPlateComplexFZ(energy, gamma, varAngle);
  G4complex Z3 = GetGasComplexFZ(energy, gamma, varAngle);
 
  G4complex R = (Z1 - Z2) * (Z1 - Z2) * (1. - H2 * H) +
                (Z2 - Z3) * (Z2 - Z3) * (1. - H3) +
                2. * (Z1 - Z2) * (Z2 - Z3) * H2 * (1. - H3);
 
  result = 2.0 * std::real(R) * (varAngle * energy / hbarc / hbarc);
 
  return result;
}

ProcessDescription()

void G4StrawTubeXTRadiator::ProcessDescription ( std::ostream & out ) const

overridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 83 of file G4StrawTubeXTRadiator.cc.

{
  out << "Simulation of forward X-ray transition radiation for the case of\n"
         "a straw tube radiator.\n";
}

Referenced by DumpInfo().

Member Data Documentation

fMatIndex3

G4int G4StrawTubeXTRadiator::fMatIndex3

protected

Definition at line 76 of file G4StrawTubeXTRadiator.hh.

Referenced by ComputeMediumPhotoAbsCof(), and G4StrawTubeXTRadiator().

fMediumPhotoAbsCof

G4SandiaTable* G4StrawTubeXTRadiator::fMediumPhotoAbsCof

protected

Definition at line 72 of file G4StrawTubeXTRadiator.hh.

Referenced by ComputeMediumPhotoAbsCof(), and GetMediumLinearPhotoAbs().

fSigma3

G4double G4StrawTubeXTRadiator::fSigma3

protected

Definition at line 74 of file G4StrawTubeXTRadiator.hh.

Referenced by G4StrawTubeXTRadiator(), and GetMediumFormationZone().

---

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

• G4StrawTubeXTRadiator.hh
• G4StrawTubeXTRadiator.cc