Definition of GGS Integrated Hit. More...

#include <GGSIntHit.h>

Inheritance diagram for GGSIntHit:

Collaboration diagram for GGSIntHit:

Public Member Functions
	GGSIntHit (const std::vector< G4double > &timeBins)
	Constructor. More...

	~GGSIntHit ()
	Destructor.

	GGSIntHit (const GGSIntHit &right)
	Copy constructor. More...

const GGSIntHit &	operator= (const GGSIntHit &right)
	Assignment operator. More...

int	operator== (const GGSIntHit &right) const
	Comparison operator. More...

void *	operator new (size_t size)
	Operator new. More...

void	operator delete (void *aHit)
	Operator delete. More...

virtual void	AddStep (const G4Step &step)
	Adds a step to the particle hit. More...

void	SetVolume (const G4VPhysicalVolume *volume)
	Setter for volume. More...

void	SetPosHitsStorage (bool flag)
	Turn on or off the storage of position hits. More...

void	SetPartHitsStorage (bool flag)
	Turn on or off the storage of particle hits. More...

void	SetAbsTranslation (const G4ThreeVector &pos)
	Set the position in global coordinates of the sensitive element. More...

void	SetID (G4int id)
	Setter for volume ID. More...

G4double	GetTotalEnergyDep () const
	Energy release getter. More...

const std::vector< G4double > &	GetEnergyDep () const
	Energy release getter. More...

G4double	GetTime () const
	Time getter. More...

const std::vector< G4double > &	GetTimeBins () const
	Time bins getter. More...

G4double	GetPathLength () const
	Path length getter. More...

const G4VPhysicalVolume *	GetVolume ()
	Getter for volume. More...

GGSPartHitsCollection *	GetPartHits ()
	Getter of container of particle hits. More...

const G4ThreeVector &	GetAbsTranslation ()
	Getter of absolute position of sensitive element. More...

G4int	GetID ()
	Getter for volume ID. More...

Protected Attributes
std::vector< G4double >	_eDep

std::vector< G4double >	_timeBins

G4double	_time

G4double	_pathLength

const G4VPhysicalVolume *	_volume

G4ThreeVector	_absTranslation

G4int	_id

bool	_isReset

GGSPartHitsCollection *	_partHits

bool	_storePosHits

Detailed Description

Definition of GGS Integrated Hit.

An integrated hits is associated to a sensitive physical volume. It carries informations about the whole energy released on that volume by all the particles traversing it during an event.

Definition at line 31 of file GGSIntHit.h.

Constructor & Destructor Documentation

GGSIntHit::GGSIntHit ( const std::vector< G4double > & timeBins )

Constructor.

Parameters

timeBins The time bins in which integrated hits will be accumulated.

Definition at line 12 of file GGSIntHit.cpp.

                                                         :
     _eDep(timeBins.size() + 1, 0.), _timeBins(timeBins), _time(-1.), _pathLength(0.), _volume(NULL), _absTranslation(0.,
         0., 0.), _id(-1), _isReset(true), _partHits(
     NULL), _storePosHits(false) {
 }

GGSIntHit::GGSIntHit ( const GGSIntHit & right )

Copy constructor.

Parameters

right The object to be copied into current one.

Definition at line 26 of file GGSIntHit.cpp.

                                            :
     G4VHit(), _timeBins(right._timeBins) {
   _eDep = right._eDep;
   _time = right._time;
   _pathLength = right._pathLength;
   _volume = right._volume;
   _absTranslation = right._absTranslation;
   _id = right._id;
   _isReset = right._isReset;
 
   if (right._partHits == NULL) {
     _partHits = NULL;
   }
   else
     _partHits = new GGSPartHitsCollection(*right._partHits);
 
   _storePosHits = right._storePosHits;
 }

Member Function Documentation

void GGSIntHit::AddStep ( const G4Step & step )

virtual

Adds a step to the particle hit.

This method adds a step to the integrated hit. The integrated hit integral quantities (scoring value and path length) will be incremented by the corresponding step quantities.

If SetPartHitsStorage has been previously called to enable the storage of particle hits, then the step will be added also to its corresponding particle hit object. The same for SetPosHitsStorage relatively to storage of position hits.

Parameters

step	The step to be added.

