#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")

virtual	~G4StrawTubeXTRadiator ()

G4double	GetMediumFormationZone (G4double, G4double, G4double)

void	ComputeMediumPhotoAbsCof ()

G4double	GetMediumLinearPhotoAbs (G4double)

G4complex	GetMediumComplexFZ (G4double, G4double, G4double)

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

Public Member Functions inherited from G4VXTRenergyLoss
	G4VXTRenergyLoss (G4LogicalVolume anEnvelope, G4Material , G4Material *, G4double, G4double, G4int, const G4String &processName="XTRenergyLoss", G4ProcessType type=fElectromagnetic)

virtual	~G4VXTRenergyLoss ()

G4bool	IsApplicable (const G4ParticleDefinition &)

G4VParticleChange *	PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)

G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)

void	BuildPhysicsTable (const G4ParticleDefinition &)

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)

G4double	GetGamma ()

G4double	GetEnergy ()

G4double	GetVarAngle ()

void	SetGamma (G4double gamma)

void	SetEnergy (G4double energy)

void	SetVarAngle (G4double varAngle)

void	SetAngleRadDistr (G4bool pAngleRadDistr)

void	SetCompton (G4bool pC)

G4PhysicsLogVector *	GetProtonVector ()

G4int	GetTotBin ()

G4PhysicsFreeVector *	GetAngleVector (G4double energy, G4int n)

Public Member Functions inherited from G4VDiscreteProcess
	G4VDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)

	G4VDiscreteProcess (G4VDiscreteProcess &)

virtual	~G4VDiscreteProcess ()

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

Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)

	G4VProcess (const G4VProcess &right)

virtual	~G4VProcess ()

G4int	operator== (const G4VProcess &right) const

G4int	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 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

virtual void	DumpInfo () 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
G4int	fMatIndex3

G4double	fSigma3

G4SandiaTable *	fMediumPhotoAbsCof

Protected Attributes inherited from G4VXTRenergyLoss
G4ParticleDefinition *	fPtrGamma

G4double *	fGammaCutInKineticEnergy

G4double	fGammaTkinCut

G4LogicalVolume *	fEnvelope

G4PhysicsTable *	fAngleDistrTable

G4PhysicsTable *	fEnergyDistrTable

G4PhysicsLogVector *	fProtonEnergyVector

G4PhysicsLogVector *	fXTREnergyVector

G4double	fTheMinEnergyTR

G4double	fTheMaxEnergyTR

G4double	fMinEnergyTR

G4double	fMaxEnergyTR

G4double	fTheMaxAngle

G4double	fTheMinAngle

G4double	fMaxThetaTR

G4int	fBinTR

G4double	fMinProtonTkin

G4double	fMaxProtonTkin

G4int	fTotBin

G4double	fGamma

G4double	fEnergy

G4double	fVarAngle

G4double	fLambda

G4double	fPlasmaCof

G4double	fCofTR

G4bool	fExitFlux

G4bool	fAngleRadDistr

G4bool	fCompton

G4double	fSigma1

G4double	fSigma2

G4int	fMatIndex1

G4int	fMatIndex2

G4int	fPlateNumber

G4double	fTotalDist

G4double	fPlateThick

G4double	fGasThick

G4double	fAlphaPlate

G4double	fAlphaGas

G4SandiaTable *	fPlatePhotoAbsCof

G4SandiaTable *	fGasPhotoAbsCof

G4ParticleChange	fParticleChange

G4PhysicsTable *	fAngleForEnergyTable

std::vector< G4PhysicsTable * >	fAngleBank

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager

G4VParticleChange *	pParticleChange

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft

G4double	currentInteractionLength

G4double	theInitialNumberOfInteractionLength

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType

G4int	theProcessSubType

G4double	thePILfactor

G4bool	enableAtRestDoIt

G4bool	enableAlongStepDoIt

G4bool	enablePostStepDoIt

G4int	verboseLevel

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

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

void	ClearNumberOfInteractionLengthLeft ()

Detailed Description

Definition at line 53 of file G4StrawTubeXTRadiator.hh.

Constructor & Destructor Documentation

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

Definition at line 40 of file G4StrawTubeXTRadiator.cc.

References ComputeMediumPhotoAbsCof(), python.hepunit::eV, G4VXTRenergyLoss::fAlphaGas, G4VXTRenergyLoss::fAlphaPlate, fMatIndex3, G4VXTRenergyLoss::fPlasmaCof, fSigma3, G4cout, G4endl, G4Material::GetElectronDensity(), G4Material::GetIndex(), G4Material::GetName(), and G4VProcess::verboseLevel.

                                                                       :
   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 = 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();
 
   // Build energy and angular integral spectra of X-ray TR photons from
   // a radiator
 
   // BuildTable();
 }

G4StrawTubeXTRadiator::~G4StrawTubeXTRadiator ( )

virtual

Definition at line 93 of file G4StrawTubeXTRadiator.cc.

94 {

95 }

Member Function Documentation

void G4StrawTubeXTRadiator::ComputeMediumPhotoAbsCof ( )

Definition at line 184 of file G4StrawTubeXTRadiator.cc.

References fMatIndex3, fMediumPhotoAbsCof, G4Material::GetMaterialTable(), and G4Material::GetSandiaTable().

Referenced by G4StrawTubeXTRadiator().

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

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

Definition at line 162 of file G4StrawTubeXTRadiator.cc.

References GetMediumFormationZone(), and GetMediumLinearPhotoAbs().

Referenced by GetStackFactor().

 {
   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;
 }

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

Definition at line 148 of file G4StrawTubeXTRadiator.cc.

References fSigma3, python.hepunit::hbarc, and G4InuclParticleNames::lambda.

Referenced by GetMediumComplexFZ().

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

G4double G4StrawTubeXTRadiator::GetMediumLinearPhotoAbs ( G4double omega )

Definition at line 196 of file G4StrawTubeXTRadiator.cc.

References fMediumPhotoAbsCof, and G4SandiaTable::GetSandiaCofForMaterial().

Referenced by GetMediumComplexFZ().

 {
   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;
 }

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

virtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 106 of file G4StrawTubeXTRadiator.cc.

References C3, G4VXTRenergyLoss::fAlphaGas, G4VXTRenergyLoss::fAlphaPlate, G4VXTRenergyLoss::fGasThick, G4VXTRenergyLoss::fPlateThick, G4VXTRenergyLoss::GetGasComplexFZ(), G4VXTRenergyLoss::GetGasFormationZone(), G4VXTRenergyLoss::GetGasLinearPhotoAbs(), GetMediumComplexFZ(), G4VXTRenergyLoss::GetPlateComplexFZ(), G4VXTRenergyLoss::GetPlateFormationZone(), G4VXTRenergyLoss::GetPlateLinearPhotoAbs(), and python.hepunit::hbarc.

 {
 
 
   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;
 
 }