G4ElasticHadrNucleusHE Class Reference

#include <G4ElasticHadrNucleusHE.hh>

Inheritance diagram for G4ElasticHadrNucleusHE:

Public Member Functions
void	ActivateFor (const G4Element *anElement)
void	ActivateFor (const G4Material *aMaterial)
G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus) override
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
G4double	ComputeMomentumCMS (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
void	DeActivateFor (const G4Element *anElement)
void	DeActivateFor (const G4Material *aMaterial)
	G4ElasticHadrNucleusHE (const G4String &name="hElasticGlauber")
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
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	GetRecoilEnergyThreshold () const
G4double	GetSlopeCof (const G4int pdg)
G4int	GetVerboseLevel () const
void	InitialiseModel () override
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
G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A) override
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
void	SetLowestEnergyLimit (G4double value)
void	SetMaxEnergy (const G4double anEnergy)
void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
void	SetMinEnergy (G4double anEnergy)
void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
void	SetRecoilEnergyThreshold (G4double val)
void	SetVerboseLevel (G4int value)
	~G4ElasticHadrNucleusHE () override

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 Member Functions
void	Binom ()
void	DefineHadronValues (G4int Z)
void	FillData (const G4ParticleDefinition *p, G4int idx, G4int Z)
G4int	FillFq2 (G4int A)
	G4ElasticHadrNucleusHE (const G4ElasticHadrNucleusHE &)
G4double	GetBinomCof (G4int n, G4int m)
G4double	GetFt (G4double Q2)
G4double	GetLightFq2 (G4int Z, G4int A, G4double Q)
G4double	GetQ2_2 (G4int N, G4int Nmax, const std::vector< G4double > &F, G4double rand)
G4double	HadrNucDifferCrSec (G4int A, G4double Q2)
G4double	HadronNucleusQ2_2 (const G4ElasticData *pElD, G4double plabGeV, G4double tmax)
G4double	HadronProtonQ2 (G4double plab, G4double tmax)
void	InFileName (std::ostringstream &, const G4ParticleDefinition *p, G4int Z)
void	InterpolateHN (G4int n, const G4double EnP[], const G4double C0P[], const G4double C1P[], const G4double B0P[], const G4double B1P[])
G4double	LineInterpol (G4double p0, G4double p2, G4double c1, G4double c2, G4double p)
G4ElasticHadrNucleusHE &	operator= (const G4ElasticHadrNucleusHE &right)
void	OutFileName (std::ostringstream &, const G4ParticleDefinition *p, G4int Z)
G4bool	ReadLine (std::ifstream &, std::vector< G4double > &)
void	WriteLine (std::ofstream &, std::vector< G4double > &)

Private Attributes
G4double	aAIm
G4double	aDIm
G4double	Aeff
G4double	BoundaryP [7]
G4double	BoundaryTG [7]
G4double	BoundaryTL [7]
G4double	Coeff0
G4double	Coeff1
G4double	Coeff2
G4double	ConstU
G4double	DDSect2
G4double	DDSect3
G4double	dQ2
G4double	Dtot11
G4double	ekinLowLimit
std::pair< G4double, G4double >	epCheckLevels
char *	fDirectory
G4double	HadrEnergy
G4double	HadrReIm
G4double	HadrSlope
G4double	HadrTot
G4double	hLabMomentum
G4double	hLabMomentum2
G4double	hMass
G4double	hMass2
G4int	iHadrCode
G4int	iHadron
G4int	iHadron1
G4bool	isMaster
G4double	lowestEnergyLimit
G4NistManager *	nistManager
G4int	nwarn
G4double	Pnucl
G4double	Q2max
G4double	R1
G4double	R2
G4double	recoilEnergyThreshold
G4HadronicInteractionRegistry *	registry
G4double	Slope0
G4double	Slope1
G4double	Slope2
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
G4ParticleDefinition *	theProton
G4double	TotP

Static Private Attributes
static G4double	fBinom [240][240] = {{0.0}}
static G4ElasticData *	fElasticData [NHADRONS][ZMAX] = {{nullptr}}
static G4double	fEnergy [NENERGY] = {0.0}
static const G4int	fHadronCode [NHADRONS]
static const G4int	fHadronType [NHADRONS]
static const G4int	fHadronType1 [NHADRONS]
static G4double	fLineF [ONQ2] = {0.0}
static G4double	fLowEdgeEnergy [NENERGY] = {0.0}
static G4bool	fRetrieveFromFile = false
static G4bool	fStoreToFile = false

Detailed Description

Definition at line 94 of file G4ElasticHadrNucleusHE.hh.

Constructor & Destructor Documentation

G4ElasticHadrNucleusHE() [1/2]

G4ElasticHadrNucleusHE::G4ElasticHadrNucleusHE ( const G4String &  name = "hElasticGlauber" )

explicit

Definition at line 244 of file G4ElasticHadrNucleusHE.cc.

  : G4HadronElastic(name), fDirectory(nullptr), isMaster(false)
{
  dQ2 = hMass = hMass2 = hLabMomentum = hLabMomentum2 = HadrEnergy 
    = R1 = R2 = Pnucl = Aeff = HadrTot = HadrSlope = HadrReIm = TotP = DDSect2
    = DDSect3 = ConstU = Slope1 = Slope2 = Coeff1 = Coeff2
    = Slope0 = Coeff0 = aAIm = aDIm = Dtot11 = Q2max = 0.0;
  iHadrCode = iHadron = iHadron1 = 0;
 
  verboseLevel = 0;
  ekinLowLimit = 400.0*CLHEP::MeV;
 
  BoundaryP[0]=9.0; BoundaryTG[0]=5.0;BoundaryTL[0]=0.;
  BoundaryP[1]=20.0;BoundaryTG[1]=1.5;BoundaryTL[1]=0.;
  BoundaryP[2]=5.0; BoundaryTG[2]=1.0;BoundaryTL[2]=1.5;
  BoundaryP[3]=8.0; BoundaryTG[3]=3.0;BoundaryTL[3]=0.;
  BoundaryP[4]=7.0; BoundaryTG[4]=3.0;BoundaryTL[4]=0.;
  BoundaryP[5]=5.0; BoundaryTG[5]=2.0;BoundaryTL[5]=0.;
  BoundaryP[6]=5.0; BoundaryTG[6]=1.5;BoundaryTL[6]=3.0;
 
  nistManager = G4NistManager::Instance();
 
  if(fEnergy[0] == 0.0) {
#ifdef G4MULTITHREADED
    G4MUTEXLOCK(&elasticMutex);
    if(fEnergy[0] == 0.0) {
#endif
      isMaster = true;
      Binom();
      // energy in GeV
      fEnergy[0] = 0.4;
      fEnergy[1] = 0.6;
      fEnergy[2] = 0.8;
      fEnergy[3] = 1.0;
      fLowEdgeEnergy[0] = 0.0;
      fLowEdgeEnergy[1] = 0.5;
      fLowEdgeEnergy[2] = 0.7;
      fLowEdgeEnergy[3] = 0.9;
      G4double f = G4Exp(G4Log(10.)*0.1);
      G4double e = f*f;
      for(G4int i=4; i<NENERGY; ++i) {
    fEnergy[i] = e;
    fLowEdgeEnergy[i] = e/f;
    e *= f*f;
      }
      if(verboseLevel > 0) {
    G4cout << "### G4ElasticHadrNucleusHE: energy points in GeV" << G4endl;
    for(G4int i=0; i<NENERGY; ++i) {
      G4cout << "  " << i << "   " << fLowEdgeEnergy[i] 
         << "  " << fEnergy[i] << G4endl;
    }
      }
#ifdef G4MULTITHREADED
    }
    G4MUTEXUNLOCK(&elasticMutex);
#endif
  }
}

References aAIm, aDIm, Aeff, Binom(), BoundaryP, BoundaryTG, BoundaryTL, Coeff0, Coeff1, Coeff2, ConstU, DDSect2, DDSect3, dQ2, Dtot11, ekinLowLimit, fEnergy, fLowEdgeEnergy, G4cout, G4endl, G4Exp(), G4Log(), G4MUTEXLOCK, G4MUTEXUNLOCK, HadrEnergy, HadrReIm, HadrSlope, HadrTot, hLabMomentum, hLabMomentum2, hMass, hMass2, iHadrCode, iHadron, iHadron1, G4NistManager::Instance(), isMaster, CLHEP::MeV, NENERGY, nistManager, Pnucl, Q2max, R1, R2, Slope0, Slope1, Slope2, TotP, and G4HadronicInteraction::verboseLevel.

~G4ElasticHadrNucleusHE()

G4ElasticHadrNucleusHE::~G4ElasticHadrNucleusHE ( )

override

Definition at line 315 of file G4ElasticHadrNucleusHE.cc.

{
  if(isMaster) {
    for(G4int j = 0; j < NHADRONS; ++j) {
      for(G4int k = 0; k < ZMAX; ++k) {
    G4ElasticData* ptr = fElasticData[j][k]; 
    if(ptr) { 
      delete ptr;
      fElasticData[j][k] = nullptr;
      for(G4int l = j+1; l < NHADRONS; ++l) {
        if(ptr == fElasticData[l][k]) { fElasticData[l][k] = nullptr; }
      }
    }
      }
    }
    delete fDirectory;
    fDirectory = nullptr;
  }
}

References fDirectory, fElasticData, isMaster, NHADRONS, and ZMAX.

G4ElasticHadrNucleusHE() [2/2]

G4ElasticHadrNucleusHE::G4ElasticHadrNucleusHE ( const G4ElasticHadrNucleusHE &  )

private

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().

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.

Binom()

void G4ElasticHadrNucleusHE::Binom ( )

private