Definition at line 104 of file GGSIntHit.cpp.

                                           {
   static const std::string routineName("GGSIntHit::AddStep");
   // Find current bin
   unsigned int currBin = 0;
   while (currBin < _eDep.size() - 1 && step.GetPreStepPoint()->GetGlobalTime() > _timeBins[currBin]) {
     currBin++;
   }
 
   // Increment energy deposit
   _eDep[currBin] += step.GetTotalEnergyDeposit();
   _pathLength += step.GetStepLength();
 
   if (_isReset) {
     // First step for this integrated hit: store initial values
     _isReset = false;
     _time = step.GetPreStepPoint()->GetGlobalTime();
   }
 }

const G4ThreeVector& GGSIntHit::GetAbsTranslation ( )

inline

Getter of absolute position of sensitive element.

Returns: The position of the sensitive element associated to the hit in global coordinates.

Definition at line 242 of file GGSIntHit.h.

                                            {
     return _absTranslation;
   }

const std::vector<G4double>& GGSIntHit::GetEnergyDep ( ) const

inline

Energy release getter.

Returns: The energy release for each time bin.

Definition at line 186 of file GGSIntHit.h.

                                                   {
     return _eDep;
   }

G4int GGSIntHit::GetID ( )

inline

Getter for volume ID.

Returns: The ID of the volume (i.e., touchable) associated to the hit.

Definition at line 250 of file GGSIntHit.h.

                 {
     return _id;
   }

GGSPartHitsCollection* GGSIntHit::GetPartHits ( )

inline

Getter of container of particle hits.

Returns: pointer to the particle hit container (NULL if storage of particle hits has not been previously turned on by calling SetPartHitsStorage).

Definition at line 234 of file GGSIntHit.h.

                                        {
     return _partHits;
   }

G4double GGSIntHit::GetPathLength ( ) const

inline

Path length getter.

Returns: The total path length of particles traversing the volume.

Definition at line 216 of file GGSIntHit.h.

                                  {
     return _pathLength;
   }

G4double GGSIntHit::GetTime ( ) const

inline

Time getter.

Returns: The time of first energy deposit.

Definition at line 194 of file GGSIntHit.h.

                            {
     return _time;
   }

const std::vector<G4double>& GGSIntHit::GetTimeBins ( ) const

inline

Time bins getter.

The returned vector contains the time bins. Element 0 is the ending time of first bin (which always starts at time 0); last element is the ending time of last bin but one, which corresponds to start time of last bin (last bin always ends at the end of event).

Returns: The time binning for energy deposit.

Definition at line 208 of file GGSIntHit.h.

                                                  {
     return _timeBins;
   }

G4double GGSIntHit::GetTotalEnergyDep ( ) const

inline

Energy release getter.

Returns: The total energy release.

Definition at line 178 of file GGSIntHit.h.

                                      {
     return accumulate(_eDep.begin(), _eDep.end(), 0.);
   }

const G4VPhysicalVolume* GGSIntHit::GetVolume ( )

inline

Getter for volume.

Returns: The physical volume associated to the hit.

Definition at line 224 of file GGSIntHit.h.

                                        {
     return _volume;
   }

void GGSIntHit::operator delete ( void * aHit )

inline

Operator delete.

This implementation uses G4Allocator.

Parameters

aHit	The current hit to be deleted.

See Also: GGSIntHitAllocator.

Definition at line 287 of file GGSIntHit.h.

                                                  {
   GGSIntHitAllocator->FreeSingle((GGSIntHit*) aHit);
 }

void * GGSIntHit::operator new ( size_t size )

inline

Operator new.

This implementation uses G4Allocator.

Parameters

size	Unused (needed by the interface).

See Also: GGSIntHitAllocator.

Definition at line 279 of file GGSIntHit.h.

                                            {
   void* aHit;
   aHit = (void*) GGSIntHitAllocator->MallocSingle();
   return aHit;
 }

const GGSIntHit & GGSIntHit::operator= ( const GGSIntHit & right )

Assignment operator.

Parameters

right The object to be assigned to current one.

Returns: The new current object.

Definition at line 47 of file GGSIntHit.cpp.

                                                             {
   if (this != &right) {
     _eDep = right._eDep;
     _timeBins = right._timeBins;
     _time = right._time;
     _pathLength = right._pathLength;
     _volume = right._volume;
     _absTranslation = right._absTranslation;
     _id = right._id;
     _isReset = right._isReset;
     if (right._partHits == NULL) {
       delete _partHits;
       _partHits = NULL;
     }
     else
       *_partHits = *(right._partHits);
     _storePosHits = right._storePosHits;
   }
   return *this;
 }

int GGSIntHit::operator== ( const GGSIntHit & right ) const

Comparison operator.

Parameters

right RHS of comparison.

Returns: true if LHS==RHS.

