G4NuclNuclDiffuseElastic Class Reference

#include <G4NuclNuclDiffuseElastic.hh>

Inheritance diagram for G4NuclNuclDiffuseElastic:

Public Member Functions
void	ActivateFor (const G4Element *anElement)
void	ActivateFor (const G4Material *aMaterial)
G4complex	Amplitude (G4double theta)
G4complex	AmplitudeFar (G4double theta)
G4complex	AmplitudeGG (G4double theta)
G4double	AmplitudeGGMod2 (G4double theta)
G4complex	AmplitudeGla (G4double theta)
G4double	AmplitudeGlaMod2 (G4double theta)
G4double	AmplitudeMod2 (G4double theta)
G4complex	AmplitudeNear (G4double theta)
G4complex	AmplitudeSim (G4double theta)
G4double	AmplitudeSimMod2 (G4double theta)
G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus) override
G4double	BesselJone (G4double z)
G4double	BesselJzero (G4double z)
G4double	BesselOneByArg (G4double z)
void	BuildAngleTable ()
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
G4double	CalcMandelstamS (const G4double mp, const G4double mt, const G4double Plab)
G4double	CalculateAm (G4double momentum, G4double n, G4double Z)
G4double	CalculateCoulombPhase (G4int n)
void	CalculateCoulombPhaseZero ()
G4double	CalculateNuclearRad (G4double A)
G4double	CalculateParticleBeta (const G4ParticleDefinition *particle, G4double momentum)
void	CalculateRutherfordAnglePar ()
G4double	CalculateZommerfeld (G4double beta, G4double Z1, G4double Z2)
G4double	ComputeMomentumCMS (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4complex	CoulombAmplitude (G4double theta)
G4double	CoulombAmplitudeMod2 (G4double theta)
G4double	DampFactor (G4double z)
void	DeActivateFor (const G4Element *anElement)
void	DeActivateFor (const G4Material *aMaterial)
	G4NuclNuclDiffuseElastic ()
G4complex	GammaLess (G4double theta)
G4complex	GammaLogarithm (G4complex xx)
G4complex	GammaLogB2n (G4complex xx)
G4complex	GammaMore (G4double theta)
G4double	GetCint (G4double x)
G4double	GetCofAlphaCoulomb ()
G4double	GetCofAlphaMax ()
G4double	GetCosHaPit2 (G4double t)
G4double	GetCoulombElasticXsc (const G4ParticleDefinition *particle, G4double theta, G4double momentum, G4double Z)
G4double	GetCoulombIntegralXsc (const G4ParticleDefinition *particle, G4double momentum, G4double Z, G4double theta1, G4double theta2)
G4double	GetCoulombTotalXsc (const G4ParticleDefinition *particle, G4double momentum, G4double Z)
G4double	GetDiffElasticProb (G4double theta)
G4double	GetDiffElasticSumProb (G4double theta)
G4double	GetDiffElasticSumProbA (G4double alpha)
G4double	GetDiffuseElasticSumXsc (const G4ParticleDefinition *particle, G4double theta, G4double momentum, G4double A, G4double Z)
G4double	GetDiffuseElasticXsc (const G4ParticleDefinition *particle, G4double theta, G4double momentum, G4double A)
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
G4double	GetErf (G4double x)
G4complex	GetErfcComp (G4complex z, G4int nMax)
G4complex	GetErfcInt (G4complex z)
G4complex	GetErfComp (G4complex z, G4int nMax)
G4complex	GetErfcSer (G4complex z, G4int nMax)
G4complex	GetErfInt (G4complex z)
G4complex	GetErfSer (G4complex z, G4int nMax)
G4double	GetExpCos (G4double x)
G4double	GetExpSin (G4double x)
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
G4double	GetFresnelDiffuseXsc (G4double theta)
G4double	GetFresnelIntegrandXsc (G4double alpha)
G4double	GetHadronNucleonXscNS (G4ParticleDefinition *pParticle, G4double pTkin, G4ParticleDefinition *tParticle)
G4double	GetIntegrandFunction (G4double theta)
G4double	GetInvCoulombElasticXsc (const G4ParticleDefinition *particle, G4double tMand, G4double momentum, G4double A, G4double Z)
G4double	GetInvElasticSumXsc (const G4ParticleDefinition *particle, G4double tMand, G4double momentum, G4double A, G4double Z)
G4double	GetInvElasticXsc (const G4ParticleDefinition *particle, G4double theta, G4double momentum, G4double A, G4double Z)
G4double	GetLegendrePol (G4int n, G4double x)
G4double	GetMaxEnergy () const
G4double	GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
G4double	GetMinEnergy () const
G4double	GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
const G4String &	GetModelName () const
G4double	GetNuclearRadius ()
G4double	GetProfileLambda ()
G4double	GetRatioGen (G4double theta)
G4double	GetRatioSim (G4double theta)
G4double	GetRecoilEnergyThreshold () const
G4double	GetRutherfordXsc (G4double theta)
G4double	GetScatteringAngle (G4int iMomentum, G4int iAngle, G4double position)
G4double	GetSinHaPit2 (G4double t)
G4double	GetSint (G4double x)
G4double	GetSlopeCof (const G4int pdg)
G4int	GetVerboseLevel () const
void	InitDynParameters (const G4ParticleDefinition *theParticle, G4double partMom)
void	Initialise ()
virtual void	InitialiseModel ()
void	InitialiseOnFly (G4double Z, G4double A)
void	InitParameters (const G4ParticleDefinition *theParticle, G4double partMom, G4double Z, G4double A)
void	InitParametersGla (const G4DynamicParticle *aParticle, G4double partMom, G4double Z, G4double A)
G4double	IntegralElasticProb (const G4ParticleDefinition *particle, G4double theta, G4double momentum, G4double A)
virtual G4bool	IsApplicable (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
G4bool	IsBlocked (const G4Element *anElement) const
G4bool	IsBlocked (const G4Material *aMaterial) const
G4double	LowestEnergyLimit () const
void	ModelDescription (std::ostream &) const override
G4bool	operator!= (const G4HadronicInteraction &right) const =delete
G4bool	operator== (const G4HadronicInteraction &right) const =delete
G4complex	PhaseFar (G4double theta)
G4complex	PhaseNear (G4double theta)
G4double	Profile (G4double theta)
G4double	ProfileFar (G4double theta)
G4double	ProfileNear (G4double theta)
G4double	SampleCoulombMuCMS (const G4ParticleDefinition *aParticle, G4double p)
virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4double	SampleT (const G4ParticleDefinition *aParticle, G4double p, G4double A)
G4double	SampleTableT (const G4ParticleDefinition *aParticle, G4double p, G4double Z, G4double A)
G4double	SampleTableThetaCMS (const G4ParticleDefinition *aParticle, G4double p, G4double Z, G4double A)
G4double	SampleThetaCMS (const G4ParticleDefinition *aParticle, G4double p, G4double A)
G4double	SampleThetaLab (const G4HadProjectile *aParticle, G4double tmass, G4double A)
void	SetCofAlpha (G4double pa)
void	SetCofAlphaCoulomb (G4double pa)
void	SetCofAlphaMax (G4double pa)
void	SetCofDelta (G4double pa)
void	SetCofFar (G4double pa)
void	SetCofLambda (G4double pa)
void	SetCofPhase (G4double pa)
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
void	SetEtaRatio (G4double pa)
void	SetHEModelLowLimit (G4double value)
void	SetLowestEnergyLimit (G4double value)
void	SetMaxEnergy (const G4double anEnergy)
void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
void	SetMaxL (G4int l)
void	SetMinEnergy (G4double anEnergy)
void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
void	SetNuclearRadiusCof (G4double r)
void	SetPlabLowLimit (G4double value)
void	SetProfileAlpha (G4double pa)
void	SetProfileDelta (G4double pd)
void	SetProfileLambda (G4double pl)
void	SetQModelLowLimit (G4double value)
void	SetRecoilEnergyThreshold (G4double val)
void	SetRecoilKinEnergyLimit (G4double value)
void	SetVerboseLevel (G4int value)
void	TestAngleTable (const G4ParticleDefinition *theParticle, G4double partMom, G4double Z, G4double A)
G4complex	TestErfcComp (G4complex z, G4int nMax)
G4complex	TestErfcInt (G4complex z)
G4complex	TestErfcSer (G4complex z, G4int nMax)
G4double	ThetaCMStoThetaLab (const G4DynamicParticle *aParticle, G4double tmass, G4double thetaCMS)
G4double	ThetaLabToThetaCMS (const G4DynamicParticle *aParticle, G4double tmass, G4double thetaLab)
virtual	~G4NuclNuclDiffuseElastic ()

Protected Member Functions
void	Block ()
G4bool	IsBlocked () const
void	SetModelName (const G4String &nam)

Protected Attributes
G4bool	isBlocked
G4double	pLocalTmax
G4int	secID
G4double	theMaxEnergy
G4double	theMinEnergy
G4HadFinalState	theParticleChange
G4int	verboseLevel

Private Attributes
std::pair< G4double, G4double >	epCheckLevels
G4bool	fAddCoulomb
G4double	fAm
std::vector< G4PhysicsTable * >	fAngleBank
G4int	fAngleBin
G4PhysicsTable *	fAngleTable
G4double	fAtomicNumber
G4double	fAtomicWeight
G4double	fBeta
G4double	fCofAlpha
G4double	fCofAlphaCoulomb
G4double	fCofAlphaMax
G4double	fCofDelta
G4double	fCofFar
G4double	fCofLambda
G4double	fCofPhase
G4double	fCoulombMuC
G4double	fCoulombPhase0
std::vector< G4String >	fElementNameVector
std::vector< G4double >	fElementNumberVector
G4int	fEnergyBin
G4PhysicsLogVector *	fEnergyVector
G4double	fEtaRatio
G4double	fHalfRutThetaTg
G4double	fHalfRutThetaTg2
G4int	fMaxL
G4double	fNuclearRadius
G4double	fNuclearRadius1
G4double	fNuclearRadius2
G4double	fNuclearRadiusCof
G4double	fNuclearRadiusSquare
const G4ParticleDefinition *	fParticle
G4double	fProfileAlpha
G4double	fProfileDelta
G4double	fProfileLambda
G4double	fReZ
G4double	fRutherfordRatio
G4double	fRutherfordTheta
G4double	fSumSigma
G4double	fWaveVector
G4double	fZommerfeld
G4double	lowEnergyLimitHE
G4double	lowEnergyLimitQ
G4double	lowEnergyRecoilLimit
G4double	lowestEnergyLimit
G4int	nwarn
G4double	plabLowLimit
G4double	recoilEnergyThreshold
G4HadronicInteractionRegistry *	registry
G4ParticleDefinition *	theAlpha
std::vector< const G4Material * >	theBlockedList
std::vector< const G4Element * >	theBlockedListElements
G4ParticleDefinition *	theDeuteron
std::vector< std::pair< G4double, const G4Material * > >	theMaxEnergyList
std::vector< std::pair< G4double, const G4Element * > >	theMaxEnergyListElements
std::vector< std::pair< G4double, const G4Material * > >	theMinEnergyList
std::vector< std::pair< G4double, const G4Element * > >	theMinEnergyListElements
G4String	theModelName
G4ParticleDefinition *	theNeutron
const G4ParticleDefinition *	thePionMinus
const G4ParticleDefinition *	thePionPlus
G4ParticleDefinition *	theProton

Detailed Description

Definition at line 61 of file G4NuclNuclDiffuseElastic.hh.

Constructor & Destructor Documentation

G4NuclNuclDiffuseElastic()

G4NuclNuclDiffuseElastic::G4NuclNuclDiffuseElastic

(

)

Definition at line 68 of file G4NuclNuclDiffuseElastic.cc.

  : G4HadronElastic("NNDiffuseElastic"), fParticle(0)
{
  SetMinEnergy( 50*MeV );
  SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
  verboseLevel = 0;
  lowEnergyRecoilLimit = 100.*keV;  
  lowEnergyLimitQ  = 0.0*GeV;  
  lowEnergyLimitHE = 0.0*GeV;  
  lowestEnergyLimit= 0.0*keV;  
  plabLowLimit     = 20.0*MeV;
 
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  theDeuteron = G4Deuteron::Deuteron();
  theAlpha    = G4Alpha::Alpha();
  thePionPlus = G4PionPlus::PionPlus();
  thePionMinus= G4PionMinus::PionMinus();
 
  fEnergyBin = 300;  // Increased from the original 200 to have no wider log-energy-bins up to 10 PeV 
  fAngleBin  = 200;
 
  fEnergyVector =  new G4PhysicsLogVector( theMinEnergy, theMaxEnergy, fEnergyBin );
  fAngleTable = 0;
 
  fParticle = 0;
  fWaveVector = 0.;
  fAtomicWeight = 0.;
  fAtomicNumber = 0.;
  fNuclearRadius = 0.;
  fBeta = 0.;
  fZommerfeld = 0.;
  fAm = 0.;
  fAddCoulomb = false;
  // Ranges of angle table relative to current Rutherford (Coulomb grazing) angle
 
  // Empirical parameters
 
  fCofAlphaMax     = 1.5;
  fCofAlphaCoulomb = 0.5;
 
  fProfileDelta  = 1.;
  fProfileAlpha  = 0.5;
 
  fCofLambda = 1.0;
  fCofDelta  = 0.04;
  fCofAlpha  = 0.095;
 
  fNuclearRadius1 = fNuclearRadius2 = fNuclearRadiusSquare 
    = fRutherfordRatio = fCoulombPhase0 = fHalfRutThetaTg = fHalfRutThetaTg2 
    = fRutherfordTheta = fProfileLambda = fCofPhase = fCofFar
    = fSumSigma = fEtaRatio = fReZ = 0.0;
  fMaxL = 0;
 
  fNuclearRadiusCof = 1.0;
  fCoulombMuC = 0.0;
}

References G4Alpha::Alpha(), G4Deuteron::Deuteron(), fAddCoulomb, fAm, fAngleBin, fAngleTable, fAtomicNumber, fAtomicWeight, fBeta, fCofAlpha, fCofAlphaCoulomb, fCofAlphaMax, fCofDelta, fCofFar, fCofLambda, fCofPhase, fCoulombMuC, fCoulombPhase0, fEnergyBin, fEnergyVector, fEtaRatio, fHalfRutThetaTg, fHalfRutThetaTg2, fMaxL, fNuclearRadius, fNuclearRadius1, fNuclearRadius2, fNuclearRadiusCof, fNuclearRadiusSquare, fParticle, fProfileAlpha, fProfileDelta, fProfileLambda, fReZ, fRutherfordRatio, fRutherfordTheta, fSumSigma, fWaveVector, fZommerfeld, G4HadronicInteraction::GetMaxEnergy(), GeV, G4HadronicParameters::Instance(), keV, lowEnergyLimitHE, lowEnergyLimitQ, lowEnergyRecoilLimit, lowestEnergyLimit, MeV, G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), plabLowLimit, G4Proton::Proton(), G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), theAlpha, theDeuteron, G4HadronicInteraction::theMaxEnergy, G4HadronicInteraction::theMinEnergy, theNeutron, thePionMinus, thePionPlus, theProton, and G4HadronicInteraction::verboseLevel.

Referenced by BuildAngleTable(), GetCint(), GetErfInt(), GetSint(), IntegralElasticProb(), SampleThetaCMS(), and TestAngleTable().

~G4NuclNuclDiffuseElastic()

G4NuclNuclDiffuseElastic::~G4NuclNuclDiffuseElastic ( )

virtual

Definition at line 130 of file G4NuclNuclDiffuseElastic.cc.

{
  if ( fEnergyVector ) {
    delete fEnergyVector;
    fEnergyVector = 0;
  }
 
  for ( std::vector<G4PhysicsTable*>::iterator it = fAngleBank.begin();
        it != fAngleBank.end(); ++it ) {
    if ( (*it) ) (*it)->clearAndDestroy();
    delete *it;
    *it = 0;
  }
  fAngleTable = 0;
}

References fAngleBank, fAngleTable, and fEnergyVector.

Member Function Documentation

ActivateFor() [1/2]

void G4HadronicInteraction::ActivateFor ( const G4Element *  anElement )

inlineinherited

Definition at line 128 of file G4HadronicInteraction.hh.

  { 
    Block(); 
    SetMaxEnergy(GetMaxEnergy(), anElement);
    SetMinEnergy(GetMinEnergy(), anElement);
  }

References G4HadronicInteraction::Block(), G4HadronicInteraction::GetMaxEnergy(), G4HadronicInteraction::GetMinEnergy(), G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().

ActivateFor() [2/2]

void G4HadronicInteraction::ActivateFor ( const G4Material *  aMaterial )

inlineinherited

Definition at line 120 of file G4HadronicInteraction.hh.

  { 
    Block(); 
    SetMaxEnergy(GetMaxEnergy(), aMaterial);
    SetMinEnergy(GetMinEnergy(), aMaterial);
  }

References G4HadronicInteraction::Block(), G4HadronicInteraction::GetMaxEnergy(), G4HadronicInteraction::GetMinEnergy(), G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().

Amplitude()

G4complex G4NuclNuclDiffuseElastic::Amplitude ( G4double  theta )

inline

Definition at line 978 of file G4NuclNuclDiffuseElastic.hh.

{
 
  G4complex out = AmplitudeNear(theta) + fCofFar*AmplitudeFar(theta);
  // G4complex out = AmplitudeNear(theta);
  // G4complex out = AmplitudeFar(theta);
  return    out;
}

References AmplitudeFar(), AmplitudeNear(), and fCofFar.

Referenced by AmplitudeMod2().

AmplitudeFar()

G4complex G4NuclNuclDiffuseElastic::AmplitudeFar ( G4double  theta )

inline

Definition at line 964 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double kappa = std::sqrt(0.5*fProfileLambda/std::sin(theta)/CLHEP::pi);
  G4complex out = G4complex(kappa/fWaveVector,0.);
  out *= ProfileFar(theta);
  out *= PhaseFar(theta);
  return out;
}

References fProfileLambda, fWaveVector, PhaseFar(), CLHEP::pi, and ProfileFar().

Referenced by Amplitude().

AmplitudeGG()

G4complex G4NuclNuclDiffuseElastic::AmplitudeGG

(

G4double

theta

)

Definition at line 1691 of file G4NuclNuclDiffuseElastic.cc.

{
  G4int n;
  G4double T12b, a, aTemp, b2, sinThetaH = std::sin(0.5*theta);
  G4double sinThetaH2 = sinThetaH*sinThetaH;
  G4complex out = G4complex(0.,0.); 
  G4complex im  = G4complex(0.,1.);
 
  a  = -fSumSigma/CLHEP::twopi/fNuclearRadiusSquare;
  b2 = fWaveVector*fWaveVector*fNuclearRadiusSquare*sinThetaH2;
 
  aTemp = a;
 
  for( n = 1; n < fMaxL; n++)
  {    
    T12b   = aTemp*G4Exp(-b2/n)/n;         
    aTemp *= a;  
    out   += T12b; 
    G4cout<<"out = "<<out<<G4endl;  
  }
  out *= -4.*im*fWaveVector/CLHEP::pi;
  out += CoulombAmplitude(theta);
  return out;
}

References CoulombAmplitude(), fMaxL, fNuclearRadiusSquare, fSumSigma, fWaveVector, G4cout, G4endl, G4Exp(), CLHEP::detail::n, CLHEP::pi, and CLHEP::twopi.

Referenced by AmplitudeGGMod2().

AmplitudeGGMod2()

G4double G4NuclNuclDiffuseElastic::AmplitudeGGMod2 ( G4double  theta )

inline

Definition at line 1069 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex out = AmplitudeGG(theta);
  G4double mod2 = out.real()*out.real() + out.imag()*out.imag();
  return   mod2;
}

References AmplitudeGG().

AmplitudeGla()

G4complex G4NuclNuclDiffuseElastic::AmplitudeGla

(

G4double

theta

)

Definition at line 1664 of file G4NuclNuclDiffuseElastic.cc.

{
  G4int n;
  G4double T12b, b, b2; // cosTheta = std::cos(theta);
  G4complex out = G4complex(0.,0.), shiftC, shiftN; 
  G4complex im  = G4complex(0.,1.);
 
  for( n = 0; n < fMaxL; n++)
  {
    shiftC = std::exp( im*2.*CalculateCoulombPhase(n) );
    // b = ( fZommerfeld + std::sqrt( fZommerfeld*fZommerfeld + n*(n+1) ) )/fWaveVector;
    b = ( std::sqrt( G4double(n*(n+1)) ) )/fWaveVector;
    b2 = b*b;
    T12b = fSumSigma*G4Exp(-b2/fNuclearRadiusSquare)/CLHEP::pi/fNuclearRadiusSquare;         
    shiftN = std::exp( -0.5*(1.-im*fEtaRatio)*T12b ) - 1.;
    out +=  (2.*n+1.)*shiftC*shiftN*GetLegendrePol(n, theta);   
  }
  out /= 2.*im*fWaveVector;
  out += CoulombAmplitude(theta);
  return out;
}

References CalculateCoulombPhase(), CoulombAmplitude(), fEtaRatio, fMaxL, fNuclearRadiusSquare, fSumSigma, fWaveVector, G4Exp(), GetLegendrePol(), CLHEP::detail::n, and CLHEP::pi.

Referenced by AmplitudeGlaMod2().

AmplitudeGlaMod2()

G4double G4NuclNuclDiffuseElastic::AmplitudeGlaMod2 ( G4double  theta )

inline

Definition at line 1058 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex out = AmplitudeGla(theta);
  G4double mod2 = out.real()*out.real() + out.imag()*out.imag();
  return   mod2;
}

References AmplitudeGla().

AmplitudeMod2()

G4double G4NuclNuclDiffuseElastic::AmplitudeMod2 ( G4double  theta )

inline

Definition at line 991 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex out = Amplitude(theta);
  G4double mod2 = out.real()*out.real() + out.imag()*out.imag();
  return   mod2;
}

References Amplitude().

AmplitudeNear()

G4complex G4NuclNuclDiffuseElastic::AmplitudeNear

(

G4double

theta

)

Definition at line 1608 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double kappa = std::sqrt(0.5*fProfileLambda/std::sin(theta)/CLHEP::pi);
  G4complex out = G4complex(kappa/fWaveVector,0.);
 
  out *= PhaseNear(theta);
 
  if( theta <= fRutherfordTheta )
  {
    out *= GammaLess(theta) + ProfileNear(theta);
    // out *= GammaMore(theta) + ProfileNear(theta);
    out += CoulombAmplitude(theta);
  }
  else
  {
    out *= GammaMore(theta) + ProfileNear(theta);
    // out *= GammaLess(theta) + ProfileNear(theta);
  }
  return out;
}

References CoulombAmplitude(), fProfileLambda, fRutherfordTheta, fWaveVector, GammaLess(), GammaMore(), PhaseNear(), CLHEP::pi, and ProfileNear().

Referenced by Amplitude().

AmplitudeSim()

G4complex G4NuclNuclDiffuseElastic::AmplitudeSim

(

G4double

theta

)

Definition at line 1633 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double sinThetaR  = 2.*fHalfRutThetaTg/(1. + fHalfRutThetaTg2);
  G4double dTheta     = 0.5*(theta - fRutherfordTheta);
  G4double sindTheta  = std::sin(dTheta);
  G4double persqrt2   = std::sqrt(0.5);
 
  G4complex order     = G4complex(persqrt2,persqrt2);
  order              *= std::sqrt(0.5*fProfileLambda/sinThetaR)*2.*sindTheta;
  // order              *= std::sqrt(0.5*fProfileLambda/sinThetaR)*2.*dTheta;
 
  G4complex out;
 
  if ( theta <= fRutherfordTheta )
  {
    out = 1. - 0.5*GetErfcInt(-order)*ProfileNear(theta);
  }
  else
  {
    out = 0.5*GetErfcInt(order)*ProfileNear(theta);
  }
 
  out *= CoulombAmplitude(theta);
 
  return out;
}

References CoulombAmplitude(), fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, GetErfcInt(), and ProfileNear().

Referenced by AmplitudeSimMod2().

AmplitudeSimMod2()

G4double G4NuclNuclDiffuseElastic::AmplitudeSimMod2 ( G4double  theta )

inline

Definition at line 1046 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex out = AmplitudeSim(theta);
  G4double mod2 = out.real()*out.real() + out.imag()*out.imag();
  return   mod2;
}

References AmplitudeSim().

ApplyYourself()

G4HadFinalState * G4HadronElastic::ApplyYourself ( const G4HadProjectile &  aTrack, G4Nucleus &  targetNucleus  )

overridevirtualinherited

Reimplemented from G4HadronicInteraction.

Reimplemented in G4NeutrinoElectronNcModel, G4NeutronElectronElModel, G4LEHadronProtonElastic, G4LEnp, and G4LEpp.

Definition at line 81 of file G4HadronElastic.cc.

{
  theParticleChange.Clear();
 
  const G4HadProjectile* aParticle = &aTrack;
  G4double ekin = aParticle->GetKineticEnergy();
 
  // no scattering below the limit
  if(ekin <= lowestEnergyLimit) {
    theParticleChange.SetEnergyChange(ekin);
    theParticleChange.SetMomentumChange(0.,0.,1.);
    return &theParticleChange;
  }
 
  G4int A = targetNucleus.GetA_asInt();
  G4int Z = targetNucleus.GetZ_asInt();
 
  // Scattered particle referred to axis of incident particle
  const G4ParticleDefinition* theParticle = aParticle->GetDefinition();
  G4double m1 = theParticle->GetPDGMass();
  G4double plab = std::sqrt(ekin*(ekin + 2.0*m1));
 
  if (verboseLevel>1) {
    G4cout << "G4HadronElastic: " 
       << aParticle->GetDefinition()->GetParticleName() 
       << " Plab(GeV/c)= " << plab/GeV  
       << " Ekin(MeV) = " << ekin/MeV 
       << " scattered off Z= " << Z 
       << " A= " << A 
       << G4endl;
  }
 
  G4double mass2 = G4NucleiProperties::GetNuclearMass(A, Z);
  G4double e1 = m1 + ekin;
  G4LorentzVector lv(0.0,0.0,plab,e1+mass2);
  G4ThreeVector bst = lv.boostVector();
  G4double momentumCMS = plab*mass2/std::sqrt(m1*m1 + mass2*mass2 + 2.*mass2*e1);
 
  pLocalTmax = 4.0*momentumCMS*momentumCMS;
 
  // Sampling in CM system
  G4double t = SampleInvariantT(theParticle, plab, Z, A);
 
  if(t < 0.0 || t > pLocalTmax) {
    // For the very rare cases where cos(theta) is greater than 1 or smaller than -1,
    // print some debugging information via a "JustWarning" exception, and resample
    // using the default algorithm
#ifdef G4VERBOSE
    if(nwarn < 2) {
      G4ExceptionDescription ed;
      ed << GetModelName() << " wrong sampling t= " << t << " tmax= " << pLocalTmax
     << " for " << aParticle->GetDefinition()->GetParticleName() 
     << " ekin=" << ekin << " MeV" 
     << " off (Z,A)=(" << Z << "," << A << ") - will be resampled" << G4endl;
      G4Exception( "G4HadronElastic::ApplyYourself", "hadEla001", JustWarning, ed);
      ++nwarn;
    }
#endif
    t = G4HadronElastic::SampleInvariantT(theParticle, plab, Z, A);
  }
 
  G4double phi  = G4UniformRand()*CLHEP::twopi;
  G4double cost = 1. - 2.0*t/pLocalTmax;
 
  if (cost > 1.0) { cost = 1.0; }
  else if(cost < -1.0) { cost = -1.0; } 
 
  G4double sint = std::sqrt((1.0-cost)*(1.0+cost));
 
  if (verboseLevel>1) {
    G4cout << " t= " << t << " tmax(GeV^2)= " << pLocalTmax/(GeV*GeV) 
       << " Pcms(GeV)= " << momentumCMS/GeV << " cos(t)=" << cost 
       << " sin(t)=" << sint << G4endl;
  }
  G4LorentzVector nlv1(momentumCMS*sint*std::cos(phi),
               momentumCMS*sint*std::sin(phi),
                       momentumCMS*cost,
               std::sqrt(momentumCMS*momentumCMS + m1*m1));
 
  nlv1.boost(bst); 
 
  G4double eFinal = nlv1.e() - m1;
  if (verboseLevel > 1) {
    G4cout <<"G4HadronElastic: m= " << m1 << " Efin(MeV)= " << eFinal 
       << " 4-M Final: " << nlv1 
       << G4endl;
  }
 
  if(eFinal <= 0.0) { 
    theParticleChange.SetMomentumChange(0.0,0.0,1.0);
    theParticleChange.SetEnergyChange(0.0);
  } else {
    theParticleChange.SetMomentumChange(nlv1.vect().unit());
    theParticleChange.SetEnergyChange(eFinal);
  }
  lv -= nlv1;
  G4double erec =  std::max(lv.e() - mass2, 0.0);
  if (verboseLevel > 1) {
    G4cout << "Recoil: " <<" m= " << mass2 << " Erec(MeV)= " << erec
       << " 4-mom: " << lv 
       << G4endl;
  }
 
  // the recoil is created if kinetic energy above the threshold
  if(erec > GetRecoilEnergyThreshold()) {
    G4ParticleDefinition * theDef = nullptr;
    if(Z == 1 && A == 1)       { theDef = theProton; }
    else if (Z == 1 && A == 2) { theDef = theDeuteron; }
    else if (Z == 1 && A == 3) { theDef = G4Triton::Triton(); }
    else if (Z == 2 && A == 3) { theDef = G4He3::He3(); }
    else if (Z == 2 && A == 4) { theDef = theAlpha; }
    else {
      theDef = 
    G4ParticleTable::GetParticleTable()->GetIonTable()->GetIon(Z,A,0.0);
    }
    G4DynamicParticle * aSec = new G4DynamicParticle(theDef, lv.vect().unit(), erec);
    theParticleChange.AddSecondary(aSec, secID);
  } else {
    theParticleChange.SetLocalEnergyDeposit(erec);
  }
 
  return &theParticleChange;
}

References A, G4HadFinalState::AddSecondary(), CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4HadFinalState::Clear(), CLHEP::HepLorentzVector::e(), e1, G4cout, G4endl, G4Exception(), G4UniformRand, G4Nucleus::GetA_asInt(), G4HadProjectile::GetDefinition(), G4IonTable::GetIon(), G4ParticleTable::GetIonTable(), G4HadProjectile::GetKineticEnergy(), G4HadronicInteraction::GetModelName(), G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetParticleName(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGMass(), G4HadronicInteraction::GetRecoilEnergyThreshold(), G4Nucleus::GetZ_asInt(), GeV, G4He3::He3(), JustWarning, G4HadronElastic::lowestEnergyLimit, G4INCL::Math::max(), MeV, G4HadronElastic::nwarn, G4HadronElastic::pLocalTmax, G4HadronElastic::SampleInvariantT(), G4HadronElastic::secID, G4HadFinalState::SetEnergyChange(), G4HadFinalState::SetLocalEnergyDeposit(), G4HadFinalState::SetMomentumChange(), G4HadronElastic::theAlpha, G4HadronElastic::theDeuteron, G4HadronicInteraction::theParticleChange, G4HadronElastic::theProton, G4Triton::Triton(), CLHEP::twopi, CLHEP::Hep3Vector::unit(), CLHEP::HepLorentzVector::vect(), G4HadronicInteraction::verboseLevel, and Z.

BesselJone()

G4double G4NuclNuclDiffuseElastic::BesselJone

(

G4double

z

)

Definition at line 2085 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double modvalue, value2, fact1, fact2, arg, shift, bessel;
 
  modvalue = std::fabs(value);
 
  if ( modvalue < 8.0 ) 
  {
    value2 = value*value;
 
    fact1  = value*(72362614232.0 + value2*(-7895059235.0 
                                  + value2*( 242396853.1
                                  + value2*(-2972611.439 
                                  + value2*( 15704.48260 
                                  + value2*(-30.16036606  ) ) ) ) ) );
 
    fact2  = 144725228442.0 + value2*(2300535178.0 
                            + value2*(18583304.74
                            + value2*(99447.43394 
                            + value2*(376.9991397 
                            + value2*1.0             ) ) ) );
    bessel = fact1/fact2;
  } 
  else 
  {
    arg    = 8.0/modvalue;
 
    value2 = arg*arg;
 
    shift  = modvalue - 2.356194491;
 
    fact1  = 1.0 + value2*( 0.183105e-2 
                 + value2*(-0.3516396496e-4
                 + value2*(0.2457520174e-5 
                 + value2*(-0.240337019e-6          ) ) ) );
 
    fact2 = 0.04687499995 + value2*(-0.2002690873e-3
                          + value2*( 0.8449199096e-5
                          + value2*(-0.88228987e-6
                          + value2*0.105787412e-6       ) ) );
 
    bessel = std::sqrt( 0.636619772/modvalue)*(std::cos(shift)*fact1 - arg*std::sin(shift)*fact2);
 
    if (value < 0.0) bessel = -bessel;
  }
  return bessel;
}

Referenced by BesselOneByArg(), GetDiffElasticProb(), GetDiffElasticSumProb(), and GetDiffElasticSumProbA().

BesselJzero()

G4double G4NuclNuclDiffuseElastic::BesselJzero

(

G4double

z

)

Definition at line 2033 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double modvalue, value2, fact1, fact2, arg, shift, bessel;
 
  modvalue = std::fabs(value);
 
  if ( value < 8.0 && value > -8.0 )
  {
    value2 = value*value;
 
    fact1  = 57568490574.0 + value2*(-13362590354.0 
                           + value2*( 651619640.7 
                           + value2*(-11214424.18 
                           + value2*( 77392.33017 
                           + value2*(-184.9052456   ) ) ) ) );
 
    fact2  = 57568490411.0 + value2*( 1029532985.0 
                           + value2*( 9494680.718
                           + value2*(59272.64853
                           + value2*(267.8532712 
                           + value2*1.0               ) ) ) );
 
    bessel = fact1/fact2;
  } 
  else 
  {
    arg    = 8.0/modvalue;
 
    value2 = arg*arg;
 
    shift  = modvalue-0.785398164;
 
    fact1  = 1.0 + value2*(-0.1098628627e-2 
                 + value2*(0.2734510407e-4
                 + value2*(-0.2073370639e-5 
                 + value2*0.2093887211e-6    ) ) );
 
    fact2  = -0.1562499995e-1 + value2*(0.1430488765e-3
                              + value2*(-0.6911147651e-5 
                              + value2*(0.7621095161e-6
                              - value2*0.934945152e-7    ) ) );
 
    bessel = std::sqrt(0.636619772/modvalue)*(std::cos(shift)*fact1 - arg*std::sin(shift)*fact2 );
  }
  return bessel;
}

Referenced by GetDiffElasticProb(), GetDiffElasticSumProb(), and GetDiffElasticSumProbA().

BesselOneByArg()

G4double G4NuclNuclDiffuseElastic::BesselOneByArg ( G4double  z )

inline

Definition at line 419 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double x2, result;
  
  if( std::fabs(x) < 0.01 )
  { 
   x     *= 0.5;
   x2     = x*x;
   result = 2. - x2 + x2*x2/6.;
  }
  else
  {
    result = BesselJone(x)/x; 
  }
  return result;
}

References BesselJone().

Referenced by GetDiffElasticProb(), GetDiffElasticSumProb(), and GetDiffElasticSumProbA().

Block()

void G4HadronicInteraction::Block ( )

inlineprotectedinherited

Definition at line 170 of file G4HadronicInteraction.hh.

{ isBlocked = true; }

References G4HadronicInteraction::isBlocked.

Referenced by G4HadronicInteraction::ActivateFor(), G4HadronicInteraction::DeActivateFor(), G4HadronicInteraction::SetMaxEnergy(), and G4HadronicInteraction::SetMinEnergy().

BuildAngleTable()

void G4NuclNuclDiffuseElastic::BuildAngleTable

(

)

Definition at line 1007 of file G4NuclNuclDiffuseElastic.cc.

{
  G4int i, j;
  G4double partMom, kinE, m1 = fParticle->GetPDGMass();
  G4double alpha1, alpha2, alphaMax, alphaCoulomb, delta = 0., sum = 0.;
 
  // G4cout<<"particle z = "<<z<<"; particle m1 = "<<m1/GeV<<" GeV"<<G4endl;
 
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
  
  fAngleTable = new G4PhysicsTable(fEnergyBin);
 
  for( i = 0; i < fEnergyBin; i++)
  {
    kinE        = fEnergyVector->GetLowEdgeEnergy(i);
 
    // G4cout<<G4endl;
    // G4cout<<"kinE = "<<kinE/MeV<<" MeV"<<G4endl;
 
    partMom     = std::sqrt( kinE*(kinE + 2*m1) );
 
    InitDynParameters(fParticle, partMom);
 
    alphaMax = fRutherfordTheta*fCofAlphaMax;
 
    if(alphaMax > pi) alphaMax = pi;
 
    // VI: Coverity complain
    //alphaMax = pi2;
 
    alphaCoulomb = fRutherfordTheta*fCofAlphaCoulomb;
 
    // G4cout<<"alphaCoulomb = "<<alphaCoulomb/degree<<"; alphaMax = "<<alphaMax/degree<<G4endl;
 
    G4PhysicsFreeVector* angleVector = new G4PhysicsFreeVector(fAngleBin-1);
 
    // G4PhysicsLogVector*  angleBins = new G4PhysicsLogVector( 0.001*alphaMax, alphaMax, fAngleBin );
 
    G4double delth = (alphaMax-alphaCoulomb)/fAngleBin;
        
    sum = 0.;
 
    // fAddCoulomb = false;
    fAddCoulomb = true;
 
    // for(j = 1; j < fAngleBin; j++)
    for(j = fAngleBin-1; j >= 1; j--)
    {
      // alpha1 = angleBins->GetLowEdgeEnergy(j-1);
      // alpha2 = angleBins->GetLowEdgeEnergy(j);
 
      // alpha1 = alphaMax*(j-1)/fAngleBin;
      // alpha2 = alphaMax*( j )/fAngleBin;
 
      alpha1 = alphaCoulomb + delth*(j-1);
      // if(alpha1 < kRlim2) alpha1 = kRlim2;
      alpha2 = alpha1 + delth;
 
      delta = integral.Legendre10(this, &G4NuclNuclDiffuseElastic::GetFresnelIntegrandXsc, alpha1, alpha2);
      // delta = integral.Legendre96(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, alpha1, alpha2);
 
      sum += delta;
      
      angleVector->PutValue( j-1 , alpha1, sum ); // alpha2
      // G4cout<<"j-1 = "<<j-1<<"; alpha2 = "<<alpha2/degree<<"; sum = "<<sum<<G4endl;
    }
    fAngleTable->insertAt(i,angleVector);
 
    // delete[] angleVector; 
    // delete[] angleBins; 
  }
  return;
}

References alpha2, fAddCoulomb, fAngleBin, fAngleTable, fCofAlphaCoulomb, fCofAlphaMax, fEnergyBin, fEnergyVector, fParticle, fRutherfordTheta, G4NuclNuclDiffuseElastic(), GetFresnelIntegrandXsc(), G4PhysicsVector::GetLowEdgeEnergy(), G4ParticleDefinition::GetPDGMass(), InitDynParameters(), G4PhysicsTable::insertAt(), pi, and G4PhysicsFreeVector::PutValue().

Referenced by Initialise(), and InitialiseOnFly().

BuildPhysicsTable()

void G4HadronicInteraction::BuildPhysicsTable ( const G4ParticleDefinition &  )

virtualinherited

Reimplemented in G4LENDorBERTModel, G4LENDCombinedModel, G4LENDGammaModel, G4LENDModel, G4ParticleHPCapture, G4ParticleHPElastic, G4ParticleHPFission, G4ParticleHPInelastic, G4ParticleHPThermalScattering, G4AblaInterface, and G4PreCompoundModel.

Definition at line 56 of file G4HadronicInteraction.cc.

{}

CalcMandelstamS()

G4double G4NuclNuclDiffuseElastic::CalcMandelstamS ( const G4double  mp, const G4double  mt, const G4double  Plab  )

inline

Definition at line 1080 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double Elab = std::sqrt ( mp * mp + Plab * Plab );
  G4double sMand  = mp*mp + mt*mt + 2*Elab*mt ;
 
  return sMand;
}

Referenced by GetHadronNucleonXscNS().

CalculateAm()

G4double G4NuclNuclDiffuseElastic::CalculateAm ( G4double  momentum, G4double  n, G4double  Z  )

inline

Definition at line 468 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double k   = momentum/CLHEP::hbarc;
  G4double ch  = 1.13 + 3.76*n*n;
  G4double zn  = 1.77*k*(1.0/G4Pow::GetInstance()->A13(Z))*CLHEP::Bohr_radius;
  G4double zn2 = zn*zn;
  fAm          = ch/zn2;
 
  return fAm;
}

References G4Pow::A13(), CLHEP::Bohr_radius, fAm, G4Pow::GetInstance(), CLHEP::hbarc, CLHEP::detail::n, and Z.

Referenced by GetCoulombElasticXsc(), GetCoulombIntegralXsc(), GetCoulombTotalXsc(), GetDiffuseElasticSumXsc(), InitDynParameters(), InitParameters(), InitParametersGla(), and TestAngleTable().

CalculateCoulombPhase()

G4double G4NuclNuclDiffuseElastic::CalculateCoulombPhase ( G4int  n )

inline

Definition at line 844 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex z        = G4complex(1. + n, fZommerfeld); 
  // G4complex gammalog = GammaLogarithm(z);
  G4complex gammalog = GammaLogB2n(z);
  return gammalog.imag();
}

References fZommerfeld, GammaLogB2n(), and CLHEP::detail::n.

Referenced by AmplitudeGla().

CalculateCoulombPhaseZero()

void G4NuclNuclDiffuseElastic::CalculateCoulombPhaseZero ( )

inline

Definition at line 831 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex z        = G4complex(1,fZommerfeld); 
  // G4complex gammalog = GammaLogarithm(z);
  G4complex gammalog = GammaLogB2n(z);
  fCoulombPhase0     = gammalog.imag();
}

References fCoulombPhase0, fZommerfeld, and GammaLogB2n().

Referenced by InitDynParameters(), InitParameters(), and InitParametersGla().

CalculateNuclearRad()

G4double G4NuclNuclDiffuseElastic::CalculateNuclearRad ( G4double  A )

inline

Definition at line 483 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double r0 = 1.*CLHEP::fermi, radius;
  // r0 *= 1.12;
  // r0 *= 1.44;
  r0 *= fNuclearRadiusCof;
 
  /*
  if( A < 50. )
  {
    if( A > 10. ) r0  = 1.16*( 1 - (1.0/G4Pow::GetInstance()->A23(A)) )*CLHEP::fermi;   // 1.08*fermi;
    else          r0  = 1.1*CLHEP::fermi;
 
    radius = r0*G4Pow::GetInstance()->A13(A);
  }
  else
  {
    r0 = 1.7*CLHEP::fermi;   // 1.7*fermi;
 
    radius = r0*G4Pow::GetInstance()->powA(A, 0.27); // 0.27);
  }
  */
  radius = r0*G4Pow::GetInstance()->A13(A);
 
  return radius;
}

References A, G4Pow::A13(), CLHEP::fermi, fNuclearRadiusCof, and G4Pow::GetInstance().

Referenced by GetDiffuseElasticSumXsc(), GetDiffuseElasticXsc(), Initialise(), InitialiseOnFly(), InitParameters(), InitParametersGla(), IntegralElasticProb(), SampleCoulombMuCMS(), SampleThetaCMS(), and TestAngleTable().

CalculateParticleBeta()

G4double G4NuclNuclDiffuseElastic::CalculateParticleBeta ( const G4ParticleDefinition *  particle, G4double  momentum  )

inline

Definition at line 441 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double mass = particle->GetPDGMass();
  G4double a    = momentum/mass;
  fBeta         = a/std::sqrt(1+a*a);
 
  return fBeta; 
}

References fBeta, and G4ParticleDefinition::GetPDGMass().

Referenced by GetCoulombElasticXsc(), GetCoulombIntegralXsc(), GetCoulombTotalXsc(), and GetDiffuseElasticSumXsc().

CalculateRutherfordAnglePar()

void G4NuclNuclDiffuseElastic::CalculateRutherfordAnglePar ( )

inline

Definition at line 858 of file G4NuclNuclDiffuseElastic.hh.

{
  fHalfRutThetaTg   = fZommerfeld/fProfileLambda;  // (fWaveVector*fNuclearRadius);
  fRutherfordTheta  = 2.*std::atan(fHalfRutThetaTg);
  fHalfRutThetaTg2  = fHalfRutThetaTg*fHalfRutThetaTg;
  // G4cout<<"fRutherfordTheta = "<<fRutherfordTheta/degree<<" degree"<<G4endl;
 
}

References fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, and fZommerfeld.

Referenced by InitDynParameters(), and InitParameters().

CalculateZommerfeld()

G4double G4NuclNuclDiffuseElastic::CalculateZommerfeld ( G4double  beta, G4double  Z1, G4double  Z2  )

inline

Definition at line 456 of file G4NuclNuclDiffuseElastic.hh.

{
  fZommerfeld = CLHEP::fine_structure_const*Z1*Z2/beta;
 
  return fZommerfeld; 
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CLHEP::fine_structure_const, fZommerfeld, and anonymous_namespace{paraMaker.cc}::Z1.

Referenced by GetCoulombElasticXsc(), GetCoulombIntegralXsc(), GetCoulombTotalXsc(), GetDiffuseElasticSumXsc(), InitDynParameters(), InitParameters(), InitParametersGla(), and TestAngleTable().

ComputeMomentumCMS()

G4double G4HadronElastic::ComputeMomentumCMS ( const G4ParticleDefinition *  p, G4double  plab, G4int  Z, G4int  A  )

inlineinherited

Definition at line 102 of file G4HadronElastic.hh.

{
  G4double m1 = p->GetPDGMass();
  G4double m12= m1*m1;
  G4double mass2 = G4NucleiProperties::GetNuclearMass(A, Z);
  return plab*mass2/std::sqrt(m12 + mass2*mass2 + 2.*mass2*std::sqrt(m12 + plab*plab));
}

References A, G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGMass(), and Z.

CoulombAmplitude()

G4complex G4NuclNuclDiffuseElastic::CoulombAmplitude ( G4double  theta )

inline

Definition at line 797 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex ca;
  
  G4double sinHalfTheta  = std::sin(0.5*theta);
  G4double sinHalfTheta2 = sinHalfTheta*sinHalfTheta; 
  sinHalfTheta2         += fAm;
 
  G4double order         = 2.*fCoulombPhase0 - fZommerfeld*G4Log(sinHalfTheta2);
  G4complex z            = G4complex(0., order);
  ca                     = std::exp(z);
 
  ca                    *= -fZommerfeld/(2.*fWaveVector*sinHalfTheta2);
 
  return ca; 
}

References fAm, fCoulombPhase0, fWaveVector, fZommerfeld, and G4Log().

Referenced by AmplitudeGG(), AmplitudeGla(), AmplitudeNear(), AmplitudeSim(), and CoulombAmplitudeMod2().

CoulombAmplitudeMod2()

G4double G4NuclNuclDiffuseElastic::CoulombAmplitudeMod2 ( G4double  theta )

inline

Definition at line 818 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex ca = CoulombAmplitude(theta);
  G4double out = ca.real()*ca.real() + ca.imag()*ca.imag();
 
  return  out;
}

References CoulombAmplitude().

DampFactor()

G4double G4NuclNuclDiffuseElastic::DampFactor ( G4double  z )

inline

Definition at line 396 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double df;
  G4double f2 = 2., f3 = 6., f4 = 24.; // first factorials
 
  // x *= pi;
 
  if( std::fabs(x) < 0.01 )
  { 
    df = 1./(1. + x/f2 + x*x/f3 + x*x*x/f4);
  }
  else
  {
    df = x/std::sinh(x); 
  }
  return df;
}

Referenced by GetDiffElasticProb(), GetDiffElasticSumProb(), and GetDiffElasticSumProbA().

DeActivateFor() [1/2]

void G4HadronicInteraction::DeActivateFor ( const G4Element *  anElement )

inherited

Definition at line 186 of file G4HadronicInteraction.cc.

{
  Block(); 
  theBlockedListElements.push_back(anElement);
}

References G4HadronicInteraction::Block(), and G4HadronicInteraction::theBlockedListElements.

DeActivateFor() [2/2]

void G4HadronicInteraction::DeActivateFor ( const G4Material *  aMaterial )

inherited

Definition at line 180 of file G4HadronicInteraction.cc.

{
  Block(); 
  theBlockedList.push_back(aMaterial);
}

References G4HadronicInteraction::Block(), and G4HadronicInteraction::theBlockedList.

Referenced by G4HadronHElasticPhysics::ConstructProcess().

GammaLess()

G4complex G4NuclNuclDiffuseElastic::GammaLess

(

G4double

theta

)

Definition at line 1553 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double sinThetaR      = 2.*fHalfRutThetaTg/(1. + fHalfRutThetaTg2);
  G4double cosHalfThetaR2 = 1./(1. + fHalfRutThetaTg2);
 
  G4double u              = std::sqrt(0.5*fProfileLambda/sinThetaR);
  G4double kappa          = u/std::sqrt(CLHEP::pi);
  G4double dTheta         = theta - fRutherfordTheta;
  u                      *= dTheta;
  G4double u2             = u*u;
  G4double u2m2p3         = u2*2./3.;
 
  G4complex im            = G4complex(0.,1.);
  G4complex order         = G4complex(u,u);
  order                  /= std::sqrt(2.);
 
  G4complex gamma         = CLHEP::pi*kappa*GetErfcInt(-order)*std::exp(im*(u*u+0.25*CLHEP::pi));
  G4complex a0            = 0.5*(1. + 4.*(1.+im*u2)*cosHalfThetaR2/3.)/sinThetaR;
  G4complex a1            = 0.5*(1. + 2.*(1.+im*u2m2p3)*cosHalfThetaR2)/sinThetaR;
  G4complex out           = gamma*(1. - a1*dTheta) - a0;
 
  return out;
}

References a0, fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, GetErfcInt(), and CLHEP::pi.

Referenced by AmplitudeNear().

GammaLogarithm()

G4complex G4NuclNuclDiffuseElastic::GammaLogarithm

(

G4complex

xx

)

Definition at line 2009 of file G4NuclNuclDiffuseElastic.cc.

{
  const G4double cof[6] = { 76.18009172947146,     -86.50532032941677,
                             24.01409824083091,      -1.231739572450155,
                              0.1208650973866179e-2, -0.5395239384953e-5  } ;
  G4int j;
  G4complex z = zz - 1.0;
  G4complex tmp = z + 5.5;
  tmp -= (z + 0.5) * std::log(tmp);
  G4complex ser = G4complex(1.000000000190015,0.);
 
  for ( j = 0; j <= 5; j++ )
  {
    z += 1.0;
    ser += cof[j]/z;
  }
  return -tmp + std::log(2.5066282746310005*ser);
}

GammaLogB2n()

G4complex G4NuclNuclDiffuseElastic::GammaLogB2n ( G4complex  xx )

inline

Definition at line 610 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex z1 = 12.*z;
  G4complex z2 = z*z;
  G4complex z3 = z2*z;
  G4complex z5 = z2*z3;
  G4complex z7 = z2*z5;
 
  z3 *= 360.;
  z5 *= 1260.;
  z7 *= 1680.;
 
  G4complex result  = (z-0.5)*std::log(z)-z+0.5*G4Log(CLHEP::twopi);
            result += 1./z1 - 1./z3 +1./z5 -1./z7;
  return result;
}

References G4Log(), and CLHEP::twopi.

Referenced by CalculateCoulombPhase(), and CalculateCoulombPhaseZero().

GammaMore()

G4complex G4NuclNuclDiffuseElastic::GammaMore

(

G4double

theta

)

Definition at line 1581 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double sinThetaR      = 2.*fHalfRutThetaTg/(1. + fHalfRutThetaTg2);
  G4double cosHalfThetaR2 = 1./(1. + fHalfRutThetaTg2);
 
  G4double u              = std::sqrt(0.5*fProfileLambda/sinThetaR);
  G4double kappa          = u/std::sqrt(CLHEP::pi);
  G4double dTheta         = theta - fRutherfordTheta;
  u                      *= dTheta;
  G4double u2             = u*u;
  G4double u2m2p3         = u2*2./3.;
 
  G4complex im            = G4complex(0.,1.);
  G4complex order         = G4complex(u,u);
  order                  /= std::sqrt(2.);
  G4complex gamma         = CLHEP::pi*kappa*GetErfcInt(order)*std::exp(im*(u*u+0.25*CLHEP::pi));
  G4complex a0            = 0.5*(1. + 4.*(1.+im*u2)*cosHalfThetaR2/3.)/sinThetaR;
  G4complex a1            = 0.5*(1. + 2.*(1.+im*u2m2p3)*cosHalfThetaR2)/sinThetaR;
  G4complex out           = -gamma*(1. - a1*dTheta) - a0;
 
  return out;
}

References a0, fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, GetErfcInt(), and CLHEP::pi.

Referenced by AmplitudeNear().

GetCint()

G4double G4NuclNuclDiffuseElastic::GetCint ( G4double  x )

inline

Definition at line 767 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double out;
 
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
 
  out= integral.Legendre96(this,&G4NuclNuclDiffuseElastic::GetCosHaPit2, 0., x );
 
  return out;
}

References G4NuclNuclDiffuseElastic(), and GetCosHaPit2().

Referenced by GetRatioGen(), and GetRatioSim().

GetCofAlphaCoulomb()

G4double G4NuclNuclDiffuseElastic::GetCofAlphaCoulomb ( )

inline

Definition at line 293 of file G4NuclNuclDiffuseElastic.hh.

{return fCofAlphaCoulomb;};

References fCofAlphaCoulomb.

GetCofAlphaMax()

G4double G4NuclNuclDiffuseElastic::GetCofAlphaMax ( )

inline

Definition at line 292 of file G4NuclNuclDiffuseElastic.hh.

{return fCofAlphaMax;};

References fCofAlphaMax.

GetCosHaPit2()

G4double G4NuclNuclDiffuseElastic::GetCosHaPit2 ( G4double  t )

inline

Definition at line 195 of file G4NuclNuclDiffuseElastic.hh.

{return std::cos(CLHEP::halfpi*t*t);};

References CLHEP::halfpi.

Referenced by GetCint().

GetCoulombElasticXsc()

G4double G4NuclNuclDiffuseElastic::GetCoulombElasticXsc ( const G4ParticleDefinition *  particle, G4double  theta, G4double  momentum, G4double  Z  )

inline

Definition at line 514 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double sinHalfTheta  = std::sin(0.5*theta);
  G4double sinHalfTheta2 = sinHalfTheta*sinHalfTheta;
  G4double beta          = CalculateParticleBeta( particle, momentum);
  G4double z             = particle->GetPDGCharge();
  G4double n             = CalculateZommerfeld( beta, z, Z );
  G4double am            = CalculateAm( momentum, n, Z);
  G4double k             = momentum/CLHEP::hbarc;
  G4double ch            = 0.5*n/k;
  G4double ch2           = ch*ch;
  G4double xsc           = ch2/((sinHalfTheta2+am)*(sinHalfTheta2+am));
 
  return xsc;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CalculateAm(), CalculateParticleBeta(), CalculateZommerfeld(), G4ParticleDefinition::GetPDGCharge(), CLHEP::hbarc, CLHEP::detail::n, and Z.

Referenced by GetInvCoulombElasticXsc().

GetCoulombIntegralXsc()

G4double G4NuclNuclDiffuseElastic::GetCoulombIntegralXsc ( const G4ParticleDefinition *  particle, G4double  momentum, G4double  Z, G4double  theta1, G4double  theta2  )

inline

Definition at line 579 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double c1 = std::cos(theta1);
  //G4cout<<"c1 = "<<c1<<G4endl;
  G4double c2 = std::cos(theta2);
  // G4cout<<"c2 = "<<c2<<G4endl;
  G4double beta          = CalculateParticleBeta( particle, momentum);
  // G4cout<<"beta = "<<beta<<G4endl;
  G4double z             = particle->GetPDGCharge();
  G4double n             = CalculateZommerfeld( beta, z, Z );
  // G4cout<<"fZomerfeld = "<<n<<G4endl;
  G4double am            = CalculateAm( momentum, n, Z);
  // G4cout<<"cof Am = "<<am<<G4endl;
  G4double k             = momentum/CLHEP::hbarc;
  // G4cout<<"k = "<<k*CLHEP::fermi<<" 1/fermi"<<G4endl;
  // G4cout<<"k*Bohr_radius = "<<k*CLHEP::Bohr_radius<<G4endl;
  G4double ch            = n/k;
  G4double ch2           = ch*ch;
  am *= 2.;
  G4double xsc = ch2*CLHEP::twopi*(c1-c2)/((1 - c1 + am)*(1 - c2 + am));
 
  return xsc;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CalculateAm(), CalculateParticleBeta(), CalculateZommerfeld(), G4ParticleDefinition::GetPDGCharge(), CLHEP::hbarc, CLHEP::detail::n, CLHEP::twopi, and Z.

GetCoulombTotalXsc()

G4double G4NuclNuclDiffuseElastic::GetCoulombTotalXsc ( const G4ParticleDefinition *  particle, G4double  momentum, G4double  Z  )

inline

Definition at line 553 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double beta          = CalculateParticleBeta( particle, momentum);
  G4cout<<"beta = "<<beta<<G4endl;
  G4double z             = particle->GetPDGCharge();
  G4double n             = CalculateZommerfeld( beta, z, Z );
  G4cout<<"fZomerfeld = "<<n<<G4endl;
  G4double am            = CalculateAm( momentum, n, Z);
  G4cout<<"cof Am = "<<am<<G4endl;
  G4double k             = momentum/CLHEP::hbarc;
  G4cout<<"k = "<<k*CLHEP::fermi<<" 1/fermi"<<G4endl;
  G4cout<<"k*Bohr_radius = "<<k*CLHEP::Bohr_radius<<G4endl;
  G4double ch            = n/k;
  G4double ch2           = ch*ch;
  G4double xsc           = ch2*CLHEP::pi/(am +am*am);
 
  return xsc;
}

References anonymous_namespace{G4PionRadiativeDecayChannel.cc}::beta, CLHEP::Bohr_radius, CalculateAm(), CalculateParticleBeta(), CalculateZommerfeld(), CLHEP::fermi, G4cout, G4endl, G4ParticleDefinition::GetPDGCharge(), CLHEP::hbarc, CLHEP::detail::n, CLHEP::pi, and Z.

GetDiffElasticProb()

G4double G4NuclNuclDiffuseElastic::GetDiffElasticProb

(

G4double

theta

)

Definition at line 397 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double sigma, bzero, bzero2, bonebyarg, bonebyarg2, damp, damp2;
  G4double delta, diffuse, gamma;
  G4double e1, e2, bone, bone2;
 
  // G4double wavek = momentum/hbarc;  // wave vector
  // G4double r0    = 1.08*fermi;
  // G4double rad   = r0*G4Pow::GetInstance()->A13(A);
 
  G4double kr    = fWaveVector*fNuclearRadius; // wavek*rad;
  G4double kr2   = kr*kr;
  G4double krt   = kr*theta;
 
  bzero      = BesselJzero(krt);
  bzero2     = bzero*bzero;    
  bone       = BesselJone(krt);
  bone2      = bone*bone;
  bonebyarg  = BesselOneByArg(krt);
  bonebyarg2 = bonebyarg*bonebyarg;  
 
  // VI - Coverity complains
  /*
  if (fParticle == theProton)
  {
    diffuse = 0.63*fermi;
    gamma   = 0.3*fermi;
    delta   = 0.1*fermi*fermi;
    e1      = 0.3*fermi;
    e2      = 0.35*fermi;
  }
  else // as proton, if were not defined 
  {
  */
    diffuse = 0.63*fermi;
    gamma   = 0.3*fermi;
    delta   = 0.1*fermi*fermi;
    e1      = 0.3*fermi;
    e2      = 0.35*fermi;
    //}
  G4double lambda = 15.; // 15 ok
 
  //  G4double kgamma    = fWaveVector*gamma;   // wavek*delta;
 
  G4double kgamma    = lambda*(1.-G4Exp(-fWaveVector*gamma/lambda));   // wavek*delta;
  G4double kgamma2   = kgamma*kgamma;
 
  // G4double dk2t  = delta*fWaveVector*fWaveVector*theta; // delta*wavek*wavek*theta;
  // G4double dk2t2 = dk2t*dk2t;
  // G4double pikdt = pi*fWaveVector*diffuse*theta;// pi*wavek*diffuse*theta;
 
  G4double pikdt    = lambda*(1.-G4Exp(-pi*fWaveVector*diffuse*theta/lambda));   // wavek*delta;
 
  damp           = DampFactor(pikdt);
  damp2          = damp*damp;
 
  G4double mode2k2 = (e1*e1+e2*e2)*fWaveVector*fWaveVector;  
  G4double e2dk3t  = -2.*e2*delta*fWaveVector*fWaveVector*fWaveVector*theta;
 
 
  sigma  = kgamma2;
  // sigma  += dk2t2;
  sigma *= bzero2;
  sigma += mode2k2*bone2 + e2dk3t*bzero*bone;
  sigma += kr2*bonebyarg2;
  sigma *= damp2;          // *rad*rad;
 
  return sigma;
}

References BesselJone(), BesselJzero(), BesselOneByArg(), DampFactor(), e1, e2, fermi, fNuclearRadius, fWaveVector, G4Exp(), G4InuclParticleNames::lambda, and pi.

Referenced by GetDiffuseElasticXsc().

GetDiffElasticSumProb()

G4double G4NuclNuclDiffuseElastic::GetDiffElasticSumProb

(

G4double

theta

)

Definition at line 477 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double sigma, bzero, bzero2, bonebyarg, bonebyarg2, damp, damp2;
  G4double delta, diffuse, gamma;
  G4double e1, e2, bone, bone2;
 
  // G4double wavek = momentum/hbarc;  // wave vector
  // G4double r0    = 1.08*fermi;
  // G4double rad   = r0*G4Pow::GetInstance()->A13(A);
 
  G4double kr    = fWaveVector*fNuclearRadius; // wavek*rad;
  G4double kr2   = kr*kr;
  G4double krt   = kr*theta;
 
  bzero      = BesselJzero(krt);
  bzero2     = bzero*bzero;    
  bone       = BesselJone(krt);
  bone2      = bone*bone;
  bonebyarg  = BesselOneByArg(krt);
  bonebyarg2 = bonebyarg*bonebyarg;  
 
  if (fParticle == theProton)
  {
    diffuse = 0.63*fermi;
    // diffuse = 0.6*fermi;
    gamma   = 0.3*fermi;
    delta   = 0.1*fermi*fermi;
    e1      = 0.3*fermi;
    e2      = 0.35*fermi;
  }
  else // as proton, if were not defined 
  {
    diffuse = 0.63*fermi;
    gamma   = 0.3*fermi;
    delta   = 0.1*fermi*fermi;
    e1      = 0.3*fermi;
    e2      = 0.35*fermi;
  }
  G4double lambda = 15.; // 15 ok
  // G4double kgamma    = fWaveVector*gamma;   // wavek*delta;
  G4double kgamma    = lambda*(1.-G4Exp(-fWaveVector*gamma/lambda));   // wavek*delta;
 
  // G4cout<<"kgamma = "<<kgamma<<G4endl;
 
  if(fAddCoulomb)  // add Coulomb correction
  {
    G4double sinHalfTheta  = std::sin(0.5*theta);
    G4double sinHalfTheta2 = sinHalfTheta*sinHalfTheta;
 
    kgamma += 0.5*fZommerfeld/kr/(sinHalfTheta2+fAm); // correction at J0()
  // kgamma += 0.65*fZommerfeld/kr/(sinHalfTheta2+fAm); // correction at J0()
  }
 
  G4double kgamma2   = kgamma*kgamma;
 
  // G4double dk2t  = delta*fWaveVector*fWaveVector*theta; // delta*wavek*wavek*theta;
  //   G4cout<<"dk2t = "<<dk2t<<G4endl;
  // G4double dk2t2 = dk2t*dk2t;
  // G4double pikdt = pi*fWaveVector*diffuse*theta;// pi*wavek*diffuse*theta;
 
  G4double pikdt    = lambda*(1.-G4Exp(-pi*fWaveVector*diffuse*theta/lambda));   // wavek*delta;
 
  // G4cout<<"pikdt = "<<pikdt<<G4endl;
 
  damp           = DampFactor(pikdt);
  damp2          = damp*damp;
 
  G4double mode2k2 = (e1*e1+e2*e2)*fWaveVector*fWaveVector;  
  G4double e2dk3t  = -2.*e2*delta*fWaveVector*fWaveVector*fWaveVector*theta;
 
  sigma  = kgamma2;
  // sigma += dk2t2;
  sigma *= bzero2;
  sigma += mode2k2*bone2; 
  sigma += e2dk3t*bzero*bone;
 
  // sigma += kr2*(1 + 8.*fZommerfeld*fZommerfeld/kr2)*bonebyarg2;  // correction at J1()/()
  sigma += kr2*bonebyarg2;  // correction at J1()/()
 
  sigma *= damp2;          // *rad*rad;
 
  return sigma;
}

References BesselJone(), BesselJzero(), BesselOneByArg(), DampFactor(), e1, e2, fAddCoulomb, fAm, fermi, fNuclearRadius, fParticle, fWaveVector, fZommerfeld, G4Exp(), G4InuclParticleNames::lambda, pi, and theProton.

Referenced by GetDiffuseElasticSumXsc().

GetDiffElasticSumProbA()

G4double G4NuclNuclDiffuseElastic::GetDiffElasticSumProbA

(

G4double

alpha

)

Definition at line 572 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double theta; 
 
  theta = std::sqrt(alpha);
 
  // theta = std::acos( 1 - alpha/2. );
 
  G4double sigma, bzero, bzero2, bonebyarg, bonebyarg2, damp, damp2;
  G4double delta, diffuse, gamma;
  G4double e1, e2, bone, bone2;
 
  // G4double wavek = momentum/hbarc;  // wave vector
  // G4double r0    = 1.08*fermi;
  // G4double rad   = r0*G4Pow::GetInstance()->A13(A);
 
  G4double kr    = fWaveVector*fNuclearRadius; // wavek*rad;
  G4double kr2   = kr*kr;
  G4double krt   = kr*theta;
 
  bzero      = BesselJzero(krt);
  bzero2     = bzero*bzero;    
  bone       = BesselJone(krt);
  bone2      = bone*bone;
  bonebyarg  = BesselOneByArg(krt);
  bonebyarg2 = bonebyarg*bonebyarg;  
 
  if (fParticle == theProton)
  {
    diffuse = 0.63*fermi;
    // diffuse = 0.6*fermi;
    gamma   = 0.3*fermi;
    delta   = 0.1*fermi*fermi;
    e1      = 0.3*fermi;
    e2      = 0.35*fermi;
  }
  else // as proton, if were not defined 
  {
    diffuse = 0.63*fermi;
    gamma   = 0.3*fermi;
    delta   = 0.1*fermi*fermi;
    e1      = 0.3*fermi;
    e2      = 0.35*fermi;
  }
  G4double lambda = 15.; // 15 ok
  // G4double kgamma    = fWaveVector*gamma;   // wavek*delta;
  G4double kgamma    = lambda*(1.-G4Exp(-fWaveVector*gamma/lambda));   // wavek*delta;
 
  // G4cout<<"kgamma = "<<kgamma<<G4endl;
 
  if(fAddCoulomb)  // add Coulomb correction
  {
    G4double sinHalfTheta  = theta*0.5; // std::sin(0.5*theta);
    G4double sinHalfTheta2 = sinHalfTheta*sinHalfTheta;
 
    kgamma += 0.5*fZommerfeld/kr/(sinHalfTheta2+fAm); // correction at J0()
  // kgamma += 0.65*fZommerfeld/kr/(sinHalfTheta2+fAm); // correction at J0()
  }
 
  G4double kgamma2   = kgamma*kgamma;
 
  // G4double dk2t  = delta*fWaveVector*fWaveVector*theta; // delta*wavek*wavek*theta;
  //   G4cout<<"dk2t = "<<dk2t<<G4endl;
  // G4double dk2t2 = dk2t*dk2t;
  // G4double pikdt = pi*fWaveVector*diffuse*theta;// pi*wavek*diffuse*theta;
 
  G4double pikdt    = lambda*(1.-G4Exp(-pi*fWaveVector*diffuse*theta/lambda));   // wavek*delta;
 
  // G4cout<<"pikdt = "<<pikdt<<G4endl;
 
  damp           = DampFactor(pikdt);
  damp2          = damp*damp;
 
  G4double mode2k2 = (e1*e1+e2*e2)*fWaveVector*fWaveVector;  
  G4double e2dk3t  = -2.*e2*delta*fWaveVector*fWaveVector*fWaveVector*theta;
 
  sigma  = kgamma2;
  // sigma += dk2t2;
  sigma *= bzero2;
  sigma += mode2k2*bone2; 
  sigma += e2dk3t*bzero*bone;
 
  // sigma += kr2*(1 + 8.*fZommerfeld*fZommerfeld/kr2)*bonebyarg2;  // correction at J1()/()
  sigma += kr2*bonebyarg2;  // correction at J1()/()
 
  sigma *= damp2;          // *rad*rad;
 
  return sigma;
}

References alpha, BesselJone(), BesselJzero(), BesselOneByArg(), DampFactor(), e1, e2, fAddCoulomb, fAm, fermi, fNuclearRadius, fParticle, fWaveVector, fZommerfeld, G4Exp(), G4InuclParticleNames::lambda, pi, and theProton.

Referenced by GetIntegrandFunction().

GetDiffuseElasticSumXsc()

G4double G4NuclNuclDiffuseElastic::GetDiffuseElasticSumXsc	(	const G4ParticleDefinition *	particle,
		G4double	theta,
		G4double	momentum,
		G4double	A,
		G4double	Z
	)

Definition at line 261 of file G4NuclNuclDiffuseElastic.cc.

{
  fParticle      = particle;
  fWaveVector    = momentum/hbarc;
  fAtomicWeight  = A;
  fAtomicNumber  = Z;
  fNuclearRadius = CalculateNuclearRad(A);
  fAddCoulomb    = false;
 
  G4double z     = particle->GetPDGCharge();
 
  G4double kRt   = fWaveVector*fNuclearRadius*theta;
  G4double kRtC  = 1.9;
 
  if( z && (kRt > kRtC) )
  {
    fAddCoulomb = true;
    fBeta       = CalculateParticleBeta( particle, momentum);
    fZommerfeld = CalculateZommerfeld( fBeta, z, fAtomicNumber);
    fAm         = CalculateAm( momentum, fZommerfeld, fAtomicNumber);
  }
  G4double sigma = fNuclearRadius*fNuclearRadius*GetDiffElasticSumProb(theta);
 
  return sigma;
}

References A, CalculateAm(), CalculateNuclearRad(), CalculateParticleBeta(), CalculateZommerfeld(), fAddCoulomb, fAm, fAtomicNumber, fAtomicWeight, fBeta, fNuclearRadius, fParticle, fWaveVector, fZommerfeld, GetDiffElasticSumProb(), G4ParticleDefinition::GetPDGCharge(), source.hepunit::hbarc, and Z.

Referenced by GetInvElasticSumXsc().

GetDiffuseElasticXsc()

G4double G4NuclNuclDiffuseElastic::GetDiffuseElasticXsc	(	const G4ParticleDefinition *	particle,
		G4double	theta,
		G4double	momentum,
		G4double	A
	)

Definition at line 190 of file G4NuclNuclDiffuseElastic.cc.

{
  fParticle      = particle;
  fWaveVector    = momentum/hbarc;
  fAtomicWeight  = A;
  fAddCoulomb    = false;
  fNuclearRadius = CalculateNuclearRad(A);
 
  G4double sigma = fNuclearRadius*fNuclearRadius*GetDiffElasticProb(theta);
 
  return sigma;
}

References A, CalculateNuclearRad(), fAddCoulomb, fAtomicWeight, fNuclearRadius, fParticle, fWaveVector, GetDiffElasticProb(), and source.hepunit::hbarc.

Referenced by GetInvElasticXsc().

GetEnergyMomentumCheckLevels()

std::pair< G4double, G4double > G4HadronicInteraction::GetEnergyMomentumCheckLevels ( ) const

virtualinherited

Reimplemented in G4TheoFSGenerator.

Definition at line 217 of file G4HadronicInteraction.cc.

{
  return epCheckLevels;
}

References G4HadronicInteraction::epCheckLevels.

Referenced by G4HadronicProcess::CheckEnergyMomentumConservation(), and G4TheoFSGenerator::GetEnergyMomentumCheckLevels().

GetErf()

G4double G4NuclNuclDiffuseElastic::GetErf ( G4double  x )

inline

Definition at line 631 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double t, z, tmp, result;
 
  z   = std::fabs(x);
  t   = 1.0/(1.0+0.5*z);
 
  tmp = t*std::exp(-z*z-1.26551223+t*(1.00002368+t*(0.37409196+t*(0.09678418+
        t*(-0.18628806+t*(0.27886807+t*(-1.13520398+t*(1.48851587+
        t*(-0.82215223+t*0.17087277)))))))));
 
  if( x >= 0.) result = 1. - tmp;
  else         result = 1. + tmp; 
    
  return result;
}

Referenced by GetErfComp(), and GetErfInt().

GetErfcComp()

G4complex G4NuclNuclDiffuseElastic::GetErfcComp ( G4complex  z, G4int  nMax  )

inline

Definition at line 652 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex erfcz = 1. - GetErfComp( z, nMax);
  return erfcz;
}

References GetErfComp().

GetErfcInt()

G4complex G4NuclNuclDiffuseElastic::GetErfcInt ( G4complex  z )

inline

Definition at line 672 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex erfcz = 1. - GetErfInt( z); // , nMax);
  return erfcz;
}

References GetErfInt().

Referenced by AmplitudeSim(), GammaLess(), and GammaMore().

GetErfComp()

G4complex G4NuclNuclDiffuseElastic::GetErfComp	(	G4complex	z,
		G4int	nMax
	)

Definition at line 1466 of file G4NuclNuclDiffuseElastic.cc.

{
  G4int n;
  G4double n2, cofn, shny, chny, fn, gn;
 
  G4double x = z.real();
  G4double y = z.imag();
 
  G4double outRe = 0., outIm = 0.;
 
  G4double twox  = 2.*x;
  G4double twoxy = twox*y;
  G4double twox2 = twox*twox;
 
  G4double cof1 = G4Exp(-x*x)/CLHEP::pi;
 
  G4double cos2xy = std::cos(twoxy);
  G4double sin2xy = std::sin(twoxy);
 
  G4double twoxcos2xy = twox*cos2xy;
  G4double twoxsin2xy = twox*sin2xy;
 
  for( n = 1; n <= nMax; n++)
  {
    n2   = n*n;
 
    cofn = G4Exp(-0.5*n2)/(n2+twox2);  // /(n2+0.5*twox2);
 
    chny = std::cosh(n*y);
    shny = std::sinh(n*y);
 
    fn  = twox - twoxcos2xy*chny + n*sin2xy*shny;
    gn  =        twoxsin2xy*chny + n*cos2xy*shny;
 
    fn *= cofn;
    gn *= cofn;
 
    outRe += fn;
    outIm += gn;
  }
  outRe *= 2*cof1;
  outIm *= 2*cof1;
 
  if(std::abs(x) < 0.0001)
  {
    outRe += GetErf(x);
    outIm += cof1*y;
  }
  else
  {
    outRe += GetErf(x) + cof1*(1-cos2xy)/twox;
    outIm += cof1*sin2xy/twox;
  }
  return G4complex(outRe, outIm);
}

References G4Exp(), GetErf(), CLHEP::detail::n, and CLHEP::pi.

Referenced by GetErfcComp(), and TestErfcComp().

GetErfcSer()

G4complex G4NuclNuclDiffuseElastic::GetErfcSer ( G4complex  z, G4int  nMax  )

inline

Definition at line 662 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex erfcz = 1. - GetErfSer( z, nMax);
  return erfcz;
}

References GetErfSer().

GetErfInt()

G4complex G4NuclNuclDiffuseElastic::GetErfInt

(

G4complex

z

)

Definition at line 1527 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double outRe, outIm;
 
  G4double x = z.real();
  G4double y = z.imag();
  fReZ       = x;
 
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
 
  outRe = integral.Legendre96(this,&G4NuclNuclDiffuseElastic::GetExpSin, 0., y );
  outIm = integral.Legendre96(this,&G4NuclNuclDiffuseElastic::GetExpCos, 0., y );
 
  outRe *= 2./std::sqrt(CLHEP::pi);
  outIm *= 2./std::sqrt(CLHEP::pi);
 
  outRe += GetErf(x);
 
  return G4complex(outRe, outIm);
}

References fReZ, G4NuclNuclDiffuseElastic(), GetErf(), GetExpCos(), GetExpSin(), and CLHEP::pi.

Referenced by GetErfcInt(), and TestErfcInt().

GetErfSer()

G4complex G4NuclNuclDiffuseElastic::GetErfSer ( G4complex  z, G4int  nMax  )

inline

Definition at line 719 of file G4NuclNuclDiffuseElastic.hh.

{
  G4int n;
  G4double a =1., b = 1., tmp;
  G4complex sum = z, d = z;
 
  for( n = 1; n <= nMax; n++)
  {
    a *= 2.;
    b *= 2.*n +1.;
    d *= z*z;
 
    tmp = a/b;
 
    sum += tmp*d;
  }
  sum *= 2.*std::exp(-z*z)/std::sqrt(CLHEP::pi);
 
  return sum;
}

References CLHEP::detail::n, and CLHEP::pi.

Referenced by GetErfcSer(), and TestErfcSer().

GetExpCos()

G4double G4NuclNuclDiffuseElastic::GetExpCos ( G4double  x )

inline

Definition at line 742 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double result;
 
  result = G4Exp(x*x-fReZ*fReZ);
  result *= std::cos(2.*x*fReZ);
  return result;
}

References fReZ, and G4Exp().

Referenced by GetErfInt().

GetExpSin()

G4double G4NuclNuclDiffuseElastic::GetExpSin ( G4double  x )

inline

Definition at line 753 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double result;
 
  result = G4Exp(x*x-fReZ*fReZ);
  result *= std::sin(2.*x*fReZ);
  return result;
}

References fReZ, and G4Exp().

Referenced by GetErfInt().

GetFatalEnergyCheckLevels()

const std::pair< G4double, G4double > G4HadronicInteraction::GetFatalEnergyCheckLevels ( ) const

virtualinherited

Reimplemented in G4FissLib, G4LFission, G4LENDFission, G4ParticleHPCapture, G4ParticleHPElastic, G4ParticleHPFission, G4ParticleHPInelastic, and G4ParticleHPThermalScattering.

Definition at line 210 of file G4HadronicInteraction.cc.

{
  // default level of Check
  return std::pair<G4double, G4double>(2.*perCent, 1. * GeV);
}

References GeV, and perCent.

Referenced by G4HadronicProcess::CheckResult().

GetFresnelDiffuseXsc()

G4double G4NuclNuclDiffuseElastic::GetFresnelDiffuseXsc ( G4double  theta )

inline

Definition at line 1023 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double ratio   = GetRatioGen(theta);
  G4double ruthXsc = GetRutherfordXsc(theta);
  G4double xsc     = ratio*ruthXsc;
  return xsc;
}

References GetRatioGen(), and GetRutherfordXsc().

Referenced by GetFresnelIntegrandXsc().

GetFresnelIntegrandXsc()

G4double G4NuclNuclDiffuseElastic::GetFresnelIntegrandXsc ( G4double  alpha )

inline

Definition at line 1035 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double theta = std::sqrt(alpha);
  G4double xsc     = GetFresnelDiffuseXsc(theta);
  return xsc;
}

References alpha, and GetFresnelDiffuseXsc().

Referenced by BuildAngleTable().

GetHadronNucleonXscNS()

G4double G4NuclNuclDiffuseElastic::GetHadronNucleonXscNS	(	G4ParticleDefinition *	pParticle,
		G4double	pTkin,
		G4ParticleDefinition *	tParticle
	)

Definition at line 1862 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double xsection(0), /*Delta,*/ A0, B0;
  G4double hpXsc(0);
  G4double hnXsc(0);
 
 
  G4double targ_mass     = tParticle->GetPDGMass();
  G4double proj_mass     = pParticle->GetPDGMass(); 
 
  G4double proj_energy   = proj_mass + pTkin; 
  G4double proj_momentum = std::sqrt(pTkin*(pTkin+2*proj_mass));
 
  G4double sMand = CalcMandelstamS ( proj_mass , targ_mass , proj_momentum );
 
  sMand         /= CLHEP::GeV*CLHEP::GeV;  // in GeV for parametrisation
  proj_momentum /= CLHEP::GeV;
  proj_energy   /= CLHEP::GeV;
  proj_mass     /= CLHEP::GeV;
  G4double logS = G4Log(sMand);
 
  // General PDG fit constants
 
 
  // fEtaRatio=Re[f(0)]/Im[f(0)]
 
  if( proj_momentum >= 1.2 )
  {
    fEtaRatio  = 0.13*(logS - 5.8579332)*G4Pow::GetInstance()->powA(sMand,-0.18);
  }       
  else if( proj_momentum >= 0.6 )
  { 
    fEtaRatio = -75.5*(G4Pow::GetInstance()->powA(proj_momentum,0.25)-0.95)/
      (G4Pow::GetInstance()->powA(3*proj_momentum,2.2)+1);     
  }
  else 
  {
    fEtaRatio = 15.5*proj_momentum/(27*proj_momentum*proj_momentum*proj_momentum+2);
  }
  G4cout<<"fEtaRatio = "<<fEtaRatio<<G4endl;
 
  // xsc
  
  if( proj_momentum >= 10. ) // high energy: pp = nn = np
    // if( proj_momentum >= 2.)
  {
    //Delta = 1.;
 
    //if( proj_energy < 40. ) Delta = 0.916+0.0021*proj_energy;
 
    //AR-12Aug2016  if( proj_momentum >= 10.)
    {
        B0 = 7.5;
        A0 = 100. - B0*G4Log(3.0e7);
 
        xsection = A0 + B0*G4Log(proj_energy) - 11
                  + 103*G4Pow::GetInstance()->powA(2*0.93827*proj_energy + proj_mass*proj_mass+
                     0.93827*0.93827,-0.165);        //  mb
    }
  }
  else // low energy pp = nn != np
  {
      if(pParticle == tParticle) // pp or nn      // nn to be pp
      {
        if( proj_momentum < 0.73 )
        {
          hnXsc = 23 + 50*( G4Pow::GetInstance()->powA( G4Log(0.73/proj_momentum), 3.5 ) );
        }
        else if( proj_momentum < 1.05  )
        {
          hnXsc = 23 + 40*(G4Log(proj_momentum/0.73))*
                         (G4Log(proj_momentum/0.73));
        }
        else  // if( proj_momentum < 10.  )
        {
          hnXsc = 39.0 + 
              75*(proj_momentum - 1.2)/(G4Pow::GetInstance()->powA(proj_momentum,3.0) + 0.15);
        }
        xsection = hnXsc;
      }
      else  // pn to be np
      {
        if( proj_momentum < 0.8 )
        {
          hpXsc = 33+30*G4Pow::GetInstance()->powA(G4Log(proj_momentum/1.3),4.0);
        }      
        else if( proj_momentum < 1.4 )
        {
          hpXsc = 33+30*G4Pow::GetInstance()->powA(G4Log(proj_momentum/0.95),2.0);
        }
        else    // if( proj_momentum < 10.  )
        {
          hpXsc = 33.3+
              20.8*(G4Pow::GetInstance()->powA(proj_momentum,2.0)-1.35)/
                 (G4Pow::GetInstance()->powA(proj_momentum,2.50)+0.95);
        }
        xsection = hpXsc;
      }
  }
  xsection *= CLHEP::millibarn; // parametrised in mb
  G4cout<<"xsection = "<<xsection/CLHEP::millibarn<<" mb"<<G4endl;
  return xsection;
}

References CalcMandelstamS(), fEtaRatio, G4cout, G4endl, G4Log(), G4Pow::GetInstance(), G4ParticleDefinition::GetPDGMass(), CLHEP::GeV, CLHEP::millibarn, and G4Pow::powA().

Referenced by InitParametersGla().

GetIntegrandFunction()

G4double G4NuclNuclDiffuseElastic::GetIntegrandFunction

(

G4double

theta

)

Definition at line 668 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double result;
 
  result  = GetDiffElasticSumProbA(alpha);
 
  // result *= 2*pi*std::sin(theta);
 
  return result;
}

References alpha, and GetDiffElasticSumProbA().

Referenced by IntegralElasticProb(), SampleThetaCMS(), and TestAngleTable().

GetInvCoulombElasticXsc()

G4double G4NuclNuclDiffuseElastic::GetInvCoulombElasticXsc	(	const G4ParticleDefinition *	particle,
		G4double	tMand,
		G4double	momentum,
		G4double	A,
		G4double	Z
	)

Definition at line 348 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double m1 = particle->GetPDGMass();
  G4LorentzVector lv1(0.,0.,plab,std::sqrt(plab*plab+m1*m1));
 
  G4int iZ = static_cast<G4int>(Z+0.5);
  G4int iA = static_cast<G4int>(A+0.5);
  G4ParticleDefinition * theDef = 0;
 
  if      (iZ == 1 && iA == 1) theDef = theProton;
  else if (iZ == 1 && iA == 2) theDef = theDeuteron;
  else if (iZ == 1 && iA == 3) theDef = G4Triton::Triton();
  else if (iZ == 2 && iA == 3) theDef = G4He3::He3();
  else if (iZ == 2 && iA == 4) theDef = theAlpha;
  else theDef = G4ParticleTable::GetParticleTable()->GetIonTable()->GetIon(iZ,iA,0);
 
  G4double tmass = theDef->GetPDGMass();
 
  G4LorentzVector lv(0.0,0.0,0.0,tmass);   
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
  lv1.boost(-bst);
 
  G4ThreeVector p1 = lv1.vect();
  G4double ptot    = p1.mag();
  G4double ptot2 = ptot*ptot;
  G4double cost = 1 - 0.5*std::fabs(tMand)/ptot2;
 
  if( cost >= 1.0 )      cost = 1.0;  
  else if( cost <= -1.0) cost = -1.0;
  
  G4double thetaCMS = std::acos(cost);
 
  G4double sigma = GetCoulombElasticXsc( particle, thetaCMS, ptot, Z );
 
  sigma *= pi/ptot2;
 
  return sigma;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), GetCoulombElasticXsc(), G4IonTable::GetIon(), G4ParticleTable::GetIonTable(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGMass(), G4He3::He3(), CLHEP::Hep3Vector::mag(), pi, theAlpha, theDeuteron, theProton, G4Triton::Triton(), CLHEP::HepLorentzVector::vect(), and Z.

GetInvElasticSumXsc()

G4double G4NuclNuclDiffuseElastic::GetInvElasticSumXsc	(	const G4ParticleDefinition *	particle,
		G4double	tMand,
		G4double	momentum,
		G4double	A,
		G4double	Z
	)

Definition at line 296 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double m1 = particle->GetPDGMass();
 
  G4LorentzVector lv1(0.,0.,plab,std::sqrt(plab*plab+m1*m1));
 
  G4int iZ = static_cast<G4int>(Z+0.5);
  G4int iA = static_cast<G4int>(A+0.5);
 
  G4ParticleDefinition* theDef = 0;
 
  if      (iZ == 1 && iA == 1) theDef = theProton;
  else if (iZ == 1 && iA == 2) theDef = theDeuteron;
  else if (iZ == 1 && iA == 3) theDef = G4Triton::Triton();
  else if (iZ == 2 && iA == 3) theDef = G4He3::He3();
  else if (iZ == 2 && iA == 4) theDef = theAlpha;
  else theDef = G4ParticleTable::GetParticleTable()->GetIonTable()->GetIon(iZ,iA,0);
 
  G4double tmass = theDef->GetPDGMass();
 
  G4LorentzVector lv(0.0,0.0,0.0,tmass);   
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
  lv1.boost(-bst);
 
  G4ThreeVector p1 = lv1.vect();
  G4double ptot    = p1.mag();
  G4double ptot2   = ptot*ptot;
  G4double cost    = 1 - 0.5*std::fabs(tMand)/ptot2;
 
  if( cost >= 1.0 )      cost = 1.0;  
  else if( cost <= -1.0) cost = -1.0;
  
  G4double thetaCMS = std::acos(cost);
 
  G4double sigma = GetDiffuseElasticSumXsc( particle, thetaCMS, ptot, A, Z );
 
  sigma *= pi/ptot2;
 
  return sigma;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), GetDiffuseElasticSumXsc(), G4IonTable::GetIon(), G4ParticleTable::GetIonTable(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGMass(), G4He3::He3(), CLHEP::Hep3Vector::mag(), pi, theAlpha, theDeuteron, theProton, G4Triton::Triton(), CLHEP::HepLorentzVector::vect(), and Z.

GetInvElasticXsc()

G4double G4NuclNuclDiffuseElastic::GetInvElasticXsc	(	const G4ParticleDefinition *	particle,
		G4double	theta,
		G4double	momentum,
		G4double	A,
		G4double	Z
	)

Definition at line 211 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double m1 = particle->GetPDGMass();
  G4LorentzVector lv1(0.,0.,plab,std::sqrt(plab*plab+m1*m1));
 
  G4int iZ = static_cast<G4int>(Z+0.5);
  G4int iA = static_cast<G4int>(A+0.5);
  G4ParticleDefinition * theDef = 0;
 
  if      (iZ == 1 && iA == 1) theDef = theProton;
  else if (iZ == 1 && iA == 2) theDef = theDeuteron;
  else if (iZ == 1 && iA == 3) theDef = G4Triton::Triton();
  else if (iZ == 2 && iA == 3) theDef = G4He3::He3();
  else if (iZ == 2 && iA == 4) theDef = theAlpha;
  else theDef = G4ParticleTable::GetParticleTable()->GetIonTable()->GetIon(iZ,iA,0);
 
  G4double tmass = theDef->GetPDGMass();
 
  G4LorentzVector lv(0.0,0.0,0.0,tmass);   
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
  lv1.boost(-bst);
 
  G4ThreeVector p1 = lv1.vect();
  G4double ptot    = p1.mag();
  G4double ptot2 = ptot*ptot;
  G4double cost = 1 - 0.5*std::fabs(tMand)/ptot2;
 
  if( cost >= 1.0 )      cost = 1.0;  
  else if( cost <= -1.0) cost = -1.0;
  
  G4double thetaCMS = std::acos(cost);
 
  G4double sigma = GetDiffuseElasticXsc( particle, thetaCMS, ptot, A);
 
  sigma *= pi/ptot2;
 
  return sigma;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), GetDiffuseElasticXsc(), G4IonTable::GetIon(), G4ParticleTable::GetIonTable(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGMass(), G4He3::He3(), CLHEP::Hep3Vector::mag(), pi, theAlpha, theDeuteron, theProton, G4Triton::Triton(), CLHEP::HepLorentzVector::vect(), and Z.

GetLegendrePol()

G4double G4NuclNuclDiffuseElastic::GetLegendrePol	(	G4int	n,
		G4double	x
	)

Definition at line 1440 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double legPol, epsilon = 1.e-6;
  G4double x = std::cos(theta);
 
  if     ( n  < 0 ) legPol = 0.;
  else if( n == 0 ) legPol = 1.;
  else if( n == 1 ) legPol = x;
  else if( n == 2 ) legPol = (3.*x*x-1.)/2.;
  else if( n == 3 ) legPol = (5.*x*x*x-3.*x)/2.;
  else if( n == 4 ) legPol = (35.*x*x*x*x-30.*x*x+3.)/8.;
  else if( n == 5 ) legPol = (63.*x*x*x*x*x-70.*x*x*x+15.*x)/8.;
  else if( n == 6 ) legPol = (231.*x*x*x*x*x*x-315.*x*x*x*x+105.*x*x-5.)/16.;
  else           
  {
    // legPol = ( (2*n-1)*x*GetLegendrePol(n-1,x) - (n-1)*GetLegendrePol(n-2,x) )/n;
 
    legPol = std::sqrt( 2./(n*CLHEP::pi*std::sin(theta+epsilon)) )*std::sin( (n+0.5)*theta+0.25*CLHEP::pi );
  }
  return legPol; 
}

References epsilon(), CLHEP::detail::n, and CLHEP::pi.

Referenced by AmplitudeGla().

GetMaxEnergy() [1/2]

G4double G4HadronicInteraction::GetMaxEnergy ( ) const

inlineinherited

Definition at line 96 of file G4HadronicInteraction.hh.

  { return theMaxEnergy; }

References G4HadronicInteraction::theMaxEnergy.

Referenced by G4HadronicInteraction::ActivateFor(), G4IonQMDPhysics::AddProcess(), G4IonINCLXXPhysics::AddProcess(), G4IonPhysics::ConstructProcess(), G4ChargeExchange::G4ChargeExchange(), G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4HadronElastic::G4HadronElastic(), G4hhElastic::G4hhElastic(), G4LowEGammaNuclearModel::G4LowEGammaNuclearModel(), G4NeutrinoElectronCcModel::G4NeutrinoElectronCcModel(), G4NeutrinoElectronNcModel::G4NeutrinoElectronNcModel(), G4NeutronRadCapture::G4NeutronRadCapture(), G4NuclNuclDiffuseElastic(), G4EnergyRangeManager::GetHadronicInteraction(), and G4HadronicProcessStore::Print().

GetMaxEnergy() [2/2]

G4double G4HadronicInteraction::GetMaxEnergy ( const G4Material *  aMaterial, const G4Element *  anElement  ) const

inherited

Definition at line 131 of file G4HadronicInteraction.cc.

{
  if(!IsBlocked()) { return theMaxEnergy; } 
  if( IsBlocked(aMaterial) || IsBlocked(anElement) ) { return 0.0; }
  if(!theMaxEnergyListElements.empty()) {
    for(auto const& elmlist : theMaxEnergyListElements) {
      if( anElement == elmlist.second )
    { return elmlist.first; }
    }
  }
  if(!theMaxEnergyList.empty()) {
    for(auto const& matlist : theMaxEnergyList) {
      if( aMaterial == matlist.second )
    { return matlist.first; }
    }
  }
  return theMaxEnergy;
}

References G4HadronicInteraction::IsBlocked(), G4HadronicInteraction::theMaxEnergy, G4HadronicInteraction::theMaxEnergyList, and G4HadronicInteraction::theMaxEnergyListElements.

GetMinEnergy() [1/2]

G4double G4HadronicInteraction::GetMinEnergy ( ) const

inlineinherited

Definition at line 83 of file G4HadronicInteraction.hh.

  { return theMinEnergy; }

References G4HadronicInteraction::theMinEnergy.

Referenced by G4HadronicInteraction::ActivateFor(), G4EnergyRangeManager::GetHadronicInteraction(), and G4HadronicProcessStore::Print().

GetMinEnergy() [2/2]

G4double G4HadronicInteraction::GetMinEnergy ( const G4Material *  aMaterial, const G4Element *  anElement  ) const

inherited

Definition at line 81 of file G4HadronicInteraction.cc.

{
  if(!IsBlocked()) { return theMinEnergy; } 
  if( IsBlocked(aMaterial) || IsBlocked(anElement) ) { return DBL_MAX; }
  if(!theMinEnergyListElements.empty()) {
    for(auto const& elmlist : theMinEnergyListElements) {
    if( anElement == elmlist.second )
      { return elmlist.first; }
    }
  }
  if(!theMinEnergyList.empty()) {
    for(auto const & matlist : theMinEnergyList) {
      if( aMaterial == matlist.second )
    { return matlist.first; }
    }
  }
  return theMinEnergy;
}

References DBL_MAX, G4HadronicInteraction::IsBlocked(), G4HadronicInteraction::theMinEnergy, G4HadronicInteraction::theMinEnergyList, and G4HadronicInteraction::theMinEnergyListElements.

GetModelName()

const G4String & G4HadronicInteraction::GetModelName ( ) const

inlineinherited

Definition at line 115 of file G4HadronicInteraction.hh.

  { return theModelName; }

References G4HadronicInteraction::theModelName.

Referenced by G4MuMinusCapturePrecompound::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4INCLXXInterface::ApplyYourself(), G4TheoFSGenerator::ApplyYourself(), G4HadronStoppingProcess::AtRestDoIt(), G4VHadronPhysics::BuildModel(), G4HadronicProcess::CheckEnergyMomentumConservation(), G4HadronicProcess::CheckResult(), G4ChargeExchangePhysics::ConstructProcess(), G4MuonicAtomDecay::DecayIt(), G4LENDModel::DumpLENDTargetInfo(), G4AblaInterface::G4AblaInterface(), G4ElectroVDNuclearModel::G4ElectroVDNuclearModel(), G4EMDissociation::G4EMDissociation(), G4ExcitedStringDecay::G4ExcitedStringDecay(), G4LEHadronProtonElastic::G4LEHadronProtonElastic(), G4LENDModel::G4LENDModel(), G4LENDorBERTModel::G4LENDorBERTModel(), G4LEnp::G4LEnp(), G4LEpp::G4LEpp(), G4LFission::G4LFission(), G4LowEGammaNuclearModel::G4LowEGammaNuclearModel(), G4LowEIonFragmentation::G4LowEIonFragmentation(), G4MuonVDNuclearModel::G4MuonVDNuclearModel(), G4NeutrinoElectronCcModel::G4NeutrinoElectronCcModel(), G4NeutrinoNucleusModel::G4NeutrinoNucleusModel(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), G4INCLXXInterface::GetDeExcitationModelName(), G4EnergyRangeManager::GetHadronicInteraction(), G4VHighEnergyGenerator::GetProjectileNucleus(), G4NeutronRadCapture::InitialiseModel(), G4BinaryCascade::ModelDescription(), G4LMsdGenerator::ModelDescription(), G4VPartonStringModel::ModelDescription(), G4TheoFSGenerator::ModelDescription(), G4VHadronPhysics::NewModel(), G4NeutrinoElectronProcess::PostStepDoIt(), G4HadronicProcess::PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), G4MuNeutrinoNucleusProcess::PostStepDoIt(), G4HadronicProcessStore::PrintModelHtml(), G4BinaryCascade::PropagateModelDescription(), G4HadronicProcessStore::RegisterInteraction(), and G4LENDModel::returnUnchanged().

GetNuclearRadius()

G4double G4NuclNuclDiffuseElastic::GetNuclearRadius ( )

inline

Definition at line 185 of file G4NuclNuclDiffuseElastic.hh.

{return fNuclearRadius;};

References fNuclearRadius.

GetProfileLambda()

G4double G4NuclNuclDiffuseElastic::GetProfileLambda ( )

inline

Definition at line 274 of file G4NuclNuclDiffuseElastic.hh.

{return fProfileLambda;};

References fProfileLambda.

GetRatioGen()

G4double G4NuclNuclDiffuseElastic::GetRatioGen

(

G4double

theta

)

Definition at line 1972 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double sinThetaR  = 2.*fHalfRutThetaTg/(1. + fHalfRutThetaTg2);
  G4double dTheta     = 0.5*(theta - fRutherfordTheta);
  G4double sindTheta  = std::sin(dTheta);
 
  G4double prof       = Profile(theta);
  G4double prof2      = prof*prof;
  // G4double profmod    = std::abs(prof);
  G4double order      = std::sqrt(fProfileLambda/sinThetaR/CLHEP::pi)*2.*sindTheta;
 
  order = std::abs(order); // since sin changes sign!
  // G4cout<<"order = "<<order<<G4endl; 
 
  G4double cosFresnel = GetCint(order);  
  G4double sinFresnel = GetSint(order);  
 
  G4double out;
 
  if ( theta <= fRutherfordTheta )
  {
    out  = 1. + 0.5*( (0.5-cosFresnel)*(0.5-cosFresnel)+(0.5-sinFresnel)*(0.5-sinFresnel) )*prof2; 
    out += ( cosFresnel + sinFresnel - 1. )*prof;
  }
  else
  {
    out = 0.5*( (0.5-cosFresnel)*(0.5-cosFresnel)+(0.5-sinFresnel)*(0.5-sinFresnel) )*prof2;
  }
 
  return out;
}

References fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, GetCint(), GetSint(), CLHEP::pi, and Profile().

Referenced by GetFresnelDiffuseXsc().

GetRatioSim()

G4double G4NuclNuclDiffuseElastic::GetRatioSim ( G4double  theta )

inline

Definition at line 1003 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double sinThetaR  = 2.*fHalfRutThetaTg/(1. + fHalfRutThetaTg2);
  G4double dTheta     = 0.5*(theta - fRutherfordTheta);
  G4double sindTheta  = std::sin(dTheta);
 
  G4double order      = std::sqrt(fProfileLambda/sinThetaR/CLHEP::pi)*2.*sindTheta;
  // G4cout<<"order = "<<order<<G4endl;  
  G4double cosFresnel = 0.5 - GetCint(order);  
  G4double sinFresnel = 0.5 - GetSint(order);  
 
  G4double out = 0.5*( cosFresnel*cosFresnel + sinFresnel*sinFresnel );
 
  return out;
}

References fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, GetCint(), GetSint(), and CLHEP::pi.

GetRecoilEnergyThreshold()

G4double G4HadronicInteraction::GetRecoilEnergyThreshold ( ) const

inlineinherited

Definition at line 141 of file G4HadronicInteraction.hh.

  { return recoilEnergyThreshold;}

References G4HadronicInteraction::recoilEnergyThreshold.

Referenced by G4ChargeExchange::ApplyYourself(), and G4HadronElastic::ApplyYourself().

GetRutherfordXsc()

G4double G4NuclNuclDiffuseElastic::GetRutherfordXsc ( G4double  theta )

inline

Definition at line 537 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double sinHalfTheta  = std::sin(0.5*theta);
  G4double sinHalfTheta2 = sinHalfTheta*sinHalfTheta;
 
  G4double ch2           = fRutherfordRatio*fRutherfordRatio;
 
  G4double xsc           = ch2/(sinHalfTheta2+fAm)/(sinHalfTheta2+fAm);
 
  return xsc;
}

References fAm, and fRutherfordRatio.

Referenced by GetFresnelDiffuseXsc().

GetScatteringAngle()

G4double G4NuclNuclDiffuseElastic::GetScatteringAngle	(	G4int	iMomentum,
		G4int	iAngle,
		G4double	position
	)

Definition at line 1086 of file G4NuclNuclDiffuseElastic.cc.

{
 G4double x1, x2, y1, y2, randAngle;
 
  if( iAngle == 0 )
  {
    randAngle = (*fAngleTable)(iMomentum)->GetLowEdgeEnergy(iAngle);
    // iAngle++;
  }
  else
  {
    if ( iAngle >= G4int((*fAngleTable)(iMomentum)->GetVectorLength()) )
    {
      iAngle = (*fAngleTable)(iMomentum)->GetVectorLength() - 1;
    }
    y1 = (*(*fAngleTable)(iMomentum))(iAngle-1);
    y2 = (*(*fAngleTable)(iMomentum))(iAngle);
 
    x1 = (*fAngleTable)(iMomentum)->GetLowEdgeEnergy(iAngle-1);
    x2 = (*fAngleTable)(iMomentum)->GetLowEdgeEnergy(iAngle);
 
    if ( x1 == x2 )   randAngle = x2;
    else
    {
      if ( y1 == y2 ) randAngle = x1 + ( x2 - x1 )*G4UniformRand();
      else
      {
        randAngle = x1 + ( position - y1 )*( x2 - x1 )/( y2 - y1 );
      }
    }
  }
  return randAngle;
}

References fAngleTable, and G4UniformRand.

Referenced by SampleTableThetaCMS().

GetSinHaPit2()

G4double G4NuclNuclDiffuseElastic::GetSinHaPit2 ( G4double  t )

inline

Definition at line 196 of file G4NuclNuclDiffuseElastic.hh.

{return std::sin(CLHEP::halfpi*t*t);};

References CLHEP::halfpi.

Referenced by GetSint().

GetSint()

G4double G4NuclNuclDiffuseElastic::GetSint ( G4double  x )

inline

Definition at line 782 of file G4NuclNuclDiffuseElastic.hh.

{
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
 
  G4double out = 
    integral.Legendre96(this,&G4NuclNuclDiffuseElastic::GetSinHaPit2, 0., x );
 
  return out;
}

References G4NuclNuclDiffuseElastic(), and GetSinHaPit2().

Referenced by GetRatioGen(), and GetRatioSim().

GetSlopeCof()

G4double G4HadronElastic::GetSlopeCof ( const G4int  pdg )

inherited

Definition at line 277 of file G4HadronElastic.cc.

{
  // The input parameter "pdg" should be the absolute value of the PDG code
  // (i.e. the same value for a particle and its antiparticle).
 
  G4double coeff = 1.0;
 
  // heavy barions
 
  static const G4double  lBarCof1S  = 0.88;
  static const G4double  lBarCof2S  = 0.76;
  static const G4double  lBarCof3S  = 0.64;
  static const G4double  lBarCof1C  = 0.784378;
  static const G4double  lBarCofSC  = 0.664378;
  static const G4double  lBarCof2SC = 0.544378;
  static const G4double  lBarCof1B  = 0.740659;
  static const G4double  lBarCofSB  = 0.620659;
  static const G4double  lBarCof2SB = 0.500659;
  
  if( pdg == 3122 || pdg == 3222 ||  pdg == 3112 || pdg == 3212  )
  {
    coeff = lBarCof1S; // Lambda, Sigma+, Sigma-, Sigma0
 
  } else if( pdg == 3322 || pdg == 3312   )
  {
    coeff = lBarCof2S; // Xi-, Xi0
  }
  else if( pdg == 3324)
  {
    coeff = lBarCof3S; // Omega
  }
  else if( pdg == 4122 ||  pdg == 4212 ||   pdg == 4222 ||   pdg == 4112   )
  {
    coeff = lBarCof1C; // LambdaC+, SigmaC+, SigmaC++, SigmaC0
  }
  else if( pdg == 4332 )
  {
    coeff = lBarCof2SC; // OmegaC
  }
  else if( pdg == 4232 || pdg == 4132 )
  {
    coeff = lBarCofSC; // XiC+, XiC0
  }
  else if( pdg == 5122 || pdg == 5222 || pdg == 5112 || pdg == 5212    )
  {
    coeff = lBarCof1B; // LambdaB, SigmaB+, SigmaB-, SigmaB0
  }
  else if( pdg == 5332 )
  {
    coeff = lBarCof2SB; // OmegaB-
  }
  else if( pdg == 5132 || pdg == 5232 ) // XiB-, XiB0
  {
    coeff = lBarCofSB;
  }
  // heavy mesons Kaons?
  static const G4double lMesCof1S = 0.82; // Kp/piP kaons?
  static const G4double llMesCof1C = 0.676568;
  static const G4double llMesCof1B = 0.610989;
  static const G4double llMesCof2C = 0.353135;
  static const G4double llMesCof2B = 0.221978;
  static const G4double llMesCofSC = 0.496568;
  static const G4double llMesCofSB = 0.430989;
  static const G4double llMesCofCB = 0.287557;
  static const G4double llMesCofEtaP = 0.88;
  static const G4double llMesCofEta = 0.76;
 
  if( pdg == 321 || pdg == 311 || pdg == 310 )
  {
    coeff = lMesCof1S; //K+-0
  }
  else if( pdg == 511 ||  pdg == 521  )
  {
    coeff = llMesCof1B; // BMeson0, BMeson+
  }
  else if(pdg == 421 ||  pdg == 411 )
  {
    coeff = llMesCof1C; // DMeson+, DMeson0
  }
  else if( pdg == 531  )
  {
    coeff = llMesCofSB; // BSMeson0
  }
  else if( pdg == 541 )
  {
    coeff = llMesCofCB; // BCMeson+-
  }
  else if(pdg == 431 ) 
  {
    coeff = llMesCofSC; // DSMeson+-
  }
  else if(pdg == 441 || pdg == 443 )
  {
    coeff = llMesCof2C; // Etac, JPsi
  }
  else if(pdg == 553 )
  {
    coeff = llMesCof2B; // Upsilon
  }
  else if(pdg == 221 )
  {
    coeff = llMesCofEta; // Eta
  }
  else if(pdg == 331 )
  {
    coeff = llMesCofEtaP; // Eta'
  } 
  return coeff;
}

GetVerboseLevel()

G4int G4HadronicInteraction::GetVerboseLevel ( ) const

inlineinherited

Definition at line 109 of file G4HadronicInteraction.hh.

  { return verboseLevel; }

References G4HadronicInteraction::verboseLevel.

InitDynParameters()

void G4NuclNuclDiffuseElastic::InitDynParameters	(	const G4ParticleDefinition *	theParticle,
		G4double	partMom
	)

Definition at line 1767 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double a = 0.;
  G4double z = theParticle->GetPDGCharge();
  G4double m1 = theParticle->GetPDGMass();
 
  fWaveVector = partMom/CLHEP::hbarc;
 
  G4double lambda = fCofLambda*fWaveVector*fNuclearRadius;
 
  if( z )
  {
    a           = partMom/m1; // beta*gamma for m1
    fBeta       = a/std::sqrt(1+a*a);
    fZommerfeld = CalculateZommerfeld( fBeta, z, fAtomicNumber);
    fRutherfordRatio = fZommerfeld/fWaveVector; 
    fAm         = CalculateAm( partMom, fZommerfeld, fAtomicNumber);
  }
  fProfileLambda = lambda; // *std::sqrt(1.-2*fZommerfeld/lambda);
  fProfileDelta  = fCofDelta*fProfileLambda;
  fProfileAlpha  = fCofAlpha*fProfileLambda;
 
  CalculateCoulombPhaseZero();
  CalculateRutherfordAnglePar();
 
  return;
}

References CalculateAm(), CalculateCoulombPhaseZero(), CalculateRutherfordAnglePar(), CalculateZommerfeld(), fAm, fAtomicNumber, fBeta, fCofAlpha, fCofDelta, fCofLambda, fNuclearRadius, fProfileAlpha, fProfileDelta, fProfileLambda, fRutherfordRatio, fWaveVector, fZommerfeld, G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), CLHEP::hbarc, and G4InuclParticleNames::lambda.

Referenced by BuildAngleTable(), and SampleCoulombMuCMS().

Initialise()

void G4NuclNuclDiffuseElastic::Initialise

(

)

Definition at line 150 of file G4NuclNuclDiffuseElastic.cc.

{
 
  // fEnergyVector = new G4PhysicsLogVector( theMinEnergy, theMaxEnergy, fEnergyBin );
 
  const G4ElementTable* theElementTable = G4Element::GetElementTable();
  size_t jEl, numOfEl = G4Element::GetNumberOfElements();
 
  // projectile radius
 
  G4double A1 = G4double( fParticle->GetBaryonNumber() );
  G4double R1 = CalculateNuclearRad(A1);
 
  for(jEl = 0 ; jEl < numOfEl; ++jEl) // application element loop
  {
    fAtomicNumber = (*theElementTable)[jEl]->GetZ();     // atomic number
    fAtomicWeight = G4NistManager::Instance()->GetAtomicMassAmu( static_cast< G4int >( fAtomicNumber ) );
 
    fNuclearRadius = CalculateNuclearRad(fAtomicWeight);
    fNuclearRadius += R1;
 
    if(verboseLevel > 0) 
    {   
      G4cout<<"G4NuclNuclDiffuseElastic::Initialise() the element: "
        <<(*theElementTable)[jEl]->GetName()<<G4endl;
    }
    fElementNumberVector.push_back(fAtomicNumber);
    fElementNameVector.push_back((*theElementTable)[jEl]->GetName());
 
    BuildAngleTable();
    fAngleBank.push_back(fAngleTable);
  }  
}

References BuildAngleTable(), CalculateNuclearRad(), fAngleBank, fAngleTable, fAtomicNumber, fAtomicWeight, fElementNameVector, fElementNumberVector, fNuclearRadius, fParticle, G4cout, G4endl, G4NistManager::GetAtomicMassAmu(), G4ParticleDefinition::GetBaryonNumber(), G4Element::GetElementTable(), G4Element::GetNumberOfElements(), G4NistManager::Instance(), and G4HadronicInteraction::verboseLevel.

InitialiseModel()

void G4HadronicInteraction::InitialiseModel ( )

virtualinherited

Reimplemented in G4ANuElNucleusCcModel, G4ANuElNucleusNcModel, G4ANuMuNucleusCcModel, G4ANuMuNucleusNcModel, G4NuElNucleusCcModel, G4NuElNucleusNcModel, G4NuMuNucleusCcModel, G4NuMuNucleusNcModel, G4AblaInterface, G4NeutronRadCapture, G4LowEGammaNuclearModel, G4PreCompoundModel, and G4ElasticHadrNucleusHE.

Definition at line 59 of file G4HadronicInteraction.cc.

{}

InitialiseOnFly()

void G4NuclNuclDiffuseElastic::InitialiseOnFly	(	G4double	Z,
		G4double	A
	)

Definition at line 978 of file G4NuclNuclDiffuseElastic.cc.

{
  fAtomicNumber  = Z;     // atomic number
  fAtomicWeight  = G4NistManager::Instance()->GetAtomicMassAmu( static_cast< G4int >( Z ) );
 
  G4double A1 = G4double( fParticle->GetBaryonNumber() );
  G4double R1 = CalculateNuclearRad(A1);
 
  fNuclearRadius = CalculateNuclearRad(fAtomicWeight) + R1;
  
  if( verboseLevel > 0 )    
  {
    G4cout<<"G4NuclNuclDiffuseElastic::Initialise() the element with Z = "
      <<Z<<"; and A = "<<A<<G4endl;
  }
  fElementNumberVector.push_back(fAtomicNumber);
 
  BuildAngleTable();
 
  fAngleBank.push_back(fAngleTable);
 
  return;
}

References A, BuildAngleTable(), CalculateNuclearRad(), fAngleBank, fAngleTable, fAtomicNumber, fAtomicWeight, fElementNumberVector, fNuclearRadius, fParticle, G4cout, G4endl, G4NistManager::GetAtomicMassAmu(), G4ParticleDefinition::GetBaryonNumber(), G4NistManager::Instance(), G4HadronicInteraction::verboseLevel, and Z.

Referenced by SampleTableThetaCMS().

InitParameters()

void G4NuclNuclDiffuseElastic::InitParameters	(	const G4ParticleDefinition *	theParticle,
		G4double	partMom,
		G4double	Z,
		G4double	A
	)

Definition at line 1722 of file G4NuclNuclDiffuseElastic.cc.

{
  fAtomicNumber  = Z;     // atomic number
  fAtomicWeight  = A;     // number of nucleons
 
  fNuclearRadius2 = CalculateNuclearRad(fAtomicWeight);
  G4double A1     = G4double( theParticle->GetBaryonNumber() );   
  fNuclearRadius1 = CalculateNuclearRad(A1);
  // fNuclearRadius = std::sqrt(fNuclearRadius1*fNuclearRadius1+fNuclearRadius2*fNuclearRadius2);
  fNuclearRadius = fNuclearRadius1 + fNuclearRadius2;
 
  G4double a = 0.;
  G4double z = theParticle->GetPDGCharge();
  G4double m1 = theParticle->GetPDGMass();
 
  fWaveVector = partMom/CLHEP::hbarc;
 
  G4double lambda = fCofLambda*fWaveVector*fNuclearRadius;
  G4cout<<"kR = "<<lambda<<G4endl;
 
  if( z )
  {
    a           = partMom/m1; // beta*gamma for m1
    fBeta       = a/std::sqrt(1+a*a);
    fZommerfeld = CalculateZommerfeld( fBeta, z, fAtomicNumber);
    fRutherfordRatio = fZommerfeld/fWaveVector; 
    fAm         = CalculateAm( partMom, fZommerfeld, fAtomicNumber);
  }
  G4cout<<"fZommerfeld = "<<fZommerfeld<<G4endl;
  fProfileLambda = lambda; // *std::sqrt(1.-2*fZommerfeld/lambda);
  G4cout<<"fProfileLambda = "<<fProfileLambda<<G4endl;
  fProfileDelta  = fCofDelta*fProfileLambda;
  fProfileAlpha  = fCofAlpha*fProfileLambda;
 
  CalculateCoulombPhaseZero();
  CalculateRutherfordAnglePar();
 
  return;
}

References A, CalculateAm(), CalculateCoulombPhaseZero(), CalculateNuclearRad(), CalculateRutherfordAnglePar(), CalculateZommerfeld(), fAm, fAtomicNumber, fAtomicWeight, fBeta, fCofAlpha, fCofDelta, fCofLambda, fNuclearRadius, fNuclearRadius1, fNuclearRadius2, fProfileAlpha, fProfileDelta, fProfileLambda, fRutherfordRatio, fWaveVector, fZommerfeld, G4cout, G4endl, G4ParticleDefinition::GetBaryonNumber(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), CLHEP::hbarc, G4InuclParticleNames::lambda, and Z.

InitParametersGla()

void G4NuclNuclDiffuseElastic::InitParametersGla	(	const G4DynamicParticle *	aParticle,
		G4double	partMom,
		G4double	Z,
		G4double	A
	)

Definition at line 1802 of file G4NuclNuclDiffuseElastic.cc.

{
  fAtomicNumber  = Z;     // target atomic number
  fAtomicWeight  = A;     // target number of nucleons
 
  fNuclearRadius2 = CalculateNuclearRad(fAtomicWeight); // target nucleus radius
  G4double A1     = G4double( aParticle->GetDefinition()->GetBaryonNumber() );   
  fNuclearRadius1 = CalculateNuclearRad(A1); // projectile nucleus radius
  fNuclearRadiusSquare = fNuclearRadius1*fNuclearRadius1+fNuclearRadius2*fNuclearRadius2;
 
 
  G4double a  = 0., kR12;
  G4double z  = aParticle->GetDefinition()->GetPDGCharge();
  G4double m1 = aParticle->GetDefinition()->GetPDGMass();
 
  fWaveVector = partMom/CLHEP::hbarc;
 
  G4double pN = A1 - z;
  if( pN < 0. ) pN = 0.;
 
  G4double tN = A - Z;
  if( tN < 0. ) tN = 0.;
 
  G4double pTkin = aParticle->GetKineticEnergy();  
  pTkin /= A1;
 
 
  fSumSigma = (Z*z+pN*tN)*GetHadronNucleonXscNS(theProton, pTkin, theProton) +
              (z*tN+pN*Z)*GetHadronNucleonXscNS(theProton, pTkin, theNeutron);
 
  G4cout<<"fSumSigma = "<<fSumSigma/CLHEP::millibarn<<" mb"<<G4endl;
  G4cout<<"pi*R2 = "<<CLHEP::pi*fNuclearRadiusSquare/CLHEP::millibarn<<" mb"<<G4endl;
  kR12 = fWaveVector*std::sqrt(fNuclearRadiusSquare);
  G4cout<<"k*sqrt(R2) = "<<kR12<<" "<<G4endl;
  fMaxL = (G4int(kR12)+1)*4;
  G4cout<<"fMaxL = "<<fMaxL<<" "<<G4endl;
 
  if( z )
  {
      a           = partMom/m1; // beta*gamma for m1
      fBeta       = a/std::sqrt(1+a*a);
      fZommerfeld = CalculateZommerfeld( fBeta, z, fAtomicNumber);
      fAm         = CalculateAm( partMom, fZommerfeld, fAtomicNumber);
  }
 
  CalculateCoulombPhaseZero();
 
 
  return;
}

References A, CalculateAm(), CalculateCoulombPhaseZero(), CalculateNuclearRad(), CalculateZommerfeld(), fAm, fAtomicNumber, fAtomicWeight, fBeta, fMaxL, fNuclearRadius1, fNuclearRadius2, fNuclearRadiusSquare, fSumSigma, fWaveVector, fZommerfeld, G4cout, G4endl, G4ParticleDefinition::GetBaryonNumber(), G4DynamicParticle::GetDefinition(), GetHadronNucleonXscNS(), G4DynamicParticle::GetKineticEnergy(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), CLHEP::hbarc, CLHEP::millibarn, CLHEP::pi, theNeutron, theProton, and Z.

IntegralElasticProb()

G4double G4NuclNuclDiffuseElastic::IntegralElasticProb	(	const G4ParticleDefinition *	particle,
		G4double	theta,
		G4double	momentum,
		G4double	A
	)

Definition at line 684 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double result;
  fParticle      = particle;
  fWaveVector    = momentum/hbarc;
  fAtomicWeight  = A;
 
  fNuclearRadius = CalculateNuclearRad(A);
 
 
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
 
  // result = integral.Legendre10(this,&G4NuclNuclDiffuseElastic::GetIntegrandFunction, 0., theta ); 
  result = integral.Legendre96(this,&G4NuclNuclDiffuseElastic::GetIntegrandFunction, 0., theta ); 
 
  return result;
}

References A, CalculateNuclearRad(), fAtomicWeight, fNuclearRadius, fParticle, fWaveVector, G4NuclNuclDiffuseElastic(), GetIntegrandFunction(), and source.hepunit::hbarc.

IsApplicable()

G4bool G4HadronicInteraction::IsApplicable ( const G4HadProjectile &  aTrack, G4Nucleus &  targetNucleus  )

virtualinherited

Reimplemented in G4DiffuseElastic, G4DiffuseElasticV2, G4hhElastic, G4NeutrinoElectronNcModel, G4NeutronElectronElModel, G4ANuElNucleusCcModel, G4ANuElNucleusNcModel, G4ANuMuNucleusCcModel, G4ANuMuNucleusNcModel, G4NeutrinoElectronCcModel, G4NeutrinoNucleusModel, G4NuElNucleusCcModel, G4NuElNucleusNcModel, G4NuMuNucleusCcModel, G4NuMuNucleusNcModel, G4CascadeInterface, and G4LMsdGenerator.

Definition at line 75 of file G4HadronicInteraction.cc.

{ 
  return true;
}

IsBlocked() [1/3]

G4bool G4HadronicInteraction::IsBlocked ( ) const

inlineprotectedinherited

Definition at line 169 of file G4HadronicInteraction.hh.

{ return isBlocked;}

References G4HadronicInteraction::isBlocked.

Referenced by G4HadronicInteraction::GetMaxEnergy(), and G4HadronicInteraction::GetMinEnergy().

IsBlocked() [2/3]

G4bool G4HadronicInteraction::IsBlocked ( const G4Element *  anElement ) const

inherited

Definition at line 202 of file G4HadronicInteraction.cc.

{
  for (auto const& elm : theBlockedListElements) {
    if (anElement == elm) return true;
  }
  return false;
}

References G4HadronicInteraction::theBlockedListElements.

IsBlocked() [3/3]

G4bool G4HadronicInteraction::IsBlocked ( const G4Material *  aMaterial ) const

inherited

Definition at line 193 of file G4HadronicInteraction.cc.

{
  for (auto const& mat : theBlockedList) {
    if (aMaterial == mat) return true;
  }
  return false;
}

References G4HadronicInteraction::theBlockedList.

LowestEnergyLimit()

G4double G4HadronElastic::LowestEnergyLimit ( ) const

inlineinherited

Definition at line 96 of file G4HadronElastic.hh.

{
  return lowestEnergyLimit;
}

References G4HadronElastic::lowestEnergyLimit.

Referenced by G4NeutrinoElectronNcModel::ApplyYourself(), and G4NeutronElectronElModel::ApplyYourself().

ModelDescription()

void G4HadronElastic::ModelDescription ( std::ostream &  outFile ) const

overridevirtualinherited

Reimplemented from G4HadronicInteraction.

Reimplemented in G4NeutrinoElectronNcModel, and G4NeutronElectronElModel.

Definition at line 70 of file G4HadronElastic.cc.

{
  outFile << "G4HadronElastic is the base class for all hadron-nucleus\n" 
          << "elastic scattering models except HP.\n" 
          << "By default it uses the Gheisha two-exponential momentum\n"
      << "transfer parameterization.  The model is fully relativistic\n"
      << "as opposed to the original Gheisha model which was not.\n"
      << "This model may be used for all long-lived hadrons at all\n"
      << "incident energies but fit the data only for relativistic scattering.\n";
}

operator!=()

G4bool G4HadronicInteraction::operator!= ( const G4HadronicInteraction &  right ) const

deleteinherited

operator==()

G4bool G4HadronicInteraction::operator== ( const G4HadronicInteraction &  right ) const

deleteinherited

PhaseFar()

G4complex G4NuclNuclDiffuseElastic::PhaseFar ( G4double  theta )

inline

Definition at line 945 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double twosigma = 2.*fCoulombPhase0; 
  twosigma -= fZommerfeld*G4Log(fHalfRutThetaTg2/(1.+fHalfRutThetaTg2));
  twosigma += fRutherfordTheta*fZommerfeld/fHalfRutThetaTg - CLHEP::halfpi;
  twosigma += fProfileLambda*theta - 0.25*CLHEP::pi;
 
  twosigma *= fCofPhase;
 
  G4complex z = G4complex(0., twosigma);
 
  return std::exp(z);
}

References fCofPhase, fCoulombPhase0, fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, fZommerfeld, G4Log(), CLHEP::halfpi, and CLHEP::pi.

Referenced by AmplitudeFar().

PhaseNear()

G4complex G4NuclNuclDiffuseElastic::PhaseNear ( G4double  theta )

inline

Definition at line 927 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double twosigma = 2.*fCoulombPhase0; 
  twosigma -= fZommerfeld*G4Log(fHalfRutThetaTg2/(1.+fHalfRutThetaTg2));
  twosigma += fRutherfordTheta*fZommerfeld/fHalfRutThetaTg - CLHEP::halfpi;
  twosigma -= fProfileLambda*theta - 0.25*CLHEP::pi;
 
  twosigma *= fCofPhase;
 
  G4complex z = G4complex(0., twosigma);
 
  return std::exp(z);
}

References fCofPhase, fCoulombPhase0, fHalfRutThetaTg, fHalfRutThetaTg2, fProfileLambda, fRutherfordTheta, fZommerfeld, G4Log(), CLHEP::halfpi, and CLHEP::pi.

Referenced by AmplitudeNear().

Profile()

G4double G4NuclNuclDiffuseElastic::Profile ( G4double  theta )

inline

Definition at line 908 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double dTheta = fRutherfordTheta - theta;
  G4double result = 0., argument = 0.;
 
  if(std::abs(dTheta) < 0.001) result = 1.;
  else
  {
    argument = fProfileDelta*dTheta;
    result   = CLHEP::pi*argument;
    result  /= std::sinh(result);
  }
  return result;
}

References fProfileDelta, fRutherfordTheta, and CLHEP::pi.

Referenced by GetRatioGen().

ProfileFar()

G4double G4NuclNuclDiffuseElastic::ProfileFar ( G4double  theta )

inline

Definition at line 892 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double dTheta   = fRutherfordTheta + theta;
  G4double argument = fProfileDelta*dTheta;
 
  G4double result   = CLHEP::pi*argument*G4Exp(fProfileAlpha*argument);
  result           /= std::sinh(CLHEP::pi*argument);
  result           /= dTheta;
 
  return result;
}

References fProfileAlpha, fProfileDelta, fRutherfordTheta, G4Exp(), and CLHEP::pi.

Referenced by AmplitudeFar().

ProfileNear()

G4double G4NuclNuclDiffuseElastic::ProfileNear ( G4double  theta )

inline

Definition at line 871 of file G4NuclNuclDiffuseElastic.hh.

{
  G4double dTheta = fRutherfordTheta - theta;
  G4double result = 0., argument = 0.;
 
  if(std::abs(dTheta) < 0.001) result = fProfileAlpha*fProfileDelta;
  else
  {
    argument = fProfileDelta*dTheta;
    result   = CLHEP::pi*argument*G4Exp(fProfileAlpha*argument);
    result  /= std::sinh(CLHEP::pi*argument);
    result  -= 1.;
    result  /= dTheta;
  }
  return result;
}

References fProfileAlpha, fProfileDelta, fRutherfordTheta, G4Exp(), and CLHEP::pi.

Referenced by AmplitudeNear(), and AmplitudeSim().

SampleCoulombMuCMS()

G4double G4NuclNuclDiffuseElastic::SampleCoulombMuCMS	(	const G4ParticleDefinition *	aParticle,
		G4double	p
	)

Definition at line 809 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double t(0.), rand(0.), mu(0.);
 
  G4double A1 = G4double( aParticle->GetBaryonNumber() );
  G4double R1 = CalculateNuclearRad(A1);
 
  fNuclearRadius  = CalculateNuclearRad(fAtomicWeight);
  fNuclearRadius += R1;
 
  InitDynParameters(fParticle, p);
 
  fCoulombMuC = fHalfRutThetaTg2/(1.+fHalfRutThetaTg2);
 
  rand = G4UniformRand();
 
  // sample (1-cosTheta) in 0 < Theta < ThetaC as Coulomb scattering
 
  mu  = fCoulombMuC*rand*fAm;
  mu /= fAm + fCoulombMuC*( 1. - rand );
 
  t = 4.*p*p*mu;
 
  return t;
}

References CalculateNuclearRad(), fAm, fAtomicWeight, fCoulombMuC, fHalfRutThetaTg2, fNuclearRadius, fParticle, G4UniformRand, G4ParticleDefinition::GetBaryonNumber(), and InitDynParameters().

Referenced by SampleInvariantT().

SampleInvariantT()

G4double G4NuclNuclDiffuseElastic::SampleInvariantT ( const G4ParticleDefinition *  p, G4double  plab, G4int  Z, G4int  A  )

virtual

Reimplemented from G4HadronElastic.

Definition at line 777 of file G4NuclNuclDiffuseElastic.cc.

{
  fParticle = aParticle;
  fAtomicNumber = G4double(Z);
  fAtomicWeight = G4double(A);
 
  G4double t(0.), m1 = fParticle->GetPDGMass();
  G4double totElab = std::sqrt(m1*m1+p*p);
  G4double mass2 = G4NucleiProperties::GetNuclearMass(A, Z);
  G4LorentzVector lv1(p,0.0,0.0,totElab);
  G4LorentzVector  lv(0.0,0.0,0.0,mass2);   
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
  lv1.boost(-bst);
 
  G4ThreeVector p1 = lv1.vect();
  G4double momentumCMS = p1.mag();
 
  // t = SampleTableT( aParticle,  momentumCMS, G4double(Z), G4double(A) ); // sample theta2 in cms
 
  t = SampleCoulombMuCMS( aParticle,  momentumCMS); // sample theta2 in cms
 
 
  return t;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), fAtomicNumber, fAtomicWeight, fParticle, G4NucleiProperties::GetNuclearMass(), G4ParticleDefinition::GetPDGMass(), CLHEP::Hep3Vector::mag(), SampleCoulombMuCMS(), CLHEP::HepLorentzVector::vect(), and Z.

SampleT()

G4double G4NuclNuclDiffuseElastic::SampleT	(	const G4ParticleDefinition *	aParticle,
		G4double	p,
		G4double	A
	)

Definition at line 709 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double theta = SampleThetaCMS( aParticle,  p, A); // sample theta in cms
  G4double t     = 2*p*p*( 1 - std::cos(theta) ); // -t !!!
  return t;
}

References A, and SampleThetaCMS().

Referenced by SampleThetaLab().

SampleTableT()

G4double G4NuclNuclDiffuseElastic::SampleTableT	(	const G4ParticleDefinition *	aParticle,
		G4double	p,
		G4double	Z,
		G4double	A
	)

Definition at line 840 of file G4NuclNuclDiffuseElastic.cc.

{
  G4double alpha = SampleTableThetaCMS( aParticle,  p, Z, A); // sample theta2 in cms
  // G4double t     = 2*p*p*( 1 - std::cos(std::sqrt(alpha)) );             // -t !!!
  G4double t     = p*p*alpha;             // -t !!!
  return t;
}

References A, alpha, SampleTableThetaCMS(), and Z.

SampleTableThetaCMS()

G4double G4NuclNuclDiffuseElastic::SampleTableThetaCMS	(	const G4ParticleDefinition *	aParticle,
		G4double	p,
		G4double	Z,
		G4double	A
	)

Definition at line 855 of file G4NuclNuclDiffuseElastic.cc.

{
  size_t iElement;
  G4int iMomentum, iAngle;  
  G4double randAngle, position, theta1, theta2, E1, E2, W1, W2, W;  
  G4double m1 = particle->GetPDGMass();
 
  for(iElement = 0; iElement < fElementNumberVector.size(); iElement++)
  {
    if( std::fabs(Z - fElementNumberVector[iElement]) < 0.5) break;
  }
  if ( iElement == fElementNumberVector.size() ) 
  {
    InitialiseOnFly(Z,A); // table preparation, if needed
 
    // iElement--;
 
    // G4cout << "G4NuclNuclDiffuseElastic: Element with atomic number " << Z
    // << " is not found, return zero angle" << G4endl;
    // return 0.; // no table for this element
  }
  // G4cout<<"iElement = "<<iElement<<G4endl;
 
  fAngleTable = fAngleBank[iElement];
 
  G4double kinE = std::sqrt(momentum*momentum + m1*m1) - m1;
 
  for( iMomentum = 0; iMomentum < fEnergyBin; iMomentum++)
  {
    // G4cout<<iMomentum<<", kinE = "<<kinE/MeV<<", vectorE = "<<fEnergyVector->GetLowEdgeEnergy(iMomentum)/MeV<<G4endl;
    
    if( kinE < fEnergyVector->GetLowEdgeEnergy(iMomentum) ) break;
  }
  // G4cout<<"iMomentum = "<<iMomentum<<", fEnergyBin -1 = "<<fEnergyBin -1<<G4endl;
 
 
  if ( iMomentum >= fEnergyBin ) iMomentum = fEnergyBin-1;   // kinE is more then theMaxEnergy
  if ( iMomentum < 0 )           iMomentum = 0; // against negative index, kinE < theMinEnergy
 
 
  if (iMomentum == fEnergyBin -1 || iMomentum == 0 )   // the table edges
  {
    position = (*(*fAngleTable)(iMomentum))(fAngleBin-2)*G4UniformRand();
 
    // G4cout<<"position = "<<position<<G4endl;
 
    for(iAngle = 0; iAngle < fAngleBin-1; iAngle++)
    {
      if( position < (*(*fAngleTable)(iMomentum))(iAngle) ) break;
    }
    if (iAngle >= fAngleBin-1) iAngle = fAngleBin-2;
 
    // G4cout<<"iAngle = "<<iAngle<<G4endl;
 
    randAngle = GetScatteringAngle(iMomentum, iAngle, position);
 
    // G4cout<<"randAngle = "<<randAngle<<G4endl;
  }
  else  // kinE inside between energy table edges
  {
    // position = (*(*fAngleTable)(iMomentum))(fAngleBin-2)*G4UniformRand();
    position = (*(*fAngleTable)(iMomentum))(0)*G4UniformRand();
 
    // G4cout<<"position = "<<position<<G4endl;
 
    for(iAngle = 0; iAngle < fAngleBin-1; iAngle++)
    {
      // if( position < (*(*fAngleTable)(iMomentum))(iAngle) ) break;
      if( position > (*(*fAngleTable)(iMomentum))(iAngle) ) break;
    }
    if (iAngle >= fAngleBin-1) iAngle = fAngleBin-2;
 
    // G4cout<<"iAngle = "<<iAngle<<G4endl;
 
    theta2  = GetScatteringAngle(iMomentum, iAngle, position);
 
    // G4cout<<"theta2 = "<<theta2<<G4endl;
 
    E2 = fEnergyVector->GetLowEdgeEnergy(iMomentum);
 
    // G4cout<<"E2 = "<<E2<<G4endl;
    
    iMomentum--;
    
    // position = (*(*fAngleTable)(iMomentum))(fAngleBin-2)*G4UniformRand();
 
    // G4cout<<"position = "<<position<<G4endl;
 
    for(iAngle = 0; iAngle < fAngleBin-1; iAngle++)
    {
      // if( position < (*(*fAngleTable)(iMomentum))(iAngle) ) break;
      if( position > (*(*fAngleTable)(iMomentum))(iAngle) ) break;
    }
    if (iAngle >= fAngleBin-1) iAngle = fAngleBin-2;
    
    theta1  = GetScatteringAngle(iMomentum, iAngle, position);
 
    // G4cout<<"theta1 = "<<theta1<<G4endl;
 
    E1 = fEnergyVector->GetLowEdgeEnergy(iMomentum);
 
    // G4cout<<"E1 = "<<E1<<G4endl;
 
    W  = 1.0/(E2 - E1);
    W1 = (E2 - kinE)*W;
    W2 = (kinE - E1)*W;
 
    randAngle = W1*theta1 + W2*theta2;
    
    // randAngle = theta2;
  }
  // G4double angle = randAngle;
  // if (randAngle > 0.) randAngle /= 2*pi*std::sin(angle);
  // G4cout<<"randAngle = "<<randAngle/degree<<G4endl;
 
  return randAngle;
}

References A, fAngleBank, fAngleBin, fAngleTable, fElementNumberVector, fEnergyBin, fEnergyVector, G4UniformRand, G4PhysicsVector::GetLowEdgeEnergy(), G4ParticleDefinition::GetPDGMass(), GetScatteringAngle(), InitialiseOnFly(), position, and Z.

Referenced by SampleTableT().

SampleThetaCMS()

G4double G4NuclNuclDiffuseElastic::SampleThetaCMS	(	const G4ParticleDefinition *	aParticle,
		G4double	p,
		G4double	A
	)

Definition at line 723 of file G4NuclNuclDiffuseElastic.cc.

{
  G4int i, iMax = 100;  
  G4double norm, result, theta1, theta2, thetaMax, sum = 0.;
 
  fParticle      = particle;
  fWaveVector    = momentum/hbarc;
  fAtomicWeight  = A;
 
  fNuclearRadius = CalculateNuclearRad(A);
  
  thetaMax = 10.174/fWaveVector/fNuclearRadius;
 
  if (thetaMax > pi) thetaMax = pi;
 
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
 
  // result = integral.Legendre10(this,&G4NuclNuclDiffuseElastic::GetIntegrandFunction, 0., theta ); 
  norm = integral.Legendre96(this,&G4NuclNuclDiffuseElastic::GetIntegrandFunction, 0., thetaMax );
 
  norm *= G4UniformRand();
 
  for(i = 1; i <= iMax; i++)
  {
    theta1 = (i-1)*thetaMax/iMax; 
    theta2 = i*thetaMax/iMax;
    sum   += integral.Legendre10(this,&G4NuclNuclDiffuseElastic::GetIntegrandFunction, theta1, theta2);
 
    if ( sum >= norm ) 
    {
      result = 0.5*(theta1 + theta2);
      break;
    }
  }
  if (i > iMax ) result = 0.5*(theta1 + theta2);
 
  G4double sigma = pi*thetaMax/iMax;
 
  result += G4RandGauss::shoot(0.,sigma);
 
  if(result < 0.) result = 0.;
  if(result > thetaMax) result = thetaMax;
 
  return result;
}

References A, CalculateNuclearRad(), fAtomicWeight, fNuclearRadius, fParticle, fWaveVector, G4NuclNuclDiffuseElastic(), G4UniformRand, GetIntegrandFunction(), source.hepunit::hbarc, pi, and G4INCL::DeJongSpin::shoot().

Referenced by SampleT().

SampleThetaLab()

G4double G4NuclNuclDiffuseElastic::SampleThetaLab	(	const G4HadProjectile *	aParticle,
		G4double	tmass,
		G4double	A
	)

Definition at line 1129 of file G4NuclNuclDiffuseElastic.cc.

{
  const G4ParticleDefinition* theParticle = aParticle->GetDefinition();
  G4double m1 = theParticle->GetPDGMass();
  G4double plab = aParticle->GetTotalMomentum();
  G4LorentzVector lv1 = aParticle->Get4Momentum();
  G4LorentzVector lv(0.0,0.0,0.0,tmass);   
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
  lv1.boost(-bst);
 
  G4ThreeVector p1 = lv1.vect();
  G4double ptot    = p1.mag();
  G4double tmax    = 4.0*ptot*ptot;
  G4double t       = 0.0;
 
 
  //
  // Sample t
  //
  
  t = SampleT( theParticle, ptot, A);
 
  // NaN finder
  if(!(t < 0.0 || t >= 0.0)) 
  {
    if (verboseLevel > 0) 
    {
      G4cout << "G4NuclNuclDiffuseElastic:WARNING: A = " << A 
         << " mom(GeV)= " << plab/GeV 
             << " S-wave will be sampled" 
         << G4endl; 
    }
    t = G4UniformRand()*tmax; 
  }
  if(verboseLevel>1)
  {
    G4cout <<" t= " << t << " tmax= " << tmax 
       << " ptot= " << ptot << G4endl;
  }
  // Sampling of angles in CM system
 
  G4double phi  = G4UniformRand()*twopi;
  G4double cost = 1. - 2.0*t/tmax;
  G4double sint;
 
  if( cost >= 1.0 ) 
  {
    cost = 1.0;
    sint = 0.0;
  }
  else if( cost <= -1.0) 
  {
    cost = -1.0;
    sint =  0.0;
  }
  else  
  {
    sint = std::sqrt((1.0-cost)*(1.0+cost));
  }    
  if (verboseLevel>1) 
  {
    G4cout << "cos(t)=" << cost << " std::sin(t)=" << sint << G4endl;
  }
  G4ThreeVector v1(sint*std::cos(phi),sint*std::sin(phi),cost);
  v1 *= ptot;
  G4LorentzVector nlv1(v1.x(),v1.y(),v1.z(),std::sqrt(ptot*ptot + m1*m1));
 
  nlv1.boost(bst); 
 
  G4ThreeVector np1 = nlv1.vect();
 
    // G4double theta = std::acos( np1.z()/np1.mag() );  // degree;
 
  G4double theta = np1.theta();
 
  return theta;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4cout, G4endl, G4UniformRand, G4HadProjectile::Get4Momentum(), G4HadProjectile::GetDefinition(), G4ParticleDefinition::GetPDGMass(), G4HadProjectile::GetTotalMomentum(), GeV, CLHEP::Hep3Vector::mag(), SampleT(), CLHEP::Hep3Vector::theta(), twopi, CLHEP::HepLorentzVector::vect(), G4HadronicInteraction::verboseLevel, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

SetCofAlpha()

void G4NuclNuclDiffuseElastic::SetCofAlpha ( G4double  pa )

inline

Definition at line 281 of file G4NuclNuclDiffuseElastic.hh.

{fCofAlpha = pa;};

References fCofAlpha.

SetCofAlphaCoulomb()

void G4NuclNuclDiffuseElastic::SetCofAlphaCoulomb ( G4double  pa )

inline

Definition at line 283 of file G4NuclNuclDiffuseElastic.hh.

{fCofAlphaCoulomb = pa;};

References fCofAlphaCoulomb.

SetCofAlphaMax()

void G4NuclNuclDiffuseElastic::SetCofAlphaMax ( G4double  pa )

inline

Definition at line 282 of file G4NuclNuclDiffuseElastic.hh.

{fCofAlphaMax = pa;};

References fCofAlphaMax.

SetCofDelta()

void G4NuclNuclDiffuseElastic::SetCofDelta ( G4double  pa )

inline

Definition at line 285 of file G4NuclNuclDiffuseElastic.hh.

{fCofDelta = pa;};

References fCofDelta.

SetCofFar()

void G4NuclNuclDiffuseElastic::SetCofFar ( G4double  pa )

inline

Definition at line 287 of file G4NuclNuclDiffuseElastic.hh.

{fCofFar = pa;};

References fCofFar.

SetCofLambda()

void G4NuclNuclDiffuseElastic::SetCofLambda ( G4double  pa )

inline

Definition at line 279 of file G4NuclNuclDiffuseElastic.hh.

{fCofLambda = pa;};

References fCofLambda.

SetCofPhase()

void G4NuclNuclDiffuseElastic::SetCofPhase ( G4double  pa )

inline

Definition at line 286 of file G4NuclNuclDiffuseElastic.hh.

{fCofPhase = pa;};

References fCofPhase.

SetEnergyMomentumCheckLevels()

void G4HadronicInteraction::SetEnergyMomentumCheckLevels ( G4double  relativeLevel, G4double  absoluteLevel  )

inlineinherited

Definition at line 149 of file G4HadronicInteraction.hh.

  { epCheckLevels.first = relativeLevel;
    epCheckLevels.second = absoluteLevel; }

References G4HadronicInteraction::epCheckLevels.

Referenced by G4BinaryCascade::G4BinaryCascade(), G4CascadeInterface::G4CascadeInterface(), and G4FTFModel::G4FTFModel().

SetEtaRatio()

void G4NuclNuclDiffuseElastic::SetEtaRatio ( G4double  pa )

inline

Definition at line 288 of file G4NuclNuclDiffuseElastic.hh.

{fEtaRatio = pa;};

References fEtaRatio.

SetHEModelLowLimit()

void G4NuclNuclDiffuseElastic::SetHEModelLowLimit ( G4double  value )

inline

Definition at line 377 of file G4NuclNuclDiffuseElastic.hh.

{
  lowEnergyLimitHE = value;
}

References lowEnergyLimitHE.

SetLowestEnergyLimit()

void G4NuclNuclDiffuseElastic::SetLowestEnergyLimit ( G4double  value )

inline

Definition at line 387 of file G4NuclNuclDiffuseElastic.hh.

{
  lowestEnergyLimit = value;
}

References lowestEnergyLimit.

SetMaxEnergy() [1/3]

void G4HadronicInteraction::SetMaxEnergy ( const G4double  anEnergy )

inlineinherited

Definition at line 102 of file G4HadronicInteraction.hh.

  { theMaxEnergy = anEnergy; }

References G4HadronicInteraction::theMaxEnergy.

Referenced by G4HadronicInteraction::ActivateFor(), G4IonINCLXXPhysics::AddProcess(), G4BertiniElectroNuclearBuilder::Build(), G4LENDBertiniGammaElectroNuclearBuilder::Build(), G4NeutronLENDBuilder::Build(), G4NeutronPHPBuilder::Build(), G4AlphaPHPBuilder::Build(), G4BertiniKaonBuilder::Build(), G4BertiniNeutronBuilder::Build(), G4BertiniPiKBuilder::Build(), G4BertiniPionBuilder::Build(), G4BertiniProtonBuilder::Build(), G4BinaryAlphaBuilder::Build(), G4BinaryDeuteronBuilder::Build(), G4BinaryHe3Builder::Build(), G4BinaryNeutronBuilder::Build(), G4BinaryPiKBuilder::Build(), G4BinaryPionBuilder::Build(), G4BinaryProtonBuilder::Build(), G4BinaryTritonBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4FTFBinaryKaonBuilder::Build(), G4FTFBinaryNeutronBuilder::Build(), G4FTFBinaryPionBuilder::Build(), G4FTFBinaryProtonBuilder::Build(), G4FTFPAntiBarionBuilder::Build(), G4FTFPKaonBuilder::Build(), G4FTFPNeutronBuilder::Build(), G4FTFPPiKBuilder::Build(), G4FTFPPionBuilder::Build(), G4FTFPProtonBuilder::Build(), G4He3PHPBuilder::Build(), G4HyperonFTFPBuilder::Build(), G4HyperonQGSPBuilder::Build(), G4INCLXXNeutronBuilder::Build(), G4INCLXXPionBuilder::Build(), G4INCLXXProtonBuilder::Build(), G4PrecoNeutronBuilder::Build(), G4PrecoProtonBuilder::Build(), G4ProtonPHPBuilder::Build(), G4QGSBinaryKaonBuilder::Build(), G4QGSBinaryNeutronBuilder::Build(), G4QGSBinaryPiKBuilder::Build(), G4QGSBinaryPionBuilder::Build(), G4QGSBinaryProtonBuilder::Build(), G4QGSPAntiBarionBuilder::Build(), G4QGSPKaonBuilder::Build(), G4QGSPLundStrFragmProtonBuilder::Build(), G4QGSPNeutronBuilder::Build(), G4QGSPPiKBuilder::Build(), G4QGSPPionBuilder::Build(), G4TritonPHPBuilder::Build(), G4QGSPProtonBuilder::Build(), G4HadronicBuilder::BuildFTFP_BERT(), G4HadronicBuilder::BuildFTFQGSP_BERT(), G4QGSBuilder::BuildModel(), G4VHadronPhysics::BuildModel(), G4HadronicBuilder::BuildQGSP_FTFP_BERT(), G4EmExtraPhysics::ConstructGammaElectroNuclear(), LBE::ConstructHad(), G4EmExtraPhysics::ConstructLENDGammaNuclear(), G4HadronDElasticPhysics::ConstructProcess(), G4HadronElasticPhysics::ConstructProcess(), G4HadronHElasticPhysics::ConstructProcess(), G4IonINCLXXPhysics::ConstructProcess(), G4IonPhysics::ConstructProcess(), G4IonPhysicsPHP::ConstructProcess(), G4IonQMDPhysics::ConstructProcess(), G4ANuElNucleusNcModel::G4ANuElNucleusNcModel(), G4ANuMuNucleusNcModel::G4ANuMuNucleusNcModel(), G4BertiniKaonBuilder::G4BertiniKaonBuilder(), G4BertiniPiKBuilder::G4BertiniPiKBuilder(), G4BertiniPionBuilder::G4BertiniPionBuilder(), G4BinaryCascade::G4BinaryCascade(), G4BinaryPiKBuilder::G4BinaryPiKBuilder(), G4BinaryPionBuilder::G4BinaryPionBuilder(), G4ChargeExchange::G4ChargeExchange(), G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4ElectroVDNuclearModel::G4ElectroVDNuclearModel(), G4EMDissociation::G4EMDissociation(), G4FissLib::G4FissLib(), G4FTFBinaryKaonBuilder::G4FTFBinaryKaonBuilder(), G4FTFBinaryNeutronBuilder::G4FTFBinaryNeutronBuilder(), G4FTFBinaryPiKBuilder::G4FTFBinaryPiKBuilder(), G4FTFBinaryPionBuilder::G4FTFBinaryPionBuilder(), G4FTFBinaryProtonBuilder::G4FTFBinaryProtonBuilder(), G4FTFPAntiBarionBuilder::G4FTFPAntiBarionBuilder(), G4FTFPKaonBuilder::G4FTFPKaonBuilder(), G4FTFPNeutronBuilder::G4FTFPNeutronBuilder(), G4FTFPPiKBuilder::G4FTFPPiKBuilder(), G4FTFPPionBuilder::G4FTFPPionBuilder(), G4FTFPProtonBuilder::G4FTFPProtonBuilder(), G4HadronElastic::G4HadronElastic(), G4HadronicAbsorptionFritiof::G4HadronicAbsorptionFritiof(), G4HadronicAbsorptionFritiofWithBinaryCascade::G4HadronicAbsorptionFritiofWithBinaryCascade(), G4hhElastic::G4hhElastic(), G4HyperonFTFPBuilder::G4HyperonFTFPBuilder(), G4HyperonQGSPBuilder::G4HyperonQGSPBuilder(), G4INCLXXPionBuilder::G4INCLXXPionBuilder(), G4LEHadronProtonElastic::G4LEHadronProtonElastic(), G4LENDModel::G4LENDModel(), G4LEnp::G4LEnp(), G4LEpp::G4LEpp(), G4LFission::G4LFission(), G4LowEGammaNuclearModel::G4LowEGammaNuclearModel(), G4MuonVDNuclearModel::G4MuonVDNuclearModel(), G4NeutrinoElectronCcModel::G4NeutrinoElectronCcModel(), G4NeutrinoElectronNcModel::G4NeutrinoElectronNcModel(), G4NeutrinoNucleusModel::G4NeutrinoNucleusModel(), G4NeutronElectronElModel::G4NeutronElectronElModel(), G4NeutronRadCapture::G4NeutronRadCapture(), G4NuclNuclDiffuseElastic(), G4NuElNucleusNcModel::G4NuElNucleusNcModel(), G4NuMuNucleusNcModel::G4NuMuNucleusNcModel(), G4ParticleHPCapture::G4ParticleHPCapture(), G4ParticleHPElastic::G4ParticleHPElastic(), G4ParticleHPFission::G4ParticleHPFission(), G4ParticleHPInelastic::G4ParticleHPInelastic(), G4ParticleHPThermalScattering::G4ParticleHPThermalScattering(), G4QGSPAntiBarionBuilder::G4QGSPAntiBarionBuilder(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), G4HadronPhysicsFTFP_BERT_HP::Neutron(), G4HadronPhysicsINCLXX::Neutron(), G4HadronPhysicsQGSP_BERT_HP::Neutron(), G4HadronPhysicsQGSP_BIC_HP::Neutron(), G4HadronPhysicsShielding::Neutron(), and G4VHadronPhysics::NewModel().

SetMaxEnergy() [2/3]

void G4HadronicInteraction::SetMaxEnergy ( G4double  anEnergy, const G4Element *  anElement  )

inherited

Definition at line 151 of file G4HadronicInteraction.cc.

{
  Block(); 
  if(!theMaxEnergyListElements.empty()) {
    for(auto & elmlist : theMaxEnergyListElements) {
      if( anElement == elmlist.second ) {
        elmlist.first = anEnergy;
        return;
      }
    }
  }
  theMaxEnergyListElements.push_back(std::pair<G4double, const G4Element *>(anEnergy, anElement));
}

References G4HadronicInteraction::Block(), and G4HadronicInteraction::theMaxEnergyListElements.

SetMaxEnergy() [3/3]

void G4HadronicInteraction::SetMaxEnergy ( G4double  anEnergy, const G4Material *  aMaterial  )

inherited

Definition at line 166 of file G4HadronicInteraction.cc.

{
  Block(); 
  if(!theMaxEnergyList.empty()) {
    for(auto & matlist: theMaxEnergyList) {
      if( aMaterial == matlist.second ) {
    matlist.first = anEnergy;
    return;
      }
    }
  }
  theMaxEnergyList.push_back(std::pair<G4double, const G4Material *>(anEnergy, aMaterial));
}

References G4HadronicInteraction::Block(), and G4HadronicInteraction::theMaxEnergyList.

SetMaxL()

void G4NuclNuclDiffuseElastic::SetMaxL ( G4int  l )

inline

Definition at line 289 of file G4NuclNuclDiffuseElastic.hh.

{fMaxL = l;};

References fMaxL.

SetMinEnergy() [1/3]

void G4HadronicInteraction::SetMinEnergy ( G4double  anEnergy )

inlineinherited

Definition at line 89 of file G4HadronicInteraction.hh.

  { theMinEnergy = anEnergy; }

References G4HadronicInteraction::theMinEnergy.

Referenced by G4HadronicInteraction::ActivateFor(), G4BertiniElectroNuclearBuilder::Build(), G4LENDBertiniGammaElectroNuclearBuilder::Build(), G4NeutronLENDBuilder::Build(), G4NeutronPHPBuilder::Build(), G4AlphaPHPBuilder::Build(), G4BertiniKaonBuilder::Build(), G4BertiniNeutronBuilder::Build(), G4BertiniPiKBuilder::Build(), G4BertiniPionBuilder::Build(), G4BertiniProtonBuilder::Build(), G4BinaryAlphaBuilder::Build(), G4BinaryDeuteronBuilder::Build(), G4BinaryHe3Builder::Build(), G4BinaryNeutronBuilder::Build(), G4BinaryPiKBuilder::Build(), G4BinaryPionBuilder::Build(), G4BinaryProtonBuilder::Build(), G4BinaryTritonBuilder::Build(), G4DeuteronPHPBuilder::Build(), G4FTFBinaryKaonBuilder::Build(), G4FTFBinaryNeutronBuilder::Build(), G4FTFBinaryPiKBuilder::Build(), G4FTFBinaryPionBuilder::Build(), G4FTFBinaryProtonBuilder::Build(), G4FTFPAntiBarionBuilder::Build(), G4FTFPKaonBuilder::Build(), G4FTFPNeutronBuilder::Build(), G4FTFPPiKBuilder::Build(), G4FTFPPionBuilder::Build(), G4FTFPProtonBuilder::Build(), G4He3PHPBuilder::Build(), G4HyperonFTFPBuilder::Build(), G4HyperonQGSPBuilder::Build(), G4INCLXXNeutronBuilder::Build(), G4INCLXXPionBuilder::Build(), G4INCLXXProtonBuilder::Build(), G4PrecoNeutronBuilder::Build(), G4PrecoProtonBuilder::Build(), G4ProtonPHPBuilder::Build(), G4QGSBinaryKaonBuilder::Build(), G4QGSBinaryNeutronBuilder::Build(), G4QGSBinaryPiKBuilder::Build(), G4QGSBinaryPionBuilder::Build(), G4QGSBinaryProtonBuilder::Build(), G4QGSPAntiBarionBuilder::Build(), G4QGSPKaonBuilder::Build(), G4QGSPLundStrFragmProtonBuilder::Build(), G4QGSPNeutronBuilder::Build(), G4QGSPPiKBuilder::Build(), G4QGSPPionBuilder::Build(), G4TritonPHPBuilder::Build(), G4QGSPProtonBuilder::Build(), G4QGSBuilder::BuildModel(), G4VHadronPhysics::BuildModel(), G4EmExtraPhysics::ConstructGammaElectroNuclear(), LBE::ConstructHad(), G4EmExtraPhysics::ConstructLENDGammaNuclear(), G4HadronElasticPhysicsHP::ConstructProcess(), G4HadronElasticPhysicsLEND::ConstructProcess(), G4HadronElasticPhysicsPHP::ConstructProcess(), G4HadronDElasticPhysics::ConstructProcess(), G4HadronElasticPhysics::ConstructProcess(), G4HadronHElasticPhysics::ConstructProcess(), G4IonElasticPhysics::ConstructProcess(), G4IonINCLXXPhysics::ConstructProcess(), G4IonPhysics::ConstructProcess(), G4IonPhysicsPHP::ConstructProcess(), G4IonQMDPhysics::ConstructProcess(), G4ANuElNucleusNcModel::G4ANuElNucleusNcModel(), G4ANuMuNucleusNcModel::G4ANuMuNucleusNcModel(), G4BertiniKaonBuilder::G4BertiniKaonBuilder(), G4BertiniPiKBuilder::G4BertiniPiKBuilder(), G4BertiniPionBuilder::G4BertiniPionBuilder(), G4BinaryCascade::G4BinaryCascade(), G4BinaryPiKBuilder::G4BinaryPiKBuilder(), G4BinaryPionBuilder::G4BinaryPionBuilder(), G4ChargeExchange::G4ChargeExchange(), G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4ElectroVDNuclearModel::G4ElectroVDNuclearModel(), G4EMDissociation::G4EMDissociation(), G4FissLib::G4FissLib(), G4FTFBinaryKaonBuilder::G4FTFBinaryKaonBuilder(), G4FTFBinaryNeutronBuilder::G4FTFBinaryNeutronBuilder(), G4FTFBinaryPiKBuilder::G4FTFBinaryPiKBuilder(), G4FTFBinaryPionBuilder::G4FTFBinaryPionBuilder(), G4FTFBinaryProtonBuilder::G4FTFBinaryProtonBuilder(), G4FTFPAntiBarionBuilder::G4FTFPAntiBarionBuilder(), G4FTFPKaonBuilder::G4FTFPKaonBuilder(), G4FTFPNeutronBuilder::G4FTFPNeutronBuilder(), G4FTFPPiKBuilder::G4FTFPPiKBuilder(), G4FTFPPionBuilder::G4FTFPPionBuilder(), G4FTFPProtonBuilder::G4FTFPProtonBuilder(), G4HadronElastic::G4HadronElastic(), G4HadronicAbsorptionBertini::G4HadronicAbsorptionBertini(), G4HadronicAbsorptionFritiof::G4HadronicAbsorptionFritiof(), G4HadronicAbsorptionFritiofWithBinaryCascade::G4HadronicAbsorptionFritiofWithBinaryCascade(), G4hhElastic::G4hhElastic(), G4HyperonFTFPBuilder::G4HyperonFTFPBuilder(), G4HyperonQGSPBuilder::G4HyperonQGSPBuilder(), G4INCLXXPionBuilder::G4INCLXXPionBuilder(), G4LEHadronProtonElastic::G4LEHadronProtonElastic(), G4LENDModel::G4LENDModel(), G4LEnp::G4LEnp(), G4LEpp::G4LEpp(), G4LFission::G4LFission(), G4LowEGammaNuclearModel::G4LowEGammaNuclearModel(), G4MuonVDNuclearModel::G4MuonVDNuclearModel(), G4NeutrinoElectronCcModel::G4NeutrinoElectronCcModel(), G4NeutrinoElectronNcModel::G4NeutrinoElectronNcModel(), G4NeutrinoNucleusModel::G4NeutrinoNucleusModel(), G4NeutronElectronElModel::G4NeutronElectronElModel(), G4NeutronRadCapture::G4NeutronRadCapture(), G4NuclNuclDiffuseElastic(), G4NuElNucleusNcModel::G4NuElNucleusNcModel(), G4NuMuNucleusNcModel::G4NuMuNucleusNcModel(), G4ParticleHPCapture::G4ParticleHPCapture(), G4ParticleHPElastic::G4ParticleHPElastic(), G4ParticleHPFission::G4ParticleHPFission(), G4ParticleHPInelastic::G4ParticleHPInelastic(), G4ParticleHPThermalScattering::G4ParticleHPThermalScattering(), G4QGSPAntiBarionBuilder::G4QGSPAntiBarionBuilder(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), G4NeutrinoElectronCcModel::IsApplicable(), G4HadronPhysicsFTFP_BERT_HP::Neutron(), G4HadronPhysicsINCLXX::Neutron(), G4HadronPhysicsQGSP_BERT_HP::Neutron(), G4HadronPhysicsQGSP_BIC_HP::Neutron(), G4HadronPhysicsShielding::Neutron(), and G4VHadronPhysics::NewModel().

SetMinEnergy() [2/3]

void G4HadronicInteraction::SetMinEnergy ( G4double  anEnergy, const G4Element *  anElement  )

inherited

Definition at line 101 of file G4HadronicInteraction.cc.

{
  Block(); 
  if(!theMinEnergyListElements.empty()) {
    for(auto & elmlist : theMinEnergyListElements) {
      if( anElement == elmlist.second ) {
    elmlist.first = anEnergy;
    return;
      }
    }
  }
  theMinEnergyListElements.push_back(std::pair<G4double, const G4Element *>(anEnergy, anElement));
}

References G4HadronicInteraction::Block(), and G4HadronicInteraction::theMinEnergyListElements.

SetMinEnergy() [3/3]

void G4HadronicInteraction::SetMinEnergy ( G4double  anEnergy, const G4Material *  aMaterial  )

inherited

Definition at line 116 of file G4HadronicInteraction.cc.

{
  Block(); 
  if(!theMinEnergyList.empty()) {
    for(auto & matlist : theMinEnergyList) {
      if( aMaterial == matlist.second ) {
    matlist.first = anEnergy;
    return;
      }
    }
  }
  theMinEnergyList.push_back(std::pair<G4double, const G4Material *>(anEnergy, aMaterial));
}

References G4HadronicInteraction::Block(), and G4HadronicInteraction::theMinEnergyList.

SetModelName()

void G4HadronicInteraction::SetModelName ( const G4String &  nam )

inlineprotectedinherited

Definition at line 166 of file G4HadronicInteraction.hh.

  { theModelName = nam; }

References G4HadronicInteraction::theModelName.

Referenced by G4ExcitedStringDecay::G4ExcitedStringDecay().

SetNuclearRadiusCof()

void G4NuclNuclDiffuseElastic::SetNuclearRadiusCof ( G4double  r )

inline

Definition at line 290 of file G4NuclNuclDiffuseElastic.hh.

{fNuclearRadiusCof = r;};

References fNuclearRadiusCof.

SetPlabLowLimit()

void G4NuclNuclDiffuseElastic::SetPlabLowLimit ( G4double  value )

inline

Definition at line 372 of file G4NuclNuclDiffuseElastic.hh.

{
  plabLowLimit = value;
}

References plabLowLimit.

SetProfileAlpha()

void G4NuclNuclDiffuseElastic::SetProfileAlpha ( G4double  pa )

inline

Definition at line 278 of file G4NuclNuclDiffuseElastic.hh.

{fProfileAlpha = pa;};

References fProfileAlpha.

SetProfileDelta()

void G4NuclNuclDiffuseElastic::SetProfileDelta ( G4double  pd )

inline

Definition at line 277 of file G4NuclNuclDiffuseElastic.hh.

{fProfileDelta = pd;};

References fProfileDelta.

SetProfileLambda()

void G4NuclNuclDiffuseElastic::SetProfileLambda ( G4double  pl )

inline

Definition at line 276 of file G4NuclNuclDiffuseElastic.hh.

{fProfileLambda = pl;};

References fProfileLambda.

SetQModelLowLimit()

void G4NuclNuclDiffuseElastic::SetQModelLowLimit ( G4double  value )

inline

Definition at line 382 of file G4NuclNuclDiffuseElastic.hh.

{
  lowEnergyLimitQ = value;
}

References lowEnergyLimitQ.

SetRecoilEnergyThreshold()

void G4HadronicInteraction::SetRecoilEnergyThreshold ( G4double  val )

inlineinherited

Definition at line 138 of file G4HadronicInteraction.hh.

  { recoilEnergyThreshold = val; }

References G4HadronicInteraction::recoilEnergyThreshold.

Referenced by G4NeutrinoElectronProcess::PostStepDoIt(), G4ElNeutrinoNucleusProcess::PostStepDoIt(), G4HadronElasticProcess::PostStepDoIt(), and G4MuNeutrinoNucleusProcess::PostStepDoIt().

SetRecoilKinEnergyLimit()

void G4NuclNuclDiffuseElastic::SetRecoilKinEnergyLimit ( G4double  value )

inline

Definition at line 367 of file G4NuclNuclDiffuseElastic.hh.

{
  lowEnergyRecoilLimit = value;
}

References lowEnergyRecoilLimit.

SetVerboseLevel()

void G4HadronicInteraction::SetVerboseLevel ( G4int  value )

inlineinherited

Definition at line 112 of file G4HadronicInteraction.hh.

  { verboseLevel = value; }

References G4HadronicInteraction::verboseLevel.

Referenced by G4CascadeInterface::SetVerboseLevel(), and G4PreCompoundDeexcitation::setVerboseLevel().

TestAngleTable()

void G4NuclNuclDiffuseElastic::TestAngleTable	(	const G4ParticleDefinition *	theParticle,
		G4double	partMom,
		G4double	Z,
		G4double	A
	)

Definition at line 1337 of file G4NuclNuclDiffuseElastic.cc.

{
  fAtomicNumber  = Z;     // atomic number
  fAtomicWeight  = A;     // number of nucleons
  fNuclearRadius = CalculateNuclearRad(fAtomicWeight);
  
     
  
  G4cout<<"G4NuclNuclDiffuseElastic::TestAngleTable() init the element with Z = "
      <<Z<<"; and A = "<<A<<G4endl;
 
  fElementNumberVector.push_back(fAtomicNumber);
 
 
 
 
  G4int i=0, j;
  G4double a = 0., z = theParticle->GetPDGCharge(),  m1 = fParticle->GetPDGMass();
  G4double alpha1=0., alpha2=0., alphaMax=0., alphaCoulomb=0.;
  G4double deltaL10 = 0., deltaL96 = 0., deltaAG = 0.;
  G4double sumL10 = 0.,sumL96 = 0.,sumAG = 0.;
  G4double epsilon = 0.001;
 
  G4Integrator<G4NuclNuclDiffuseElastic,G4double(G4NuclNuclDiffuseElastic::*)(G4double)> integral;
  
  fAngleTable = new G4PhysicsTable(fEnergyBin);
 
  fWaveVector = partMom/hbarc;
 
  G4double kR     = fWaveVector*fNuclearRadius;
  G4double kR2    = kR*kR;
  G4double kRmax  = 10.6; // 10.6, 18, 10.174; ~ 3 maxima of J1 or 15., 25.
  G4double kRcoul = 1.2; // 1.4, 2.5; // on the first slope of J1
 
  alphaMax = kRmax*kRmax/kR2;
 
  if (alphaMax > 4.) alphaMax = 4.;  // vmg05-02-09: was pi2 
 
  alphaCoulomb = kRcoul*kRcoul/kR2;
 
  if( z )
  {
      a           = partMom/m1; // beta*gamma for m1
      fBeta       = a/std::sqrt(1+a*a);
      fZommerfeld = CalculateZommerfeld( fBeta, z, fAtomicNumber);
      fAm         = CalculateAm( partMom, fZommerfeld, fAtomicNumber);
  }
  G4PhysicsFreeVector* angleVector = new G4PhysicsFreeVector(fAngleBin-1);
 
  // G4PhysicsLogVector*  angleBins = new G4PhysicsLogVector( 0.001*alphaMax, alphaMax, fAngleBin );
    
  
  fAddCoulomb = false;
 
  for(j = 1; j < fAngleBin; j++)
  {
      // alpha1 = angleBins->GetLowEdgeEnergy(j-1);
      // alpha2 = angleBins->GetLowEdgeEnergy(j);
 
    alpha1 = alphaMax*(j-1)/fAngleBin;
    alpha2 = alphaMax*( j )/fAngleBin;
 
    if( ( alpha2 > alphaCoulomb ) && z ) fAddCoulomb = true;
 
    deltaL10 = integral.Legendre10(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, alpha1, alpha2);
    deltaL96 = integral.Legendre96(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, alpha1, alpha2);
    deltaAG  = integral.AdaptiveGauss(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, 
                                       alpha1, alpha2,epsilon);
 
      // G4cout<<alpha1<<"\t"<<std::sqrt(alpha1)/degree<<"\t"
      //     <<deltaL10<<"\t"<<deltaL96<<"\t"<<deltaAG<<G4endl;
 
    sumL10 += deltaL10;
    sumL96 += deltaL96;
    sumAG  += deltaAG;
 
    G4cout<<alpha1<<"\t"<<std::sqrt(alpha1)/degree<<"\t"
            <<sumL10<<"\t"<<sumL96<<"\t"<<sumAG<<G4endl;
      
    angleVector->PutValue( j-1 , alpha1, sumL10 ); // alpha2
  }
  fAngleTable->insertAt(i,angleVector);
  fAngleBank.push_back(fAngleTable);
 
  /*
  // Integral over all angle range - Bad accuracy !!!
 
  sumL10 = integral.Legendre10(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, 0., alpha2);
  sumL96 = integral.Legendre96(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, 0., alpha2);
  sumAG  = integral.AdaptiveGauss(this, &G4NuclNuclDiffuseElastic::GetIntegrandFunction, 
                                       0., alpha2,epsilon);
  G4cout<<G4endl;
  G4cout<<alpha2<<"\t"<<std::sqrt(alpha2)/degree<<"\t"
            <<sumL10<<"\t"<<sumL96<<"\t"<<sumAG<<G4endl;
  */
  return;
}

References A, alpha2, CalculateAm(), CalculateNuclearRad(), CalculateZommerfeld(), degree, epsilon(), fAddCoulomb, fAm, fAngleBank, fAngleBin, fAngleTable, fAtomicNumber, fAtomicWeight, fBeta, fElementNumberVector, fEnergyBin, fNuclearRadius, fParticle, fWaveVector, fZommerfeld, G4cout, G4endl, G4NuclNuclDiffuseElastic(), GetIntegrandFunction(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), source.hepunit::hbarc, G4PhysicsTable::insertAt(), G4PhysicsFreeVector::PutValue(), and Z.

TestErfcComp()

G4complex G4NuclNuclDiffuseElastic::TestErfcComp ( G4complex  z, G4int  nMax  )

inline

Definition at line 683 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex miz = G4complex( z.imag(), -z.real() ); 
  G4complex erfcz = 1. - GetErfComp( miz, nMax);
  G4complex w = std::exp(-z*z)*erfcz;
  return w;
}

References GetErfComp().

TestErfcInt()

G4complex G4NuclNuclDiffuseElastic::TestErfcInt ( G4complex  z )

inline

Definition at line 707 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex miz = G4complex( z.imag(), -z.real() ); 
  G4complex erfcz = 1. - GetErfInt( miz); // , nMax);
  G4complex w = std::exp(-z*z)*erfcz;
  return w;
}

References GetErfInt().

TestErfcSer()

G4complex G4NuclNuclDiffuseElastic::TestErfcSer ( G4complex  z, G4int  nMax  )

inline

Definition at line 695 of file G4NuclNuclDiffuseElastic.hh.

{
  G4complex miz = G4complex( z.imag(), -z.real() ); 
  G4complex erfcz = 1. - GetErfSer( miz, nMax);
  G4complex w = std::exp(-z*z)*erfcz;
  return w;
}

References GetErfSer().

ThetaCMStoThetaLab()

G4double G4NuclNuclDiffuseElastic::ThetaCMStoThetaLab	(	const G4DynamicParticle *	aParticle,
		G4double	tmass,
		G4double	thetaCMS
	)

Definition at line 1217 of file G4NuclNuclDiffuseElastic.cc.

{
  const G4ParticleDefinition* theParticle = aParticle->GetDefinition();
  G4double m1 = theParticle->GetPDGMass();
  // G4double plab = aParticle->GetTotalMomentum();
  G4LorentzVector lv1 = aParticle->Get4Momentum();
  G4LorentzVector lv(0.0,0.0,0.0,tmass);   
 
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
 
  lv1.boost(-bst);
 
  G4ThreeVector p1 = lv1.vect();
  G4double ptot    = p1.mag();
 
  G4double phi  = G4UniformRand()*twopi;
  G4double cost = std::cos(thetaCMS);
  G4double sint;
 
  if( cost >= 1.0 ) 
  {
    cost = 1.0;
    sint = 0.0;
  }
  else if( cost <= -1.0) 
  {
    cost = -1.0;
    sint =  0.0;
  }
  else  
  {
    sint = std::sqrt((1.0-cost)*(1.0+cost));
  }    
  if (verboseLevel>1) 
  {
    G4cout << "cos(tcms)=" << cost << " std::sin(tcms)=" << sint << G4endl;
  }
  G4ThreeVector v1(sint*std::cos(phi),sint*std::sin(phi),cost);
  v1 *= ptot;
  G4LorentzVector nlv1(v1.x(),v1.y(),v1.z(),std::sqrt(ptot*ptot + m1*m1));
 
  nlv1.boost(bst); 
 
  G4ThreeVector np1 = nlv1.vect();
 
 
  G4double thetaLab = np1.theta();
 
  return thetaLab;
}

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4cout, G4endl, G4UniformRand, G4DynamicParticle::Get4Momentum(), G4DynamicParticle::GetDefinition(), G4ParticleDefinition::GetPDGMass(), CLHEP::Hep3Vector::mag(), CLHEP::Hep3Vector::theta(), twopi, CLHEP::HepLorentzVector::vect(), G4HadronicInteraction::verboseLevel, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

ThetaLabToThetaCMS()

G4double G4NuclNuclDiffuseElastic::ThetaLabToThetaCMS	(	const G4DynamicParticle *	aParticle,
		G4double	tmass,
		G4double	thetaLab
	)

Definition at line 1278 of file G4NuclNuclDiffuseElastic.cc.

