GGS(GenericGEANT4Simulation)Software  2.99.0
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Macros
Public Member Functions
GGSIntHitSD Class Reference

Sensitive detector class for integrated hits. More...

#include <GGSIntHitSD.h>

Inheritance diagram for GGSIntHitSD:
Inheritance graph
[legend]
Collaboration diagram for GGSIntHitSD:
Collaboration graph
[legend]

Public Member Functions

 GGSIntHitSD (G4String name)
 Constructor. More...
 
 ~GGSIntHitSD ()
 Destructor.
 
void Initialize (G4HCofThisEvent *hitCollection)
 Initializes the sensitive detector. More...
 
G4bool ProcessHits (G4Step *aStep, G4TouchableHistory *ROHist)
 The hit processing method. More...
 
void EndOfEvent (G4HCofThisEvent *hitCollection)
 Sets the hits collection, as required by specifications of mother class. More...
 
void SetTimeBin (G4double time)
 Time bins setter. More...
 
const std::vector< G4double > & GetTimeBins ()
 Time bins getter. More...
 
void SetPartHitsStorage (bool store)
 Turn on or off the storage (i.e. persistency) of particle hits. More...
 
void SetTouchableIDComputer (const std::string &tidcClassName)
 Sets the touchable ID computer class. More...
 
bool GetPartHitsStorage ()
 Return current status of particle hits storage (i.e. persistence). More...
 
void SetPosHitsStorage (bool store)
 Turn on or off the storage (i.e. persistence) of position hits. More...
 
bool GetPosHitsStorage ()
 Return current status of particle hits storage (i.e. persistence). More...
 
void SetIntHitThreshold (float minEnergy)
 Sets the minimum energy for integrated hits. More...
 
void SetPartHitThreshold (float minEnergy)
 Sets the minimum energy for particle hits. More...
 
void SetPosHitThreshold (float minEnergy)
 Sets the minimum energy for position hits. More...
 

Detailed Description

Sensitive detector class for integrated hits.

This sensitive detector takes care of summing up all energy releases inside the sensitive volume it manages, so that a single hit is produced for each single volume with non-zero energy release for each event.

Definition at line 29 of file GGSIntHitSD.h.

Constructor & Destructor Documentation

GGSIntHitSD::GGSIntHitSD ( G4String  name)

Constructor.

Parameters
nameThe sensitive detector's name.

Definition at line 21 of file GGSIntHitSD.cpp.

