#include <G4LowEnergyBremsstrahlung.hh>

Inheritance diagram for G4LowEnergyBremsstrahlung:

Collaboration diagram for G4LowEnergyBremsstrahlung:

Public Member Functions
	G4LowEnergyBremsstrahlung (const G4String &processName="LowEnBrem")

	~G4LowEnergyBremsstrahlung ()

G4bool	IsApplicable (const G4ParticleDefinition &)

void	BuildPhysicsTable (const G4ParticleDefinition &particleType)

G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &step)

void	SetCutForLowEnSecPhotons (G4double cut)

void	SetAngularGenerator (G4RDVBremAngularDistribution *distribution)

void	SetAngularGenerator (const G4String &name)

void	PrintInfoDefinition ()

Public Member Functions inherited from G4eLowEnergyLoss
	G4eLowEnergyLoss (const G4String &)

	~G4eLowEnergyLoss ()

void	BuildDEDXTable (const G4ParticleDefinition &aParticleType)

G4double	GetContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)

G4VParticleChange *	AlongStepDoIt (const G4Track &track, const G4Step &Step)

void	ActivateFluorescence (G4bool val)

G4bool	Fluorescence () const

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

	G4RDVeLowEnergyLoss (G4RDVeLowEnergyLoss &)

virtual	~G4RDVeLowEnergyLoss ()

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

	G4VContinuousDiscreteProcess (G4VContinuousDiscreteProcess &)

virtual	~G4VContinuousDiscreteProcess ()

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety, G4GPILSelection *selection)

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)