{
  const G4ParticleDefinition* theParticle = aParticle->GetDefinition();
  G4double m1 = theParticle->GetPDGMass();
  G4double plab = aParticle->GetTotalMomentum();
  G4LorentzVector lv1 = aParticle->Get4Momentum();
  G4LorentzVector lv(0.0,0.0,0.0,tmass);   
 
  lv += lv1;
 
  G4ThreeVector bst = lv.boostVector();
 
  // lv1.boost(-bst);
 
  // G4ThreeVector p1 = lv1.vect();
  // G4double ptot    = p1.mag();
 
  G4double phi  = G4UniformRand()*twopi;
  G4double cost = std::cos(thetaLab);
  G4double sint;
 
  if( cost >= 1.0 ) 
  {
    cost = 1.0;
    sint = 0.0;
  }
  else if( cost <= -1.0) 
  {
    cost = -1.0;
    sint =  0.0;
  }
  else  
  {
    sint = std::sqrt((1.0-cost)*(1.0+cost));
  }    
  if (verboseLevel>1) 
  {
    G4cout << "cos(tlab)=" << cost << " std::sin(tlab)=" << sint << G4endl;
  }
  G4ThreeVector v1(sint*std::cos(phi),sint*std::sin(phi),cost);
  v1 *= plab;
  G4LorentzVector nlv1(v1.x(),v1.y(),v1.z(),std::sqrt(plab*plab + m1*m1));
 
  nlv1.boost(-bst); 
 
  G4ThreeVector np1 = nlv1.vect();
 
 
  G4double thetaCMS = np1.theta();
 
  return thetaCMS;
}

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4cout, G4endl, G4UniformRand, G4DynamicParticle::Get4Momentum(), G4DynamicParticle::GetDefinition(), G4ParticleDefinition::GetPDGMass(), G4DynamicParticle::GetTotalMomentum(), CLHEP::Hep3Vector::theta(), twopi, CLHEP::HepLorentzVector::vect(), G4HadronicInteraction::verboseLevel, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Field Documentation

epCheckLevels

std::pair<G4double, G4double> G4HadronicInteraction::epCheckLevels

privateinherited

Definition at line 198 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetEnergyMomentumCheckLevels(), and G4HadronicInteraction::SetEnergyMomentumCheckLevels().

fAddCoulomb

G4bool G4NuclNuclDiffuseElastic::fAddCoulomb

private

Definition at line 337 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetDiffuseElasticSumXsc(), GetDiffuseElasticXsc(), and TestAngleTable().

fAm

G4double G4NuclNuclDiffuseElastic::fAm

private

Definition at line 336 of file G4NuclNuclDiffuseElastic.hh.

Referenced by CalculateAm(), CoulombAmplitude(), G4NuclNuclDiffuseElastic(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetDiffuseElasticSumXsc(), GetRutherfordXsc(), InitDynParameters(), InitParameters(), InitParametersGla(), SampleCoulombMuCMS(), and TestAngleTable().

fAngleBank

std::vector<G4PhysicsTable*> G4NuclNuclDiffuseElastic::fAngleBank

private

Definition at line 316 of file G4NuclNuclDiffuseElastic.hh.

Referenced by Initialise(), InitialiseOnFly(), SampleTableThetaCMS(), TestAngleTable(), and ~G4NuclNuclDiffuseElastic().

fAngleBin

G4int G4NuclNuclDiffuseElastic::fAngleBin

private

Definition at line 312 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), SampleTableThetaCMS(), and TestAngleTable().

fAngleTable

G4PhysicsTable* G4NuclNuclDiffuseElastic::fAngleTable

private

Definition at line 315 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), GetScatteringAngle(), Initialise(), InitialiseOnFly(), SampleTableThetaCMS(), TestAngleTable(), and ~G4NuclNuclDiffuseElastic().

fAtomicNumber

G4double G4NuclNuclDiffuseElastic::fAtomicNumber

private

Definition at line 325 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetDiffuseElasticSumXsc(), InitDynParameters(), Initialise(), InitialiseOnFly(), InitParameters(), InitParametersGla(), SampleInvariantT(), and TestAngleTable().

fAtomicWeight

G4double G4NuclNuclDiffuseElastic::fAtomicWeight

private

Definition at line 324 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetDiffuseElasticSumXsc(), GetDiffuseElasticXsc(), Initialise(), InitialiseOnFly(), InitParameters(), InitParametersGla(), IntegralElasticProb(), SampleCoulombMuCMS(), SampleInvariantT(), SampleThetaCMS(), and TestAngleTable().

fBeta

G4double G4NuclNuclDiffuseElastic::fBeta

private

Definition at line 333 of file G4NuclNuclDiffuseElastic.hh.

Referenced by CalculateParticleBeta(), G4NuclNuclDiffuseElastic(), GetDiffuseElasticSumXsc(), InitDynParameters(), InitParameters(), InitParametersGla(), and TestAngleTable().

fCofAlpha

G4double G4NuclNuclDiffuseElastic::fCofAlpha

private

Definition at line 349 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitDynParameters(), InitParameters(), and SetCofAlpha().

fCofAlphaCoulomb

G4double G4NuclNuclDiffuseElastic::fCofAlphaCoulomb

private

Definition at line 355 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), GetCofAlphaCoulomb(), and SetCofAlphaCoulomb().

fCofAlphaMax

G4double G4NuclNuclDiffuseElastic::fCofAlphaMax

private

Definition at line 354 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), GetCofAlphaMax(), and SetCofAlphaMax().

fCofDelta

G4double G4NuclNuclDiffuseElastic::fCofDelta

private

Definition at line 350 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitDynParameters(), InitParameters(), and SetCofDelta().

fCofFar

G4double G4NuclNuclDiffuseElastic::fCofFar

private

Definition at line 352 of file G4NuclNuclDiffuseElastic.hh.

Referenced by Amplitude(), G4NuclNuclDiffuseElastic(), and SetCofFar().

fCofLambda

G4double G4NuclNuclDiffuseElastic::fCofLambda

private

Definition at line 348 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitDynParameters(), InitParameters(), and SetCofLambda().

fCofPhase

G4double G4NuclNuclDiffuseElastic::fCofPhase

private

Definition at line 351 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), PhaseFar(), PhaseNear(), and SetCofPhase().

fCoulombMuC

G4double G4NuclNuclDiffuseElastic::fCoulombMuC

private

Definition at line 362 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and SampleCoulombMuCMS().

fCoulombPhase0

G4double G4NuclNuclDiffuseElastic::fCoulombPhase0

private

Definition at line 339 of file G4NuclNuclDiffuseElastic.hh.

Referenced by CalculateCoulombPhaseZero(), CoulombAmplitude(), G4NuclNuclDiffuseElastic(), PhaseFar(), and PhaseNear().

fElementNameVector

std::vector<G4String> G4NuclNuclDiffuseElastic::fElementNameVector

private

Definition at line 319 of file G4NuclNuclDiffuseElastic.hh.

Referenced by Initialise().

fElementNumberVector

std::vector<G4double> G4NuclNuclDiffuseElastic::fElementNumberVector

private

Definition at line 318 of file G4NuclNuclDiffuseElastic.hh.

Referenced by Initialise(), InitialiseOnFly(), SampleTableThetaCMS(), and TestAngleTable().

fEnergyBin

G4int G4NuclNuclDiffuseElastic::fEnergyBin

private

Definition at line 311 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), SampleTableThetaCMS(), and TestAngleTable().

fEnergyVector

G4PhysicsLogVector* G4NuclNuclDiffuseElastic::fEnergyVector

private

Definition at line 314 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), SampleTableThetaCMS(), and ~G4NuclNuclDiffuseElastic().

fEtaRatio

G4double G4NuclNuclDiffuseElastic::fEtaRatio

private

Definition at line 359 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeGla(), G4NuclNuclDiffuseElastic(), GetHadronNucleonXscNS(), and SetEtaRatio().

fHalfRutThetaTg

G4double G4NuclNuclDiffuseElastic::fHalfRutThetaTg

private

Definition at line 340 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeSim(), CalculateRutherfordAnglePar(), G4NuclNuclDiffuseElastic(), GammaLess(), GammaMore(), GetRatioGen(), GetRatioSim(), PhaseFar(), and PhaseNear().

fHalfRutThetaTg2

G4double G4NuclNuclDiffuseElastic::fHalfRutThetaTg2

private

Definition at line 341 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeSim(), CalculateRutherfordAnglePar(), G4NuclNuclDiffuseElastic(), GammaLess(), GammaMore(), GetRatioGen(), GetRatioSim(), PhaseFar(), PhaseNear(), and SampleCoulombMuCMS().

fMaxL

G4int G4NuclNuclDiffuseElastic::fMaxL

private

Definition at line 357 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeGG(), AmplitudeGla(), G4NuclNuclDiffuseElastic(), InitParametersGla(), and SetMaxL().

fNuclearRadius

G4double G4NuclNuclDiffuseElastic::fNuclearRadius

private

Definition at line 329 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetDiffElasticProb(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetDiffuseElasticSumXsc(), GetDiffuseElasticXsc(), GetNuclearRadius(), InitDynParameters(), Initialise(), InitialiseOnFly(), InitParameters(), IntegralElasticProb(), SampleCoulombMuCMS(), SampleThetaCMS(), and TestAngleTable().

fNuclearRadius1

G4double G4NuclNuclDiffuseElastic::fNuclearRadius1

private

Definition at line 327 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitParameters(), and InitParametersGla().

fNuclearRadius2

G4double G4NuclNuclDiffuseElastic::fNuclearRadius2

private

Definition at line 328 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitParameters(), and InitParametersGla().

fNuclearRadiusCof

G4double G4NuclNuclDiffuseElastic::fNuclearRadiusCof

private

Definition at line 331 of file G4NuclNuclDiffuseElastic.hh.

Referenced by CalculateNuclearRad(), G4NuclNuclDiffuseElastic(), and SetNuclearRadiusCof().

fNuclearRadiusSquare

G4double G4NuclNuclDiffuseElastic::fNuclearRadiusSquare

private

Definition at line 330 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeGG(), AmplitudeGla(), G4NuclNuclDiffuseElastic(), and InitParametersGla().

fParticle

const G4ParticleDefinition* G4NuclNuclDiffuseElastic::fParticle

private

Definition at line 321 of file G4NuclNuclDiffuseElastic.hh.

Referenced by BuildAngleTable(), G4NuclNuclDiffuseElastic(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetDiffuseElasticSumXsc(), GetDiffuseElasticXsc(), Initialise(), InitialiseOnFly(), IntegralElasticProb(), SampleCoulombMuCMS(), SampleInvariantT(), SampleThetaCMS(), and TestAngleTable().

fProfileAlpha

G4double G4NuclNuclDiffuseElastic::fProfileAlpha

private

Definition at line 346 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitDynParameters(), InitParameters(), ProfileFar(), ProfileNear(), and SetProfileAlpha().

fProfileDelta

G4double G4NuclNuclDiffuseElastic::fProfileDelta

private

Definition at line 345 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), InitDynParameters(), InitParameters(), Profile(), ProfileFar(), ProfileNear(), and SetProfileDelta().

fProfileLambda

G4double G4NuclNuclDiffuseElastic::fProfileLambda

private

Definition at line 344 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeFar(), AmplitudeNear(), AmplitudeSim(), CalculateRutherfordAnglePar(), G4NuclNuclDiffuseElastic(), GammaLess(), GammaMore(), GetProfileLambda(), GetRatioGen(), GetRatioSim(), InitDynParameters(), InitParameters(), PhaseFar(), PhaseNear(), and SetProfileLambda().

fReZ

G4double G4NuclNuclDiffuseElastic::fReZ

private

Definition at line 361 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetErfInt(), GetExpCos(), and GetExpSin().

fRutherfordRatio

G4double G4NuclNuclDiffuseElastic::fRutherfordRatio

private

Definition at line 335 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetRutherfordXsc(), InitDynParameters(), and InitParameters().

fRutherfordTheta

G4double G4NuclNuclDiffuseElastic::fRutherfordTheta

private

Definition at line 342 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeNear(), AmplitudeSim(), BuildAngleTable(), CalculateRutherfordAnglePar(), G4NuclNuclDiffuseElastic(), GammaLess(), GammaMore(), GetRatioGen(), GetRatioSim(), PhaseFar(), PhaseNear(), Profile(), ProfileFar(), and ProfileNear().

fSumSigma

G4double G4NuclNuclDiffuseElastic::fSumSigma

private

Definition at line 358 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeGG(), AmplitudeGla(), G4NuclNuclDiffuseElastic(), and InitParametersGla().

fWaveVector

G4double G4NuclNuclDiffuseElastic::fWaveVector

private

Definition at line 323 of file G4NuclNuclDiffuseElastic.hh.

Referenced by AmplitudeFar(), AmplitudeGG(), AmplitudeGla(), AmplitudeNear(), CoulombAmplitude(), G4NuclNuclDiffuseElastic(), GetDiffElasticProb(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetDiffuseElasticSumXsc(), GetDiffuseElasticXsc(), InitDynParameters(), InitParameters(), InitParametersGla(), IntegralElasticProb(), SampleThetaCMS(), and TestAngleTable().

fZommerfeld

G4double G4NuclNuclDiffuseElastic::fZommerfeld

private

Definition at line 334 of file G4NuclNuclDiffuseElastic.hh.

Referenced by CalculateCoulombPhase(), CalculateCoulombPhaseZero(), CalculateRutherfordAnglePar(), CalculateZommerfeld(), CoulombAmplitude(), G4NuclNuclDiffuseElastic(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetDiffuseElasticSumXsc(), InitDynParameters(), InitParameters(), InitParametersGla(), PhaseFar(), PhaseNear(), and TestAngleTable().

isBlocked

G4bool G4HadronicInteraction::isBlocked

protectedinherited

Definition at line 188 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::Block(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), and G4HadronicInteraction::IsBlocked().

lowEnergyLimitHE

G4double G4NuclNuclDiffuseElastic::lowEnergyLimitHE

private

Definition at line 306 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and SetHEModelLowLimit().

lowEnergyLimitQ

G4double G4NuclNuclDiffuseElastic::lowEnergyLimitQ

private

Definition at line 307 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and SetQModelLowLimit().

lowEnergyRecoilLimit

G4double G4NuclNuclDiffuseElastic::lowEnergyRecoilLimit

private

Definition at line 305 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and SetRecoilKinEnergyLimit().

lowestEnergyLimit

G4double G4NuclNuclDiffuseElastic::lowestEnergyLimit

private

Definition at line 308 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and SetLowestEnergyLimit().

nwarn

G4int G4HadronElastic::nwarn

privateinherited

Definition at line 88 of file G4HadronElastic.hh.

Referenced by G4HadronElastic::ApplyYourself(), and G4HadronElastic::G4HadronElastic().

plabLowLimit

G4double G4NuclNuclDiffuseElastic::plabLowLimit

private

Definition at line 309 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and SetPlabLowLimit().

pLocalTmax

G4double G4HadronElastic::pLocalTmax

protectedinherited

Definition at line 77 of file G4HadronElastic.hh.

Referenced by G4HadronElastic::ApplyYourself(), G4HadronElastic::G4HadronElastic(), G4ElasticHadrNucleusHE::SampleInvariantT(), G4HadronElastic::SampleInvariantT(), and G4LowEHadronElastic::SampleInvariantT().

recoilEnergyThreshold

G4double G4HadronicInteraction::recoilEnergyThreshold

privateinherited

Definition at line 194 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetRecoilEnergyThreshold(), and G4HadronicInteraction::SetRecoilEnergyThreshold().

registry

G4HadronicInteractionRegistry* G4HadronicInteraction::registry

privateinherited

Definition at line 192 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::G4HadronicInteraction(), and G4HadronicInteraction::~G4HadronicInteraction().

secID

G4int G4HadronElastic::secID

protectedinherited

Definition at line 78 of file G4HadronElastic.hh.

Referenced by G4NeutrinoElectronNcModel::ApplyYourself(), G4NeutronElectronElModel::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4LEHadronProtonElastic::ApplyYourself(), G4LEnp::ApplyYourself(), G4LEpp::ApplyYourself(), G4HadronElastic::G4HadronElastic(), G4LEHadronProtonElastic::G4LEHadronProtonElastic(), G4LEnp::G4LEnp(), G4LEpp::G4LEpp(), G4NeutrinoElectronNcModel::G4NeutrinoElectronNcModel(), and G4NeutronElectronElModel::G4NeutronElectronElModel().

theAlpha

G4ParticleDefinition* G4NuclNuclDiffuseElastic::theAlpha

private

Definition at line 300 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetInvCoulombElasticXsc(), GetInvElasticSumXsc(), and GetInvElasticXsc().

theBlockedList

std::vector<const G4Material *> G4HadronicInteraction::theBlockedList

privateinherited

Definition at line 204 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::DeActivateFor(), and G4HadronicInteraction::IsBlocked().

theBlockedListElements

std::vector<const G4Element *> G4HadronicInteraction::theBlockedListElements

privateinherited

Definition at line 205 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::DeActivateFor(), and G4HadronicInteraction::IsBlocked().

theDeuteron

G4ParticleDefinition* G4NuclNuclDiffuseElastic::theDeuteron

private

Definition at line 299 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetInvCoulombElasticXsc(), GetInvElasticSumXsc(), and GetInvElasticXsc().

theMaxEnergy

G4double G4HadronicInteraction::theMaxEnergy

protectedinherited

Definition at line 186 of file G4HadronicInteraction.hh.

Referenced by G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4HadronicInteraction::G4HadronicInteraction(), G4hhElastic::G4hhElastic(), G4NuclNuclDiffuseElastic(), G4HadronicInteraction::GetMaxEnergy(), and G4HadronicInteraction::SetMaxEnergy().

theMaxEnergyList

std::vector<std::pair<G4double, const G4Material *> > G4HadronicInteraction::theMaxEnergyList

privateinherited

Definition at line 201 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetMaxEnergy(), and G4HadronicInteraction::SetMaxEnergy().

theMaxEnergyListElements

std::vector<std::pair<G4double, const G4Element *> > G4HadronicInteraction::theMaxEnergyListElements

privateinherited

Definition at line 203 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetMaxEnergy(), and G4HadronicInteraction::SetMaxEnergy().

theMinEnergy

G4double G4HadronicInteraction::theMinEnergy

protectedinherited

Definition at line 185 of file G4HadronicInteraction.hh.

Referenced by G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4hhElastic::G4hhElastic(), G4NuclNuclDiffuseElastic(), G4HadronicInteraction::GetMinEnergy(), and G4HadronicInteraction::SetMinEnergy().

theMinEnergyList

std::vector<std::pair<G4double, const G4Material *> > G4HadronicInteraction::theMinEnergyList

privateinherited

Definition at line 200 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetMinEnergy(), and G4HadronicInteraction::SetMinEnergy().

theMinEnergyListElements

std::vector<std::pair<G4double, const G4Element *> > G4HadronicInteraction::theMinEnergyListElements

privateinherited

Definition at line 202 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetMinEnergy(), and G4HadronicInteraction::SetMinEnergy().

theModelName

G4String G4HadronicInteraction::theModelName

privateinherited

Definition at line 196 of file G4HadronicInteraction.hh.

Referenced by G4HadronicInteraction::GetModelName(), and G4HadronicInteraction::SetModelName().

theNeutron

G4ParticleDefinition* G4NuclNuclDiffuseElastic::theNeutron

private

Definition at line 298 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), and InitParametersGla().

theParticleChange

G4HadFinalState G4HadronicInteraction::theParticleChange

protectedinherited

Definition at line 172 of file G4HadronicInteraction.hh.

Referenced by G4WilsonAbrasionModel::ApplyYourself(), G4EMDissociation::ApplyYourself(), G4LENDCapture::ApplyYourself(), G4LENDElastic::ApplyYourself(), G4LENDFission::ApplyYourself(), G4LENDInelastic::ApplyYourself(), G4ElectroVDNuclearModel::ApplyYourself(), G4ParticleHPThermalScattering::ApplyYourself(), G4NeutrinoElectronNcModel::ApplyYourself(), G4NeutronElectronElModel::ApplyYourself(), G4LFission::ApplyYourself(), G4ANuElNucleusCcModel::ApplyYourself(), G4ANuElNucleusNcModel::ApplyYourself(), G4ANuMuNucleusCcModel::ApplyYourself(), G4ANuMuNucleusNcModel::ApplyYourself(), G4MuonVDNuclearModel::ApplyYourself(), G4NeutrinoElectronCcModel::ApplyYourself(), G4NuElNucleusCcModel::ApplyYourself(), G4NuElNucleusNcModel::ApplyYourself(), G4NuMuNucleusCcModel::ApplyYourself(), G4NuMuNucleusNcModel::ApplyYourself(), G4QMDReaction::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4LowEGammaNuclearModel::ApplyYourself(), G4ChargeExchange::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4LEHadronProtonElastic::ApplyYourself(), G4LEnp::ApplyYourself(), G4LEpp::ApplyYourself(), G4BinaryCascade::ApplyYourself(), G4CascadeInterface::ApplyYourself(), G4LMsdGenerator::ApplyYourself(), G4ElectroVDNuclearModel::CalculateEMVertex(), G4MuonVDNuclearModel::CalculateEMVertex(), G4ElectroVDNuclearModel::CalculateHadronicVertex(), G4MuonVDNuclearModel::CalculateHadronicVertex(), G4NeutrinoNucleusModel::CoherentPion(), G4CascadeInterface::copyOutputToHadronicResult(), G4BinaryCascade::DebugEpConservation(), G4BinaryCascade::DebugFinalEpConservation(), G4NeutrinoNucleusModel::FinalBarion(), G4NeutrinoNucleusModel::FinalMeson(), G4WilsonAbrasionModel::GetAbradedNucleons(), G4CascadeInterface::NoInteraction(), G4CascadeInterface::Propagate(), G4NeutrinoNucleusModel::RecoilDeexcitation(), G4LEHadronProtonElastic::~G4LEHadronProtonElastic(), G4LEnp::~G4LEnp(), and G4LFission::~G4LFission().

thePionMinus

const G4ParticleDefinition* G4NuclNuclDiffuseElastic::thePionMinus

private

Definition at line 303 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic().

thePionPlus

const G4ParticleDefinition* G4NuclNuclDiffuseElastic::thePionPlus

private

Definition at line 302 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic().

theProton

G4ParticleDefinition* G4NuclNuclDiffuseElastic::theProton

private

Definition at line 297 of file G4NuclNuclDiffuseElastic.hh.

Referenced by G4NuclNuclDiffuseElastic(), GetDiffElasticSumProb(), GetDiffElasticSumProbA(), GetInvCoulombElasticXsc(), GetInvElasticSumXsc(), GetInvElasticXsc(), and InitParametersGla().

verboseLevel

G4int G4HadronicInteraction::verboseLevel

protectedinherited

Definition at line 177 of file G4HadronicInteraction.hh.

Referenced by G4WilsonAbrasionModel::ApplyYourself(), G4EMDissociation::ApplyYourself(), G4LFission::ApplyYourself(), G4MuMinusCapturePrecompound::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4LowEGammaNuclearModel::ApplyYourself(), G4ChargeExchange::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4LEHadronProtonElastic::ApplyYourself(), G4LEnp::ApplyYourself(), G4LEpp::ApplyYourself(), G4CascadeInterface::ApplyYourself(), G4CascadeInterface::checkFinalResult(), G4CascadeInterface::copyOutputToHadronicResult(), G4CascadeInterface::copyOutputToReactionProducts(), G4LENDModel::create_used_target_map(), G4CascadeInterface::createBullet(), G4CascadeInterface::createTarget(), G4ElasticHadrNucleusHE::DefineHadronValues(), G4ElasticHadrNucleusHE::FillData(), G4ElasticHadrNucleusHE::FillFq2(), G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4ElasticHadrNucleusHE::G4ElasticHadrNucleusHE(), G4EMDissociation::G4EMDissociation(), G4hhElastic::G4hhElastic(), G4NuclNuclDiffuseElastic(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), G4ElasticHadrNucleusHE::GetFt(), G4ElasticHadrNucleusHE::GetLightFq2(), G4ElasticHadrNucleusHE::GetQ2_2(), G4HadronicInteraction::GetVerboseLevel(), G4ElasticHadrNucleusHE::HadronNucleusQ2_2(), G4ElasticHadrNucleusHE::HadronProtonQ2(), G4LFission::init(), G4DiffuseElastic::Initialise(), G4DiffuseElasticV2::Initialise(), Initialise(), G4DiffuseElastic::InitialiseOnFly(), G4DiffuseElasticV2::InitialiseOnFly(), InitialiseOnFly(), G4CascadeInterface::makeDynamicParticle(), G4CascadeInterface::NoInteraction(), G4CascadeInterface::Propagate(), G4ElasticHadrNucleusHE::SampleInvariantT(), G4AntiNuclElastic::SampleThetaCMS(), G4DiffuseElastic::SampleThetaLab(), SampleThetaLab(), G4AntiNuclElastic::SampleThetaLab(), G4WilsonAbrasionModel::SetUseAblation(), G4HadronicInteraction::SetVerboseLevel(), G4WilsonAbrasionModel::SetVerboseLevel(), G4DiffuseElastic::ThetaCMStoThetaLab(), G4DiffuseElasticV2::ThetaCMStoThetaLab(), ThetaCMStoThetaLab(), G4DiffuseElastic::ThetaLabToThetaCMS(), G4DiffuseElasticV2::ThetaLabToThetaCMS(), and ThetaLabToThetaCMS().

---

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

• source/processes/hadronic/models/coherent_elastic/include/G4NuclNuclDiffuseElastic.hh
• source/processes/hadronic/models/coherent_elastic/src/G4NuclNuclDiffuseElastic.cc