#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;
}

Referenced by G4NeutrinoPhysics(), and operator=().

◆ ~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;
}

◆ 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;
}

◆ SetNuEleCcBias()

void G4NeutrinoPhysics::SetNuEleCcBias ( G4double bf )

Definition at line 111 of file G4NeutrinoPhysics.cc.

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

◆ SetNuEleNcBias()

void G4NeutrinoPhysics::SetNuEleNcBias ( G4double bf )

Definition at line 116 of file G4NeutrinoPhysics.cc.

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

◆ SetNuNucleusBias()

void G4NeutrinoPhysics::SetNuNucleusBias ( G4double bf )

Definition at line 121 of file G4NeutrinoPhysics.cc.

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

◆ SetNuOscDistanceBias()

void G4NeutrinoPhysics::SetNuOscDistanceBias ( G4double bf )

Definition at line 126 of file G4NeutrinoPhysics.cc.

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

◆ SetNuOscDistanceName()

void G4NeutrinoPhysics::SetNuOscDistanceName ( const G4String & dn )

Definition at line 136 of file G4NeutrinoPhysics.cc.

{
  fNuOscDistanceName = dn;
}

◆ 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:

Public Member Functions

Additional Inherited Members