#include <G4NeutrinoElectronNcXsc.hh>

Inheritance diagram for G4NeutrinoElectronNcXsc:

Public Member Functions
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

virtual void	CrossSectionDescription (std::ostream &) const

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

bool	ForAllAtomsAndEnergies () const

	G4NeutrinoElectronNcXsc ()

G4double	GetBiasingFactor ()

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=nullptr)

G4double	GetCutEnergy ()

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material )

virtual G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=nullptr, const G4Element elm=nullptr, const G4Material mat=nullptr)

G4double	GetMaxKinEnergy () const

G4double	GetMinKinEnergy () const

const G4String &	GetName () const

virtual G4int	GetVerboseLevel () const

virtual G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material )

virtual G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element elm=nullptr, const G4Material *mat=nullptr)

virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)

void	SetBiasingFactor (G4double bf)

void	SetCutEnergy (G4double ec)

void	SetForAllAtomsAndEnergies (G4bool val)

void	SetMaxKinEnergy (G4double value)

void	SetMinKinEnergy (G4double value)

void	SetName (const G4String &nam)

virtual void	SetVerboseLevel (G4int value)

	~G4NeutrinoElectronNcXsc ()

Protected Attributes
G4double	fBiasingFactor

G4double	fCofXsc

G4double	fCutEnergy

G4double	fSin2tW

G4String	name

G4int	verboseLevel

Private Attributes
G4bool	isForAllAtomsAndEnergies

G4double	maxKinEnergy

G4double	minKinEnergy

G4CrossSectionDataSetRegistry *	registry

Detailed Description

Definition at line 39 of file G4NeutrinoElectronNcXsc.hh.

Constructor & Destructor Documentation

◆ G4NeutrinoElectronNcXsc()

G4NeutrinoElectronNcXsc::G4NeutrinoElectronNcXsc ( )

Definition at line 41 of file G4NeutrinoElectronNcXsc.cc.

 : G4VCrossSectionDataSet("NuElectronNcXsc")
{
  // PDG2016: Gf=1.1663787(6)e-5*(hc)^3/GeV^2
  // fCofXsc  = Gf*Gf*MeC2*2/pi
 
  fCofXsc  = 1.36044e-22;
  fCofXsc *= hbarc*hbarc*electron_mass_c2;
  fCofXsc /= halfpi;
 
  // G4cout<<"hbarc = "<<hbarc/MeV/fermi<<" MeV*fermi"<<G4endl;
 
  // PDG2016: sin^2 theta Weinberg
 
  fSin2tW = 0.23129; // 0.2312;
 
  fCutEnergy = 0.; // default value
  fBiasingFactor = 1.; 
}

References source.hepunit::electron_mass_c2, fBiasingFactor, fCofXsc, fCutEnergy, fSin2tW, halfpi, and source.hepunit::hbarc.

◆ ~G4NeutrinoElectronNcXsc()

G4NeutrinoElectronNcXsc::~G4NeutrinoElectronNcXsc ( )

Definition at line 61 of file G4NeutrinoElectronNcXsc.cc.

62{}

Member Function Documentation

◆ BuildPhysicsTable()

void G4VCrossSectionDataSet::BuildPhysicsTable ( const G4ParticleDefinition & )

virtualinherited

◆ ComputeCrossSection()

G4double G4VCrossSectionDataSet::ComputeCrossSection	(	const G4DynamicParticle *	part,
		const G4Element *	elm,
		const G4Material *	mat = `nullptr`
	)

inherited

Definition at line 81 of file G4VCrossSectionDataSet.cc.

