G4ModifiedMephi Class Reference

#include <G4ModifiedMephi.hh>

Inheritance diagram for G4ModifiedMephi:

Public Member Functions
	G4ModifiedMephi (const G4String &name="")
	~G4ModifiedMephi () override
G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double gEnergy, G4int Z, const G4Material *mat=nullptr) final
void	SamplePairDirections (const G4DynamicParticle *dp, G4double elecKinEnergy, G4double posiKinEnergy, G4ThreeVector &dirElectron, G4ThreeVector &dirPositron, G4int Z=0, const G4Material *mat=nullptr) final
void	PrintGeneratorInformation () const override
G4ModifiedMephi &	operator= (const G4ModifiedMephi &right)=delete
	G4ModifiedMephi (const G4ModifiedMephi &)=delete
Public Member Functions inherited from G4VEmAngularDistribution
	G4VEmAngularDistribution (const G4String &name)
virtual	~G4VEmAngularDistribution ()
virtual G4ThreeVector &	SampleDirection (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
virtual G4ThreeVector &	SampleDirectionForShell (const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, G4int shellID, const G4Material *)
virtual void	SamplePairDirections (const G4DynamicParticle *dp, G4double elecKinEnergy, G4double posiKinEnergy, G4ThreeVector &dirElectron, G4ThreeVector &dirPositron, G4int Z=0, const G4Material *mat=nullptr)
virtual void	PrintGeneratorInformation () const
const G4String &	GetName () const
G4VEmAngularDistribution &	operator= (const G4VEmAngularDistribution &right)=delete
	G4VEmAngularDistribution (const G4VEmAngularDistribution &)=delete

Additional Inherited Members
Protected Attributes inherited from G4VEmAngularDistribution
G4ThreeVector	fLocalDirection
G4bool	fPolarisation

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.

  : G4VEmAngularDistribution("ModifiedMephi")
{}    

~G4ModifiedMephi()

G4ModifiedMephi::~G4ModifiedMephi ( )

override

Definition at line 57 of file G4ModifiedMephi.cc.

{}

G4ModifiedMephi() [2/2]

G4ModifiedMephi::G4ModifiedMephi ( const G4ModifiedMephi &  )

delete

Member Function Documentation

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

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

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

---

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

• G4ModifiedMephi.hh
• G4ModifiedMephi.cc