Par04SensitiveDetector Class Reference

Sensitive detector. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par04SensitiveDetector:

Public Member Functions
	Par04SensitiveDetector (G4String aName)
	Par04SensitiveDetector (G4String aName, G4ThreeVector aNbOfCells, G4ThreeVector aNSizeOfCells)
virtual	~Par04SensitiveDetector ()
virtual void	Initialize (G4HCofThisEvent *HCE) final
	Create hit collection.
virtual G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *aROhist) final
	Process energy deposit from the full simulation.
virtual G4bool	ProcessHits (const G4FastHit *aHit, const G4FastTrack *aTrack, G4TouchableHistory *aROhist) final
	Process energy deposit from the fast simulation.
Par04Hit *	RetrieveAndSetupHit (G4ThreeVector aPosition)
	Process energy deposit - common part for full and fast simulation It is invoked from ProcessHits() methods, and sets basic hit properties (position, etc.), common for hit from fast and full simulation.
virtual void	EndOfEvent (G4HCofThisEvent *aHC) final
	Rewrite hits map to a vector.
Public Member Functions inherited from G4VFastSimSensitiveDetector
G4bool	Hit (const G4FastHit *aHit, const G4FastTrack *aTrack, G4TouchableHandle *aTouchable)

Private Attributes
Par04HitsCollection *	fHitsCollection = nullptr
	Collection of hits.
std::unordered_map< G4int, std::unique_ptr< Par04Hit > >	fHitsMap
	Map of hits to be used in runtime.
G4int	fHitCollectionID = -1
	ID of collection of hits.
G4ThreeVector	fMeshNbOfCells = { 10, 10, 10 }
	Number of mesh readout cells in cylindrical coordinates.
G4ThreeVector	fMeshSizeOfCells = { 1 * m, 2 * CLHEP::pi / 10., 1 * m }
	Size of mesh readout cells in cylindrical coordinates.
G4ThreeVector	fEntrancePosition = { -1, -1, -1 }
	Retrieved once per event: position of entering particle.
G4ThreeVector	fEntranceDirection = { -1, -1, -1 }
	Retrieved once per event: direction of entering particle.

Detailed Description

Sensitive detector.

Describes how to store the energy deposited within the detector. It derives from two classes: G4VSensitiveDetector and G4VFastSimSensitiveDetector. Addition of G4VFastSimSensitiveDetector is necessary in order to handle the energy deposits from the fast simulation.

Two ProcessHits() methods are introduced to handle energy deposited from full (detailed) simulation, and from fast simulation. The common part is handled by RetrieveAdnSetupHit() method.

Definition at line 92 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

Par04SensitiveDetector() [1/2]

Par04SensitiveDetector::Par04SensitiveDetector

(

G4String

aName

)

Definition at line 52 of file Par04SensitiveDetector.cc.

  : G4VSensitiveDetector(aName)
{
  collectionName.insert("hits");
}

Par04SensitiveDetector() [2/2]

Par04SensitiveDetector::Par04SensitiveDetector	(	G4String	aName,
		G4ThreeVector	aNbOfCells,
		G4ThreeVector	aNSizeOfCells
	)

Definition at line 59 of file Par04SensitiveDetector.cc.

  : G4VSensitiveDetector(aName)
  , fMeshNbOfCells(aNb)
  , fMeshSizeOfCells(aSize)
{
  collectionName.insert("hits");
}

~Par04SensitiveDetector()

Par04SensitiveDetector::~Par04SensitiveDetector ( )

virtualdefault

Member Function Documentation

Initialize()

void Par04SensitiveDetector::Initialize ( G4HCofThisEvent *  HCE )

finalvirtual

Create hit collection.

Definition at line 74 of file Par04SensitiveDetector.cc.

{
  fHitsCollection = new Par04HitsCollection(SensitiveDetectorName, collectionName[0]);
  if(fHitCollectionID < 0)
  {
    fHitCollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(fHitsCollection);
  }
  aHCE->AddHitsCollection(fHitCollectionID, fHitsCollection);
 
  // reset entrance position
  fEntrancePosition.set(-1, -1, -1);
  fEntranceDirection.set(-1, -1, -1);
}

ProcessHits() [1/2]

G4bool Par04SensitiveDetector::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  aROhist  )

finalvirtual

Process energy deposit from the full simulation.

Definition at line 90 of file Par04SensitiveDetector.cc.

{
  G4double edep = aStep->GetTotalEnergyDeposit();
  if(edep == 0.)
    return true;
 
  auto hit = RetrieveAndSetupHit(aStep->GetPostStepPoint()->GetPosition());
  if(hit == nullptr)
    return true;
 
  // Add energy deposit from G4Step
  hit->AddEdep(edep);
  // Increment the counter
  hit->AddNdep(1);
 
  // Fill time information from G4Step
  // If it's already filled, choose hit with earliest global time
  if(hit->GetTime() == -1 || hit->GetTime() > aStep->GetTrack()->GetGlobalTime())
    hit->SetTime(aStep->GetTrack()->GetGlobalTime());
 
  // Set hit type to full simulation (only if hit is not already marked as fast
  // sim)
  if(hit->GetType() != 1)
    hit->SetType(0);
 
  return true;
}

ProcessHits() [2/2]

G4bool Par04SensitiveDetector::ProcessHits ( const G4FastHit *  aHit, const G4FastTrack *  aTrack, G4TouchableHistory *  aROhist  )

finalvirtual

Process energy deposit from the fast simulation.

Implements G4VFastSimSensitiveDetector.

Definition at line 120 of file Par04SensitiveDetector.cc.

