G4hhElastic Class Reference

#include <G4hhElastic.hh>

Inheritance diagram for G4hhElastic:

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 &)
void	BuildTableT (G4ParticleDefinition *target, G4ParticleDefinition *projectile)
void	BuildTableTest (G4ParticleDefinition *target, G4ParticleDefinition *projectile, G4double plab)
void	CalculateBQ (G4double b)
void	CalculateBqQ12 (G4double b)
void	CalculateBqQ123 (G4double b)
void	CalculateBqQ13 (G4double b)
G4double	ComputeMomentumCMS (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
void	DeActivateFor (const G4Element *anElement)
void	DeActivateFor (const G4Material *aMaterial)
	G4hhElastic ()
	G4hhElastic (G4ParticleDefinition *target, G4ParticleDefinition *projectile)
	G4hhElastic (G4ParticleDefinition *target, G4ParticleDefinition *projectile, G4double plab)
G4complex	GetAqq ()
G4complex	GetAQQ ()
G4complex	GetAqQ ()
G4double	GetBq ()
G4double	GetBQ ()
G4double	GetBqQ ()
G4double	GetCofF2 ()
G4double	GetCofF3 ()
G4double	GetCofS1 ()
G4double	GetCofS2 ()
G4double	GetCofS3 ()
G4double	GetdsdtF1 (G4double s, G4double q)
G4double	GetdsdtF12 (G4double s, G4double q)
G4double	GetdsdtF123 (G4double q)
G4double	GetdsdtF123qQgG (G4double q)
G4double	GetdsdtF12qQgG (G4double s, G4double q)
G4double	GetdsdtF13qQG (G4double s, G4double q)
G4double	GetdsdtF1qQgG (G4double s, G4double q)
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
G4double	GetExpRatioF123 (G4double s, G4double q)
G4complex	GetF1 (G4double qp)
G4complex	GetF1qQgG (G4double qp)
G4complex	GetF2 (G4double qp)
G4complex	GetF2qQgG (G4double qp)
G4complex	GetF3 (G4double qp)
G4complex	GetF3qQgG (G4double qp)
virtual const std::pair< G4double, G4double >	GetFatalEnergyCheckLevels () const
G4double	GetImCof ()
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	GetOpticalRatio ()
G4double	GetRA ()
G4double	GetRB ()
G4double	GetRecoilEnergyThreshold () const
G4double	GetRg ()
G4double	GetRG ()
G4double	GetRhoReIm ()
G4double	GetRq ()
G4double	GetRQ ()
G4double	GetSlopeCof (const G4int pdg)
G4double	GetSpp ()
G4double	GetTransfer (G4int iMomentum, G4int iTransfer, G4double position)
G4int	GetVerboseLevel () const
void	Initialise ()
virtual void	InitialiseModel ()
virtual G4bool	IsApplicable (const G4HadProjectile &, G4Nucleus &)
G4bool	IsBlocked (const G4Element *anElement) const
G4bool	IsBlocked (const G4Material *aMaterial) const
G4double	LowestEnergyLimit () const
void	ModelDescription (std::ostream &) const override
G4bool	operator!= (const G4HadronicInteraction &right) const =delete
G4bool	operator== (const G4HadronicInteraction &right) const =delete
G4complex	Phi13 ()
G4complex	Phi14 ()
G4complex	Phi23 ()
G4complex	Phi24 ()
G4complex	Pomeron ()
G4double	SampleBisectionalT (const G4ParticleDefinition *p, G4double plab)
G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int, G4int)
G4double	SampleTest (G4double tMin)
void	SetAlphaP (G4double a)
void	SetBq (G4double b)
void	SetBQ (G4double b)
void	SetBqQ (G4double b)
void	SetCofF2 (G4double f)
void	SetCofF3 (G4double f)
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
void	SetEta (G4double E)
void	SetImCof (G4double a)
void	SetLambda (G4double L)
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	SetParameters ()
void	SetParametersCMS (G4double plab)
void	SetRA (G4double rn, G4double pq, G4double pQ)
void	SetRB (G4double rn, G4double pq, G4double pQ)
void	SetRecoilEnergyThreshold (G4double val)
void	SetSigmaTot (G4double stot)
void	SetSpp (G4double spp)
void	SetVerboseLevel (G4int value)
virtual	~G4hhElastic ()

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

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

Private Attributes
std::pair< G4double, G4double >	epCheckLevels
G4double	fAlpha
G4double	fAlphaP
std::vector< G4PhysicsTable * >	fBankT
G4double	fBeta
G4int	fBinT
G4double	fBq
G4double	fBQ
G4double	fBqQ
G4double	fCofF2
G4double	fCofF3
G4double	fDelta
G4int	fEnergyBin
G4PhysicsLogVector *	fEnergyVector
G4double	fEta
G4double	fExpSlope
G4double	fGamma
G4HadronNucleonXsc *	fHadrNuclXsc
G4double	fImCof
G4int	fInTkin
G4double	fLambda
G4double	fLambdaFF
G4double	fMassDif2
G4double	fMassProj
G4double	fMassSum2
G4double	fMassTarg
G4double	fMff2
G4double	fMQ
G4double	fMq
G4double	fOldTkin
G4double	fOptRatio
G4double	fPcms
G4ParticleDefinition *	fProjectile
G4double	fQcof
G4double	fRA
G4double	fRB
G4double	fRG
G4double	fRg
G4double	fRhoReIm
G4double	fRQ
G4double	fRq
G4double	fSigmaTot
G4double	fSo
G4double	fSpp
G4PhysicsTable *	fTableT
G4ParticleDefinition *	fTarget
G4double	lowEnergyLimitHE
G4double	lowEnergyLimitQ
G4double	lowEnergyRecoilLimit
G4double	lowestEnergyLimit
G4int	nwarn
G4double	plabLowLimit
G4double	recoilEnergyThreshold
G4HadronicInteractionRegistry *	registry
G4ParticleDefinition *	theAlpha
std::vector< const G4Material * >	theBlockedList
std::vector< const G4Element * >	theBlockedListElements
G4ParticleDefinition *	theDeuteron
std::vector< std::pair< G4double, const G4Material * > >	theMaxEnergyList
std::vector< std::pair< G4double, const G4Element * > >	theMaxEnergyListElements
std::vector< std::pair< G4double, const G4Material * > >	theMinEnergyList
std::vector< std::pair< G4double, const G4Element * > >	theMinEnergyListElements
G4String	theModelName
G4ParticleDefinition *	theNeutron
G4ParticleDefinition *	thePionMinus
G4ParticleDefinition *	thePionPlus
G4ParticleDefinition *	theProton

Static Private Attributes
static const G4double	theNuclNuclData [19][6]
static const G4double	thePiKaNuclData [8][6]

Detailed Description

Definition at line 61 of file G4hhElastic.hh.

Constructor & Destructor Documentation

G4hhElastic() [1/3]

G4hhElastic::G4hhElastic

(

)

Definition at line 73 of file G4hhElastic.cc.

  : G4HadronElastic("HadrHadrElastic")
{
  SetMinEnergy( 1.*GeV );
  SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
  verboseLevel = 0;
  lowEnergyRecoilLimit = 100.*keV;  
  lowEnergyLimitQ  = 0.0*GeV;  
  lowEnergyLimitHE = 0.0*GeV;  
  lowestEnergyLimit= 0.0*keV;  
  plabLowLimit     = 20.0*MeV;
 
  fRhoReIm=fSigmaTot=fOptRatio=fSpp=fPcms=0.0;
  fInTkin=0;
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  thePionPlus = G4PionPlus::PionPlus();
  thePionMinus= G4PionMinus::PionMinus();
 
  fTarget  = G4Proton::Proton();
  fProjectile  = 0;
  fHadrNuclXsc = new G4HadronNucleonXsc();
 
  fEnergyBin = 200;
  fBinT  = 514; // 514; // 500; // 200;
 
  fEnergyVector =  new G4PhysicsLogVector( theMinEnergy, theMaxEnergy, fEnergyBin );
 
  fTableT = 0;
  fOldTkin = 0.;
  SetParameters();
 
  Initialise();
}

References fBinT, fEnergyBin, fEnergyVector, fHadrNuclXsc, fInTkin, fOldTkin, fOptRatio, fPcms, fProjectile, fRhoReIm, fSigmaTot, fSpp, fTableT, fTarget, G4HadronicInteraction::GetMaxEnergy(), GeV, Initialise(), G4HadronicParameters::Instance(), keV, lowEnergyLimitHE, lowEnergyLimitQ, lowEnergyRecoilLimit, lowestEnergyLimit, MeV, G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), plabLowLimit, G4Proton::Proton(), G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), SetParameters(), G4HadronicInteraction::theMaxEnergy, G4HadronicInteraction::theMinEnergy, theNeutron, thePionMinus, thePionPlus, theProton, and G4HadronicInteraction::verboseLevel.

Referenced by BuildTableT(), and BuildTableTest().

G4hhElastic() [2/3]

G4hhElastic::G4hhElastic	(	G4ParticleDefinition *	target,
		G4ParticleDefinition *	projectile,
		G4double	plab
	)

Definition at line 114 of file G4hhElastic.cc.

  : G4HadronElastic("HadrHadrElastic")
{
  SetMinEnergy( 1.*GeV );
  SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
  verboseLevel         = 0;
  lowEnergyRecoilLimit = 100.*keV;  
  lowEnergyLimitQ      = 0.0*GeV;  
  lowEnergyLimitHE     = 0.0*GeV;  
  lowestEnergyLimit    = 0.0*keV;  
  plabLowLimit         = 20.0*MeV;
 
  fRhoReIm=fSigmaTot=fOptRatio=fSpp=fPcms=0.0;
  fInTkin=0;
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  thePionPlus = G4PionPlus::PionPlus();
  thePionMinus= G4PionMinus::PionMinus();
 
  fTarget      = target;
  fProjectile  = projectile;
  fMassTarg   = fTarget->GetPDGMass();
  fMassProj   = fProjectile->GetPDGMass();
  fMassSum2   = (fMassTarg+fMassProj)*(fMassTarg+fMassProj);
  fMassDif2   = (fMassTarg-fMassProj)*(fMassTarg-fMassProj);
  fHadrNuclXsc = new G4HadronNucleonXsc();
 
  fEnergyBin = 200;
  fBinT      = 514; // 200;
 
  fEnergyVector =  new G4PhysicsLogVector( theMinEnergy, theMaxEnergy, fEnergyBin );
  fTableT       = 0;
  fOldTkin = 0.;
 
 
  SetParameters();
  SetParametersCMS( plab);
}

References fBinT, fEnergyBin, fEnergyVector, fHadrNuclXsc, fInTkin, fMassDif2, fMassProj, fMassSum2, fMassTarg, fOldTkin, fOptRatio, fPcms, fProjectile, fRhoReIm, fSigmaTot, fSpp, fTableT, fTarget, G4HadronicInteraction::GetMaxEnergy(), G4ParticleDefinition::GetPDGMass(), GeV, G4HadronicParameters::Instance(), keV, lowEnergyLimitHE, lowEnergyLimitQ, lowEnergyRecoilLimit, lowestEnergyLimit, MeV, G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), plabLowLimit, G4Proton::Proton(), G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), SetParameters(), SetParametersCMS(), G4HadronicInteraction::theMaxEnergy, G4HadronicInteraction::theMinEnergy, theNeutron, thePionMinus, thePionPlus, theProton, and G4HadronicInteraction::verboseLevel.

G4hhElastic() [3/3]

G4hhElastic::G4hhElastic	(	G4ParticleDefinition *	target,
		G4ParticleDefinition *	projectile
	)

Definition at line 159 of file G4hhElastic.cc.

  : G4HadronElastic("HadrHadrElastic")
{
  SetMinEnergy( 1.*GeV );
  SetMaxEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
  verboseLevel = 0;
  lowEnergyRecoilLimit = 100.*keV;  
  lowEnergyLimitQ  = 0.0*GeV;  
  lowEnergyLimitHE = 0.0*GeV;  
  lowestEnergyLimit= 0.0*keV;  
  plabLowLimit     = 20.0*MeV;
 
  fRhoReIm=fSigmaTot=fOptRatio=fSpp=fPcms=0.0;
  fInTkin=0;
 
  fTarget = target; // later vmg
  fProjectile = projectile;
  theProton   = G4Proton::Proton();
  theNeutron  = G4Neutron::Neutron();
  thePionPlus = G4PionPlus::PionPlus();
  thePionMinus= G4PionMinus::PionMinus();
 
  fTarget  = G4Proton::Proton(); // later vmg
  fMassTarg   = fTarget->GetPDGMass();
  fMassProj   = fProjectile->GetPDGMass();
  fMassSum2   = (fMassTarg+fMassProj)*(fMassTarg+fMassProj);
  fMassDif2   = (fMassTarg-fMassProj)*(fMassTarg-fMassProj);
  fHadrNuclXsc = new G4HadronNucleonXsc();
 
  fEnergyBin = 200;
  fBinT  = 514; // 514; // 500; // 200;
 
  fEnergyVector =  new G4PhysicsLogVector( theMinEnergy, theMaxEnergy, fEnergyBin );
 
  fTableT = 0;
  fOldTkin = 0.;
 
  SetParameters();
}

