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.

{
  // 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 66 of file G4SPSAngDistribution.cc.

{
    G4MUTEXDESTROY(mutex);
}

Member Function Documentation

DefineAngRefAxes()

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

Definition at line 94 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 760 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 227 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSource::ListSource().

GetDistType()

G4String G4SPSAngDistribution::GetDistType

(

)

Definition at line 197 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSource::ListSource().

GetMaxPhi()

G4double G4SPSAngDistribution::GetMaxPhi

(

)

Definition at line 221 of file G4SPSAngDistribution.cc.

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

GetMaxTheta()

G4double G4SPSAngDistribution::GetMaxTheta

(

)

Definition at line 209 of file G4SPSAngDistribution.cc.

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

GetMinPhi()

G4double G4SPSAngDistribution::GetMinPhi

(

)

Definition at line 215 of file G4SPSAngDistribution.cc.

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

GetMinTheta()

G4double G4SPSAngDistribution::GetMinTheta

(

)

Definition at line 203 of file G4SPSAngDistribution.cc.

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

ReSetHist()

void G4SPSAngDistribution::ReSetHist

(

const G4String &

atype

)

Definition at line 741 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 71 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 142 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBeamSigmaInAngX()

void G4SPSAngDistribution::SetBeamSigmaInAngX

(

G4double

r

)

Definition at line 148 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBeamSigmaInAngY()

void G4SPSAngDistribution::SetBeamSigmaInAngY

(

G4double

r

)

Definition at line 154 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetBiasRndm()

void G4SPSAngDistribution::SetBiasRndm

(

G4SPSRandomGenerator *

a

)

Definition at line 173 of file G4SPSAngDistribution.cc.

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

Referenced by G4SingleParticleSource::G4SingleParticleSource().

SetFocusPoint()

void G4SPSAngDistribution::SetFocusPoint

(

const G4ThreeVector &

input

)

Definition at line 245 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMaxPhi()

void G4SPSAngDistribution::SetMaxPhi

(

G4double

maxp

)

Definition at line 136 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMaxTheta()

void G4SPSAngDistribution::SetMaxTheta

(

G4double

maxt

)

Definition at line 130 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 124 of file G4SPSAngDistribution.cc.

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

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

SetMinTheta()

void G4SPSAngDistribution::SetMinTheta

(

G4double

mint

)

Definition at line 118 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 160 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 167 of file G4SPSAngDistribution.cc.

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

Referenced by G4SingleParticleSource::G4SingleParticleSource().

SetUserWRTSurface()

void G4SPSAngDistribution::SetUserWRTSurface

(

G4bool

wrtSurf

)

Definition at line 251 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 263 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 179 of file G4SPSAngDistribution.cc.

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

Referenced by G4SingleParticleSource::SetVerbosity().

UserDefAngPhi()

void G4SPSAngDistribution::UserDefAngPhi

(

const G4ThreeVector &

input

)

Definition at line 233 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 185 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