virtual G4VParticleChange *	AtRestDoIt (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 Member Functions
G4double	GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

Protected Member Functions inherited from G4eLowEnergyLoss
virtual std::vector < G4DynamicParticle * > *	DeexciteAtom (const G4MaterialCutsCouple *, G4double, G4double)

Protected Member Functions inherited from G4RDVeLowEnergyLoss
G4double	GetLossWithFluct (const G4DynamicParticle aParticle, const G4MaterialCutsCouple couple, G4double MeanLoss, G4double step)

Protected Member Functions inherited from G4VContinuousDiscreteProcess
void	SetGPILSelection (G4GPILSelection selection)

G4GPILSelection	GetGPILSelection () const

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

void	ClearNumberOfInteractionLengthLeft ()

Additional Inherited Members
Static Public Member Functions inherited from G4eLowEnergyLoss
static void	SetNbOfProcesses (G4int nb)

static void	PlusNbOfProcesses ()

static void	MinusNbOfProcesses ()

static G4int	GetNbOfProcesses ()

static void	SetLowerBoundEloss (G4double val)

static void	SetUpperBoundEloss (G4double val)

static void	SetNbinEloss (G4int nb)

static G4double	GetLowerBoundEloss ()

static G4double	GetUpperBoundEloss ()

static G4int	GetNbinEloss ()

Static Public Member Functions inherited from G4RDVeLowEnergyLoss
static void	SetRndmStep (G4bool value)

static void	SetEnlossFluc (G4bool value)

static void	SetStepFunction (G4double c1, G4double c2)

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

Static Protected Member Functions inherited from G4RDVeLowEnergyLoss
static G4PhysicsTable *	BuildRangeTable (G4PhysicsTable theDEDXTable, G4PhysicsTable theRangeTable, G4double Tmin, G4double Tmax, G4int nbin)

static G4PhysicsTable *	BuildLabTimeTable (G4PhysicsTable theDEDXTable, G4PhysicsTable theLabTimeTable, G4double Tmin, G4double Tmax, G4int nbin)

static G4PhysicsTable *	BuildProperTimeTable (G4PhysicsTable theDEDXTable, G4PhysicsTable ProperTimeTable, G4double Tmin, G4double Tmax, G4int nbin)

static G4PhysicsTable *	BuildRangeCoeffATable (G4PhysicsTable theRangeTable, G4PhysicsTable theCoeffATable, G4double Tmin, G4double Tmax, G4int nbin)

static G4PhysicsTable *	BuildRangeCoeffBTable (G4PhysicsTable theRangeTable, G4PhysicsTable theCoeffBTable, G4double Tmin, G4double Tmax, G4int nbin)

static G4PhysicsTable *	BuildRangeCoeffCTable (G4PhysicsTable theRangeTable, G4PhysicsTable theCoeffCTable, G4double Tmin, G4double Tmax, G4int nbin)

static G4PhysicsTable *	BuildInverseRangeTable (G4PhysicsTable theRangeTable, G4PhysicsTable theRangeCoeffATable, G4PhysicsTable theRangeCoeffBTable, G4PhysicsTable theRangeCoeffCTable, G4PhysicsTable *theInverseRangeTable, G4double Tmin, G4double Tmax, G4int nbin)

Protected Attributes inherited from G4eLowEnergyLoss
G4PhysicsTable *	theLossTable

G4double	MinKineticEnergy

G4double	Charge

G4double	lastCharge

Protected Attributes inherited from G4RDVeLowEnergyLoss
const G4Material *	lastMaterial

G4int	imat

G4double	f1Fluct

G4double	f2Fluct

G4double	e1Fluct

G4double	e2Fluct

G4double	rateFluct

G4double	ipotFluct

G4double	e1LogFluct

G4double	e2LogFluct

G4double	ipotLogFluct

const G4int	nmaxCont1

const G4int	nmaxCont2

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 inherited from G4eLowEnergyLoss
static G4PhysicsTable *	theDEDXElectronTable = 0

static G4PhysicsTable *	theDEDXPositronTable = 0

static G4PhysicsTable *	theRangeElectronTable = 0

static G4PhysicsTable *	theRangePositronTable = 0

static G4PhysicsTable *	theInverseRangeElectronTable = 0

static G4PhysicsTable *	theInverseRangePositronTable = 0

static G4PhysicsTable *	theLabTimeElectronTable = 0

static G4PhysicsTable *	theLabTimePositronTable = 0

static G4PhysicsTable *	theProperTimeElectronTable = 0

static G4PhysicsTable *	theProperTimePositronTable = 0

static G4int	NbOfProcesses = 2

static G4int	CounterOfElectronProcess = 0

static G4int	CounterOfPositronProcess = 0

static G4PhysicsTable **	RecorderOfElectronProcess

static G4PhysicsTable **	RecorderOfPositronProcess

Static Protected Attributes inherited from G4RDVeLowEnergyLoss
static G4double	ParticleMass

static G4double	taulow

static G4double	tauhigh

static G4double	ltaulow

static G4double	ltauhigh

static G4double	dRoverRange = 20*perCent

static G4double	finalRange = 200*micrometer

static G4double	c1lim = dRoverRange

static G4double	c2lim = 2.(1.-dRoverRange)finalRange

static G4double	c3lim = -(1.-dRoverRange)finalRangefinalRange

static G4bool	rndmStepFlag = false

static G4bool	EnlossFlucFlag = true

Detailed Description

Definition at line 70 of file G4LowEnergyBremsstrahlung.hh.

Constructor & Destructor Documentation

G4LowEnergyBremsstrahlung::G4LowEnergyBremsstrahlung ( const G4String & processName = "LowEnBrem" )

Definition at line 85 of file G4LowEnergyBremsstrahlung.cc.

   : G4eLowEnergyLoss(nam),
   crossSectionHandler(0),
   theMeanFreePath(0),
   energySpectrum(0)
 {
   cutForPhotons = 0.;
   verboseLevel = 0;
   generatorName = "tsai";
   angularDistribution = new G4RDModifiedTsai("TsaiGenerator"); // default generator
 //  angularDistribution->PrintGeneratorInformation();
   TsaiAngularDistribution = new G4RDModifiedTsai("TsaiGenerator");
 }

G4LowEnergyBremsstrahlung::~G4LowEnergyBremsstrahlung ( )

Definition at line 116 of file G4LowEnergyBremsstrahlung.cc.

 {
   if(crossSectionHandler) delete crossSectionHandler;
   if(energySpectrum) delete energySpectrum;
   if(theMeanFreePath) delete theMeanFreePath;
   delete angularDistribution;
   delete TsaiAngularDistribution;
   energyBins.clear();
 }

Member Function Documentation

void G4LowEnergyBremsstrahlung::BuildPhysicsTable ( const G4ParticleDefinition & particleType )

virtual

Reimplemented from G4VProcess.

Definition at line 127 of file G4LowEnergyBremsstrahlung.cc.

 {
   if(verboseLevel > 0) {
     G4cout << "G4LowEnergyBremsstrahlung::BuildPhysicsTable start"
            << G4endl;
       }
 
   cutForSecondaryPhotons.clear();
 
   // Create and fill BremsstrahlungParameters once
   if( energySpectrum != 0 ) delete energySpectrum;
   energyBins.clear();
   for(size_t i=0; i<15; i++) {
     G4double x = 0.1*((G4double)i);
     if(i == 0)  x = 0.01;
     if(i == 10) x = 0.95;
     if(i == 11) x = 0.97;
     if(i == 12) x = 0.99;
     if(i == 13) x = 0.995;
     if(i == 14) x = 1.0;
     energyBins.push_back(x);
   }
   const G4String dataName("/brem/br-sp.dat");
   energySpectrum = new G4RDeBremsstrahlungSpectrum(energyBins,dataName);
 
   if(verboseLevel > 0) {
     G4cout << "G4LowEnergyBremsstrahlungSpectrum is initialized"
            << G4endl;
       }
 
   // Create and fill G4RDCrossSectionHandler once
 
   if( crossSectionHandler != 0 ) delete crossSectionHandler;
   G4RDVDataSetAlgorithm* interpolation = new G4RDLogLogInterpolation();
   G4double lowKineticEnergy  = GetLowerBoundEloss();
   G4double highKineticEnergy = GetUpperBoundEloss();
   G4int    totBin = GetNbinEloss();
   crossSectionHandler = new G4RDBremsstrahlungCrossSectionHandler(energySpectrum, interpolation);
   crossSectionHandler->Initialise(0,lowKineticEnergy, highKineticEnergy, totBin);
   crossSectionHandler->LoadShellData("brem/br-cs-");
 
   if (verboseLevel > 0) {
     G4cout << GetProcessName()
            << " is created; Cross section data: "
            << G4endl;
     crossSectionHandler->PrintData();
     G4cout << "Parameters: "
            << G4endl;
     energySpectrum->PrintData();
   }
 
   // Build loss table for Bremsstrahlung
 
   BuildLossTable(aParticleType);
 
   if(verboseLevel > 0) {
     G4cout << "The loss table is built"
            << G4endl;
       }
 
   if (&aParticleType==G4Electron::Electron()) {
 
     RecorderOfElectronProcess[CounterOfElectronProcess] = (*this).theLossTable;
     CounterOfElectronProcess++;
     PrintInfoDefinition();  
 
   } else {
 
     RecorderOfPositronProcess[CounterOfPositronProcess] = (*this).theLossTable;
     CounterOfPositronProcess++;
   }
 
   // Build mean free path data using cut values
 
   if( theMeanFreePath != 0 ) delete theMeanFreePath;
   theMeanFreePath = crossSectionHandler->
                     BuildMeanFreePathForMaterials(&cutForSecondaryPhotons);
 
   if(verboseLevel > 0) {
     G4cout << "The MeanFreePath table is built"
            << G4endl;
       }
 
   // Build common DEDX table for all ionisation processes
 
   BuildDEDXTable(aParticleType);
 
   if(verboseLevel > 0) {
     G4cout << "G4LowEnergyBremsstrahlung::BuildPhysicsTable end"
            << G4endl;
       }
  
 }

Here is the call graph for this function:

G4double G4LowEnergyBremsstrahlung::GetMeanFreePath	(	const G4Track &	track,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)

protectedvirtual

Implements G4eLowEnergyLoss.

Definition at line 387 of file G4LowEnergyBremsstrahlung.cc.

 {
   *cond = NotForced;
   G4int index = (track.GetMaterialCutsCouple())->GetIndex();
   const G4RDVEMDataSet* data = theMeanFreePath->GetComponent(index);
   G4double meanFreePath = data->FindValue(track.GetKineticEnergy());
   return meanFreePath;
 }

Here is the call graph for this function:

G4bool G4LowEnergyBremsstrahlung::IsApplicable ( const G4ParticleDefinition & particle )

virtual

Reimplemented from G4eLowEnergyLoss.

Definition at line 381 of file G4LowEnergyBremsstrahlung.cc.

 {
   return (  (&particle == G4Electron::Electron())  );
 }

Here is the call graph for this function:

G4VParticleChange * G4LowEnergyBremsstrahlung::PostStepDoIt	(	const G4Track &	track,
		const G4Step &	step
	)

virtual

Implements G4eLowEnergyLoss.

Definition at line 302 of file G4LowEnergyBremsstrahlung.cc.

 {
   aParticleChange.Initialize(track);
 
   const G4MaterialCutsCouple* couple = track.GetMaterialCutsCouple();
   G4double kineticEnergy = track.GetKineticEnergy();
   G4int index = couple->GetIndex();
   G4double tCut = cutForSecondaryPhotons[index];
 
   // Control limits
   if(tCut >= kineticEnergy)
      return G4VContinuousDiscreteProcess::PostStepDoIt(track, step);
 
   G4int Z = crossSectionHandler->SelectRandomAtom(couple, kineticEnergy);
 
   G4double tGamma = energySpectrum->SampleEnergy(Z, tCut, kineticEnergy, kineticEnergy);
   G4double totalEnergy = kineticEnergy + electron_mass_c2;
   G4double finalEnergy = kineticEnergy - tGamma; // electron/positron final energy  
   G4double theta = 0;
 
   if((kineticEnergy < 1*MeV && kineticEnergy > 1*keV && generatorName == "2bn")){
       theta = angularDistribution->PolarAngle(kineticEnergy,finalEnergy,Z);
   }else{
       theta = TsaiAngularDistribution->PolarAngle(kineticEnergy,finalEnergy,Z);
   }
 
   G4double phi   = twopi * G4UniformRand();
   G4double dirZ  = std::cos(theta);
   G4double sinTheta  = std::sqrt(1. - dirZ*dirZ);
   G4double dirX  = sinTheta*std::cos(phi);
   G4double dirY  = sinTheta*std::sin(phi);
 
   G4ThreeVector gammaDirection (dirX, dirY, dirZ);
   G4ThreeVector electronDirection = track.GetMomentumDirection();
 
   //
   // Update the incident particle 
   //
   gammaDirection.rotateUz(electronDirection);   
     
   // Kinematic problem
   if (finalEnergy < 0.) {
     tGamma += finalEnergy;
     finalEnergy = 0.0;
   }
 
   G4double momentum = std::sqrt((totalEnergy + electron_mass_c2)*kineticEnergy);
 
   G4double finalX = momentum*electronDirection.x() - tGamma*gammaDirection.x();
   G4double finalY = momentum*electronDirection.y() - tGamma*gammaDirection.y();
   G4double finalZ = momentum*electronDirection.z() - tGamma*gammaDirection.z();
       
   aParticleChange.SetNumberOfSecondaries(1);
   G4double norm = 1./std::sqrt(finalX*finalX + finalY*finalY + finalZ*finalZ);
   aParticleChange.ProposeMomentumDirection(finalX*norm, finalY*norm, finalZ*norm);
   aParticleChange.ProposeEnergy( finalEnergy );
 
   // create G4DynamicParticle object for the gamma 
   G4DynamicParticle* aGamma= new G4DynamicParticle (G4Gamma::Gamma(),
                             gammaDirection, tGamma);
   aParticleChange.AddSecondary(aGamma);
 
   return G4VContinuousDiscreteProcess::PostStepDoIt(track, step);
 }

Here is the call graph for this function:

void G4LowEnergyBremsstrahlung::PrintInfoDefinition ( )

Definition at line 369 of file G4LowEnergyBremsstrahlung.cc.

 {
   G4String comments = "Total cross sections from EEDL database.";
   comments += "\n      Gamma energy sampled from a parameterised formula.";
   comments += "\n      Implementation of the continuous dE/dx part.";  
   comments += "\n      At present it can be used for electrons ";
   comments += "in the energy range [250eV,100GeV].";
   comments += "\n      The process must work with G4LowEnergyIonisation.";
   
   G4cout << G4endl << GetProcessName() << ":  " << comments << G4endl;
 }         

Here is the call graph for this function:

Here is the caller graph for this function:

void G4LowEnergyBremsstrahlung::SetAngularGenerator ( G4RDVBremAngularDistribution * distribution )

Definition at line 403 of file G4LowEnergyBremsstrahlung.cc.

 {
   angularDistribution = distribution;
   angularDistribution->PrintGeneratorInformation();
 }

Here is the call graph for this function:

void G4LowEnergyBremsstrahlung::SetAngularGenerator ( const G4String & name )

Definition at line 409 of file G4LowEnergyBremsstrahlung.cc.

 {
   if (name == "tsai") 
     {
       delete angularDistribution;
       angularDistribution = new G4RDModifiedTsai("TsaiGenerator");
       generatorName = name;
     }
   else if (name == "2bn")
     {
       delete angularDistribution;
       angularDistribution = new G4RDGenerator2BN("2BNGenerator");
       generatorName = name;
     }
   else if (name == "2bs")
     {
        delete angularDistribution;
        angularDistribution = new G4RDGenerator2BS("2BSGenerator");
        generatorName = name;
     }
   else
     {
       G4Exception("G4LowEnergyBremsstrahlung::SetAngularGenerator()",
                   "InvalidSetup", FatalException, "Generator does not exist!");
     }
 
   angularDistribution->PrintGeneratorInformation();
 }

Here is the call graph for this function:

void G4LowEnergyBremsstrahlung::SetCutForLowEnSecPhotons ( G4double cut )

Definition at line 398 of file G4LowEnergyBremsstrahlung.cc.

 {
   cutForPhotons = cut;
 }

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

source/geant4.10.03.p03/examples/advanced/eRosita/physics/include/G4LowEnergyBremsstrahlung.hh
source/geant4.10.03.p03/examples/advanced/eRosita/physics/src/G4LowEnergyBremsstrahlung.cc

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation