Geant4-11
G4PSTrackCounter.cc
Go to the documentation of this file.
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26//
27//
28// G4PSTrackCounter
29#include "G4PSTrackCounter.hh"
30#include "G4UnitsTable.hh"
31#include "G4VScoreHistFiller.hh"
32
34// (Description)
35// This is a primitive scorer class for scoring number of tracks in a cell.
36//
37// Created: 2007-02-02 Tsukasa ASO, Akinori Kimura.
38// 2010-07-22 Introduce Unit specification.
39//
41
44 , HCID(-1)
45 , fDirection(direction)
46 , EvtMap(0)
47 , weighted(false)
48{
49 SetUnit("");
50}
51
53
55{
56 G4StepPoint* preStep = aStep->GetPreStepPoint();
57 G4StepPoint* posStep = aStep->GetPostStepPoint();
58 G4bool IsEnter = preStep->GetStepStatus() == fGeomBoundary;
59 G4bool IsExit = posStep->GetStepStatus() == fGeomBoundary;
60
61 // Regular : count in prestep volume.
62 G4int index = GetIndex(aStep);
63
64 // G4cout << " trk " << aStep->GetTrack()->GetTrackID()
65 // << " index " << index << " In " << IsEnter << " Out " <<IsExit
66 // << " PreStatus " << preStep->GetStepStatus()
67 // << " PostStatus " << posStep->GetStepStatus()
68 // << " w " << aStep->GetPreStepPoint()->GetWeight()
69 // << G4endl;
70
71 G4bool flag = FALSE;
72
73 if(IsEnter && fDirection == fCurrent_In)
74 flag = TRUE;
75 else if(IsExit && fDirection == fCurrent_Out)
76 flag = TRUE;
77 else if((IsExit || IsEnter) && fDirection == fCurrent_InOut)
78 flag = TRUE;
79
80 if(flag)
81 {
82 // G4cout << " Count " << G4endl;
83 G4double val = 1.0;
84 if(weighted)
85 val *= aStep->GetPreStepPoint()->GetWeight();
86 EvtMap->add(index, val);
87
88 if(hitIDMap.size() > 0 && hitIDMap.find(index) != hitIDMap.end())
89 {
90 auto filler = G4VScoreHistFiller::Instance();
91 if(!filler)
92 {
94 "G4PSTrackCounter::ProcessHits", "SCORER0123", JustWarning,
95 "G4TScoreHistFiller is not instantiated!! Histogram is not filled.");
96 }
97 else
98 {
99 filler->FillH1(hitIDMap[index],
100 aStep->GetPreStepPoint()->GetKineticEnergy(), val);
101 }
102 }
103 }
104
105 return TRUE;
106}
107
109{
111 if(HCID < 0)
112 {
114 }
116 // G4cout << "----- PS Initialize " << G4endl;
117}
118
120
122
124
126{
127 G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
128 G4cout << " PrimitiveScorer " << GetName() << G4endl;
129 G4cout << " Number of entries " << EvtMap->entries() << G4endl;
130 std::map<G4int, G4double*>::iterator itr = EvtMap->GetMap()->begin();
131 for(; itr != EvtMap->GetMap()->end(); itr++)
132 {
133 G4cout << " copy no.: " << itr->first
134 << " track count: " << *(itr->second) << " [tracks] " << G4endl;
135 }
136}
137
139{
140 if(unit == "")
141 {
142 unitName = unit;
143 unitValue = 1.0;
144 }
145 else
146 {
147 G4String msg = "Invalid unit [" + unit + "] (Current unit is [" +
148 GetUnit() + "] ) for " + GetName();
149 G4Exception("G4PSTrackCounter::SetUnit", "DetPS0018", JustWarning, msg);
150 }
151}
@ JustWarning
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
@ fCurrent_Out
@ fCurrent_In
@ fCurrent_InOut
@ fGeomBoundary
Definition: G4StepStatus.hh:43
double G4double
Definition: G4Types.hh:83
bool G4bool
Definition: G4Types.hh:86
int G4int
Definition: G4Types.hh:85
#define G4endl
Definition: G4ios.hh:57
G4GLOB_DLL std::ostream G4cout
#define TRUE
Definition: Globals.hh:27
#define FALSE
Definition: Globals.hh:23
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)
virtual void SetUnit(const G4String &unit)
G4THitsMap< G4double > * EvtMap
virtual G4bool ProcessHits(G4Step *, G4TouchableHistory *)
virtual void EndOfEvent(G4HCofThisEvent *)
virtual void DrawAll()
G4PSTrackCounter(G4String name, G4int direction, G4int depth=0)
virtual void PrintAll()
virtual ~G4PSTrackCounter()
virtual void clear()
virtual void Initialize(G4HCofThisEvent *)
G4StepStatus GetStepStatus() const
G4double GetKineticEnergy() const
G4double GetWeight() const
Definition: G4Step.hh:62
G4StepPoint * GetPreStepPoint() const
G4StepPoint * GetPostStepPoint() const
std::map< G4int, G4int > hitIDMap
virtual G4int GetIndex(G4Step *)
G4String GetName() const
const G4String & GetUnit() const
G4MultiFunctionalDetector * detector
G4int GetCollectionID(G4int)
static G4VScoreHistFiller * Instance()
Map_t * GetMap() const
Definition: G4THitsMap.hh:155
size_t entries() const
Definition: G4THitsMap.hh:158
void clear()
Definition: G4THitsMap.hh:524
size_t add(const G4int &key, U *&aHit) const
Definition: G4THitsMap.hh:177
const char * name(G4int ptype)