G4SPSAngDistribution Class Reference

#include <G4SPSAngDistribution.hh>

Public Member Functions
	G4SPSAngDistribution ()
	~G4SPSAngDistribution ()
void	SetAngDistType (const G4String &)
void	DefineAngRefAxes (const G4String &, const G4ThreeVector &)
void	SetMinTheta (G4double)
void	SetMinPhi (G4double)
void	SetMaxTheta (G4double)
void	SetMaxPhi (G4double)
void	SetBeamSigmaInAngR (G4double)
void	SetBeamSigmaInAngX (G4double)
void	SetBeamSigmaInAngY (G4double)
void	UserDefAngTheta (const G4ThreeVector &)
void	UserDefAngPhi (const G4ThreeVector &)
void	SetFocusPoint (const G4ThreeVector &)
void	SetParticleMomentumDirection (const G4ParticleMomentum &aMomDirection)
void	SetUseUserAngAxis (G4bool)
void	SetUserWRTSurface (G4bool)
void	SetPosDistribution (G4SPSPosDistribution *a)
void	SetBiasRndm (G4SPSRandomGenerator *a)
G4ParticleMomentum	GenerateOne ()
void	ReSetHist (const G4String &)
void	SetVerbosity (G4int a)
G4String	GetDistType ()
G4double	GetMinTheta ()
G4double	GetMaxTheta ()
G4double	GetMinPhi ()
G4double	GetMaxPhi ()
G4ThreeVector	GetDirection ()

Detailed Description

Definition at line 66 of file G4SPSAngDistribution.hh.

Constructor & Destructor Documentation

G4SPSAngDistribution()

G4SPSAngDistribution::G4SPSAngDistribution

(

)

Definition at line 38 of file G4SPSAngDistribution.cc.

  : Theta(0.), Phi(0.)
{
  // Angular distribution Variables
  G4ThreeVector zero;
  particle_momentum_direction = G4ParticleMomentum(0,0,-1);
 
  AngDistType = "planar"; 
  AngRef1 = CLHEP::HepXHat;
  AngRef2 = CLHEP::HepYHat;
  AngRef3 = CLHEP::HepZHat;
  MinTheta = 0.;
  MaxTheta = pi;
  MinPhi = 0.;
  MaxPhi = twopi;
  DR = 0.;
  DX = 0.;
  DY = 0.;
  FocusPoint = G4ThreeVector(0., 0., 0.);
  UserDistType = "NULL";
  UserWRTSurface = true;
  UserAngRef = false;
  IPDFThetaExist = false;
  IPDFPhiExist = false;
  verbosityLevel = 0;
 
  G4MUTEXINIT(mutex);
}

~G4SPSAngDistribution()

G4SPSAngDistribution::~G4SPSAngDistribution

(

)

Definition at line 67 of file G4SPSAngDistribution.cc.

{
    G4MUTEXDESTROY(mutex);
}

Member Function Documentation

DefineAngRefAxes()

void G4SPSAngDistribution::DefineAngRefAxes	(	const G4String &	refname,
		const G4ThreeVector &	ref
	)

Definition at line 95 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  if (refname == "angref1")
    AngRef1 = ref.unit(); // x'
  else if (refname == "angref2")
    AngRef2 = ref.unit(); // vector in x'y' plane
 
  // User defines x' (AngRef1) and a vector in the x'y'
  // plane (AngRef2). Then, AngRef1 x AngRef2 = AngRef3
  // the z' vector. Then, AngRef3 x AngRef1 = AngRef2
  // which will now be y'.
 
  AngRef3 = AngRef1.cross(AngRef2); // z'
  AngRef2 = AngRef3.cross(AngRef1); // y'
  UserAngRef = true ;
  if(verbosityLevel == 2)
  {
    G4cout << "Angular distribution rotation axes " << AngRef1
           << " " << AngRef2 << " " << AngRef3 << G4endl;
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

GenerateOne()

G4ParticleMomentum G4SPSAngDistribution::GenerateOne

(

)

Definition at line 766 of file G4SPSAngDistribution.cc.

{
  // Local copy for thread safety
  //
  G4ParticleMomentum localM = particle_momentum_direction;
 
  // Angular stuff
  //
  if(AngDistType == "iso")
    GenerateIsotropicFlux(localM);
  else if(AngDistType == "cos")
    GenerateCosineLawFlux(localM);
  else if(AngDistType == "planar")
    GeneratePlanarFlux(localM);
  else if(AngDistType == "beam1d" || AngDistType == "beam2d" )
    GenerateBeamFlux(localM);
  else if(AngDistType == "user")
    GenerateUserDefFlux(localM);
  else if(AngDistType == "focused")
    GenerateFocusedFlux(localM);
  else
    G4cout << "Error: AngDistType has unusual value" << G4endl;
  return localM;
}

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

GetDirection()

G4ThreeVector G4SPSAngDistribution::GetDirection

(

)

Definition at line 228 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  return particle_momentum_direction;
}