Definition at line 1362 of file G4ElasticHadrNucleusHE.cc.

{
  for(G4int N = 0; N < 240; ++N) {
    G4double J = 1.0;
    for(G4int M = 0; M <= N; ++M) {
      G4double Fact1 = 1.0;
      if (N > 0 && N > M && M > 0 ) {
    J *= (G4double)(N-M+1)/(G4double)M;
    Fact1 = J;
      }
      fBinom[N][M] = Fact1;
    }
  }
}

References fBinom, and M.

Referenced by G4ElasticHadrNucleusHE().

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().

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.

{}

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.

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().

DefineHadronValues()

void G4ElasticHadrNucleusHE::DefineHadronValues ( G4int  Z )

private

Definition at line 962 of file G4ElasticHadrNucleusHE.cc.

{
  G4double sHadr = 2.*HadrEnergy*protonM+protonM2+hMass2;
  G4double sqrS  = std::sqrt(sHadr);
  
  if(verboseLevel>2) {
    G4cout << "GetHadrValues: Z= " << Z << " iHadr= " << iHadron 
       << " E(GeV)= " << HadrEnergy << " sqrS= " << sqrS
       << " plab= " << hLabMomentum   
       <<"  E - m  "<<HadrEnergy - hMass<< G4endl;
  }
  G4double TotN = 0.0;
  G4double logE = G4Log(HadrEnergy);
  G4double logS = G4Log(sHadr);
           TotP = 0.0;
 
  switch (iHadron) {
  case 0:                  //  proton, neutron
  case 6:
 
    if(hLabMomentum > 10) {
      TotP = TotN = 7.5*logE - 40.12525 + 103*G4Exp(-logS*0.165);// mb
 
    } else {
      // ==================  neutron  ================
 
      if( hLabMomentum > 1.4 ) {
    TotN = 33.3+15.2*(hLabMomentum2-1.35)/
      (G4Exp(G4Log(hLabMomentum)*2.37)+0.95);
        
      } else if(hLabMomentum > 0.8) {
    G4double A0 = logE + 0.0513;
    TotN = 33.0 + 25.5*A0*A0;  
      } else {
    G4double A0 = logE - 0.2634;  // log(1.3)
    TotN = 33.0 + 30.*A0*A0*A0*A0;
      }
      //  =================  proton  ===============
 
      if(hLabMomentum >= 1.05) {
    TotP = 39.0+75.*(hLabMomentum-1.2)/(hLabMomentum2*hLabMomentum+0.15);
      } else if(hLabMomentum >= 0.7) {
    G4double A0 = logE + 0.3147;
    TotP = 23.0 + 40.*A0*A0;
      } else {
    TotP = 23.+50.*G4Exp(G4Log(G4Log(0.73/hLabMomentum))*3.5);
      }
    }
    HadrTot = 0.5*(TotP+TotN);
    //  ...................................................
    //  Proton slope
    if(hLabMomentum >= 2.)       { HadrSlope = 5.44 + 0.88*logS; }
    else if(hLabMomentum >= 0.5) { HadrSlope = 3.73*hLabMomentum-0.37; }
    else                         { HadrSlope = 1.5; }
 
    //  ...................................................
    if(hLabMomentum >= 1.2) {
      HadrReIm  = 0.13*(logS - 5.8579332)*G4Exp(-logS*0.18);
    } else if(hLabMomentum >= 0.6) { 
      HadrReIm = -75.5*(G4Exp(G4Log(hLabMomentum)*0.25)-0.95)/
    (G4Exp(G4Log(3*hLabMomentum)*2.2)+1);     
    } else {
      HadrReIm = 15.5*hLabMomentum/(27*hLabMomentum2*hLabMomentum+2);
    }
    //  ...................................................
    DDSect2   = 2.2;                              //mb*GeV-2
    DDSect3   = 0.6;                               //mb*GeV-2
    //  ================== lambda  ==================
    if( iHadrCode == 3122) {
      HadrTot   *= 0.88;
      HadrSlope *=0.85;
      //  ================== sigma +  ==================
    } else if( iHadrCode == 3222) {
      HadrTot   *=0.81;
      HadrSlope *=0.85;
      //  ================== sigma 0,-  ==================
    } else if(iHadrCode == 3112 || iHadrCode == 3212 ) {
      HadrTot   *=0.88;
      HadrSlope *=0.85;
      //  ===================  xi  =================
    } else if( iHadrCode == 3312 || iHadrCode == 3322 ) {
      HadrTot   *=0.77;
      HadrSlope *=0.75;
      //  =================  omega  =================
    } else if( iHadrCode == 3334) {
      HadrTot   *=0.78;
      HadrSlope *=0.7;
    }
    break;
    //  ===========================================================
  case 1:              //   antiproton
  case 7:              //   antineutron
 
    HadrTot   = 5.2+5.2*logE + 123.2/sqrS;     //  mb
    HadrSlope = 8.32+0.57*logS;                //(GeV/c)^-2
 
    if( HadrEnergy < 1000 ) {
      HadrReIm  = 0.06*(sqrS-2.236)*(sqrS-14.14)*G4Exp(-logS*0.8);
    } else {
      HadrReIm  = 0.6*(logS - 5.8579332)*G4Exp(-logS*0.25);
    }
    DDSect2   = 11;                            //mb*(GeV/c)^-2
    DDSect3   = 3;                             //mb*(GeV/c)^-2
    //  ================== lambda  ==================
    if( iHadrCode == -3122) {
      HadrTot   *= 0.88;
      HadrSlope *=0.85;
      //  ================== sigma +  ==================
    } else if( iHadrCode == -3222) {
      HadrTot   *=0.81;
      HadrSlope *=0.85;
      //  ================== sigma 0,-  ==================
    } else if(iHadrCode == -3112 || iHadrCode == -3212 ) {
      HadrTot   *=0.88;
      HadrSlope *=0.85;
    //  ===================  xi  =================
    } else if( iHadrCode == -3312 || iHadrCode == -3322 ) {
      HadrTot   *=0.77;
      HadrSlope *=0.75;
      //  =================  omega  =================
    } else if( iHadrCode == -3334) {
      HadrTot   *=0.78;
      HadrSlope *=0.7;
    }
    break;
    //  -------------------------------------------
  case 2:             //   pi plus, pi minus
  case 3:
 
    if(hLabMomentum >= 3.5) {
      TotP = 10.6+2.*logE + 25.*G4Exp(-logE*0.43); // mb
      //  =========================================
    } else if(hLabMomentum >= 1.15) {
      G4double x = (hLabMomentum - 2.55)/0.55; 
      G4double y = (hLabMomentum - 1.47)/0.225;
      TotP = 3.2*G4Exp(-x*x) + 12.*G4Exp(-y*y) + 27.5;
      //  =========================================
    } else if(hLabMomentum >= 0.4) {
      TotP  = 88*(logE+0.2877)*(logE+0.2877)+14.0;
    //  =========================================
    } else {
      G4double x = (hLabMomentum - 0.29)/0.085;
      TotP = 20. + 180.*G4Exp(-x*x);
    }
    //  -------------------------------------------
 
    if(hLabMomentum >= 3.0 ) {
      TotN = 10.6 + 2.*logE + 30.*G4Exp(-logE*0.43); // mb
    } else if(hLabMomentum >= 1.3) {
      G4double x = (hLabMomentum - 2.1)/0.4;
      G4double y = (hLabMomentum - 1.4)/0.12;
      TotN = 36.1+0.079 - 4.313*logE + 3.*G4Exp(-x*x) + 1.5*G4Exp(-y*y);
    } else if(hLabMomentum >= 0.65) {
      G4double x = (hLabMomentum - 0.72)/0.06;
      G4double y = (hLabMomentum - 1.015)/0.075;
      TotN = 36.1 + 10.*G4Exp(-x*x) + 24*G4Exp(-y*y);
    } else if(hLabMomentum >= 0.37) {
      G4double x = G4Log(hLabMomentum/0.48);
      TotN = 26. + 110.*x*x;
    } else {
      G4double x = (hLabMomentum - 0.29)/0.07;
      TotN = 28.0 + 40.*G4Exp(-x*x);
    }
    HadrTot = (TotP+TotN)*0.5;
    //  ........................................
    HadrSlope = 7.28+0.245*logS;        // GeV-2
    HadrReIm  = 0.2*(logS - 4.6051702)*G4Exp(-logS*0.15);
 
    DDSect2   = 0.7;                               //mb*GeV-2
    DDSect3   = 0.27;                              //mb*GeV-2
 
    break;
    //  ==========================================================
  case 4:            //  K plus
 
    HadrTot   = 10.6+1.8*logE + 9.0*G4Exp(-logE*0.55);  // mb
    if(HadrEnergy>100) { HadrSlope = 15.0; }
    else { HadrSlope = 1.0+1.76*logS - 2.84/sqrS; }   // GeV-2
 
    HadrReIm  = 0.4*(sHadr-20)*(sHadr-150)*G4Exp(-G4Log(sHadr+50)*2.1);
    DDSect2   = 0.7;                             //mb*GeV-2
    DDSect3   = 0.21;                            //mb*GeV-2
    break;
    //  =========================================================
  case 5:              //   K minus
 
    HadrTot   = 10+1.8*logE + 25./sqrS; // mb
    HadrSlope = 6.98+0.127*logS;         // GeV-2
    HadrReIm  = 0.4*(sHadr-20)*(sHadr-20)*G4Exp(-G4Log(sHadr+50)*2.1);
    DDSect2   = 0.7;                             //mb*GeV-2
    DDSect3   = 0.27;                            //mb*GeV-2
    break;
  }   
  //  =========================================================
  if(verboseLevel>2) {
    G4cout << "HadrTot= " << HadrTot << " HadrSlope= " << HadrSlope
       << " HadrReIm= " << HadrReIm << " DDSect2= " << DDSect2
       << " DDSect3= " << DDSect3 << G4endl;
  }
  if(Z != 1) return;
 
  // Scattering of protons
 
  Coeff0 = Coeff1 = Coeff2 = 0.0;
  Slope0 = Slope1 = 1.0;
  Slope2 = 5.0;
 
  // data for iHadron=0
  static const G4double EnP0[6]={1.5,3.0,5.0,9.0,14.0,19.0};
  static const G4double C0P0[6]={0.15,0.02,0.06,0.08,0.0003,0.0002};
  static const G4double C1P0[6]={0.05,0.02,0.03,0.025,0.0,0.0};
  static const G4double B0P0[6]={1.5,2.5,3.0,4.5,1.4,1.25};
  static const G4double B1P0[6]={5.0,1.0,3.5,4.0,4.8,4.8};
      
  // data for iHadron=6,7
  static const G4double EnN[5]={1.5,5.0,10.0,14.0,20.0};
  static const G4double C0N[5]={0.0,0.0,0.02,0.02,0.01};
  static const G4double C1N[5]={0.06,0.008,0.0015,0.001,0.0003};
  static const G4double B0N[5]={1.5,2.5,3.8,3.8,3.5};
  static const G4double B1N[5]={1.5,2.2,3.6,4.5,4.8};
 
  // data for iHadron=1
  static const G4double EnP[2]={1.5,4.0};
  static const G4double C0P[2]={0.001,0.0005};
  static const G4double C1P[2]={0.003,0.001};
  static const G4double B0P[2]={2.5,4.5};
  static const G4double B1P[2]={1.0,4.0};
 
  // data for iHadron=2
  static const G4double EnPP[4]={1.0,2.0,3.0,4.0};
  static const G4double C0PP[4]={0.0,0.0,0.0,0.0};
  static const G4double C1PP[4]={0.15,0.08,0.02,0.01};
  static const G4double B0PP[4]={1.5,2.8,3.8,3.8};
  static const G4double B1PP[4]={0.8,1.6,3.6,4.6};
 
  // data for iHadron=3
  static const G4double EnPPN[4]={1.0,2.0,3.0,4.0};
  static const G4double C0PPN[4]={0.0,0.0,0.0,0.0};
  static const G4double C1PPN[4]={0.0,0.0,0.0,0.0};
  static const G4double B0PPN[4]={1.5,2.8,3.8,3.8};
  static const G4double B1PPN[4]={0.8,1.6,3.6,4.6};
 
  // data for iHadron=4
  static const G4double EnK[4]={1.4,2.33,3.0,5.0};
  static const G4double C0K[4]={0.0,0.0,0.0,0.0};
  static const G4double C1K[4]={0.01,0.007,0.005,0.003};
  static const G4double B0K[4]={1.5,2.0,3.8,3.8};
  static const G4double B1K[4]={1.6,1.6,1.6,1.6};
 
  // data for iHadron=5
  static const G4double EnKM[2]={1.4,4.0};
  static const G4double C0KM[2]={0.006,0.002};
  static const G4double C1KM[2]={0.00,0.00};
  static const G4double B0KM[2]={2.5,3.5};
  static const G4double B1KM[2]={1.6,1.6};
 
  switch(iHadron) {
  case 0:
 
    if(hLabMomentum <BoundaryP[0]) {
      InterpolateHN(6,EnP0,C0P0,C1P0,B0P0,B1P0);
    }
    Coeff2 = 0.8/hLabMomentum2;
    break; 
 
  case 6:
 
    if(hLabMomentum < BoundaryP[1]) {
      InterpolateHN(5,EnN,C0N,C1N,B0N,B1N);
    }
    Coeff2 = 0.8/hLabMomentum2;
    break; 
 
  case 1:
  case 7:
    if(hLabMomentum <  BoundaryP[2]) {
      InterpolateHN(2,EnP,C0P,C1P,B0P,B1P);
    }
    break; 
 
  case 2:
 
    if(hLabMomentum < BoundaryP[3]) {
      InterpolateHN(4,EnPP,C0PP,C1PP,B0PP,B1PP);
    }
    Coeff2 = 0.02/hLabMomentum;
    break; 
 
  case 3:
 
    if(hLabMomentum < BoundaryP[4]) {
      InterpolateHN(4,EnPPN,C0PPN,C1PPN,B0PPN,B1PPN);
    }
    Coeff2 = 0.02/hLabMomentum;
    break;
 
  case 4:
 
    if(hLabMomentum < BoundaryP[5]) {
      InterpolateHN(4,EnK,C0K,C1K,B0K,B1K);
    }
    if(hLabMomentum < 1) { Coeff2 = 0.34; }
    else  { Coeff2 = 0.34/(hLabMomentum2*hLabMomentum); }
    break; 
 
  case 5:
    if(hLabMomentum < BoundaryP[6]) {
      InterpolateHN(2,EnKM,C0KM,C1KM,B0KM,B1KM);
    }
    if(hLabMomentum < 1) { Coeff2 = 0.01; }
    else  { Coeff2 = 0.01/(hLabMomentum2*hLabMomentum); }
    break; 
  }
 
  if(verboseLevel > 2) {
    G4cout<<"  HadrVal : Plasb  "<<hLabMomentum
      <<"  iHadron  "<<iHadron<<"  HadrTot  "<<HadrTot<<G4endl;
  }
}

