Ex03CalorimeterSD Class Reference

The calorimeter sensitive detector. More...

#include <Ex03CalorimeterSD.h>

Inheritance diagram for Ex03CalorimeterSD:

Public Member Functions
	Ex03CalorimeterSD (const char *name, Ex03DetectorConstruction *detector)
	Ex03CalorimeterSD (const Ex03CalorimeterSD &origin, Ex03DetectorConstruction *detector)
	Ex03CalorimeterSD ()
virtual	~Ex03CalorimeterSD ()
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.
Ex03DetectorConstruction *	fDetector
	Detector construction.
TClonesArray *	fCalCollection
	Hits collection.
Int_t	fAbsorberVolId
	The absorber volume Id.
Int_t	fGapVolId
	The gap volume Id.
Bool_t	fAssemblyHierarchyPrinted
	Whether assembly levels were printed.
Int_t	fVerboseLevel
	Verbosity level.
Ex03CalorHit	fReflectedHit
	Reflected contribution to the layer hits.

Detailed Description

The calorimeter sensitive detector.

Date

06/03/2003

Author

I. Hrivnacova; IPN, Orsay

Definition at line 36 of file Ex03CalorimeterSD.h.

Constructor & Destructor Documentation

Ex03CalorimeterSD() [1/3]

Ex03CalorimeterSD::Ex03CalorimeterSD	(	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 Ex03CalorimeterSD.cxx.

  : TNamed(name, ""),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    fAbsorberVolId(0),
    fGapVolId(0),
    fAssemblyHierarchyPrinted(kFALSE),
    fVerboseLevel(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();
}

Ex03CalorimeterSD() [2/3]

Ex03CalorimeterSD::Ex03CalorimeterSD	(	const Ex03CalorimeterSD &	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 Ex03CalorimeterSD.cxx.

  : TNamed(origin),
    fMC(0),
    fDetector(detector),
    fCalCollection(0),
    fAbsorberVolId(origin.fAbsorberVolId),
    fGapVolId(origin.fGapVolId),
    fAssemblyHierarchyPrinted(kFALSE),
    fVerboseLevel(origin.fVerboseLevel),
    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();
}

Ex03CalorimeterSD() [3/3]

Ex03CalorimeterSD::Ex03CalorimeterSD

(

)

Default constructor

Definition at line 87 of file Ex03CalorimeterSD.cxx.

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

~Ex03CalorimeterSD()

Ex03CalorimeterSD::~Ex03CalorimeterSD ( )

virtual

Destructor

Definition at line 101 of file Ex03CalorimeterSD.cxx.

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

Member Function Documentation

Initialize()

void Ex03CalorimeterSD::Initialize

(

)

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

Definition at line 136 of file Ex03CalorimeterSD.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()

Bool_t Ex03CalorimeterSD::ProcessHits

(

)

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

Definition at line 153 of file Ex03CalorimeterSD.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;
 
  if (fDetector->GetUseAssemblies() && !fAssemblyHierarchyPrinted) {
    fAssemblyHierarchyPrinted = kTRUE;
    cout << "Assembly volume hierarchy:" << endl;
    for (Int_t off = 0; off <= 6; ++off) {
      Int_t offCopyNo = -1;
      TString offName = fMC->CurrentVolOffName(off);
      Int_t offId = fMC->CurrentVolOffID(off, offCopyNo);
      cout << "  off=" << off << " name=" << offName << " id=" << offId
           << " copyNo=" << offCopyNo << endl;
    }
  }
 
  const Bool_t reflected = fDetector->GetUseReflection() && copyNo == 2;
  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 (reflected) fReflectedHit.AddAbs(edep, step);
  }
 
  if (id == fGapVolId) {
    GetHit(copyNo)->AddGap(edep, step);
    if (reflected) fReflectedHit.AddGap(edep, step);
  }
 
  return true;
}

EndOfEvent()

void Ex03CalorimeterSD::EndOfEvent

(

)

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

Definition at line 202 of file Ex03CalorimeterSD.cxx.

{
  /// Print hits collection (if verbose) and reset hits afterwards.
 
  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 Ex03CalorimeterSD::Register

(

)

Register the hits collection in Root manager.

Definition at line 226 of file Ex03CalorimeterSD.cxx.

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

Print()

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

virtual

Print the hits collection.

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

(

)

const

Print the total values for all layers.

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

inline

Set verbose level

Parameters

level

The new verbose level value

Definition at line 79 of file Ex03CalorimeterSD.h.

{
  fVerboseLevel = level;
}

GetHit()

Ex03CalorHit * Ex03CalorimeterSD::GetHit

(

Int_t

i

)

const

Returns

The hit for the specified layer.

Parameters

i	The layer number

Definition at line 114 of file Ex03CalorimeterSD.cxx.

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

ResetHits()

void Ex03CalorimeterSD::ResetHits ( )

private

Reset all hits in the hits collection.

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

private

The VMC implementation.

Definition at line 64 of file Ex03CalorimeterSD.h.

fDetector

Ex03DetectorConstruction* Ex03CalorimeterSD::fDetector

private

Detector construction.

Definition at line 65 of file Ex03CalorimeterSD.h.

fCalCollection

TClonesArray* Ex03CalorimeterSD::fCalCollection

private

Hits collection.

Definition at line 66 of file Ex03CalorimeterSD.h.

fAbsorberVolId

Int_t Ex03CalorimeterSD::fAbsorberVolId

private

The absorber volume Id.

Definition at line 67 of file Ex03CalorimeterSD.h.

fGapVolId

Int_t Ex03CalorimeterSD::fGapVolId

private

The gap volume Id.

Definition at line 68 of file Ex03CalorimeterSD.h.

fAssemblyHierarchyPrinted

Bool_t Ex03CalorimeterSD::fAssemblyHierarchyPrinted

private

Whether assembly levels were printed.

Definition at line 69 of file Ex03CalorimeterSD.h.

fVerboseLevel

Int_t Ex03CalorimeterSD::fVerboseLevel

private

Verbosity level.

Definition at line 70 of file Ex03CalorimeterSD.h.

fReflectedHit

Ex03CalorHit Ex03CalorimeterSD::fReflectedHit

private

Reflected contribution to the layer hits.

Definition at line 72 of file Ex03CalorimeterSD.h.

---

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

• examples/E03/E03a/include/Ex03CalorimeterSD.h
• examples/E03/E03a/src/Ex03CalorimeterSD.cxx