Ex03dCalorimeterSD Class Reference

The calorimeter sensitive detector. More...

#include <Ex03dCalorimeterSD.h>

Inheritance diagram for Ex03dCalorimeterSD:

Public Member Functions
	Ex03dCalorimeterSD (const char *name, Ex03DetectorConstruction *detector)
	Ex03dCalorimeterSD (const Ex03dCalorimeterSD &origin, Ex03DetectorConstruction *detector)
	Ex03dCalorimeterSD ()
virtual	~Ex03dCalorimeterSD ()
void	Initialize ()
Bool_t	ProcessHits ()
void	EndOfEvent ()
void	Register ()
virtual void	Print (Option_t *option="") const
void	PrintTotal () const
void	SetVerboseLevel (Int_t level)
Ex03CalorHit *	GetHit (Int_t i) const

Static Public Member Functions
static void	PrintTotal (std::vector< Ex03CalorHit > *collection)

Private Member Functions
void	ResetHits ()

Private Attributes
TVirtualMC *	fMC
	The VMC implementation.
Ex03DetectorConstruction *	fDetector
	Detector construction.
std::vector< Ex03CalorHit > *	fCalCollection
Int_t	fAbsorberVolId
	< The vector of particle (persistent)
Int_t	fGapVolId
	The gap volume Id.
Int_t	fVerboseLevel
	Verbosity level.

Detailed Description

The calorimeter sensitive detector.

A variant of the Ex03CalorimeterSD class updated for RNTuple output.

Date

07/07/2026

Author

Radoslaw Karabowicz; GSI

Definition at line 39 of file Ex03dCalorimeterSD.h.

Constructor & Destructor Documentation

Ex03dCalorimeterSD() [1/3]

Ex03dCalorimeterSD::Ex03dCalorimeterSD	(	const char *	name,
		Ex03DetectorConstruction *	detector )

Standard constructor. Create hits collection and an empty hit for each layer As the copy numbers may start from 0 or 1 (depending on geometry model, we create one more layer for this case.)

Parameters

name	The calorimeter hits collection name
detector	The detector construction

Definition at line 65 of file Ex03dCalorimeterSD.cxx.

  : TNamed(name, ""),
    fMC(0),
    fDetector(detector),
    fAbsorberVolId(0),
    fGapVolId(0),
    fVerboseLevel(1)
{
  /// Standard constructor.
  /// Create hits collection and an empty hit for each layer
  /// As the copy numbers may start from 0 or 1 (depending on
  /// geometry model, we create one more layer for this case.)
  /// \param name      The calorimeter hits collection name
  /// \param detector  The detector construction
 
  for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
    fCalCollection->push_back(Ex03CalorHit());
}

Ex03dCalorimeterSD() [2/3]

Ex03dCalorimeterSD::Ex03dCalorimeterSD	(	const Ex03dCalorimeterSD &	origin,
		Ex03DetectorConstruction *	detector )

Copy constructor (for clonig on worker thread in MT mode). Create hits collection and an empty hit for each layer As the copy numbers may start from 0 or 1 (depending on geometry model, we create one more layer for this case.)

Parameters

origin	The source object (on master).
detector	The detector construction

Definition at line 86 of file Ex03dCalorimeterSD.cxx.

  : TNamed(origin),
    fMC(0),
    fDetector(detector),
    fAbsorberVolId(origin.fAbsorberVolId),
    fGapVolId(origin.fGapVolId),
    fVerboseLevel(origin.fVerboseLevel)
{
  /// Copy constructor (for clonig on worker thread in MT mode).
  /// Create hits collection and an empty hit for each layer
  /// As the copy numbers may start from 0 or 1 (depending on
  /// geometry model, we create one more layer for this case.)
  /// \param origin    The source object (on master).
  /// \param detector  The detector construction
 
  for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
    fCalCollection->push_back(Ex03CalorHit());
}

Ex03dCalorimeterSD() [3/3]

Ex03dCalorimeterSD::Ex03dCalorimeterSD

(

)

Default constructor

Definition at line 107 of file Ex03dCalorimeterSD.cxx.

  : TNamed(), fDetector(0), fAbsorberVolId(0), fGapVolId(0), fVerboseLevel(1)
{
  /// Default constructor
}

~Ex03dCalorimeterSD()

Ex03dCalorimeterSD::~Ex03dCalorimeterSD ( )

virtual

Destructor

Definition at line 114 of file Ex03dCalorimeterSD.cxx.

{
  /// Destructor
 
  delete fCalCollection;
}

Member Function Documentation

PrintTotal() [1/2]

void Ex03dCalorimeterSD::PrintTotal ( std::vector< Ex03CalorHit > * collection )

static

Print the total values for all layers for the given collection

Definition at line 39 of file Ex03dCalorimeterSD.cxx.

{
  /// Print the total values for all layers for the given collection
 
  Double_t totEAbs = 0.;
  Double_t totLAbs = 0.;
  Double_t totEGap = 0.;
  Double_t totLGap = 0.;
 
  Int_t nofHits = collection->size();
  for (auto& hit : *collection) {
    totEAbs += hit.GetEdepAbs();
    totLAbs += hit.GetTrakAbs();
    totEGap += hit.GetEdepGap();
    totLGap += hit.GetTrakGap();
  }
 
  cout << "   Absorber: total energy (MeV): " << setw(7) << totEAbs * 1.0e03
       << "       total track length (cm):  " << setw(7) << totLAbs << endl
       << "   Gap:      total energy (MeV): " << setw(7) << totEGap * 1.0e03
       << "       total track length (cm):  " << setw(7) << totLGap << endl;
}

