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

#include <SensitiveDetector.h>

Inheritance diagram for SensitiveDetector:

Collaboration diagram for SensitiveDetector:

[legend]

Public Member Functions
	SensitiveDetector ()
	Constructor.

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...

Private Member Functions
bool	finishTrack ()
	finish a track

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
StoreArray< MCParticle >	m_mcParticles
	store array of the MCParticles

StoreArray< FANGSSimHit >	m_simHits
	store array of the SimHits

RelationArray	m_relMCSimHit {m_mcParticles, m_simHits}
	relation array of the MCParticle -> SimHit relation

std::stack< SensorTraversal >	m_tracks
	Stack of tracks to keep track of particles.

Const::EDetector	m_subdetector
	Subdetector the class belongs to.

Static Private Attributes
static std::map< std::string, RelationArray::EConsolidationAction >	s_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 FANGS detector.

Definition at line 28 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

relation	RelationArray 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

name	Name 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

activeStatus bool 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

step	the 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 45 of file SensitiveDetector.cc.

     {
       //Get scintilator ladder and sensor number
       const G4TouchableHistory* touchable = dynamic_cast<const G4TouchableHistory*>(step->GetPreStepPoint()->GetTouchable());
       int ladderID = touchable->GetVolume(2)->GetCopyNo();
       int sensorID = touchable->GetVolume(1)->GetCopyNo();
  
       //Get Track information
       const G4Track& track = *step->GetTrack();
       const int trackID = track.GetTrackID();
       const int pdgCode = step->GetTrack()->GetDefinition()->GetPDGEncoding();
       // deposited Energy in Geant4 Units
       double depEnergy  = step->GetTotalEnergyDeposit();
       // convert into Belle2 Units
       depEnergy *= Unit::GeV / CLHEP::GeV;
  
       // get pre and post step point
       const G4StepPoint& preStep = *step->GetPreStepPoint();
       const G4StepPoint& postStep = *step->GetPostStepPoint();
  
       // check if this is the same sensor traversal, otherwise add one to the stack
       if (m_tracks.empty() || (!m_tracks.top().check(trackID, ladderID, sensorID))) {
         m_tracks.push(SensorTraversal());
       }
       // get the top of the stack
       SensorTraversal& traversal = m_tracks.top();
  
       //If new track, remember values
       if (traversal.getTrackID() == 0) {
         bool isPrimary = Simulation::TrackInfo::getInfo(track).hasStatus(MCParticle::c_PrimaryParticle);
         //Add start position
         const G4ThreeVector preStepPos = preStep.GetPosition() / CLHEP::mm * Unit::mm;
         const G4ThreeVector preStepMom = preStep.GetMomentum() / CLHEP::MeV * Unit::MeV;
         //const G4AffineTransform& localToGlobalTransform = preStep.GetTouchableHandle()->GetHistory()->GetTopTransform().Inverse();
         const G4AffineTransform& localToGlobalTransform = preStep.GetTouchableHandle()->GetHistory()->GetTopTransform();
         const G4ThreeVector localpreStepPos = localToGlobalTransform.TransformPoint(preStep.GetPosition()) / CLHEP::mm * Unit::mm;
         const double time = preStep.GetGlobalTime() / CLHEP::ns * Unit::ns;
         traversal.setInitial(trackID, ladderID, sensorID, pdgCode, isPrimary, preStepPos, localpreStepPos, preStepMom, time);
         //Remember if the track came from the outside
         if (preStep.GetStepStatus() == fGeomBoundary) traversal.hasEntered();
       }
       //Add current step
       const G4ThreeVector postStepPos = postStep.GetPosition() / CLHEP::mm * Unit::mm;
       const double length = step->GetStepLength() / CLHEP::cm * Unit::cm;
       traversal.add(postStepPos, depEnergy, length);
  
       //check if we are leaving the volume
       bool isLeaving = (postStep.GetStepStatus() == fGeomBoundary);
       if (isLeaving) traversal.hasLeft();
  
       //if we are leaving or the track is stopped, finish it
       if (isLeaving || track.GetTrackStatus() >= fStopAndKill) {
         bool contained = traversal.isContained();
         bool saved = finishTrack();
         //we mark all particles as important which created a hit and entered or left the volume
         if (saved && !contained) {
           Simulation::TrackInfo::getInfo(track).setIgnore(false);
         }
         return saved;
       }
       // Track not finished, return false right now
       return false;
     }

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

beast/fangs/simulation/include/SensitiveDetector.h
beast/fangs/simulation/src/SensitiveDetector.cc

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Function Documentation

◆ ProcessHits()

◆ registerMCParticleRelation() [1/2]

◆ registerMCParticleRelation() [2/2]

◆ setActive()

◆ step()