#include <G4XTRGaussRadModel.hh>

Inheritance diagram for G4XTRGaussRadModel:

Public Member Functions
	G4XTRGaussRadModel (G4LogicalVolume anEnvelope, G4double, G4double, G4Material , G4Material *, G4double, G4double, G4int, const G4String &processName="XTRGaussRadModel")

	~G4XTRGaussRadModel () override

G4double	SpectralXTRdEdx (G4double energy) override

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)

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

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)

void	ClearNumberOfInteractionLengthLeft ()

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

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 48 of file G4XTRGaussRadModel.hh.

Constructor & Destructor Documentation

◆ G4XTRGaussRadModel()

G4XTRGaussRadModel::G4XTRGaussRadModel	(	G4LogicalVolume *	anEnvelope,
		G4double	alphaPlate,
		G4double	alphaGas,
		G4Material *	foilMat,
		G4Material *	gasMat,
		G4double	aa,
		G4double	b,
		G4int	n,
		const G4String &	processName = "XTRGaussRadModel" )

explicit

Definition at line 42 of file G4XTRGaussRadModel.cc.

  : G4VXTRenergyLoss(anEnvelope, foilMat, gasMat, aa, b, n, processName)
{
  if( verboseLevel > 0 )
    G4cout << "G4XTRGaussRadModel EM process is called"
           << G4endl;
 
  fAlphaPlate = alphaPlate; // ~40
  fAlphaGas   = alphaGas;   // ~10 
  fExitFlux   = true;       // XTR photons are moved to the end of radiator
}

◆ ~G4XTRGaussRadModel()

G4XTRGaussRadModel::~G4XTRGaussRadModel ( )

overridedefault

Member Function Documentation

◆ DumpInfo()

void G4XTRGaussRadModel::DumpInfo ( ) const

inlineoverridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 64 of file G4XTRGaussRadModel.hh.

64{ ProcessDescription(G4cout); };

G4XTRGaussRadModel::ProcessDescription

void ProcessDescription(std::ostream &) const override

Definition G4XTRGaussRadModel.cc:62

◆ GetStackFactor()

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

overridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 143 of file G4XTRGaussRadModel.cc.

{
  G4double result(0.); 
 
  G4double Ma, Mb, aMa, bMb, sigma;
 
  G4double sa = fPlateThick/fAlphaPlate;
  G4double sb = fGasThick/fAlphaGas;
 
  Ma =  GetPlateLinearPhotoAbs(energy);  
  aMa = fPlateThick * Ma;
  Mb = GetGasLinearPhotoAbs(energy);  
  bMb   = fGasThick * Mb; 
  sigma = aMa + bMb; // m1*t1+m2*t2 dimensionless
  G4double nn = G4double( fPlateNumber );
 
  // Gauss fluctuation of foil and gas gaps according to
  // RMS = sa = a/fAlphaPlate and RMS = sb = b/fAlphaGas
 
  G4complex med(0.,1.);
  G4complex Z1   = GetPlateComplexFZ( energy, gamma, varAngle);
  G4complex por1 = -0.5*med*fPlateThick/Z1 - 0.125*sa*sa/Z1/Z1;
  
  G4complex Z2   = GetGasComplexFZ( energy, gamma, varAngle);
  G4complex por2 = -0.5*med*fGasThick/Z2 - 0.125*sb*sb/Z2/Z2;
  
  G4complex npor = (por1+por2)*nn;
  
  G4complex Ha = exp(por1); 
  G4complex Hb = exp(por2);
  
  G4complex H  = Ha*Hb;
  G4complex Hn = exp(npor); 
 
  G4complex A  = ( 1.0 - Ha ) * ( 1.0 - Hb ) / ( 1. - H );
  G4complex B1 = ( 1.0 - Ha ) * ( 1.0 - Ha ) * Hb / ( 1. - H );
  
  G4complex R  = B1*( exp(-sigma*nn) - Hn )/( exp(-sigma) - H );
 
  R += A * ( 1 - exp(-sigma*nn) ) / ( 1. - exp(-sigma) ); // -> A*nn
  
  R *= OneInterfaceXTRdEdx(energy, gamma, varAngle);
  
  result      = 2.0 * std::real(R);
 
  return result;
}

◆ ProcessDescription()

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

overridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 62 of file G4XTRGaussRadModel.cc.

{
  out << "Simulation of forward X-ray transition radiation generated by\n"
         "relativistic charged particles crossing the interface between\n"
         "two materials.\n";
}

Referenced by DumpInfo().

◆ SpectralXTRdEdx()

G4double G4XTRGaussRadModel::SpectralXTRdEdx ( G4double energy )

overridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 71 of file G4XTRGaussRadModel.cc.

{
  static constexpr G4double cofPHC = 4. * pi * hbarc;
  G4double result, sum = 0., tmp, cof1, cof2, cofMin, theta2, theta2k;
  G4double aMa, bMb, sigma, dump;
  G4int k, kMax, kMin;
 
  aMa   = fPlateThick * GetPlateLinearPhotoAbs(energy);
  bMb   = fGasThick * GetGasLinearPhotoAbs(energy);
  sigma = 0.5 * (aMa + bMb);
  dump  = std::exp(-fPlateNumber * sigma);
  if(verboseLevel > 2)
    G4cout << " dump = " << dump << G4endl;
  tmp  = (fSigma1 - fSigma2) / cofPHC / energy;
  cof1 = fPlateThick * tmp;
  cof2 = fGasThick * tmp;
 
  cofMin = energy * (fPlateThick + fGasThick) / fGamma / fGamma;
  cofMin += (fPlateThick * fSigma1 + fGasThick * fSigma2) / energy;
  cofMin /= cofPHC;
 
  theta2 = cofPHC / (energy * (fPlateThick + fGasThick));
 
  kMin = G4int(cofMin);
  if(cofMin > kMin)
    kMin++;
 
  kMax = kMin + 200; // 99; // 49; //
 
  if(verboseLevel > 2)
  {
    G4cout << cof1 << "     " << cof2 << "        " << cofMin << G4endl;
    G4cout << "kMin = " << kMin << ";    kMax = " << kMax << G4endl;
  }
  for(k = kMin; k <= kMax; ++k)
  {
    tmp    = pi * fPlateThick * (k + cof2) / (fPlateThick + fGasThick);
    result = (k - cof1) * (k - cof1) * (k + cof2) * (k + cof2);
    if(k == kMin && kMin == G4int(cofMin))
    {
      sum +=
        0.5 * std::sin(tmp) * std::sin(tmp) * std::abs(k - cofMin) / result;
    }
    else
    {
      sum += std::sin(tmp) * std::sin(tmp) * std::abs(k - cofMin) / result;
    }
    theta2k = std::sqrt(theta2 * std::abs(k - cofMin));
 
    if(verboseLevel > 2)
    {
      G4cout << k << "   " << theta2k << "     "
             << std::sin(tmp) * std::sin(tmp) * std::abs(k - cofMin) / result
             << "      " << sum << G4endl;
    }
  }
  result = 2 * (cof1 + cof2) * (cof1 + cof2) * sum / energy;
  result *= dump * (-1 + dump + 2 * fPlateNumber);
  return result;
}

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4XTRGaussRadModel()

◆ ~G4XTRGaussRadModel()

Member Function Documentation

◆ DumpInfo()

◆ GetStackFactor()

◆ ProcessDescription()

◆ SpectralXTRdEdx()