Belle II Software  release-08-01-10
SensitiveDetector Class Reference

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

#include <SensitiveDetector.h>

Inheritance diagram for SensitiveDetector:
Collaboration diagram for SensitiveDetector:

Public Member Functions

 SensitiveDetector ()
 Constructor.
 
int saveG4TrackInfo (int trackID, int PDG, float Mass, float Energy, float vtx[3], float mom[3])
 Save saveG4TrackInfo into datastore.
 

Static Public Member Functions

static const std::map< std::string, RelationArray::EConsolidationAction > & getMCParticleRelations ()
 Return a list of all registered Relations with MCParticles.
 
static void setActive (bool activeStatus)
 Enable/Disable all Sensitive Detectors. More...
 
static void registerMCParticleRelation (const std::string &name, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
 Register an relation involving MCParticles. More...
 
static void registerMCParticleRelation (const RelationArray &relation, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
 Overload to make it easer to register MCParticle relations. More...
 

Protected Member Functions

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

Protected Attributes

int m_trackID
 track id
 

Private Member Functions

virtual bool ProcessHits (G4Step *aStep, G4TouchableHistory *aROhist)
 Check if recording hits is enabled and if so call step() and set the correct MCParticle flag. More...
 

Private Attributes

Const::EDetector m_subdetector
 Subdetector the class belongs to.
 

Static Private Attributes

static std::map< std::string, RelationArray::EConsolidationActions_mcRelations
 Static set holding all relations which have to be updated at the end of the Event.
 
static bool s_active
 Static bool which indicates wether recording of hits is enabled.
 

Detailed Description

Sensitive Detector implementation of the HE3TUBE detector.

Definition at line 22 of file SensitiveDetector.h.

Member Function Documentation

◆ ProcessHits()

bool ProcessHits ( G4Step *  aStep,
G4TouchableHistory *  aROhist 
)
inlineprivatevirtualinherited

Check if recording hits is enabled and if so call step() and set the correct MCParticle flag.

Called by Geant4 for each step inside the sensitive volumes attached

Definition at line 94 of file SensitiveDetectorBase.h.

◆ registerMCParticleRelation() [1/2]

static void registerMCParticleRelation ( const RelationArray relation,
RelationArray::EConsolidationAction  ignoreAction = RelationArray::c_negativeWeight 
)
inlinestaticinherited

Overload to make it easer to register MCParticle relations.

Parameters
relationRelationArray to register
ignoreAction

Definition at line 66 of file SensitiveDetectorBase.h.

◆ registerMCParticleRelation() [2/2]

void registerMCParticleRelation ( const std::string &  name,
RelationArray::EConsolidationAction  ignoreAction = RelationArray::c_negativeWeight 
)
staticinherited

Register an relation involving MCParticles.

All Relations which point from an MCParticle to something have to be registered with addMCParticleRelation() because the index of the MCParticles might change at the end of the event. During simulation, the TrackID should be used as index of the MCParticle

Parameters
nameName of the relation to register
ignoreAction

Definition at line 22 of file SensitiveDetectorBase.cc.

◆ setActive()

static void setActive ( bool  activeStatus)
inlinestaticinherited

Enable/Disable all Sensitive Detectors.

By default, all sensitive detectors won't create hits to make it possible to use the Geant4 Navigator for non-simulation purposes. Only during simulation the sensitive detectors will be enabled to record hits

Parameters
activeStatusbool to indicate wether hits should be recorded

Definition at line 50 of file SensitiveDetectorBase.h.

◆ step()

bool step ( G4Step *  step,
G4TouchableHistory *   
)
overrideprotectedvirtual

Step processing method.

Parameters
stepthe G4Step with the current step information
Returns
true if a Hit has been created, false if the hit was ignored

Implements SensitiveDetectorBase.

Definition at line 51 of file SensitiveDetector.cc.

52  {
53  //Get Track information
54  const G4Track& track = *step->GetTrack();
55  const int trackID = track.GetTrackID();
56  const double depEnergy = step->GetTotalEnergyDeposit() * CLHEP::MeV;
57  const double nielEnergy = step->GetNonIonizingEnergyDeposit() * CLHEP::MeV;
58  const G4ThreeVector G4tkPos = step->GetTrack()->GetPosition();
59  float tkPos[3];
60  tkPos[0] = G4tkPos.x() * CLHEP::cm;
61  tkPos[1] = G4tkPos.y() * CLHEP::cm;
62  tkPos[2] = G4tkPos.z() * CLHEP::cm;
63  const G4ThreeVector G4tkMom = step->GetTrack()->GetMomentum();
64  float tkMom[3];
65  tkMom[0] = G4tkMom.x() * CLHEP::MeV;
66  tkMom[1] = G4tkMom.y() * CLHEP::MeV;
67  tkMom[2] = G4tkMom.z() * CLHEP::MeV;
68  const G4ThreeVector G4tkMomDir = step->GetTrack()->GetMomentumDirection();
69  float tkMomDir[3];
70  tkMomDir[0] = G4tkMomDir.x() * CLHEP::MeV;
71  tkMomDir[1] = G4tkMomDir.y() * CLHEP::MeV;
72  tkMomDir[2] = G4tkMomDir.z() * CLHEP::MeV;
73  const int tkPDG = step->GetTrack()->GetDefinition()->GetPDGEncoding();
74  const double tkKEnergy = step->GetTrack()->GetKineticEnergy();
75  const int detNb = step->GetTrack()->GetVolume()->GetCopyNo();
76  const double GlTime = step->GetPreStepPoint()->GetGlobalTime();
77  //Ignore everything below 1eV
78  //if (depEnergy < CLHEP::eV) return false;
79 
80  //Get the datastore arrays
81  StoreArray<MCParticle> mcParticles;
82  StoreArray<He3tubeSimHit> simHits;
83  RelationArray relMCSimHit(mcParticles, simHits);
84 
85  //find out if the process that created the particle was a neutron process
86  bool neuProc = false;
87  G4String CPName;
88  if (step->GetTrack()->GetCreatorProcess() != 0) {
89  const G4VProcess* creator = step->GetTrack()->GetCreatorProcess();
90  CPName = creator->GetProcessName();
91  if (G4StrUtil::contains(CPName, "Neutron")) neuProc = true;
92  }
93 
94  //set TrackID
95  m_trackID = track.GetTrackID();
96  //Save Hit if track leaves volume or is killed
97  if (track.GetNextVolume() != track.GetVolume() || track.GetTrackStatus() >= fStopAndKill) {
98  auto& mcparticle = Simulation::TrackInfo::getInfo(track);
99  int PDG = mcparticle.getPDG();
100  float Mass = mcparticle.getMass();
101  float Energy = mcparticle.getEnergy();
102  float vtx[3];
103  vtx[0] = mcparticle.getProductionVertex().X();
104  vtx[1] = mcparticle.getProductionVertex().Y();
105  vtx[2] = mcparticle.getProductionVertex().Z();
106  float mom[3];
107  mom[0] = mcparticle.getMomentum().X();
108  mom[1] = mcparticle.getMomentum().Y();
109  mom[2] = mcparticle.getMomentum().Z();
110  saveG4TrackInfo(m_trackID, PDG, Mass, Energy, vtx, mom);
111  //Reset TrackID
112  m_trackID = 0;
113  }
114 
115  StoreArray<He3tubeSimHit> He3tubeHits;
116  He3tubeSimHit* hit = He3tubeHits.appendNew(
117  depEnergy,
118  nielEnergy,
119  tkPDG,
120  tkKEnergy,
121  detNb,
122  GlTime,
123  tkPos,
124  tkMom,
125  tkMomDir,
126  neuProc
127  );
128 
129  //Add Relation between SimHit and MCParticle with a weight of 1. Since
130  //the MCParticle index is not yet defined we use the trackID from Geant4
131  relMCSimHit.add(trackID, hit->getArrayIndex(), 1.0);
132 
133  return true;
134  }
static Payload getInfo(Carrier &obj)
Static function to just return UserInformation attached to the obj of type Carrier.
Definition: UserInfo.h:100
int saveG4TrackInfo(int trackID, int PDG, float Mass, float Energy, float vtx[3], float mom[3])
Save saveG4TrackInfo into datastore.
bool step(G4Step *step, G4TouchableHistory *) override
Step processing method.

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