G4AdjointPosOnPhysVolGenerator Class Reference

#include <G4AdjointPosOnPhysVolGenerator.hh>

Collaboration diagram for G4AdjointPosOnPhysVolGenerator:

Public Member Functions
G4VPhysicalVolume *	DefinePhysicalVolume (const G4String &aName)
void	DefinePhysicalVolume1 (const G4String &aName)
G4double	ComputeAreaOfExtSurface ()
G4double	ComputeAreaOfExtSurface (G4int NStat)
G4double	ComputeAreaOfExtSurface (G4double epsilon)
G4double	ComputeAreaOfExtSurface (G4VSolid *aSolid)
G4double	ComputeAreaOfExtSurface (G4VSolid *aSolid, G4int NStat)
G4double	ComputeAreaOfExtSurface (G4VSolid *aSolid, G4double epsilon)
void	GenerateAPositionOnTheExtSurfaceOfASolid (G4VSolid *aSolid, G4ThreeVector &p, G4ThreeVector &direction)
void	GenerateAPositionOnTheExtSurfaceOfTheSolid (G4ThreeVector &p, G4ThreeVector &direction)
void	GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume (G4ThreeVector &p, G4ThreeVector &direction)
void	GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume (G4ThreeVector &p, G4ThreeVector &direction, G4double &costh_to_normal)
void	SetSolid (G4VSolid *aSolid)
G4double	GetAreaOfExtSurfaceOfThePhysicalVolume ()
G4double	GetCosThDirComparedToNormal ()

Static Public Member Functions
static G4AdjointPosOnPhysVolGenerator *	GetInstance ()

Private Member Functions
	G4AdjointPosOnPhysVolGenerator ()
	~G4AdjointPosOnPhysVolGenerator ()
G4double	ComputeAreaOfExtSurfaceStartingFromSphere (G4VSolid *aSolid, G4int NStat)
G4double	ComputeAreaOfExtSurfaceStartingFromBox (G4VSolid *aSolid, G4int NStat)
void	GenerateAPositionOnASolidBoundary (G4VSolid *aSolid, G4ThreeVector &p, G4ThreeVector &direction)
G4double	GenerateAPositionOnASphereBoundary (G4VSolid *aSolid, G4ThreeVector &p, G4ThreeVector &direction)
G4double	GenerateAPositionOnABoxBoundary (G4VSolid *aSolid, G4ThreeVector &p, G4ThreeVector &direction)
void	ComputeTransformationFromPhysVolToWorld ()

Private Attributes
G4VSolid *	theSolid
G4VPhysicalVolume *	thePhysicalVolume
G4bool	UseSphere
G4String	ModelOfSurfaceSource
G4AffineTransform	theTransformationFromPhysVolToWorld
G4double	AreaOfExtSurfaceOfThePhysicalVolume
G4double	CosThDirComparedToNormal

Static Private Attributes
static G4ThreadLocal G4AdjointPosOnPhysVolGenerator *	theInstance = 0

Detailed Description

Definition at line 62 of file G4AdjointPosOnPhysVolGenerator.hh.

Constructor & Destructor Documentation

G4AdjointPosOnPhysVolGenerator()

G4AdjointPosOnPhysVolGenerator::G4AdjointPosOnPhysVolGenerator ( )

private

Definition at line 66 of file G4AdjointPosOnPhysVolGenerator.cc.

   : theSolid(0), thePhysicalVolume(0),
     UseSphere(true), ModelOfSurfaceSource("OnSolid"),
     AreaOfExtSurfaceOfThePhysicalVolume(0.), CosThDirComparedToNormal(0.)
{ 
}

Here is the caller graph for this function:

~G4AdjointPosOnPhysVolGenerator()

G4AdjointPosOnPhysVolGenerator::~G4AdjointPosOnPhysVolGenerator ( )

private

Definition at line 60 of file G4AdjointPosOnPhysVolGenerator.cc.

{ 
}

Here is the caller graph for this function:

Member Function Documentation

ComputeAreaOfExtSurface() [1/6]

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurface

(

)

Definition at line 109 of file G4AdjointPosOnPhysVolGenerator.cc.

{
   return ComputeAreaOfExtSurface(theSolid); 
}

Here is the caller graph for this function:

ComputeAreaOfExtSurface() [2/6]

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurface

(

G4int

NStat

)

Definition at line 115 of file G4AdjointPosOnPhysVolGenerator.cc.

