#include <G4HESigmaZeroInelastic.hh>

Inheritance diagram for G4HESigmaZeroInelastic:

Public Member Functions
	G4HESigmaZeroInelastic ()

	~G4HESigmaZeroInelastic ()

virtual void	ModelDescription (std::ostream &) const

void	SetMaxNumberOfSecondaries (G4int maxnumber)

void	SetVerboseLevel (G4int verbose)

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

G4int	GetNumberOfSecondaries ()

Public Member Functions inherited from G4HEInelastic
	G4HEInelastic (const G4String &modelName="HEInelastic")

	~G4HEInelastic ()

void	SetMaxNumberOfSecondaries (const G4int maxnumber)

void	SetVerboseLevel (const G4int level)

void	ForceEnergyConservation (G4bool energyConservation)

G4bool	EnergyConservation (void)

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

G4double	Amin (G4double a, G4double b)

G4double	Amax (G4double a, G4double b)

G4int	Imin (G4int a, G4int b)

G4int	Imax (G4int a, G4int b)

void	FillParticleChange (G4HEVector pv[], G4int aVecLength)

G4double	pmltpc (G4int np, G4int nm, G4int nz, G4int n, G4double b, G4double c)

G4int	Factorial (G4int n)

G4double	NuclearInelasticity (G4double incidentKineticEnergy, G4double atomicWeight, G4double atomicNumber)

G4double	NuclearExcitation (G4double incidentKineticEnergy, G4double atomicWeight, G4double atomicNumber, G4double &excitationEnergyCascade, G4double &excitationEnergyEvaporation)

void	HighEnergyCascading (G4bool &successful, G4HEVector pv[], G4int &vecLen, G4double &excitationEnergyGNP, G4double &excitationEnergyDTA, const G4HEVector &incidentParticle, const G4HEVector &targetParticle, G4double atomicWeight, G4double atomicNumber)

void	HighEnergyClusterProduction (G4bool &successful, G4HEVector pv[], G4int &vecLen, G4double &excitationEnergyGNP, G4double &excitationEnergyDTA, const G4HEVector &incidentParticle, const G4HEVector &targetParticle, G4double atomicWeight, G4double atomicNumber)

void	TuningOfHighEnergyCascading (G4HEVector pv[], G4int &vecLen, const G4HEVector &incidentParticle, const G4HEVector &targetParticle, G4double atomicWeight, G4double atomicNumber)

void	MediumEnergyCascading (G4bool &successful, G4HEVector pv[], G4int &vecLen, G4double &excitationEnergyGNP, G4double &excitationEnergyDTA, const G4HEVector &incidentParticle, const G4HEVector &targetParticle, G4double atomicWeight, G4double atomicNumber)

void	MediumEnergyClusterProduction (G4bool &successful, G4HEVector pv[], G4int &vecLen, G4double &excitationEnergyGNP, G4double &excitationEnergyDTA, const G4HEVector &incidentParticle, const G4HEVector &targetParticle, G4double atomicWeight, G4double atomicNumber)

void	QuasiElasticScattering (G4bool &successful, G4HEVector pv[], G4int &vecLen, G4double &excitationEnergyGNP, G4double &excitationEnergyDTA, const G4HEVector &incidentParticle, const G4HEVector &targetParticle, G4double atomicWeight, G4double atomicNumber)

void	ElasticScattering (G4bool &successful, G4HEVector pv[], G4int &vecLen, const G4HEVector &incidentParticle, G4double atomicWeight, G4double atomicNumber)

G4int	rtmi (G4double *x, G4double xli, G4double xri, G4double eps, G4int iend, G4double aa, G4double bb, G4double cc, G4double dd, G4double rr)

G4double	fctcos (G4double t, G4double aa, G4double bb, G4double cc, G4double dd, G4double rr)

void	StrangeParticlePairProduction (const G4double availableEnergy, const G4double centerOfMassEnergy, G4HEVector pv[], G4int &vecLen, const G4HEVector &incidentParticle, const G4HEVector &targetParticle)

G4double	NBodyPhaseSpace (const G4double totalEnergy, const G4bool constantCrossSection, G4HEVector pv[], G4int &vecLen)

G4double	NBodyPhaseSpace (G4int npart, G4HEVector pv[], G4double wmax, G4double wfcn, G4int maxtrial, G4int ntrial)

G4double	gpdk (G4double a, G4double b, G4double c)

void	QuickSort (G4double arr[], const G4int lidx, const G4int ridx)

G4double	Alam (G4double a, G4double b, G4double c)

G4double	CalculatePhaseSpaceWeight (G4int npart)

G4double	normal (void)

G4double	GammaRand (G4double avalue)

G4double	Erlang (G4int mvalue)

G4int	Poisson (G4double x)

void	SetParticles (void)

Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")

virtual	~G4HadronicInteraction ()

virtual G4HadFinalState *	ApplyYourself (const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)=0

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

virtual G4bool	IsApplicable (const G4HadProjectile &, G4Nucleus &)

G4double	GetMinEnergy () const

G4double	GetMinEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMinEnergy (G4double anEnergy)

void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)

G4double	GetMaxEnergy () const

G4double	GetMaxEnergy (const G4Material aMaterial, const G4Element anElement) const

void	SetMaxEnergy (const G4double anEnergy)

void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)

void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)

const G4HadronicInteraction *	GetMyPointer () const

G4int	GetVerboseLevel () const

void	SetVerboseLevel (G4int value)

const G4String &	GetModelName () const

void	DeActivateFor (const G4Material *aMaterial)

void	ActivateFor (const G4Material *aMaterial)

void	DeActivateFor (const G4Element *anElement)

void	ActivateFor (const G4Element *anElement)

G4bool	IsBlocked (const G4Material *aMaterial) const

G4bool	IsBlocked (const G4Element *anElement) const

void	SetRecoilEnergyThreshold (G4double val)

G4double	GetRecoilEnergyThreshold () const

G4bool	operator== (const G4HadronicInteraction &right) const

G4bool	operator!= (const G4HadronicInteraction &right) const

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

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

void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)

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

Public Attributes
G4int	verboseLevel

G4int	MAXPART

G4int	vecLength

Public Attributes inherited from G4HEInelastic
G4int	verboseLevel

G4int	MAXPART

G4bool	conserveEnergy

G4HEVector	PionPlus

G4HEVector	PionZero

G4HEVector	PionMinus

G4HEVector	KaonPlus

G4HEVector	KaonZero

G4HEVector	AntiKaonZero

G4HEVector	KaonMinus

G4HEVector	KaonZeroShort

G4HEVector	KaonZeroLong

G4HEVector	Proton

G4HEVector	AntiProton

G4HEVector	Neutron

G4HEVector	AntiNeutron

G4HEVector	Lambda

G4HEVector	AntiLambda

G4HEVector	SigmaPlus

G4HEVector	SigmaZero

G4HEVector	SigmaMinus

G4HEVector	AntiSigmaPlus

G4HEVector	AntiSigmaZero

G4HEVector	AntiSigmaMinus

G4HEVector	XiZero

G4HEVector	XiMinus

G4HEVector	AntiXiZero

G4HEVector	AntiXiMinus

G4HEVector	OmegaMinus

G4HEVector	AntiOmegaMinus

G4HEVector	Deuteron

G4HEVector	Triton

G4HEVector	Alpha

G4HEVector	Gamma

Additional Inherited Members
Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)

G4bool	IsBlocked () const

void	Block ()

Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange

G4int	verboseLevel

G4double	theMinEnergy

G4double	theMaxEnergy

G4bool	isBlocked

Detailed Description

Definition at line 52 of file G4HESigmaZeroInelastic.hh.

Constructor & Destructor Documentation

◆ G4HESigmaZeroInelastic()

G4HESigmaZeroInelastic::G4HESigmaZeroInelastic ( )

inline

Definition at line 55 of file G4HESigmaZeroInelastic.hh.

                             : G4HEInelastic("G4HESigmaZeroInelastic")
    {
      theMinEnergy = 20*CLHEP::GeV;
      theMaxEnergy = 10*CLHEP::TeV;
      MAXPART      = 2048;
      verboseLevel = 0; 
      G4cout << "WARNING: model G4HESigmaZeroInelastic is being deprecated and will\n"
             << "disappear in Geant4 version 10.0"  << G4endl; 
    }

◆ ~G4HESigmaZeroInelastic()

G4HESigmaZeroInelastic::~G4HESigmaZeroInelastic ( )

inline

Definition at line 65 of file G4HESigmaZeroInelastic.hh.

65{};

Member Function Documentation

◆ ApplyYourself()

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

virtual

Implements G4HadronicInteraction.

Definition at line 56 of file G4HESigmaZeroInelastic.cc.

{
  G4HEVector* pv = new G4HEVector[MAXPART];
  const G4HadProjectile* aParticle = &aTrack;
  G4HEVector incidentParticle(aParticle);
     
  G4HELambdaInelastic theLambdaInelastic;
  theLambdaInelastic.SetMaxNumberOfSecondaries(MAXPART);
  theLambdaInelastic.SetVerboseLevel(verboseLevel);
    
  G4double incidentTotalMomentum = incidentParticle.getTotalMomentum();
  G4double pgam = G4UniformRand()*incidentTotalMomentum*0.75;
  G4HEVector incidentLambda; 
  incidentLambda.SmulAndUpdate(incidentParticle, 
                      (incidentTotalMomentum - pgam)/incidentTotalMomentum);                    
  G4DynamicParticle* aLambda = new G4DynamicParticle();
  aLambda->SetDefinition(G4Lambda::Lambda());
  aLambda->SetMomentum(incidentLambda.getMomentum());
  G4HadProjectile aLambdaTrack(*aLambda);
  G4HadFinalState* result =
          theLambdaInelastic.ApplyYourself(aLambdaTrack, targetNucleus);            
  vecLength = theLambdaInelastic.GetNumberOfSecondaries();
    
  pv[vecLength] = Gamma;
  pv[vecLength].setMomentum(incidentParticle.getMomentum());
  pv[vecLength].SmulAndUpdate( pv[vecLength],pgam/incidentTotalMomentum);
  G4DynamicParticle* aPhoton = new G4DynamicParticle();
  aPhoton->SetDefinition(G4Gamma::Gamma());
  aPhoton->SetMomentum(pv[vecLength].getMomentum());
  result->AddSecondary(aPhoton);
  delete [] pv;
  return result;
} 

◆ GetNumberOfSecondaries()

G4int G4HESigmaZeroInelastic::GetNumberOfSecondaries ( )

inline

Definition at line 80 of file G4HESigmaZeroInelastic.hh.

81 { return vecLength; }

◆ ModelDescription()

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

virtual

Reimplemented from G4HadronicInteraction.

Definition at line 41 of file G4HESigmaZeroInelastic.cc.

{
  outFile << "G4HESigmaZeroInelastic is one of the High Energy\n"
          << "Parameterized (HEP) models used to implement inelastic\n"
          << "Sigma0 scattering from nuclei.  It is a re-engineered\n"
          << "version of the GHEISHA code of H. Fesefeldt.  It divides the\n"
          << "initial collision products into backward- and forward-going\n"
          << "clusters which are then decayed into final state hadrons.\n"
          << "The model does not conserve energy on an event-by-event\n"
          << "basis.  It may be applied to Sigma0 with initial energies\n"
          << "above 20 GeV.\n";
}

◆ SetMaxNumberOfSecondaries()

void G4HESigmaZeroInelastic::SetMaxNumberOfSecondaries ( G4int maxnumber )

inline

Definition at line 73 of file G4HESigmaZeroInelastic.hh.

74 { MAXPART = maxnumber; };

◆ SetVerboseLevel()

void G4HESigmaZeroInelastic::SetVerboseLevel ( G4int verbose )

inline

Definition at line 75 of file G4HESigmaZeroInelastic.hh.

76 { verboseLevel = verbose;};

Member Data Documentation

◆ MAXPART

G4int G4HESigmaZeroInelastic::MAXPART

Definition at line 70 of file G4HESigmaZeroInelastic.hh.

Referenced by ApplyYourself(), G4HESigmaZeroInelastic(), and SetMaxNumberOfSecondaries().

◆ vecLength

G4int G4HESigmaZeroInelastic::vecLength

Definition at line 71 of file G4HESigmaZeroInelastic.hh.

Referenced by ApplyYourself(), and GetNumberOfSecondaries().

◆ verboseLevel

G4int G4HESigmaZeroInelastic::verboseLevel

Definition at line 69 of file G4HESigmaZeroInelastic.hh.

Referenced by ApplyYourself(), G4HESigmaZeroInelastic(), and SetVerboseLevel().

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

Public Member Functions

Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4HESigmaZeroInelastic()

◆ ~G4HESigmaZeroInelastic()

Member Function Documentation

◆ ApplyYourself()

◆ GetNumberOfSecondaries()

◆ ModelDescription()

◆ SetMaxNumberOfSecondaries()

◆ SetVerboseLevel()

Member Data Documentation

◆ MAXPART

◆ vecLength

◆ verboseLevel