#include <G4SPSPosDistribution.hh>

Public Member Functions
	G4SPSPosDistribution ()

	~G4SPSPosDistribution ()

void	SetPosDisType (G4String)

void	SetPosDisShape (G4String)

void	SetCentreCoords (G4ThreeVector)

void	SetPosRot1 (G4ThreeVector)

void	SetPosRot2 (G4ThreeVector)

void	SetHalfX (G4double)

void	SetHalfY (G4double)

void	SetHalfZ (G4double)

void	SetRadius (G4double)

void	SetRadius0 (G4double)

void	SetBeamSigmaInR (G4double)

void	SetBeamSigmaInX (G4double)

void	SetBeamSigmaInY (G4double)

void	SetParAlpha (G4double)

void	SetParTheta (G4double)

void	SetParPhi (G4double)

void	ConfineSourceToVolume (G4String)

void	SetBiasRndm (G4SPSRandomGenerator *a)

void	SetVerbosity (G4int a)

G4ThreeVector	GenerateOne ()

G4String	GetPosDisType () const

G4String	GetPosDisShape () const

G4ThreeVector	GetCentreCoords () const

G4double	GetHalfX () const

G4double	GetHalfY () const

G4double	GetHalfZ () const

G4double	GetRadius () const

G4ThreeVector	GetSideRefVec1 () const

G4ThreeVector	GetSideRefVec2 () const

G4ThreeVector	GetSideRefVec3 () const

G4String	GetSourcePosType () const

G4ThreeVector	GetParticlePos () const

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 156 of file G4SPSPosDistribution.hh.

Constructor & Destructor Documentation

G4SPSPosDistribution::G4SPSPosDistribution ( )

Definition at line 67 of file G4SPSPosDistribution.cc.

                                            : PosRndm(0)
 {
   SourcePosType = "Point";
   Shape = "NULL";
   CentreCoords = G4ThreeVector(0,0,0);
   Rotx = CLHEP::HepXHat;
   Roty = CLHEP::HepYHat;
   Rotz = CLHEP::HepZHat;
   halfx = 0.;
   halfy = 0.;
   halfz = 0.;
   Radius = 0.;
   Radius0 = 0.;
   SR = 0.;
   SX = 0.;
   SY = 0.;
   ParAlpha = 0.;
   ParTheta = 0.;
   ParPhi = 0.;
   Confine = false; //If true confines source distribution to VolName
   VolName = "NULL";
   verbosityLevel = 0 ;
   G4MUTEXINIT(a_mutex);
 }

G4SPSPosDistribution::~G4SPSPosDistribution ( )

Definition at line 92 of file G4SPSPosDistribution.cc.

 {
   G4MUTEXDESTROY(a_mutex);
 }

Member Function Documentation

void G4SPSPosDistribution::ConfineSourceToVolume ( G4String Vname )

