VMC Examples Version 6.6
Loading...
Searching...
No Matches
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)
 
Ex03CalorHitGetHit (Int_t i) const
 

Private Member Functions

void ResetHits ()
 

Private Attributes

TVirtualMC * fMC
 The VMC implementation.
 
Ex03DetectorConstructionfDetector
 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.

Date
06/03/2003
Author
I. Hrivnacova; IPN, Orsay

Definition at line 35 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
nameThe calorimeter hits collection name
detectorThe detector construction

Definition at line 37 of file Ex03CalorimeterSD.cxx.

39 : TNamed(name, ""),
40 fMC(0),
41 fDetector(detector),
44 fGapVolId(0),
46{
47 /// Standard constructor.
48 /// Create hits collection and an empty hit for each layer
49 /// As the copy numbers may start from 0 or 1 (depending on
50 /// geometry model, we create one more layer for this case.)
51 /// \param name The calorimeter hits collection name
52 /// \param detector The detector construction
53
54 fCalCollection = new TClonesArray("Ex03CalorHit", 500);
55 for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
56 new ((*fCalCollection)[i]) Ex03CalorHit();
57}
The calorimeter hit.
Ex03DetectorConstruction * fDetector
Detector construction.
TClonesArray * fCalCollection
Hits collection.
Int_t fVerboseLevel
Verbosity level.
Int_t fGapVolId
The gap volume Id.
Int_t fAbsorberVolId
The absorber volume Id.
TVirtualMC * fMC
The VMC implementation.

◆ 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
originThe source object (on master).
detectorThe detector construction

Definition at line 60 of file Ex03CalorimeterSD.cxx.

62 : TNamed(origin),
63 fMC(0),
64 fDetector(detector),
67 fGapVolId(origin.fGapVolId),
69{
70 /// Copy constructor (for clonig on worker thread in MT mode).
71 /// Create hits collection and an empty hit for each layer
72 /// As the copy numbers may start from 0 or 1 (depending on
73 /// geometry model, we create one more layer for this case.)
74 /// \param origin The source object (on master).
75 /// \param detector The detector construction
76
77 fCalCollection = new TClonesArray("Ex03CalorHit", 500);
78 for (Int_t i = 0; i < fDetector->GetNbOfLayers() + 1; i++)
79 new ((*fCalCollection)[i]) Ex03CalorHit();
80}

◆ Ex03CalorimeterSD() [3/3]

Ex03CalorimeterSD::Ex03CalorimeterSD ( )

Default constructor

Definition at line 83 of file Ex03CalorimeterSD.cxx.

84 : TNamed(),
85 fDetector(0),
88 fGapVolId(0),
90{
91 /// Default constructor
92}

◆ ~Ex03CalorimeterSD()

Ex03CalorimeterSD::~Ex03CalorimeterSD ( )
virtual

Destructor

Definition at line 95 of file Ex03CalorimeterSD.cxx.

96{
97 /// Destructor
98
99 if (fCalCollection) fCalCollection->Delete();
100 delete fCalCollection;
101}

Member Function Documentation

◆ Initialize()

void Ex03CalorimeterSD::Initialize ( )

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

Definition at line 130 of file Ex03CalorimeterSD.cxx.

131{
132 /// Register hits collection in the Root manager;
133 /// set sensitive volumes.
134
135 if (TMCRootManager::Instance()) Register();
136
137 // Keep the pointer to TVirtualMC object as a data member
138 // to avoid a possible performance penalty due to a frequent retrieval
139 // from the thread-local storage
140 fMC = gMC;
141
142 fAbsorberVolId = fMC->VolId("ABSO");
143 fGapVolId = fMC->VolId("GAPX");
144}

◆ ProcessHits()

Bool_t Ex03CalorimeterSD::ProcessHits ( )

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

Definition at line 147 of file Ex03CalorimeterSD.cxx.

148{
149 /// Account energy deposit and track lengths for each layer in its hit.
150
151 Int_t copyNo;
152 Int_t id = fMC->CurrentVolID(copyNo);
153
154 if (id != fAbsorberVolId && id != fGapVolId) return false;
155
156 fMC->CurrentVolOffID(2, copyNo);
157 // cout << "Got copyNo "<< copyNo << " " << fMC->CurrentVolPath() << endl;
158
159 Double_t edep = fMC->Edep();
160
161 Double_t step = 0.;
162 if (fMC->TrackCharge() != 0.) step = fMC->TrackStep();
163
164 if (!GetHit(copyNo)) {
165 std::cerr << "No hit found for layer with copyNo = " << copyNo << endl;
166 return false;
167 }
168
169 if (id == fAbsorberVolId) {
170 GetHit(copyNo)->AddAbs(edep, step);
171 }
172
173 if (id == fGapVolId) {
174 GetHit(copyNo)->AddGap(edep, step);
175 }
176
177 return true;
178}
void AddAbs(Double_t de, Double_t dl)
void AddGap(Double_t de, Double_t dl)
Ex03CalorHit * GetHit(Int_t i) const

