Par03SensitiveDetector Class Reference

Sensitive detector. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par03SensitiveDetector:

Public Member Functions
	Par03SensitiveDetector (G4String aName)
	Par03SensitiveDetector (G4String aName, G4int aNumLayers, G4int aNumPhi, G4int aNumRho)
virtual	~Par03SensitiveDetector ()
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.
Par03Hit *	RetrieveAndSetupHit (G4TouchableHistory *aTouchable)
	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.
Public Member Functions inherited from G4VFastSimSensitiveDetector
G4bool	Hit (const G4FastHit *aHit, const G4FastTrack *aTrack, G4TouchableHandle *aTouchable)

Private Attributes
Par03HitsCollection *	fHitsCollection = nullptr
	Collection of hits.
G4int	fHitCollectionID = -1
	ID of collection of hits.
G4int	fCellNoZ = 10
	Number of readout cells along z axis.
G4int	fCellNoRho = 10
	Number of readout cells along radius of cylinder.
G4int	fCellNoPhi = 10
	Number of readout cells along azimuthal angle.

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 68 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

Par03SensitiveDetector() [1/2]

Par03SensitiveDetector::Par03SensitiveDetector

(

G4String

aName

)

Definition at line 35 of file Par03SensitiveDetector.cc.

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

Par03SensitiveDetector() [2/2]

Par03SensitiveDetector::Par03SensitiveDetector	(	G4String	aName,
		G4int	aNumLayers,
		G4int	aNumPhi,
		G4int	aNumRho
	)

Definition at line 42 of file Par03SensitiveDetector.cc.

  : G4VSensitiveDetector(aName)
  , fCellNoZ(aNumLayers)
  , fCellNoRho(aNumRho)
  , fCellNoPhi(aNumPhi)
{
  collectionName.insert("hits");
}

~Par03SensitiveDetector()

Par03SensitiveDetector::~Par03SensitiveDetector ( )

virtualdefault

Member Function Documentation

Initialize()

void Par03SensitiveDetector::Initialize ( G4HCofThisEvent *  HCE )

finalvirtual

Create hit collection.

Definition at line 58 of file Par03SensitiveDetector.cc.

{
  fHitsCollection =
    new Par03HitsCollection(SensitiveDetectorName, collectionName[0]);
  if(fHitCollectionID < 0)
  {
    fHitCollectionID =
      G4SDManager::GetSDMpointer()->GetCollectionID(fHitsCollection);
  }
  aHCE->AddHitsCollection(fHitCollectionID, fHitsCollection);
 
  // fill calorimeter hits with zero energy deposition
  for(G4int iphi = 0; iphi < fCellNoPhi; iphi++)
    for(G4int irho = 0; irho < fCellNoRho; irho++)
      for(G4int iz = 0; iz < fCellNoZ; iz++)
      {
        auto  hit = new Par03Hit();
        fHitsCollection->insert(hit);
      }
}

ProcessHits() [1/2]

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

finalvirtual

Process energy deposit from the full simulation.

Definition at line 81 of file Par03SensitiveDetector.cc.

{
  G4double edep = aStep->GetTotalEnergyDeposit();
  if(edep == 0.)
    return true;
 
  auto  aTouchable =
    (G4TouchableHistory*) (aStep->GetPreStepPoint()->GetTouchable());
 
  auto hit = RetrieveAndSetupHit(aTouchable);
 
  // Add energy deposit from G4Step
  hit->AddEdep(edep);
 
  // 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 Par03SensitiveDetector::ProcessHits ( const G4FastHit *  aHit, const G4FastTrack *  aTrack, G4TouchableHistory *  aROhist  )

finalvirtual

Process energy deposit from the fast simulation.

Implements G4VFastSimSensitiveDetector.

Definition at line 111 of file Par03SensitiveDetector.cc.

{
  G4double edep = aHit->GetEnergy();
  if(edep == 0.)
    return true;
 
  auto hit = RetrieveAndSetupHit(aTouchable);
 
  // Add energy deposit from G4FastHit
  hit->AddEdep(edep);
 
  // 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()

Par03Hit * Par03SensitiveDetector::RetrieveAndSetupHit

(

G4TouchableHistory *

aTouchable

)

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 141 of file Par03SensitiveDetector.cc.

{
  G4int rhoNo = aTouchable->GetCopyNumber(0);  // cell
  G4int phiNo = aTouchable->GetCopyNumber(1);  // segment
  G4int zNo   = aTouchable->GetCopyNumber(2);  // layer
 
  std::size_t hitID = fCellNoRho * fCellNoZ * phiNo + fCellNoZ * rhoNo + zNo;
 
  if(hitID >= fHitsCollection->entries())
  {
    G4Exception(
      "Par03SensitiveDetector::RetrieveAndSetupHit()", "InvalidSetup",
      FatalException,
      "Size of hit collection in Par03SensitiveDetector is smaller than the "
      "number of cells created in Par03DetectorConstruction!");
  }
  Par03Hit* hit = (*fHitsCollection)[hitID];
 
  if(hit->GetRhoId() < 0)
  {
    hit->SetRhoId(rhoNo);
    hit->SetPhiId(phiNo);
    hit->SetZid(zNo);
    hit->SetLogV(aTouchable->GetVolume(0)->GetLogicalVolume());
    G4AffineTransform transform = aTouchable->GetHistory()->GetTopTransform();
    hit->SetRot(transform.NetRotation());
    transform.Invert();
    hit->SetPos(transform.NetTranslation());
  }
  return hit;
}

Member Data Documentation

fHitsCollection

Par03HitsCollection* Par03SensitiveDetector::fHitsCollection = nullptr

private

Collection of hits.

Definition at line 74 of file Par03SensitiveDetector.hh.

fHitCollectionID

G4int Par03SensitiveDetector::fHitCollectionID = -1

private

ID of collection of hits.

Definition at line 76 of file Par03SensitiveDetector.hh.

fCellNoZ

G4int Par03SensitiveDetector::fCellNoZ = 10

private

Number of readout cells along z axis.

Definition at line 78 of file Par03SensitiveDetector.hh.

fCellNoRho

G4int Par03SensitiveDetector::fCellNoRho = 10

private

Number of readout cells along radius of cylinder.

Definition at line 80 of file Par03SensitiveDetector.hh.

fCellNoPhi

G4int Par03SensitiveDetector::fCellNoPhi = 10

private

Number of readout cells along azimuthal angle.

Definition at line 82 of file Par03SensitiveDetector.hh.

---

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

• .doxygen/Doxymodules_parameterisations.h
• extended/parameterisations/Par03/include/Par03SensitiveDetector.hh
• extended/parameterisations/Par03/src/Par03SensitiveDetector.cc