G4SPSAngDistribution Class Reference

#include <G4SPSAngDistribution.hh>

Collaboration diagram for G4SPSAngDistribution:

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

Private Member Functions
void	GenerateFocusedFlux (G4ParticleMomentum &outputMom)
void	GenerateIsotropicFlux (G4ParticleMomentum &outputMom)
void	GenerateCosineLawFlux (G4ParticleMomentum &outputMom)
void	GenerateBeamFlux (G4ParticleMomentum &outputMom)
void	GeneratePlanarFlux (G4ParticleMomentum &outputMom)
void	GenerateUserDefFlux (G4ParticleMomentum &outputMom)
G4double	GenerateUserDefTheta ()
G4double	GenerateUserDefPhi ()

Private Attributes
G4String	AngDistType
G4ThreeVector	AngRef1
G4ThreeVector	AngRef2
G4ThreeVector	AngRef3
G4double	MinTheta
G4double	MaxTheta
G4double	MinPhi
G4double	MaxPhi
G4double	DR
G4double	DX
G4double	DY
G4double	Theta
G4double	Phi
G4ThreeVector	FocusPoint
G4bool	IPDFThetaExist
G4bool	IPDFPhiExist
G4PhysicsOrderedFreeVector	UDefThetaH
G4PhysicsOrderedFreeVector	IPDFThetaH
G4PhysicsOrderedFreeVector	UDefPhiH
G4PhysicsOrderedFreeVector	IPDFPhiH
G4String	UserDistType
G4bool	UserWRTSurface
G4bool	UserAngRef
G4ParticleMomentum	particle_momentum_direction
G4SPSPosDistribution *	posDist
G4SPSRandomGenerator *	angRndm
G4int	verbosityLevel
G4PhysicsOrderedFreeVector	ZeroPhysVector
G4Mutex	mutex

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::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::~G4SPSAngDistribution

(

)

Definition at line 84 of file G4SPSAngDistribution.cc.

{
    G4MUTEXDESTROY(mutex);
}

Member Function Documentation

DefineAngRefAxes()

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:

GenerateBeamFlux()

void G4SPSAngDistribution::GenerateBeamFlux ( G4ParticleMomentum &  outputMom )

private

Definition at line 255 of file G4SPSAngDistribution.cc.

