#include <G4NeutrinoElectronNcXsc.hh>

Inheritance diagram for G4NeutrinoElectronNcXsc:

Public Member Functions
	G4NeutrinoElectronNcXsc ()

	~G4NeutrinoElectronNcXsc ()

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

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

void	SetCutEnergy (G4double ec)

G4double	GetCutEnergy ()

void	SetBiasingFactor (G4double bf)

G4double	GetBiasingFactor ()

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

virtual G4bool	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material mat=nullptr)

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

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

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

virtual G4double	GetElementCrossSection (const G4DynamicParticle , G4int Z, const G4Material mat=nullptr)

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

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

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) const

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

bool	ForAllAtomsAndEnergies () const

void	SetForAllAtomsAndEnergies (G4bool val)

const G4String &	GetName () const

Protected Attributes
G4double	fCofXsc

G4double	fSin2tW

G4double	fCutEnergy

G4double	fBiasingFactor

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

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

◆ ~G4NeutrinoElectronNcXsc()

G4NeutrinoElectronNcXsc::~G4NeutrinoElectronNcXsc ( )

Definition at line 61 of file G4NeutrinoElectronNcXsc.cc.

62{}

Member Function Documentation

◆ GetBiasingFactor()

G4double G4NeutrinoElectronNcXsc::GetBiasingFactor ( )

inline

Definition at line 58 of file G4NeutrinoElectronNcXsc.hh.

58{return fBiasingFactor;};

◆ GetCutEnergy()

G4double G4NeutrinoElectronNcXsc::GetCutEnergy ( )

inline

Definition at line 55 of file G4NeutrinoElectronNcXsc.hh.

55{return 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;
}

Referenced by G4NeutrinoElectronTotXsc::GetElementCrossSection().

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

Referenced by G4NeutrinoElectronTotXsc::IsElementApplicable().

◆ SetBiasingFactor()

void G4NeutrinoElectronNcXsc::SetBiasingFactor ( G4double bf )

inline

Definition at line 57 of file G4NeutrinoElectronNcXsc.hh.

57{fBiasingFactor=bf;};

Referenced by G4NeutrinoElectronTotXsc::SetBiasingFactor(), and G4NeutrinoElectronTotXsc::SetBiasingFactors().

◆ SetCutEnergy()

void G4NeutrinoElectronNcXsc::SetCutEnergy ( G4double ec )

inline

Definition at line 54 of file G4NeutrinoElectronNcXsc.hh.

54{fCutEnergy=ec;};

Member Data Documentation

◆ fBiasingFactor

G4double G4NeutrinoElectronNcXsc::fBiasingFactor

protected

Definition at line 65 of file G4NeutrinoElectronNcXsc.hh.

Referenced by G4NeutrinoElectronNcXsc(), GetBiasingFactor(), GetElementCrossSection(), and SetBiasingFactor().

◆ fCofXsc

G4double G4NeutrinoElectronNcXsc::fCofXsc

protected

Definition at line 62 of file G4NeutrinoElectronNcXsc.hh.

Referenced by G4NeutrinoElectronNcXsc(), and GetElementCrossSection().

◆ fCutEnergy

G4double G4NeutrinoElectronNcXsc::fCutEnergy

protected

Definition at line 64 of file G4NeutrinoElectronNcXsc.hh.

Referenced by G4NeutrinoElectronNcXsc(), GetCutEnergy(), GetElementCrossSection(), IsElementApplicable(), and SetCutEnergy().

◆ fSin2tW

G4double G4NeutrinoElectronNcXsc::fSin2tW

protected

Definition at line 63 of file G4NeutrinoElectronNcXsc.hh.

Referenced by G4NeutrinoElectronNcXsc(), and GetElementCrossSection().

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

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4NeutrinoElectronNcXsc()

◆ ~G4NeutrinoElectronNcXsc()

Member Function Documentation

◆ GetBiasingFactor()

◆ GetCutEnergy()

◆ GetElementCrossSection()

◆ IsElementApplicable()

◆ SetBiasingFactor()

◆ SetCutEnergy()

Member Data Documentation

◆ fBiasingFactor

◆ fCofXsc

◆ fCutEnergy

◆ fSin2tW