G4ElNucleusSFcs Class Reference

#include <G4ElNucleusSFcs.hh>

Inheritance diagram for G4ElNucleusSFcs:

Public Member Functions
	G4ElNucleusSFcs ()
virtual	~G4ElNucleusSFcs ()
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4bool	IsElementApplicable (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	ThresholdEnergy ()
G4double	GetRatio (G4int Z, G4int A)
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	GetElementCrossSection (const G4DynamicParticle *, G4int Z, 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	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

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel
G4String	name

Detailed Description

Definition at line 51 of file G4ElNucleusSFcs.hh.

Constructor & Destructor Documentation

G4ElNucleusSFcs()

G4ElNucleusSFcs::G4ElNucleusSFcs

(

)

Definition at line 42 of file G4ElNucleusSFcs.cc.

                                :G4VCrossSectionDataSet(Default_Name())
{
  fCHIPScs = new G4ElectroNuclearCrossSection();  
}

~G4ElNucleusSFcs()

G4ElNucleusSFcs::~G4ElNucleusSFcs ( )

virtual

Definition at line 47 of file G4ElNucleusSFcs.cc.

{
  if( fCHIPScs != nullptr ) delete fCHIPScs;
}

Member Function Documentation

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 53 of file G4ElNucleusSFcs.cc.

{
    outFile << "G4ElNucleusSFcs provides the inelastic\n"
    << "cross section for e- and e+ interactions with nuclei according to the structure function approach."
    << "all energies.\n";
}

Default_Name()

static const char * G4ElNucleusSFcs::Default_Name ( )

inlinestatic

Definition at line 58 of file G4ElNucleusSFcs.hh.

{return "ElectronNucleusSFcs";}

GetIsoCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 70 of file G4ElNucleusSFcs.cc.

{
  G4double xsc(0.), ratio(1.);
 
  xsc = fCHIPScs->GetElementCrossSection( aPart, ZZ, nullptr ); //, const G4Material*);
 
  ratio = GetRatio(ZZ,AA);
 
  if( ratio > 0. ) xsc /= ratio;
  
  return xsc;
}

GetRatio()

G4double G4ElNucleusSFcs::GetRatio	(	G4int	Z,
		G4int	A )

Definition at line 88 of file G4ElNucleusSFcs.cc.

{
  G4double ratio(1.);
 
  if      ( Z == 1   && A == 1 )   return 1.51;
  else if ( Z == 1   && A == 2 )   return 0.33;
  else if ( Z == 1   && A == 3 )   return 0.27;
  else if ( Z == 2   && A == 4 )   return 1.81;
  else if ( Z == 6   && A == 12 )  return 2.26;
  else if ( Z == 7   && A == 14 )  return 2.47;
  else if ( Z == 8   && A == 16 )  return 2.61;
  else if ( Z == 13  && A == 27 )  return 2.57;
  else if ( Z == 14  && A == 28 )  return 2.49;
  else if ( Z == 18  && A == 40 )  return 2.72;
  else if ( Z == 22  && A == 48 )  return 2.71;
  else if ( Z == 26  && A == 56 )  return 2.79;
  else if ( Z == 29  && A == 64 )  return 2.78;
  else if ( Z == 32  && A == 73 )  return 2.87;
  else if ( Z == 42  && A == 96 )  return 3.02;
  else if ( Z == 46  && A == 106 ) return 3.02;
  else if ( Z == 47  && A == 108 ) return 2.99;
  else if ( Z == 48  && A == 112 ) return 3.00;
  else if ( Z == 74  && A == 184 ) return 3.44;
  else if ( Z == 79  && A == 200 ) return 3.49;
  else if ( Z == 82  && A == 207 ) return 3.48;
  else if ( Z == 92  && A == 238 ) return 3.88;
  else
  {
    G4int it(0), iMax(19);
    G4double zz = G4double(Z);
 
    for ( it = 0; it < iMax; ++it ) if ( zz <= fZZ[it] ) break;
 
    if     ( it == 0 )    return fRR[0];
    else if( it == iMax ) return fRR[iMax-1];
    else
    {
      G4double x1 = fZZ[it-1];
      G4double x2 = fZZ[it];
      G4double y1 = fRR[it-1];
      G4double y2 = fRR[it];
      
      if( x1 >= x2 ) return fRR[it];
      else
      {
        G4double angle = (y2-y1)/(x2-x1);
    ratio = y1 + ( zz - x1 )*angle;
      }
    } 
  }
  return ratio;
}

Referenced by GetIsoCrossSection().

IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 60 of file G4ElNucleusSFcs.cc.

{
 
  G4double eKin   = aPart->GetKineticEnergy();
  G4double thTkin = ThresholdEnergy();
  
  return ( Z > 0 && Z < 120 && eKin > thTkin );
}

ThresholdEnergy()

G4double G4ElNucleusSFcs::ThresholdEnergy

(

)

Definition at line 143 of file G4ElNucleusSFcs.cc.

{
  G4double thTkin = 134.9766*CLHEP::MeV; // Pi0 threshold for the nucleon
  
  thTkin *= 1.3; // nucleon motion
  
  return thTkin;
}

Referenced by IsElementApplicable().

---

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

• G4ElNucleusSFcs.hh
• G4ElNucleusSFcs.cc