G4EmExtraPhysics Class Reference

#include <G4EmExtraPhysics.hh>

Inheritance diagram for G4EmExtraPhysics:

Public Member Functions
	G4EmExtraPhysics (G4int ver=1)
	G4EmExtraPhysics (const G4String &name)
	~G4EmExtraPhysics () override
void	ConstructParticle () override
void	ConstructProcess () override
void	Synch (G4bool val)
void	SynchAll (G4bool val)
void	GammaNuclear (G4bool val)
void	LENDGammaNuclear (G4bool val)
void	ElectroNuclear (G4bool val)
void	MuonNuclear (G4bool val)
void	GammaToMuMu (G4bool val)
void	MuonToMuMu (G4bool val)
void	PositronToMuMu (G4bool val)
void	PositronToHadrons (G4bool val)
void	GammaToMuMuFactor (G4double val)
void	PositronToMuMuFactor (G4double val)
void	PositronToHadronsFactor (G4double val)
void	GammaNuclearLEModelLimit (G4double val)
void	NeutrinoActivated (G4bool val)
void	NuETotXscActivated (G4bool val)
void	SetUseGammaNuclearXS (G4bool val)
void	SetNuEleCcBias (G4double bf)
void	SetNuEleNcBias (G4double bf)
void	SetNuNucleusBias (G4double bf)
void	SetNuDetectorName (const G4String &dn)
Public Member Functions inherited from G4VPhysicsConstructor
	G4VPhysicsConstructor (const G4String &="")
	G4VPhysicsConstructor (const G4String &name, G4int physics_type)
virtual	~G4VPhysicsConstructor ()
virtual void	ConstructParticle ()=0
virtual void	ConstructProcess ()=0
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 53 of file G4EmExtraPhysics.hh.

Constructor & Destructor Documentation

G4EmExtraPhysics() [1/2]

G4EmExtraPhysics::G4EmExtraPhysics

(

G4int

ver = 1

)

Definition at line 140 of file G4EmExtraPhysics.cc.

                                           : 
  G4VPhysicsConstructor("G4GammaLeptoNuclearPhys"),
  fGNLowEnergyLimit(200*CLHEP::MeV),
  verbose(ver)
{
  theMessenger = new G4EmMessenger(this);
  SetPhysicsType(bEmExtra);
  if(verbose > 1) G4cout << "### G4EmExtraPhysics" << G4endl;
}

G4EmExtraPhysics() [2/2]

G4EmExtraPhysics::G4EmExtraPhysics

(

const G4String &

name

)

Definition at line 150 of file G4EmExtraPhysics.cc.

  : G4EmExtraPhysics(1)
{}

~G4EmExtraPhysics()

G4EmExtraPhysics::~G4EmExtraPhysics ( )

override

Definition at line 154 of file G4EmExtraPhysics.cc.

{
  delete theMessenger;
  theMessenger = nullptr;
}

Member Function Documentation

ConstructParticle()

void G4EmExtraPhysics::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 279 of file G4EmExtraPhysics.cc.

{
  G4Gamma::Gamma();
  G4Electron::Electron();
  G4Positron::Positron();
  G4MuonPlus::MuonPlus();
  G4MuonMinus::MuonMinus();
 
  G4AntiNeutrinoE::AntiNeutrinoE();
  G4NeutrinoE::NeutrinoE();
  G4AntiNeutrinoMu::AntiNeutrinoMu();
  G4NeutrinoMu::NeutrinoMu();
  G4AntiNeutrinoTau::AntiNeutrinoTau();
  G4NeutrinoTau::NeutrinoTau();
}

ConstructProcess()

void G4EmExtraPhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 295 of file G4EmExtraPhysics.cc.

{
  G4ParticleDefinition* gamma     = G4Gamma::Gamma();
  G4ParticleDefinition* electron  = G4Electron::Electron();
  G4ParticleDefinition* positron  = G4Positron::Positron();
  G4ParticleDefinition* muonplus  = G4MuonPlus::MuonPlus();
  G4ParticleDefinition* muonminus = G4MuonMinus::MuonMinus();
 
  G4PhysicsListHelper* ph = G4PhysicsListHelper::GetPhysicsListHelper();
  G4LossTableManager* emManager = G4LossTableManager::Instance();
 
  if(gnActivated) { ConstructGammaElectroNuclear(); }
 
  if(munActivated) {
    G4MuonNuclearProcess* muNucProcess = new G4MuonNuclearProcess();
    G4MuonVDNuclearModel* muNucModel = new G4MuonVDNuclearModel();
    muNucProcess->RegisterMe(muNucModel);
    ph->RegisterProcess( muNucProcess, muonplus);
    ph->RegisterProcess( muNucProcess, muonminus);
  }
  if(gmumuActivated) {
    G4GammaConversionToMuons* theGammaToMuMu = new G4GammaConversionToMuons();
    theGammaToMuMu->SetCrossSecFactor(gmumuFactor);
    G4GammaGeneralProcess* sp = 
      static_cast<G4GammaGeneralProcess*>(emManager->GetGammaGeneralProcess());
    if(nullptr != sp) {
      sp->AddMMProcess(theGammaToMuMu);
    } else {
      ph->RegisterProcess(theGammaToMuMu, gamma);
    }
  }  
  if(mmumuActivated) {
    auto proc = new G4MuonToMuonPairProduction();
    ph->RegisterProcess(proc, muonplus);
    ph->RegisterProcess(proc, muonminus);
  }
  if(pmumuActivated) {
    G4AnnihiToMuPair* thePosiToMuMu = new G4AnnihiToMuPair();
    thePosiToMuMu->SetCrossSecFactor(pmumuFactor);
    ph->RegisterProcess(thePosiToMuMu, positron);
    G4AnnihiToMuPair* thePosiToTauTau = new G4AnnihiToMuPair("AnnihiToTauPair");
    thePosiToTauTau->SetCrossSecFactor(pmumuFactor);
    ph->RegisterProcess(thePosiToTauTau, positron);
  }  
  if(phadActivated) {
    G4eeToHadrons* thePosiToHadrons = new G4eeToHadrons();
    thePosiToHadrons->SetCrossSecFactor(phadFactor);
    ph->RegisterProcess(thePosiToHadrons, positron);
  }  
  if(synActivated) {
    G4SynchrotronRadiation* theSynchRad = new G4SynchrotronRadiation();
    ph->RegisterProcess( theSynchRad, electron);
    ph->RegisterProcess( theSynchRad, positron);
    if(synActivatedForAll) {
      ph->RegisterProcess( theSynchRad, muonplus);
      ph->RegisterProcess( theSynchRad, muonminus);
 
      ph->RegisterProcess( theSynchRad, G4Proton::Proton());
      ph->RegisterProcess( theSynchRad, G4AntiProton::AntiProton());
      ph->RegisterProcess( theSynchRad, G4PionPlus::PionPlus());
      ph->RegisterProcess( theSynchRad, G4PionMinus::PionMinus());
      ph->RegisterProcess( theSynchRad, G4GenericIon::GenericIon());
    }
  }
  if( fNuActivated )
  {
    G4ParticleDefinition* anuelectron = G4AntiNeutrinoE::AntiNeutrinoE();
    G4ParticleDefinition* nuelectron  = G4NeutrinoE::NeutrinoE(); 
    G4ParticleDefinition* anumuon     = G4AntiNeutrinoMu::AntiNeutrinoMu();
    G4ParticleDefinition* numuon      = G4NeutrinoMu::NeutrinoMu();
    G4ParticleDefinition* anutau      = G4AntiNeutrinoTau::AntiNeutrinoTau();
    G4ParticleDefinition* nutau       = G4NeutrinoTau::NeutrinoTau();
 
    G4NeutrinoElectronProcess* theNuEleProcess = 
      new G4NeutrinoElectronProcess(fNuDetectorName);
    G4NeutrinoElectronTotXsc* theNuEleTotXsc = new G4NeutrinoElectronTotXsc();
 
    if(fNuETotXscActivated)
    {
      G4double bftot = std::max(fNuEleCcBias,fNuEleNcBias);
      theNuEleProcess->SetBiasingFactor(bftot);
    }
    else
    {
      theNuEleProcess->SetBiasingFactors(fNuEleCcBias,fNuEleNcBias);
      theNuEleTotXsc->SetBiasingFactors(fNuEleCcBias,fNuEleNcBias);
    }
    theNuEleProcess->AddDataSet(theNuEleTotXsc);
 
    G4NeutrinoElectronCcModel* ccModel = new G4NeutrinoElectronCcModel();
    G4NeutrinoElectronNcModel* ncModel = new G4NeutrinoElectronNcModel();
    theNuEleProcess->RegisterMe(ccModel);
    theNuEleProcess->RegisterMe(ncModel);
 
    ph->RegisterProcess(theNuEleProcess, anuelectron);
    ph->RegisterProcess(theNuEleProcess, nuelectron);
    ph->RegisterProcess(theNuEleProcess, anumuon);
    ph->RegisterProcess(theNuEleProcess, numuon);
    ph->RegisterProcess(theNuEleProcess, anutau);
    ph->RegisterProcess(theNuEleProcess, nutau);
 
    // nu_mu nucleus interactions
 
    G4MuNeutrinoNucleusProcess* theNuMuNucleusProcess = new G4MuNeutrinoNucleusProcess(fNuDetectorName);
    G4MuNeutrinoNucleusTotXsc* theNuMuNucleusTotXsc = new G4MuNeutrinoNucleusTotXsc();
    
    if(fNuETotXscActivated)
    {
      theNuMuNucleusProcess->SetBiasingFactor(fNuNucleusBias);
    }
    theNuMuNucleusProcess->AddDataSet(theNuMuNucleusTotXsc);
 
    G4NuMuNucleusCcModel* numunuclcc = new G4NuMuNucleusCcModel();
    G4NuMuNucleusNcModel* numunuclnc = new G4NuMuNucleusNcModel();
    G4ANuMuNucleusCcModel* anumunuclcc = new G4ANuMuNucleusCcModel();
    G4ANuMuNucleusNcModel* anumunuclnc = new G4ANuMuNucleusNcModel();
    
    theNuMuNucleusProcess->RegisterMe(numunuclcc);
    theNuMuNucleusProcess->RegisterMe(numunuclnc);
    theNuMuNucleusProcess->RegisterMe(anumunuclcc);
    theNuMuNucleusProcess->RegisterMe(anumunuclnc);
 
    ph->RegisterProcess(theNuMuNucleusProcess, anumuon);
    ph->RegisterProcess(theNuMuNucleusProcess, numuon);    
 
    // nu_tau nucleus interactions
 
    G4TauNeutrinoNucleusProcess* theNuTauNucleusProcess = new G4TauNeutrinoNucleusProcess(fNuDetectorName);
    G4TauNeutrinoNucleusTotXsc* theNuTauNucleusTotXsc = new G4TauNeutrinoNucleusTotXsc();
    
    if(fNuETotXscActivated)
    {
      theNuTauNucleusProcess->SetBiasingFactor(fNuNucleusBias);
    }
    theNuTauNucleusProcess->AddDataSet(theNuTauNucleusTotXsc);
 
    G4NuTauNucleusCcModel* nutaunuclcc = new G4NuTauNucleusCcModel();
    G4NuTauNucleusNcModel* nutaunuclnc = new G4NuTauNucleusNcModel();
    G4ANuTauNucleusCcModel* anutaunuclcc = new G4ANuTauNucleusCcModel();
    G4ANuTauNucleusNcModel* anutaunuclnc = new G4ANuTauNucleusNcModel();
    
    theNuTauNucleusProcess->RegisterMe(nutaunuclcc);
    theNuTauNucleusProcess->RegisterMe(nutaunuclnc);
    theNuTauNucleusProcess->RegisterMe(anutaunuclcc);
    theNuTauNucleusProcess->RegisterMe(anutaunuclnc);
 
    ph->RegisterProcess(theNuTauNucleusProcess, anutau);
    ph->RegisterProcess(theNuTauNucleusProcess, nutau);    
 
    // nu_e nucleus interactions
 
    G4ElNeutrinoNucleusProcess* theNuElNucleusProcess = new G4ElNeutrinoNucleusProcess(fNuDetectorName);
    G4ElNeutrinoNucleusTotXsc* theNuElNucleusTotXsc = new G4ElNeutrinoNucleusTotXsc();
    
    if(fNuETotXscActivated)
    {
      theNuElNucleusProcess->SetBiasingFactor(fNuNucleusBias);
    }
    theNuElNucleusProcess->AddDataSet(theNuElNucleusTotXsc);
 
    G4NuElNucleusCcModel* nuelnuclcc = new G4NuElNucleusCcModel();
    G4NuElNucleusNcModel* nuelnuclnc = new G4NuElNucleusNcModel();
    G4ANuElNucleusCcModel* anuelnuclcc = new G4ANuElNucleusCcModel();
    G4ANuElNucleusNcModel* anuelnuclnc = new G4ANuElNucleusNcModel();
    
    theNuElNucleusProcess->RegisterMe(nuelnuclcc);
    theNuElNucleusProcess->RegisterMe(nuelnuclnc);
    theNuElNucleusProcess->RegisterMe(anuelnuclcc);
    theNuElNucleusProcess->RegisterMe(anuelnuclnc);
 
    ph->RegisterProcess(theNuElNucleusProcess, anuelectron);
    ph->RegisterProcess(theNuElNucleusProcess, nuelectron);    
  }
}

ElectroNuclear()

void G4EmExtraPhysics::ElectroNuclear

(

G4bool

val

)

Definition at line 183 of file G4EmExtraPhysics.cc.

{
  eActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

GammaNuclear()

void G4EmExtraPhysics::GammaNuclear

(

G4bool

val

)

Definition at line 171 of file G4EmExtraPhysics.cc.

{
  gnActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

GammaNuclearLEModelLimit()

void G4EmExtraPhysics::GammaNuclearLEModelLimit

(

G4double

val

)

Definition at line 260 of file G4EmExtraPhysics.cc.

{
  if(val <= CLHEP::MeV) { 
    fGNLowEnergyLimit = 0.0;
 
    // lowenergy model should not be applied at high energy
  } else if(val <= CLHEP::GeV) {
    fGNLowEnergyLimit = val;
    gLENDActivated = false;
  }
}

Referenced by G4EmMessenger::SetNewValue().

GammaToMuMu()

void G4EmExtraPhysics::GammaToMuMu

(

G4bool

val

)

Definition at line 193 of file G4EmExtraPhysics.cc.

{
  gmumuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

GammaToMuMuFactor()

void G4EmExtraPhysics::GammaToMuMuFactor

(

G4double

val

)

Definition at line 213 of file G4EmExtraPhysics.cc.

{
  if(val > 0.0) gmumuFactor = val;
}

Referenced by G4EmMessenger::SetNewValue().

LENDGammaNuclear()

void G4EmExtraPhysics::LENDGammaNuclear

(

G4bool

val

)

Definition at line 176 of file G4EmExtraPhysics.cc.

{
  gLENDActivated = val;
  // LEND cannot be used with low-energy model
  if(val) { fGNLowEnergyLimit = 0.0; }
}

Referenced by G4EmMessenger::SetNewValue(), and Shielding::Shielding().

MuonNuclear()

void G4EmExtraPhysics::MuonNuclear

(

G4bool

val

)

Definition at line 188 of file G4EmExtraPhysics.cc.

{
  munActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

MuonToMuMu()

void G4EmExtraPhysics::MuonToMuMu

(

G4bool

val

)

Definition at line 198 of file G4EmExtraPhysics.cc.

{
  mmumuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

NeutrinoActivated()

void G4EmExtraPhysics::NeutrinoActivated

(

G4bool

val

)

Definition at line 230 of file G4EmExtraPhysics.cc.

{
  fNuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

NuETotXscActivated()

void G4EmExtraPhysics::NuETotXscActivated

(

G4bool

val

)

Definition at line 235 of file G4EmExtraPhysics.cc.

{
  fNuETotXscActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToHadrons()

void G4EmExtraPhysics::PositronToHadrons

(

G4bool

val

)

Definition at line 208 of file G4EmExtraPhysics.cc.

{
  phadActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToHadronsFactor()

void G4EmExtraPhysics::PositronToHadronsFactor

(

G4double

val

)

Definition at line 223 of file G4EmExtraPhysics.cc.

{
  if(val > 0.0) phadFactor = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToMuMu()

void G4EmExtraPhysics::PositronToMuMu

(

G4bool

val

)

Definition at line 203 of file G4EmExtraPhysics.cc.

{
  pmumuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToMuMuFactor()

void G4EmExtraPhysics::PositronToMuMuFactor

(

G4double

val

)

Definition at line 218 of file G4EmExtraPhysics.cc.

{
  if(val > 0.0) pmumuFactor = val;
}

Referenced by G4EmMessenger::SetNewValue().

SetNuDetectorName()

void G4EmExtraPhysics::SetNuDetectorName

(

const G4String &

dn

)

Definition at line 272 of file G4EmExtraPhysics.cc.

{
  fNuDetectorName = dn;
}

Referenced by G4EmMessenger::SetNewValue().

SetNuEleCcBias()

void G4EmExtraPhysics::SetNuEleCcBias

(

G4double

bf

)

Definition at line 245 of file G4EmExtraPhysics.cc.

{
  if(bf > 0.0) fNuEleCcBias = bf;
}

Referenced by G4EmMessenger::SetNewValue().

SetNuEleNcBias()

void G4EmExtraPhysics::SetNuEleNcBias

(

G4double

bf

)

Definition at line 250 of file G4EmExtraPhysics.cc.

{
  if(bf > 0.0) fNuEleNcBias = bf;
}

Referenced by G4EmMessenger::SetNewValue().

SetNuNucleusBias()

void G4EmExtraPhysics::SetNuNucleusBias

(

G4double

bf

)

Definition at line 255 of file G4EmExtraPhysics.cc.

{
  if(bf > 0.0) fNuNucleusBias = bf;
}

Referenced by G4EmMessenger::SetNewValue().

SetUseGammaNuclearXS()

void G4EmExtraPhysics::SetUseGammaNuclearXS

(

G4bool

val

)

Definition at line 240 of file G4EmExtraPhysics.cc.

{
  fUseGammaNuclearXS = val;
}

Referenced by G4EmMessenger::SetNewValue().

Synch()

void G4EmExtraPhysics::Synch

(

G4bool

val

)

Definition at line 160 of file G4EmExtraPhysics.cc.

{
  synActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

SynchAll()

void G4EmExtraPhysics::SynchAll

(

G4bool

val

)

Definition at line 165 of file G4EmExtraPhysics.cc.

{
  synActivatedForAll = val;
  if(synActivatedForAll) { synActivated = true; }
}

Referenced by G4EmMessenger::SetNewValue().

---

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

• G4EmExtraPhysics.hh
• G4EmExtraPhysics.cc