{
  G4int Z = elm->GetZasInt();
 
  if (IsElementApplicable(part, Z, mat)) { 
    return GetElementCrossSection(part, Z, mat);
  }
 
  // isotope-wise cross section making sum over available
  // isotope cross sections, which may be incomplete, so
  // the result is corrected 
  size_t nIso = elm->GetNumberOfIsotopes();    
  G4double fact = 0.0;
  G4double xsec = 0.0;
 
  // user-defined isotope abundances        
  const G4IsotopeVector* isoVector = elm->GetIsotopeVector();
  const G4double* abundVector = elm->GetRelativeAbundanceVector();
 
  for (size_t j=0; j<nIso; ++j) {
    const G4Isotope* iso = (*isoVector)[j];
    G4int A = iso->GetN();
    if(abundVector[j] > 0.0 && IsIsoApplicable(part, Z, A, elm, mat)) {
      fact += abundVector[j];
      xsec += abundVector[j]*GetIsoCrossSection(part, Z, A, iso, elm, mat);
    }
  }
  return (fact > 0.0) ? xsec/fact : 0.0; 
}

References A, G4VCrossSectionDataSet::GetElementCrossSection(), G4VCrossSectionDataSet::GetIsoCrossSection(), G4Element::GetIsotopeVector(), G4Isotope::GetN(), G4Element::GetNumberOfIsotopes(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZasInt(), G4VCrossSectionDataSet::IsElementApplicable(), G4VCrossSectionDataSet::IsIsoApplicable(), and Z.

Referenced by G4VCrossSectionDataSet::GetCrossSection().

◆ CrossSectionDescription()

void G4VCrossSectionDataSet::CrossSectionDescription ( std::ostream & outFile ) const

virtualinherited

Definition at line 177 of file G4VCrossSectionDataSet.cc.

{
  outFile << "The description for this cross section data set has not been written yet.\n";
}

Referenced by G4CrossSectionDataStore::PrintCrossSectionHtml().

◆ DumpPhysicsTable()

void G4VCrossSectionDataSet::DumpPhysicsTable ( const G4ParticleDefinition & )

virtualinherited

Reimplemented in G4CrossSectionElastic, G4CrossSectionInelastic, G4LENDCrossSection, G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, G4ParticleHPInelasticData, G4ParticleHPJENDLHEData, G4ParticleHPThermalScatteringData, and G4UPiNuclearCrossSection.

Definition at line 174 of file G4VCrossSectionDataSet.cc.

175{}

◆ ForAllAtomsAndEnergies()

bool G4VCrossSectionDataSet::ForAllAtomsAndEnergies ( ) const

inlineinherited

Definition at line 230 of file G4VCrossSectionDataSet.hh.

{
  return isForAllAtomsAndEnergies;
}

References G4VCrossSectionDataSet::isForAllAtomsAndEnergies.

Referenced by G4CrossSectionDataStore::AddDataSet().

◆ GetBiasingFactor()

G4double G4NeutrinoElectronNcXsc::GetBiasingFactor ( )

inline

Definition at line 58 of file G4NeutrinoElectronNcXsc.hh.

58{return fBiasingFactor;};

References fBiasingFactor.

◆ GetCrossSection()

G4double G4VCrossSectionDataSet::GetCrossSection	(	const G4DynamicParticle *	dp,
		const G4Element *	elm,
		const G4Material *	mat = `nullptr`
	)

inlineinherited

Definition at line 187 of file G4VCrossSectionDataSet.hh.

{
  return ComputeCrossSection(dp, elm, mat);
}

References G4VCrossSectionDataSet::ComputeCrossSection().

◆ GetCutEnergy()

G4double G4NeutrinoElectronNcXsc::GetCutEnergy ( )

inline

Definition at line 55 of file G4NeutrinoElectronNcXsc.hh.

55{return fCutEnergy;};

References fCutEnergy.

◆ GetElementCrossSection()

G4double G4NeutrinoElectronNcXsc::GetElementCrossSection	(	const G4DynamicParticle *	aPart,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 of file G4NeutrinoElectronNcXsc.cc.

{
  G4double result = 0., cofL, cofR, cofL2, cofR2, cofLR;
 
  G4double energy = aPart->GetTotalEnergy();
  G4String pName   = aPart->GetDefinition()->GetParticleName();
 
  if( pName == "nu_e")
  {
    cofL = 0.5 + fSin2tW;
    cofR = fSin2tW;
  }
  else if( pName == "anti_nu_e")
  {
    cofL = fSin2tW;
    cofR = 0.5 + fSin2tW;
  }
  else if( pName == "nu_mu")
  {
    cofL = -0.5 + fSin2tW;
    cofR = fSin2tW;
  }
  else if( pName == "anti_nu_mu")
  {
    cofL = fSin2tW;
    cofR = -0.5 + fSin2tW;
  }
  else if( pName == "nu_tau") // vmg: nu_tau as nu_mu ???
  {
    cofL = -0.5 + fSin2tW;
    cofR = fSin2tW;
  }
  else if( pName == "anti_nu_tau")
  {
    cofL = fSin2tW;
    cofR = -0.5 + fSin2tW;
  }
  else
  {
    return result;
  }
  // if( energy <= electron_mass_c2 ) return result;
 
  cofL2 = cofL*cofL;
  cofR2 = cofR*cofR;
  cofLR = cofL*cofR;
 
  if( fCutEnergy > 0. )
  {
    G4double tM  = 2.*energy*energy/(electron_mass_c2 + 2.*energy);
    G4double tM2 = tM*tM;
    G4double tM3 = tM*tM2;
    G4double tC  = fCutEnergy;
    G4double tC2 = tC*tC;
    G4double tC3 = tC*tC2;
 
    result  = (cofL2+cofR2)*(tM-tC);
    result -= (cofR2+cofLR*0.5*electron_mass_c2/energy)*(tM2-tC2)/energy;
    result += cofR2*(tM3-tC3)/energy/energy/3.;
  }
  else
  {
    G4double rtM  = 2.*energy/(electron_mass_c2 + 2.*energy);
    G4double rtM2 = rtM*rtM;
    G4double rtM3 = rtM*rtM2;
 
    result  = (cofL2+cofR2)*rtM*energy;
    result -= (cofR2*energy+cofLR*0.5*electron_mass_c2)*rtM2;
    result += cofR2*rtM3*energy/3.;
  }
  // result = cofL*cofL + cofR*cofR/3.;
  // G4cout<<"cofL2 + cofR2/3. = "<<result<<G4endl;
  // result -= 0.5*cofL*cofR*electron_mass_c2/energy;
 
  G4double aa = 1.;
  G4double bb = 1.7;
  G4double gw = 2.141*GeV;
  G4double dd = 5000.;
  G4double mw = 80.385*GeV;
  G4double mz = 91.1876*GeV;
 
  G4double emass = electron_mass_c2;
  G4double totS  = 2.*energy*emass + emass*emass;
 
  if( energy > 50.*GeV )
  {
    result *= bb;
    result /= 1.+ aa*totS/mz/mz;
 
    if( pName == "anti_nu_e")
    {
      result *= 1. + dd*gw*gw*totS/( (totS-mw*mw)*(totS-mw*mw)+gw*gw*mw*mw );
    }
  }
 
  result *= fCofXsc; //*energy;
 
  result *= ZZ;  // incoherent sum over  all element electrons
 
  result *= fBiasingFactor;
 
  return result;
}

References source.hepunit::electron_mass_c2, G4INCL::KinematicsUtils::energy(), fBiasingFactor, fCofXsc, fCutEnergy, fSin2tW, G4DynamicParticle::GetDefinition(), G4ParticleDefinition::GetParticleName(), G4DynamicParticle::GetTotalEnergy(), and GeV.

Referenced by G4NeutrinoElectronTotXsc::GetElementCrossSection().

◆ GetIsoCrossSection()

G4double G4VCrossSectionDataSet::GetIsoCrossSection	(	const G4DynamicParticle *	dynPart,
		G4int	Z,
		G4int	A,
		const G4Isotope *	iso = `nullptr`,
		const G4Element *	elm = `nullptr`,
		const G4Material *	mat = `nullptr`
	)

virtualinherited

Definition at line 130 of file G4VCrossSectionDataSet.cc.

{
  G4ExceptionDescription ed;
  ed << "GetIsoCrossSection is not implemented in <" << name << ">\n"
     << "Particle: " << dynPart->GetDefinition()->GetParticleName()
     << "  Ekin(MeV)= "  << dynPart->GetKineticEnergy()/MeV;
  if(mat) { ed << "  material: " << mat->GetName(); }
  if(elm) { ed << " element: " << elm->GetName(); }
  ed << " target Z= " << Z << " A= " << A << G4endl;
  G4Exception("G4VCrossSectionDataSet::GetIsoCrossSection", "had001", 
              FatalException, ed);
  return 0.0;
}

References A, FatalException, G4endl, G4Exception(), G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4Element::GetName(), G4Material::GetName(), G4ParticleDefinition::GetParticleName(), MeV, G4VCrossSectionDataSet::name, and Z.

Referenced by G4VCrossSectionDataSet::ComputeCrossSection(), G4GammaNuclearXS::GetIsoCrossSection(), and G4GammaNuclearXS::Initialise().

◆ GetMaxKinEnergy()

G4double G4VCrossSectionDataSet::GetMaxKinEnergy ( ) const

inlineinherited

Definition at line 220 of file G4VCrossSectionDataSet.hh.

{
  return maxKinEnergy;
}

References G4VCrossSectionDataSet::maxKinEnergy.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4LENDCrossSection::IsIsoApplicable(), G4ParticleHPCaptureData::IsIsoApplicable(), G4ParticleHPElasticData::IsIsoApplicable(), G4ParticleHPFissionData::IsIsoApplicable(), and G4ParticleHPInelasticData::IsIsoApplicable().

◆ GetMinKinEnergy()

G4double G4VCrossSectionDataSet::GetMinKinEnergy ( ) const

inlineinherited

Definition at line 210 of file G4VCrossSectionDataSet.hh.

{
  return minKinEnergy;
}

References G4VCrossSectionDataSet::minKinEnergy.

Referenced by G4CrossSectionElastic::IsElementApplicable(), G4CrossSectionInelastic::IsElementApplicable(), G4LENDCrossSection::IsIsoApplicable(), G4ParticleHPCaptureData::IsIsoApplicable(), G4ParticleHPElasticData::IsIsoApplicable(), G4ParticleHPFissionData::IsIsoApplicable(), and G4ParticleHPInelasticData::IsIsoApplicable().

◆ GetName()

const G4String & G4VCrossSectionDataSet::GetName ( ) const

inlineinherited

Definition at line 225 of file G4VCrossSectionDataSet.hh.

{
  return name;
}

References G4VCrossSectionDataSet::name.

Referenced by G4LENDCrossSection::DumpLENDTargetInfo(), G4LENDCrossSection::DumpPhysicsTable(), G4LENDCrossSection::GetIsoCrossSection(), and G4CrossSectionDataStore::PrintCrossSectionHtml().

◆ GetVerboseLevel()

G4int G4VCrossSectionDataSet::GetVerboseLevel ( ) const

inlinevirtualinherited

Reimplemented in G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, and G4ParticleHPInelasticData.

Definition at line 195 of file G4VCrossSectionDataSet.hh.

{
  return verboseLevel;
}

References G4VCrossSectionDataSet::verboseLevel.

◆ IsElementApplicable()

G4bool G4NeutrinoElectronNcXsc::IsElementApplicable	(	const G4DynamicParticle *	aPart,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 67 of file G4NeutrinoElectronNcXsc.cc.

{
  G4bool result  = false;
  G4String pName = aPart->GetDefinition()->GetParticleName();
  G4double minEnergy = 0., energy = aPart->GetTotalEnergy();
  // Z *= 1;
  if( fCutEnergy > 0. ) // min detected recoil electron energy
  {
    minEnergy = 0.5*(fCutEnergy+sqrt(fCutEnergy*(fCutEnergy+2.*electron_mass_c2)));
  }
  if( ( pName == "nu_e"   || pName == "anti_nu_e"   || 
        pName == "nu_mu"  || pName == "anti_nu_mu"  || 
        pName == "nu_tau" || pName == "anti_nu_tau"   ) &&
        energy > minEnergy                                 )
  {
    result = true;
  }
  return result;
}

References source.hepunit::electron_mass_c2, G4INCL::KinematicsUtils::energy(), fCutEnergy, G4DynamicParticle::GetDefinition(), G4ParticleDefinition::GetParticleName(), and G4DynamicParticle::GetTotalEnergy().

Referenced by G4NeutrinoElectronTotXsc::IsElementApplicable().

◆ IsIsoApplicable()

G4bool G4VCrossSectionDataSet::IsIsoApplicable	(	const G4DynamicParticle *	,
		G4int	Z,
		G4int	A,
		const G4Element *	elm = `nullptr`,
		const G4Material *	mat = `nullptr`
	)

virtualinherited

Reimplemented in G4MuNeutrinoNucleusTotXsc, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticXS, G4GammaNuclearXS, G4BGGNucleonElasticXS, G4BGGPionElasticXS, G4BGGPionInelasticXS, G4BGGNucleonInelasticXS, G4ParticleHPBGGNucleonInelasticXS, G4ParticleInelasticXS, G4LENDCrossSection, G4ParticleHPCaptureData, G4ParticleHPElasticData, G4ParticleHPFissionData, G4ParticleHPInelasticData, G4ParticleHPThermalScatteringData, G4ElNeutrinoNucleusTotXsc, G4PhotoNuclearCrossSection, G4ChipsAntiBaryonElasticXS, G4ChipsAntiBaryonInelasticXS, G4ChipsHyperonElasticXS, G4ChipsHyperonInelasticXS, G4ChipsKaonMinusElasticXS, G4ChipsKaonMinusInelasticXS, G4ChipsKaonPlusElasticXS, G4ChipsKaonPlusInelasticXS, G4ChipsKaonZeroElasticXS, G4ChipsKaonZeroInelasticXS, G4ChipsNeutronElasticXS, G4ChipsNeutronInelasticXS, G4ChipsPionMinusElasticXS, G4ChipsPionMinusInelasticXS, G4ChipsPionPlusElasticXS, G4ChipsPionPlusInelasticXS, G4ChipsProtonElasticXS, and G4ChipsProtonInelasticXS.

Definition at line 72 of file G4VCrossSectionDataSet.cc.

{
  return false;
}

Referenced by G4VCrossSectionDataSet::ComputeCrossSection().

◆ SelectIsotope()

const G4Isotope * G4VCrossSectionDataSet::SelectIsotope	(	const G4Element *	anElement,
		G4double	kinEnergy,
		G4double	logE
	)

virtualinherited

Reimplemented in G4GammaNuclearXS, G4NeutronCaptureXS, G4NeutronElasticXS, G4NeutronInelasticXS, and G4ParticleInelasticXS.

Definition at line 149 of file G4VCrossSectionDataSet.cc.

{
  size_t nIso = anElement->GetNumberOfIsotopes();
  const G4Isotope* iso = anElement->GetIsotope(0);
 
  // more than 1 isotope
  if(1 < nIso) {
    const G4double* abundVector = anElement->GetRelativeAbundanceVector();
    G4double sum = 0.0;
    G4double q = G4UniformRand();
    for (size_t j=0; j<nIso; ++j) {
      sum += abundVector[j];
      if(q <= sum) {
    iso = anElement->GetIsotope(j);
    break;
      }
    }
  }
  return iso;
}