References BoundaryP, Coeff0, Coeff1, Coeff2, DDSect2, DDSect3, G4cout, G4endl, G4Exp(), G4Log(), HadrEnergy, HadrReIm, HadrSlope, HadrTot, hLabMomentum, hLabMomentum2, hMass, hMass2, iHadrCode, iHadron, InterpolateHN(), protonM, protonM2, Slope0, Slope1, Slope2, TotP, G4HadronicInteraction::verboseLevel, and Z.

Referenced by FillData(), and HadronProtonQ2().

FillData()

void G4ElasticHadrNucleusHE::FillData ( const G4ParticleDefinition *  p, G4int  idx, G4int  Z  )

private

Definition at line 441 of file G4ElasticHadrNucleusHE.cc.

{
#ifdef G4MULTITHREADED
  G4MUTEXLOCK(&elasticMutex);
  if(!fElasticData[idx][Z]) {
#endif
    G4int A = G4lrint(nistManager->GetAtomicMassAmu(Z));
    G4ElasticData* pElD = new G4ElasticData(p, Z, A, fEnergy);
    if(fRetrieveFromFile) { 
      std::ostringstream ss;
      InFileName(ss, p, Z); 
      std::ifstream infile(ss.str(), std::ios::in);
      for(G4int i=0; i<NENERGY; ++i) {
    if(ReadLine(infile, pElD->fCumProb[i])) {
      continue;
    } else {
      fRetrieveFromFile = false;
          break;
    }
      }
      infile.close();
    }
    R1     = pElD->R1;
    R2     = pElD->R2;
    Aeff   = pElD->Aeff;
    Pnucl  = pElD->Pnucl;
    dQ2    = pElD->dQ2;
    if(verboseLevel > 0) {
      G4cout<<"### FillData for " << p->GetParticleName() 
        << " Z= " << Z << " idx= " << idx << " iHadron= " << iHadron 
        <<" iHadron1= " << iHadron1 << " iHadrCode= " << iHadrCode
            <<"\n   R1= " << R1 << " R2= " << R2 << " Aeff= " << Aeff 
        <<" Pnucl= " << Pnucl << G4endl;
    }
 
    if(!fRetrieveFromFile) {  
      for(G4int i=0; i<NENERGY; ++i) {
    G4double T = fEnergy[i];
    hLabMomentum2 = T*(T + 2.*hMass);
    hLabMomentum  = std::sqrt(hLabMomentum2);
    HadrEnergy = hMass + T;
    DefineHadronValues(Z);
    Q2max = pElD->maxQ2[i];
 
    G4int length  = FillFq2(A); 
    (pElD->fCumProb[i]).reserve(length);
    G4double norm = 1.0/fLineF[length-1];
 
    if(verboseLevel > 0) {
      G4cout << "### i= " << i << " Z= " << Z << " A= " << A 
         << " length= " << length << " Q2max= " << Q2max << G4endl;
    }
 
    (pElD->fCumProb[i]).push_back(0.0);
    for(G4int ii=1; ii<length-1; ++ii) {
      (pElD->fCumProb[i]).push_back(fLineF[ii]*norm);
      if(verboseLevel > 2) {
        G4cout << "    ii= " << ii << " val= " 
           << (pElD->fCumProb[i])[ii] << G4endl;
      }
    }
    (pElD->fCumProb[i]).push_back(1.0);
      }
    }
 
    if(fStoreToFile) {
      std::ostringstream ss;
      OutFileName(ss, p, Z); 
      std::ofstream fileout(ss.str());
      for(G4int i=0; i<NENERGY; ++i) {
    WriteLine(fileout, pElD->fCumProb[i]);
      }
      fileout.close();
    }
    
    if(verboseLevel > 0) {
      G4cout << " G4ElasticHadrNucleusHE::FillData done for idx= " << idx
         << " for " << p->GetParticleName() << " Z= " << Z
         << " A= " << A << G4endl;
    }
    fElasticData[idx][Z] = pElD;
 
#ifdef G4MULTITHREADED
  }
  G4MUTEXUNLOCK(&elasticMutex);
#endif  
}