{
  G4double edep = aHit->GetEnergy();
  if(edep == 0.)
    return true;
 
  auto hit = RetrieveAndSetupHit(aHit->GetPosition());
  if(hit == nullptr)
    return true;
 
  // Add energy deposit from G4FastHit
  hit->AddEdep(edep);
  // Increment the counter
  hit->AddNdep(1);
 
  // Fill time information from G4FastTrack
  // If it's already filled, choose hit with earliest global time
  if(hit->GetTime() == -1 || hit->GetTime() > aTrack->GetPrimaryTrack()->GetGlobalTime())
  {
    hit->SetTime(aTrack->GetPrimaryTrack()->GetGlobalTime());
  }
 
  // Set hit type to fast simulation (even if hit was already marked as full
  // sim, overwrite it)
  hit->SetType(1);
 
  return true;
}

RetrieveAndSetupHit()

Par04Hit * Par04SensitiveDetector::RetrieveAndSetupHit

(

G4ThreeVector

aPosition

)

Process energy deposit - common part for full and fast simulation It is invoked from ProcessHits() methods, and sets basic hit properties (position, etc.), common for hit from fast and full simulation.

Definition at line 152 of file Par04SensitiveDetector.cc.

{
  if(fEntrancePosition.x() == -1)
  {
    auto  info = dynamic_cast<Par04EventInformation*>(
      G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetUserInformation());
    if(info == nullptr)
      return nullptr;
    fEntrancePosition  = info->GetPosition();
    fEntranceDirection = info->GetDirection();
  }
 
  auto delta = aGlobalPosition - fEntrancePosition;
 
  // Calculate rotation matrix along the particle momentum direction
  // It will rotate the shower axes to match the incoming particle direction
  G4RotationMatrix rotMatrix = G4RotationMatrix();
  double particleTheta       = fEntranceDirection.theta();
  double particlePhi         = fEntranceDirection.phi();
  rotMatrix.rotateZ(-particlePhi);
  rotMatrix.rotateY(-particleTheta);
  G4RotationMatrix rotMatrixInv = CLHEP::inverseOf(rotMatrix);
 
  delta = rotMatrix * delta;
 
  G4int rhoNo = std::floor(delta.perp() / fMeshSizeOfCells.x());
  G4int phiNo = std::floor((CLHEP::pi + delta.phi()) / fMeshSizeOfCells.y());
  G4int zNo   = std::floor(delta.z() / fMeshSizeOfCells.z());
 
  std::size_t hitID =
    fMeshNbOfCells.x() * fMeshNbOfCells.z() * phiNo + fMeshNbOfCells.z() * rhoNo + zNo;
  
  if(zNo >= fMeshNbOfCells.z() || rhoNo >= fMeshNbOfCells.x() || zNo < 0)
  {
    return nullptr;
  }
 
  auto hit = fHitsMap[hitID].get();
  if (hit==nullptr) {
    fHitsMap[hitID] = std::unique_ptr<Par04Hit>(new Par04Hit());
    hit = fHitsMap[hitID].get();
    hit->SetPhiId(phiNo);
    hit->SetRhoId(rhoNo);
    hit->SetZid(zNo);
    hit->SetRot(rotMatrixInv);
    hit->SetPos(fEntrancePosition +
                rotMatrixInv * G4ThreeVector(0, 0, (zNo + 0.5) * fMeshSizeOfCells.z()));
  
  }
  return hit;
}

EndOfEvent()

void Par04SensitiveDetector::EndOfEvent ( G4HCofThisEvent *  aHC )

finalvirtual

Rewrite hits map to a vector.

Definition at line 206 of file Par04SensitiveDetector.cc.

{
  for(const auto& hits: fHitsMap){
    fHitsCollection->insert(new Par04Hit(*hits.second.get()));
  }
  fHitsMap.clear();
}

Member Data Documentation

fHitsCollection

Par04HitsCollection* Par04SensitiveDetector::fHitsCollection = nullptr

private

Collection of hits.

Definition at line 82 of file Par04SensitiveDetector.hh.

fHitsMap

std::unordered_map<G4int, std::unique_ptr<Par04Hit> > Par04SensitiveDetector::fHitsMap

private

Map of hits to be used in runtime.

Definition at line 84 of file Par04SensitiveDetector.hh.

fHitCollectionID

G4int Par04SensitiveDetector::fHitCollectionID = -1

private

ID of collection of hits.

Definition at line 86 of file Par04SensitiveDetector.hh.

fMeshNbOfCells

G4ThreeVector Par04SensitiveDetector::fMeshNbOfCells = { 10, 10, 10 }

private

Number of mesh readout cells in cylindrical coordinates.

Definition at line 88 of file Par04SensitiveDetector.hh.

{ 10, 10, 10 };

fMeshSizeOfCells

G4ThreeVector Par04SensitiveDetector::fMeshSizeOfCells = { 1 * m, 2 * CLHEP::pi / 10., 1 * m }

private

Size of mesh readout cells in cylindrical coordinates.

Definition at line 90 of file Par04SensitiveDetector.hh.

{ 1 * m, 2 * CLHEP::pi / 10., 1 * m };

fEntrancePosition

G4ThreeVector Par04SensitiveDetector::fEntrancePosition = { -1, -1, -1 }

private

Retrieved once per event: position of entering particle.

Definition at line 92 of file Par04SensitiveDetector.hh.

{ -1, -1, -1 };

fEntranceDirection

G4ThreeVector Par04SensitiveDetector::fEntranceDirection = { -1, -1, -1 }

private

Retrieved once per event: direction of entering particle.

Definition at line 94 of file Par04SensitiveDetector.hh.

{ -1, -1, -1 };

---

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

• .doxygen/Doxymodules_parameterisations.h
• extended/parameterisations/Par04/include/Par04SensitiveDetector.hh
• extended/parameterisations/Par04/src/Par04SensitiveDetector.cc