References fBinT, fEnergyBin, fEnergyVector, fHadrNuclXsc, fInTkin, fMassDif2, fMassProj, fMassSum2, fMassTarg, fOldTkin, fOptRatio, fPcms, fProjectile, fRhoReIm, fSigmaTot, fSpp, fTableT, fTarget, G4HadronicInteraction::GetMaxEnergy(), G4ParticleDefinition::GetPDGMass(), GeV, G4HadronicParameters::Instance(), keV, lowEnergyLimitHE, lowEnergyLimitQ, lowEnergyRecoilLimit, lowestEnergyLimit, MeV, G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), plabLowLimit, G4Proton::Proton(), G4HadronicInteraction::SetMaxEnergy(), G4HadronicInteraction::SetMinEnergy(), SetParameters(), G4HadronicInteraction::theMaxEnergy, G4HadronicInteraction::theMinEnergy, theNeutron, thePionMinus, thePionPlus, theProton, and G4HadronicInteraction::verboseLevel.

~G4hhElastic()

G4hhElastic::~G4hhElastic ( )

virtual

Definition at line 205 of file G4hhElastic.cc.

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

References fBankT, fEnergyVector, fHadrNuclXsc, and fTableT.

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.

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.

{}

BuildTableT()

void G4hhElastic::BuildTableT	(	G4ParticleDefinition *	target,
		G4ParticleDefinition *	projectile
	)

Definition at line 254 of file G4hhElastic.cc.

{
  G4int iTkin, jTransfer;
  G4double plab, Tkin, tMax;
  G4double t1, t2, dt, delta = 0., sum = 0.;
 
  fTarget     = target;
  fProjectile = projectile;
  fMassTarg   = fTarget->GetPDGMass();
  fMassProj   = fProjectile->GetPDGMass();
  fMassSum2   = (fMassTarg+fMassProj)*(fMassTarg+fMassProj);
  fMassDif2   = (fMassTarg-fMassProj)*(fMassTarg-fMassProj);
 
  G4Integrator<G4hhElastic,G4double(G4hhElastic::*)(G4double)> integral;
  // G4HadronNucleonXsc*          hnXsc = new G4HadronNucleonXsc();
  fTableT = new G4PhysicsTable(fEnergyBin);
 
  for( iTkin = 0; iTkin < fEnergyBin; iTkin++)
  {
    Tkin  = fEnergyVector->GetLowEdgeEnergy(iTkin);
    plab  = std::sqrt( Tkin*( Tkin + 2*fMassProj ) );
    // G4DynamicParticle*  theDynamicParticle = new G4DynamicParticle(projectile,
    //                                           G4ParticleMomentum(0.,0.,1.),
    //                                           Tkin);
    // fSigmaTot = fHadrNuclXsc->GetHadronNucleonXscNS( theDynamicParticle, target );
 
    SetParametersCMS( plab );
 
    tMax  = 4.*fPcms*fPcms;
    if( tMax > 15.*GeV*GeV ) tMax = 15.*GeV*GeV; // Check vs. energy ???
 
    G4PhysicsFreeVector* vectorT = new G4PhysicsFreeVector(fBinT-1);
    sum = 0.;
    dt  = tMax/fBinT;
 
    // for(j = 1; j < fBinT; j++)
 
    for( jTransfer = fBinT-1; jTransfer >= 1; jTransfer--)
    {
      t1 = dt*(jTransfer-1);
      t2 = t1 + dt;
 
      if( fMassProj > 900.*MeV ) // pp, pn
      {
        delta = integral.Legendre10(this, &G4hhElastic::GetdsdtF123, t1, t2);
        // delta = integral.Legendre96(this, &G4hhElastic::GetdsdtF123, t1, t2);
      }
      else // pi+-p, K+-p
      {
        delta = integral.Legendre10(this, &G4hhElastic::GetdsdtF123qQgG, t1, t2);
        // delta = integral.Legendre96(this, &G4hhElastic::GetdsdtF123qQgG, t1, t2);
      }
      sum += delta;
      vectorT->PutValue( jTransfer-1, t1, sum ); // t2
    }
    // vectorT->PutValue( fBinT-1, dt*(fBinT-1), 0. ); // t2
    fTableT->insertAt( iTkin, vectorT );
    // delete theDynamicParticle;
  }
  // delete hnXsc;
 
  return;
}

References fBinT, fEnergyBin, fEnergyVector, fMassDif2, fMassProj, fMassSum2, fMassTarg, fPcms, fProjectile, fTableT, fTarget, G4hhElastic(), GetdsdtF123(), GetdsdtF123qQgG(), G4PhysicsVector::GetLowEdgeEnergy(), G4ParticleDefinition::GetPDGMass(), GeV, G4PhysicsTable::insertAt(), MeV, G4PhysicsFreeVector::PutValue(), and SetParametersCMS().

Referenced by Initialise().

BuildTableTest()

void G4hhElastic::BuildTableTest	(	G4ParticleDefinition *	target,
		G4ParticleDefinition *	projectile,
		G4double	plab
	)

Definition at line 520 of file G4hhElastic.cc.