{
   return ComputeAreaOfExtSurface(theSolid,NStats); 
}

Here is the call graph for this function:

ComputeAreaOfExtSurface() [3/6]

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurface

(

G4double

epsilon

)

Definition at line 121 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  return ComputeAreaOfExtSurface(theSolid,eps); 
}

Here is the call graph for this function:

ComputeAreaOfExtSurface() [4/6]

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurface

(

G4VSolid *

aSolid

)

Definition at line 127 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  return ComputeAreaOfExtSurface(aSolid,1.e-3); 
}

Here is the call graph for this function:

ComputeAreaOfExtSurface() [5/6]

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurface	(	G4VSolid *	aSolid,
		G4int	NStat
	)

Definition at line 133 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  if (ModelOfSurfaceSource == "OnSolid" ){
    if (UseSphere){
        return ComputeAreaOfExtSurfaceStartingFromSphere(aSolid,NStats);    
    }
    else {
        return ComputeAreaOfExtSurfaceStartingFromBox(aSolid,NStats);
    }
  }
  else {
    G4ThreeVector p,dir;
    if (ModelOfSurfaceSource == "ExternalSphere" ) return GenerateAPositionOnASphereBoundary(aSolid, p,dir);
    return GenerateAPositionOnABoxBoundary(aSolid, p,dir);
  }
}

Here is the call graph for this function:

ComputeAreaOfExtSurface() [6/6]

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurface	(	G4VSolid *	aSolid,
		G4double	epsilon
	)

Definition at line 151 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  G4int Nstats = G4int(1./(eps*eps));
  return ComputeAreaOfExtSurface(aSolid,Nstats);
}

Here is the call graph for this function:

ComputeAreaOfExtSurfaceStartingFromBox()

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurfaceStartingFromBox ( G4VSolid *  aSolid, G4int  NStat  )

private

Definition at line 177 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  if ( Nstat <= 0 ) return 0.;
  G4double area=1.;
  G4int i=0;
  G4int j=0;
  while (i<Nstat){
    G4ThreeVector p, direction;
    area = GenerateAPositionOnABoxBoundary( aSolid,p, direction);
    G4double dist_to_in = aSolid->DistanceToIn(p,direction);
    if (dist_to_in<kInfinity/2.) i++;
    j++;
 }
 area=area*double(i)/double(j);
 return area;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeAreaOfExtSurfaceStartingFromSphere()

G4double G4AdjointPosOnPhysVolGenerator::ComputeAreaOfExtSurfaceStartingFromSphere ( G4VSolid *  aSolid, G4int  NStat  )

private

Definition at line 195 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  if ( Nstat <= 0 ) return 0.;
  G4double area=1.;
  G4int i=0;
  G4int j=0;
  while (i<Nstat){
    G4ThreeVector p, direction;
    area = GenerateAPositionOnASphereBoundary( aSolid,p, direction);
    G4double dist_to_in = aSolid->DistanceToIn(p,direction);
    if (dist_to_in<kInfinity/2.) i++;
    j++;
 }
 area=area*double(i)/double(j);
 return area;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ComputeTransformationFromPhysVolToWorld()

void G4AdjointPosOnPhysVolGenerator::ComputeTransformationFromPhysVolToWorld ( )

private