References A, G4ElasticData::Aeff, Aeff, DefineHadronValues(), G4ElasticData::dQ2, dQ2, G4ElasticData::fCumProb, fElasticData, fEnergy, FillFq2(), fLineF, fRetrieveFromFile, fStoreToFile, G4cout, G4endl, G4lrint(), G4MUTEXLOCK, G4MUTEXUNLOCK, G4NistManager::GetAtomicMassAmu(), G4ParticleDefinition::GetParticleName(), HadrEnergy, hLabMomentum, hLabMomentum2, hMass, iHadrCode, iHadron, iHadron1, InFileName(), G4ElasticData::maxQ2, NENERGY, nistManager, OutFileName(), G4ElasticData::Pnucl, Pnucl, Q2max, G4ElasticData::R1, R1, G4ElasticData::R2, R2, ReadLine(), G4HadronicInteraction::verboseLevel, WriteLine(), and Z.

Referenced by InitialiseModel(), and SampleInvariantT().

FillFq2()

G4int G4ElasticHadrNucleusHE::FillFq2 ( G4int  A )

private

Definition at line 661 of file G4ElasticHadrNucleusHE.cc.

{
  G4double curQ2, curSec;
  G4double curSum = 0.0;
  G4double totSum = 0.0;
 
  G4double ddQ2 = dQ2*0.1;
  G4double Q2l  = 0.0;
 
  G4int ii = 0;
  for(ii=1; ii<ONQ2-1; ++ii) {
    curSum = curSec = 0.0;
 
    for(G4int jj=0; jj<10; ++jj) {
      curQ2 = Q2l+(jj + 0.5)*ddQ2;
      if(curQ2 >= Q2max) { break; }
      curSec = HadrNucDifferCrSec(A, curQ2);
      curSum += curSec;
    }
    G4double del = (curQ2 >= Q2max) ? Q2max - Q2l : dQ2;
    Q2l    += del;
    curSum *= del*0.1;
    totSum += curSum;
    fLineF[ii] = totSum;
    if (verboseLevel>2) {
      G4cout<<ii << ". FillFq2: A= " << A << " Q2= "<<Q2l<<" dQ2= "
        <<dQ2<<" Tot= "<<totSum << " dTot " <<curSum
        <<" curSec= " <<curSec<<G4endl;
    }
    if(totSum*1.e-4 > curSum || Q2l >= Q2max) { break; }
  }
  ii = std::min(ii, ONQ2-2);
  curQ2 = Q2l;
  G4double xx = R1*(Q2max - curQ2);
  if(xx > 0.0) {
    xx = (xx > 20.) ? 0.0 : G4Exp(-xx);
    curSec = HadrNucDifferCrSec(A, curQ2);
    totSum += curSec*(1.0 - xx)/R1;
  }
  fLineF[ii + 1] = totSum;
  if (verboseLevel>1) {
    G4cout << "### FillFq2 done curQ2= " << curQ2 << " Q2max= "<< Q2max 
           << " sumG= " << fLineF[ONQ2-2] << "  totSum= " << totSum 
       << " Nbins= " << ii + 1 << G4endl;
  }
  return ii + 2;
}

References A, dQ2, fLineF, G4cout, G4endl, G4Exp(), HadrNucDifferCrSec(), G4INCL::Math::min(), ONQ2, Q2max, R1, and G4HadronicInteraction::verboseLevel.

Referenced by FillData().

GetBinomCof()

G4double G4ElasticHadrNucleusHE::GetBinomCof ( G4int  n, G4int  m  )

inlineprivate

Definition at line 221 of file G4ElasticHadrNucleusHE.hh.

{
  return (numN >= numM && numN < 240) ? fBinom[numN][numM] : 0.0;
}

References fBinom.

Referenced by GetLightFq2().

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().

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().

GetFt()

G4double G4ElasticHadrNucleusHE::GetFt ( G4double  Q2 )

private

Definition at line 1284 of file G4ElasticHadrNucleusHE.cc.

{
  G4double Fdistr=0;
  G4double SqrQ2 = std::sqrt(Q2);
 
  Fdistr = (1-Coeff1-Coeff0) / HadrSlope*(1-G4Exp(-HadrSlope*Q2))
    + Coeff0*(1-G4Exp(-Slope0*Q2))
    + Coeff2/Slope2*G4Exp(Slope2*ConstU)*(G4Exp(Slope2*Q2)-1)
    + 2*Coeff1/Slope1*(1/Slope1-(1/Slope1+SqrQ2)*G4Exp(-Slope1*SqrQ2));
 
  if (verboseLevel>1) {
    G4cout<<"Old:  Coeff0 Coeff1 Coeff2 "<<Coeff0<<"  "
          <<Coeff1<<"  "<<Coeff2<<"  Slope Slope0 Slope1 Slope2 "
          <<HadrSlope<<"  "<<Slope0<<"  "<<Slope1<<"  "<<Slope2
          <<"  Fdistr "<<Fdistr<<G4endl;
  }
  return Fdistr;
}

References Coeff0, Coeff1, Coeff2, ConstU, G4cout, G4endl, G4Exp(), HadrSlope, Slope0, Slope1, Slope2, and G4HadronicInteraction::verboseLevel.

Referenced by HadronProtonQ2().

GetLightFq2()

G4double G4ElasticHadrNucleusHE::GetLightFq2 ( G4int  Z, G4int  A, G4double  Q  )

private

Definition at line 711 of file G4ElasticHadrNucleusHE.cc.

{
  // Scattering off proton
  if(Z == 1) 
  {
    G4double SqrQ2  = std::sqrt(Q2);
    G4double valueConstU = 2.*(hMass2 + protonM2) - Q2;
 
    G4double y = (1.-Coeff1-Coeff0)/HadrSlope*(1.-G4Exp(-HadrSlope*Q2))
      + Coeff0*(1.-G4Exp(-Slope0*Q2))
      + Coeff2/Slope2*G4Exp(Slope2*valueConstU)*(G4Exp(Slope2*Q2)-1.)
      + 2.*Coeff1/Slope1*(1./Slope1-(1./Slope1+SqrQ2)*G4Exp(-Slope1*SqrQ2));
 
    return y;
  }
 
  // The preparing of probability function  
 
  G4double prec = A > 208  ?  1.0e-7 : 1.0e-6;
 
  G4double    Stot     = HadrTot*MbToGeV2;     //  Gev^-2
  G4double    Bhad     = HadrSlope;         //  GeV^-2
  G4double    Asq      = 1+HadrReIm*HadrReIm;
  G4double    Rho2     = std::sqrt(Asq);
 
  if(verboseLevel >1) {
    G4cout<<" Fq2 Before for i Tot B Im "<<HadrTot<<"  "<<HadrSlope<<"  "
      <<HadrReIm<<G4endl;
  }
  if(verboseLevel > 1) {
    G4cout << "GetFq2: Stot= " << Stot << " Bhad= " << Bhad 
           <<"  Im "<<HadrReIm 
           << " Asq= " << Asq << G4endl;
    G4cout << "R1= " << R1 << " R2= " << R2 << " Pnucl= " << Pnucl <<G4endl;
  }
  G4double    R12      = R1*R1;
  G4double    R22      = R2*R2;
  G4double    R12B     = R12+2*Bhad;
  G4double    R22B     = R22+2*Bhad;
 
  G4double    Norm     = (R12*R1-Pnucl*R22*R2); // HP->Aeff;
 
  G4double    R13      = R12*R1/R12B;
  G4double    R23      = Pnucl*R22*R2/R22B;
  G4double    Unucl    = Stot/twopi*R13/Norm;
  G4double    UnucRho2 = -Unucl*Rho2;
 
  G4double    FiH      = std::asin(HadrReIm/Rho2);
  G4double    NN2      = R23/R13;
 
  if(verboseLevel > 2) {
    G4cout << "UnucRho2= " << UnucRho2 << " FiH= " << FiH << " NN2= " << NN2 
       << " Norm= " << Norm << G4endl;
  }
  G4double    Prod0 = 0.;
  G4double    N1    = -1.0;
 
  for(G4int i1 = 1; i1<= A; ++i1) 
    {
      N1 *= (-Unucl*Rho2*(A-i1+1)/(G4double)i1);
      G4double Prod1 = 0.;
      G4double N2    = -1.;
 
      for(G4int i2 = 1; i2<=A; ++i2) 
        {
          N2 *= (-Unucl*Rho2*(A-i2+1)/(G4double)i2);
          G4double Prod2 = 0; 
          G4double N5    = -1/NN2;
      for(G4int j2=0; j2<= i2; ++j2) 
            {
              G4double Prod3 = 0;
              G4double exp2  = 1./((G4double)j2/R22B+(G4double)(i2-j2)/R12B);
              N5 *= (-NN2);
              G4double N4 = -1./NN2;
          for(G4int j1=0; j1<=i1; ++j1) 
        {
          G4double exp1  = 1./((G4double)j1/R22B+(G4double)(i1-j1)/R12B);
          G4double dddd  = 0.25*(exp1+exp2);
          N4    *= (-NN2);
          Prod3 += 
                    N4*exp1*exp2*(1.-G4Exp(-Q2*dddd))*GetBinomCof(i1,j1)/dddd;
        }                                   // j1
          Prod2 += Prod3*N5*GetBinomCof(i2,j2);
        }                                      // j2
      Prod1 += Prod2*N2*std::cos(FiH*(i1-i2));
 
      if (std::abs(Prod2*N2/Prod1)<prec) break;
        }                                         // i2
      Prod0 += Prod1*N1;
      if(std::abs(N1*Prod1/Prod0) < prec) break;
    }                                           // i1
 
  const G4double fact = 0.25*CLHEP::pi/MbToGeV2; 
  Prod0 *= fact;  //  This is in mb
 
  if(verboseLevel>1) {
    G4cout << "GetLightFq2 Z= " << Z << " A= " << A 
       <<" Q2= " << Q2 << " Res= " << Prod0 << G4endl;
  }
  return Prod0;
}

