G4URRNeutrons Class Reference

#include <G4URRNeutrons.hh>

Inheritance diagram for G4URRNeutrons:

Public Member Functions
	G4URRNeutrons (G4int ver=1)
virtual	~G4URRNeutrons ()
void	ConstructProcess () override
Public Member Functions inherited from G4VHadronPhysics
	G4VHadronPhysics (const G4String &name="hInelastic", G4int verbose=0)
virtual	~G4VHadronPhysics ()
void	ConstructParticle () override
G4HadronicInteraction *	BuildModel (G4VHadronModelBuilder *, G4double emin, G4double emax)
G4HadronicInteraction *	NewModel (G4HadronicInteraction *, G4double emin, G4double emax)
	G4VHadronPhysics (G4VHadronPhysics &)=delete
G4VHadronPhysics &	operator= (const G4VHadronPhysics &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 55 of file G4URRNeutrons.hh.

Constructor & Destructor Documentation

G4URRNeutrons()

G4URRNeutrons::G4URRNeutrons ( G4int ver = 1 )

explicit

Definition at line 68 of file G4URRNeutrons.cc.

: G4VHadronPhysics( "URRNeutrons", ver ) {}

~G4URRNeutrons()

G4URRNeutrons::~G4URRNeutrons ( )

virtual

Definition at line 71 of file G4URRNeutrons.cc.

{}

Member Function Documentation

ConstructProcess()

void G4URRNeutrons::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 74 of file G4URRNeutrons.cc.

                                     {
  // Find elastic, capture, fission and inelastic processes of neutron
  // (from the physics list on which this constructor is applied on top);
  // then look at their hadronic final-state models in order to disable those that
  // would otherwise fully overlap with the URR hadronic final-state models.
  // Note that for disabling a hadronic model - not being defined the "DeRegister"
  // method - it is enough to set to 0.0 the max energy of the model. 
  
  if ( G4HadronicParameters::Instance()->GetVerboseLevel() > 1 ) {
    G4cout << "### " << GetPhysicsName() << " Construct Processes " << G4endl;
  }
  G4Neutron* part = G4Neutron::Neutron();
 
  // Elastic (including, eventually, thermal scattering)
  G4HadronicProcess* elasticProcess = G4PhysListUtil::FindElasticProcess( part );
  if ( elasticProcess == nullptr ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron elastic treatment: "
       << "NOT found elastic process => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int niElastic = static_cast< G4int >( (elasticProcess->GetHadronicInteractionList()).size() );
  if ( niElastic < 1 ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron elastic treatment: "
       << "NOT found any elastic model => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int indexHPelastic = -1;
  G4int indexHPthermalScattering = -1;
  for ( G4int index = 0; index < niElastic; ++index ) {
    if ( (elasticProcess->GetHadronicInteractionList())[index]->GetModelName() == "NeutronHPElastic" ) {
      indexHPelastic = index;
    } else if ( (elasticProcess->GetHadronicInteractionList())[index]->GetModelName() == "NeutronHPThermalScattering" ) {
      indexHPthermalScattering = index;
    }
  }
  if ( indexHPelastic >= 0 ) {
    (elasticProcess->GetHadronicInteractionList())[indexHPelastic]->SetMaxEnergy( 0.0 );  // Disabled
    G4cout << G4endl << " G4URRNeutrons::ConstructProcess() : found NeutronHPElastic => Disabled !" << G4endl;
  } else {
    G4cout << "### " << GetPhysicsName()
       << " WARNING: NOT found NeutronHPElastic => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4bool isThermalScatteringOn = false;
  if ( indexHPthermalScattering > 0 ) isThermalScatteringOn = true;
  elasticProcess->AddDataSet( new G4ParticleHPElasticDataPT );
  elasticProcess->RegisterMe( new G4ParticleHPElasticURR( isThermalScatteringOn ) );
 
  // Capture
  G4HadronicProcess* captureProcess = G4PhysListUtil::FindCaptureProcess( part );
  if ( captureProcess == nullptr ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron capture treatment: "
       << "NOT found capture process => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int niCapture = static_cast< G4int >( (captureProcess->GetHadronicInteractionList()).size() );
  if ( niCapture < 1 ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron capture treatment: "
       << "NOT found any capture model => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int indexHPcapture = -1;
  for ( G4int index = 0; index < niCapture; ++index ) {
    // For QGSP_BERT_HP, since G4 11.2, the neutron capture model is called either "nRadCaptureHP"
    // or "nuDEX_neutronCapture" (the latter when NuDEX is used).
    G4String nameNeutronCaptureModel = (captureProcess->GetHadronicInteractionList())[index]->GetModelName();
    if ( nameNeutronCaptureModel == "NeutronHPCapture"  ||
     nameNeutronCaptureModel == "nRadCaptureHP"     ||
     nameNeutronCaptureModel == "nuDEX_neutronCapture" ) {
      indexHPcapture = index;
    }
  }
  if ( indexHPcapture >= 0 ) {
    (captureProcess->GetHadronicInteractionList())[indexHPcapture]->SetMaxEnergy( 0.0 );  // Disabled
    G4cout << G4endl << " G4URRNeutrons::ConstructProcess() : found "
       << (captureProcess->GetHadronicInteractionList())[indexHPcapture]->GetModelName()
           << " => Disabled !" << G4endl;
  } else {
    G4cout << "### " << GetPhysicsName()
       << " WARNING: NOT found any expected neutron capture model => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  captureProcess->AddDataSet( new G4ParticleHPCaptureDataPT );
  captureProcess->RegisterMe( new G4ParticleHPCaptureURR );
  
  // Fission
  G4HadronicProcess* fissionProcess  = G4PhysListUtil::FindFissionProcess( part );
  if ( fissionProcess == nullptr ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron fission treatment: "
       << "NOT found fission process => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int niFission = static_cast< G4int >( (fissionProcess->GetHadronicInteractionList()).size() );
  if ( niFission < 1 ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron fission treatment: "
       << "NOT found any fission model => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int indexHPfission = -1;
  for ( G4int index = 0; index < niFission; ++index ) {
    // For QGSP_BERT_HP, since G4 11.2, the neutron fission model is called "nFissionVI"
    G4String nameNeutronFissionModel = (fissionProcess->GetHadronicInteractionList())[index]->GetModelName();
    if ( nameNeutronFissionModel == "NeutronHPFission"  ||  nameNeutronFissionModel == "nFissionVI" ) {
      indexHPfission = index;
    }
  }
  if ( indexHPfission >= 0 ) {
    (fissionProcess->GetHadronicInteractionList())[indexHPfission]->SetMaxEnergy( 0.0 );  // Disabled
    G4cout << G4endl << " G4URRNeutrons::ConstructProcess() : found "
           << (fissionProcess->GetHadronicInteractionList())[indexHPfission]->GetModelName()  // Disabled
           << " => Disabled !" << G4endl;
  } else {
    G4cout << "### " << GetPhysicsName()
       << " WARNING: NOT found any expected neutron fission model => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  fissionProcess->RegisterMe( new G4ParticleHPFissionURR );
  fissionProcess->AddDataSet( new G4ParticleHPFissionDataPT );
 
  // Inelastic
  G4HadronicProcess* inelasticProcess = G4PhysListUtil::FindInelasticProcess( part );
  if ( inelasticProcess == nullptr ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron inelastic treatment: "
       << "NOT found inelastic process => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int niInelastic = static_cast< G4int >( (inelasticProcess->GetHadronicInteractionList()).size() );
  if ( niInelastic < 1 ) {
    G4cout << "### " << GetPhysicsName() << " WARNING: Fail to add URR neutron inelastic treatment: "
       << "NOT found any inelastic model => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  G4int indexHPinelastic = -1;
  for ( G4int index = 0; index < niInelastic; ++index ) {
    if ( (inelasticProcess->GetHadronicInteractionList())[index]->GetModelName() == "NeutronHPInelastic" ) {
      indexHPinelastic = index;
    }
  }
  if ( indexHPinelastic >= 0 ) {
    (inelasticProcess->GetHadronicInteractionList())[indexHPinelastic]->SetMaxEnergy( 0.0 );  // Disabled
    G4cout << G4endl << " G4URRNeutrons::ConstructProcess() : found NeutronHPInelastic => Disabled !" << G4endl;
  } else {
    G4cout << "### " << GetPhysicsName()
       << " WARNING: NOT found NeutronHPInelastic => G4URRNeutrons returns without doing anything !" << G4endl;
    return;
  }
  inelasticProcess->AddDataSet( new G4ParticleHPInelasticDataPT );
  inelasticProcess->RegisterMe( new G4ParticleHPInelasticURR );
}

---

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

• G4URRNeutrons.hh
• G4URRNeutrons.cc