Referenced by G4GeneralParticleSource::ListSource().

GetDistType()

G4String G4SPSAngDistribution::GetDistType

(

)

Definition at line 198 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  return AngDistType;
}

Referenced by G4GeneralParticleSource::ListSource().

GetMaxPhi()

G4double G4SPSAngDistribution::GetMaxPhi

(

)

Definition at line 222 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  return MaxPhi;
}

GetMaxTheta()

G4double G4SPSAngDistribution::GetMaxTheta

(

)

Definition at line 210 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  return MaxTheta;
}

GetMinPhi()

G4double G4SPSAngDistribution::GetMinPhi

(

)

Definition at line 216 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  return MinPhi;
}

GetMinTheta()

G4double G4SPSAngDistribution::GetMinTheta

(

)

Definition at line 204 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  return MinTheta;
}

ReSetHist()

void G4SPSAngDistribution::ReSetHist

(

const G4String &

atype

)

Definition at line 747 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  if (atype == "theta")
  {
    UDefThetaH = IPDFThetaH = ZeroPhysVector ;
    IPDFThetaExist = false ;
  }
  else if (atype == "phi")
  {    
    UDefPhiH = IPDFPhiH = ZeroPhysVector ;
    IPDFPhiExist = false ;
  } 
  else
  {
    G4cout << "Error, histtype not accepted " << G4endl;
  }
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetAngDistType()

void G4SPSAngDistribution::SetAngDistType

(

const G4String &

atype

)

Definition at line 72 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  if(atype != "iso" && atype != "cos" && atype != "user" && atype != "planar"
     && atype != "beam1d" && atype != "beam2d"  && atype != "focused")
  {
    G4cout << "Error, distribution must be iso, cos, planar, beam1d, beam2d, focused or user"
           << G4endl;
  }
  else
  {
    AngDistType = atype;
  }
  if (AngDistType == "cos")  { MaxTheta = pi/2.; }
  if (AngDistType == "user")
  {
    UDefThetaH = IPDFThetaH = ZeroPhysVector;
    IPDFThetaExist = false;
    UDefPhiH = IPDFPhiH = ZeroPhysVector;
    IPDFPhiExist = false;
  }
}

Referenced by G4AdjointPrimaryGenerator::G4AdjointPrimaryGenerator(), G4AdjointPrimaryGenerator::SetAdjointPrimarySourceOnAnExtSurfaceOfAVolume(), G4GeneralParticleSourceMessenger::SetNewValue(), and G4AdjointPrimaryGenerator::SetSphericalAdjointPrimarySource().

SetBeamSigmaInAngR()

void G4SPSAngDistribution::SetBeamSigmaInAngR

(

G4double

r

)

Definition at line 143 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  DR = r;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBeamSigmaInAngX()

void G4SPSAngDistribution::SetBeamSigmaInAngX

(

G4double

r

)

Definition at line 149 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  DX = r;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBeamSigmaInAngY()

void G4SPSAngDistribution::SetBeamSigmaInAngY

(

G4double

r

)

Definition at line 155 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  DY = r;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBiasRndm()

void G4SPSAngDistribution::SetBiasRndm

