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 ()
G4ThreeVector	GetDirection ()
G4ParticleMomentum	GenerateOne ()

Detailed Description

Andrea Dotti Feb 2015 Important: This is a shared class between threads. Only one thread should use the set-methods here. Note that this is exactly what is achieved using UI commands. If you use the set methods to set defaults in your application take care that only one thread is executing them. In addition take care of calling these methods before the run is started Do not use these setters during the event loop

Definition at line 174 of file G4SPSAngDistribution.hh.

Constructor & Destructor Documentation

G4SPSAngDistribution::G4SPSAngDistribution

(

)

Definition at line 55 of file G4SPSAngDistribution.cc.

  : Theta(0.), Phi(0.), posDist(0),angRndm(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

(

)

Definition at line 84 of file G4SPSAngDistribution.cc.

{
    G4MUTEXDESTROY(mutex);
}

Member Function Documentation

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

Definition at line 107 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;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4ParticleMomentum G4SPSAngDistribution::GenerateOne

(

)

Definition at line 656 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;
}

Here is the caller graph for this function:

G4ThreeVector G4SPSAngDistribution::GetDirection

(

)

Definition at line 215 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

G4String G4SPSAngDistribution::GetDistType

(

)

Definition at line 210 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

G4double G4SPSAngDistribution::GetMaxPhi

(

)

Definition at line 214 of file G4SPSAngDistribution.cc.

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

G4double G4SPSAngDistribution::GetMaxTheta

(

)

Definition at line 212 of file G4SPSAngDistribution.cc.

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

G4double G4SPSAngDistribution::GetMinPhi

(

)

Definition at line 213 of file G4SPSAngDistribution.cc.

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

G4double G4SPSAngDistribution::GetMinTheta

(

)

Definition at line 211 of file G4SPSAngDistribution.cc.

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

void G4SPSAngDistribution::ReSetHist

(

G4String

atype

)

Definition at line 641 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;
  }
}

Here is the caller graph for this function:

void G4SPSAngDistribution::SetAngDistType

(

G4String

atype

)

Definition at line 90 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 ;
  }
}

Here is the caller graph for this function:

void G4SPSAngDistribution::SetBeamSigmaInAngR

(

G4double

r

)

Definition at line 153 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetBeamSigmaInAngX

(

G4double

r

)

Definition at line 159 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetBeamSigmaInAngY

(

G4double

r

)

Definition at line 165 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetBiasRndm

(

G4SPSRandomGenerator *

a

)

Definition at line 184 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetFocusPoint

(

G4ThreeVector

input

)

Definition at line 230 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetMaxPhi

(

G4double

maxp

)

Definition at line 147 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetMaxTheta

(

G4double

maxt

)

Definition at line 141 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetMinPhi

(

G4double

minp

)

Definition at line 135 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetMinTheta

(

G4double

mint

)

Definition at line 129 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetParticleMomentumDirection

(

G4ParticleMomentum

aMomentumDirection

)

Definition at line 171 of file G4SPSAngDistribution.cc.

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

}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4SPSAngDistribution::SetPosDistribution

(

G4SPSPosDistribution *

a

)

Definition at line 178 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetUserWRTSurface

(

G4bool

wrtSurf

)

Definition at line 236 of file G4SPSAngDistribution.cc.

{
    G4AutoLock l(&mutex);
  // 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;
}

Here is the caller graph for this function:

void G4SPSAngDistribution::SetUseUserAngAxis

(

G4bool

userang

)

Definition at line 247 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::SetVerbosity

(

G4int

a

)

Definition at line 190 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

void G4SPSAngDistribution::UserDefAngPhi

(

G4ThreeVector

input

)

Definition at line 217 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); 
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4SPSAngDistribution::UserDefAngTheta

(

G4ThreeVector

input

)

Definition at line 197 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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p02/source/event/include/G4SPSAngDistribution.hh
• source/geant4.10.03.p02/source/event/src/G4SPSAngDistribution.cc