◆ EndOfEvent()

void Ex03CalorimeterSD::EndOfEvent ( )

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

Definition at line 181 of file Ex03CalorimeterSD.cxx.

182{
183 /// Print hits collection (if verbose) and reset hits afterwards.
184
185 if (fVerboseLevel > 1) Print();
186
187 // Reset hits collection
188 ResetHits();
189}
virtual void Print(Option_t *option="") const

◆ Register()

void Ex03CalorimeterSD::Register ( )

Register the hits collection in Root manager.

Definition at line 192 of file Ex03CalorimeterSD.cxx.

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

◆ Print()

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

Print the hits collection.

Definition at line 200 of file Ex03CalorimeterSD.cxx.

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

◆ PrintTotal()

void Ex03CalorimeterSD::PrintTotal ( ) const

Print the total values for all layers.

Definition at line 212 of file Ex03CalorimeterSD.cxx.

213{
214 /// Print the total values for all layers.
215
216 Double_t totEAbs = 0.;
217 Double_t totLAbs = 0.;
218 Double_t totEGap = 0.;
219 Double_t totLGap = 0.;
220
221 Int_t nofHits = fCalCollection->GetEntriesFast();
222 for (Int_t i = 0; i < nofHits; i++) {
223 totEAbs += GetHit(i)->GetEdepAbs();
224 totLAbs += GetHit(i)->GetTrakAbs();
225 totEGap += GetHit(i)->GetEdepGap();
226 totLGap += GetHit(i)->GetTrakGap();
227 }
228
229 cout << " Absorber: total energy (MeV): " << setw(7) << totEAbs * 1.0e03
230 << " total track length (cm): " << setw(7) << totLAbs << endl
231 << " Gap: total energy (MeV): " << setw(7) << totEGap * 1.0e03
232 << " total track length (cm): " << setw(7) << totLGap << endl;
233}
Double_t GetEdepGap()
Double_t GetTrakGap()
Double_t GetEdepAbs()
Double_t GetTrakAbs()

◆ SetVerboseLevel()

void Ex03CalorimeterSD::SetVerboseLevel ( Int_t level)
inline

Set verbose level

Parameters
levelThe new verbose level value

Definition at line 75 of file Ex03CalorimeterSD.h.

76{
77 fVerboseLevel = level;
78}

◆ GetHit()

Ex03CalorHit * Ex03CalorimeterSD::GetHit ( Int_t i) const
Returns
The hit for the specified layer.
Parameters
iThe layer number

Definition at line 108 of file Ex03CalorimeterSD.cxx.

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

◆ ResetHits()

void Ex03CalorimeterSD::ResetHits ( )
private

Reset all hits in the hits collection.

Definition at line 117 of file Ex03CalorimeterSD.cxx.

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

Member Data Documentation

◆ fMC

TVirtualMC* Ex03CalorimeterSD::fMC
private

The VMC implementation.

Definition at line 63 of file Ex03CalorimeterSD.h.

◆ fDetector

Ex03DetectorConstruction* Ex03CalorimeterSD::fDetector
private

Detector construction.

Definition at line 64 of file Ex03CalorimeterSD.h.

◆ fCalCollection

TClonesArray* Ex03CalorimeterSD::fCalCollection
private

Hits collection.

Definition at line 65 of file Ex03CalorimeterSD.h.

◆ fAbsorberVolId

Int_t Ex03CalorimeterSD::fAbsorberVolId
private

The absorber volume Id.

Definition at line 66 of file Ex03CalorimeterSD.h.

◆ fGapVolId

Int_t Ex03CalorimeterSD::fGapVolId
private

The gap volume Id.

Definition at line 67 of file Ex03CalorimeterSD.h.

◆ fVerboseLevel

Int_t Ex03CalorimeterSD::fVerboseLevel
private

Verbosity level.

Definition at line 68 of file Ex03CalorimeterSD.h.


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