References A, Coeff0, Coeff1, Coeff2, G4cout, G4endl, G4Exp(), GetBinomCof(), HadrReIm, HadrSlope, HadrTot, hMass2, MbToGeV2, CLHEP::pi, Pnucl, CLHEP::prec, protonM2, R1, R2, Slope0, Slope1, Slope2, twopi, G4HadronicInteraction::verboseLevel, and Z.

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::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().

GetQ2_2()

G4double G4ElasticHadrNucleusHE::GetQ2_2 ( G4int  N, G4int  Nmax, const std::vector< G4double > &  F, G4double  rand  )

private

Definition at line 599 of file G4ElasticHadrNucleusHE.cc.

{
  //G4cout << "GetQ2_2 kk= " << kk << " kmax= " << kmax << "  size= " 
  //     << F.size() << "  rand= " << ranUni << G4endl;
  if(kk == kmax-1) {
    G4double X1 = dQ2*kk;
    G4double F1 = F[kk-1];
    G4double X2 = Q2max;
    G4double xx = R1*(X2 - X1);
    xx = (xx > 20.) ? 0.0 : G4Exp(-xx);
    G4double Y = X1 - G4Log(1.0 - (ranUni - F1)*(1.0 - xx)/(1.0 - F1))/R1;
    return Y;
  } 
  G4double F1, F2, F3, X1, X2, X3;
 
  if(kk == 1 || kk == 0) {
    F1 = F[0]; 
    F2 = F[1];
    F3 = F[2];
    X1 = 0.0;
    X2 = dQ2;
    X3 = dQ2*2;
  } else {
    F1 = F[kk-2];
    F2 = F[kk-1];
    F3 = F[kk];
    X1 = dQ2*(kk-2);
    X2 = dQ2*(kk-1);
    X3 = dQ2*kk;
  }
  if(verboseLevel > 1) {
    G4cout << "GetQ2_2 kk= " << kk << " X2= " << X2 << " X3= " << X3 
       << " F2= " << F2 << " F3= " << F3 << " Rndm= " << ranUni << G4endl;
  }
 
  G4double F12 = F1*F1;
  G4double F22 = F2*F2;
  G4double F32 = F3*F3;
 
  G4double D0  = F12*F2+F1*F32+F3*F22-F32*F2-F22*F1-F12*F3;
 
  if(verboseLevel > 2) {
    G4cout << "       X1= " << X1 << " F1= " << F1 << "  D0= " 
           << D0 << G4endl; 
  }
  G4double Y;
  if(std::abs(D0) < 1.e-9) { 
    Y = X2 + (ranUni - F2)*(X3 - X2)/(F3 - F2);
  } else {
    G4double DA = X1*F2+X3*F1+X2*F3-X3*F2-X1*F3-X2*F1;
    G4double DB = X2*F12+X1*F32+X3*F22-X2*F32-X3*F12-X1*F22;
    G4double DC = X3*F2*F12+X2*F1*F32+X1*F3*F22
               -X1*F2*F32-X2*F3*F12-X3*F1*F22;
    Y = (DA*ranUni*ranUni + DB*ranUni + DC)/D0;
  }
  return Y;
}

References dQ2, F12, F22, F32, G4cout, G4endl, G4Exp(), G4Log(), Q2max, R1, G4HadronicInteraction::verboseLevel, and Y().

Referenced by HadronNucleusQ2_2().

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().

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.

HadrNucDifferCrSec()

G4double G4ElasticHadrNucleusHE::HadrNucDifferCrSec ( G4int  A, G4double  Q2  )

private

Definition at line 816 of file G4ElasticHadrNucleusHE.cc.

{
  //   ------ All external kinematical variables are in MeV -------
  //            ------ but internal in GeV !!!!  ------
 
  // Scattering of proton
  if(A == 1) 
  {
    G4double SqrQ2  = std::sqrt(aQ2);
    G4double valueConstU = hMass2 + protonM2-2*protonM*HadrEnergy - aQ2;
    
    BoundaryTL[0] = Q2max;
    BoundaryTL[1] = Q2max;
    BoundaryTL[3] = Q2max;
    BoundaryTL[4] = Q2max;
    BoundaryTL[5] = Q2max;
 
    G4double dSigPodT = HadrTot*HadrTot*(1+HadrReIm*HadrReIm)*
      ( Coeff1*G4Exp(-Slope1*SqrQ2)+
    Coeff2*G4Exp( Slope2*(valueConstU)+aQ2)+
    (1-Coeff1-Coeff0)*G4Exp(-HadrSlope*aQ2)+
    Coeff0*G4Exp(-Slope0*aQ2) )*2.568/(16*pi);
 
    return dSigPodT;
  }
 
  G4double    Stot     = HadrTot*MbToGeV2; 
  G4double    Bhad     = HadrSlope; 
  G4double    Asq      = 1+HadrReIm*HadrReIm;
  G4double    Rho2     = std::sqrt(Asq);
  G4double    R12      = R1*R1;
  G4double    R22      = R2*R2;
  G4double    R12B     = R12+2*Bhad;
  G4double    R22B     = R22+2*Bhad;
  G4double    R12Ap    = R12+20;
  G4double    R22Ap    = R22+20;
  G4double    R13Ap    = R12*R1/R12Ap;
  G4double    R23Ap    = R22*R2*Pnucl/R22Ap;
  G4double    R23dR13  = R23Ap/R13Ap;
  G4double    R12Apd   = 2/R12Ap;
  G4double    R22Apd   = 2/R22Ap;
  G4double R12ApdR22Ap = 0.5*(R12Apd+R22Apd);
 
  G4double DDSec1p  = (DDSect2+DDSect3*G4Log(0.53*HadrEnergy/R1));
  G4double DDSec2p  = (DDSect2+DDSect3*G4Log(0.53*HadrEnergy/
                             std::sqrt((R12+R22)*0.5)));
  G4double DDSec3p  = (DDSect2+DDSect3*G4Log(0.53*HadrEnergy/R2));
 
  G4double    Norm     = (R12*R1-Pnucl*R22*R2)*Aeff;
  G4double    R13      = R12*R1/R12B;
  G4double    R23      = Pnucl*R22*R2/R22B;
  G4double    Unucl    = Stot/(twopi*Norm)*R13;
  G4double    UnuclScr = Stot/(twopi*Norm)*R13Ap;
  G4double    SinFi    = HadrReIm/Rho2;
  G4double    FiH      = std::asin(SinFi);
  G4double    N        = -1;
  G4double    N2       = R23/R13;
 
  G4double ImElasticAmpl0 = 0;
  G4double ReElasticAmpl0 = 0;
  G4double exp1;
 
  for(G4int i=1; i<=A; ++i) {
    N  *= (-Unucl*Rho2*(A-i+1)/(G4double)i);
    G4double N4 = 1;
    G4double medTot = R12B/(G4double)i;
    G4double Prod1  = G4Exp(-aQ2*R12B/(G4double)(4*i))*medTot;
 
    for(G4int l=1; l<=i; ++l) {
      exp1 = l/R22B+(i-l)/R12B;
      N4 *= (-N2*(i-l+1)/(G4double)l);
      G4double expn4 = N4/exp1;
      Prod1  += expn4*G4Exp(-aQ2/(exp1*4));
      medTot += expn4;
    }  // end l
 
    G4double dcos = N*std::cos(FiH*i);
    ReElasticAmpl0  += Prod1*N*std::sin(FiH*i);
    ImElasticAmpl0  += Prod1*dcos;
    if(std::abs(Prod1*N/ImElasticAmpl0) < 0.000001) break;
  }      // i
 
  static const G4double pi25 = CLHEP::pi/2.568;
  ImElasticAmpl0 *= pi25;   // The amplitude in mB
  ReElasticAmpl0 *= pi25;   // The amplitude in mB
 
  G4double C1 = R13Ap*R13Ap*0.5*DDSec1p;
  G4double C2 = 2*R23Ap*R13Ap*0.5*DDSec2p;
  G4double C3 = R23Ap*R23Ap*0.5*DDSec3p;
  
  G4double N1p  = 1;
  G4double Din1 = 0.5*(C1*G4Exp(-aQ2/8*R12Ap)/2*R12Ap-
               C2/R12ApdR22Ap*G4Exp(-aQ2/(4*R12ApdR22Ap))+
               C3*R22Ap/2*G4Exp(-aQ2/8*R22Ap));
 
  G4double DTot1 = 0.5*(C1*0.5*R12Ap-C2/R12ApdR22Ap+C3*R22Ap*0.5);
 
  for(G4int i=1; i<= A-2; ++i) {
    N1p *= (-UnuclScr*Rho2*(A-i-1)/(G4double)i);
    G4double N2p  = 1;
    G4double Din2 = 0;
    G4double DmedTot = 0;
    G4double BinCoeff = 1.0;
    for(G4int l=0; l<=i; ++l) {
      if(l > 0) { BinCoeff *= (i-l+1)/(G4double)l; }
 
      exp1  = l/R22B+(i-l)/R12B;
      G4double exp1p = exp1+R12Apd;
      G4double exp2p = exp1+R12ApdR22Ap;
      G4double exp3p = exp1+R22Apd;
 
      Din2 += N2p*BinCoeff*(C1/exp1p*G4Exp(-aQ2/(4*exp1p))-
                C2/exp2p*G4Exp(-aQ2/(4*exp2p))+
                C3/exp3p*G4Exp(-aQ2/(4*exp3p)));
 
      DmedTot += N2p*BinCoeff*(C1/exp1p-C2/exp2p+C3/exp3p);
 
      N2p *= -R23dR13;
    }     // l
 
    G4double dcos = N1p*std::cos(FiH*i)/(G4double)((i+2)*(i+1));
    Din1  += Din2*dcos;
    DTot1 += DmedTot*dcos;
    
    if(std::abs(Din2*N1p/Din1) < 0.000001) break;
  }           //  i
  G4double gg = (G4double)(A*(A-1)*4)/(Norm*Norm); 
 
  Din1  *= (-gg);
  DTot1 *= 5*gg;
 
  //  ----------------  dSigma/d|-t|,  mb/(GeV/c)^-2  -----------------
 
  G4double DiffCrSec2 = (ReElasticAmpl0*ReElasticAmpl0+
             (ImElasticAmpl0+Din1)*
             (ImElasticAmpl0+Din1))/twopi;
 
  Dtot11 = DTot1;
  aAIm   = ImElasticAmpl0;
  aDIm   = Din1;
 
  return DiffCrSec2;  //  dSig/d|-t|,  mb/(GeV/c)^-2
} 