Definition at line 390 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  G4VPhysicalVolume* daughter =thePhysicalVolume;
  G4LogicalVolume* mother = thePhysicalVolume->GetMotherLogical();
  theTransformationFromPhysVolToWorld = G4AffineTransform();
  G4PhysicalVolumeStore* thePhysVolStore =G4PhysicalVolumeStore::GetInstance();
  while (mother){
    theTransformationFromPhysVolToWorld *=
    G4AffineTransform(daughter->GetFrameRotation(),daughter->GetObjectTranslation());
    for ( unsigned int i=0; i< thePhysVolStore->size();i++){
        if ((*thePhysVolStore)[i]->GetLogicalVolume() == mother){
            daughter = (*thePhysVolStore)[i];
            mother =daughter->GetMotherLogical();
            break;
        };      
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefinePhysicalVolume()

G4VPhysicalVolume * G4AdjointPosOnPhysVolGenerator::DefinePhysicalVolume

(

const G4String &

aName

)

Definition at line 75 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  thePhysicalVolume = 0;
  theSolid =0;
  G4PhysicalVolumeStore* thePhysVolStore =G4PhysicalVolumeStore::GetInstance();
  for ( unsigned int i=0; i< thePhysVolStore->size();i++){
    G4String vol_name =(*thePhysVolStore)[i]->GetName();
    if (vol_name == ""){
        vol_name = (*thePhysVolStore)[i]->GetLogicalVolume()->GetName();
    }
    if (vol_name == aName){
        thePhysicalVolume = (*thePhysVolStore)[i];
    }
  }
  if (thePhysicalVolume){
    theSolid = thePhysicalVolume->GetLogicalVolume()->GetSolid();
    ComputeTransformationFromPhysVolToWorld();
    /*AreaOfExtSurfaceOfThePhysicalVolume=ComputeAreaOfExtSurface(1.e-3);
    G4cout<<"Monte Carlo  Estimate of the  area of the external surface :"<<AreaOfExtSurfaceOfThePhysicalVolume/m/m<<" m2"<<std::endl;*/
  }
  else {
    G4cout<<"The physical volume with name "<<aName<<" does not exist!!"<<std::endl;
    G4cout<<"Before generating a source on an external surface of a volume you should select another physical volume"<<std::endl; 
  }
  return thePhysicalVolume;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DefinePhysicalVolume1()

void G4AdjointPosOnPhysVolGenerator::DefinePhysicalVolume1

(

const G4String &

aName

)

Definition at line 103 of file G4AdjointPosOnPhysVolGenerator.cc.

{
   thePhysicalVolume = DefinePhysicalVolume(aName);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateAPositionOnABoxBoundary()

G4double G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnABoxBoundary ( G4VSolid *  aSolid, G4ThreeVector &  p, G4ThreeVector &  direction  )

private

Definition at line 268 of file G4AdjointPosOnPhysVolGenerator.cc.

{

  G4double ran_var,px,py,pz,minX,maxX,minY,maxY,minZ,maxZ;
  
  // values needed for CalculateExtent signature

  G4VoxelLimits limit;                // Unlimited
  G4AffineTransform origin;

  // min max extents of pSolid along X,Y,Z

  aSolid->CalculateExtent(kXAxis,limit,origin,minX,maxX);
  aSolid->CalculateExtent(kYAxis,limit,origin,minY,maxY);
  aSolid->CalculateExtent(kZAxis,limit,origin,minZ,maxZ);
  
  G4double scale=.1;
  minX-=scale*std::abs(minX);
  minY-=scale*std::abs(minY);
  minZ-=scale*std::abs(minZ);
  maxX+=scale*std::abs(maxX);
  maxY+=scale*std::abs(maxY);
  maxZ+=scale*std::abs(maxZ);
  
  G4double dX=(maxX-minX);
  G4double dY=(maxY-minY);
  G4double dZ=(maxZ-minZ);

  G4double XY_prob=2.*dX*dY;
  G4double YZ_prob=2.*dY*dZ;
  G4double ZX_prob=2.*dZ*dX;
  G4double area=XY_prob+YZ_prob+ZX_prob;
  XY_prob/=area;
  YZ_prob/=area;
  ZX_prob/=area;
  
  ran_var=G4UniformRand();
  G4double cos_th2 = G4UniformRand();
  G4double sth = std::sqrt(1.-cos_th2);
  G4double cth = std::sqrt(cos_th2);
  G4double phi=G4UniformRand()*3.1415926*2;
  G4double dirX = sth*std::cos(phi);
  G4double dirY = sth*std::sin(phi);
  G4double dirZ = cth;
  if (ran_var <=XY_prob){ //on the XY faces
    G4double ran_var1=ran_var/XY_prob;
    G4double ranX=ran_var1;
    if (ran_var1<=0.5){
        pz=minZ;
        direction=G4ThreeVector(dirX,dirY,dirZ);
        ranX=ran_var1*2.;
    } 
    else{
        pz=maxZ;
        direction=-G4ThreeVector(dirX,dirY,dirZ);
        ranX=(ran_var1-0.5)*2.;
    }
    G4double ranY=G4UniformRand();
    px=minX+(maxX-minX)*ranX;
    py=minY+(maxY-minY)*ranY;
  }
  else if (ran_var <=(XY_prob+YZ_prob)){ //on the YZ faces 
    G4double ran_var1=(ran_var-XY_prob)/YZ_prob;
    G4double ranY=ran_var1;
    if (ran_var1<=0.5){
        px=minX;
        direction=G4ThreeVector(dirZ,dirX,dirY);
        ranY=ran_var1*2.;
    } 
    else{
        px=maxX;
        direction=-G4ThreeVector(dirZ,dirX,dirY);
        ranY=(ran_var1-0.5)*2.;
    }
    G4double ranZ=G4UniformRand();
    py=minY+(maxY-minY)*ranY;
    pz=minZ+(maxZ-minZ)*ranZ;
        
  }
  else{ //on the ZX faces 
    G4double ran_var1=(ran_var-XY_prob-YZ_prob)/ZX_prob;
    G4double ranZ=ran_var1;
    if (ran_var1<=0.5){
        py=minY;
        direction=G4ThreeVector(dirY,dirZ,dirX);
        ranZ=ran_var1*2.;
    } 
    else{
        py=maxY;
        direction=-G4ThreeVector(dirY,dirZ,dirX);
        ranZ=(ran_var1-0.5)*2.;
    }
    G4double ranX=G4UniformRand();
    px=minX+(maxX-minX)*ranX;
    pz=minZ+(maxZ-minZ)*ranZ;
  }
  
  p=G4ThreeVector(px,py,pz);
  return area;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateAPositionOnASolidBoundary()

void G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnASolidBoundary ( G4VSolid *  aSolid, G4ThreeVector &  p, G4ThreeVector &  direction  )

private

Definition at line 213 of file G4AdjointPosOnPhysVolGenerator.cc.

{ 
  G4bool find_pos =false;
  while (!find_pos){
    if (UseSphere) GenerateAPositionOnASphereBoundary( aSolid,p, direction);
    else  GenerateAPositionOnABoxBoundary( aSolid,p, direction);
    G4double dist_to_in = aSolid->DistanceToIn(p,direction);
    if (dist_to_in<kInfinity/2.) {
      find_pos =true;
      p+= 0.999999*direction*dist_to_in;
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateAPositionOnASphereBoundary()

G4double G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnASphereBoundary ( G4VSolid *  aSolid, G4ThreeVector &  p, G4ThreeVector &  direction  )

private

Definition at line 228 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  G4double minX,maxX,minY,maxY,minZ,maxZ;

  // values needed for CalculateExtent signature

  G4VoxelLimits limit;                // Unlimited
  G4AffineTransform origin;

  // min max extents of pSolid along X,Y,Z

  aSolid->CalculateExtent(kXAxis,limit,origin,minX,maxX);
  aSolid->CalculateExtent(kYAxis,limit,origin,minY,maxY);
  aSolid->CalculateExtent(kZAxis,limit,origin,minZ,maxZ);

  G4ThreeVector center = G4ThreeVector((minX+maxX)/2.,(minY+maxY)/2.,(minZ+maxZ)/2.);
  
  G4double dX=(maxX-minX)/2.;
  G4double dY=(maxY-minY)/2.;
  G4double dZ=(maxZ-minZ)/2.;
  G4double scale=1.01;
  G4double r=scale*std::sqrt(dX*dX+dY*dY+dZ*dZ);

  G4double cos_th2 = G4UniformRand();
  G4double theta = std::acos(std::sqrt(cos_th2));
  G4double phi=G4UniformRand()*3.1415926*2;
  direction.setRThetaPhi(1.,theta,phi);
  direction=-direction;
  G4double cos_th = (1.-2.*G4UniformRand());
  theta = std::acos(cos_th);
  if (G4UniformRand() <0.5) theta=3.1415926-theta;
  phi=G4UniformRand()*3.1415926*2;
  p.setRThetaPhi(r,theta,phi);
  p+=center;
  direction.rotateY(theta);
  direction.rotateZ(phi);
  return 4.*3.1415926*r*r;;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateAPositionOnTheExtSurfaceOfASolid()

void G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnTheExtSurfaceOfASolid	(	G4VSolid *	aSolid,
		G4ThreeVector &	p,
		G4ThreeVector &	direction
	)

Definition at line 157 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  if (ModelOfSurfaceSource == "OnSolid" ){
    GenerateAPositionOnASolidBoundary(aSolid, p,direction);
        return;
  }
  if (ModelOfSurfaceSource == "ExternalSphere" ) {
    GenerateAPositionOnASphereBoundary(aSolid, p, direction);
    return;
  } 
    GenerateAPositionOnABoxBoundary(aSolid, p, direction);
    return;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume() [1/2]

void G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume	(	G4ThreeVector &	p,
		G4ThreeVector &	direction
	)

Definition at line 370 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  if (!thePhysicalVolume) {
    G4cout<<"Before generating a source on an external surface of volume you should select a physical volume"<<std::endl; 
    return;
  };
  GenerateAPositionOnTheExtSurfaceOfTheSolid(p,direction);
  p = theTransformationFromPhysVolToWorld.TransformPoint(p);
  direction = theTransformationFromPhysVolToWorld.TransformAxis(direction);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume() [2/2]

void G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume	(	G4ThreeVector &	p,
		G4ThreeVector &	direction,
		G4double &	costh_to_normal
	)

Definition at line 382 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  GenerateAPositionOnTheExtSurfaceOfThePhysicalVolume(p, direction);
  costh_to_normal = CosThDirComparedToNormal;
}                                       

Here is the call graph for this function:

GenerateAPositionOnTheExtSurfaceOfTheSolid()

void G4AdjointPosOnPhysVolGenerator::GenerateAPositionOnTheExtSurfaceOfTheSolid	(	G4ThreeVector &	p,
		G4ThreeVector &	direction
	)

Definition at line 171 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  GenerateAPositionOnTheExtSurfaceOfASolid(theSolid,p,direction);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetAreaOfExtSurfaceOfThePhysicalVolume()

G4double G4AdjointPosOnPhysVolGenerator::GetAreaOfExtSurfaceOfThePhysicalVolume ( )

inline

Definition at line 93 of file G4AdjointPosOnPhysVolGenerator.hh.

{return AreaOfExtSurfaceOfThePhysicalVolume;}

GetCosThDirComparedToNormal()

G4double G4AdjointPosOnPhysVolGenerator::GetCosThDirComparedToNormal ( )

inline

Definition at line 94 of file G4AdjointPosOnPhysVolGenerator.hh.

{return CosThDirComparedToNormal;}

Here is the call graph for this function:

GetInstance()

G4AdjointPosOnPhysVolGenerator * G4AdjointPosOnPhysVolGenerator::GetInstance ( void  )

static

Definition at line 49 of file G4AdjointPosOnPhysVolGenerator.cc.

{
  if(!theInstance)
  {
    theInstance = new  G4AdjointPosOnPhysVolGenerator;
  }
  return theInstance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetSolid()

void G4AdjointPosOnPhysVolGenerator::SetSolid ( G4VSolid *  aSolid )

inline

Definition at line 92 of file G4AdjointPosOnPhysVolGenerator.hh.

{theSolid=aSolid;}

Member Data Documentation

AreaOfExtSurfaceOfThePhysicalVolume

G4double G4AdjointPosOnPhysVolGenerator::AreaOfExtSurfaceOfThePhysicalVolume

private

Definition at line 118 of file G4AdjointPosOnPhysVolGenerator.hh.

CosThDirComparedToNormal

G4double G4AdjointPosOnPhysVolGenerator::CosThDirComparedToNormal

private

Definition at line 119 of file G4AdjointPosOnPhysVolGenerator.hh.

ModelOfSurfaceSource

G4String G4AdjointPosOnPhysVolGenerator::ModelOfSurfaceSource

private

Definition at line 116 of file G4AdjointPosOnPhysVolGenerator.hh.

theInstance

G4ThreadLocal G4AdjointPosOnPhysVolGenerator * G4AdjointPosOnPhysVolGenerator::theInstance = 0

staticprivate

Definition at line 111 of file G4AdjointPosOnPhysVolGenerator.hh.

thePhysicalVolume

G4VPhysicalVolume* G4AdjointPosOnPhysVolGenerator::thePhysicalVolume

private

Definition at line 113 of file G4AdjointPosOnPhysVolGenerator.hh.

theSolid

G4VSolid* G4AdjointPosOnPhysVolGenerator::theSolid

private

Definition at line 112 of file G4AdjointPosOnPhysVolGenerator.hh.

theTransformationFromPhysVolToWorld

G4AffineTransform G4AdjointPosOnPhysVolGenerator::theTransformationFromPhysVolToWorld

private

Definition at line 117 of file G4AdjointPosOnPhysVolGenerator.hh.

UseSphere

G4bool G4AdjointPosOnPhysVolGenerator::UseSphere

private

Definition at line 115 of file G4AdjointPosOnPhysVolGenerator.hh.

---

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

• G4AdjointPosOnPhysVolGenerator.hh
• G4AdjointPosOnPhysVolGenerator.cc