{
  G4double theta, phi;
  G4double px, py, pz;
  if (AngDistType == "beam1d")
    { 
      theta = G4RandGauss::shoot(0.0,DR);
      phi = twopi * G4UniformRand();
    }
  else 
    {
      px = G4RandGauss::shoot(0.0,DX);
      py = G4RandGauss::shoot(0.0,DY);
      theta = std::sqrt (px*px + py*py);
      if (theta != 0.) { 
            phi = std::acos(px/theta);
            if ( py < 0.) phi = -phi;
        } else {
            phi = 0.0;
        }
    }
  px = -std::sin(theta) * std::cos(phi);
  py = -std::sin(theta) * std::sin(phi);
  pz = -std::cos(theta);
  G4double finx, finy,  finz ;
  finx = px, finy =py, finz =pz;
  if (UserAngRef){
    // Apply Angular Rotation Matrix
    // x * AngRef1, y * AngRef2 and z * AngRef3
    finx = (px * AngRef1.x()) + (py * AngRef2.x()) + (pz * AngRef3.x());
    finy = (px * AngRef1.y()) + (py * AngRef2.y()) + (pz * AngRef3.y());
    finz = (px * AngRef1.z()) + (py * AngRef2.z()) + (pz * AngRef3.z());
    G4double ResMag = std::sqrt((finx*finx) + (finy*finy) + (finz*finz));
    finx = finx/ResMag;
    finy = finy/ResMag;
    finz = finz/ResMag;
  }
  mom.setX(finx);
  mom.setY(finy);
  mom.setZ(finz);

  // particle_momentum_direction now holds unit momentum vector.
  if(verbosityLevel >= 1)
    G4cout << "Generating beam vector: " << mom << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateCosineLawFlux()

void G4SPSAngDistribution::GenerateCosineLawFlux ( G4ParticleMomentum &  outputMom )

private

Definition at line 375 of file G4SPSAngDistribution.cc.

{
  // Method to generate flux distributed with a cosine law
  G4double px, py, pz;
  G4double rndm, rndm2;
  //
  G4double sintheta, sinphi,costheta,cosphi;
  rndm = angRndm->GenRandTheta();
  sintheta = std::sqrt( rndm * (std::sin(MaxTheta)*std::sin(MaxTheta) - std::sin(MinTheta)*std::sin(MinTheta) ) 
           +std::sin(MinTheta)*std::sin(MinTheta) );
  costheta = std::sqrt(1. -sintheta*sintheta);
  
  rndm2 = angRndm->GenRandPhi();
  Phi = MinPhi + (MaxPhi - MinPhi) * rndm2; 
  sinphi = std::sin(Phi);
  cosphi = std::cos(Phi);

  px = -sintheta * cosphi;
  py = -sintheta * sinphi;
  pz = -costheta;

  // for volume and ponit source use mother or user defined co-ordinates
  // for plane and surface source user surface-normal or userdefined co-ordinates
  //
  G4double finx, finy, finz;
  if (posDist->GetSourcePosType() == "Point" || posDist->GetSourcePosType() == "Volume") {
    if (UserAngRef){
      // Apply Rotation Matrix
      finx = (px * AngRef1.x()) + (py * AngRef2.x()) + (pz * AngRef3.x());
      finy = (px * AngRef1.y()) + (py * AngRef2.y()) + (pz * AngRef3.y());
      finz = (px * AngRef1.z()) + (py * AngRef2.z()) + (pz * AngRef3.z());
    } else {
      finx = px;
      finy = py;
      finz = pz;
    }
  } else {    // for plane and surface source   
    if (UserAngRef){
      // Apply Rotation Matrix
      finx = (px * AngRef1.x()) + (py * AngRef2.x()) + (pz * AngRef3.x());
      finy = (px * AngRef1.y()) + (py * AngRef2.y()) + (pz * AngRef3.y());
      finz = (px * AngRef1.z()) + (py * AngRef2.z()) + (pz * AngRef3.z());
    } else {
      finx = (px*posDist->GetSideRefVec1().x()) + (py*posDist->GetSideRefVec2().x()) + (pz*posDist->GetSideRefVec3().x());
      finy = (px*posDist->GetSideRefVec1().y()) + (py*posDist->GetSideRefVec2().y()) + (pz*posDist->GetSideRefVec3().y());
      finz = (px*posDist->GetSideRefVec1().z()) + (py*posDist->GetSideRefVec2().z()) + (pz*posDist->GetSideRefVec3().z());
    }
  }
  G4double ResMag = std::sqrt((finx*finx) + (finy*finy) + (finz*finz));
  finx = finx/ResMag;
  finy = finy/ResMag;
  finz = finz/ResMag;

  mom.setX(finx);
  mom.setY(finy);
  mom.setZ(finz);

  // particle_momentum_direction now contains unit momentum vector.
  if(verbosityLevel >= 1)
    {
      G4cout << "Resultant cosine-law unit momentum vector " << mom << G4endl;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateFocusedFlux()

void G4SPSAngDistribution::GenerateFocusedFlux ( G4ParticleMomentum &  outputMom )

private

Definition at line 301 of file G4SPSAngDistribution.cc.

{
  mom = (FocusPoint - posDist->GetParticlePos()).unit();
  // 
  // particle_momentum_direction now holds unit momentum vector.
  if(verbosityLevel >= 1)
    G4cout << "Generating focused vector: " << mom << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateIsotropicFlux()

void G4SPSAngDistribution::GenerateIsotropicFlux ( G4ParticleMomentum &  outputMom )

private

Definition at line 310 of file G4SPSAngDistribution.cc.

{
  // generates isotropic flux.
  // No vectors are needed.
  G4double rndm, rndm2;
  G4double px, py, pz;

  //
  G4double sintheta, sinphi,costheta,cosphi;
  rndm = angRndm->GenRandTheta();
  costheta = std::cos(MinTheta) - rndm * (std::cos(MinTheta) - std::cos(MaxTheta));
  sintheta = std::sqrt(1. - costheta*costheta);
  
  rndm2 = angRndm->GenRandPhi();
  Phi = MinPhi + (MaxPhi - MinPhi) * rndm2; 
  sinphi = std::sin(Phi);
  cosphi = std::cos(Phi);

  px = -sintheta * cosphi;
  py = -sintheta * sinphi;
  pz = -costheta;

  // for volume and ponit source use mother or user defined co-ordinates
  // for plane and surface source user surface-normal or userdefined co-ordinates
  //
  G4double finx, finy, finz;
  if (posDist->GetSourcePosType() == "Point" || posDist->GetSourcePosType() == "Volume") {
    if (UserAngRef){
      // Apply Rotation Matrix
      // x * AngRef1, y * AngRef2 and z * AngRef3
      finx = (px * AngRef1.x()) + (py * AngRef2.x()) + (pz * AngRef3.x());
      finy = (px * AngRef1.y()) + (py * AngRef2.y()) + (pz * AngRef3.y());
      finz = (px * AngRef1.z()) + (py * AngRef2.z()) + (pz * AngRef3.z());
    } else {
      finx = px;
      finy = py;
      finz = pz;
    }
  } else {    // for plane and surface source   
    if (UserAngRef){
      // Apply Rotation Matrix
      // x * AngRef1, y * AngRef2 and z * AngRef3
      finx = (px * AngRef1.x()) + (py * AngRef2.x()) + (pz * AngRef3.x());
      finy = (px * AngRef1.y()) + (py * AngRef2.y()) + (pz * AngRef3.y());
      finz = (px * AngRef1.z()) + (py * AngRef2.z()) + (pz * AngRef3.z());
    } else {
      finx = (px*posDist->GetSideRefVec1().x()) + (py*posDist->GetSideRefVec2().x()) + (pz*posDist->GetSideRefVec3().x());
      finy = (px*posDist->GetSideRefVec1().y()) + (py*posDist->GetSideRefVec2().y()) + (pz*posDist->GetSideRefVec3().y());
      finz = (px*posDist->GetSideRefVec1().z()) + (py*posDist->GetSideRefVec2().z()) + (pz*posDist->GetSideRefVec3().z());
    }
  }
  G4double ResMag = std::sqrt((finx*finx) + (finy*finy) + (finz*finz));
  finx = finx/ResMag;
  finy = finy/ResMag;
  finz = finz/ResMag;

  mom.setX(finx);
  mom.setY(finy);
  mom.setZ(finz);

  // particle_momentum_direction now holds unit momentum vector.
  if(verbosityLevel >= 1)
    G4cout << "Generating isotropic vector: " << mom << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateOne()

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 call graph for this function:

Here is the caller graph for this function:

GeneratePlanarFlux()

void G4SPSAngDistribution::GeneratePlanarFlux ( G4ParticleMomentum &  outputMom )

private

Definition at line 439 of file G4SPSAngDistribution.cc.

{
  // particle_momentum_direction now contains unit momentum vector.
  // nothing need be done here as the m-directions have been set directly
  // under this option
  if(verbosityLevel >= 1)
    {
      G4cout << "Resultant Planar wave  momentum vector " << mom << G4endl;
    }
}

Here is the caller graph for this function:

GenerateUserDefFlux()

void G4SPSAngDistribution::GenerateUserDefFlux ( G4ParticleMomentum &  outputMom )

private

Definition at line 450 of file G4SPSAngDistribution.cc.

{
  G4double rndm, px, py, pz, pmag;

  if(UserDistType == "NULL")
    G4cout << "Error: UserDistType undefined" << G4endl;
  else if(UserDistType == "theta") {
    Theta = 10.;
    while(Theta > MaxTheta || Theta < MinTheta)
      Theta = GenerateUserDefTheta();
    Phi = 10.;
    while(Phi > MaxPhi || Phi < MinPhi) {
      rndm = angRndm->GenRandPhi();
      Phi = twopi * rndm;
    }
  }
  else if(UserDistType == "phi") {
    Theta = 10.;
    while(Theta > MaxTheta || Theta < MinTheta)
      {
    rndm = angRndm->GenRandTheta();
    Theta = std::acos(1. - (2. * rndm));
      }
    Phi = 10.;
    while(Phi > MaxPhi || Phi < MinPhi)
      Phi = GenerateUserDefPhi();
  }
  else if(UserDistType == "both")
    {
      Theta = 10.;
      while(Theta > MaxTheta || Theta < MinTheta)      
    Theta = GenerateUserDefTheta();
      Phi = 10.;
      while(Phi > MaxPhi || Phi < MinPhi)
    Phi = GenerateUserDefPhi();
    }
  px = -std::sin(Theta) * std::cos(Phi);
  py = -std::sin(Theta) * std::sin(Phi);
  pz = -std::cos(Theta);

  pmag = std::sqrt((px*px) + (py*py) + (pz*pz));

  if(!UserWRTSurface) {
    G4double finx, finy, finz;      
    if (UserAngRef) {
      // Apply Rotation Matrix
      // x * AngRef1, y * AngRef2 and z * AngRef3
      finx = (px * AngRef1.x()) + (py * AngRef2.x()) + (pz * AngRef3.x());
      finy = (px * AngRef1.y()) + (py * AngRef2.y()) + (pz * AngRef3.y());
      finz = (px * AngRef1.z()) + (py * AngRef2.z()) + (pz * AngRef3.z());
    } else {  // use mother co-ordinates
      finx = px;
      finy = py;
      finz = pz;
    }
    G4double ResMag = std::sqrt((finx*finx) + (finy*finy) + (finz*finz));
    finx = finx/ResMag;
    finy = finy/ResMag;
    finz = finz/ResMag;
    
    mom.setX(finx);
    mom.setY(finy);
    mom.setZ(finz);
  } 
  else  {  // UserWRTSurface = true
    G4double pxh = px/pmag;
    G4double pyh = py/pmag;
    G4double pzh = pz/pmag;
    if(verbosityLevel > 1) {
      G4cout <<"SideRefVecs " <<posDist->GetSideRefVec1()<<posDist->GetSideRefVec2()<<posDist->GetSideRefVec3()<<G4endl;
      G4cout <<"Raw Unit vector "<<pxh<<","<<pyh<<","<<pzh<<G4endl;
    }
    G4double resultx = (pxh*posDist->GetSideRefVec1().x()) + (pyh*posDist->GetSideRefVec2().x()) +
      (pzh*posDist->GetSideRefVec3().x());
    
    G4double resulty = (pxh*posDist->GetSideRefVec1().y()) + (pyh*posDist->GetSideRefVec2().y()) +
      (pzh*posDist->GetSideRefVec3().y());
    
    G4double resultz = (pxh*posDist->GetSideRefVec1().z()) + (pyh*posDist->GetSideRefVec2().z()) +
      (pzh*posDist->GetSideRefVec3().z());
    
    G4double ResMag = std::sqrt((resultx*resultx) + (resulty*resulty) + (resultz*resultz));
    resultx = resultx/ResMag;
    resulty = resulty/ResMag;
    resultz = resultz/ResMag;
    
    mom.setX(resultx);
    mom.setY(resulty);
    mom.setZ(resultz);
  }
  
  // particle_momentum_direction now contains unit momentum vector.
  if(verbosityLevel > 0 )
    {
      G4cout << "Final User Defined momentum vector " << particle_momentum_direction << G4endl;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateUserDefPhi()

G4double G4SPSAngDistribution::GenerateUserDefPhi ( )

private

Definition at line 594 of file G4SPSAngDistribution.cc.

{
  // Create cumulative histogram if not already done so. Then use RandFlat
  //::shoot to generate the output Theta value.

  if(UserDistType == "NULL" || UserDistType == "theta")
    {
      // No user defined phi distribution
      G4cout << "Error ***********************" << G4endl;
      G4cout << "UserDistType = " << UserDistType << G4endl;
      return(0.);
    }
  else
    {
      // UserDistType = phi or both and so a phi distribution
      // is defined. This should be integrated if not already done.
      G4AutoLock l(&mutex);
      if(IPDFPhiExist == false)
      {
          // IPDF has not been created, so create it
          G4double bins[1024],vals[1024], sum;
          G4int ii;
          G4int maxbin = G4int(UDefPhiH.GetVectorLength());
          bins[0] = UDefPhiH.GetLowEdgeEnergy(size_t(0));
          vals[0] = UDefPhiH(size_t(0));
          sum = vals[0];
          for(ii=1;ii<maxbin;ii++)
          {
              bins[ii] = UDefPhiH.GetLowEdgeEnergy(size_t(ii));
              vals[ii] = UDefPhiH(size_t(ii)) + vals[ii-1];
              sum = sum + UDefPhiH(size_t(ii));
          }
          for(ii=0;ii<maxbin;ii++)
          {
              vals[ii] = vals[ii]/sum;
              IPDFPhiH.InsertValues(bins[ii], vals[ii]);
          }
          // Make IPDFPhiExist = true
          IPDFPhiExist = true;
      }
      l.unlock();
      // IPDF has been create so carry on
      G4double rndm = G4UniformRand();
      return(IPDFPhiH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateUserDefTheta()

G4double G4SPSAngDistribution::GenerateUserDefTheta ( )

private

Definition at line 548 of file G4SPSAngDistribution.cc.

{
  // Create cumulative histogram if not already done so. Then use RandFlat
  //::shoot to generate the output Theta value.
  if(UserDistType == "NULL" || UserDistType == "phi")
    {
      // No user defined theta distribution
      G4cout << "Error ***********************" << G4endl;
      G4cout << "UserDistType = " << UserDistType << G4endl;
      return (0.);
    }
  else
    {
      // UserDistType = theta or both and so a theta distribution
      // is defined. This should be integrated if not already done.
      G4AutoLock l(&mutex);
      if(IPDFThetaExist == false)
      {
          // IPDF has not been created, so create it
          G4double bins[1024],vals[1024], sum;
          G4int ii;
          G4int maxbin = G4int(UDefThetaH.GetVectorLength());
          bins[0] = UDefThetaH.GetLowEdgeEnergy(size_t(0));
          vals[0] = UDefThetaH(size_t(0));
          sum = vals[0];
          for(ii=1;ii<maxbin;ii++)
          {
              bins[ii] = UDefThetaH.GetLowEdgeEnergy(size_t(ii));
              vals[ii] = UDefThetaH(size_t(ii)) + vals[ii-1];
              sum = sum + UDefThetaH(size_t(ii));
          }
          for(ii=0;ii<maxbin;ii++)
          {
              vals[ii] = vals[ii]/sum;
              IPDFThetaH.InsertValues(bins[ii], vals[ii]);
          }
          // Make IPDFThetaExist = true
          IPDFThetaExist = true;
      }
      l.unlock();
      // IPDF has been create so carry on
      G4double rndm = G4UniformRand();
      return(IPDFThetaH.GetEnergy(rndm));
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDirection()

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:

GetDistType()

G4String G4SPSAngDistribution::GetDistType

(

)

Definition at line 210 of file G4SPSAngDistribution.cc.

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

Here is the caller graph for this function:

GetMaxPhi()

G4double G4SPSAngDistribution::GetMaxPhi

(

)

Definition at line 214 of file G4SPSAngDistribution.cc.

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

GetMaxTheta()

G4double G4SPSAngDistribution::GetMaxTheta

(

)

Definition at line 212 of file G4SPSAngDistribution.cc.

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

GetMinPhi()

G4double G4SPSAngDistribution::GetMinPhi

(

)

Definition at line 213 of file G4SPSAngDistribution.cc.

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

GetMinTheta()

G4double G4SPSAngDistribution::GetMinTheta

(

)

Definition at line 211 of file G4SPSAngDistribution.cc.

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

ReSetHist()

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:

SetAngDistType()

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:

SetBeamSigmaInAngR()

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:

SetBeamSigmaInAngX()

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:

SetBeamSigmaInAngY()

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:

SetBiasRndm()

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:

SetFocusPoint()

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:

SetMaxPhi()

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:

SetMaxTheta()

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:

SetMinPhi()

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:

SetMinTheta()

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:

SetParticleMomentumDirection()

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:

SetPosDistribution()

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:

SetUserWRTSurface()

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:

SetUseUserAngAxis()

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:

SetVerbosity()

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:

UserDefAngPhi()

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:

UserDefAngTheta()

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:

Member Data Documentation

AngDistType

G4String G4SPSAngDistribution::AngDistType

private

Definition at line 228 of file G4SPSAngDistribution.hh.

AngRef1

G4ThreeVector G4SPSAngDistribution::AngRef1

private

Definition at line 229 of file G4SPSAngDistribution.hh.

AngRef2

G4ThreeVector G4SPSAngDistribution::AngRef2

private

Definition at line 229 of file G4SPSAngDistribution.hh.

AngRef3

G4ThreeVector G4SPSAngDistribution::AngRef3

private

Definition at line 229 of file G4SPSAngDistribution.hh.

angRndm

G4SPSRandomGenerator* G4SPSAngDistribution::angRndm

private

Definition at line 247 of file G4SPSAngDistribution.hh.

DR

G4double G4SPSAngDistribution::DR

private

Definition at line 231 of file G4SPSAngDistribution.hh.

DX

G4double G4SPSAngDistribution::DX

private

Definition at line 231 of file G4SPSAngDistribution.hh.

DY

G4double G4SPSAngDistribution::DY

private

Definition at line 231 of file G4SPSAngDistribution.hh.

FocusPoint

G4ThreeVector G4SPSAngDistribution::FocusPoint

private

Definition at line 233 of file G4SPSAngDistribution.hh.

IPDFPhiExist

G4bool G4SPSAngDistribution::IPDFPhiExist

private

Definition at line 234 of file G4SPSAngDistribution.hh.

IPDFPhiH

G4PhysicsOrderedFreeVector G4SPSAngDistribution::IPDFPhiH

private

Definition at line 238 of file G4SPSAngDistribution.hh.

IPDFThetaExist

G4bool G4SPSAngDistribution::IPDFThetaExist

private

Definition at line 234 of file G4SPSAngDistribution.hh.

IPDFThetaH

G4PhysicsOrderedFreeVector G4SPSAngDistribution::IPDFThetaH

private

Definition at line 236 of file G4SPSAngDistribution.hh.

MaxPhi

G4double G4SPSAngDistribution::MaxPhi

private

Definition at line 230 of file G4SPSAngDistribution.hh.

MaxTheta

G4double G4SPSAngDistribution::MaxTheta

private

Definition at line 230 of file G4SPSAngDistribution.hh.

MinPhi

G4double G4SPSAngDistribution::MinPhi

private

Definition at line 230 of file G4SPSAngDistribution.hh.

MinTheta

G4double G4SPSAngDistribution::MinTheta

private

Definition at line 230 of file G4SPSAngDistribution.hh.

mutex

G4Mutex G4SPSAngDistribution::mutex

private

Definition at line 254 of file G4SPSAngDistribution.hh.

particle_momentum_direction

G4ParticleMomentum G4SPSAngDistribution::particle_momentum_direction

private

Definition at line 244 of file G4SPSAngDistribution.hh.

Phi

G4double G4SPSAngDistribution::Phi

private

Definition at line 232 of file G4SPSAngDistribution.hh.

posDist

G4SPSPosDistribution* G4SPSAngDistribution::posDist

private

Definition at line 246 of file G4SPSAngDistribution.hh.

Theta

G4double G4SPSAngDistribution::Theta

private

Definition at line 232 of file G4SPSAngDistribution.hh.

UDefPhiH

G4PhysicsOrderedFreeVector G4SPSAngDistribution::UDefPhiH

private

Definition at line 237 of file G4SPSAngDistribution.hh.

UDefThetaH

G4PhysicsOrderedFreeVector G4SPSAngDistribution::UDefThetaH

private

Definition at line 235 of file G4SPSAngDistribution.hh.

UserAngRef

G4bool G4SPSAngDistribution::UserAngRef

private

Definition at line 242 of file G4SPSAngDistribution.hh.

UserDistType

G4String G4SPSAngDistribution::UserDistType

private

Definition at line 239 of file G4SPSAngDistribution.hh.

UserWRTSurface

G4bool G4SPSAngDistribution::UserWRTSurface

private

Definition at line 240 of file G4SPSAngDistribution.hh.

verbosityLevel

G4int G4SPSAngDistribution::verbosityLevel

private

Definition at line 250 of file G4SPSAngDistribution.hh.

ZeroPhysVector

G4PhysicsOrderedFreeVector G4SPSAngDistribution::ZeroPhysVector

private

Definition at line 252 of file G4SPSAngDistribution.hh.

---

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

• G4SPSAngDistribution.hh
• G4SPSAngDistribution.cc