G4XTRGammaRadModel Class Reference

#include <G4XTRGammaRadModel.hh>

Inheritance diagram for G4XTRGammaRadModel:

Public Member Functions
	G4XTRGammaRadModel (G4LogicalVolume *anEnvelope, G4double, G4double, G4Material *, G4Material *, G4double, G4double, G4int, const G4String &processName="XTRgammaRadiator")
virtual	~G4XTRGammaRadModel ()
G4double	GetStackFactor (G4double energy, G4double gamma, G4double varAngle)

---

Detailed Description

Definition at line 56 of file G4XTRGammaRadModel.hh.

---

Constructor & Destructor Documentation

G4XTRGammaRadModel::G4XTRGammaRadModel	(	G4LogicalVolume *	anEnvelope,
		G4double	,
		G4double	,
		G4Material *	,
		G4Material *	,
		G4double	,
		G4double	,
		G4int	,
		const G4String &	processName = "XTRgammaRadiator"
	)

Definition at line 42 of file G4XTRGammaRadModel.cc.

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

                                                                  :
  G4VXTRenergyLoss(anEnvelope,foilMat,gasMat,a,b,n,processName)
{
  G4cout<<"Gammma distributed X-ray TR radiator model is called"<<G4endl ;

  // Build energy and angular integral spectra of X-ray TR photons from
  // a radiator

  fAlphaPlate = alphaPlate ;
  fAlphaGas   = alphaGas   ;
  G4cout<<"fAlphaPlate = "<<fAlphaPlate<<" ; fAlphaGas = "<<fAlphaGas<<G4endl ;
  fExitFlux = true;
  //  BuildTable() ;
}

G4XTRGammaRadModel::~G4XTRGammaRadModel

(

)

[virtual]

Definition at line 64 of file G4XTRGammaRadModel.cc.

{
  ;
}

---

Member Function Documentation

G4double G4XTRGammaRadModel::GetStackFactor	(	G4double	energy,
		G4double	gamma,
		G4double	varAngle
	)			[virtual]

Reimplemented from G4VXTRenergyLoss.

Definition at line 80 of file G4XTRGammaRadModel.cc.

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

{
  G4double result, Qa, Qb, Q, Za, Zb, Ma, Mb ;
  
  Za = GetPlateFormationZone(energy,gamma,varAngle) ;
  Zb = GetGasFormationZone(energy,gamma,varAngle) ;

  Ma = GetPlateLinearPhotoAbs(energy) ;
  Mb = GetGasLinearPhotoAbs(energy) ;

  Qa = ( 1.0 + fPlateThick*Ma/fAlphaPlate ) ;
  Qa = std::pow(Qa,-fAlphaPlate) ;
  Qb = ( 1.0 + fGasThick*Mb/fAlphaGas ) ;
  Qb = std::pow(Qb,-fAlphaGas) ;
  Q  = Qa*Qb ;

  G4complex Ca(1.0+0.5*fPlateThick*Ma/fAlphaPlate,fPlateThick/Za/fAlphaPlate) ; 
  G4complex Cb(1.0+0.5*fGasThick*Mb/fAlphaGas,fGasThick/Zb/fAlphaGas) ; 

  G4complex Ha = std::pow(Ca,-fAlphaPlate) ;  
  G4complex Hb = std::pow(Cb,-fAlphaGas) ;
  G4complex H  = Ha*Hb ;

  G4complex F1 = ( 0.5*(1+Qa)*(1.0+H) - Ha - Qa*Hb )/(1.0-H) ;

  G4complex F2 = (1.0-Ha)*(Qa-Ha)*Hb/(1.0-H)/(Q-H) ;

  F2          *= std::pow(Q,G4double(fPlateNumber)) - std::pow(H,fPlateNumber) ;

  result      = ( 1 - std::pow(Q,G4double(fPlateNumber)) )/( 1 - Q ) ;

  G4complex stack  = result*F1;
            stack += F2;
            stack *= 2.0*OneInterfaceXTRdEdx(energy,gamma,varAngle);

            result = std::real(stack);

            // result     *= 2.0*std::real(F1);
            // result     += 2.0*std::real(F2);

  return      result ;
}

---

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

• G4XTRGammaRadModel.hh
• G4XTRGammaRadModel.cc

---