Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions
G4RegularXTRadiator Class Reference

#include <G4RegularXTRadiator.hh>

Inheritance diagram for G4RegularXTRadiator:
G4VXTRenergyLoss G4VDiscreteProcess G4VProcess

Public Member Functions

 G4RegularXTRadiator (G4LogicalVolume *anEnvelope, G4Material *, G4Material *, G4double, G4double, G4int, const G4String &processName="XTRegularRadiator")
 
 ~G4RegularXTRadiator ()
 
G4double SpectralXTRdEdx (G4double energy)
 
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)
 
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 &)
 

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

Detailed Description

Definition at line 50 of file G4RegularXTRadiator.hh.

Constructor & Destructor Documentation

G4RegularXTRadiator::G4RegularXTRadiator ( G4LogicalVolume *  anEnvelope,
G4Material *  foilMat,
G4Material *  gasMat,
G4double  a,
G4double  b,
G4int  n,
const G4String &  processName = "XTRegularRadiator" 
)

Definition at line 42 of file G4RegularXTRadiator.cc.

References G4VXTRenergyLoss::fAlphaGas, G4VXTRenergyLoss::fAlphaPlate, G4cout, and G4endl.

45  :
46  G4VXTRenergyLoss(anEnvelope,foilMat,gasMat,a,b,n,processName)
47 {
48  G4cout<<"Regular X-ray TR radiator EM process is called"<<G4endl ;
49 
50  // Build energy and angular integral spectra of X-ray TR photons from
51  // a radiator
52 
53  fAlphaPlate = 10000;
54  fAlphaGas = 1000;
55  G4cout<<"fAlphaPlate = "<<fAlphaPlate<<" ; fAlphaGas = "<<fAlphaGas<<G4endl ;
56 
57  // BuildTable() ;
58 }
G4cout
G4GLOB_DLL std::ostream G4cout
n
const G4int n
Definition: G4UrQMD1_3Interface.hh:144
G4VXTRenergyLoss::G4VXTRenergyLoss
G4VXTRenergyLoss(G4LogicalVolume *anEnvelope, G4Material *, G4Material *, G4double, G4double, G4int, const G4String &processName="XTRenergyLoss", G4ProcessType type=fElectromagnetic)
Definition: G4VXTRenergyLoss.cc:63
G4endl
#define G4endl
Definition: G4ios.hh:61
G4RegularXTRadiator::~G4RegularXTRadiator ( )

Definition at line 62 of file G4RegularXTRadiator.cc.

63 {
64  ;
65 }

Member Function Documentation

G4double G4RegularXTRadiator::GetStackFactor ( G4double  energy,
G4double  gamma,
G4double  varAngle 
)
virtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 171 of file G4RegularXTRadiator.cc.

References G4VXTRenergyLoss::fAlphaGas, G4VXTRenergyLoss::fAlphaPlate, G4VXTRenergyLoss::fGasThick, G4VXTRenergyLoss::fPlateNumber, G4VXTRenergyLoss::fPlateThick, G4VXTRenergyLoss::GetGasFormationZone(), G4VXTRenergyLoss::GetGasLinearPhotoAbs(), G4VXTRenergyLoss::GetPlateFormationZone(), G4VXTRenergyLoss::GetPlateLinearPhotoAbs(), and G4VXTRenergyLoss::OneInterfaceXTRdEdx().

173 {
174 
175  // some gamma (10000/1000) like algorithm
176 
177  G4double result, Za, Zb, Ma, Mb;
178 
179  Za = GetPlateFormationZone(energy,gamma,varAngle);
180  Zb = GetGasFormationZone(energy,gamma,varAngle);
181 
182  Ma = GetPlateLinearPhotoAbs(energy);
183  Mb = GetGasLinearPhotoAbs(energy);
184 
185 
186  G4complex Ca(1.0+0.5*fPlateThick*Ma/fAlphaPlate,fPlateThick/Za/fAlphaPlate);
187  G4complex Cb(1.0+0.5*fGasThick*Mb/fAlphaGas,fGasThick/Zb/fAlphaGas);
188 
189  G4complex Ha = std::pow(Ca,-fAlphaPlate);
190  G4complex Hb = std::pow(Cb,-fAlphaGas);
191  G4complex H = Ha*Hb;
192 
193  G4complex F1 = (1.0 - Ha)*(1.0 - Hb )/(1.0 - H)
194  * G4double(fPlateNumber);
195 
196  G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
197  * (1.0 - std::pow(H,fPlateNumber));
198 
199  G4complex R = (F1 + F2)*OneInterfaceXTRdEdx(energy,gamma,varAngle);
200 
201  result = 2.0*std::real(R);
202 
203  return result;
204 
205  /*
206  // numerically stable but slow algorithm
207 
208  G4double result, Qa, Qb, Q, aZa, bZb, aMa, bMb; // , D;
209 
210  aZa = fPlateThick/GetPlateFormationZone(energy,gamma,varAngle);
211  bZb = fGasThick/GetGasFormationZone(energy,gamma,varAngle);
212  aMa = fPlateThick*GetPlateLinearPhotoAbs(energy);
213  bMb = fGasThick*GetGasLinearPhotoAbs(energy);
214  Qa = std::exp(-aMa);
215  Qb = std::exp(-bMb);
216  Q = Qa*Qb;
217  G4complex Ha( std::exp(-0.5*aMa)*std::cos(aZa),
218  -std::exp(-0.5*aMa)*std::sin(aZa) );
219  G4complex Hb( std::exp(-0.5*bMb)*std::cos(bZb),
220  -std::exp(-0.5*bMb)*std::sin(bZb) );
221  G4complex H = Ha*Hb;
222 
223  G4complex Hs = conj(H);
224  D = 1.0 /( (1 - std::sqrt(Q))*(1 - std::sqrt(Q)) +
225  4*std::sqrt(Q)*std::sin(0.5*(aZa+bZb))*std::sin(0.5*(aZa+bZb)) );
226  G4complex F1 = (1.0 - Ha)*(1.0 - Hb)*(1.0 - Hs)
227  * G4double(fPlateNumber)*D;
228  G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb*(1.0-Hs)*(1.0-Hs)
229  * (1.0 - std::pow(H,fPlateNumber)) * D*D;
230  G4complex R = (F1 + F2)*OneInterfaceXTRdEdx(energy,gamma,varAngle);
231 
232 
233  G4complex S(0.,0.), c(1.,0.);
234  G4int k;
235  for(k = 1; k < fPlateNumber; k++)
236  {
237  c *= H;
238  S += ( G4double(fPlateNumber) - G4double(k) )*c;
239  }
240  G4complex R = (2.- Ha - 1./Ha)*S + (1. - Ha)*G4double(fPlateNumber);
241  R *= OneInterfaceXTRdEdx(energy,gamma,varAngle);
242  result = 2.0*std::real(R);
243  return result;
244  */
245 }
G4VXTRenergyLoss::GetGasFormationZone
G4double GetGasFormationZone(G4double, G4double, G4double)
Definition: G4VXTRenergyLoss.cc:1105
G4VXTRenergyLoss::GetPlateLinearPhotoAbs
G4double GetPlateLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1085
G4VXTRenergyLoss::GetGasLinearPhotoAbs
G4double GetGasLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1157
G4VXTRenergyLoss::GetPlateFormationZone
G4double GetPlateFormationZone(G4double, G4double, G4double)
Definition: G4VXTRenergyLoss.cc:1033
energy
double precision function energy(A, Z)
Definition: dpm25nuc6.f:4106
G4complex
std::complex< G4double > G4complex
Definition: G4Types.hh:81
G4VXTRenergyLoss::OneInterfaceXTRdEdx
G4complex OneInterfaceXTRdEdx(G4double energy, G4double gamma, G4double varAngle)
Definition: G4VXTRenergyLoss.cc:858
G4double
double G4double
Definition: G4Types.hh:76
G4double G4RegularXTRadiator::SpectralXTRdEdx ( G4double  energy)
virtual

Reimplemented from G4VXTRenergyLoss.

Definition at line 71 of file G4RegularXTRadiator.cc.

References energy(), G4VXTRenergyLoss::fGamma, G4VXTRenergyLoss::fGasThick, G4VXTRenergyLoss::fPlateNumber, G4VXTRenergyLoss::fPlateThick, G4VXTRenergyLoss::fSigma1, G4VXTRenergyLoss::fSigma2, G4cout, G4endl, G4VXTRenergyLoss::GetGasLinearPhotoAbs(), G4VXTRenergyLoss::GetPlateLinearPhotoAbs(), python.hepunit::hbarc, python.hepunit::pi, and G4VProcess::verboseLevel.

72 {
73  G4double result, sum = 0., tmp, cof1, cof2, cofMin, cofPHC, theta2, theta2k;
74  G4double aMa, bMb ,sigma, dump;
75  G4int k, kMax, kMin;
76 
77  aMa = fPlateThick*GetPlateLinearPhotoAbs(energy);
78  bMb = fGasThick*GetGasLinearPhotoAbs(energy);
79  sigma = 0.5*(aMa + bMb);
80  dump = std::exp(-fPlateNumber*sigma);
81  if(verboseLevel > 2) G4cout<<" dump = "<<dump<<G4endl;
82  cofPHC = 4*pi*hbarc;
83  tmp = (fSigma1 - fSigma2)/cofPHC/energy;
84  cof1 = fPlateThick*tmp;
85  cof2 = fGasThick*tmp;
86 
87  cofMin = energy*(fPlateThick + fGasThick)/fGamma/fGamma;
88  cofMin += (fPlateThick*fSigma1 + fGasThick*fSigma2)/energy;
89  cofMin /= cofPHC;
90 
91  theta2 = cofPHC/(energy*(fPlateThick + fGasThick));
92 
93  // if (fGamma < 1200) kMin = G4int(cofMin); // 1200 ?
94  // else kMin = 1;
95 
96 
97  kMin = G4int(cofMin);
98  if (cofMin > kMin) kMin++;
99 
100  // tmp = (fPlateThick + fGasThick)*energy*fMaxThetaTR;
101  // tmp /= cofPHC;
102  // kMax = G4int(tmp);
103  // if(kMax < 0) kMax = 0;
104  // kMax += kMin;
105 
106 
107  kMax = kMin + 49; // 19; // kMin + G4int(tmp);
108 
109  // tmp /= fGamma;
110  // if( G4int(tmp) < kMin ) kMin = G4int(tmp);
111 
112  if(verboseLevel > 2)
113  {
114  G4cout<<cof1<<" "<<cof2<<" "<<cofMin<<G4endl;
115  G4cout<<"kMin = "<<kMin<<"; kMax = "<<kMax<<G4endl;
116  }
117  for( k = kMin; k <= kMax; k++ )
118  {
119  tmp = pi*fPlateThick*(k + cof2)/(fPlateThick + fGasThick);
120  result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
121  // tmp = std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
122  if( k == kMin && kMin == G4int(cofMin) )
123  {
124  sum += 0.5*std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
125  }
126  else
127  {
128  sum += std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
129  }
130  theta2k = std::sqrt(theta2*std::abs(k-cofMin));
131 
132  if(verboseLevel > 2)
133  {
134  // G4cout<<"k = "<<k<<"; sqrt(theta2k) = "<<theta2k<<"; tmp = "<<std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result
135  // <<"; sum = "<<sum<<G4endl;
136  G4cout<<k<<" "<<theta2k<<" "<<std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result
137  <<" "<<sum<<G4endl;
138  }
139  }
140  result = 2*( cof1 + cof2 )*( cof1 + cof2 )*sum/energy;
141  // result *= ( 1 - std::exp(-0.5*fPlateNumber*sigma) )/( 1 - std::exp(-0.5*sigma) );
142  // fPlateNumber;
143  result *= ( 1 - dump + 2*dump*fPlateNumber );
144  /*
145  fEnergy = energy;
146  // G4Integrator<G4VXTRenergyLoss,G4double(G4VXTRenergyLoss::*)(G4double)> integral;
147  G4Integrator<G4TransparentRegXTRadiator,G4double(G4VXTRenergyLoss::*)(G4double)> integral;
148 
149  tmp = integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
150  0.0,0.3*fMaxThetaTR) +
151  integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
152  0.3*fMaxThetaTR,0.6*fMaxThetaTR) +
153  integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
154  0.6*fMaxThetaTR,fMaxThetaTR) ;
155  result += tmp;
156  */
157  return result;
158 }
G4VProcess::verboseLevel
G4int verboseLevel
Definition: G4VProcess.hh:368
G4VXTRenergyLoss::GetPlateLinearPhotoAbs
G4double GetPlateLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1085
G4VXTRenergyLoss::GetGasLinearPhotoAbs
G4double GetGasLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1157
G4int
int G4int
Definition: G4Types.hh:78
energy
double precision function energy(A, Z)
Definition: dpm25nuc6.f:4106
G4cout
G4GLOB_DLL std::ostream G4cout
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76
The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:57:00 for Geant4.10 by   doxygen 1.8.7