Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Attributes | Static Protected Attributes
G4TransitionRadiation Class Referenceabstract

#include <G4TransitionRadiation.hh>

Inheritance diagram for G4TransitionRadiation:
G4VDiscreteProcess G4VProcess G4ForwardXrayTR

Public Member Functions

 G4TransitionRadiation (const G4String &processName="TR", G4ProcessType type=fElectromagnetic)
 
virtual ~G4TransitionRadiation ()
 
G4bool IsApplicable (const G4ParticleDefinition &aParticleType)
 
G4double GetMeanFreePath (const G4Track &, G4double, G4ForceCondition *condition)
 
G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double SpectralAngleTRdensity (G4double energy, G4double varAngle) const =0
 
G4double IntegralOverEnergy (G4double energy1, G4double energy2, G4double varAngle) const
 
G4double IntegralOverAngle (G4double energy, G4double varAngle1, G4double varAngle2) const
 
G4double AngleIntegralDistribution (G4double varAngle1, G4double varAngle2) const
 
G4double EnergyIntegralDistribution (G4double energy1, G4double energy2) const
 
- 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 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 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 fMatIndex1
 
G4int fMatIndex2
 
G4double fGamma
 
G4double fEnergy
 
G4double fVarAngle
 
G4double fMinEnergy
 
G4double fMaxEnergy
 
G4double fMaxTheta
 
G4double fSigma1
 
G4double fSigma2
 
- 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
 

Static Protected Attributes

static const G4int fSympsonNumber = 100
 
static const G4int fGammaNumber = 15
 
static const G4int fPointNumber = 100
 

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 57 of file G4TransitionRadiation.hh.

Constructor & Destructor Documentation

G4TransitionRadiation::G4TransitionRadiation ( const G4String &  processName = "TR",
G4ProcessType  type = fElectromagnetic 
)

Definition at line 60 of file G4TransitionRadiation.cc.

References fEnergy, fGamma, fMatIndex1, fMatIndex2, fMaxEnergy, fMaxTheta, fMinEnergy, fSigma1, fSigma2, fTransitionRadiation, fVarAngle, and G4VProcess::SetProcessSubType().

61  : G4VDiscreteProcess(processName, type)
62 {
63  SetProcessSubType(fTransitionRadiation);
64  fMatIndex1 = fMatIndex2 = 0;
65 
66  fGamma = fEnergy = fVarAngle = fMinEnergy = fMaxEnergy = fMaxTheta = fSigma1 = fSigma2 = 0.0;
67 }
fTransitionRadiation
Definition: G4EmProcessSubType.hh:66
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
G4TransitionRadiation::~G4TransitionRadiation ( )
virtual

Definition at line 74 of file G4TransitionRadiation.cc.

75 {}

Member Function Documentation

G4double G4TransitionRadiation::AngleIntegralDistribution ( G4double  varAngle1,
G4double  varAngle2 
) const

Definition at line 156 of file G4TransitionRadiation.cc.

References fMaxEnergy, fMinEnergy, fSympsonNumber, and IntegralOverEnergy().

158 {
159  G4int i ;
160  G4double h , sumEven = 0.0 , sumOdd = 0.0 ;
161  h = 0.5*(varAngle2 - varAngle1)/fSympsonNumber ;
162  for(i=1;i<fSympsonNumber;i++)
163  {
164  sumEven += IntegralOverEnergy(fMinEnergy,
165  fMinEnergy +0.3*(fMaxEnergy-fMinEnergy),
166  varAngle1 + 2*i*h)
167  + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
168  fMaxEnergy,
169  varAngle1 + 2*i*h);
170  sumOdd += IntegralOverEnergy(fMinEnergy,
171  fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
172  varAngle1 + (2*i - 1)*h)
173  + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
174  fMaxEnergy,
175  varAngle1 + (2*i - 1)*h) ;
176  }
177  sumOdd += IntegralOverEnergy(fMinEnergy,
178  fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
179  varAngle1 + (2*fSympsonNumber - 1)*h)
180  + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
181  fMaxEnergy,
182  varAngle1 + (2*fSympsonNumber - 1)*h) ;
183 
184  return h*(IntegralOverEnergy(fMinEnergy,
185  fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
186  varAngle1)
187  + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
188  fMaxEnergy,
189  varAngle1)
190  + IntegralOverEnergy(fMinEnergy,
191  fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
192  varAngle2)
193  + IntegralOverEnergy(fMinEnergy + 0.3*(fMaxEnergy - fMinEnergy),
194  fMaxEnergy,
195  varAngle2)
196  + 4.0*sumOdd + 2.0*sumEven )/3.0 ;
197 }
G4TransitionRadiation::fSympsonNumber
static const G4int fSympsonNumber
Definition: G4TransitionRadiation.hh:109
G4int
int G4int
Definition: G4Types.hh:78
G4double
double G4double
Definition: G4Types.hh:76
G4TransitionRadiation::IntegralOverEnergy
G4double IntegralOverEnergy(G4double energy1, G4double energy2, G4double varAngle) const
Definition: G4TransitionRadiation.cc:104
G4double G4TransitionRadiation::EnergyIntegralDistribution ( G4double  energy1,
G4double  energy2 
) const

Definition at line 206 of file G4TransitionRadiation.cc.

References fMaxTheta, fSympsonNumber, and IntegralOverAngle().

208 {
209  G4int i ;
210  G4double h , sumEven = 0.0 , sumOdd = 0.0 ;
211  h = 0.5*(energy2 - energy1)/fSympsonNumber ;
212  for(i=1;i<fSympsonNumber;i++)
213  {
214  sumEven += IntegralOverAngle(energy1 + 2*i*h,0.0,0.01*fMaxTheta )
215  + IntegralOverAngle(energy1 + 2*i*h,0.01*fMaxTheta,fMaxTheta);
216  sumOdd += IntegralOverAngle(energy1 + (2*i - 1)*h,0.0,0.01*fMaxTheta)
217  + IntegralOverAngle(energy1 + (2*i - 1)*h,0.01*fMaxTheta,fMaxTheta) ;
218  }
219  sumOdd += IntegralOverAngle(energy1 + (2*fSympsonNumber - 1)*h,
220  0.0,0.01*fMaxTheta)
221  + IntegralOverAngle(energy1 + (2*fSympsonNumber - 1)*h,
222  0.01*fMaxTheta,fMaxTheta) ;
223 
224  return h*(IntegralOverAngle(energy1,0.0,0.01*fMaxTheta)
225  + IntegralOverAngle(energy1,0.01*fMaxTheta,fMaxTheta)
226  + IntegralOverAngle(energy2,0.0,0.01*fMaxTheta)
227  + IntegralOverAngle(energy2,0.01*fMaxTheta,fMaxTheta)
228  + 4.0*sumOdd + 2.0*sumEven )/3.0 ;
229 }
G4TransitionRadiation::fSympsonNumber
static const G4int fSympsonNumber
Definition: G4TransitionRadiation.hh:109
G4int
int G4int
Definition: G4Types.hh:78
G4TransitionRadiation::IntegralOverAngle
G4double IntegralOverAngle(G4double energy, G4double varAngle1, G4double varAngle2) const
Definition: G4TransitionRadiation.cc:130
G4double
double G4double
Definition: G4Types.hh:76
G4double G4TransitionRadiation::GetMeanFreePath ( const G4Track &  ,
G4double  ,
G4ForceCondition *  condition 
)
virtual

Implements G4VDiscreteProcess.

Definition at line 83 of file G4TransitionRadiation.cc.

References DBL_MAX, and Forced.

86 {
87  *condition = Forced;
88  return DBL_MAX; // so TR doesn't limit mean free path
89 }
condition
G4double condition(const G4ErrorSymMatrix &m)
Forced
Definition: G4ForceCondition.hh:54
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4double G4TransitionRadiation::IntegralOverAngle ( G4double  energy,
G4double  varAngle1,
G4double  varAngle2 
) const

Definition at line 130 of file G4TransitionRadiation.cc.

References fSympsonNumber, and SpectralAngleTRdensity().

Referenced by EnergyIntegralDistribution().

133 {
134  G4int i ;
135  G4double h , sumEven = 0.0 , sumOdd = 0.0 ;
136  h = 0.5*(varAngle2 - varAngle1)/fSympsonNumber ;
137  for(i=1;i<fSympsonNumber;i++)
138  {
139  sumEven += SpectralAngleTRdensity(energy,varAngle1 + 2*i*h) ;
140  sumOdd += SpectralAngleTRdensity(energy,varAngle1 + (2*i - 1)*h) ;
141  }
142  sumOdd += SpectralAngleTRdensity(energy,varAngle1 + (2*fSympsonNumber - 1)*h) ;
143 
144  return h*( SpectralAngleTRdensity(energy,varAngle1)
145  + SpectralAngleTRdensity(energy,varAngle2)
146  + 4.0*sumOdd + 2.0*sumEven )/3.0 ;
147 }
G4TransitionRadiation::fSympsonNumber
static const G4int fSympsonNumber
Definition: G4TransitionRadiation.hh:109
G4int
int G4int
Definition: G4Types.hh:78
energy
double precision function energy(A, Z)
Definition: dpm25nuc6.f:4106
G4TransitionRadiation::SpectralAngleTRdensity
virtual G4double SpectralAngleTRdensity(G4double energy, G4double varAngle) const =0
G4double
double G4double
Definition: G4Types.hh:76
G4double G4TransitionRadiation::IntegralOverEnergy ( G4double  energy1,
G4double  energy2,
G4double  varAngle 
) const

Definition at line 104 of file G4TransitionRadiation.cc.

References fSympsonNumber, and SpectralAngleTRdensity().

Referenced by AngleIntegralDistribution().

