Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
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 ()
 
virtual void ProcessDescription (std::ostream &) const override
 
virtual void DumpInfo () const override
 
 G4VXTRenergyLoss (G4VXTRenergyLoss &)=delete
 
G4VXTRenergyLossoperator= (const G4VXTRenergyLoss &right)=delete
 
virtual G4double GetStackFactor (G4double energy, G4double gamma, G4double varAngle)
 
virtual G4bool IsApplicable (const G4ParticleDefinition &) override
 
virtual G4VParticleChangePostStepDoIt (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 ()
 
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 ()
 
G4PhysicsLogVectorGetProtonVector ()
 
G4int GetTotBin ()
 
G4PhysicsFreeVectorGetAngleVector (G4double energy, G4int n)
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &aName, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
G4VDiscreteProcessoperator= (const G4VDiscreteProcess &)=delete
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4VParticleChangePostStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChangeAtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChangeAlongStepDoIt (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 ()
 
G4VProcessoperator= (const G4VProcess &)=delete
 
G4bool operator== (const G4VProcess &right) const
 
G4bool operator!= (const G4VProcess &right) const
 
virtual G4VParticleChangePostStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChangeAlongStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChangeAtRestDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0
 
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 G4bool IsApplicable (const G4ParticleDefinition &)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void PreparePhysicsTable (const G4ParticleDefinition &)
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
const G4StringGetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
 
const G4StringGetProcessName () const
 
G4ProcessType GetProcessType () const
 
void SetProcessType (G4ProcessType)
 
G4int GetProcessSubType () const
 
void SetProcessSubType (G4int)
 
virtual const G4VProcessGetCreatorProcess () const
 
virtual void StartTracking (G4Track *)
 
virtual void EndTracking ()
 
virtual void SetProcessManager (const G4ProcessManager *)
 
virtual const G4ProcessManagerGetProcessManager ()
 
virtual void ResetNumberOfInteractionLengthLeft ()
 
G4double GetNumberOfInteractionLengthLeft () const
 
G4double GetTotalNumberOfInteractionLengthTraversed () const
 
G4bool isAtRestDoItIsEnabled () const
 
G4bool isAlongStepDoItIsEnabled () const
 
G4bool isPostStepDoItIsEnabled () const
 
virtual void DumpInfo () const
 
virtual void ProcessDescription (std::ostream &outfile) const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcessGetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Protected Attributes

G4SandiaTablefMediumPhotoAbsCof
 
G4double fSigma3
 
G4int fMatIndex3
 
- Protected Attributes inherited from G4VXTRenergyLoss
G4double fTheMinEnergyTR
 
G4double fTheMaxEnergyTR
 
G4double fTheMinAngle
 
G4double fTheMaxAngle
 
G4int fTotBin
 
G4int fBinTR
 
G4ParticleDefinitionfPtrGamma
 
G4doublefGammaCutInKineticEnergy
 
G4LogicalVolumefEnvelope
 
G4PhysicsTablefAngleDistrTable
 
G4PhysicsTablefEnergyDistrTable
 
G4PhysicsTablefAngleForEnergyTable
 
G4PhysicsLogVectorfProtonEnergyVector
 
G4PhysicsLogVectorfXTREnergyVector
 
G4SandiaTablefPlatePhotoAbsCof
 
G4SandiaTablefGasPhotoAbsCof
 
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 G4ProcessManageraProcessManager = nullptr
 
G4VParticleChangepParticleChange = 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 G4StringGetProcessTypeName (G4ProcessType)
 
virtual G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0
 
- 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.

41 : G4VXTRenergyLoss(anEnvelope, foilMat, gasMat, a, b, 1, processName)
42{
43 if(verboseLevel > 0)
44 G4cout << "Straw tube X-ray TR radiator EM process is called" << G4endl;
45
46 if(unishut)
47 {
48 fAlphaPlate = 1. / 3.;
49 fAlphaGas = 12.4;
50 if(verboseLevel > 0)
51 G4cout << "straw uniform shooting: "
52 << "fAlphaPlate = " << fAlphaPlate
53 << " ; fAlphaGas = " << fAlphaGas << G4endl;
54 }
55 else
56 {
57 fAlphaPlate = 0.5;
58 fAlphaGas = 5.;
59 if(verboseLevel > 0)
60 G4cout << "straw isotropical shooting: "
61 << "fAlphaPlate = " << fAlphaPlate
62 << " ; fAlphaGas = " << fAlphaGas << G4endl;
63 }
64
65 // index of medium material
66 fMatIndex3 = (G4int)mediumMat->GetIndex();
67 if(verboseLevel > 0)
68 G4cout << "medium material = " << mediumMat->GetName() << G4endl;
69
70 // plasma energy squared for plate material
71 fSigma3 = fPlasmaCof * mediumMat->GetElectronDensity();
72 if(verboseLevel > 0)
73 G4cout << "medium plasma energy = " << std::sqrt(fSigma3) / eV << " eV"
74 << G4endl;
75
76 // Compute cofs for preparation of linear photo absorption in external medium
78}
int G4int
Definition: G4Types.hh:85
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
G4double GetElectronDensity() const
Definition: G4Material.hh:212
const G4String & GetName() const
Definition: G4Material.hh:172
size_t GetIndex() const
Definition: G4Material.hh:255
G4int verboseLevel
Definition: G4VProcess.hh:360
static constexpr G4double fPlasmaCof

◆ ~G4StrawTubeXTRadiator()

G4StrawTubeXTRadiator::~G4StrawTubeXTRadiator ( )
default

Member Function Documentation

◆ ComputeMediumPhotoAbsCof()

void G4StrawTubeXTRadiator::ComputeMediumPhotoAbsCof ( )

Definition at line 162 of file G4StrawTubeXTRadiator.cc.

163{
164 const G4MaterialTable* theMaterialTable = G4Material::GetMaterialTable();
165 const G4Material* mat = (*theMaterialTable)[fMatIndex3];
167}
std::vector< G4Material * > G4MaterialTable
G4SandiaTable * GetSandiaTable() const
Definition: G4Material.hh:224
static G4MaterialTable * GetMaterialTable()
Definition: G4Material.cc:677
G4SandiaTable * fMediumPhotoAbsCof

Referenced by G4StrawTubeXTRadiator().

◆ DumpInfo()

void G4StrawTubeXTRadiator::DumpInfo ( ) const
inlineoverridevirtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 69 of file G4StrawTubeXTRadiator.hh.

void ProcessDescription(std::ostream &) const override

◆ GetMediumComplexFZ()

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

Definition at line 142 of file G4StrawTubeXTRadiator.cc.

145{
146 G4double cof, length, delta, real_v, image_v;
147
148 length = 0.5 * GetMediumFormationZone(omega, gamma, varAngle);
149 delta = length * GetMediumLinearPhotoAbs(omega);
150 cof = 1.0 / (1.0 + delta * delta);
151
152 real_v = length * cof;
153 image_v = real_v * delta;
154
155 G4complex zone(real_v, image_v);
156 return zone;
157}
double G4double
Definition: G4Types.hh:83
std::complex< G4double > G4complex
Definition: G4Types.hh:88
G4double GetMediumFormationZone(G4double, G4double, G4double)
G4double GetMediumLinearPhotoAbs(G4double)

Referenced by GetStackFactor().

◆ GetMediumFormationZone()

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

Definition at line 130 of file G4StrawTubeXTRadiator.cc.

133{
134 G4double cof, lambda;
135 lambda = 1.0 / gamma / gamma + varAngle + fSigma3 / omega / omega;
136 cof = 2.0 * hbarc / omega / lambda;
137 return cof;
138}

Referenced by GetMediumComplexFZ().

◆ GetMediumLinearPhotoAbs()

G4double G4StrawTubeXTRadiator::GetMediumLinearPhotoAbs ( G4double  omega)

Definition at line 172 of file G4StrawTubeXTRadiator.cc.

173{
174 G4double omega2, omega3, omega4;
175
176 omega2 = omega * omega;
177 omega3 = omega2 * omega;
178 omega4 = omega2 * omega2;
179
180 const G4double* SandiaCof =
182
183 G4double cross = SandiaCof[0] / omega + SandiaCof[1] / omega2 +
184 SandiaCof[2] / omega3 + SandiaCof[3] / omega4;
185 return cross;
186}
G4double GetSandiaCofForMaterial(G4int, G4int) const

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.

97{
98 G4double result, L2, L3, M2, M3;
99
100 L2 = GetPlateFormationZone(energy, gamma, varAngle);
101 L3 = GetGasFormationZone(energy, gamma, varAngle);
102
103 M2 = GetPlateLinearPhotoAbs(energy);
104 M3 = GetGasLinearPhotoAbs(energy);
105
106 G4complex C2(1.0 + 0.5 * fPlateThick * M2 / fAlphaPlate,
107 fPlateThick / L2 / fAlphaPlate);
108 G4complex C3(1.0 + 0.5 * fGasThick * M3 / fAlphaGas,
109 fGasThick / L3 / fAlphaGas);
110
111 G4complex H2 = std::pow(C2, -fAlphaPlate);
112 G4complex H3 = std::pow(C3, -fAlphaGas);
113 G4complex H = H2 * H3;
114
115 G4complex Z1 = GetMediumComplexFZ(energy, gamma, varAngle);
116 G4complex Z2 = GetPlateComplexFZ(energy, gamma, varAngle);
117 G4complex Z3 = GetGasComplexFZ(energy, gamma, varAngle);
118
119 G4complex R = (Z1 - Z2) * (Z1 - Z2) * (1. - H2 * H) +
120 (Z2 - Z3) * (Z2 - Z3) * (1. - H3) +
121 2. * (Z1 - Z2) * (Z2 - Z3) * H2 * (1. - H3);
122
123 result = 2.0 * std::real(R) * (varAngle * energy / hbarc / hbarc);
124
125 return result;
126}
const double C2
#define C3
G4complex GetMediumComplexFZ(G4double, G4double, G4double)
G4double GetPlateLinearPhotoAbs(G4double)
G4double GetGasFormationZone(G4double, G4double, G4double)
G4complex GetPlateComplexFZ(G4double, G4double, G4double)
G4double GetPlateFormationZone(G4double, G4double, G4double)
G4double GetGasLinearPhotoAbs(G4double)
G4complex GetGasComplexFZ(G4double, G4double, G4double)
G4double energy(const ThreeVector &p, const G4double m)

◆ ProcessDescription()

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

Reimplemented from G4VXTRenergyLoss.

Definition at line 83 of file G4StrawTubeXTRadiator.cc.

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

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: