G4EmExtraPhysics Class Reference

#include <G4EmExtraPhysics.hh>

Inheritance diagram for G4EmExtraPhysics:

Public Member Functions
	G4EmExtraPhysics (G4int ver=1)
	G4EmExtraPhysics (const G4String &name)
virtual	~G4EmExtraPhysics ()
void	ConstructParticle ()
void	ConstructProcess ()
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	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
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4int	GetInstanceID () const
virtual void	TerminateWorker ()

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
G4String	namePhysics
G4int	typePhysics
G4ParticleTable *	theParticleTable
G4int	g4vpcInstanceID
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 123 of file G4EmExtraPhysics.cc.

                                           : 
  G4VPhysicsConstructor("G4GammaLeptoNuclearPhys"),
  gnActivated (true),
  eActivated  (true),
  gLENDActivated(false),
  munActivated(true),
  synActivated(false),
  synActivatedForAll(false),
  gmumuActivated(false),
  pmumuActivated(false),
  phadActivated (false),
  fNuActivated (false),
  fNuETotXscActivated (false),
  fUseGammaNuclearXS(false),
  gmumuFactor (1.0),
  pmumuFactor (1.0),
  phadFactor  (1.0),
  fNuEleCcBias(1.0),
  fNuEleNcBias(1.0),
  fNuNucleusBias(1.0),
  fGNLowEnergyLimit(200*CLHEP::MeV),
  fNuDetectorName("0"),
  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 152 of file G4EmExtraPhysics.cc.

  : G4EmExtraPhysics(1)
{}

~G4EmExtraPhysics()

G4EmExtraPhysics::~G4EmExtraPhysics ( )

virtual

Definition at line 156 of file G4EmExtraPhysics.cc.

{
  delete theMessenger;
  theMessenger = nullptr;
}

Member Function Documentation

ConstructParticle()

void G4EmExtraPhysics::ConstructParticle ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 276 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 ( )

virtual

Implements G4VPhysicsConstructor.

Definition at line 292 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 = 
      (G4GammaGeneralProcess*)emManager->GetGammaGeneralProcess();
    if(sp) {
      sp->AddMMProcess(theGammaToMuMu);
    } else {
      ph->RegisterProcess(theGammaToMuMu, gamma);
    }
  }  
  if(pmumuActivated) {
    G4AnnihiToMuPair* thePosiToMuMu = new G4AnnihiToMuPair();
    thePosiToMuMu->SetCrossSecFactor(pmumuFactor);
    ph->RegisterProcess(thePosiToMuMu, 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) {
      auto myParticleIterator=GetParticleIterator();
      myParticleIterator->reset();
      G4ParticleDefinition* particle = nullptr;
 
      while( (*myParticleIterator)() ) {
    particle = myParticleIterator->value();
    if( particle->GetPDGStable() && particle->GetPDGCharge() != 0.0) { 
      if(verbose > 1) {
        G4cout << "### G4SynchrotronRadiation for " 
           << particle->GetParticleName() << G4endl;
      }
      ph->RegisterProcess( theSynchRad, particle);
    }
      }
    }
  }
  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_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 185 of file G4EmExtraPhysics.cc.

{
  eActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

GammaNuclear()

void G4EmExtraPhysics::GammaNuclear

(

G4bool

val

)

Definition at line 173 of file G4EmExtraPhysics.cc.

{
  gnActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

GammaNuclearLEModelLimit()

void G4EmExtraPhysics::GammaNuclearLEModelLimit

(

G4double

val

)

Definition at line 257 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 195 of file G4EmExtraPhysics.cc.

{
  gmumuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

GammaToMuMuFactor()

void G4EmExtraPhysics::GammaToMuMuFactor

(

G4double

val

)

Definition at line 210 of file G4EmExtraPhysics.cc.

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

Referenced by G4EmMessenger::SetNewValue().

LENDGammaNuclear()

void G4EmExtraPhysics::LENDGammaNuclear

(

G4bool

val

)

Definition at line 178 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 190 of file G4EmExtraPhysics.cc.

{
  munActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

NeutrinoActivated()

void G4EmExtraPhysics::NeutrinoActivated

(

G4bool

val

)

Definition at line 227 of file G4EmExtraPhysics.cc.

{
  fNuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

NuETotXscActivated()

void G4EmExtraPhysics::NuETotXscActivated

(

G4bool

val

)

Definition at line 232 of file G4EmExtraPhysics.cc.

{
  fNuETotXscActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToHadrons()

void G4EmExtraPhysics::PositronToHadrons

(

G4bool

val

)

Definition at line 205 of file G4EmExtraPhysics.cc.

{
  phadActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToHadronsFactor()

void G4EmExtraPhysics::PositronToHadronsFactor

(

G4double

val

)

Definition at line 220 of file G4EmExtraPhysics.cc.

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

Referenced by G4EmMessenger::SetNewValue().

PositronToMuMu()

void G4EmExtraPhysics::PositronToMuMu

(

G4bool

val

)

Definition at line 200 of file G4EmExtraPhysics.cc.

{
  pmumuActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

PositronToMuMuFactor()

void G4EmExtraPhysics::PositronToMuMuFactor

(

G4double

val

)

Definition at line 215 of file G4EmExtraPhysics.cc.

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

Referenced by G4EmMessenger::SetNewValue().

SetNuDetectorName()

void G4EmExtraPhysics::SetNuDetectorName

(

const G4String &

dn

)

Definition at line 269 of file G4EmExtraPhysics.cc.

{
  fNuDetectorName = dn;
}

Referenced by G4EmMessenger::SetNewValue().

SetNuEleCcBias()

void G4EmExtraPhysics::SetNuEleCcBias

(

G4double

bf

)

Definition at line 242 of file G4EmExtraPhysics.cc.

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

Referenced by G4EmMessenger::SetNewValue().

SetNuEleNcBias()

void G4EmExtraPhysics::SetNuEleNcBias

(

G4double

bf

)

Definition at line 247 of file G4EmExtraPhysics.cc.

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

Referenced by G4EmMessenger::SetNewValue().

SetNuNucleusBias()

void G4EmExtraPhysics::SetNuNucleusBias

(

G4double

bf

)

Definition at line 252 of file G4EmExtraPhysics.cc.

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

Referenced by G4EmMessenger::SetNewValue().

SetUseGammaNuclearXS()

void G4EmExtraPhysics::SetUseGammaNuclearXS

(

G4bool

val

)

Definition at line 237 of file G4EmExtraPhysics.cc.

{
  fUseGammaNuclearXS = val;
}

Referenced by G4EmMessenger::SetNewValue().

Synch()

void G4EmExtraPhysics::Synch

(

G4bool

val

)

Definition at line 162 of file G4EmExtraPhysics.cc.

{
  synActivated = val;
}

Referenced by G4EmMessenger::SetNewValue().

SynchAll()

void G4EmExtraPhysics::SynchAll

(

G4bool

val

)

Definition at line 167 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