{
  G4int jTransfer;
  G4double tMax; // , sQq, sQG;
  G4double t1, t2, dt, delta = 0., sum = 0. ; // , threshold;
 
  fTarget     = target;
  fProjectile = projectile;
  fMassTarg =  fTarget->GetPDGMass();
  fMassProj =  fProjectile->GetPDGMass();
  fMassSum2 = (fMassTarg+fMassProj)*(fMassTarg+fMassProj);
  fMassDif2 = (fMassTarg-fMassProj)*(fMassTarg-fMassProj);
  fSpp  = fMassProj*fMassProj + fMassTarg*fMassTarg + 2.*fMassTarg*std::sqrt(plab*plab + fMassProj*fMassProj);
  fPcms = std::sqrt( (fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
 
  G4cout<<"fMassTarg = "<<fMassTarg<<" MeV; fMassProj = "<<fMassProj<<" MeV"<<G4endl;
  tMax = 4.*fPcms*fPcms;
  if( tMax > 15.*GeV*GeV ) tMax = 15.*GeV*GeV; // Check vs. energy ???
 
 
  G4Integrator<G4hhElastic,G4double(G4hhElastic::*)(G4double)> integral;
  fTableT = new G4PhysicsTable(1);
  G4PhysicsFreeVector* vectorT = new G4PhysicsFreeVector(fBinT-1);
 
  sum = 0.;
  dt = tMax/G4double(fBinT);
  G4cout<<"s = "<<std::sqrt(fSpp)/GeV<<" GeV; fPcms = "<<fPcms/GeV
    <<" GeV; qMax = "<<tMax/GeV/GeV<<" GeV2; dt = "<<dt/GeV/GeV<<" GeV2"<<G4endl;
 
  // G4cout<<"fRA = "<<fRA*GeV<<"; fRB = "<<fRB*GeV<<G4endl;
 
    // for(jTransfer = 1; jTransfer < fBinT; jTransfer++)
  for( jTransfer = fBinT-1; jTransfer >= 1; jTransfer-- )
  {
    t1 = dt*(jTransfer-1);
    t2 = t1 + dt;
 
    if( fMassProj > 900.*MeV ) // pp, pn
    {
      delta = integral.Legendre10(this, &G4hhElastic::GetdsdtF123, t1, t2);
      // threshold = integral.Legendre96(this, &G4hhElastic::GetdsdtF123, t1, tMax);
    }
    else // pi+-p, K+-p
    {
      delta = integral.Legendre10(this, &G4hhElastic::GetdsdtF123qQgG, t1, t2);
      // threshold = integral.Legendre96(this, &G4hhElastic::GetdsdtF123qQgG, t1, tMax);
      // delta = integral.Legendre96(this, &G4hhElastic::GetdsdtF123, t1, t2);
    }
    sum += delta;
    // G4cout<<delta<<"\t"<<sum<<"\t"<<threshold<<G4endl;   
  
    // sQq = GetdsdtF123(q1);
    // sQG = GetdsdtF123qQgG(q1);
    // G4cout<<q1/GeV<<"\t"<<sQG*GeV*GeV/millibarn<<"\t"<<sQq*GeV*GeV/millibarn<<G4endl;
    // G4cout<<"sum = "<<sum<<", "; 
    
    vectorT->PutValue( jTransfer-1, t1, sum );  // t2
  }
  // vectorT->PutValue( fBinT-1, dt*(fBinT-1), 0. ); // t2
  fTableT->insertAt( 0, vectorT );
  fBankT.push_back( fTableT );  // 0
 
  // for(jTransfer = 0; jTransfer < fBinT-1; jTransfer++) 
  //   G4cout<<(*(*fTableT)(0))(jTransfer)/sum<<"\t\t"<<G4Pow::GetInstance()->powN(2.,-jTransfer)<<G4endl;
  
  return;
}

References fBankT, fBinT, fMassDif2, fMassProj, fMassSum2, fMassTarg, fPcms, fProjectile, fSpp, fTableT, fTarget, G4cout, G4endl, G4hhElastic(), GetdsdtF123(), GetdsdtF123qQgG(), G4ParticleDefinition::GetPDGMass(), GeV, G4PhysicsTable::insertAt(), MeV, and G4PhysicsFreeVector::PutValue().

CalculateBQ()

void G4hhElastic::CalculateBQ ( G4double  b )

inline

Definition at line 957 of file G4hhElastic.hh.

{
  fBq = b1;
  G4double s1 = GetCofS1();
  G4double s2 = GetCofS2();
  G4double s3 = GetCofS3();
  G4double sqrtBq = std::sqrt(fBq);
 
  // cofs of the fBQ 3rd equation
 
  G4double a = s3*sqrtBq;
  G4double b = s1*fBq - 1.;
  G4double c = (s2*fBq - 2.)*sqrtBq;
  G4double d = 1. - fBq;
 
  // cofs of the incomplete 3rd equation
 
  G4double p  = c/a;
           p -= b*b/a/a/3.;
  G4double q  = d/a;
           q -= b*c/a/a/3.;
           q += 2*b*b*b/a/a/a/27.;
 
  // cofs for the incomplete colutions
 
  G4double D  = p*p*p/3./3./3.;
           D += q*q/2./2.;
  G4complex A1 = G4complex(- q/2., std::sqrt(-D) );
  G4complex A  = std::pow(A1,1./3.);
 
  G4complex B1 = G4complex(- q/2., -std::sqrt(-D) );
  G4complex B  = std::pow(B1,1./3.);
 
  // roots of the incomplete 3rd equation
 
  G4complex y1 = A + B;
  G4complex y2 = -0.5*(A + B) + 0.5*std::sqrt(3.)*(A - B)*G4complex(0.,1.);
  G4complex y3 = -0.5*(A + B) - 0.5*std::sqrt(3.)*(A - B)*G4complex(0.,1.);
 
  G4complex x1 = y1 - b/a/3.;
  G4complex x2 = y2 - b/a/3.;
  G4complex x3 = y3 - b/a/3.;
 
  G4cout<<"re_x1 = "<<real(x1)<<"; re_x2 = "<<real(x2)<<"; re_x3 = "<<real(x3)<<G4endl;
  G4cout<<"im_x1 = "<<imag(x1)<<"; im_x2 = "<<imag(x2)<<"; im_x3 = "<<imag(x3)<<G4endl;
 
  G4double r1 = real(x1)*real(x1);
  G4double r2 = real(x2)*real(x2);
  G4double r3 = real(x3)*real(x3);
 
  if( r1 <= 1. && r1 >= 0. )      fBQ = r1;
  else if( r2 <= 1. && r2 >= 0. ) fBQ = r2;
  else if( r3 <= 1. && r3 >= 0. )      fBQ = r3;
  else fBQ = 1.;
  // fBQ = real(x3)*real(x3);
  G4double sqrtBqBQ = std::sqrt(fBq*fBQ);
  fOptRatio  = fBq + fBQ + 2.*sqrtBqBQ - 1.;
  fOptRatio /= sqrtBqBQ*( GetCofS1()*sqrtBqBQ + GetCofS2()*fBq + GetCofS3()*fBQ );
  G4cout<<"F123, fOptRatio = "<<fOptRatio<<G4endl;
  
  return ;
}

References A, B(), D(), fBq, fBQ, fOptRatio, G4cout, G4endl, GetCofS1(), GetCofS2(), and GetCofS3().

CalculateBqQ12()

void G4hhElastic::CalculateBqQ12 ( G4double  b )

inline

Definition at line 786 of file G4hhElastic.hh.

{
  fBq = b1;
 
  G4complex z1324 = G4complex(1./8./CLHEP::pi,0.);
            z1324 /= Phi13() + Phi24() + fLambda + fEta;
        G4double c1324 = real(z1324)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex z1423 = G4complex(1./8./CLHEP::pi,0.);
            z1423 /= Phi14() + Phi23() + fLambda + fEta;
        G4double c1423 = real(z1423)/(CLHEP::hbarc*CLHEP::hbarc);
 
  fBQ  = 1. - 2.*fBq;
  fBQ /= 2. - fSigmaTot*fBq*(c1324+1423); 
 
  G4double ratio  = 2.*(fBq + fBQ) - fSigmaTot*fBq*fBQ*(c1324 + c1423); 
  G4cout<<"ratio = "<<ratio<<G4endl;
 
  return ;
}

References fBq, fBQ, fEta, fLambda, fSigmaTot, G4cout, G4endl, CLHEP::hbarc, Phi13(), Phi14(), Phi23(), Phi24(), and CLHEP::pi.

CalculateBqQ123()

void G4hhElastic::CalculateBqQ123 ( G4double  b )

inline

Definition at line 812 of file G4hhElastic.hh.

{
  fBq = b1;
 
  G4complex z1324 = fCofF2*G4complex(1./8./CLHEP::pi,0.);
            z1324 /= Phi13() + Phi24() + fLambda + fEta;
        G4double c1324 = real(z1324)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex z1423 = fCofF2*G4complex(1./8./CLHEP::pi,0.);
            z1423 /= Phi14() + Phi23() + fLambda + fEta;
        G4double c1423 = real(z1423)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex z1314 = fCofF3*G4complex(1./8./CLHEP::pi,0.);
            z1314 /= Phi13() + Phi14() + fEta;
        G4double c1314 = real(z1314)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex z2324 = fCofF3*G4complex(1./8./CLHEP::pi,0.);
            z2324 /= Phi23() + Phi24() + fEta;
        G4double c2324 = real(z2324)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex z1323 = fCofF3*G4complex(1./8./CLHEP::pi,0.);
            z1323 /= Phi13() + Phi23() + fLambda;
        G4double c1323 = real(z1323)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex z1424 = fCofF3*G4complex(1./8./CLHEP::pi,0.);
            z1424 /= Phi14() + Phi24() + fLambda;
        G4double c1424 = real(z1424)/(CLHEP::hbarc*CLHEP::hbarc);
 
  G4double A = fSigmaTot*c2324;
  G4double B = fSigmaTot*fBq*(c1324 + c1423 + c1323 + c1424) - 2.;
  G4double C = 1. + fSigmaTot*fBq*fBq*c1314 - 2*fBq;
  G4cout<<"A = "<<A<<"; B = "<<B<<"; C = "<<C<<G4endl;
  G4cout<<"determinant = "<<B*B-4.*A*C<<G4endl;
 
  G4double x1 = ( -B - std::sqrt(B*B-4.*A*C) )/2./A;
  G4double x2 = ( -B + std::sqrt(B*B-4.*A*C) )/2./A;
  G4cout<<"x1 = "<<x1<<"; x2 = "<<x2<<G4endl;
 
  if( B*B-4.*A*C < 1.e-6 ) fBQ = std::abs(-B/2./A);
  else if ( B < 0.)        fBQ = std::abs( ( -B - std::sqrt(B*B-4.*A*C) )/2./A);
  else                     fBQ  = std::abs( ( -B + std::sqrt(B*B-4.*A*C) )/2./A);
 
  fOptRatio  = 2*(fBq+fBQ) - fSigmaTot*fBq*fBQ*(c1324 + c1423 + c1323 + c1424);
  fOptRatio -= fSigmaTot*fBq*fBq*c1314 + fSigmaTot*c2324*fBQ*fBQ;
  G4cout<<"BqQ123, fOptRatio = "<<fOptRatio<<G4endl;
 
  return ;
}

References A, B(), C(), fBq, fBQ, fCofF2, fCofF3, fEta, fLambda, fOptRatio, fSigmaTot, G4cout, G4endl, CLHEP::hbarc, Phi13(), Phi14(), Phi23(), Phi24(), and CLHEP::pi.

CalculateBqQ13()

void G4hhElastic::CalculateBqQ13 ( G4double  b )

inline

Definition at line 763 of file G4hhElastic.hh.

{
  fBQ = b2;
 
  G4complex z1424 = G4complex(1./8./CLHEP::pi,0.);
            z1424 /= Phi14() + Phi24() + fAlpha;
        G4double c1424 = real(z1424)/(CLHEP::hbarc*CLHEP::hbarc);
 
  fBqQ  = 1. - fBQ;
  fBQ /= 1. - fSigmaTot*fBQ*c1424; 
 
  G4cout<<"fSigmaTot*fBQ*c1424 = "<<fSigmaTot*fBQ*c1424<<G4endl;
 
  G4double ratio  = fBqQ + fBQ - fSigmaTot*fBqQ*fBQ*c1424; 
  G4cout<<"ratio = "<<ratio<<G4endl;
 
  return ;
}

References fAlpha, fBQ, fBqQ, fSigmaTot, G4cout, G4endl, CLHEP::hbarc, Phi14(), Phi24(), and CLHEP::pi.

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

GetAqq()

G4complex G4hhElastic::GetAqq ( )

inline

Definition at line 867 of file G4hhElastic.hh.

{
  G4double re, im;
 
  re = fRq*fRq/8. + fAlphaP*G4Log(fSpp/fSo) + 8.*fLambda/9.;
  im = -0.5*fAlphaP*fImCof*CLHEP::pi;
  return G4complex(re,im);
}

References fAlphaP, fImCof, fLambda, fRq, fSo, fSpp, G4Log(), and CLHEP::pi.

Referenced by GetAqQ(), GetCofS2(), GetF1(), GetF2(), and GetF3().

GetAQQ()

G4complex G4hhElastic::GetAQQ ( )

inline

Definition at line 880 of file G4hhElastic.hh.

{
  G4double re, im;
 
  re = fRQ*fRQ/8. + fAlphaP*G4Log(fSpp/fSo) + 2.*fLambda/9.;
  im = -0.5*fAlphaP*fImCof*CLHEP::pi;
  return G4complex(re,im);
}

References fAlphaP, fImCof, fLambda, fRQ, fSo, fSpp, G4Log(), and CLHEP::pi.

Referenced by GetAqQ(), GetCofS3(), GetF1(), GetF2(), and GetF3().

GetAqQ()

G4complex G4hhElastic::GetAqQ ( )

inline

Definition at line 893 of file G4hhElastic.hh.

{
  G4complex z = 0.5*( GetAqq() + GetAQQ() );
  return z;
}

References GetAqq(), and GetAQQ().

Referenced by GetCofS1(), GetCofS2(), GetCofS3(), GetF1(), GetF2(), and GetF3().

GetBq()

G4double G4hhElastic::GetBq ( )

inline

Definition at line 169 of file G4hhElastic.hh.

{ return fBq;};

References fBq.

GetBQ()

G4double G4hhElastic::GetBQ ( )

inline

Definition at line 170 of file G4hhElastic.hh.

{ return fBQ;};

References fBQ.

GetBqQ()

G4double G4hhElastic::GetBqQ ( )

inline

Definition at line 171 of file G4hhElastic.hh.

{ return fBqQ;};

References fBqQ.

GetCofF2()

G4double G4hhElastic::GetCofF2 ( )

inline

Definition at line 191 of file G4hhElastic.hh.

{return fCofF2;};

References fCofF2.

GetCofF3()

G4double G4hhElastic::GetCofF3 ( )

inline

Definition at line 192 of file G4hhElastic.hh.

{return fCofF3;};

References fCofF3.

GetCofS1()

G4double G4hhElastic::GetCofS1 ( )

inline

Definition at line 903 of file G4hhElastic.hh.

{
  G4complex z = 1./( GetAqQ() + 4.*fLambda/9. );
  G4double result = real(z);
  result /= 4.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc;
  result *= fSigmaTot*fCofF2;
  return result;
}

References fCofF2, fLambda, fSigmaTot, GetAqQ(), CLHEP::hbarc, and CLHEP::pi.

Referenced by CalculateBQ().

GetCofS2()

G4double G4hhElastic::GetCofS2 ( )

inline

Definition at line 916 of file G4hhElastic.hh.

{
  G4complex z = 1./( GetAqq() + GetAqQ() - 4.*fLambda/9. );
  G4double result = real(z);
  result /= 4.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc;
  result *= fSigmaTot*fCofF3;
  return result;
}

References fCofF3, fLambda, fSigmaTot, GetAqq(), GetAqQ(), CLHEP::hbarc, and CLHEP::pi.

Referenced by CalculateBQ().

GetCofS3()

G4double G4hhElastic::GetCofS3 ( )

inline

Definition at line 929 of file G4hhElastic.hh.

{
  G4complex z = 1./( GetAQQ() + GetAqQ() + 2.*fLambda/9. );
  G4double result = real(z);
  result /= 4.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc;
  result *= fSigmaTot*fCofF3;
  return result;
}

References fCofF3, fLambda, fSigmaTot, GetAQQ(), GetAqQ(), CLHEP::hbarc, and CLHEP::pi.

Referenced by CalculateBQ().

GetdsdtF1()

G4double G4hhElastic::GetdsdtF1 ( G4double  s, G4double  q  )

inline

Definition at line 1042 of file G4hhElastic.hh.

{
  fSpp = spp;
  G4double p = std::sqrt(0.25*spp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
  G4complex F1 = GetF1(t);
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F1)*real(F1) + imag(F1)*imag(F1);  
  return dsdt;
}

References fSpp, GetF1(), CLHEP::pi, and CLHEP::proton_mass_c2.

GetdsdtF12()

G4double G4hhElastic::GetdsdtF12 ( G4double  s, G4double  q  )

inline

Definition at line 1085 of file G4hhElastic.hh.

{
  fSpp = spp;
  G4double p = std::sqrt(0.25*spp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
  G4complex F1 = GetF1(t) - GetF2(t);
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F1)*real(F1) + imag(F1)*imag(F1);  
  return dsdt;
}

References fSpp, GetF1(), GetF2(), CLHEP::pi, and CLHEP::proton_mass_c2.

GetdsdtF123()

G4double G4hhElastic::GetdsdtF123 ( G4double  q )

inline

Definition at line 1128 of file G4hhElastic.hh.

{
  G4double  p   = std::sqrt(0.25*fSpp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
  G4complex F1  = GetF1(t);
            F1 -= fCofF2*GetF2(t);
            F1 -= fCofF3*GetF3(t);
  G4double dsdt = CLHEP::pi/p/p;
  dsdt         *= real(F1)*real(F1) + imag(F1)*imag(F1);  
  return dsdt;
}

References fCofF2, fCofF3, fSpp, GetF1(), GetF2(), GetF3(), CLHEP::pi, and CLHEP::proton_mass_c2.

Referenced by BuildTableT(), and BuildTableTest().

GetdsdtF123qQgG()

G4double G4hhElastic::GetdsdtF123qQgG ( G4double  q )

inline

Definition at line 746 of file G4hhElastic.hh.

{
  G4double p = std::sqrt( (fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp );
 
  G4complex F123  = GetF1qQgG(t); // - fCofF2*GetF2qQgG(t) - fCofF3*GetF3qQgG(t);
            F123 -= fCofF2*GetF2qQgG(t);
            F123 -= fCofF3*GetF3qQgG(t);
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F123)*real(F123) + imag(F123)*imag(F123);  
  return dsdt;
}

References fCofF2, fCofF3, fMassDif2, fMassSum2, fSpp, GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), and CLHEP::pi.

Referenced by BuildTableT(), and BuildTableTest().

GetdsdtF12qQgG()

G4double G4hhElastic::GetdsdtF12qQgG ( G4double  s, G4double  q  )

inline

Definition at line 676 of file G4hhElastic.hh.

{
  fSpp = spp;
  G4double p = std::sqrt((fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
 
  G4complex F12 = GetF1qQgG(t) - GetF2qQgG(t);
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F12)*real(F12) + imag(F12)*imag(F12);  
  return dsdt;
}

References F12, fMassDif2, fMassSum2, fSpp, GetF1qQgG(), GetF2qQgG(), and CLHEP::pi.

GetdsdtF13qQG()

G4double G4hhElastic::GetdsdtF13qQG ( G4double  s, G4double  q  )

inline

Definition at line 584 of file G4hhElastic.hh.

{
  fSpp = spp;
  G4double p = std::sqrt((fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
  G4double k = p/CLHEP::hbarc;
 
  G4complex exp14 = fBqQ*std::exp(-(Phi14() + fBeta*fBeta*fLambda + fGamma*fGamma*fEta)*t);
  G4complex exp24 = fBQ*std::exp(-(Phi24() + fAlpha*fAlpha*fLambda + fGamma*fGamma*fEta)*t);
 
  G4complex F1  = exp14 + exp24;
 
  F1 *=  0.25*k*fSigmaTot/CLHEP::pi;
  F1 *= G4complex(0.,1.);
 
  // 1424
 
  G4complex z1424  = -(Phi24() + fAlpha*fLambda)*(Phi24() + fAlpha*fLambda);
            z1424 /= Phi14() + Phi24() + fLambda;
            z1424 += Phi24() + fAlpha*fAlpha*fLambda + fGamma*fGamma*fEta;
 
 
  G4complex exp1424  = std::exp(-z1424*t);
            exp1424 /= Phi14() + Phi24() + fLambda;
 
  G4complex F3  = fBqQ*fBQ*exp1424;
 
 
  F3 *=  0.25*k/CLHEP::pi;
  F3 *= G4complex(0.,1.);
  F3 *= fSigmaTot*fSigmaTot/(8.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc);
 
  G4complex F13  = F1 - F3;
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F13)*real(F13) + imag(F13)*imag(F13);  
 
  return dsdt;
}

References F13, fAlpha, fBeta, fBQ, fBqQ, fEta, fGamma, fLambda, fMassDif2, fMassSum2, fSigmaTot, fSpp, CLHEP::hbarc, Phi14(), Phi24(), and CLHEP::pi.

GetdsdtF1qQgG()

G4double G4hhElastic::GetdsdtF1qQgG ( G4double  s, G4double  q  )

inline

Definition at line 627 of file G4hhElastic.hh.

{
  fSpp = spp;
  G4double p = std::sqrt((fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
 
  G4complex F1 = GetF1qQgG(t);
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F1)*real(F1) + imag(F1)*imag(F1);  
  return dsdt;
}

References fMassDif2, fMassSum2, fSpp, GetF1qQgG(), and CLHEP::pi.

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

GetExpRatioF123()

G4double G4hhElastic::GetExpRatioF123 ( G4double  s, G4double  q  )

inline

Definition at line 1143 of file G4hhElastic.hh.

{
  fSpp = spp;
  G4double p = std::sqrt(0.25*spp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
 
  // qQ-ds/dt
 
  G4complex F1 = GetF1(t) - fCofF2*GetF2(t) - fCofF3*GetF3(t);
 
  G4double dsdt = CLHEP::pi/p/p;
  dsdt *= real(F1)*real(F1) + imag(F1)*imag(F1); 
 
  // exponent ds/dt
 
  G4complex F10 = GetF1(0.) - fCofF2*GetF2(0.) - fCofF3*GetF3(0.);
 
  fRhoReIm = real(F10)/imag(F10);
 
  G4double dsdt0 = CLHEP::pi/p/p;
  dsdt0 *= real(F10)*real(F10) + imag(F10)*imag(F10); 
 
  dsdt0 *= G4Exp(-fExpSlope*t);
 
  G4double ratio = dsdt/dsdt0;
 
  return ratio;
}

References F10, fCofF2, fCofF3, fExpSlope, fRhoReIm, fSpp, G4Exp(), GetF1(), GetF2(), GetF3(), CLHEP::pi, and CLHEP::proton_mass_c2.

GetF1()

G4complex G4hhElastic::GetF1 ( G4double  qp )

inline

Definition at line 1024 of file G4hhElastic.hh.

{
  G4double p = std::sqrt(0.25*fSpp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
  G4double  k    = p/CLHEP::hbarc;
  G4complex exp1 = fBq*std::exp(-GetAqq()*t);
  G4complex exp2 = fBQ*std::exp(-GetAQQ()*t);
  G4complex exp3 = 2.*std::sqrt(fBq*fBQ)*std::exp(-GetAqQ()*t);
 
  G4complex res  = exp1 + exp2 + exp3;
  res *=  0.25*k*fSigmaTot/CLHEP::pi;
  res *= G4complex(0.,1.);
  return res;
}

References fBq, fBQ, fSigmaTot, fSpp, GetAqq(), GetAQQ(), GetAqQ(), CLHEP::hbarc, CLHEP::pi, and CLHEP::proton_mass_c2.

Referenced by GetdsdtF1(), GetdsdtF12(), GetdsdtF123(), and GetExpRatioF123().

GetF1qQgG()

G4complex G4hhElastic::GetF1qQgG ( G4double  qp )

inline

Definition at line 562 of file G4hhElastic.hh.

{
  G4double p = std::sqrt((fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
  G4double k = p/CLHEP::hbarc;
 
  G4complex exp13 = fBq*std::exp(-(Phi13() + fBeta*fBeta*fLambda + fDelta*fDelta*fEta)*t);
  G4complex exp14 = fBq*std::exp(-(Phi14() + fBeta*fBeta*fLambda + fGamma*fGamma*fEta)*t);
  G4complex exp23 = fBQ*std::exp(-(Phi23() + fAlpha*fAlpha*fLambda + fDelta*fDelta*fEta)*t);
  G4complex exp24 = fBQ*std::exp(-(Phi24() + fAlpha*fAlpha*fLambda + fGamma*fGamma*fEta)*t);
 
  G4complex res  = exp13 + exp14 + exp23 + exp24;
 
  res *=  0.25*k*fSigmaTot/CLHEP::pi;
  res *= G4complex(0.,1.);
 
  return res;
}

References fAlpha, fBeta, fBq, fBQ, fDelta, fEta, fGamma, fLambda, fMassDif2, fMassSum2, fSigmaTot, fSpp, CLHEP::hbarc, Phi13(), Phi14(), Phi23(), Phi24(), and CLHEP::pi.

Referenced by GetdsdtF123qQgG(), GetdsdtF12qQgG(), and GetdsdtF1qQgG().

GetF2()

G4complex G4hhElastic::GetF2 ( G4double  qp )

inline

Definition at line 1057 of file G4hhElastic.hh.

{
  G4double p = std::sqrt(0.25*fSpp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
  G4double  k   = p/CLHEP::hbarc;
  G4complex z1  = GetAqq()*GetAQQ() - 16.*fLambda*fLambda/81.;
            z1 /= 2.*(GetAqQ() + 4.*fLambda/9.);
  G4complex exp1 = std::exp(-z1*t);
 
  G4complex z2 = 0.5*( GetAqQ() - 4.*fLambda/9.);
 
  G4complex exp2 = std::exp(-z2*t);
 
  G4complex res  = exp1 + exp2;
 
  G4complex z3 = GetAqQ() + 4.*fLambda/9.;
 
  res *=  0.25*k/CLHEP::pi;
  res *= G4complex(0.,1.);
  res /= z3;
  res *= fBq*fBQ*fSigmaTot*fSigmaTot/(8.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc);
 
  return res;
}

References fBq, fBQ, fLambda, fSigmaTot, fSpp, GetAqq(), GetAQQ(), GetAqQ(), CLHEP::hbarc, CLHEP::pi, and CLHEP::proton_mass_c2.

Referenced by GetdsdtF12(), GetdsdtF123(), and GetExpRatioF123().

GetF2qQgG()

G4complex G4hhElastic::GetF2qQgG ( G4double  qp )

inline

Definition at line 643 of file G4hhElastic.hh.

{
  G4double p = std::sqrt((fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
  G4double k = p/CLHEP::hbarc;
 
  G4complex z1324  = -(Phi24() + fAlpha*fLambda + fGamma*fEta)*(Phi24() + fAlpha*fLambda + fGamma*fEta);
            z1324 /= Phi13() + Phi24() + fLambda + fEta;
            z1324 += Phi24() + fAlpha*fAlpha*fLambda + fGamma*fGamma*fEta;
 
  G4complex exp1324  = std::exp(-z1324*t);
            exp1324 /= Phi13() + Phi24() + fLambda + fEta;
 
  G4complex z1423 = -(Phi23() + fAlpha*fLambda + fDelta*fEta)*(Phi24() + fAlpha*fLambda + fDelta*fEta);;
            z1423 /= Phi14() + Phi23() + fLambda + fEta;
            z1423 += Phi23() + fAlpha*fAlpha*fLambda + fDelta*fDelta*fEta;
 
  G4complex exp1423 = std::exp(-z1423*t);
            exp1423 /= Phi14() + Phi23() + fLambda + fEta; 
 
  G4complex res  = exp1324 + exp1423;
 
 
  res *=  0.25*k/CLHEP::pi;
  res *= G4complex(0.,1.);
  res *= fBq*fBQ*fSigmaTot*fSigmaTot/(8.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc); // or 4. ???
 
  return res;
}

References fAlpha, fBq, fBQ, fDelta, fEta, fGamma, fLambda, fMassDif2, fMassSum2, fSigmaTot, fSpp, CLHEP::hbarc, Phi13(), Phi14(), Phi23(), Phi24(), and CLHEP::pi.

Referenced by GetdsdtF123qQgG(), and GetdsdtF12qQgG().

GetF3()

G4complex G4hhElastic::GetF3 ( G4double  qp )

inline

Definition at line 1100 of file G4hhElastic.hh.

{
  G4double p = std::sqrt(0.25*fSpp - CLHEP::proton_mass_c2*CLHEP::proton_mass_c2);
  G4double  k   = p/CLHEP::hbarc;
  G4complex z1  = GetAqq()*GetAqQ() - 4.*fLambda*fLambda/81.;
            z1 /= GetAqq() + GetAqQ() - 4.*fLambda/9.;
 
  G4complex exp1 = std::exp(-z1*t)*fBq/(GetAqq() + GetAqQ() - 4.*fLambda/9.);
 
  G4complex z2 = GetAqQ()*GetAQQ() - 1.*fLambda*fLambda/81.;
            z2 /= GetAQQ() + GetAqQ() + 2.*fLambda/9.;
 
  G4complex exp2 = std::exp(-z2*t)*fBQ/(GetAQQ() + GetAqQ() + 2.*fLambda/9.);
 
  G4complex res  = exp1 + exp2;
 
 
  res *=  0.25*k/CLHEP::pi;
  res *= G4complex(0.,1.);
  res *= std::sqrt(fBq*fBQ)*fSigmaTot*fSigmaTot/(4.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc);
 
  return res;
}

References fBq, fBQ, fLambda, fSigmaTot, fSpp, GetAqq(), GetAQQ(), GetAqQ(), CLHEP::hbarc, CLHEP::pi, and CLHEP::proton_mass_c2.

Referenced by GetdsdtF123(), and GetExpRatioF123().

GetF3qQgG()

G4complex G4hhElastic::GetF3qQgG ( G4double  qp )

inline

Definition at line 692 of file G4hhElastic.hh.

{
  G4double p = std::sqrt( (fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
  G4double k = p/CLHEP::hbarc;
 
            // 1314
 
  G4complex z1314  = -(Phi14() + fGamma*fEta)*(Phi14() + fGamma*fEta);
            z1314 /= Phi13() + Phi14() + fEta;
            z1314 += Phi14() + fBeta*fBeta*fLambda + fGamma*fGamma*fEta;
 
  G4complex exp1314  = std::exp(-z1314*t);
            exp1314 /= Phi13() + Phi14() + fEta;
 
        // 2324
 
  G4complex z2324 = -(Phi24() + fGamma*fEta)*(Phi24() + fGamma*fEta);;
            z2324 /= Phi24() + Phi23() + fEta;
            z2324 += Phi24() + fAlpha*fAlpha*fLambda + fGamma*fGamma*fEta;
 
  G4complex exp2324 = std::exp(-z2324*t);
            exp2324 /= Phi24() + Phi23() + fEta; 
 
        // 1323
 
  G4complex z1323  = -(Phi23() + fAlpha*fLambda)*(Phi23() + fAlpha*fLambda);
            z1323 /= Phi13() + Phi23() + fLambda;
            z1323 += Phi23() + fAlpha*fAlpha*fLambda + fDelta*fDelta*fEta;
 
  G4complex exp1323  = std::exp(-z1323*t);
            exp1323 /= Phi13() + Phi23() + fLambda;
 
        // 1424
 
  G4complex z1424  = -(Phi24() + fAlpha*fLambda)*(Phi24() + fAlpha*fLambda);
            z1424 /= Phi14() + Phi24() + fLambda;
            z1424 += Phi24() + fAlpha*fAlpha*fLambda + fGamma*fGamma*fEta;
 
  G4complex exp1424  = std::exp(-z1424*t);
            exp1424 /= Phi14() + Phi24() + fLambda;
 
        G4complex res  = fBq*fBq*exp1314 + fBQ*fBQ*exp2324 + fBq*fBQ*exp1323 + fBq*fBQ*exp1424;
 
  res *=  0.25*k/CLHEP::pi;
  res *= G4complex(0.,1.);
  res *= fSigmaTot*fSigmaTot/(8.*CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc);
 
  return res;
}

References fAlpha, fBeta, fBq, fBQ, fDelta, fEta, fGamma, fLambda, fMassDif2, fMassSum2, fSigmaTot, fSpp, CLHEP::hbarc, Phi13(), Phi14(), Phi23(), Phi24(), and CLHEP::pi.

Referenced by GetdsdtF123qQgG().

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

GetImCof()

G4double G4hhElastic::GetImCof ( )

inline

Definition at line 186 of file G4hhElastic.hh.

{return fImCof;};

References fImCof.

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

GetOpticalRatio()

G4double G4hhElastic::GetOpticalRatio ( )

inline

Definition at line 942 of file G4hhElastic.hh.

{
  return fOptRatio;
  // G4double sqrtBqBQ = std::sqrt(fBq*fBQ);
  // G4double result = fBq + fBQ + 2.*sqrtBqBQ - 1.;
  // result /= sqrtBqBQ*( GetCofS1()*sqrtBqBQ + GetCofS2()*fBq + GetCofS3()*fBQ );
  // return result;
}

References fOptRatio.

GetRA()

G4double G4hhElastic::GetRA ( )

inline

Definition at line 194 of file G4hhElastic.hh.

{ return fRA;};

References fRA.

Referenced by SetParametersCMS().

GetRB()

G4double G4hhElastic::GetRB ( )

inline

Definition at line 198 of file G4hhElastic.hh.

{ return fRB;};

References fRB.

Referenced by SetParametersCMS().

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

GetRg()

G4double G4hhElastic::GetRg ( )

inline

Definition at line 199 of file G4hhElastic.hh.

{ return fRg;};

References fRg.

GetRG()

G4double G4hhElastic::GetRG ( )

inline

Definition at line 200 of file G4hhElastic.hh.

{ return fRG;};

References fRG.

GetRhoReIm()

G4double G4hhElastic::GetRhoReIm ( )

inline

Definition at line 175 of file G4hhElastic.hh.

{ return fRhoReIm;};

References fRhoReIm.

GetRq()

G4double G4hhElastic::GetRq ( )

inline

Definition at line 195 of file G4hhElastic.hh.

{ return fRq;};

References fRq.

GetRQ()

G4double G4hhElastic::GetRQ ( )

inline

Definition at line 196 of file G4hhElastic.hh.

{ return fRQ;};

References fRQ.

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

GetSpp()

G4double G4hhElastic::GetSpp ( )

inline

Definition at line 166 of file G4hhElastic.hh.

{return fSpp;};

References fSpp.

GetTransfer()

G4double G4hhElastic::GetTransfer	(	G4int	iMomentum,
		G4int	iTransfer,
		G4double	position
	)

Definition at line 629 of file G4hhElastic.cc.

{
  G4double x1, x2, y1, y2, randTransfer, delta, mean, epsilon = 1.e-6;
 
  if( iTransfer == 0 )
  {
    randTransfer = (*fTableT)(iTkin)->GetLowEdgeEnergy(iTransfer);
    // iTransfer++;
  }
  else
  {
    if ( iTransfer >= G4int((*fTableT)(iTkin)->GetVectorLength()) )
    {
      iTransfer = (*fTableT)(iTkin)->GetVectorLength() - 1;
    }
    y1 = (*(*fTableT)(iTkin))(iTransfer-1);
    y2 = (*(*fTableT)(iTkin))(iTransfer);
 
    x1 = (*fTableT)(iTkin)->GetLowEdgeEnergy(iTransfer-1);
    x2 = (*fTableT)(iTkin)->GetLowEdgeEnergy(iTransfer);
 
    delta = y2 - y1;
    mean  = y2 + y1;
 
    if ( x1 == x2 ) randTransfer = x2;
    else
    {
      // if ( y1 == y2 ) 
      if ( delta < epsilon*mean ) 
           randTransfer = x1 + ( x2 - x1 )*G4UniformRand();
      else randTransfer = x1 + ( position - y1 )*( x2 - x1 )/delta; // ( y2 - y1 );
    }
  }
  return randTransfer;
}

References epsilon(), fTableT, and G4UniformRand.

Referenced by SampleBisectionalT(), SampleInvariantT(), and SampleTest().

GetVerboseLevel()

G4int G4HadronicInteraction::GetVerboseLevel ( ) const

inlineinherited

Definition at line 109 of file G4HadronicInteraction.hh.

  { return verboseLevel; }

References G4HadronicInteraction::verboseLevel.

Initialise()

void G4hhElastic::Initialise

(

)

Definition at line 230 of file G4hhElastic.cc.

{
  // pp,pn
 
  fProjectile = G4Proton::Proton();
  BuildTableT(fTarget, fProjectile);
  fBankT.push_back(fTableT); // 0
 
  // pi+-p
  
  fProjectile = G4PionPlus::PionPlus();
  BuildTableT(fTarget, fProjectile);
  fBankT.push_back(fTableT);  // 1
  //K+-p
  fProjectile = G4KaonPlus::KaonPlus();
  BuildTableT(fTarget, fProjectile);
  fBankT.push_back(fTableT);  // 2
  
}

References BuildTableT(), fBankT, fProjectile, fTableT, fTarget, G4KaonPlus::KaonPlus(), G4PionPlus::PionPlus(), and G4Proton::Proton().

Referenced by G4hhElastic().

InitialiseModel()

void G4HadronicInteraction::InitialiseModel ( )

virtualinherited

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

Definition at line 59 of file G4HadronicInteraction.cc.

{}

IsApplicable()

G4bool G4hhElastic::IsApplicable ( const G4HadProjectile &  projectile, G4Nucleus &  nucleus  )

inlinevirtual

Reimplemented from G4HadronicInteraction.

Definition at line 256 of file G4hhElastic.hh.

{
  if( ( projectile.GetDefinition() == G4Proton::Proton() ||
        projectile.GetDefinition() == G4Neutron::Neutron() ||
        projectile.GetDefinition() == G4PionPlus::PionPlus() ||
        projectile.GetDefinition() == G4PionMinus::PionMinus() ||
        projectile.GetDefinition() == G4KaonPlus::KaonPlus() ||
        projectile.GetDefinition() == G4KaonMinus::KaonMinus() ) &&
 
        nucleus.GetZ_asInt() < 2 ) return true;
  else                              return false;
}

References G4HadProjectile::GetDefinition(), G4Nucleus::GetZ_asInt(), G4KaonMinus::KaonMinus(), G4KaonPlus::KaonPlus(), G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), and G4Proton::Proton().

IsBlocked() [1/3]

G4bool G4HadronicInteraction::IsBlocked ( ) const

inlineprotectedinherited

Definition at line 169 of file G4HadronicInteraction.hh.

{ return isBlocked;}

References G4HadronicInteraction::isBlocked.

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

IsBlocked() [2/3]

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

inherited

Definition at line 202 of file G4HadronicInteraction.cc.

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

References G4HadronicInteraction::theBlockedListElements.

IsBlocked() [3/3]

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

inherited

Definition at line 193 of file G4HadronicInteraction.cc.

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

References G4HadronicInteraction::theBlockedList.

LowestEnergyLimit()

G4double G4HadronElastic::LowestEnergyLimit ( ) const

inlineinherited

Definition at line 96 of file G4HadronElastic.hh.

{
  return lowestEnergyLimit;
}

References G4HadronElastic::lowestEnergyLimit.

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

ModelDescription()

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

overridevirtualinherited

Reimplemented from G4HadronicInteraction.

Reimplemented in G4NeutrinoElectronNcModel, and G4NeutronElectronElModel.

Definition at line 70 of file G4HadronElastic.cc.

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

operator!=()

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

deleteinherited

operator==()

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

deleteinherited

Phi13()

G4complex G4hhElastic::Phi13 ( )

inline

Definition at line 520 of file G4hhElastic.hh.

{
  G4double re = (fRq*fRq + fRg*fRg)/16.;
  G4complex result(re,0.);
  result += Pomeron();
  return result;
}

References fRg, fRq, and Pomeron().

Referenced by CalculateBqQ12(), CalculateBqQ123(), GetF1qQgG(), GetF2qQgG(), and GetF3qQgG().

Phi14()

G4complex G4hhElastic::Phi14 ( )

inline

Definition at line 530 of file G4hhElastic.hh.

{
  G4double re = (fRq*fRq + fRG*fRG)/16.;
  G4complex result(re,0.);
  result += Pomeron();
  return result;
}

References fRG, fRq, and Pomeron().

Referenced by CalculateBqQ12(), CalculateBqQ123(), CalculateBqQ13(), GetdsdtF13qQG(), GetF1qQgG(), GetF2qQgG(), and GetF3qQgG().

Phi23()

G4complex G4hhElastic::Phi23 ( )

inline

Definition at line 540 of file G4hhElastic.hh.

{
  G4double re = (fRQ*fRQ + fRg*fRg)/16.;
  G4complex result(re,0.);
  result += Pomeron();
  return result;
}

References fRg, fRQ, and Pomeron().

Referenced by CalculateBqQ12(), CalculateBqQ123(), GetF1qQgG(), GetF2qQgG(), and GetF3qQgG().

Phi24()

G4complex G4hhElastic::Phi24 ( )

inline

Definition at line 550 of file G4hhElastic.hh.

{
  G4double re = (fRQ*fRQ + fRG*fRG)/16.;
  G4complex result(re,0.);
  result += Pomeron();
  return result;
}

References fRG, fRQ, and Pomeron().

Referenced by CalculateBqQ12(), CalculateBqQ123(), CalculateBqQ13(), GetdsdtF13qQG(), GetF1qQgG(), GetF2qQgG(), and GetF3qQgG().

Pomeron()

G4complex G4hhElastic::Pomeron ( )

inline

Definition at line 509 of file G4hhElastic.hh.

{
  G4double re, im;
 
  re = fAlphaP*G4Log(fSpp/fSo);
  im = -0.5*fAlphaP*fImCof*CLHEP::pi;
  return G4complex(re,im);
} 

References fAlphaP, fImCof, fSo, fSpp, G4Log(), and CLHEP::pi.

Referenced by Phi13(), Phi14(), Phi23(), and Phi24().

SampleBisectionalT()

G4double G4hhElastic::SampleBisectionalT	(	const G4ParticleDefinition *	p,
		G4double	plab
	)

Definition at line 439 of file G4hhElastic.cc.

{
  G4int iTkin, iTransfer;
  G4double t,  position, m1 = aParticle->GetPDGMass();
  G4double Tkin = std::sqrt(m1*m1+p*p) - m1;
 
  if( aParticle == G4Proton::Proton() || aParticle == G4Neutron::Neutron() )
  {
    fTableT = fBankT[0];
  }
  if( aParticle == G4PionPlus::PionPlus() || aParticle == G4PionMinus::PionMinus() )
  {
    fTableT = fBankT[1];
  }
  if( aParticle == G4KaonPlus::KaonPlus() || aParticle == G4KaonMinus::KaonMinus() )
  {
    fTableT = fBankT[2];
  }
  G4double delta    = std::abs(Tkin - fOldTkin)/(Tkin + fOldTkin);
  G4double deltaMax = 1.e-2;
 
  if ( delta < deltaMax ) iTkin = fInTkin; 
  else
  {  
    for( iTkin = 0; iTkin < fEnergyBin; iTkin++ )
    {
      if( Tkin < fEnergyVector->GetLowEdgeEnergy(iTkin) ) break;
    }
  }
  if ( iTkin >= fEnergyBin ) iTkin = fEnergyBin-1;   // Tkin is more then theMaxEnergy
  if ( iTkin < 0 )           iTkin = 0; // against negative index, Tkin < theMinEnergy
 
  fOldTkin = Tkin;
  fInTkin = iTkin;
 
  if (iTkin == fEnergyBin -1 || iTkin == 0 )   // the table edges
  {
    position = (*(*fTableT)(iTkin))(0)*G4UniformRand();
 
    for(iTransfer = 0; iTransfer < fBinT-1; iTransfer++)
    {
      if( position >= (*(*fTableT)(iTkin))(iTransfer) ) break;
    }
    if (iTransfer >= fBinT-1) iTransfer = fBinT-2;
 
    t = GetTransfer(iTkin, iTransfer, position);
 
 
  }
  else  // Tkin inside between energy table edges
  {
    G4double rand = G4UniformRand();
    position = (*(*fTableT)(iTkin))(0)*rand;
 
    //
    // (*fTableT)(iTkin)->GetLowEdgeEnergy(fBinT-2);
    G4int sTransfer = 0, fTransfer =  fBinT - 2, dTransfer = fTransfer - sTransfer;
    G4double y2;
 
    for( iTransfer = 0; iTransfer < fBinT - 1; iTransfer++ )
    {
      // dTransfer %= 2;
      dTransfer /= 2;
      // dTransfer *= 0.5;
      y2 = (*(*fTableT)(iTkin))( sTransfer + dTransfer );
 
      if( y2 > position ) sTransfer += dTransfer;
 
      // if( dTransfer <= 1 ) break;
      if( dTransfer < 1 ) break;
    }
    t = (*fTableT)(iTkin)->GetLowEdgeEnergy(sTransfer); // +(-0.5+rand)*(*fTableT)(iTkin)->GetLowEdgeEnergy(3);
  }
  return t;
}

References fBankT, fBinT, fEnergyBin, fInTkin, fOldTkin, fTableT, G4UniformRand, G4ParticleDefinition::GetPDGMass(), GetTransfer(), G4KaonMinus::KaonMinus(), G4KaonPlus::KaonPlus(), G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), position, and G4Proton::Proton().

SampleInvariantT()

G4double G4hhElastic::SampleInvariantT ( const G4ParticleDefinition *  p, G4double  plab, G4int  , G4int    )

virtual

Reimplemented from G4HadronElastic.

Definition at line 322 of file G4hhElastic.cc.

{
  G4int iTkin, iTransfer;
  G4double t, t2, position, m1 = aParticle->GetPDGMass();
  G4double Tkin = std::sqrt(m1*m1+p*p) - m1;
 
  if( aParticle == G4Proton::Proton() || aParticle == G4Neutron::Neutron() )
  {
    fTableT = fBankT[0];
  }
  if( aParticle == G4PionPlus::PionPlus() || aParticle == G4PionMinus::PionMinus() )
  {
    fTableT = fBankT[1];
  }
  if( aParticle == G4KaonPlus::KaonPlus() || aParticle == G4KaonMinus::KaonMinus() )
  {
    fTableT = fBankT[2];
  }
 
  G4double delta    = std::abs(Tkin - fOldTkin)/(Tkin + fOldTkin);
  G4double deltaMax = 1.e-2;
 
  if ( delta < deltaMax ) iTkin = fInTkin; 
  else
  {  
    for( iTkin = 0; iTkin < fEnergyBin; iTkin++)
    {
      if( Tkin < fEnergyVector->GetLowEdgeEnergy(iTkin) ) break;
    }
  }
  if ( iTkin >= fEnergyBin ) iTkin = fEnergyBin-1;   // Tkin is more then theMaxEnergy
  if ( iTkin < 0 )           iTkin = 0; // against negative index, Tkin < theMinEnergy
 
  fOldTkin = Tkin;
  fInTkin = iTkin;
 
  if (iTkin == fEnergyBin -1 || iTkin == 0 )   // the table edges
  {
    position = (*(*fTableT)(iTkin))(0)*G4UniformRand();
 
    // G4cout<<"position = "<<position<<G4endl;
 
    for(iTransfer = 0; iTransfer < fBinT-1; iTransfer++)
    {
      if( position >= (*(*fTableT)(iTkin))(iTransfer) ) break;
    }
    if (iTransfer >= fBinT-1) iTransfer = fBinT-2;
 
    // G4cout<<"iTransfer = "<<iTransfer<<G4endl;
 
    t = GetTransfer(iTkin, iTransfer, position);
 
    // G4cout<<"t = "<<t<<G4endl;
  }
  else  // Tkin inside between energy table edges
  {
    // position = (*(*fTableT)(iTkin))(fBinT-2)*G4UniformRand();
    position = (*(*fTableT)(iTkin))(0)*G4UniformRand();
 
    // G4cout<<"position = "<<position<<G4endl;
 
    for(iTransfer = 0; iTransfer < fBinT-1; iTransfer++)
    {
      // if( position < (*(*fTableT)(iTkin))(iTransfer) ) break;
      if( position >= (*(*fTableT)(iTkin))(iTransfer) ) break;
    }
    if (iTransfer >= fBinT-1) iTransfer = fBinT-2;
 
    // G4cout<<"iTransfer = "<<iTransfer<<G4endl;
 
    t2  = GetTransfer(iTkin, iTransfer, position);
    return t2;
    /*
    G4double t1, E1, E2, W, W1, W2;
    // G4cout<<"t2 = "<<t2<<G4endl;
 
    E2 = fEnergyVector->GetLowEdgeEnergy(iTkin);
 
    // G4cout<<"E2 = "<<E2<<G4endl;
    
    iTkin--;
    
    // position = (*(*fTableT)(iTkin))(fBinT-2)*G4UniformRand();
 
    // G4cout<<"position = "<<position<<G4endl;
 
    for(iTransfer = 0; iTransfer < fBinT-1; iTransfer++)
    {
      // if( position < (*(*fTableT)(iTkin))(iTransfer) ) break;
      if( position >= (*(*fTableT)(iTkin))(iTransfer) ) break;
    }
    if (iTransfer >= fBinT-1) iTransfer = fBinT-2;
    
    t1  = GetTransfer(iTkin, iTransfer, position);
 
    // G4cout<<"t1 = "<<t1<<G4endl;
 
    E1 = fEnergyVector->GetLowEdgeEnergy(iTkin);
 
    // G4cout<<"E1 = "<<E1<<G4endl;
 
    W  = 1.0/(E2 - E1);
    W1 = (E2 - Tkin)*W;
    W2 = (Tkin - E1)*W;
 
    t = W1*t1 + W2*t2;
    */
  }
  return t;
}

References fBankT, fBinT, fEnergyBin, fInTkin, fOldTkin, fTableT, G4UniformRand, G4ParticleDefinition::GetPDGMass(), GetTransfer(), G4KaonMinus::KaonMinus(), G4KaonPlus::KaonPlus(), G4Neutron::Neutron(), G4PionMinus::PionMinus(), G4PionPlus::PionPlus(), position, and G4Proton::Proton().

SampleTest()

G4double G4hhElastic::SampleTest

(

G4double

tMin

)

Definition at line 593 of file G4hhElastic.cc.

{
  G4int iTkin, iTransfer, iTmin;
  G4double t, position;
  // G4double qMin = std::sqrt(tMin);
 
  fTableT = fBankT[0];
  iTkin = 0;
 
  for(iTransfer = 0; iTransfer < fBinT-1; iTransfer++)
  {
    // if( qMin <= (*fTableT)(iTkin)->GetLowEdgeEnergy(iTransfer) ) break;
    if( tMin <= (*fTableT)(iTkin)->GetLowEdgeEnergy(iTransfer) ) break;
  }
  iTmin = iTransfer-1;
  if(iTmin < 0 ) iTmin = 0;
 
  position = (*(*fTableT)(iTkin))(iTmin)*G4UniformRand();
 
  for( iTmin = 0; iTransfer < fBinT-1; iTransfer++)
  {
    if( position > (*(*fTableT)(iTkin))(iTransfer) ) break;
  }
  if (iTransfer >= fBinT-1) iTransfer = fBinT-2;
 
  t = GetTransfer(iTkin, iTransfer, position);
 
  return t;
}

References fBankT, fBinT, fTableT, G4UniformRand, GetTransfer(), and position.

SetAlphaP()

void G4hhElastic::SetAlphaP ( G4double  a )

inline

Definition at line 184 of file G4hhElastic.hh.

{fAlphaP = a;};

References fAlphaP.

SetBq()

void G4hhElastic::SetBq ( G4double  b )

inline

Definition at line 172 of file G4hhElastic.hh.

{fBq = b;};

References fBq.

Referenced by SetParametersCMS().

SetBQ()

void G4hhElastic::SetBQ ( G4double  b )

inline

Definition at line 173 of file G4hhElastic.hh.

{fBQ = b;};

References fBQ.

Referenced by SetParametersCMS().

SetBqQ()

void G4hhElastic::SetBqQ ( G4double  b )

inline

Definition at line 174 of file G4hhElastic.hh.

{fBqQ = b;};

References fBqQ.

SetCofF2()

void G4hhElastic::SetCofF2 ( G4double  f )

inline

Definition at line 189 of file G4hhElastic.hh.

{fCofF2 = f;};

References fCofF2.

Referenced by SetParametersCMS().

SetCofF3()

void G4hhElastic::SetCofF3 ( G4double  f )

inline

Definition at line 190 of file G4hhElastic.hh.

{fCofF3 = f;};

References fCofF3.

Referenced by SetParametersCMS().

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

SetEta()

void G4hhElastic::SetEta ( G4double  E )

inline

Definition at line 188 of file G4hhElastic.hh.

{fEta = E;};

References fEta.

Referenced by SetParametersCMS().

SetImCof()

void G4hhElastic::SetImCof ( G4double  a )

inline

Definition at line 185 of file G4hhElastic.hh.

{fImCof = a;};

References fImCof.

Referenced by SetParametersCMS().

SetLambda()

void G4hhElastic::SetLambda ( G4double  L )

inline

Definition at line 187 of file G4hhElastic.hh.

{fLambda = L;};

References fLambda, and L.

Referenced by SetParametersCMS().

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

SetParameters()

void G4hhElastic::SetParameters ( )

inline

Definition at line 271 of file G4hhElastic.hh.

{
  // masses
 
  fMq     = 0.36*CLHEP::GeV; // 0.441*GeV; // 0.36*GeV;
  fMQ     = 0.441*CLHEP::GeV;
  fMff2   = 0.26*CLHEP::GeV*CLHEP::GeV; // 0.25*GeV*GeV; // 0.5*GeV*GeV;
 
  fAlpha = 1./3.;
  fBeta  = 1. - fAlpha;
 
  fGamma = 1./2.; // 1./3.; // 
  fDelta = 1. - fGamma; // 1./2.;
 
  // radii and exp cof
 
  fRA       = 6.5/CLHEP::GeV; // 7.3/GeV;  //  3.25/GeV; // 7./GeV; // 2./GeV; // 1./GeV;
  fRq       = 0.173*fRA; // 2.4/GeV;
  fRQ       = 0.316*fRA; // 1./GeV; // 2./GeV; // 1./GeV;
  fRB       = 6.5/CLHEP::GeV; // 7.3/GeV;  //  3.25/GeV; // 7./GeV; // 2./GeV; // 1./GeV;
  fRg       = 0.173*fRA; // 2.4/GeV;
  fRG       = 0.173*fRA; // 1./GeV; // 2./GeV; // 1./GeV;
 
  fAlphaP   = 0.15/CLHEP::GeV/CLHEP::GeV; // 0.15/GeV/GeV;
  fLambda   = 0.25*fRA*fRA; // 0.25
  fEta      = 0.25*fRB*fRB; // 0.25
  fImCof    = 6.5;
  fCofF2 = 1.;
  fCofF3 = 1.;
 
  fBq = 0.02; // 0.21; // 1./3.;
  fBQ = 1. + fBq - 2*std::sqrt(fBq); // 1 - fBq; // 2./3.;
  fBqQ = std::sqrt(fBq*fBQ);
 
  fLambdaFF = 1.5/CLHEP::GeV/CLHEP::GeV; // 0.15/GeV/GeV;
  fSo     = 1.*CLHEP::GeV*CLHEP::GeV;
  fQcof   = 0.009*CLHEP::GeV;
  fExpSlope = 19.9/CLHEP::GeV/CLHEP::GeV;
}

References fAlpha, fAlphaP, fBeta, fBq, fBQ, fBqQ, fCofF2, fCofF3, fDelta, fEta, fExpSlope, fGamma, fImCof, fLambda, fLambdaFF, fMff2, fMQ, fMq, fQcof, fRA, fRB, fRG, fRg, fRQ, fRq, fSo, and CLHEP::GeV.

Referenced by G4hhElastic().

SetParametersCMS()

void G4hhElastic::SetParametersCMS ( G4double  plab )

inline

Definition at line 316 of file G4hhElastic.hh.

{
  G4int i;
  G4double trMass = 900.*CLHEP::MeV, Tkin;
  G4double sl, sh, ds, rAl, rAh, drA, rBl, rBh, drB, bql, bqh, dbq, bQl, bQh, dbQ, cIl, cIh, dcI;
 
  Tkin = std::sqrt(fMassProj*fMassProj + plab*plab) - fMassProj;
 
  G4DynamicParticle*  theDynamicParticle = new G4DynamicParticle(fProjectile,
                                              G4ParticleMomentum(0.,0.,1.),
                                              Tkin);
  fSigmaTot = fHadrNuclXsc->GetHadronNucleonXscNS( theDynamicParticle, fTarget );
 
  delete theDynamicParticle;
 
  fSpp  = fMassProj*fMassProj + fMassTarg*fMassTarg + 2.*fMassTarg*std::sqrt(plab*plab + fMassProj*fMassProj);
  fPcms = std::sqrt( (fSpp - fMassSum2)*(fSpp - fMassDif2)/4./fSpp);
 
  G4double sCMS = std::sqrt(fSpp);
 
  if( fMassProj > trMass ) // p,n,pb on p
  {
    this->SetCofF2(1.);
    this->SetCofF3(1.);
    fGamma = 1./3.; // 1./3.; // 
    fDelta = 1. - fGamma; // 1./2.;
 
    if( sCMS <= theNuclNuclData[0][0]*CLHEP::GeV ) // low edge, as s=2.76754
    {
      this->SetRA(theNuclNuclData[0][1]/CLHEP::GeV,0.173,0.316);
      this->SetRB(theNuclNuclData[0][2]/CLHEP::GeV,0.173,0.316); 
 
      this->SetBq(theNuclNuclData[0][3]);
      this->SetBQ(theNuclNuclData[0][4]);
      this->SetImCof(theNuclNuclData[0][5]);
 
      this->SetLambda(0.25*this->GetRA()*this->GetRA());
      this->SetEta(0.25*this->GetRB()*this->GetRB());
    }
    else if( sCMS >= theNuclNuclData[17][0]*CLHEP::GeV ) // high edge, as s=7000 ???
    {
      this->SetRA(theNuclNuclData[17][1]/CLHEP::GeV,0.173,0.316);
      this->SetRB(theNuclNuclData[17][2]/CLHEP::GeV,0.173,0.316); 
 
      this->SetBq(theNuclNuclData[17][3]);
      this->SetBQ(theNuclNuclData[17][4]);
      this->SetImCof(theNuclNuclData[17][5]);
 
      this->SetLambda(0.25*this->GetRA()*this->GetRA());
      this->SetEta(0.25*this->GetRB()*this->GetRB());
    }
    else // in approximation between array points
    {
      for( i = 0; i < 19; i++ ) if( sCMS <= theNuclNuclData[i][0]*CLHEP::GeV ) break;
      if( i == 0 ) i++;
      if( i == 19 ) i--;
 
      sl = theNuclNuclData[i-1][0]*CLHEP::GeV;
      sh = theNuclNuclData[i][0]*CLHEP::GeV;
      ds = (sCMS - sl)/(sh - sl);
 
      rAl = theNuclNuclData[i-1][1]/CLHEP::GeV;
      rAh = theNuclNuclData[i][1]/CLHEP::GeV;
      drA = rAh - rAl;
 
      rBl = theNuclNuclData[i-1][2]/CLHEP::GeV;
      rBh = theNuclNuclData[i][2]/CLHEP::GeV;
      drB = rBh - rBl;
 
      bql = theNuclNuclData[i-1][3];
      bqh = theNuclNuclData[i][3];
      dbq = bqh - bql;
 
      bQl = theNuclNuclData[i-1][4];
      bQh = theNuclNuclData[i][4];
      dbQ = bQh - bQl;
 
      cIl = theNuclNuclData[i-1][5];
      cIh = theNuclNuclData[i][5];
      dcI = cIh - cIl;
 
      this->SetRA(rAl+drA*ds,0.173,0.316);
      this->SetRB(rBl+drB*ds,0.173,0.316); 
 
      this->SetBq(bql+dbq*ds);
      this->SetBQ(bQl+dbQ*ds);
      this->SetImCof(cIl+dcI*ds);
 
      this->SetLambda(0.25*this->GetRA()*this->GetRA());
      this->SetEta(0.25*this->GetRB()*this->GetRB());
    }
  }
  else // pi, K
  {
    this->SetCofF2(1.);
    this->SetCofF3(-1.);
    fGamma = 1./2.; // 1./3.; // 
    fDelta = 1. - fGamma; // 1./2.;
  
    if( sCMS <= thePiKaNuclData[0][0]*CLHEP::GeV ) // low edge, as s=2.76754
    {
      this->SetRA(thePiKaNuclData[0][1]/CLHEP::GeV,0.173,0.316);
      this->SetRB(thePiKaNuclData[0][2]/CLHEP::GeV,0.173,0.173); 
 
      this->SetBq(thePiKaNuclData[0][3]);
      this->SetBQ(thePiKaNuclData[0][4]);
      this->SetImCof(thePiKaNuclData[0][5]);
 
      this->SetLambda(0.25*this->GetRA()*this->GetRA());
      this->SetEta(this->GetRB()*this->GetRB()/6.);
    }
    else if( sCMS >= thePiKaNuclData[7][0]*CLHEP::GeV ) // high edge, as s=7000 ???
    {
      this->SetRA(thePiKaNuclData[7][1]/CLHEP::GeV,0.173,0.316);
      this->SetRB(thePiKaNuclData[7][2]/CLHEP::GeV,0.173,0.173); 
 
      this->SetBq(thePiKaNuclData[7][3]);
      this->SetBQ(thePiKaNuclData[7][4]);
      this->SetImCof(thePiKaNuclData[7][5]);
 
      this->SetLambda(0.25*this->GetRA()*this->GetRA());
      this->SetEta(this->GetRB()*this->GetRB()/6.);
    }
    else // in approximation between array points
    {
      for( i = 0; i < 8; i++ ) if( sCMS <= thePiKaNuclData[i][0]*CLHEP::GeV ) break;
      if( i == 0 ) i++;
      if( i == 8 ) i--;
 
      sl = thePiKaNuclData[i-1][0]*CLHEP::GeV;
      sh = thePiKaNuclData[i][0]*CLHEP::GeV;
      ds = (sCMS - sl)/(sh - sl);
 
      rAl = thePiKaNuclData[i-1][1]/CLHEP::GeV;
      rAh = thePiKaNuclData[i][1]/CLHEP::GeV;
      drA = rAh - rAl;
 
      rBl = thePiKaNuclData[i-1][2]/CLHEP::GeV;
      rBh = thePiKaNuclData[i][2]/CLHEP::GeV;
      drB = rBh - rBl;
 
      bql = thePiKaNuclData[i-1][3];
      bqh = thePiKaNuclData[i][3];
      dbq = bqh - bql;
 
      bQl = thePiKaNuclData[i-1][4];
      bQh = thePiKaNuclData[i][4];
      dbQ = bQh - bQl;
 
      cIl = thePiKaNuclData[i-1][5];
      cIh = thePiKaNuclData[i][5];
      dcI = cIh - cIl;
 
      this->SetRA(rAl+drA*ds,0.173,0.316);
      this->SetRB(rBl+drB*ds,0.173,0.173); 
 
      this->SetBq(bql+dbq*ds);
      this->SetBQ(bQl+dbQ*ds);
      this->SetImCof(cIl+dcI*ds);
 
      this->SetLambda(0.25*this->GetRA()*this->GetRA());
      this->SetEta(this->GetRB()*this->GetRB()/6.);
    }
   }
  return;
}

References anonymous_namespace{G4QuasiElRatios.cc}::ds, fDelta, fGamma, fHadrNuclXsc, fMassDif2, fMassProj, fMassSum2, fMassTarg, fPcms, fProjectile, fSigmaTot, fSpp, fTarget, G4HadronNucleonXsc::GetHadronNucleonXscNS(), GetRA(), GetRB(), CLHEP::GeV, CLHEP::MeV, SetBq(), SetBQ(), SetCofF2(), SetCofF3(), SetEta(), SetImCof(), SetLambda(), SetRA(), SetRB(), theNuclNuclData, and thePiKaNuclData.

Referenced by BuildTableT(), and G4hhElastic().

SetRA()

void G4hhElastic::SetRA ( G4double  rn, G4double  pq, G4double  pQ  )

inline

Definition at line 487 of file G4hhElastic.hh.

{
  fRA = rA;
  fRq = fRA*pq;
  fRQ = fRA*pQ;
}

References fRA, fRQ, and fRq.

Referenced by SetParametersCMS().

SetRB()

void G4hhElastic::SetRB ( G4double  rn, G4double  pq, G4double  pQ  )

inline

Definition at line 498 of file G4hhElastic.hh.

{
  fRB = rB;
  fRg = fRB*pg;
  fRG = fRB*pG;
}

References fRB, fRG, fRg, and demo::pg.

Referenced by SetParametersCMS().

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

SetSigmaTot()

void G4hhElastic::SetSigmaTot ( G4double  stot )

inline

Definition at line 164 of file G4hhElastic.hh.

{fSigmaTot = stot;};

References fSigmaTot.

SetSpp()

void G4hhElastic::SetSpp ( G4double  spp )

inline

Definition at line 165 of file G4hhElastic.hh.

{fSpp = spp;};

References fSpp.

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

Field Documentation

epCheckLevels

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

privateinherited

Definition at line 198 of file G4HadronicInteraction.hh.

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

fAlpha

G4double G4hhElastic::fAlpha

private

Definition at line 132 of file G4hhElastic.hh.

Referenced by CalculateBqQ13(), GetdsdtF13qQG(), GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), and SetParameters().

fAlphaP

G4double G4hhElastic::fAlphaP

private

Definition at line 141 of file G4hhElastic.hh.

Referenced by GetAqq(), GetAQQ(), Pomeron(), SetAlphaP(), and SetParameters().

fBankT

std::vector<G4PhysicsTable*> G4hhElastic::fBankT

private

Definition at line 113 of file G4hhElastic.hh.

Referenced by BuildTableTest(), Initialise(), SampleBisectionalT(), SampleInvariantT(), SampleTest(), and ~G4hhElastic().

fBeta

G4double G4hhElastic::fBeta

private

Definition at line 133 of file G4hhElastic.hh.

Referenced by GetdsdtF13qQG(), GetF1qQgG(), GetF3qQgG(), and SetParameters().

fBinT

G4int G4hhElastic::fBinT

private

Definition at line 109 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), SampleBisectionalT(), SampleInvariantT(), and SampleTest().

fBq

G4double G4hhElastic::fBq

private

Definition at line 153 of file G4hhElastic.hh.

Referenced by CalculateBQ(), CalculateBqQ12(), CalculateBqQ123(), GetBq(), GetF1(), GetF1qQgG(), GetF2(), GetF2qQgG(), GetF3(), GetF3qQgG(), SetBq(), and SetParameters().

fBQ

G4double G4hhElastic::fBQ

private

Definition at line 154 of file G4hhElastic.hh.

Referenced by CalculateBQ(), CalculateBqQ12(), CalculateBqQ123(), CalculateBqQ13(), GetBQ(), GetdsdtF13qQG(), GetF1(), GetF1qQgG(), GetF2(), GetF2qQgG(), GetF3(), GetF3qQgG(), SetBQ(), and SetParameters().

fBqQ

G4double G4hhElastic::fBqQ

private

Definition at line 155 of file G4hhElastic.hh.

Referenced by CalculateBqQ13(), GetBqQ(), GetdsdtF13qQG(), SetBqQ(), and SetParameters().

fCofF2

G4double G4hhElastic::fCofF2

private

Definition at line 146 of file G4hhElastic.hh.

Referenced by CalculateBqQ123(), GetCofF2(), GetCofS1(), GetdsdtF123(), GetdsdtF123qQgG(), GetExpRatioF123(), SetCofF2(), and SetParameters().

fCofF3

G4double G4hhElastic::fCofF3

private

Definition at line 147 of file G4hhElastic.hh.

Referenced by CalculateBqQ123(), GetCofF3(), GetCofS2(), GetCofS3(), GetdsdtF123(), GetdsdtF123qQgG(), GetExpRatioF123(), SetCofF3(), and SetParameters().

fDelta

G4double G4hhElastic::fDelta

private

Definition at line 139 of file G4hhElastic.hh.

Referenced by GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), SetParameters(), and SetParametersCMS().

fEnergyBin

G4int G4hhElastic::fEnergyBin

private

Definition at line 108 of file G4hhElastic.hh.

Referenced by BuildTableT(), G4hhElastic(), SampleBisectionalT(), and SampleInvariantT().

fEnergyVector

G4PhysicsLogVector* G4hhElastic::fEnergyVector

private

Definition at line 111 of file G4hhElastic.hh.

Referenced by BuildTableT(), G4hhElastic(), and ~G4hhElastic().

fEta

G4double G4hhElastic::fEta

private

Definition at line 144 of file G4hhElastic.hh.

Referenced by CalculateBqQ12(), CalculateBqQ123(), GetdsdtF13qQG(), GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), SetEta(), and SetParameters().

fExpSlope

G4double G4hhElastic::fExpSlope

private

Definition at line 149 of file G4hhElastic.hh.

Referenced by GetExpRatioF123(), and SetParameters().

fGamma

G4double G4hhElastic::fGamma

private

Definition at line 138 of file G4hhElastic.hh.

Referenced by GetdsdtF13qQG(), GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), SetParameters(), and SetParametersCMS().

fHadrNuclXsc

G4HadronNucleonXsc* G4hhElastic::fHadrNuclXsc

private

Definition at line 247 of file G4hhElastic.hh.

Referenced by G4hhElastic(), SetParametersCMS(), and ~G4hhElastic().

fImCof

G4double G4hhElastic::fImCof

private

Definition at line 145 of file G4hhElastic.hh.

Referenced by GetAqq(), GetAQQ(), GetImCof(), Pomeron(), SetImCof(), and SetParameters().

fInTkin

G4int G4hhElastic::fInTkin

private

Definition at line 243 of file G4hhElastic.hh.

Referenced by G4hhElastic(), SampleBisectionalT(), and SampleInvariantT().

fLambda

G4double G4hhElastic::fLambda

private

Definition at line 143 of file G4hhElastic.hh.

Referenced by CalculateBqQ12(), CalculateBqQ123(), GetAqq(), GetAQQ(), GetCofS1(), GetCofS2(), GetCofS3(), GetdsdtF13qQG(), GetF1qQgG(), GetF2(), GetF2qQgG(), GetF3(), GetF3qQgG(), SetLambda(), and SetParameters().

fLambdaFF

G4double G4hhElastic::fLambdaFF

private

Definition at line 142 of file G4hhElastic.hh.

Referenced by SetParameters().

fMassDif2

G4double G4hhElastic::fMassDif2

private

Definition at line 126 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), GetdsdtF123qQgG(), GetdsdtF12qQgG(), GetdsdtF13qQG(), GetdsdtF1qQgG(), GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), and SetParametersCMS().

fMassProj

G4double G4hhElastic::fMassProj

private

Definition at line 124 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), and SetParametersCMS().

fMassSum2

G4double G4hhElastic::fMassSum2

private

Definition at line 125 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), GetdsdtF123qQgG(), GetdsdtF12qQgG(), GetdsdtF13qQG(), GetdsdtF1qQgG(), GetF1qQgG(), GetF2qQgG(), GetF3qQgG(), and SetParametersCMS().

fMassTarg

G4double G4hhElastic::fMassTarg

private

Definition at line 123 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), and SetParametersCMS().

fMff2

G4double G4hhElastic::fMff2

private

Definition at line 119 of file G4hhElastic.hh.

Referenced by SetParameters().

fMQ

G4double G4hhElastic::fMQ

private

Definition at line 120 of file G4hhElastic.hh.

Referenced by SetParameters().

fMq

G4double G4hhElastic::fMq

private

Definition at line 121 of file G4hhElastic.hh.

Referenced by SetParameters().

fOldTkin

G4double G4hhElastic::fOldTkin

private

Definition at line 244 of file G4hhElastic.hh.

Referenced by G4hhElastic(), SampleBisectionalT(), and SampleInvariantT().

fOptRatio

G4double G4hhElastic::fOptRatio

private

Definition at line 156 of file G4hhElastic.hh.

Referenced by CalculateBQ(), CalculateBqQ123(), G4hhElastic(), and GetOpticalRatio().

fPcms

G4double G4hhElastic::fPcms

private

Definition at line 158 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), and SetParametersCMS().

fProjectile

G4ParticleDefinition* G4hhElastic::fProjectile

private

Definition at line 95 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), Initialise(), and SetParametersCMS().