Definition at line 70 of file GGSIntHit.cpp.

                                                       {
   if ((_eDep == right._eDep) && //
       (_timeBins == right._timeBins) && //
       (_time == right._time) && //
       (_pathLength == right._pathLength) && //
       (_volume == right._volume) && //
       (_absTranslation == right._absTranslation) && //
       (_id == right._id) && //
       (_isReset == right._isReset)) {
     if (_partHits != NULL) {
       if (right._partHits != NULL) {
         return (*_partHits == *(right._partHits));
       }
 
       else {
         return false;
       }
     }
     else {
       if (right._partHits != NULL) {
         return false;
       }
       else {
         return true;
       }
     }
   }
   else
     return false;
 
 }

void GGSIntHit::SetAbsTranslation ( const G4ThreeVector & pos )

inline

Set the position in global coordinates of the sensitive element.

Parameters

pos	Absolute position of the sensitive element associated to the hit.

Definition at line 136 of file GGSIntHit.h.

                                                    {
     _absTranslation = pos;
   }

void GGSIntHit::SetID ( G4int id )

inline

Setter for volume ID.

An integer identification code for the sensitive volume (i.e. touchable). For output on Root files with GGSRootFileService it is necessary that each volume have a unique (name, ID) pair. This ID must then be set appropriately in order to have a unique (name, ID) for each sensitive volume.

Replicated volumes as G4PVReplica produce many touchables (sensitive volumes), each with a different replica number. However, if a G4PVReplica is placed multiple times (via multiple placements of a G4PVPlacement whose associated logical volume is filled by the G4PVReplica) then all the touchables will have the name of the G4PVReplica, but they will not have different IDs since touchable 0 will have replica number equal to 0 in all its placements. Since it is not possible to manually set a different replica number for each touchable, the absolute ID must be computed in order to be unique for each touchable. As an example, the current implementation in GGSIntHitSD::ProcessHits uses the copy number of the mother G4PVPlacement in order to compute a unique ID for each touchable as:

ID = touchableCopyNumber + touchableMultiplicity*motherCopyNumber

In this way, a replicate volume with 10 touchables placed many times will generate ID=0 for the first placement of the first touchable and ID=10 for the second placement of the first touchable, and so on (as long as each G4PVPlacement have different copy number; this must be ensured by the user in the implementation of the geometry). This numbering scheme for replicated volumes may be reused by other sensitive detectors, although it is not mandatory as long as each sensitive touchable has a unique (name, ID) pair.

It is responsibility of the sensitive detector class to properly fill this field.

Parameters

ID	The ID of the volume (i.e., touchable) associated to the hit.

Definition at line 170 of file GGSIntHit.h.

                        {
     _id = id;
   }

void GGSIntHit::SetPartHitsStorage ( bool flag )

Turn on or off the storage of particle hits.

Steps added to integrated hit by AddStep are by default not stored. This method enables or disables step storage as GGSPartHit objects (i.e., integration of steps over the sensitive volume for each particle).

Calling this method with "false" parameter will automatically disable position hits storage, i.e., it will call SetPosHitsStorage(false).

Parameters

flag	if true, steps will be stored as particle hits.

Definition at line 144 of file GGSIntHit.cpp.

                                             {
 
   if (flag == true) {
     if (!_partHits) {
       _partHits = new GGSPartHitsCollection;
     }
   }
   else {
     _storePosHits = false;
     delete _partHits;
     _partHits = NULL;
   }
 }

void GGSIntHit::SetPosHitsStorage ( bool flag )

Turn on or off the storage of position hits.

Steps can be stored as position hits. If this method is called with "true" parameter, each particle hit associated to the integrated hit will contain its related position hits.

Calling this method with "true" parameter will automatically enable particle hits storage, i.e., it will call SetPartHitsStorage(true).

Parameters

flag	if true, steps will be stored as position hits.

See Also: SetPartHitsStorage; GGSPartHit::SetPosHitsStorage

Definition at line 125 of file GGSIntHit.cpp.

                                            {
 
   if (flag == true) {
     _storePosHits = true;
     SetPartHitsStorage(true);
 
   }
   else {
     _storePosHits = false;
     if (_partHits) {
       for (G4int i = 0; i < _partHits->entries(); i++) {
         ((*_partHits)[i])->SetPosHitsStorage(false);
       }
     }
   }
 }

void GGSIntHit::SetVolume ( const G4VPhysicalVolume * volume )

inline

Setter for volume.

Parameters

name	The physical volume associated to the hit.

Definition at line 99 of file GGSIntHit.h.

                                                   {
     _volume = volume;
   }

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

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation