Belle II Software  light-2403-persian
AWESOMESensitiveDetector Class Reference

Sensitive Detector implementation of the AWESOME detector. More...

#include <AWESOMESensitiveDetector.h>

Inheritance diagram for AWESOMESensitiveDetector:
Collaboration diagram for AWESOMESensitiveDetector:

Public Member Functions

 AWESOMESensitiveDetector ()
 Constructor.
 
bool step (G4Step *step, G4TouchableHistory *) override
 Step processing method. More...
 

Private Attributes

StoreArray< MCParticlem_MCParticles
 MC particles.
 
StoreArray< AWESOMESimHitm_SimHits
 AWESOME simulated hits.
 
RelationArray m_MCParticlesToSimHits {m_MCParticles, m_SimHits}
 Relation array between MCParticles and AWESOMESimHits.
 

Detailed Description

Sensitive Detector implementation of the AWESOME detector.

Definition at line 27 of file AWESOMESensitiveDetector.h.

Member Function Documentation

◆ step()

bool step ( G4Step *  step,
G4TouchableHistory *   
)
override

Step processing method.

Parameters
stepThe G4Step with the current step information.
Returns
true if a hit was created, false otherwise.

Definition at line 44 of file AWESOMESensitiveDetector.cc.

45 {
46  /* Get some basic information from Geant4. */
47  const G4Track& track = *step->GetTrack();
48  const int trackID = track.GetTrackID();
49  const double energyDep = step->GetTotalEnergyDeposit() * Unit::MeV;
50  G4StepPoint* preStep = step->GetPreStepPoint();
51  G4StepPoint* postStep = step->GetPostStepPoint();
52  /* Ignore everything below 1 eV. */
53  if (energyDep < Unit::eV)
54  return false;
55  /*
56  * Compute other useful quantities to be stored.
57  * We must be sure that the quantities are stored using the proper Belle2 units:
58  * positions in cm, time in ns, etc.
59  */
60  const CLHEP::Hep3Vector position = 0.5 * (preStep->GetPosition() + postStep->GetPosition()) / CLHEP::cm; // Now in cm
61  const double time = 0.5 * (preStep->GetGlobalTime() + postStep->GetGlobalTime()); // Already in ns
62  /* Store the simulated hit. */
63  AWESOMESimHit* simHit = m_SimHits.appendNew();
64  simHit->setEnergyDep(energyDep);
65  simHit->setPosition(B2Vector3{position.x(), position.y(), position.z()});
66  simHit->setTime(time);
67  /*
68  * Add a relation between the current MCParticle and the simulated hit.
69  * Since the MCParticle index is not yet defined we use the trackID from Geant4.
70  */
71  m_MCParticlesToSimHits.add(trackID, simHit->getArrayIndex());
72  return true;
73 }
A Geant4 simulated hit for the AWESOME detector.
Definition: AWESOMESimHit.h:33
void setEnergyDep(float energyDep)
Set the deposited energy.
Definition: AWESOMESimHit.h:71
void setTime(float time)
Set the time.
void setPosition(B2Vector3< float > position)
Set the vector for position.
Definition: AWESOMESimHit.h:89
StoreArray< AWESOMESimHit > m_SimHits
AWESOME simulated hits.
RelationArray m_MCParticlesToSimHits
Relation array between MCParticles and AWESOMESimHits.
bool step(G4Step *step, G4TouchableHistory *) override
Step processing method.
A fast and root compatible alternative to TVector3.
Definition: B2Vector3.h:42
DataType x() const
access variable X (= .at(0) without boundary check)
Definition: B2Vector3.h:425
void add(index_type from, index_type to, weight_type weight=1.0)
Add a new element to the relation.
int getArrayIndex() const
Returns this object's array index (in StoreArray), or -1 if not found.
static const double eV
[electronvolt]
Definition: Unit.h:112
static const double MeV
[megaelectronvolt]
Definition: Unit.h:114

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