fQcof

G4double G4hhElastic::fQcof

private

Definition at line 159 of file G4hhElastic.hh.

Referenced by SetParameters().

fRA

G4double G4hhElastic::fRA

private

Definition at line 129 of file G4hhElastic.hh.

Referenced by GetRA(), SetParameters(), and SetRA().

fRB

G4double G4hhElastic::fRB

private

Definition at line 135 of file G4hhElastic.hh.

Referenced by GetRB(), SetParameters(), and SetRB().

fRG

G4double G4hhElastic::fRG

private

Definition at line 136 of file G4hhElastic.hh.

Referenced by GetRG(), Phi14(), Phi24(), SetParameters(), and SetRB().

fRg

G4double G4hhElastic::fRg

private

Definition at line 137 of file G4hhElastic.hh.

Referenced by GetRg(), Phi13(), Phi23(), SetParameters(), and SetRB().

fRhoReIm

G4double G4hhElastic::fRhoReIm

private

Definition at line 148 of file G4hhElastic.hh.

Referenced by G4hhElastic(), GetExpRatioF123(), and GetRhoReIm().

fRQ

G4double G4hhElastic::fRQ

private

Definition at line 130 of file G4hhElastic.hh.

Referenced by GetAQQ(), GetRQ(), Phi23(), Phi24(), SetParameters(), and SetRA().