Definition at line 275 of file G4SPSPosDistribution.cc.

 {
   VolName = Vname;
   if(verbosityLevel == 2)
     G4cout << VolName << G4endl;
   G4VPhysicalVolume *tempPV      = NULL;
   G4PhysicalVolumeStore *PVStore = 0;
   G4String theRequiredVolumeName = VolName;
   PVStore      = G4PhysicalVolumeStore::GetInstance();
   G4int      i = 0;
   G4bool found = false;
   if(verbosityLevel == 2)
     G4cout << PVStore->size() << G4endl;
   while (!found && i<G4int(PVStore->size())) {
     tempPV = (*PVStore)[i];
     found  = tempPV->GetName() == theRequiredVolumeName;
     if(verbosityLevel == 2)
       G4cout << i << " " << " " << tempPV->GetName() << " " << theRequiredVolumeName << " " << found << G4endl;
     if (!found)
       {i++;}
   }
   // found = true then the volume exists else it doesnt.
   if(found == true)
     {
       if(verbosityLevel >= 1)
     G4cout << "Volume " << VolName << " exists" << G4endl;
       Confine = true;
     }
   else
     {
       G4cout << " **** Error: Volume does not exist **** " << G4endl;
       G4cout << " Ignoring confine condition" << G4endl;
       Confine = false;
       VolName = "NULL";
     }
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector G4SPSPosDistribution::GenerateOne ( )

Definition at line 1049 of file G4SPSPosDistribution.cc.

 {
   //
   G4ThreeVector localP;
   G4bool srcconf = false;
   G4int LoopCount = 0;
   while(srcconf == false)
     {
       if(SourcePosType == "Point")
           GeneratePointSource(localP);
       else if(SourcePosType == "Beam")
           GeneratePointsInBeam(localP);
       else if(SourcePosType == "Plane")
           GeneratePointsInPlane(localP);
       else if(SourcePosType == "Surface")
           GeneratePointsOnSurface(localP);
       else if(SourcePosType == "Volume")
           GeneratePointsInVolume(localP);
       else
       {
           G4ExceptionDescription msg;
           msg << "Error: SourcePosType undefined\n";
           msg << "Generating point source\n";
           G4Exception("G4SPSPosDistribution::GenerateOne()","G4GPS001",JustWarning,msg);
           GeneratePointSource(localP);
       }
       if(Confine == true)
       {
           srcconf = IsSourceConfined(localP);
           // if source in confined srcconf = true terminating the loop
           // if source isnt confined srcconf = false and loop continues
       }
       else if(Confine == false)
           srcconf = true; // terminate loop
       LoopCount++;
       if(LoopCount == 100000)
       {
           G4ExceptionDescription msg;
           msg << "LoopCount = 100000\n";
           msg << "Either the source distribution >> confinement\n";
           msg << "or any confining volume may not overlap with\n";
           msg << "the source distribution or any confining volumes\n";
           msg << "may not exist\n"<< G4endl;
           msg << "If you have set confine then this will be ignored\n";
           msg << "for this event.\n" << G4endl;
           G4Exception("G4SPSPosDistribution::GenerateOne()","G4GPS001",JustWarning,msg);
           srcconf = true; //Avoids an infinite loop
       }
     }
   ThreadData.Get().CParticlePos=localP;
   return localP;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector G4SPSPosDistribution::GetCentreCoords ( ) const

Definition at line 201 of file G4SPSPosDistribution.cc.

 {
     return CentreCoords;
 }

Here is the caller graph for this function:

G4double G4SPSPosDistribution::GetHalfX ( ) const

Definition at line 206 of file G4SPSPosDistribution.cc.

 {
     return halfx;
 }

G4double G4SPSPosDistribution::GetHalfY ( ) const

Definition at line 211 of file G4SPSPosDistribution.cc.

 {
     return halfy;
 }

G4double G4SPSPosDistribution::GetHalfZ ( ) const

Definition at line 216 of file G4SPSPosDistribution.cc.

 {
     return halfz;
 }

G4ThreeVector G4SPSPosDistribution::GetParticlePos ( ) const

Definition at line 241 of file G4SPSPosDistribution.cc.

                                                          {
     return ThreadData.Get().CParticlePos;
 }

Here is the call graph for this function:

G4String G4SPSPosDistribution::GetPosDisShape ( ) const

Definition at line 196 of file G4SPSPosDistribution.cc.

 {
     return Shape;
 }

Here is the caller graph for this function:

G4String G4SPSPosDistribution::GetPosDisType ( ) const

Definition at line 191 of file G4SPSPosDistribution.cc.

 {
     return SourcePosType;
 }

Here is the caller graph for this function:

G4double G4SPSPosDistribution::GetRadius ( ) const

Definition at line 221 of file G4SPSPosDistribution.cc.

 {
     return   Radius;
 }

G4ThreeVector G4SPSPosDistribution::GetSideRefVec1 ( ) const

Definition at line 245 of file G4SPSPosDistribution.cc.

                                                          {
     return ThreadData.Get().CSideRefVec1;
 }

Here is the call graph for this function:

G4ThreeVector G4SPSPosDistribution::GetSideRefVec2 ( ) const

Definition at line 249 of file G4SPSPosDistribution.cc.

                                                          {
     return ThreadData.Get().CSideRefVec2;
 }

Here is the call graph for this function:

G4ThreeVector G4SPSPosDistribution::GetSideRefVec3 ( ) const

Definition at line 253 of file G4SPSPosDistribution.cc.

                                                          {
   return ThreadData.Get().CSideRefVec3;
 }

Here is the call graph for this function:

G4String G4SPSPosDistribution::GetSourcePosType ( ) const

Definition at line 237 of file G4SPSPosDistribution.cc.

                                                       {
     return SourcePosType;
 }

void G4SPSPosDistribution::SetBeamSigmaInR ( G4double r )

Definition at line 160 of file G4SPSPosDistribution.cc.

 {
   SX = SY = r;
   SR = r;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetBeamSigmaInX ( G4double r )

Definition at line 166 of file G4SPSPosDistribution.cc.

 {
   SX = r;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetBeamSigmaInY ( G4double r )

Definition at line 171 of file G4SPSPosDistribution.cc.

 {
   SY = r;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetBiasRndm ( G4SPSRandomGenerator * a )

Definition at line 226 of file G4SPSPosDistribution.cc.

 {
     G4AutoLock l(&a_mutex);
     PosRndm = a ;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetCentreCoords ( G4ThreeVector coordsOfCentre )

Definition at line 107 of file G4SPSPosDistribution.cc.

 {
   CentreCoords = coordsOfCentre;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetHalfX ( G4double xhalf )

Definition at line 135 of file G4SPSPosDistribution.cc.

 {
   halfx = xhalf;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetHalfY ( G4double yhalf )

Definition at line 140 of file G4SPSPosDistribution.cc.

 {
   halfy = yhalf;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetHalfZ ( G4double zhalf )

Definition at line 145 of file G4SPSPosDistribution.cc.

 {
   halfz = zhalf;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetParAlpha ( G4double paralp )

Definition at line 176 of file G4SPSPosDistribution.cc.

 {
   ParAlpha = paralp;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetParPhi ( G4double parphi )

Definition at line 186 of file G4SPSPosDistribution.cc.

 {
   ParPhi = parphi;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetParTheta ( G4double parthe )

Definition at line 181 of file G4SPSPosDistribution.cc.

 {
   ParTheta = parthe;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetPosDisShape ( G4String shapeType )

Definition at line 102 of file G4SPSPosDistribution.cc.

 {
      Shape = shapeType;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetPosDisType ( G4String PosType )

Important: This is a shared class between threads. Only one thread should use the set-methods here. Note that this is achieved by UI commands. If you use these 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 97 of file G4SPSPosDistribution.cc.

 {
     SourcePosType = PosType;
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetPosRot1 ( G4ThreeVector posrot1 )

Definition at line 112 of file G4SPSPosDistribution.cc.

 {
   // This should be x'
   Rotx = posrot1;
   if(verbosityLevel == 2)
     {
       G4cout << "Vector x' " << Rotx << G4endl;
     }
   GenerateRotationMatrices();
 }

Here is the caller graph for this function:

void G4SPSPosDistribution::SetPosRot2 ( G4ThreeVector posrot2 )

Definition at line 123 of file G4SPSPosDistribution.cc.

 {
   // This is a vector in the plane x'y' but need not
   // be y'
   Roty = posrot2;
   if(verbosityLevel == 2)
     {
       G4cout << "The vector in the x'-y' plane " << Roty << G4endl;
     }
   GenerateRotationMatrices();
 }