G4NeutrinoPhysics Class Reference

#include <G4NeutrinoPhysics.hh>

Inheritance diagram for G4NeutrinoPhysics:

Public Member Functions
	G4NeutrinoPhysics (G4int ver=1)
	~G4NeutrinoPhysics () override
void	ConstructParticle () override
void	ConstructProcess () override
void	NuETotXscActivated (G4bool val)
void	SetNuOscillation (G4bool val)
void	SetNuEleCcBias (G4double bf)
void	SetNuEleNcBias (G4double bf)
void	SetNuNucleusBias (G4double bf)
void	SetNuOscDistanceBias (G4double bf)
void	SetNuDetectorName (const G4String &dn)
void	SetNuOscDistanceName (const G4String &dn)
G4NeutrinoPhysics &	operator= (const G4NeutrinoPhysics &right)=delete
	G4NeutrinoPhysics (const G4NeutrinoPhysics &)=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 45 of file G4NeutrinoPhysics.hh.

Constructor & Destructor Documentation

G4NeutrinoPhysics() [1/2]

G4NeutrinoPhysics::G4NeutrinoPhysics

(

G4int

ver = 1

)

Definition at line 88 of file G4NeutrinoPhysics.cc.

                                             : 
  G4VPhysicsConstructor("NeutrinoPhys"),
  verbose(ver)
{
  theMessenger = new G4NeutrinoPhysicsMessenger(this);
  if(verbose > 1) G4cout << "### G4NeutrinoPhysics" << G4endl;
}

~G4NeutrinoPhysics()

G4NeutrinoPhysics::~G4NeutrinoPhysics ( )

override

Definition at line 96 of file G4NeutrinoPhysics.cc.

{
  delete theMessenger;
}

G4NeutrinoPhysics() [2/2]

G4NeutrinoPhysics::G4NeutrinoPhysics ( const G4NeutrinoPhysics & )

delete

Member Function Documentation

ConstructParticle()

void G4NeutrinoPhysics::ConstructParticle ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 143 of file G4NeutrinoPhysics.cc.

{
  G4Electron::Electron();
  G4AntiNeutrinoE::AntiNeutrinoE();
  G4NeutrinoE::NeutrinoE();
  G4AntiNeutrinoMu::AntiNeutrinoMu();
  G4NeutrinoMu::NeutrinoMu();
  G4AntiNeutrinoTau::AntiNeutrinoTau();
  G4NeutrinoTau::NeutrinoTau();
}

ConstructProcess()

void G4NeutrinoPhysics::ConstructProcess ( )

overridevirtual

Implements G4VPhysicsConstructor.

Definition at line 154 of file G4NeutrinoPhysics.cc.

{
  const G4ParticleDefinition* p[6] = {
    G4AntiNeutrinoE::AntiNeutrinoE(),
    G4NeutrinoE::NeutrinoE(),
    G4AntiNeutrinoMu::AntiNeutrinoMu(),
    G4NeutrinoMu::NeutrinoMu(),
    G4AntiNeutrinoTau::AntiNeutrinoTau(),
    G4NeutrinoTau::NeutrinoTau()
  };
 
  // neutrino vacuum oscillation process
  if (fNuOscillation) {
    auto theNuVacOscProcess = new G4NuVacOscProcess(fNuOscDistanceName);
    theNuVacOscProcess->SetBiasingFactor(fNuOscDistanceBias);
    
 
    for (G4int i=0; i<6; ++i) {
      p[i]->GetProcessManager()->AddDiscreteProcess(theNuVacOscProcess);
    }
  }
 
  // neutrino-electron process
  auto 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);
 
  for (G4int i=0; i<6; ++i) {
    p[i]->GetProcessManager()->AddDiscreteProcess(theNuEleProcess);
  }
 
  // nu_mu nucleus interactions
  auto theNuMuNucleusProcess = new G4MuNeutrinoNucleusProcess(fNuDetectorName);
  auto 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);
 
  for (G4int i=2; i<=3; ++i) {
    p[i]->GetProcessManager()->AddDiscreteProcess(theNuMuNucleusProcess);
  }
 
  // nu_tau nucleus interactions
  auto theNuTauNucleusProcess = new G4TauNeutrinoNucleusProcess(fNuDetectorName);
  auto 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);
 
  for (G4int i=4; i<=5; ++i) {
    p[i]->GetProcessManager()->AddDiscreteProcess(theNuTauNucleusProcess);
  }
 
  // nu_e nucleus interactions
  auto theNuElNucleusProcess = new G4ElNeutrinoNucleusProcess(fNuDetectorName);
  auto 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);
 
  for (G4int i=0; i<=1; ++i) {
    p[i]->GetProcessManager()->AddDiscreteProcess(theNuElNucleusProcess);
  }
}

NuETotXscActivated()

void G4NeutrinoPhysics::NuETotXscActivated

(

G4bool

val

)

Definition at line 101 of file G4NeutrinoPhysics.cc.

{
  fNuETotXscActivated = val;
}

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

operator=()

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

delete

SetNuDetectorName()

void G4NeutrinoPhysics::SetNuDetectorName

(

const G4String &

dn

)

Definition at line 131 of file G4NeutrinoPhysics.cc.

{
  fNuDetectorName = dn;
}

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

SetNuEleCcBias()

void G4NeutrinoPhysics::SetNuEleCcBias

(

G4double

bf

)

Definition at line 111 of file G4NeutrinoPhysics.cc.

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

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

SetNuEleNcBias()

void G4NeutrinoPhysics::SetNuEleNcBias

(

G4double

bf

)

Definition at line 116 of file G4NeutrinoPhysics.cc.

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

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

SetNuNucleusBias()

void G4NeutrinoPhysics::SetNuNucleusBias

(

G4double

bf

)

Definition at line 121 of file G4NeutrinoPhysics.cc.

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

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

SetNuOscDistanceBias()

void G4NeutrinoPhysics::SetNuOscDistanceBias

(

G4double

bf

)

Definition at line 126 of file G4NeutrinoPhysics.cc.

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

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

SetNuOscDistanceName()

void G4NeutrinoPhysics::SetNuOscDistanceName

(

const G4String &

dn

)

Definition at line 136 of file G4NeutrinoPhysics.cc.

{
  fNuOscDistanceName = dn;
}

Referenced by G4NeutrinoPhysicsMessenger::SetNewValue().

SetNuOscillation()

void G4NeutrinoPhysics::SetNuOscillation

(

G4bool

val

)

Definition at line 106 of file G4NeutrinoPhysics.cc.

{
  fNuOscillation = val;
}

---

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

• G4NeutrinoPhysics.hh
• G4NeutrinoPhysics.cc