References A, aAIm, aDIm, Aeff, BoundaryTL, C1, C2, C3, Coeff0, Coeff1, Coeff2, DDSect2, DDSect3, Dtot11, G4Exp(), G4Log(), HadrEnergy, HadrReIm, HadrSlope, HadrTot, hMass2, MbToGeV2, CLHEP::pi, pi, Pnucl, protonM, protonM2, Q2max, R1, R2, Slope0, Slope1, Slope2, and twopi.

Referenced by FillFq2().

HadronNucleusQ2_2()

G4double G4ElasticHadrNucleusHE::HadronNucleusQ2_2 ( const G4ElasticData *  pElD, G4double  plabGeV, G4double  tmax  )

private

Definition at line 553 of file G4ElasticHadrNucleusHE.cc.

{
  G4double ekin  = std::sqrt(hMass2 + plab*plab) - hMass;
 
  if(verboseLevel > 1) {
    G4cout<<"Q2_2: ekin(GeV)= " << ekin << "  plab(GeV/c)= " << plab 
          <<"  tmax(GeV2)= " << tmax <<G4endl;
  }
  // Find closest energy bin
  G4int idx; 
  for(idx=0; idx<NENERGY-1; ++idx) {
    if(ekin <= fLowEdgeEnergy[idx+1]) { break; }
  }
  //G4cout << "   idx= " << idx << G4endl;
 
  // Select kinematics for node energy
  R1    = pElD->R1;
  dQ2   = pElD->dQ2;
  Q2max = pElD->maxQ2[idx];
  G4int length = (pElD->fCumProb[idx]).size();
 
  G4double Rand = G4UniformRand();
 
  G4int iNumbQ2 = 0;
  for(iNumbQ2=1; iNumbQ2<length; ++iNumbQ2) {
    if(Rand <= (pElD->fCumProb[idx])[iNumbQ2]) { break; }
  }
  iNumbQ2 = std::min(iNumbQ2, length - 1);
  G4double Q2 = GetQ2_2(iNumbQ2, length, pElD->fCumProb[idx], Rand);
  Q2 = std::min(Q2, Q2max);
  Q2 *= tmax/Q2max;
 
  if(verboseLevel > 1) {
    G4cout<<" HadrNucleusQ2_2(2): Q2= "<<Q2<<" iNumbQ2= " << iNumbQ2 
      << " rand= " << Rand << " Q2max= " << Q2max 
          << " tmax= " << tmax << G4endl;
  }
  return Q2;
}       

References G4ElasticData::dQ2, dQ2, G4ElasticData::fCumProb, fLowEdgeEnergy, G4cout, G4endl, G4UniformRand, GetQ2_2(), hMass, hMass2, G4ElasticData::maxQ2, G4INCL::Math::min(), NENERGY, Q2max, G4ElasticData::R1, R1, and G4HadronicInteraction::verboseLevel.

Referenced by SampleInvariantT().

HadronProtonQ2()

G4double G4ElasticHadrNucleusHE::HadronProtonQ2 ( G4double  plab, G4double  tmax  )

private

Definition at line 1306 of file G4ElasticHadrNucleusHE.cc.

{
  hLabMomentum  = plab;
  hLabMomentum2 = hLabMomentum*hLabMomentum;
  HadrEnergy    = std::sqrt(hMass2 + hLabMomentum2);
  DefineHadronValues(1);
 
  G4double Sh = 2.0*protonM*HadrEnergy+protonM2+hMass2; // GeV
  ConstU = 2*protonM2+2*hMass2-Sh;
 
  BoundaryTL[0] = tmax;
  BoundaryTL[1] = tmax;
  BoundaryTL[3] = tmax;
  BoundaryTL[4] = tmax;
  BoundaryTL[5] = tmax;
 
  G4double MaxTR = (plab < BoundaryP[iHadron1]) ? 
    BoundaryTL[iHadron1] : BoundaryTG[iHadron1]; 
 
  if (verboseLevel>1) {
    G4cout<<"3  GetKin. : iHadron1  "<<iHadron1
      <<"  Bound.P[iHadron1] "<<BoundaryP[iHadron1]
      <<"  Bound.TL[iHadron1] "<<BoundaryTL[iHadron1]
      <<"  Bound.TG[iHadron1] "<<BoundaryTG[iHadron1]
      <<"  MaxT MaxTR "<<tmax<<"  "<<MaxTR<<G4endl;
  }
 
  G4double rand = G4UniformRand();
 
  G4double DDD0=MaxTR*0.5, DDD1=0.0, DDD2=MaxTR;
 
  G4double norm  = 1.0/GetFt(MaxTR);
  G4double delta = GetFt(DDD0)*norm - rand;
 
  static const G4int maxNumberOfLoops = 10000;
  G4int loopCounter = -1;
  while ( (std::abs(delta) > 0.0001) && 
          ++loopCounter < maxNumberOfLoops )  /* Loop checking, 10.08.2015, A.Ribon */
    {
      if(delta>0) 
      {
        DDD2 = DDD0;
        DDD0 = (DDD0+DDD1)*0.5;
      }
      else if(delta<0.0)
      {
        DDD1 = DDD0; 
        DDD0 = (DDD0+DDD2)*0.5;
      }
      delta = GetFt(DDD0)*norm - rand;
    }
  return (loopCounter >= maxNumberOfLoops) ? 0.0 : DDD0;
}

References BoundaryP, BoundaryTG, BoundaryTL, ConstU, DefineHadronValues(), G4cout, G4endl, G4UniformRand, GetFt(), HadrEnergy, hLabMomentum, hLabMomentum2, hMass2, iHadron1, protonM, protonM2, and G4HadronicInteraction::verboseLevel.

Referenced by SampleInvariantT().

InFileName()

void G4ElasticHadrNucleusHE::InFileName ( std::ostringstream &  ss, const G4ParticleDefinition *  p, G4int  Z  )

private

Definition at line 1380 of file G4ElasticHadrNucleusHE.cc.

{
  if(!fDirectory) {
    fDirectory = std::getenv("G4LEDATA");
    if (fDirectory) { 
      ss << fDirectory << "/";
    }
  }
  OutFileName(ss, p, Z);
}

References fDirectory, OutFileName(), and Z.

Referenced by FillData().

InitialiseModel()

void G4ElasticHadrNucleusHE::InitialiseModel ( )

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 337 of file G4ElasticHadrNucleusHE.cc.

{
  if(!isMaster) { return; }
  G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();
  
  for(G4int i=0; i<2; ++i) {
    const G4ParticleDefinition* p = G4PionPlus::PionPlus();
    if(1 == i) { p = G4PionMinus::PionMinus(); } 
    iHadrCode = fHadronCode[i]; 
    iHadron   = fHadronType[i];
    iHadron1  = fHadronType1[i];
    hMass     = p->GetPDGMass()*invGeV;
    hMass2    = hMass*hMass;
    for(size_t j=0; j<numOfCouples; ++j) {
      auto mat = theCoupleTable->GetMaterialCutsCouple(j)->GetMaterial();
      auto elmVec = mat->GetElementVector();
      size_t numOfElem = mat->GetNumberOfElements();
      for(size_t k=0; k<numOfElem; ++k) {
    G4int Z = std::min((*elmVec)[k]->GetZasInt(), ZMAX-1);
    if(!fElasticData[i][Z]) { 
          if(1 == i && Z > 1) { 
        fElasticData[1][Z] = fElasticData[0][Z]; 
      } else {
        FillData(p, i, Z);
      } 
    }
      }
    }
  }
}

References fElasticData, fHadronCode, fHadronType, fHadronType1, FillData(), G4Material::GetElementVector(), G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4ParticleDefinition::GetPDGMass(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), hMass, hMass2, iHadrCode, iHadron, iHadron1, invGeV, isMaster, G4INCL::Math::min(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), Z, and ZMAX.

InterpolateHN()

void G4ElasticHadrNucleusHE::InterpolateHN ( G4int  n, const G4double  EnP[], const G4double  C0P[], const G4double  C1P[], const G4double  B0P[], const G4double  B1P[]  )

private

Definition at line 532 of file G4ElasticHadrNucleusHE.cc.

{
  G4int i; 
 
  for(i=1; i<n; ++i) { if(hLabMomentum <= EnP[i]) { break; } }
  if(i == n) { i = n - 1; }
 
  Coeff0 = LineInterpol(EnP[i], EnP[i-1], C0P[i], C0P[i-1], hLabMomentum);
  Coeff1 = LineInterpol(EnP[i], EnP[i-1], C1P[i], C1P[i-1], hLabMomentum);
  Slope0 = LineInterpol(EnP[i], EnP[i-1], B0P[i], B0P[i-1], hLabMomentum);
  Slope1 = LineInterpol(EnP[i], EnP[i-1], B1P[i], B1P[i-1], hLabMomentum);
 
//  G4cout<<"  InterpolHN:  n i "<<n<<"  "<<i<<"  Mom "
//        <<hLabMomentum<<G4endl;
}