fRq

G4double G4hhElastic::fRq

private

Definition at line 131 of file G4hhElastic.hh.

Referenced by GetAqq(), GetRq(), Phi13(), Phi14(), SetParameters(), and SetRA().

fSigmaTot

G4double G4hhElastic::fSigmaTot

private

Definition at line 152 of file G4hhElastic.hh.

Referenced by CalculateBqQ12(), CalculateBqQ123(), CalculateBqQ13(), G4hhElastic(), GetCofS1(), GetCofS2(), GetCofS3(), GetdsdtF13qQG(), GetF1(), GetF1qQgG(), GetF2(), GetF2qQgG(), GetF3(), GetF3qQgG(), SetParametersCMS(), and SetSigmaTot().

fSo

G4double G4hhElastic::fSo

private

Definition at line 150 of file G4hhElastic.hh.

Referenced by GetAqq(), GetAQQ(), Pomeron(), and SetParameters().

fSpp

G4double G4hhElastic::fSpp

private

Definition at line 157 of file G4hhElastic.hh.

Referenced by BuildTableTest(), G4hhElastic(), GetAqq(), GetAQQ(), GetdsdtF1(), GetdsdtF12(), GetdsdtF123(), GetdsdtF123qQgG(), GetdsdtF12qQgG(), GetdsdtF13qQG(), GetdsdtF1qQgG(), GetExpRatioF123(), GetF1(), GetF1qQgG(), GetF2(), GetF2qQgG(), GetF3(), GetF3qQgG(), GetSpp(), Pomeron(), SetParametersCMS(), and SetSpp().

