G4HadronPhysicsLEND Class Reference

#include <G4HadronPhysicsLEND.hh>

Inheritance diagram for G4HadronPhysicsLEND:

Public Member Functions
	G4HadronPhysicsLEND (G4int verbose=1, const G4String &eval="")
	~G4HadronPhysicsLEND () override=default
void	ConstructProcess () override
void	ConstructParticle () override
	G4HadronPhysicsLEND (G4HadronPhysicsLEND &)=delete
G4HadronPhysicsLEND &	operator= (const G4HadronPhysicsLEND &right)=delete
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 44 of file G4HadronPhysicsLEND.hh.

Constructor & Destructor Documentation

G4HadronPhysicsLEND() [1/2]

G4HadronPhysicsLEND::G4HadronPhysicsLEND	(	G4int	verbose = 1,
		const G4String &	eval = "" )

Definition at line 56 of file G4HadronPhysicsLEND.cc.

 : G4VPhysicsConstructor("hadron inelastic LEND")
{
    verbose = verb;
    evaluation = eva;
}

Referenced by G4HadronPhysicsLEND(), and operator=().

~G4HadronPhysicsLEND()

G4HadronPhysicsLEND::~G4HadronPhysicsLEND ( )

overridedefault

G4HadronPhysicsLEND() [2/2]

G4HadronPhysicsLEND::G4HadronPhysicsLEND ( G4HadronPhysicsLEND & )

delete

Member Function Documentation

ConstructParticle()

void G4HadronPhysicsLEND::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 126 of file G4HadronPhysicsLEND.cc.

{
    G4Neutron::Neutron();
    G4Gamma::Gamma();
}

ConstructProcess()

void G4HadronPhysicsLEND::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 63 of file G4HadronPhysicsLEND.cc.

{
    //....neutron-induced reactions
    
    //....create neutron processes if they do not exist
    if( G4PhysListUtil::FindInelasticProcess(G4Neutron::Neutron()) == nullptr ){
        auto neutron_processes = new G4NeutronBuilder( true ); // Fission on
        AddBuilder(neutron_processes);
        neutron_processes->Build();
    }
 
    //....get pointers to neutron processes so can add LEND model
    auto* neutron_inelastic = (G4HadronInelasticProcess*)G4PhysListUtil::FindInelasticProcess( G4Neutron::Neutron() );
    auto* neutron_capture   = (G4NeutronCaptureProcess*)G4PhysListUtil::FindCaptureProcess( G4Neutron::Neutron() );
    auto* neutron_fission   = (G4NeutronFissionProcess*)G4PhysListUtil::FindFissionProcess( G4Neutron::Neutron() );
        
    auto* neutronLENDBuilder = new G4NeutronLENDBuilder(evaluation);
    neutronLENDBuilder->Build(neutron_inelastic);
    neutronLENDBuilder->Build(neutron_capture);
    neutronLENDBuilder->Build(neutron_fission);
 
    //....photon-induced reactions
    
    //....use gamma hadronic process if it exists
    //          check Inelastic process using G4PhysListUtil
    //          then check for GammaNuclear in G4GammaGeneralProcess
    auto* gamma_inelastic = (G4HadronInelasticProcess*)G4PhysListUtil::FindInelasticProcess( G4Gamma::Gamma() );
    
    if( gamma_inelastic == nullptr){
        G4LossTableManager* emManager = G4LossTableManager::Instance();
        auto gproc = dynamic_cast<G4GammaGeneralProcess*>(emManager->GetGammaGeneralProcess());
        if (gproc != nullptr) {
            G4HadronicProcess* gnuc = gproc->GetGammaNuclear();
            if (gnuc != nullptr )
                gamma_inelastic = (G4HadronInelasticProcess*)gnuc;
        }
    }
    //....did not find existing process, so create photonuclear process
    if( gamma_inelastic == nullptr){
        //std::cout << "DMW: HadPhysLEND = gamma inelastic pointer not found\n";
        gamma_inelastic = new G4HadronInelasticProcess( "photonNuclear", G4Gamma::Gamma() );
        G4PhysicsListHelper* plHelper = G4PhysicsListHelper::GetPhysicsListHelper();
        plHelper->RegisterProcess(gamma_inelastic, G4Gamma::Gamma());
 
        //....need these XS or LENDonly crashes
        auto xsreg = G4CrossSectionDataSetRegistry::Instance();
        G4VCrossSectionDataSet* xs = nullptr;
        xs = xsreg->GetCrossSectionDataSet("GammaNuclearXS");
        if(nullptr == xs) xs = new G4GammaNuclearXS();
        gamma_inelastic->AddDataSet(xs);
    }
 
    //....add LEND photonuclear models
    //auto* theGammaReactionLowE = new G4LENDorBERTModel( G4Gamma::Gamma() );
    auto* theGammaReactionLowE = new G4LENDCombinedModel( G4Gamma::Gamma() ); // use combined LEND models
    theGammaReactionLowE->SetMaxEnergy(maxLEND_Energy);
    theGammaReactionLowE->DumpLENDTargetInfo(true);
    gamma_inelastic->RegisterMe(theGammaReactionLowE);
    
    G4LENDCombinedCrossSection* theGammaCrossSectionLowE = new G4LENDCombinedCrossSection( G4Gamma::Gamma() );
    gamma_inelastic->AddDataSet(theGammaCrossSectionLowE);
}

operator=()

G4HadronPhysicsLEND & G4HadronPhysicsLEND::operator= ( const G4HadronPhysicsLEND & right )

delete

---

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

• G4HadronPhysicsLEND.hh
• G4HadronPhysicsLEND.cc