Initialize()

void Ex03dCalorimeterSD::Initialize

(

)

Definition at line 147 of file Ex03dCalorimeterSD.cxx.

{
  // Keep the pointer to TVirtualMC object as a data member
  // to avoid a possible performance penalty due to a frequent retrieval
  // from the thread-local storage
  fMC = gMC;
 
  fAbsorberVolId = fMC->VolId("ABSO");
  fGapVolId = fMC->VolId("GAPX");
}

ProcessHits()

Bool_t Ex03dCalorimeterSD::ProcessHits

(

)

Account energy deposit and track lengths for each layer in its hit.

Definition at line 159 of file Ex03dCalorimeterSD.cxx.

{
  /// Account energy deposit and track lengths for each layer in its hit.
 
  Int_t copyNo;
  Int_t id = fMC->CurrentVolID(copyNo);
 
  if (id != fAbsorberVolId && id != fGapVolId) return false;
 
  fMC->CurrentVolOffID(2, copyNo);
  // cout << "Got copyNo "<< copyNo << " " << fMC->CurrentVolPath() << endl;
 
  Double_t edep = fMC->Edep();
 
  Double_t step = 0.;
  if (fMC->TrackCharge() != 0.) step = fMC->TrackStep();
 
  if (!GetHit(copyNo)) {
    std::cerr << "No hit found for layer with copyNo = " << copyNo << endl;
    return false;
  }
 
  if (id == fAbsorberVolId) {
    GetHit(copyNo)->AddAbs(edep, step);
  }
 
  if (id == fGapVolId) {
    GetHit(copyNo)->AddGap(edep, step);
  }
 
  return true;
}

EndOfEvent()

void Ex03dCalorimeterSD::EndOfEvent

(

)

Print hits collection (if verbose) and reset hits afterwards.

Definition at line 193 of file Ex03dCalorimeterSD.cxx.

{
  /// Print hits collection (if verbose) and reset hits afterwards.
 
  if (fVerboseLevel > 1) Print();
 
  // Reset hits collection
  ResetHits();
}

Register()

void Ex03dCalorimeterSD::Register

(

)

Register the hits collection in Root manager.

Definition at line 204 of file Ex03dCalorimeterSD.cxx.

{
  /// Register the hits collection in Root manager.
  TMCRootManager::Instance()->Register("hits", fCalCollection);
}

Print()

void Ex03dCalorimeterSD::Print ( Option_t * option = "" ) const

virtual

Print the hits collection.

Definition at line 211 of file Ex03dCalorimeterSD.cxx.

{
  /// Print the hits collection.
 
  Int_t nofHits = fCalCollection->size();
 
  cout << "\n-------->Hits Collection: in this event: " << endl;
 
  for (Int_t i = 0; i < nofHits; i++) (*fCalCollection)[i].Print();
}

PrintTotal() [2/2]

void Ex03dCalorimeterSD::PrintTotal

(

)

const

Print the total values for all layers.

Definition at line 223 of file Ex03dCalorimeterSD.cxx.

{
  /// Print the total values for all layers.
 
  PrintTotal(fCalCollection);
}

SetVerboseLevel()

void Ex03dCalorimeterSD::SetVerboseLevel ( Int_t level )

inline

Set verbose level

Parameters

level

The new verbose level value

Definition at line 84 of file Ex03dCalorimeterSD.h.

{
  fVerboseLevel = level;
}

GetHit()

Ex03CalorHit * Ex03dCalorimeterSD::GetHit

(

Int_t

i

)

const

Returns

The hit for the specified layer.

Parameters

i	The layer number

Definition at line 126 of file Ex03dCalorimeterSD.cxx.

{
  /// \return   The hit for the specified layer.
  /// \param i  The layer number
 
  return &(*fCalCollection)[i];
}

ResetHits()

void Ex03dCalorimeterSD::ResetHits ( )

private

Reset all hits in the hits collection.

Definition at line 135 of file Ex03dCalorimeterSD.cxx.

{
  /// Reset all hits in the hits collection.
 
  for (Int_t i = 0; i < fCalCollection->size(); i++) GetHit(i)->Reset();
}

Member Data Documentation

fMC

TVirtualMC* Ex03dCalorimeterSD::fMC

private

The VMC implementation.

Definition at line 70 of file Ex03dCalorimeterSD.h.

fDetector

Ex03DetectorConstruction* Ex03dCalorimeterSD::fDetector

private

Detector construction.

Definition at line 71 of file Ex03dCalorimeterSD.h.

fCalCollection

std::vector<Ex03CalorHit>* Ex03dCalorimeterSD::fCalCollection

private

Initial value:

{
    new std::vector<Ex03CalorHit>
  }

Definition at line 72 of file Ex03dCalorimeterSD.h.

                                         {
    new std::vector<Ex03CalorHit>
  };                    //! < The vector of particle (persistent)

fAbsorberVolId

Int_t Ex03dCalorimeterSD::fAbsorberVolId

private

< The vector of particle (persistent)

The absorber volume Id

Definition at line 75 of file Ex03dCalorimeterSD.h.

fGapVolId

Int_t Ex03dCalorimeterSD::fGapVolId

private

The gap volume Id.

Definition at line 76 of file Ex03dCalorimeterSD.h.

fVerboseLevel

Int_t Ex03dCalorimeterSD::fVerboseLevel

private

Verbosity level.

Definition at line 77 of file Ex03dCalorimeterSD.h.

---

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

• examples/E03/E03d/include/Ex03dCalorimeterSD.h
• examples/E03/E03d/src/Ex03dCalorimeterSD.cxx