fTableT

G4PhysicsTable* G4hhElastic::fTableT

private

Definition at line 112 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), GetTransfer(), Initialise(), SampleBisectionalT(), SampleInvariantT(), SampleTest(), and ~G4hhElastic().

fTarget

G4ParticleDefinition* G4hhElastic::fTarget

private

Definition at line 94 of file G4hhElastic.hh.

Referenced by BuildTableT(), BuildTableTest(), G4hhElastic(), Initialise(), and SetParametersCMS().

isBlocked

G4bool G4HadronicInteraction::isBlocked

protectedinherited

Definition at line 188 of file G4HadronicInteraction.hh.

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

lowEnergyLimitHE

G4double G4hhElastic::lowEnergyLimitHE

private

Definition at line 103 of file G4hhElastic.hh.

Referenced by G4hhElastic().

lowEnergyLimitQ

G4double G4hhElastic::lowEnergyLimitQ

private

Definition at line 104 of file G4hhElastic.hh.

Referenced by G4hhElastic().

lowEnergyRecoilLimit

G4double G4hhElastic::lowEnergyRecoilLimit

private

Definition at line 102 of file G4hhElastic.hh.

Referenced by G4hhElastic().

lowestEnergyLimit

G4double G4hhElastic::lowestEnergyLimit

private

Definition at line 105 of file G4hhElastic.hh.

Referenced by G4hhElastic().

nwarn

G4int G4HadronElastic::nwarn

privateinherited

Definition at line 88 of file G4HadronElastic.hh.

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

plabLowLimit

G4double G4hhElastic::plabLowLimit

private

Definition at line 106 of file G4hhElastic.hh.

Referenced by G4hhElastic().

pLocalTmax

G4double G4HadronElastic::pLocalTmax

protectedinherited

Definition at line 77 of file G4HadronElastic.hh.

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

recoilEnergyThreshold

G4double G4HadronicInteraction::recoilEnergyThreshold

privateinherited

Definition at line 194 of file G4HadronicInteraction.hh.

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

registry

G4HadronicInteractionRegistry* G4HadronicInteraction::registry

privateinherited

Definition at line 192 of file G4HadronicInteraction.hh.

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

secID

G4int G4HadronElastic::secID

protectedinherited

Definition at line 78 of file G4HadronElastic.hh.

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

theAlpha

G4ParticleDefinition* 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(), 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(), 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* G4hhElastic::theNeutron

private

Definition at line 97 of file G4hhElastic.hh.

Referenced by G4hhElastic().

theNuclNuclData

const G4double G4hhElastic::theNuclNuclData

staticprivate

Initial value:

= 
{
  
 
  { 2.76754,    4.8,    4.8,    0.05,   0.742441,   10.5 }, 
  { 3.07744,    5.4,    5.4,    0.02,   0.83818,    6.5  }, 
  { 3.36305,    5.2,    5.2,    0.02,   0.838893,   7.5  }, 
  { 4.32941,    6,  6,  0.03,   0.769389,   7.5  }, 
  { 4.62126,    6,  6,  0.03,   0.770111,   6.5  }, 
 
  { 5.47416,    4.5,    4.5,    0.03,   0.813185,   7.5  }, 
  { 6.15088,    6.5,    6.5,    0.02,   0.799539,   6.5  }, 
  { 6.77474,    5.2,    5.2,    0.03,   0.784901,   7.5  }, 
  { 9.77775,    7,  7,  0.03,   0.742531,   6.5  }, 
                
  { 10.4728,    5.2,    5.2,    0.03,   0.780439,   7.5  }, 
 
  { 13.7631,    7,  7,  0.008,  0.8664,         5.0  }, 
  { 19.4184,    6.8,    6.8,    0.009,  0.861337,   2.5  }, 
  { 23.5,   6.8,    6.8,    0.007,  0.878112,   1.5  }, 
                
  { 24.1362,    7.2,    7.2,    0.008,  0.864745,   5.5  },
  { 52.8,   6.8,    6.8,    0.008,  0.871929,   1.5  }, 
 
  { 546,        7.4,    7.4,    0.013,  0.845877,   5.5  }, 
  { 1960,   7.8,    7.8,    0.022,  0.809062,   7.5  }, 
  { 7000,   8,  8,  0.024,  0.820441,   5.5  },  
  { 13000,  8.5,    8.5,    0.03,   0.796721,   10.5  }  
 
}

Definition at line 245 of file G4hhElastic.hh.

Referenced by SetParametersCMS().

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

thePiKaNuclData

const G4double G4hhElastic::thePiKaNuclData

staticprivate

Initial value:

= 
{
  
 
  { 2.5627,     3.8,    3.3,    0.22,   0.222,          1.5 }, 
  { 2.93928,    4.3,    3.8,    0.2,    0.250601,       1.3 }, 
  { 3.22326,    4.8,    4.3,    0.13,   0.32751,    2.5 }, 
  { 7.80704,    5.5,    5,  0.13,   0.340631,   2.5 }, 
  { 9.7328, 5,  4.5,    0.05,   0.416319,   5.5  }, 
 
  { 13.7315,    5.3,    4.8,    0.05,   0.418426,   5.5  }, 
  { 16.6359,    6.3,    5.8,    0.05,   0.423817,   5.5  }, 
  { 19.3961,    5,  4.5,    0.05,   0.413477,   3.5  } 
 
}

Definition at line 246 of file G4hhElastic.hh.

Referenced by SetParametersCMS().

thePionMinus

G4ParticleDefinition* G4hhElastic::thePionMinus

private

Definition at line 99 of file G4hhElastic.hh.

Referenced by G4hhElastic().

thePionPlus

G4ParticleDefinition* G4hhElastic::thePionPlus

private

Definition at line 98 of file G4hhElastic.hh.

Referenced by G4hhElastic().

theProton

G4ParticleDefinition* G4hhElastic::theProton

private

Definition at line 96 of file G4hhElastic.hh.

Referenced by G4hhElastic().

verboseLevel

G4int G4HadronicInteraction::verboseLevel

protectedinherited

Definition at line 177 of file G4HadronicInteraction.hh.

Referenced by G4WilsonAbrasionModel::ApplyYourself(), G4EMDissociation::ApplyYourself(), G4LFission::ApplyYourself(), G4MuMinusCapturePrecompound::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4LowEGammaNuclearModel::ApplyYourself(), G4ChargeExchange::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4LEHadronProtonElastic::ApplyYourself(), G4LEnp::ApplyYourself(), G4LEpp::ApplyYourself(), G4CascadeInterface::ApplyYourself(), G4CascadeInterface::checkFinalResult(), G4CascadeInterface::copyOutputToHadronicResult(), G4CascadeInterface::copyOutputToReactionProducts(), G4LENDModel::create_used_target_map(), G4CascadeInterface::createBullet(), G4CascadeInterface::createTarget(), G4ElasticHadrNucleusHE::DefineHadronValues(), G4ElasticHadrNucleusHE::FillData(), G4ElasticHadrNucleusHE::FillFq2(), G4DiffuseElastic::G4DiffuseElastic(), G4DiffuseElasticV2::G4DiffuseElasticV2(), G4ElasticHadrNucleusHE::G4ElasticHadrNucleusHE(), G4EMDissociation::G4EMDissociation(), G4hhElastic(), G4NuclNuclDiffuseElastic::G4NuclNuclDiffuseElastic(), G4WilsonAbrasionModel::G4WilsonAbrasionModel(), G4ElasticHadrNucleusHE::GetFt(), G4ElasticHadrNucleusHE::GetLightFq2(), G4ElasticHadrNucleusHE::GetQ2_2(), G4HadronicInteraction::GetVerboseLevel(), G4ElasticHadrNucleusHE::HadronNucleusQ2_2(), G4ElasticHadrNucleusHE::HadronProtonQ2(), G4LFission::init(), G4DiffuseElastic::Initialise(), G4DiffuseElasticV2::Initialise(), G4NuclNuclDiffuseElastic::Initialise(), G4DiffuseElastic::InitialiseOnFly(), G4DiffuseElasticV2::InitialiseOnFly(), G4NuclNuclDiffuseElastic::InitialiseOnFly(), G4CascadeInterface::makeDynamicParticle(), G4CascadeInterface::NoInteraction(), G4CascadeInterface::Propagate(), G4ElasticHadrNucleusHE::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/G4hhElastic.hh
• source/processes/hadronic/models/coherent_elastic/src/G4hhElastic.cc