EDepSim::DensityPositionGenerator Class Reference

#include <EDepSimDensityPositionGenerator.hh>

Inheritance diagram for EDepSim::DensityPositionGenerator:

Public Member Functions
	DensityPositionGenerator (const G4String &name)
virtual	~DensityPositionGenerator ()
virtual G4LorentzVector	GetPosition ()
	Return a candidate vertex. More...
virtual bool	ForcePosition ()
Public Member Functions inherited from EDepSim::VConstrainedPositionGenerator
	VConstrainedPositionGenerator (const G4String &name)
virtual	~VConstrainedPositionGenerator ()
void	SetVolumeName (const G4String &volume)
	Set the name of the volume to be sampled for a vertex. More...
void	ClearPositionTests (void)
	Clear the current set of vertex checks. More...
void	CheckVolumeName (const G4String &name)
void	CheckNotVolumeName (const G4String &name)
void	CheckVolumeMaterial (const G4String &name)
	Check that the vertex is inside of a material specified by name. More...
void	CheckNotVolumeMaterial (const G4String &name)
	Check that the vertex is not inside of a material specified by name. More...
void	CheckMinX (double x)
	Check that the vertex X position is greater than some value. More...
void	CheckMaxX (double x)
	Check that the vertex X position is less than some value. More...
void	CheckMinY (double y)
	Check that the vertex Y position is greater than some value. More...
void	CheckMaxY (double y)
	Check that the vertex Y position is less than some value. More...
void	CheckMinZ (double z)
	Check that the vertex Z position is greater than some value. More...
void	CheckMaxZ (double z)
	Check that the vertex Z position is less than some value. More...
Public Member Functions inherited from EDepSim::VPositionGenerator
	VPositionGenerator (const G4String &name)
virtual	~VPositionGenerator ()
G4String	GetName () const
	Return the name of the generator. More...

Private Attributes
double	fMaximumDensity
	The maximum density in the detector. More...

Additional Inherited Members
Public Types inherited from EDepSim::VConstrainedPositionGenerator
typedef std::vector< PositionTest * >	PositionTests
Protected Member Functions inherited from EDepSim::VConstrainedPositionGenerator
virtual bool	ValidPosition (const G4LorentzVector &vtx)
	Return true if the vertex is valid. This is used in the derived class. More...
G4LorentzVector	TrialPosition ()
	Generate a trial position uniformly in the sample box. More...

Detailed Description

Definition at line 8 of file EDepSimDensityPositionGenerator.hh.

Constructor & Destructor Documentation

EDepSim::DensityPositionGenerator::DensityPositionGenerator

(

const G4String &

name

)

Definition at line 13 of file EDepSimDensityPositionGenerator.cc.

    : EDepSim::VConstrainedPositionGenerator(name), fMaximumDensity(-1) {}

EDepSim::DensityPositionGenerator::~DensityPositionGenerator ( )

virtual

Definition at line 16 of file EDepSimDensityPositionGenerator.cc.

{}

Member Function Documentation

bool EDepSim::DensityPositionGenerator::ForcePosition ( )

virtual

Flag if the vertex should be forced to the candidate vertex returned by GetPosition().

Reimplemented from EDepSim::VConstrainedPositionGenerator.

Definition at line 55 of file EDepSimDensityPositionGenerator.cc.

                                                    {
    return true;
}

G4LorentzVector EDepSim::DensityPositionGenerator::GetPosition ( )

virtual

Return a candidate vertex.

Implements EDepSim::VPositionGenerator.

Definition at line 18 of file EDepSimDensityPositionGenerator.cc.

                                                             {
    if (fMaximumDensity < 0) {
        const G4MaterialTable* mTable = G4Material::GetMaterialTable();
        for (G4MaterialTable::const_iterator m = mTable->begin();
             m != mTable->end();
             ++m) {
            fMaximumDensity = std::max(fMaximumDensity, (*m)->GetDensity());
        }
    }
    // Get the navigator.
    G4Navigator* navigator 
        = G4TransportationManager::GetTransportationManager()
        ->GetNavigatorForTracking();
    for (int finiteLoop = 0; finiteLoop<100000; ++finiteLoop) {
        G4LorentzVector vtx = TrialPosition();

        // Get the volume that contains the point.
        G4VPhysicalVolume* volume = navigator->LocateGlobalPointAndSetup(vtx);

        // Skip this vertex if it's not in a legal volume.
        if (!volume) continue;
            
        // Get the density.
        double density
            = volume->GetLogicalVolume()->GetMaterial()->GetDensity();

        // Skip this vertex if it misses the sampling.
        if (density < fMaximumDensity*G4UniformRand()) continue;

        // The vertex meets the density sampling, so check if it is valid.
        if (ValidPosition(vtx)) return vtx;
    }
    EDepSimError("EDepSim::DensityPositionGenerator::GetPosition:"
              << " No valid position found");
    throw;
}

Member Data Documentation

double EDepSim::DensityPositionGenerator::fMaximumDensity

private

The maximum density in the detector.

Definition at line 22 of file EDepSimDensityPositionGenerator.hh.

---

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

• edep-sim/src/kinem/EDepSimDensityPositionGenerator.hh
• edep-sim/src/kinem/EDepSimDensityPositionGenerator.cc