This is optional (temporary) class that provides information on track parameters on aerogel plane, until tracking is not prepared. More...

#include <SensitiveAero.h>

Inheritance diagram for SensitiveAero:

Public Member Functions
	SensitiveAero ()
	Constructor.

bool	step (G4Step aStep, G4TouchableHistory ) override
	Process each step and calculate variables defined in PXDSimHit.

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.

static void	registerMCParticleRelation (const std::string &name, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
	Register an relation involving MCParticles.

static void	registerMCParticleRelation (const RelationArray &relation, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
	Overload to make it easer to register MCParticle relations.

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.

Private Attributes
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

This is optional (temporary) class that provides information on track parameters on aerogel plane, until tracking is not prepared.

Definition at line 22 of file SensitiveAero.h.

Constructor & Destructor Documentation

◆ SensitiveAero()

SensitiveAero ( )

Constructor.

Definition at line 30 of file SensitiveAero.cc.

                                :
      Simulation::SensitiveDetectorBase("ARICH", Const::ARICH)
    {
 
      // registration of store arrays and relations
 
      StoreArray<MCParticle> mcParticles;
      StoreArray<ARICHAeroHit> aeroHits;
      aeroHits.registerInDataStore();
      mcParticles.registerRelationTo(aeroHits);
 
      // additional registration of MCParticle relation (required for correct relations)
 
      RelationArray relation(mcParticles, aeroHits);
      registerMCParticleRelation(relation, RelationArray::c_deleteElement);
 
    }

Member Function Documentation

◆ getMCParticleRelations()

static const std::map< std::string, RelationArray::EConsolidationAction > & getMCParticleRelations ( )

inlinestaticinherited

Return a list of all registered Relations with MCParticles.

Definition at line 42 of file SensitiveDetectorBase.h.

42{ return s_mcRelations; }

Belle2::Simulation::SensitiveDetectorBase::s_mcRelations

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.

Definition: SensitiveDetectorBase.h:87

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

    {
      if (!s_active) return false;
      bool result = step(aStep, aROhist);
      // Do not include hits from invalid detector (beast,teastbeam, etc.)
      if (result && (m_subdetector != Const::invalidDetector)) TrackInfo::getInfo(*aStep->GetTrack()).addSeenInDetector(m_subdetector);
      return result;
    }

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

67 { registerMCParticleRelation(relation.getName(), ignoreAction); }

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

    {
      std::pair<std::map<std::string, RelationArray::EConsolidationAction>::iterator, bool> insert = s_mcRelations.insert(std::make_pair(
            name, ignoreAction));
      //If the relation already exists and the ignoreAction is different we do have a problem
      if (!insert.second && insert.first->second != ignoreAction) {
        B2FATAL("MCParticle Relation " << name << " already registered with different ignore action.");
      }
    }

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

50{ s_active = activeStatus; }

◆ step()

bool step	(	G4Step *	aStep,
		G4TouchableHistory *
	)

overridevirtual

Process each step and calculate variables defined in PXDSimHit.

Parameters

aStep Current Geant4 step in the sensitive medium.

Returns: true if a hit was stored, o.w. false.

Implements SensitiveDetectorBase.

Definition at line 48 of file SensitiveAero.cc.

    {
      // Get track parameters
 
      G4Track* aTrack = aStep->GetTrack();
 
      G4StepPoint* PostPosition = aStep->GetPostStepPoint();
      G4ThreeVector worldPosition = PostPosition->GetPosition();
      G4ParticleDefinition* particle = aTrack->GetDefinition();
      G4double  PDGCharge = particle->GetPDGCharge();
      G4ThreeVector momentum = PostPosition->GetMomentum();
 
      // Save only tracks of charged particles
      if (PDGCharge == 0) return (true);
 
      // Track parameters are saved at the entrance in aerogel
 
      if ((PostPosition->GetStepStatus() == fGeomBoundary) && (momentum.z() > 0.0)) {
 
        /*       B2INFO ("SensAero: " << aTrack->GetDefinition()->GetParticleName()
           << " " << aTrack->GetTrackID()
           << " " << aTrack->GetParentID()
           << " " << G4BestUnit(worldPosition,"Length")
           << " " << G4BestUnit(aTrack->GetMomentum(), "Energy")
           << " " << G4BestUnit(aTrack->GetGlobalTime(), "Time")
           << " Edep is " << G4BestUnit(aStep->GetTotalEnergyDeposit(),"Energy"));
        */
 
        int trackID = aTrack->GetTrackID();
        int PDGEncoding = particle->GetPDGEncoding();
 
        ROOT::Math::XYZVector TPosition(worldPosition.x() * Unit::mm, worldPosition.y() * Unit::mm, worldPosition.z() * Unit::mm);
        ROOT::Math::XYZVector TMomentum(momentum.x() * Unit::MeV, momentum.y() * Unit::MeV, momentum.z() * Unit::MeV);
 
        // write the hit in datastore"
        StoreArray<ARICHAeroHit> aeroHits;
        ARICHAeroHit* aeroHit = aeroHits.appendNew(PDGEncoding, TPosition, TMomentum);
 
        // Create relation to MCParticle
        StoreArray<MCParticle> mcParticles;
        RelationArray rel(mcParticles, aeroHits);
        rel.add(trackID, aeroHit->getArrayIndex());
 
      }
 
      return true;
 
    }

Member Data Documentation

◆ m_subdetector

Const::EDetector m_subdetector

privateinherited

Subdetector the class belongs to.

Definition at line 91 of file SensitiveDetectorBase.h.

◆ s_active

bool s_active

staticprivateinherited

Static bool which indicates wether recording of hits is enabled.

Definition at line 89 of file SensitiveDetectorBase.h.

◆ s_mcRelations

map< string, RelationArray::EConsolidationAction > s_mcRelations

staticprivateinherited

Static set holding all relations which have to be updated at the end of the Event.

Definition at line 87 of file SensitiveDetectorBase.h.

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

arich/simulation/include/SensitiveAero.h
arich/simulation/src/SensitiveAero.cc

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ SensitiveAero()

Member Function Documentation

◆ getMCParticleRelations()

◆ ProcessHits()

◆ registerMCParticleRelation() [1/2]

◆ registerMCParticleRelation() [2/2]

◆ setActive()

◆ step()

Member Data Documentation

◆ m_subdetector

◆ s_active

◆ s_mcRelations