Ex03cCalorimeterSD Class Reference

The calorimeter sensitive detector. More...

#include <Ex03cCalorimeterSD.h>

Inheritance diagram for Ex03cCalorimeterSD:

Public Member Functions
	Ex03cCalorimeterSD (const char *name, Ex03cDetectorConstruction *detector)
	Ex03cCalorimeterSD (const Ex03cCalorimeterSD &origin, Ex03cDetectorConstruction *detector)
	Ex03cCalorimeterSD ()
virtual	~Ex03cCalorimeterSD ()
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

Private Member Functions
void	ResetHits ()

Private Attributes
TVirtualMC *	fMC
	The VMC implementation.
Ex03cDetectorConstruction *	fDetector
	Detector construction.
TClonesArray *	fCalCollection
	Hits collection.
Int_t	fAbsorberVolId
	The absorber volume Id.
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 multiple engine runs.

Date

21/08/2019

Author

Benedikt Volkel, CERN

Definition at line 38 of file Ex03cCalorimeterSD.h.

Constructor & Destructor Documentation

Ex03cCalorimeterSD() [1/3]

Ex03cCalorimeterSD::Ex03cCalorimeterSD	(	const char *	name,
		Ex03cDetectorConstruction *	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 38 of file Ex03cCalorimeterSD.cxx.

  : TNamed(name, ""),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    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
 
  fCalCollection = new TClonesArray("Ex03CalorHit", 500);
  for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
    new ((*fCalCollection)[i]) Ex03CalorHit();
}

Ex03cCalorimeterSD() [2/3]

Ex03cCalorimeterSD::Ex03cCalorimeterSD	(	const Ex03cCalorimeterSD &	origin,
		Ex03cDetectorConstruction *	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 61 of file Ex03cCalorimeterSD.cxx.

  : TNamed(origin),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    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
 
  fCalCollection = new TClonesArray("Ex03CalorHit", 500);
  for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
    new ((*fCalCollection)[i]) Ex03CalorHit();
}

Ex03cCalorimeterSD() [3/3]

Ex03cCalorimeterSD::Ex03cCalorimeterSD

(

)

Default constructor

Definition at line 84 of file Ex03cCalorimeterSD.cxx.

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

~Ex03cCalorimeterSD()

Ex03cCalorimeterSD::~Ex03cCalorimeterSD ( )

virtual

Destructor

Definition at line 96 of file Ex03cCalorimeterSD.cxx.

{
  /// Destructor
 
  if (fCalCollection) fCalCollection->Delete();
  delete fCalCollection;
}

Member Function Documentation

Initialize()

void Ex03cCalorimeterSD::Initialize

(

)

Register hits collection in the Root manager; set sensitive volumes.

Definition at line 131 of file Ex03cCalorimeterSD.cxx.

{
  /// Register hits collection in the Root manager;
  /// set sensitive volumes.
 
  if (TMCRootManager::Instance()) Register();
 
  // 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
  if (TMCManager::Instance()) {
    TMCManager::Instance()->ConnectEnginePointer(fMC);
  }
  else {
    fMC = gMC;
  }
 
  fAbsorberVolId = fMC->VolId("ABSO");
  fGapVolId = fMC->VolId("GAPX");
}

ProcessHits()

Bool_t Ex03cCalorimeterSD::ProcessHits

(

)

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

Definition at line 153 of file Ex03cCalorimeterSD.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 Ex03cCalorimeterSD::EndOfEvent

(

)

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

Definition at line 187 of file Ex03cCalorimeterSD.cxx.

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

Register()

void Ex03cCalorimeterSD::Register

(

)

Register the hits collection in Root manager.

Definition at line 198 of file Ex03cCalorimeterSD.cxx.

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

Print()

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

virtual

Print the hits collection.

Definition at line 206 of file Ex03cCalorimeterSD.cxx.

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

PrintTotal()

void Ex03cCalorimeterSD::PrintTotal

(

)

const

Print the total values for all layers.

Definition at line 218 of file Ex03cCalorimeterSD.cxx.

{
  /// Print the total values for all layers.
 
  Double_t totEAbs = 0.;
  Double_t totLAbs = 0.;
  Double_t totEGap = 0.;
  Double_t totLGap = 0.;
 
  Int_t nofHits = fCalCollection->GetEntriesFast();
  for (Int_t i = 0; i < nofHits; i++) {
    totEAbs += GetHit(i)->GetEdepAbs();
    totLAbs += GetHit(i)->GetTrakAbs();
    totEGap += GetHit(i)->GetEdepGap();
    totLGap += GetHit(i)->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;
}

SetVerboseLevel()

void Ex03cCalorimeterSD::SetVerboseLevel ( Int_t level )

inline

Set verbose level

Parameters

level

The new verbose level value

Definition at line 78 of file Ex03cCalorimeterSD.h.

{
  fVerboseLevel = level;
}

GetHit()

Ex03CalorHit * Ex03cCalorimeterSD::GetHit

(

Int_t

i

)

const

Returns

The hit for the specified layer.

Parameters

i	The layer number

Definition at line 109 of file Ex03cCalorimeterSD.cxx.

{
  /// \return   The hit for the specified layer.
  /// \param i  The layer number
 
  return (Ex03CalorHit*)fCalCollection->At(i);
}

ResetHits()

void Ex03cCalorimeterSD::ResetHits ( )

private

Reset all hits in the hits collection.

Definition at line 118 of file Ex03cCalorimeterSD.cxx.

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

Member Data Documentation

fMC

TVirtualMC* Ex03cCalorimeterSD::fMC

private

The VMC implementation.

Definition at line 66 of file Ex03cCalorimeterSD.h.

fDetector

Ex03cDetectorConstruction* Ex03cCalorimeterSD::fDetector

private

Detector construction.

Definition at line 67 of file Ex03cCalorimeterSD.h.

fCalCollection

TClonesArray* Ex03cCalorimeterSD::fCalCollection

private

Hits collection.

Definition at line 68 of file Ex03cCalorimeterSD.h.

fAbsorberVolId

Int_t Ex03cCalorimeterSD::fAbsorberVolId

private

The absorber volume Id.

Definition at line 69 of file Ex03cCalorimeterSD.h.

fGapVolId

Int_t Ex03cCalorimeterSD::fGapVolId

private

The gap volume Id.

Definition at line 70 of file Ex03cCalorimeterSD.h.

fVerboseLevel

Int_t Ex03cCalorimeterSD::fVerboseLevel

private

Verbosity level.

Definition at line 71 of file Ex03cCalorimeterSD.h.

---

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

• examples/E03/E03c/include/Ex03cCalorimeterSD.h
• examples/E03/E03c/src/Ex03cCalorimeterSD.cxx