#include <G4ModifiedMephi.hh>

Inheritance diagram for G4ModifiedMephi:

Public Member Functions
	G4ModifiedMephi (const G4ModifiedMephi &)=delete

	G4ModifiedMephi (const G4String &name="")

const G4String &	GetName () const

G4ModifiedMephi &	operator= (const G4ModifiedMephi &right)=delete

void	PrintGeneratorInformation () const override

G4ThreeVector &	SampleDirection (const G4DynamicParticle dp, G4double gEnergy, G4int Z, const G4Material mat=nullptr) final

virtual G4ThreeVector &	SampleDirectionForShell (const G4DynamicParticle dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material )

void	SamplePairDirections (const G4DynamicParticle dp, G4double elecKinEnergy, G4double posiKinEnergy, G4ThreeVector &dirElectron, G4ThreeVector &dirPositron, G4int Z=0, const G4Material mat=nullptr) final

	~G4ModifiedMephi () override

Protected Attributes
G4ThreeVector	fLocalDirection

G4bool	fPolarisation

Private Member Functions
G4double	SampleCosTheta (G4double primKinEnergy, G4double gEnergy, G4double mass)

Private Attributes
G4String	fName

Detailed Description

Definition at line 55 of file G4ModifiedMephi.hh.

Constructor & Destructor Documentation

◆ G4ModifiedMephi() [1/2]

G4ModifiedMephi::G4ModifiedMephi ( const G4String & name = "" )

explicit

Definition at line 51 of file G4ModifiedMephi.cc.

52 : G4VEmAngularDistribution("ModifiedMephi")

53{}

G4VEmAngularDistribution::G4VEmAngularDistribution

G4VEmAngularDistribution(const G4String &name)

Definition: G4VEmAngularDistribution.cc:55

◆ ~G4ModifiedMephi()

G4ModifiedMephi::~G4ModifiedMephi ( )

override

Definition at line 57 of file G4ModifiedMephi.cc.

58{}

◆ G4ModifiedMephi() [2/2]

G4ModifiedMephi::G4ModifiedMephi ( const G4ModifiedMephi & )

delete

Member Function Documentation

◆ GetName()

const G4String & G4VEmAngularDistribution::GetName ( ) const

inlineinherited

Definition at line 111 of file G4VEmAngularDistribution.hh.

{
  return fName;
}

References G4VEmAngularDistribution::fName.

◆ operator=()

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

delete

◆ PrintGeneratorInformation()

void G4ModifiedMephi::PrintGeneratorInformation ( ) const

overridevirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 125 of file G4ModifiedMephi.cc.

{
  G4cout << "\n" << G4endl;
  G4cout << "Angular Generator is Modified Mephi" << G4endl;
} 

References G4cout, and G4endl.

◆ SampleCosTheta()

G4double G4ModifiedMephi::SampleCosTheta	(	G4double	primKinEnergy,
		G4double	gEnergy,
		G4double	mass
	)

private

Definition at line 81 of file G4ModifiedMephi.cc.

{
  G4double gam  = 1.0 + primKinEnergy/mass;
  G4double rmax = gam*CLHEP::halfpi*std::min(1.0, gam*mass/gEnergy - 1.0);
  G4double rmax2= rmax*rmax;
  G4double x = G4UniformRand()*rmax2/(1.0 + rmax2);
 
  return std::cos(std::sqrt(x/(1.0 - x))/gam);
}

References G4UniformRand, G4InuclParticleNames::gam, CLHEP::halfpi, and G4INCL::Math::min().

Referenced by SampleDirection(), and SamplePairDirections().

◆ SampleDirection()

G4ThreeVector & G4ModifiedMephi::SampleDirection	(	const G4DynamicParticle *	dp,
		G4double	gEnergy,
		G4int	Z,
		const G4Material *	mat = `nullptr`
	)

finalvirtual

Implements G4VEmAngularDistribution.

Definition at line 63 of file G4ModifiedMephi.cc.

{
  // Sample gamma angle (Z - axis along the parent particle).
  G4double cost = SampleCosTheta(dp->GetKineticEnergy(), gEnergy, 
                                 dp->GetDefinition()->GetPDGMass());
  G4double sint = std::sqrt((1.0 - cost)*(1.0 + cost));
  G4double phi  = CLHEP::twopi*G4UniformRand(); 
 
  fLocalDirection.set(sint*std::cos(phi), sint*std::sin(phi), cost);
  fLocalDirection.rotateUz(dp->GetMomentumDirection());
 
  return fLocalDirection;
}

References G4VEmAngularDistribution::fLocalDirection, G4UniformRand, G4DynamicParticle::GetDefinition(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), G4ParticleDefinition::GetPDGMass(), CLHEP::Hep3Vector::rotateUz(), SampleCosTheta(), CLHEP::Hep3Vector::set(), and CLHEP::twopi.

◆ SampleDirectionForShell()

G4ThreeVector & G4VEmAngularDistribution::SampleDirectionForShell	(	const G4DynamicParticle *	dp,
		G4double	finalTotalEnergy,
		G4int	Z,
		G4int	shellID,
		const G4Material *	mat
	)

virtualinherited

Reimplemented in G4DeltaAngle, G4DNABornAngle, and G4DNARuddAngle.

Definition at line 69 of file G4VEmAngularDistribution.cc.

{
  return SampleDirection(dp, finalTotalEnergy, Z, mat); 
}

References G4VEmAngularDistribution::SampleDirection(), and Z.

Referenced by G4DNABornIonisationModel1::SampleSecondaries(), G4DNABornIonisationModel2::SampleSecondaries(), G4DNAEmfietzoglouIonisationModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), G4DNARuddIonisationModel::SampleSecondaries(), G4MicroElecInelasticModel::SampleSecondaries(), and G4MicroElecInelasticModel_new::SampleSecondaries().

◆ SamplePairDirections()

void G4ModifiedMephi::SamplePairDirections	(	const G4DynamicParticle *	dp,
		G4double	elecKinEnergy,
		G4double	posiKinEnergy,
		G4ThreeVector &	dirElectron,
		G4ThreeVector &	dirPositron,
		G4int	Z = `0`,
		const G4Material *	mat = `nullptr`
	)

finalvirtual

Reimplemented from G4VEmAngularDistribution.

Definition at line 95 of file G4ModifiedMephi.cc.

{
  G4double phi  = CLHEP::twopi * G4UniformRand();
  G4double sinp = std::sin(phi);
  G4double cosp = std::cos(phi);
 
  G4double ekin = dp->GetKineticEnergy();
  G4double etwo = elecKinEnergy + posiKinEnergy;
  G4double mass = dp->GetDefinition()->GetPDGMass();
 
  G4double cost = SampleCosTheta(ekin, etwo, mass);
  G4double sint = std::sqrt((1.0 - cost)*(1.0 + cost));
 
  dirElectron.set(sint*cosp, sint*sinp, cost);
  dirElectron.rotateUz(dp->GetMomentumDirection());
 
  cost = SampleCosTheta(ekin, etwo, mass);
  sint = std::sqrt((1.0 - cost)*(1.0 + cost));
 
  dirPositron.set(-sint*cosp, -sint*sinp, cost);
  dirPositron.rotateUz(dp->GetMomentumDirection());
}