References Coeff0, Coeff1, hLabMomentum, LineInterpol(), CLHEP::detail::n, Slope0, and Slope1.

Referenced by DefineHadronValues().

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.

LineInterpol()

G4double G4ElasticHadrNucleusHE::LineInterpol ( G4double  p0, G4double  p2, G4double  c1, G4double  c2, G4double  p  )

inlineprivate

Definition at line 211 of file G4ElasticHadrNucleusHE.hh.

{
  return c1+(p-p1)*(c2-c1)/(p2-p1);
}

Referenced by InterpolateHN().

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 G4ElasticHadrNucleusHE::ModelDescription ( std::ostream &  outFile ) const

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 305 of file G4ElasticHadrNucleusHE.cc.

{
  outFile << "G4ElasticHadrNucleusHE is a hadron-nucleus elastic scattering\n"
      << "model developed by N. Starkov which uses a Glauber model\n"
      << "parameterization to calculate the final state.  It is valid\n"
      << "for all hadrons with incident momentum above 0.4 GeV/c.\n";
}

operator!=()

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

deleteinherited

operator=()

G4ElasticHadrNucleusHE & G4ElasticHadrNucleusHE::operator= ( const G4ElasticHadrNucleusHE &  right )

private

operator==()

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

deleteinherited

OutFileName()

void G4ElasticHadrNucleusHE::OutFileName ( std::ostringstream &  ss, const G4ParticleDefinition *  p, G4int  Z  )

private

Definition at line 1395 of file G4ElasticHadrNucleusHE.cc.

{
  ss << "hedata/" << p->GetParticleName() << Z << ".dat";
}

References G4ParticleDefinition::GetParticleName(), and Z.

Referenced by FillData(), and InFileName().

ReadLine()

G4bool G4ElasticHadrNucleusHE::ReadLine ( std::ifstream &  infile, std::vector< G4double > &  v  )

private

Definition at line 1403 of file G4ElasticHadrNucleusHE.cc.

{
  G4int n(0);
  infile >> n;
  if (infile.fail()) { return false; }
  if(n > 0) {
    v.reserve(n);
    G4double x(0.0);
    for(G4int i=0; i<n; ++i) {
      infile >> x;
      if (infile.fail()) { return false; }
      v.emplace_back(x);
    }
  }
  return true;
}

References CLHEP::detail::n.

Referenced by FillData().

SampleInvariantT()

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

overridevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 373 of file G4ElasticHadrNucleusHE.cc.

{
  G4double mass = p->GetPDGMass();
  G4double kine = sqrt(inLabMom*inLabMom + mass*mass) - mass;
  if(kine <= ekinLowLimit) {
    return G4HadronElastic::SampleInvariantT(p,inLabMom,iZ,A);
  }
  G4int Z = std::min(iZ,ZMAX-1);
  G4double Q2 = 0.0;
  iHadrCode = p->GetPDGEncoding();
 
  // below computations in GeV/c
  hMass  = mass*invGeV;
  hMass2 = hMass*hMass;
  G4double plab = inLabMom*invGeV;
  G4double tmax = pLocalTmax*invGeV2;
 
  if(verboseLevel > 1) {
    G4cout<< "G4ElasticHadrNucleusHE::SampleT: " 
      << " for " << p->GetParticleName() 
      << " at Z= " << Z << " A= " << A
      << " plab(GeV)= " << plab
      << " hadrCode= " << iHadrCode
      << G4endl;
  }
  iHadron = -1;
  G4int idx;
  for(idx=0; idx<NHADRONS; ++idx) {
    if(iHadrCode == fHadronCode[idx]) { 
      iHadron = fHadronType[idx];
      iHadron1 = fHadronType1[idx];
      break; 
    }
  }
  // Hadron is not in the list
  if(0 > iHadron) { return 0.0; }
 
  if(Z==1) {
    Q2 = HadronProtonQ2(plab, tmax);
 
    if (verboseLevel>1) {
      G4cout<<"  Proton : Q2  "<<Q2<<G4endl;
    }
  } else {
    const G4ElasticData* ElD1 = fElasticData[idx][Z];
 
    // Construct elastic data
    if(!ElD1) { 
      FillData(p, idx, Z); 
      ElD1 = fElasticData[idx][Z];
      if(!ElD1) { return 0.0; }
    }
 
    // sample scattering
    Q2 = HadronNucleusQ2_2(ElD1, plab, tmax);
 
    if(verboseLevel > 1) {
      G4cout<<" SampleT: Q2(GeV^2)= "<<Q2<< "  t/tmax= " 
            << Q2/tmax <<G4endl;
    }
  }
  return Q2*GeV2;
}

References A, ekinLowLimit, fElasticData, fHadronCode, fHadronType, fHadronType1, FillData(), G4cout, G4endl, G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGEncoding(), G4ParticleDefinition::GetPDGMass(), GeV2, HadronNucleusQ2_2(), HadronProtonQ2(), hMass, hMass2, iHadrCode, iHadron, iHadron1, invGeV, invGeV2, G4INCL::Math::min(), NHADRONS, G4HadronElastic::pLocalTmax, G4HadronElastic::SampleInvariantT(), G4HadronicInteraction::verboseLevel, Z, and ZMAX.

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().

SetLowestEnergyLimit()

void G4HadronElastic::SetLowestEnergyLimit ( G4double  value )

inlineinherited

Definition at line 91 of file G4HadronElastic.hh.

{
  lowestEnergyLimit = value;
}

References G4HadronElastic::lowestEnergyLimit.

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

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::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.

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::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().

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().

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().

WriteLine()

void G4ElasticHadrNucleusHE::WriteLine ( std::ofstream &  outfile, std::vector< G4double > &  v  )

private

Definition at line 1423 of file G4ElasticHadrNucleusHE.cc.

{
  G4int n = v.size();
  outfile << n << G4endl;
  if(n > 0) {
    for(G4int i=0; i<n; ++i) {
      outfile << v[i] << " ";
    }
    outfile << G4endl;
  }
}

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

Referenced by FillData().

Field Documentation

aAIm

G4double G4ElasticHadrNucleusHE::aAIm

private

Definition at line 186 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), and HadrNucDifferCrSec().

aDIm

G4double G4ElasticHadrNucleusHE::aDIm

private

Definition at line 186 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), and HadrNucDifferCrSec().

Aeff

G4double G4ElasticHadrNucleusHE::Aeff

private

Definition at line 189 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), G4ElasticHadrNucleusHE(), and HadrNucDifferCrSec().

BoundaryP

G4double G4ElasticHadrNucleusHE::BoundaryP[7]

private

Definition at line 180 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), and HadronProtonQ2().

BoundaryTG

G4double G4ElasticHadrNucleusHE::BoundaryTG[7]

private

Definition at line 180 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), and HadronProtonQ2().

BoundaryTL

G4double G4ElasticHadrNucleusHE::BoundaryTL[7]

private

Definition at line 180 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), HadrNucDifferCrSec(), and HadronProtonQ2().

Coeff0

G4double G4ElasticHadrNucleusHE::Coeff0

private

Definition at line 184 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), HadrNucDifferCrSec(), and InterpolateHN().

Coeff1

G4double G4ElasticHadrNucleusHE::Coeff1

private

Definition at line 183 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), HadrNucDifferCrSec(), and InterpolateHN().

Coeff2

G4double G4ElasticHadrNucleusHE::Coeff2

private

Definition at line 183 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), and HadrNucDifferCrSec().

ConstU

G4double G4ElasticHadrNucleusHE::ConstU

private

Definition at line 177 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), GetFt(), and HadronProtonQ2().

DDSect2

G4double G4ElasticHadrNucleusHE::DDSect2

private

Definition at line 177 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), and HadrNucDifferCrSec().

DDSect3

G4double G4ElasticHadrNucleusHE::DDSect3

private

Definition at line 177 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), and HadrNucDifferCrSec().

dQ2

G4double G4ElasticHadrNucleusHE::dQ2

private

Definition at line 166 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), FillFq2(), G4ElasticHadrNucleusHE(), GetQ2_2(), and HadronNucleusQ2_2().

Dtot11

G4double G4ElasticHadrNucleusHE::Dtot11

private

Definition at line 186 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), and HadrNucDifferCrSec().

ekinLowLimit

G4double G4ElasticHadrNucleusHE::ekinLowLimit

private

Definition at line 165 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), and SampleInvariantT().

epCheckLevels

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

privateinherited

Definition at line 198 of file G4HadronicInteraction.hh.

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

fBinom

G4double G4ElasticHadrNucleusHE::fBinom = {{0.0}}

staticprivate

Definition at line 194 of file G4ElasticHadrNucleusHE.hh.

Referenced by Binom(), and GetBinomCof().

fDirectory

char* G4ElasticHadrNucleusHE::fDirectory

private

Definition at line 198 of file G4ElasticHadrNucleusHE.hh.

Referenced by InFileName(), and ~G4ElasticHadrNucleusHE().

fElasticData

G4ElasticData * G4ElasticHadrNucleusHE::fElasticData = {{nullptr}}

staticprivate

Definition at line 196 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), InitialiseModel(), SampleInvariantT(), and ~G4ElasticHadrNucleusHE().

fEnergy

G4double G4ElasticHadrNucleusHE::fEnergy = {0.0}

staticprivate

Definition at line 192 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), and G4ElasticHadrNucleusHE().

fHadronCode

const G4int G4ElasticHadrNucleusHE::fHadronCode

staticprivate

Initial value:

= 
{211,-211,2112,2212,321,-321,130,310,311,-311,
 3122,3222,3112,3212,3312,3322,3334,
 -2212,-2112,-3122,-3222,-3112,-3212,-3312,-3322,-3334}

Definition at line 157 of file G4ElasticHadrNucleusHE.hh.

Referenced by InitialiseModel(), and SampleInvariantT().

fHadronType

const G4int G4ElasticHadrNucleusHE::fHadronType

staticprivate

Initial value:

= 
{2,3,6,0,4,5,4,4,4,5,
 0,0,0,0,0,0,0,
 1,7,1,1,1,1,1,1,1}

Definition at line 158 of file G4ElasticHadrNucleusHE.hh.

Referenced by InitialiseModel(), and SampleInvariantT().

fHadronType1

const G4int G4ElasticHadrNucleusHE::fHadronType1

staticprivate

Initial value:

= 
{3,4,1,0,5,6,5,5,5,6,
 0,0,0,0,0,0,0,
 2,2,2,2,2,2,2,2,2}

Definition at line 159 of file G4ElasticHadrNucleusHE.hh.

Referenced by InitialiseModel(), and SampleInvariantT().

fLineF

G4double G4ElasticHadrNucleusHE::fLineF = {0.0}

staticprivate

Definition at line 191 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), and FillFq2().

fLowEdgeEnergy

G4double G4ElasticHadrNucleusHE::fLowEdgeEnergy = {0.0}

staticprivate

Definition at line 193 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), and HadronNucleusQ2_2().

fRetrieveFromFile

G4bool G4ElasticHadrNucleusHE::fRetrieveFromFile = false

staticprivate

Definition at line 162 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData().

fStoreToFile

G4bool G4ElasticHadrNucleusHE::fStoreToFile = false

staticprivate

Definition at line 161 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData().

HadrEnergy

G4double G4ElasticHadrNucleusHE::HadrEnergy

private

Definition at line 173 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), FillData(), G4ElasticHadrNucleusHE(), HadrNucDifferCrSec(), and HadronProtonQ2().

HadrReIm

G4double G4ElasticHadrNucleusHE::HadrReIm

private

Definition at line 176 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetLightFq2(), and HadrNucDifferCrSec().

HadrSlope

G4double G4ElasticHadrNucleusHE::HadrSlope

private

Definition at line 176 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), and HadrNucDifferCrSec().

HadrTot

G4double G4ElasticHadrNucleusHE::HadrTot

private

Definition at line 176 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetLightFq2(), and HadrNucDifferCrSec().

hLabMomentum

G4double G4ElasticHadrNucleusHE::hLabMomentum

private

Definition at line 171 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), FillData(), G4ElasticHadrNucleusHE(), HadronProtonQ2(), and InterpolateHN().

hLabMomentum2

G4double G4ElasticHadrNucleusHE::hLabMomentum2

private

Definition at line 172 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), FillData(), G4ElasticHadrNucleusHE(), and HadronProtonQ2().

hMass

G4double G4ElasticHadrNucleusHE::hMass

private

Definition at line 169 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), FillData(), G4ElasticHadrNucleusHE(), HadronNucleusQ2_2(), InitialiseModel(), and SampleInvariantT().

hMass2

G4double G4ElasticHadrNucleusHE::hMass2

private

Definition at line 170 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetLightFq2(), HadrNucDifferCrSec(), HadronNucleusQ2_2(), HadronProtonQ2(), InitialiseModel(), and SampleInvariantT().

iHadrCode

G4int G4ElasticHadrNucleusHE::iHadrCode

private

Definition at line 154 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), FillData(), G4ElasticHadrNucleusHE(), InitialiseModel(), and SampleInvariantT().

iHadron

G4int G4ElasticHadrNucleusHE::iHadron

private

Definition at line 155 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), FillData(), G4ElasticHadrNucleusHE(), InitialiseModel(), and SampleInvariantT().

iHadron1

G4int G4ElasticHadrNucleusHE::iHadron1

private

Definition at line 156 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), G4ElasticHadrNucleusHE(), HadronProtonQ2(), InitialiseModel(), and SampleInvariantT().

isBlocked

G4bool G4HadronicInteraction::isBlocked

protectedinherited

Definition at line 188 of file G4HadronicInteraction.hh.

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

isMaster

G4bool G4ElasticHadrNucleusHE::isMaster

private

Definition at line 200 of file G4ElasticHadrNucleusHE.hh.

Referenced by G4ElasticHadrNucleusHE(), InitialiseModel(), and ~G4ElasticHadrNucleusHE().

lowestEnergyLimit

G4double G4HadronElastic::lowestEnergyLimit

privateinherited

Definition at line 87 of file G4HadronElastic.hh.

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

nistManager

G4NistManager* G4ElasticHadrNucleusHE::nistManager

private

Definition at line 197 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), and G4ElasticHadrNucleusHE().

nwarn

G4int G4HadronElastic::nwarn

privateinherited

Definition at line 88 of file G4HadronElastic.hh.

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

pLocalTmax

G4double G4HadronElastic::pLocalTmax

protectedinherited

Definition at line 77 of file G4HadronElastic.hh.

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

Pnucl

G4double G4ElasticHadrNucleusHE::Pnucl

private

Definition at line 189 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), G4ElasticHadrNucleusHE(), GetLightFq2(), and HadrNucDifferCrSec().

Q2max

G4double G4ElasticHadrNucleusHE::Q2max

private

Definition at line 189 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), FillFq2(), G4ElasticHadrNucleusHE(), GetQ2_2(), HadrNucDifferCrSec(), and HadronNucleusQ2_2().

R1

G4double G4ElasticHadrNucleusHE::R1

private

Definition at line 189 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), FillFq2(), G4ElasticHadrNucleusHE(), GetLightFq2(), GetQ2_2(), HadrNucDifferCrSec(), and HadronNucleusQ2_2().

R2

G4double G4ElasticHadrNucleusHE::R2

private

Definition at line 189 of file G4ElasticHadrNucleusHE.hh.

Referenced by FillData(), G4ElasticHadrNucleusHE(), GetLightFq2(), and HadrNucDifferCrSec().

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().

Slope0

G4double G4ElasticHadrNucleusHE::Slope0

private

Definition at line 184 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), HadrNucDifferCrSec(), and InterpolateHN().

Slope1

G4double G4ElasticHadrNucleusHE::Slope1

private

Definition at line 183 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), HadrNucDifferCrSec(), and InterpolateHN().

Slope2

G4double G4ElasticHadrNucleusHE::Slope2

private

Definition at line 183 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), G4ElasticHadrNucleusHE(), GetFt(), GetLightFq2(), and HadrNucDifferCrSec().

theAlpha

G4ParticleDefinition* G4HadronElastic::theAlpha

privateinherited

Definition at line 85 of file G4HadronElastic.hh.

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

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* G4HadronElastic::theDeuteron

privateinherited

Definition at line 84 of file G4HadronElastic.hh.

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

theMaxEnergy

G4double G4HadronicInteraction::theMaxEnergy

protectedinherited

Definition at line 186 of file G4HadronicInteraction.hh.

Referenced by G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4HadronicInteraction::G4HadronicInteraction(), G4hhElastic::G4hhElastic(), G4NuclNuclDiffuseElastic::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::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* G4HadronElastic::theNeutron

privateinherited

Definition at line 83 of file G4HadronElastic.hh.

Referenced by G4HadronElastic::G4HadronElastic().

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().

theProton

G4ParticleDefinition* G4HadronElastic::theProton

privateinherited

Definition at line 82 of file G4HadronElastic.hh.

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

TotP

G4double G4ElasticHadrNucleusHE::TotP

private

Definition at line 176 of file G4ElasticHadrNucleusHE.hh.

Referenced by DefineHadronValues(), and G4ElasticHadrNucleusHE().

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(), DefineHadronValues(), FillData(), FillFq2(), G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4ElasticHadrNucleusHE(), G4EMDissociation::G4EMDissociation(), G4hhElastic::G4hhElastic(), G4NuclNuclDiffuseElastic::G4NuclNuclDiffuseElastic(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), GetFt(), GetLightFq2(), GetQ2_2(), G4HadronicInteraction::GetVerboseLevel(), HadronNucleusQ2_2(), HadronProtonQ2(), G4LFission::init(), G4DiffuseElastic::Initialise(), G4DiffuseElasticV2::Initialise(), G4NuclNuclDiffuseElastic::Initialise(), G4DiffuseElastic::InitialiseOnFly(), G4DiffuseElasticV2::InitialiseOnFly(), G4NuclNuclDiffuseElastic::InitialiseOnFly(), G4CascadeInterface::makeDynamicParticle(), G4CascadeInterface::NoInteraction(), G4CascadeInterface::Propagate(), SampleInvariantT(), G4AntiNuclElastic::SampleThetaCMS(), G4DiffuseElastic::SampleThetaLab(), G4NuclNuclDiffuseElastic::SampleThetaLab(), G4AntiNuclElastic::SampleThetaLab(), G4WilsonAbrasionModel::SetUseAblation(), G4HadronicInteraction::SetVerboseLevel(), G4WilsonAbrasionModel::SetVerboseLevel(), G4DiffuseElastic::ThetaCMStoThetaLab(), G4DiffuseElasticV2::ThetaCMStoThetaLab(), G4NuclNuclDiffuseElastic::ThetaCMStoThetaLab(), G4DiffuseElastic::ThetaLabToThetaCMS(), G4DiffuseElasticV2::ThetaLabToThetaCMS(), and G4NuclNuclDiffuseElastic::ThetaLabToThetaCMS().

---

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

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