22  : G4VSensitiveDetector(name), _intHitCollection(NULL), _timeBins(0), _intHitClass("GGSIntHit"),
23  _partHitClass("GGSPartHit"), _posHitClass("GGSPosHit"), _storePartHits(false), _storePosHits(false),
24  _idComputer(nullptr), _intHitThreshold{0.}, _partHitThreshold{0.}, _posHitThreshold{0.}, _currPartHit{nullptr} {
25 
26  collectionName.insert(name);
27  _messenger = new GGSIntHitSDMessenger(this);
28 
29  auto params = name.substr(name.find_first_of('.') + 1); // Remove volume name
30  auto secondPointPos = params.find_first_of('.');
31  if (secondPointPos != std::string::npos) {
32  // Parameters have been specified so let's handle them
33  params = params.substr(secondPointPos + 1); // Remove GGSIntHitSD
34  if (params.size() > 0) {
35  // Retrieve hit classes
36  std::string intHitClassName, partHitClassName, posHitClassName;
37  auto start = params.find_first_of('[');
38  auto end = params.find_first_of(']');
39  if (start < std::string::npos && end != std::string::npos) {
40  intHitClassName = params.substr(start + 1, end - start - 1);
41  } else {
42  throw std::runtime_error(
43  std::string("[GGSIntHitSD::GGSIntHitSD] Missing definition of custom integrated hit class in \"") + params +
44  ("\""));
45  }
46 
47  start = params.find_first_of('[', end + 1);
48  end = params.find_first_of(']', start + 1);
49  if (start != std::string::npos && end != std::string::npos) {
50  partHitClassName = params.substr(start + 1, end - start - 1);
51  } else {
52  throw std::runtime_error(
53  std::string("[GGSIntHitSD::GGSIntHitSD] Missing definition of custom particle hit class in \"") + params +
54  ("\""));
55  }
56 
57  start = params.find_first_of('[', end + 1);
58  end = params.find_first_of(']', start + 1);
59  if (start != std::string::npos && end != std::string::npos) {
60  posHitClassName = params.substr(start + 1, end - start - 1);
61  } else {
62  throw std::runtime_error(
63  std::string("[GGSIntHitSD::GGSIntHitSD] Missing definition of custom position hit class in \"") + params +
64  ("\""));
65  }
66 
67  if (intHitClassName.size() > 0) {
68  auto &buildersList = GGSIntHitFactory::GetInstance().GetListOfRegisteredBuilders();
69  if (std::find(buildersList.begin(), buildersList.end(), intHitClassName) == buildersList.end()) {
70  throw std::runtime_error(std::string("[GGSIntHitSD::GGSIntHitSD] Integrated hit class \"") + intHitClassName +
71  "\" not found.");
72  }
73  _intHitClass = intHitClassName;
74  }
75  if (partHitClassName.size() > 0) {
76  auto &buildersList = GGSPartHitFactory::GetInstance().GetListOfRegisteredBuilders();
77  if (std::find(buildersList.begin(), buildersList.end(), partHitClassName) == buildersList.end()) {
78  throw std::runtime_error(std::string("[GGSIntHitSD::GGSIntHitSD] Particle hit class \"") + partHitClassName +
79  "\" not found.");
80  }
81  _partHitClass = partHitClassName;
82  }
83  if (posHitClassName.size() > 0) {
84  auto &buildersList = GGSPosHitFactory::GetInstance().GetListOfRegisteredBuilders();
85  if (std::find(buildersList.begin(), buildersList.end(), posHitClassName) == buildersList.end()) {
86  throw std::runtime_error(std::string("[GGSIntHitSD::GGSIntHitSD] Position hit class \"") + posHitClassName +
87  "\" not found.");
88  }
89  _posHitClass = posHitClassName;
90  }
91  }
92  }
93 
94  // Create the default ID computer
95  SetTouchableIDComputer("GGSUniqueTouchableIDComputer");
96 }
GGSFactory::GetListOfRegisteredBuilders
const std::vector< std::string > & GetListOfRegisteredBuilders()
Gets a vector containing the names of available registered builders.
Definition: GGSFactory.h:111
GGSFactory::GetInstance
static GGSFactory & GetInstance()
Getter method for singleton pointer.
Definition: GGSFactory.h:72
GGSIntHitSDMessenger
The integrated hit SD messenger class.
Definition: GGSIntHitSDMessenger.h:26
GGSIntHitSD::SetTouchableIDComputer
void SetTouchableIDComputer(const std::string &tidcClassName)
Sets the touchable ID computer class.
Definition: GGSIntHitSD.cpp:100

Member Function Documentation

void GGSIntHitSD::EndOfEvent ( G4HCofThisEvent *  hitCollection)

Sets the hits collection, as required by specifications of mother class.

Parameters
hitCollectionThe hits collection.

Definition at line 247 of file GGSIntHitSD.cpp.

247  {
248 
249  G4int IHID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
250  hitCollection->AddHitsCollection(IHID, _intHitCollection);
251 }
bool GGSIntHitSD::GetPartHitsStorage ( )
inline

Return current status of particle hits storage (i.e. persistence).

Returns
true if particle hits are going to be stored.

Definition at line 117 of file GGSIntHitSD.h.

117 { return _storePartHits; }
bool GGSIntHitSD::GetPosHitsStorage ( )
inline

Return current status of particle hits storage (i.e. persistence).

Returns
true if particle hits are going to be stored.

Definition at line 142 of file GGSIntHitSD.h.

142 { return _storePosHits; }
const std::vector<G4double>& GGSIntHitSD::GetTimeBins ( )
inline

