#include <G4BertiniElectroNuclearBuilder.hh>

Inheritance diagram for G4BertiniElectroNuclearBuilder:

Public Member Functions
	G4BertiniElectroNuclearBuilder (G4bool eNucl)

virtual	~G4BertiniElectroNuclearBuilder ()=default

virtual void	Build ()

G4BertiniElectroNuclearBuilder &	operator= (const G4BertiniElectroNuclearBuilder &right)=delete

	G4BertiniElectroNuclearBuilder (const G4BertiniElectroNuclearBuilder &)=delete

Protected Attributes
G4HadronInelasticProcess *	thePhotoNuclearProcess {nullptr}

G4CascadeInterface *	theGammaReaction {nullptr}

G4bool	eActivated

Detailed Description

Definition at line 33 of file G4BertiniElectroNuclearBuilder.hh.

Constructor & Destructor Documentation

◆ G4BertiniElectroNuclearBuilder() [1/2]

G4BertiniElectroNuclearBuilder::G4BertiniElectroNuclearBuilder ( G4bool eNucl )

explicit

Definition at line 71 of file G4BertiniElectroNuclearBuilder.cc.

72 : eActivated(eNucl)

73{}

G4BertiniElectroNuclearBuilder::eActivated

G4bool eActivated

Definition G4BertiniElectroNuclearBuilder.hh:50

◆ ~G4BertiniElectroNuclearBuilder()

virtual G4BertiniElectroNuclearBuilder::~G4BertiniElectroNuclearBuilder ( )

virtualdefault

◆ G4BertiniElectroNuclearBuilder() [2/2]

G4BertiniElectroNuclearBuilder::G4BertiniElectroNuclearBuilder ( const G4BertiniElectroNuclearBuilder & )

delete

Member Function Documentation

◆ Build()

void G4BertiniElectroNuclearBuilder::Build ( )

virtual

Reimplemented in G4LENDBertiniGammaElectroNuclearBuilder.

Definition at line 75 of file G4BertiniElectroNuclearBuilder.cc.

{
  // gamma
  thePhotoNuclearProcess = new G4HadronInelasticProcess( "photonNuclear", G4Gamma::Definition() );
  thePhotoNuclearProcess->AddDataSet( new G4PhotoNuclearCrossSection() );
  theGammaReaction = new G4CascadeInterface();
 
  auto theModel = new G4TheoFSGenerator;
 
  auto theStringModel = new G4QGSModel< G4GammaParticipants >;
  auto theStringDecay = new G4ExcitedStringDecay( new G4QGSMFragmentation() );
  theStringModel->SetFragmentationModel(theStringDecay);
 
  auto theCascade = new G4GeneratorPrecompoundInterface();
 
  theModel->SetTransport(theCascade);
  theModel->SetHighEnergyGenerator(theStringModel);
 
  G4ProcessManager* aProcMan = nullptr;
 
  theGammaReaction->SetMaxEnergy(3.5*GeV);
  thePhotoNuclearProcess->RegisterMe(theGammaReaction);
  theModel->SetMinEnergy(3.*GeV);
  theModel->SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
  thePhotoNuclearProcess->RegisterMe(theModel);
  
  G4GammaGeneralProcess* sp = 
    dynamic_cast<G4GammaGeneralProcess*>(G4LossTableManager::Instance()->GetGammaGeneralProcess());
  if ( nullptr != sp ) {
    sp->AddHadProcess(thePhotoNuclearProcess);
  } else {
    aProcMan = G4Gamma::Gamma()->GetProcessManager();
    aProcMan->AddDiscreteProcess(thePhotoNuclearProcess);
  }
 
  // e+, e-
  if (eActivated) {  
 
    auto theElectronNuclearProcess = new G4ElectronNuclearProcess();
    auto thePositronNuclearProcess = new G4PositronNuclearProcess();
    auto theElectroReaction = new G4ElectroVDNuclearModel();
 
    aProcMan = G4Electron::Electron()->GetProcessManager();
    theElectronNuclearProcess->RegisterMe(theElectroReaction);
    aProcMan->AddDiscreteProcess(theElectronNuclearProcess);
  
    aProcMan = G4Positron::Positron()->GetProcessManager();
    thePositronNuclearProcess->RegisterMe(theElectroReaction);
    aProcMan->AddDiscreteProcess(thePositronNuclearProcess);
  }
}