(

G4SPSRandomGenerator *

a

)

Definition at line 174 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  angRndm = a;
}

Referenced by G4SingleParticleSource::G4SingleParticleSource().

SetFocusPoint()

void G4SPSAngDistribution::SetFocusPoint

(

const G4ThreeVector &

input

)

Definition at line 246 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  FocusPoint = input;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMaxPhi()

void G4SPSAngDistribution::SetMaxPhi

(

G4double

maxp

)

Definition at line 137 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  MaxPhi = maxp;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMaxTheta()

void G4SPSAngDistribution::SetMaxTheta

(

G4double

maxt

)

Definition at line 131 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  MaxTheta = maxt;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue(), and G4AdjointPrimaryGenerator::SetSphericalAdjointPrimarySource().

SetMinPhi()

void G4SPSAngDistribution::SetMinPhi

(

G4double

minp

)

Definition at line 125 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  MinPhi = minp;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMinTheta()

void G4SPSAngDistribution::SetMinTheta

(

G4double

mint

)

Definition at line 119 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  MinTheta = mint;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue(), and G4AdjointPrimaryGenerator::SetSphericalAdjointPrimarySource().

SetParticleMomentumDirection()

void G4SPSAngDistribution::SetParticleMomentumDirection

(

const G4ParticleMomentum &

aMomDirection

)

Definition at line 161 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  particle_momentum_direction = aMomentumDirection.unit();
}

Referenced by G4AdjointPrimaryGenerator::GenerateAdjointPrimaryVertex(), G4AdjointPrimaryGenerator::GenerateFwdPrimaryVertex(), and G4GeneralParticleSourceMessenger::SetNewValue().

SetPosDistribution()

void G4SPSAngDistribution::SetPosDistribution

(

G4SPSPosDistribution *

a

)

Definition at line 168 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  posDist = a;
}

Referenced by G4SingleParticleSource::G4SingleParticleSource().

SetUserWRTSurface()

void G4SPSAngDistribution::SetUserWRTSurface

(

G4bool

wrtSurf

)

Definition at line 252 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
 
  // if UserWRTSurface = true then the user wants momenta with respect
  // to the surface normals.
  // When doing this theta has to be 0-90 only otherwise there will be
  // errors, which currently are flagged anywhere.
  //
  UserWRTSurface = wrtSurf;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetUseUserAngAxis()

void G4SPSAngDistribution::SetUseUserAngAxis

(

G4bool

userang

)

Definition at line 264 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
 
  // if UserAngRef = true  the angular distribution is defined wrt 
  // the user defined coordinates
  //
  UserAngRef = userang;
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetVerbosity()

void G4SPSAngDistribution::SetVerbosity

(

G4int

a

)

Definition at line 180 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  verbosityLevel = a;
}

Referenced by G4SingleParticleSource::SetVerbosity().

UserDefAngPhi()

void G4SPSAngDistribution::UserDefAngPhi

(

const G4ThreeVector &

input

)

Definition at line 234 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  if(UserDistType == "NULL") UserDistType = "phi";
  if(UserDistType == "theta") UserDistType = "both";  
  G4double phhi, val;
  phhi = input.x();
  val = input.y();
  if(verbosityLevel >= 1) G4cout << "In UserDefAngPhi" << G4endl;
  UDefPhiH.InsertValues(phhi, val); 
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

UserDefAngTheta()

void G4SPSAngDistribution::UserDefAngTheta

(

const G4ThreeVector &

input

)

Definition at line 186 of file G4SPSAngDistribution.cc.

{
  G4AutoLock l(&mutex);
  if(UserDistType == "NULL") UserDistType = "theta";
  if(UserDistType == "phi") UserDistType = "both";  
  G4double thi, val;
  thi = input.x();
  val = input.y();
  if(verbosityLevel >= 1) G4cout << "In UserDefAngTheta" << G4endl;
  UDefThetaH.InsertValues(thi, val);
}

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

---

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

• G4SPSAngDistribution.hh
• G4SPSAngDistribution.cc