Time bins getter.

First element is the end time of first bin (which always starts at t=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
Vector of time bins.

Definition at line 83 of file GGSIntHitSD.h.

83 { return _timeBins; }
void GGSIntHitSD::Initialize ( G4HCofThisEvent *  hitCollection)

Initializes the sensitive detector.

Parameters
hitCollectionUnused (needed by the interface).

Definition at line 107 of file GGSIntHitSD.cpp.

107  {
108 
109  // The hits collection is added to G4HCofThisEvent in EndOfEvent.
110  // It is automatically deleted at the end of the event so it must be recreated every time.
111  // However, G4Allocator should provide fast memory allocation for this purpose.
112  _intHitCollection = new GGSIntHitsCollection(SensitiveDetectorName, collectionName[0]);
113 }
GGSIntHitsCollection
G4THitsCollection< GGSIntHit > GGSIntHitsCollection
Alias for G4 template hits collection for GGSIntHit.
Definition: GGSIntHit.h:256
G4bool GGSIntHitSD::ProcessHits ( G4Step *  aStep,
G4TouchableHistory *  ROHist 
)

The hit processing method.

For each physical sensitive volume, this method sums up the total energy release during this step to the total releases of previous steps.

Parameters
aStepThe current step
ROHistUnused (needed by the interface).
Returns
true if energy deposit of current step is greater than 0.

Definition at line 115 of file GGSIntHitSD.cpp.

115  {
116  static const std::string routineName("GGSIntHitSD::ProcessHits");
117 
118  // The G4TouchableHistory argument of this method is non-null only if a readout geometry has been
119  // defined. Since this feature is not supported, it is currently always null.
120  // TODO: investigate the possibility to define a readout geometry
121 
122  if (aStep->GetTotalEnergyDeposit() == 0. && !_storePartHits)
123  return false;
124 
125  int touchableID = _idComputer->ComputeTouchableID(aStep);
126  G4TouchableHandle touchable = aStep->GetPreStepPoint()->GetTouchableHandle();
127  G4ThreeVector volPos = touchable->GetTranslation(); // Absolute volume translation
128  G4VPhysicalVolume *physVol = touchable->GetVolume();
129 
130  /* *** Add step to integrated hit *** */
131  GGSIntHit *intHit = nullptr;
132  bool intFound = false;
133  // Search for existing integrated hit for this volume
134  auto nIntHits = _intHitCollection->entries();
135  for (int iHit = nIntHits - 1; iHit >= 0; iHit--) {
136  intHit = (*(_intHitCollection))[iHit];
137  if (physVol == intHit->GetVolume() && touchableID == intHit->GetID() && volPos == intHit->GetAbsTranslation()) {
138  intFound = true;
139  break;
140  }
141  }
142  if (!intFound) {
143  // If we reach this point, then this is a new hit
144  intHit = GGSIntHitFactory::GetInstance().CreateObject(_intHitClass, _timeBins).release();
145  if (!intHit) {
146  throw std::runtime_error(
147  std::string("[GGSIntHitSD::ProcessHits] Can't create integrated hit of class ").append(_intHitClass));
148  }
149  // Set global hit properties
150  intHit->SetAbsTranslation(volPos);
151  intHit->SetID(touchableID);
152  intHit->SetVolume(physVol);
153  intHit->SetPartHitsStorage(_storePartHits);
154  intHit->SetPosHitsStorage(_storePosHits);
155  // Call the user initialization routine
156  intHit->UserInit(aStep);
157  // Insert the new integrated hit into the hit collection
158  _intHitCollection->insert(intHit);
159  }
160  // Add the current step to the hit
161  intHit->AddStep(*aStep);
162 
163  if (_storePartHits) {
164  /* *** Add step to particle hit *** */
165  // Create the particle hit if it does not exists, i.e. on first step inside the sensitive volume.
166  // If a particle entered in the volume in previous step and immediately got stopped but not killed then
167  // IsFirstStepInVolume for the current step will be true, so check also that kinetic energy is non-zero for pre-step
168  if (aStep->IsFirstStepInVolume() && aStep->GetPreStepPoint()->GetKineticEnergy() != 0) {
169  if (_currPartHit != nullptr) {
170  // On first step there must be no current hit
171  throw std::runtime_error("GGSIntHitSD::ProcessHit: existing particle hit on first step in sensitive volume");
172  }
173  _currPartHit = GGSPartHitFactory::GetInstance().CreateObject(_partHitClass).release();
174  if (!_currPartHit) {
175  throw std::runtime_error(
176  std::string("[GGSIntHitSD::ProcessHits] Can't create particle hit of class ").append(_partHitClass));
177  }
178  // Set global hit properties
179  _currPartHit->SetPosHitsStorage(_storePosHits);
180  }
181 
182  if (_currPartHit != nullptr) {
183  // Add the current step to the hit
184  _currPartHit->AddStep(*aStep);
185 
186  // Check if this is the last step in the volume
187  // Do not use aStep->IsLastStepInVolume(); see https://geant4-forum.web.cern.ch/t/last-step-in-a-volume/856
188  if ((aStep->GetPostStepPoint()->GetStepStatus() == G4StepStatus::fGeomBoundary &&
189  aStep->GetTrack()->GetTrackStatus() != G4TrackStatus::fStopButAlive) ||
190  aStep->GetTrack()->GetTrackStatus() == G4TrackStatus::fStopAndKill ||
191  aStep->GetTrack()->GetTrackStatus() == G4TrackStatus::fKillTrackAndSecondaries) {
192  if (_currPartHit->GetEnergyDeposit() > _partHitThreshold) {
193  // Hit above threshold: insert it into the hit collection
194  intHit->GetPartHits()->insert(_currPartHit);
195  } else {
196  // Hit below threshold: delete it
197  delete _currPartHit;
198  }
199  // Do not point at the hit anymore (at this point, it's either stored in the collection or deleted)
200  _currPartHit = nullptr;
201  }
202 
203  if (_storePosHits && aStep->GetTotalEnergyDeposit() > _posHitThreshold) {
204  /* *** Create position hit for current step *** */
205  GGSPosHit *posHit = GGSPosHitFactory::GetInstance().CreateObject(_posHitClass).release();
206  if (!posHit) {
207  throw std::runtime_error(
208  std::string("[GGSIntHitSD::ProcessHits] Can't create position hit of class ").append(_posHitClass));
209  }
210  // Insert the new position hit into the hit collection
211  auto *partHit =
212  (_currPartHit ? _currPartHit : (*(intHit->GetPartHits()))[intHit->GetPartHits()->GetSize() - 1]);
213  partHit->GetPosHits()->insert(posHit);
214  // Set members of the position hit
215  posHit->SetStep(*aStep);
216  }
217  } else {
218  std::vector<std::string> statusNames{
219  "fAlive", "fStopButAlive", "fStopAndKill", "fKillTrackAndSecondaries", "fSuspend", "fPostponeToNextEvent"};
220  GGSCOUT(WARNING) << "No particle hit has been created for particle:\n";
221  GGSCCOUT(WARNING) << " - type: "
222  << aStep->GetTrack()->GetDynamicParticle()->GetParticleDefinition()->GetParticleName() << "\n";
223  GGSCCOUT(WARNING) << " - creator process: " << aStep->GetTrack()->GetCreatorProcess()->GetProcessName() << "\n";
224  GGSCCOUT(WARNING) << " - trackID: " << aStep->GetTrack()->GetTrackID() << "\n";
225  GGSCCOUT(WARNING) << " - status: " << statusNames[static_cast<int>(aStep->GetTrack()->GetTrackStatus())] << "\n";
226  GGSCCOUT(WARNING) << " - eDep: " << aStep->GetTotalEnergyDeposit() / CLHEP::MeV << " MeV \n";
227  GGSCCOUT(WARNING) << " - kin. energy (pre) : " << aStep->GetPreStepPoint()->GetKineticEnergy() / CLHEP::MeV
228  << " MeV \n";
229  GGSCCOUT(WARNING) << " - kin. energy (post): " << aStep->GetPostStepPoint()->GetKineticEnergy() / CLHEP::MeV
230  << " MeV \n";
231  GGSCCOUT(WARNING) << " - process (pre) : "
232  << (aStep->GetPreStepPoint()->GetProcessDefinedStep()
233  ? aStep->GetPreStepPoint()->GetProcessDefinedStep()->GetProcessName()
234  : "none")
235  << "\n";
236  GGSCCOUT(WARNING) << " - process (post): "
237  << (aStep->GetPostStepPoint()->GetProcessDefinedStep()
238  ? aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName()
239  : "none")
240  << "\n";
241  }
242  }
243 
244  return true;
245 }
GGSCCOUT
#define GGSCCOUT(level)
Smart log utility which prints no header at the beginning of the line.
Definition: GGSSmartLog.h:126
GGSIntHit::AddStep
virtual void AddStep(const G4Step &step)
Adds a step to the particle hit.
Definition: GGSIntHit.cpp:93
GGSIntHit::GetAbsTranslation
const G4ThreeVector & GetAbsTranslation()
Getter of absolute position of sensitive element.
Definition: GGSIntHit.h:232
GGSPartHit::AddStep
virtual void AddStep(const G4Step &step)
Adds a step to the particle hit.
Definition: GGSPartHit.cpp:101
GGSPartHit::GetEnergyDeposit
G4double GetEnergyDeposit() const
Energy deposit getter.
Definition: GGSPartHit.h:148
GGSIntHit::GetPartHits
GGSPartHitsCollection * GetPartHits()
Getter of container of particle hits.
Definition: GGSIntHit.h:226
GGSFactory::CreateObject
std::unique_ptr< T > CreateObject(const std::string &name, ConstructorArgs...arguments)
Create an object.
Definition: GGSFactory.h:116
GGSIntHit::SetAbsTranslation
void SetAbsTranslation(const G4ThreeVector &pos)
Set the position in global coordinates of the sensitive element.
Definition: GGSIntHit.h:144
GGSIntHit::SetPosHitsStorage
void SetPosHitsStorage(bool flag)
Turn on or off the storage of position hits.
Definition: GGSIntHit.cpp:114
GGSIntHit::GetVolume
const G4VPhysicalVolume * GetVolume()
Getter for volume.
Definition: GGSIntHit.h:218
GGSPosHit
Definition of GGS Position Hit.
Definition: GGSPosHit.h:33
GGSIntHit::UserInit
virtual void UserInit(G4Step *aStep)
User initialization of hit global properties.
Definition: GGSIntHit.h:89
GGSPosHit::SetStep
virtual void SetStep(const G4Step &step)
Definition: GGSPosHit.cpp:60
GGSIntHit::GetID
G4int GetID()
Getter for volume ID.
Definition: GGSIntHit.h:238
GGSIntHit::SetID
void SetID(G4int id)
Setter for volume ID.
Definition: GGSIntHit.h:176
GGSPartHit::SetPosHitsStorage
void SetPosHitsStorage(bool flag)
Turn on or off the storage of position hits.
Definition: GGSPartHit.cpp:121
GGSIntHit::SetVolume
void SetVolume(const G4VPhysicalVolume *volume)
Setter for volume.
Definition: GGSIntHit.h:109
GGSIntHit
Definition of GGS Integrated Hit.
Definition: GGSIntHit.h:32
GGSFactory::GetInstance
static GGSFactory & GetInstance()
Getter method for singleton pointer.
Definition: GGSFactory.h:72
GGSIntHit::SetPartHitsStorage
void SetPartHitsStorage(bool flag)
Turn on or off the storage of particle hits.
Definition: GGSIntHit.cpp:132
void GGSIntHitSD::SetIntHitThreshold ( float  minEnergy)
inline

Sets the minimum energy for integrated hits.

Hits below this threshold will be removed from the hits array.

Parameters
minEnergyThe minimum energy for integrated hits.

Definition at line 150 of file GGSIntHitSD.h.

150 { _intHitThreshold = minEnergy; }
void GGSIntHitSD::SetPartHitsStorage ( bool  store)
inline

Turn on or off the storage (i.e. persistency) of particle hits.

This method tells to the SD whether to activate the storage of particle hits or not. By default, storage is not activated.

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

Parameters
storeif true, steps will be stored as position hits.
See Also
SetPosHitsStorage
GGSIntHit::SetPosHitsStorage

Definition at line 97 of file GGSIntHitSD.h.

97  {
98  _storePartHits = store;
99  if (store == false) {
100  SetPosHitsStorage(false);
101  }
102  }
GGSIntHitSD::SetPosHitsStorage
void SetPosHitsStorage(bool store)
Turn on or off the storage (i.e. persistence) of position hits.
Definition: GGSIntHitSD.h:131
void GGSIntHitSD::SetPartHitThreshold ( float  minEnergy)
inline

Sets the minimum energy for particle hits.

Hits below this threshold will be removed from the hits array.

Parameters
minEnergyThe minimum energy for particle hits.

Definition at line 158 of file GGSIntHitSD.h.

158 { _partHitThreshold = minEnergy; }
void GGSIntHitSD::SetPosHitsStorage ( bool  store)
inline

Turn on or off the storage (i.e. persistence) of position hits.

This method tells to the SD whether to activate the storage of position hits or not. By default, storage is not activated.

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

Parameters
storeif true, steps will be stored as position hits.
See Also
SetPartHitsStorage
GGSIntHit::SetPosHitsStorage

Definition at line 131 of file GGSIntHitSD.h.

131  {
132  _storePosHits = store;
133  if (store == true) {
134  SetPartHitsStorage(true);
135  }
136  }
GGSIntHitSD::SetPartHitsStorage
void SetPartHitsStorage(bool store)
Turn on or off the storage (i.e. persistency) of particle hits.
Definition: GGSIntHitSD.h:97
void GGSIntHitSD::SetPosHitThreshold ( float  minEnergy)
inline

Sets the minimum energy for position hits.

Hits below this threshold will be removed from the hits array.

Parameters
minEnergyThe minimum energy for position hits.

Definition at line 166 of file GGSIntHitSD.h.

166 { _posHitThreshold = minEnergy; }
void GGSIntHitSD::SetTimeBin ( G4double  time)
inline

Time bins setter.

This method adds a division in the timeline, effectively adding a time bin.

Parameters
timeThe time of first energy deposit.

Definition at line 70 of file GGSIntHitSD.h.

70  {
71  _timeBins.push_back(time);
72  std::sort(_timeBins.begin(), _timeBins.end());
73  }
void GGSIntHitSD::SetTouchableIDComputer ( const std::string &  tidcClassName)

Sets the touchable ID computer class.

This method sets the touchable ID computer class that will be used by the sensitive detector.

Parameters
tidcClassNamethe name of the touchable ID computer class.

Definition at line 100 of file GGSIntHitSD.cpp.

100  {
101  _idComputer = GGSFactory<GGSTouchableIDComputer>::GetInstance().CreateObject(tidcClassName);
102  if (!_idComputer) {
103  throw std::runtime_error("[GGSIntHitSD::GGSIntHitSD] Touchable ID computer class " + tidcClassName + " not found");
104  }
105 }
GGSFactory::CreateObject
std::unique_ptr< T > CreateObject(const std::string &name, ConstructorArgs...arguments)
Create an object.
Definition: GGSFactory.h:116
GGSFactory::GetInstance
static GGSFactory & GetInstance()
Getter method for singleton pointer.
Definition: GGSFactory.h:72
The documentation for this class was generated from the following files:
Generated on Mon May 13 2024 08:14:00 for GGS(GenericGEANT4Simulation)Software by   doxygen 1.8.5