107 {
108  G4int i ;
109  G4double h , sumEven = 0.0 , sumOdd = 0.0 ;
110  h = 0.5*(energy2 - energy1)/fSympsonNumber ;
111  for(i=1;i<fSympsonNumber;i++)
112  {
113  sumEven += SpectralAngleTRdensity(energy1 + 2*i*h,varAngle) ;
114  sumOdd += SpectralAngleTRdensity(energy1 + (2*i - 1)*h,varAngle) ;
115  }
116  sumOdd += SpectralAngleTRdensity(energy1 + (2*fSympsonNumber - 1)*h,varAngle) ;
117  return h*( SpectralAngleTRdensity(energy1,varAngle)
118  + SpectralAngleTRdensity(energy2,varAngle)
119  + 4.0*sumOdd + 2.0*sumEven )/3.0 ;
120 }
G4TransitionRadiation::fSympsonNumber
static const G4int fSympsonNumber
Definition: G4TransitionRadiation.hh:109
G4int
int G4int
Definition: G4Types.hh:78
G4TransitionRadiation::SpectralAngleTRdensity
virtual G4double SpectralAngleTRdensity(G4double energy, G4double varAngle) const =0
G4double
double G4double
Definition: G4Types.hh:76
G4bool G4TransitionRadiation::IsApplicable ( const G4ParticleDefinition &  aParticleType)
virtual

Reimplemented from G4VProcess.

Definition at line 78 of file G4TransitionRadiation.cc.

References G4ParticleDefinition::GetPDGCharge().

79 {
80  return ( aParticleType.GetPDGCharge() != 0.0 );
81 }
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:163
G4VParticleChange * G4TransitionRadiation::PostStepDoIt ( const G4Track &  ,
const G4Step &   
)
virtual

Reimplemented from G4VDiscreteProcess.

Definition at line 91 of file G4TransitionRadiation.cc.

References G4VProcess::aParticleChange, and G4VProcess::ClearNumberOfInteractionLengthLeft().

93 {
94  ClearNumberOfInteractionLengthLeft();
95  return &aParticleChange;
96 }
G4VProcess::ClearNumberOfInteractionLengthLeft
void ClearNumberOfInteractionLengthLeft()
Definition: G4VProcess.hh:447
G4VProcess::aParticleChange
G4ParticleChange aParticleChange
Definition: G4VProcess.hh:289
virtual G4double G4TransitionRadiation::SpectralAngleTRdensity ( G4double  energy,
G4double  varAngle 
) const
pure virtual

Implemented in G4ForwardXrayTR.

Referenced by IntegralOverAngle(), and IntegralOverEnergy().

Field Documentation

G4double G4TransitionRadiation::fEnergy
protected

Definition at line 105 of file G4TransitionRadiation.hh.

Referenced by G4TransitionRadiation().

G4double G4TransitionRadiation::fGamma
protected

Definition at line 104 of file G4TransitionRadiation.hh.

Referenced by G4TransitionRadiation().

const G4int G4TransitionRadiation::fGammaNumber = 15
staticprotected

Definition at line 110 of file G4TransitionRadiation.hh.

G4int G4TransitionRadiation::fMatIndex1
protected

Definition at line 99 of file G4TransitionRadiation.hh.

Referenced by G4ForwardXrayTR::BuildXrayTRtables(), G4ForwardXrayTR::G4ForwardXrayTR(), G4TransitionRadiation(), and G4ForwardXrayTR::PostStepDoIt().

G4int G4TransitionRadiation::fMatIndex2
protected

Definition at line 100 of file G4TransitionRadiation.hh.

Referenced by G4ForwardXrayTR::BuildXrayTRtables(), G4ForwardXrayTR::G4ForwardXrayTR(), G4TransitionRadiation(), and G4ForwardXrayTR::PostStepDoIt().

G4double G4TransitionRadiation::fMaxEnergy
protected

Definition at line 114 of file G4TransitionRadiation.hh.

Referenced by AngleIntegralDistribution(), and G4TransitionRadiation().

G4double G4TransitionRadiation::fMaxTheta
protected

Definition at line 115 of file G4TransitionRadiation.hh.

Referenced by EnergyIntegralDistribution(), and G4TransitionRadiation().

G4double G4TransitionRadiation::fMinEnergy
protected

Definition at line 113 of file G4TransitionRadiation.hh.

Referenced by AngleIntegralDistribution(), and G4TransitionRadiation().

const G4int G4TransitionRadiation::fPointNumber = 100
staticprotected

Definition at line 111 of file G4TransitionRadiation.hh.

G4double G4TransitionRadiation::fSigma1
protected

Definition at line 117 of file G4TransitionRadiation.hh.

Referenced by G4TransitionRadiation().

G4double G4TransitionRadiation::fSigma2
protected

Definition at line 118 of file G4TransitionRadiation.hh.

Referenced by G4TransitionRadiation().

const G4int G4TransitionRadiation::fSympsonNumber = 100
staticprotected

Definition at line 109 of file G4TransitionRadiation.hh.

Referenced by AngleIntegralDistribution(), EnergyIntegralDistribution(), IntegralOverAngle(), and IntegralOverEnergy().

G4double G4TransitionRadiation::fVarAngle
protected

Definition at line 106 of file G4TransitionRadiation.hh.

Referenced by G4TransitionRadiation().

The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:57:07 for Geant4.10 by   doxygen 1.8.7