#include <G4NeutrinoElectronCcXsc.hh>

Inheritance diagram for G4NeutrinoElectronCcXsc:

Public Member Functions
	G4NeutrinoElectronCcXsc ()

	~G4NeutrinoElectronCcXsc ()

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	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	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition , const G4Element , 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 G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition , 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 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

void	SetName (const G4String &nam)

G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete

	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Protected Attributes
G4double	fCofXsc

G4double	fSin2tW

G4double	fCutEnergy

G4double	fBiasingFactor

G4ParticleDefinition *	theMuonMinus

G4ParticleDefinition *	theTauMinus

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel {0}

G4String	name

Detailed Description

Definition at line 41 of file G4NeutrinoElectronCcXsc.hh.

Constructor & Destructor Documentation

◆ G4NeutrinoElectronCcXsc()

G4NeutrinoElectronCcXsc::G4NeutrinoElectronCcXsc ( )

Definition at line 43 of file G4NeutrinoElectronCcXsc.cc.

 : G4VCrossSectionDataSet("NuElectronCcXsc")
{
  // 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<<"fCofXsc = "<<fCofXsc*GeV/cm2<<" cm2/GeV"<<G4endl;
 
  // G4cout<<"hbarc = "<<hbarc/MeV/fermi<<" MeV*fermi"<<G4endl;
 
  // PDG2016: sin^2 theta Weinberg
 
  fSin2tW = 0.23129; // 0.2312;
 
  fCutEnergy = 0.; // default value
 
  fBiasingFactor = 1.; // default as physics
 
  theMuonMinus = G4MuonMinus::MuonMinus(); 
  theTauMinus  = G4TauMinus::TauMinus(); 
}

◆ ~G4NeutrinoElectronCcXsc()

G4NeutrinoElectronCcXsc::~G4NeutrinoElectronCcXsc ( )

Definition at line 69 of file G4NeutrinoElectronCcXsc.cc.

70{}

Member Function Documentation

◆ GetBiasingFactor()

G4double G4NeutrinoElectronCcXsc::GetBiasingFactor ( )

inline

Definition at line 60 of file G4NeutrinoElectronCcXsc.hh.

60{return fBiasingFactor;};

◆ GetCutEnergy()

G4double G4NeutrinoElectronCcXsc::GetCutEnergy ( )

inline

Definition at line 57 of file G4NeutrinoElectronCcXsc.hh.

57{return fCutEnergy;};

◆ GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 99 of file G4NeutrinoElectronCcXsc.cc.

