#include <G4CoherentPairProductionPhysics.hh>

Inheritance diagram for G4CoherentPairProductionPhysics:

Public Member Functions
	G4CoherentPairProductionPhysics (const G4String &name="Coherent Pair Production Physics")

	~G4CoherentPairProductionPhysics ()=default

void	ConstructParticle () override

void	ConstructProcess () override

void	ActivateIncoherentScattering ()

void	SetNameChannelingModel (const G4String &nameChannelingModel)
	set functions

void	SetNameG4Region (const G4String &nameG4Region)
	set name of G4Region to where the G4ChannelingFastSimModel is active

void	SetLowEnergyLimit (G4double energy)

void	SetHighAngleLimit (G4double angle)

void	SetPPKineticEnergyCut (G4double kineticEnergyCut)

void	SetSamplingPairsNumber (G4int nPairs)

void	SetChargeParticleAngleFactor (G4double chargeParticleAngleFactor)

void	SetNTrajectorySteps (G4int nTrajectorySteps)
	set number of trajectory steps of a single particle (e- or e+)

Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")

	G4VPhysicsConstructor (const G4String &name, G4int physics_type)

virtual	~G4VPhysicsConstructor ()

void	SetPhysicsName (const G4String &="")

const G4String &	GetPhysicsName () const

void	SetPhysicsType (G4int)

G4int	GetPhysicsType () const

G4int	GetInstanceID () const

virtual void	TerminateWorker ()

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

Additional Inherited Members
Static Public Member Functions inherited from G4VPhysicsConstructor
static const G4VPCManager &	GetSubInstanceManager ()

Protected Types inherited from G4VPhysicsConstructor
using	PhysicsBuilder_V = G4VPCData::PhysicsBuilders_V

Protected Member Functions inherited from G4VPhysicsConstructor
G4bool	RegisterProcess (G4VProcess process, G4ParticleDefinition particle)

G4ParticleTable::G4PTblDicIterator *	GetParticleIterator () const

PhysicsBuilder_V	GetBuilders () const

void	AddBuilder (G4PhysicsBuilderInterface *bld)

Protected Attributes inherited from G4VPhysicsConstructor
G4int	verboseLevel = 0

G4String	namePhysics = ""

G4int	typePhysics = 0

G4ParticleTable *	theParticleTable = nullptr

G4int	g4vpcInstanceID = 0

Static Protected Attributes inherited from G4VPhysicsConstructor
static G4RUN_DLL G4VPCManager	subInstanceManager

Detailed Description

Definition at line 36 of file G4CoherentPairProductionPhysics.hh.

Constructor & Destructor Documentation

◆ G4CoherentPairProductionPhysics()

G4CoherentPairProductionPhysics::G4CoherentPairProductionPhysics ( const G4String & name = "Coherent Pair Production Physics" )

Definition at line 38 of file G4CoherentPairProductionPhysics.cc.

                                                                                    :
    G4VPhysicsConstructor(name)
{}

◆ ~G4CoherentPairProductionPhysics()

G4CoherentPairProductionPhysics::~G4CoherentPairProductionPhysics ( )

default

Member Function Documentation

◆ ActivateIncoherentScattering()

void G4CoherentPairProductionPhysics::ActivateIncoherentScattering ( )

inline

activate incoherent scattering (standard gamma conversion should be switched off in physics list)

Definition at line 49 of file G4CoherentPairProductionPhysics.hh.

49{fIncoherentScattering = true;}

◆ ConstructParticle()

void G4CoherentPairProductionPhysics::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 44 of file G4CoherentPairProductionPhysics.cc.

{
    G4Gamma::GammaDefinition();  // Define the gamma particle
}

◆ ConstructProcess()

void G4CoherentPairProductionPhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 51 of file G4CoherentPairProductionPhysics.cc.

{
    //create the gamma process
    G4CoherentPairProduction* gammaProcess = new G4CoherentPairProduction();
 
    if(verboseLevel > 0) {
        G4cout << "G4CoherentPairProductionPhysics::ConstructProcess" << G4endl;
    }
 
    G4RegionStore* regionStore = G4RegionStore::GetInstance();
    G4Region* RegionCh = regionStore->GetRegion(fNameRegion);
 
    //if the region is not found
        if(RegionCh==0)
        {
            G4Exception("GetRegion",// Origin of the exception
                        "001",                   // Unique error code
                        FatalException,          // Terminate the program
                        "Region is not found! The program will terminate.");
        }
        else
        {
            //get channeling model
            G4bool someflag=false;
            G4ChannelingFastSimModel* Channeling =
                static_cast<G4ChannelingFastSimModel*>
                (RegionCh->GetFastSimulationManager()->
                 GetFastSimulationModel(fNameChannelingModel,0,someflag));
 
            //if channeling model is not found
            if(Channeling==0)
            {
                G4Exception("GetFastSimulationModel",// Origin of the exception
                            "001",                   // Unique error code
                            FatalException,          // Terminate the program
                            "Input channeling model is not found! The program will terminate."
                            );
            }
            else
            {
                gammaProcess->Input(Channeling->GetCrystalData());
            }
        }
 
    //set functions
    if(fIncoherentScattering){gammaProcess->ActivateIncoherentScattering();}
    gammaProcess->SetLowEnergyLimit(fLowEnergyLimit);
    gammaProcess->SetHighAngleLimit(fHighAngleLimit);
    gammaProcess->SetPPKineticEnergyCut(fPPKineticEnergyCut);
    gammaProcess->SetSamplingPairsNumber(fNMCPairs);
    gammaProcess->SetChargeParticleAngleFactor(fChargeParticleAngleFactor);
    gammaProcess->SetNTrajectorySteps(fNTrajectorySteps);
 
    //set the name of G4Region in which the model is applicable
    gammaProcess->SetG4RegionName(fNameRegion);
 
    //get process manager for gamma
    G4ProcessManager* pManager = G4Gamma::Gamma()->GetProcessManager();
 
    //register the G4CoherentPairProduction process
    pManager->AddDiscreteProcess(gammaProcess);
}

◆ SetChargeParticleAngleFactor()

void G4CoherentPairProductionPhysics::SetChargeParticleAngleFactor ( G4double chargeParticleAngleFactor )

inline

set the number of particle angles 1/gamma in pair production defining the width of the angular distribution of pair sampling in the Baier-Katkov Integral

Definition at line 72 of file G4CoherentPairProductionPhysics.hh.

73 {fChargeParticleAngleFactor = chargeParticleAngleFactor;}