G4SPSAngDistribution Class Reference

#include <G4SPSAngDistribution.hh>

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

Detailed Description

Definition at line 155 of file G4SPSAngDistribution.hh.

Constructor & Destructor Documentation

G4SPSAngDistribution::G4SPSAngDistribution

(

)

Definition at line 54 of file G4SPSAngDistribution.cc.

References CLHEP::HepXHat, CLHEP::HepYHat, CLHEP::HepZHat, python.hepunit::pi, and python.hepunit::twopi.

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

G4SPSAngDistribution::~G4SPSAngDistribution

(

)

Definition at line 81 of file G4SPSAngDistribution.cc.

{}

Member Function Documentation

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

Definition at line 101 of file G4SPSAngDistribution.cc.

References CLHEP::Hep3Vector::cross(), G4cout, G4endl, and CLHEP::Hep3Vector::unit().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

G4ParticleMomentum G4SPSAngDistribution::GenerateOne

(

)

Definition at line 602 of file G4SPSAngDistribution.cc.

References G4cout, and G4endl.

Referenced by G4SingleParticleSource::GeneratePrimaryVertex().

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

G4String G4SPSAngDistribution::GetDistType ( )

inline

Definition at line 188 of file G4SPSAngDistribution.hh.

{ return AngDistType;}

G4double G4SPSAngDistribution::GetMaxPhi ( )

inline

Definition at line 192 of file G4SPSAngDistribution.hh.

{ return MaxPhi; }

G4double G4SPSAngDistribution::GetMaxTheta ( )

inline

Definition at line 190 of file G4SPSAngDistribution.hh.

{ return MaxTheta; }

G4double G4SPSAngDistribution::GetMinPhi ( )

inline

Definition at line 191 of file G4SPSAngDistribution.hh.

{ return MinPhi; }

G4double G4SPSAngDistribution::GetMinTheta ( )

inline

Definition at line 189 of file G4SPSAngDistribution.hh.

{ return MinTheta; }

void G4SPSAngDistribution::ReSetHist

(

G4String

atype

)

Definition at line 588 of file G4SPSAngDistribution.cc.

References G4cout, and G4endl.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

void G4SPSAngDistribution::SetAngDistType

(

G4String

atype

)

Definition at line 85 of file G4SPSAngDistribution.cc.

References G4cout, G4endl, and python.hepunit::pi.

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

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

void G4SPSAngDistribution::SetBeamSigmaInAngR

(

G4double

r

)

Definition at line 142 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  DR = r;
}

void G4SPSAngDistribution::SetBeamSigmaInAngX

(

G4double

r

)

Definition at line 147 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  DX = r;
}

void G4SPSAngDistribution::SetBeamSigmaInAngY

(

G4double

r

)

Definition at line 152 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  DY = r;
}

void G4SPSAngDistribution::SetBiasRndm ( G4SPSRandomGenerator *  a )

inline

Definition at line 181 of file G4SPSAngDistribution.hh.

References test::a.

Referenced by G4SingleParticleSource::G4SingleParticleSource().

{angRndm = a;}

void G4SPSAngDistribution::SetFocusPoint

(

G4ThreeVector

input

)

Definition at line 181 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  FocusPoint = input;
}

void G4SPSAngDistribution::SetMaxPhi

(

G4double

maxp

)

Definition at line 137 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  MaxPhi = maxp;
}

void G4SPSAngDistribution::SetMaxTheta

(

G4double

maxt

)

Definition at line 132 of file G4SPSAngDistribution.cc.

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

{
  MaxTheta = maxt;
}

void G4SPSAngDistribution::SetMinPhi

(

G4double

minp

)

Definition at line 127 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  MinPhi = minp;
}

void G4SPSAngDistribution::SetMinTheta

(

G4double

mint

)

Definition at line 122 of file G4SPSAngDistribution.cc.

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

{
  MinTheta = mint;
}

void G4SPSAngDistribution::SetParticleMomentumDirection ( G4ParticleMomentum  aMomentumDirection )

inline

Definition at line 175 of file G4SPSAngDistribution.hh.

References CLHEP::Hep3Vector::unit().

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

  { particle_momentum_direction =  aMomentumDirection.unit(); }

void G4SPSAngDistribution::SetPosDistribution ( G4SPSPosDistribution *  a )

inline

Definition at line 180 of file G4SPSAngDistribution.hh.

References test::a.

Referenced by G4SingleParticleSource::G4SingleParticleSource().

{posDist = a; }

void G4SPSAngDistribution::SetUserWRTSurface

(

G4bool

wrtSurf

)

Definition at line 186 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  // This is only applied in user mode?
  // 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;
}

void G4SPSAngDistribution::SetUseUserAngAxis

(

G4bool

userang

)

Definition at line 196 of file G4SPSAngDistribution.cc.

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

{
  // if UserAngRef = true  the angular distribution is defined wrt 
  // the user defined co-ordinates
  UserAngRef = userang;
}

void G4SPSAngDistribution::SetVerbosity ( G4int  a )

inline

Definition at line 186 of file G4SPSAngDistribution.hh.

References test::a.

Referenced by G4SingleParticleSource::SetVerbosity().

{verbosityLevel = a; }

void G4SPSAngDistribution::UserDefAngPhi

(

G4ThreeVector

input

)

Definition at line 169 of file G4SPSAngDistribution.cc.

References G4cout, G4endl, G4PhysicsOrderedFreeVector::InsertValues(), CLHEP::Hep3Vector::x(), and CLHEP::Hep3Vector::y().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

void G4SPSAngDistribution::UserDefAngTheta

(

G4ThreeVector

input

)

Definition at line 157 of file G4SPSAngDistribution.cc.

References G4cout, G4endl, G4PhysicsOrderedFreeVector::InsertValues(), CLHEP::Hep3Vector::x(), and CLHEP::Hep3Vector::y().

Referenced by G4GeneralParticleSourceMessenger::SetNewValue().

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

---

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

• G4SPSAngDistribution.hh
• G4SPSAngDistribution.cc