Ex03bCalorimeterSD Class Reference

The calorimeter sensitive detector. More...

#include <Ex03bCalorimeterSD.h>

Inheritance diagram for Ex03bCalorimeterSD:

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

Private Member Functions
void	ResetHits ()

Private Attributes
TVirtualMC *	fMC
	The VMC implementation.
Ex03DetectorConstruction *	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.
Int_t	fPrintModulo
	The event modulus number to be printed.
Ex03CalorHit	fReflectedHit
	Reflected contribution to the layer hits.

Detailed Description

The calorimeter sensitive detector.

A variant of the Ex03CalorimeterSD: a sensitive detector class derived fron the new TVirtualMCSensitiveDetector interface

Author

I. Hrivnacova; IPN, Orsay

Definition at line 39 of file Ex03bCalorimeterSD.h.

Constructor & Destructor Documentation

Ex03bCalorimeterSD() [1/3]

Ex03bCalorimeterSD::Ex03bCalorimeterSD	(	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 37 of file Ex03bCalorimeterSD.cxx.

  : TVirtualMCSensitiveDetector(name, ""),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    fAbsorberVolId(0),
    fGapVolId(0),
    fVerboseLevel(1),
    fPrintModulo(1),
    fReflectedHit()
{
  /// 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();
}

Ex03bCalorimeterSD() [2/3]

Ex03bCalorimeterSD::Ex03bCalorimeterSD	(	const Ex03bCalorimeterSD &	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 62 of file Ex03bCalorimeterSD.cxx.

  : TVirtualMCSensitiveDetector(origin),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    fAbsorberVolId(origin.fAbsorberVolId),
    fGapVolId(origin.fGapVolId),
    fVerboseLevel(origin.fVerboseLevel),
    fPrintModulo(origin.fPrintModulo),
    fReflectedHit()
{
  /// 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();
}

Ex03bCalorimeterSD() [3/3]

Ex03bCalorimeterSD::Ex03bCalorimeterSD

(

)

Default constructor

Definition at line 87 of file Ex03bCalorimeterSD.cxx.

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

~Ex03bCalorimeterSD()

Ex03bCalorimeterSD::~Ex03bCalorimeterSD ( )

virtual

Destructor

Definition at line 100 of file Ex03bCalorimeterSD.cxx.

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

Member Function Documentation

Initialize()

void Ex03bCalorimeterSD::Initialize

(

)

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

Definition at line 135 of file Ex03bCalorimeterSD.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
  fMC = gMC;
 
  fAbsorberVolId = fMC->VolId("ABSO");
  fGapVolId = fMC->VolId("GAPX");
}

ProcessHits()

void Ex03bCalorimeterSD::ProcessHits ( )

virtual

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

Definition at line 152 of file Ex03bCalorimeterSD.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;
 
  const Bool_t reflected = fDetector->GetUseReflection() && copyNo == 2;
  fMC->CurrentVolOffID(2, copyNo);
 
  Double_t edep = fMC->Edep();
 
  Double_t step = 0.;
  if (fMC->TrackCharge() != 0.) step = fMC->TrackStep();
 
  if (!GetHit(copyNo)) {
    return;
  }
 
  if (id == fAbsorberVolId) {
    GetHit(copyNo)->AddAbs(edep, step);
    if (reflected) fReflectedHit.AddAbs(edep, step);
  }
 
  if (id == fGapVolId) {
    GetHit(copyNo)->AddGap(edep, step);
    if (reflected) fReflectedHit.AddGap(edep, step);
  }
}

EndOfEvent()

void Ex03bCalorimeterSD::EndOfEvent

(

)

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

Definition at line 185 of file Ex03bCalorimeterSD.cxx.

{
  /// Print hits collection (if verbose) and reset hits afterwards.
 
  if (gMC->CurrentEvent() % fPrintModulo == 0) {
    PrintTotal();
  }
 
  if (fVerboseLevel > 1) Print();
 
  if (fDetector->GetUseReflection()) {
    cout << "   Reflected placements (included in layer totals):" << endl;
    fReflectedHit.Print();
  }
  if (fDetector->GetRequireReflectedHits()) {
    if (fReflectedHit.GetEdepAbs() <= 0. || fReflectedHit.GetTrakAbs() <= 0. ||
        fReflectedHit.GetEdepGap() <= 0. || fReflectedHit.GetTrakGap() <= 0.) {
      Fatal("EndOfEvent", "Missing hits in reflected absorber or gap");
    }
    cout << "   Reflected hits test passed" << endl;
  }
  fReflectedHit.Reset();
 
  // Reset hits collection
  ResetHits();
}

Register()

void Ex03bCalorimeterSD::Register

(

)

Register the hits collection in Root manager.

Definition at line 213 of file Ex03bCalorimeterSD.cxx.

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

Print()

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

virtual

Print the hits collection.

Definition at line 221 of file Ex03bCalorimeterSD.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 Ex03bCalorimeterSD::PrintTotal

(

)

const

Print the total values for all layers.

Definition at line 233 of file Ex03bCalorimeterSD.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 Ex03bCalorimeterSD::SetVerboseLevel ( Int_t level )

inline

Set verbose level

Parameters

level

The new verbose level value

Definition at line 83 of file Ex03bCalorimeterSD.h.

{
  fVerboseLevel = level;
}

SetPrintModulo()

void Ex03bCalorimeterSD::SetPrintModulo ( Int_t value )

inline

Set the event modulus number to be printed

Parameters

value

The new event modulus number value

Definition at line 90 of file Ex03bCalorimeterSD.h.

{
  fPrintModulo = value;
}

GetHit()

Ex03CalorHit * Ex03bCalorimeterSD::GetHit

(

Int_t

i

)

const

Returns

The hit for the specified layer.

Parameters

i	The layer number

Definition at line 113 of file Ex03bCalorimeterSD.cxx.

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

ResetHits()

void Ex03bCalorimeterSD::ResetHits ( )

private

Reset all hits in the hits collection.

Definition at line 122 of file Ex03bCalorimeterSD.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* Ex03bCalorimeterSD::fMC

private

The VMC implementation.

Definition at line 68 of file Ex03bCalorimeterSD.h.

fDetector

Ex03DetectorConstruction* Ex03bCalorimeterSD::fDetector

private

Detector construction.

Definition at line 69 of file Ex03bCalorimeterSD.h.

fCalCollection

TClonesArray* Ex03bCalorimeterSD::fCalCollection

private

Hits collection.

Definition at line 70 of file Ex03bCalorimeterSD.h.

fAbsorberVolId

Int_t Ex03bCalorimeterSD::fAbsorberVolId

private

The absorber volume Id.

Definition at line 71 of file Ex03bCalorimeterSD.h.

fGapVolId

Int_t Ex03bCalorimeterSD::fGapVolId

private

The gap volume Id.

Definition at line 72 of file Ex03bCalorimeterSD.h.

fVerboseLevel

Int_t Ex03bCalorimeterSD::fVerboseLevel

private

Verbosity level.

Definition at line 73 of file Ex03bCalorimeterSD.h.

fPrintModulo

Int_t Ex03bCalorimeterSD::fPrintModulo

private

The event modulus number to be printed.

Definition at line 74 of file Ex03bCalorimeterSD.h.

fReflectedHit

Ex03CalorHit Ex03bCalorimeterSD::fReflectedHit

private

Reflected contribution to the layer hits.

Definition at line 76 of file Ex03bCalorimeterSD.h.

---

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

• examples/E03/E03b/include/Ex03bCalorimeterSD.h
• examples/E03/E03b/src/Ex03bCalorimeterSD.cxx