{
  G4double result = 0., totS, fmass, fmass2, emass=electron_mass_c2, emass2;
 
  G4double energy = aPart->GetTotalEnergy();
  G4String pName   = aPart->GetDefinition()->GetParticleName();
 
  emass2 = emass*emass;
  totS   = 2.*energy*emass + emass2;
 
  if( pName == "anti_nu_e" || pName == "nu_mu")
  {
    fmass  = theMuonMinus->GetPDGMass();
    fmass2 = fmass*fmass;
    result = (1. - fmass2/totS)*(1. - fmass2/totS);
  }
  else if( pName == "anti_nu_mu")
  {
    fmass  = theMuonMinus->GetPDGMass();
    fmass2 = fmass*fmass;
 
    result  = (1.+ emass2/totS)*(1.+ fmass2/totS);
    result += (1.- emass2/totS)*(1.- fmass2/totS)/3.;
    result *= 0.25*(1. - fmass2/totS)*(1. - fmass2/totS);
  }
  else if( pName == "nu_tau") 
  {
    fmass  = theTauMinus->GetPDGMass();
    fmass2 = fmass*fmass;
    result = (1. - fmass2/totS)*(1. - fmass2/totS);
  }
  else if( pName == "anti_nu_tau")
  {
    fmass  = theTauMinus->GetPDGMass();
    fmass2 = fmass*fmass;
 
    result  = (1.+ emass2/totS)*(1.+ fmass2/totS);
    result += (1.- emass2/totS)*(1.- fmass2/totS)/3.;
    result *= 0.25*(1. - fmass2/totS)*(1. - fmass2/totS);
  }
  else
  {
    return result;
  }
  // if( energy <= electron_mass_c2 ) return result;
 
  G4double aa = 1.;
  G4double bb = 1.7;
  G4double gw = 2.141*GeV;
  G4double dd = 5000.;
  G4double mw = 80.385*GeV;
 
  if( energy > 50.*GeV )
  {
    result *= bb;
    result /= 1.+ aa*totS/mw/mw;
 
    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 *= energy + 0.5*emass;
  result *= ZZ;  // incoherent sum over  all element electrons
 
  result *= fBiasingFactor; // biasing up, if set >1
 
  return result;
}

◆ IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 75 of file G4NeutrinoElectronCcXsc.cc.

{
  G4bool result  = false;
  G4String pName = aPart->GetDefinition()->GetParticleName();
  G4double minEnergy = 0., energy = aPart->GetTotalEnergy();
  G4double fmass, emass = electron_mass_c2;
 
  if(    pName == "anti_nu_e" ||   pName == "nu_mu"   || pName == "anti_nu_mu"  ) fmass = theMuonMinus->GetPDGMass(); 
  else if( pName == "nu_tau"  || pName == "anti_nu_tau" ) fmass = theTauMinus->GetPDGMass(); 
  else fmass = emass;
 
  minEnergy = (fmass-emass)*(fmass+emass)/emass;
  
  if( ( pName == "nu_mu"  || pName == "anti_nu_mu"  || 
        pName == "nu_tau" || pName == "anti_nu_tau"   ) &&
        energy > minEnergy                                 )
  {
    result = true;
  }
  return result;
}

◆ SetBiasingFactor()

void G4NeutrinoElectronCcXsc::SetBiasingFactor ( G4double bf )

inline

Definition at line 59 of file G4NeutrinoElectronCcXsc.hh.

59{fBiasingFactor=bf;};

◆ SetCutEnergy()

void G4NeutrinoElectronCcXsc::SetCutEnergy ( G4double ec )

inline

Definition at line 56 of file G4NeutrinoElectronCcXsc.hh.

56{fCutEnergy=ec;};

Member Data Documentation

◆ fBiasingFactor

G4double G4NeutrinoElectronCcXsc::fBiasingFactor

protected

Definition at line 67 of file G4NeutrinoElectronCcXsc.hh.

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

◆ fCofXsc

G4double G4NeutrinoElectronCcXsc::fCofXsc

protected

Definition at line 64 of file G4NeutrinoElectronCcXsc.hh.

Referenced by G4NeutrinoElectronCcXsc(), and GetElementCrossSection().

◆ fCutEnergy

G4double G4NeutrinoElectronCcXsc::fCutEnergy

protected

Definition at line 66 of file G4NeutrinoElectronCcXsc.hh.

Referenced by G4NeutrinoElectronCcXsc(), GetCutEnergy(), and SetCutEnergy().

◆ fSin2tW

G4double G4NeutrinoElectronCcXsc::fSin2tW

protected

Definition at line 65 of file G4NeutrinoElectronCcXsc.hh.

Referenced by G4NeutrinoElectronCcXsc().

◆ theMuonMinus

G4ParticleDefinition* G4NeutrinoElectronCcXsc::theMuonMinus

protected

Definition at line 69 of file G4NeutrinoElectronCcXsc.hh.

Referenced by G4NeutrinoElectronCcXsc(), GetElementCrossSection(), and IsElementApplicable().

◆ theTauMinus

G4ParticleDefinition* G4NeutrinoElectronCcXsc::theTauMinus

protected

Definition at line 70 of file G4NeutrinoElectronCcXsc.hh.

Referenced by G4NeutrinoElectronCcXsc(), GetElementCrossSection(), and IsElementApplicable().

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

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4NeutrinoElectronCcXsc()

◆ ~G4NeutrinoElectronCcXsc()

Member Function Documentation

◆ GetBiasingFactor()

◆ GetCutEnergy()

◆ GetElementCrossSection()

◆ IsElementApplicable()

◆ SetBiasingFactor()

◆ SetCutEnergy()

Member Data Documentation

◆ fBiasingFactor

◆ fCofXsc

◆ fCutEnergy

◆ fSin2tW

◆ theMuonMinus

◆ theTauMinus