Geant4 11.4.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4PSCellFlux Class Reference

#include <G4PSCellFlux.hh>

Inheritance diagram for G4PSCellFlux:

Public Member Functions

 G4PSCellFlux (const G4String &name, G4int depth=0)
 G4PSCellFlux (const G4String &name, const G4String &unit, G4int depth=0)
 ~G4PSCellFlux () override=default
void Weighted (G4bool flg=true)
void Initialize (G4HCofThisEvent *) override
void clear () override
void PrintAll () override
virtual void SetUnit (const G4String &unit)
Public Member Functions inherited from G4VPrimitivePlotter
 ~G4VPrimitivePlotter () override=default
void Plot (G4int copyNo, G4int histID)
G4int GetNumberOfHist () const
 G4VPrimitiveScorer (const G4String &name, G4int depth=0)
Public Member Functions inherited from G4VPrimitiveScorer
 G4VPrimitiveScorer (const G4String &name, G4int depth=0)
virtual ~G4VPrimitiveScorer ()=default
G4int GetCollectionID (G4int)
virtual void EndOfEvent (G4HCofThisEvent *)
virtual void DrawAll ()
void SetUnit (const G4String &unit)
const G4StringGetUnit () const
G4double GetUnitValue () const
void ScoreWeighted (G4bool flg=false)
G4bool IsScoreWeighted () const
void SetMultiFunctionalDetector (G4MultiFunctionalDetector *d)
G4MultiFunctionalDetectorGetMultiFunctionalDetector () const
const G4StringGetName () const
void SetFilter (G4VSDFilter *f)
G4VSDFilterGetFilter () const
void SetScoreWeightCalculator (G4ScoreWeightCalculator calculator)
void SetVerboseLevel (G4int vl)
G4int GetVerboseLevel () const
void SetNijk (G4int i, G4int j, G4int k)

Protected Member Functions

G4bool ProcessHits (G4Step *, G4TouchableHistory *) override
virtual G4double ComputeVolume (G4Step *, G4int idx)
virtual void DefineUnitAndCategory ()
Protected Member Functions inherited from G4VPrimitiveScorer
G4VSolidComputeSolid (G4Step *aStep, G4int replicaIdx)
G4VSolidComputeCurrentSolid (G4Step *aStep)
virtual G4int GetIndex (G4Step *)
void CheckAndSetUnit (const G4String &unit, const G4String &category)

Additional Inherited Members

Protected Attributes inherited from G4VPrimitivePlotter
std::map< G4int, G4inthitIDMap
Protected Attributes inherited from G4VPrimitiveScorer
G4String primitiveName
G4MultiFunctionalDetectordetector {nullptr}
G4VSDFilterfilter {nullptr}
G4int verboseLevel {0}
G4int indexDepth
G4String unitName {"NoUnit"}
G4double unitValue {1.0}
G4int fNi {0}
G4int fNj {0}
G4int fNk {0}
G4bool scoreWeighted {false}
G4ScoreWeightCalculator fScoreWeightCalculator

Detailed Description

Definition at line 57 of file G4PSCellFlux.hh.

Constructor & Destructor Documentation

◆ G4PSCellFlux() [1/2]

G4PSCellFlux::G4PSCellFlux ( const G4String & name,
G4int depth = 0 )

Definition at line 58 of file G4PSCellFlux.cc.

59 : G4PSCellFlux(name, "percm2", depth)
60{}
G4PSCellFlux(const G4String &name, G4int depth=0)

Referenced by G4PSCellFlux(), and G4PSCellFlux3D::G4PSCellFlux3D().

◆ G4PSCellFlux() [2/2]

G4PSCellFlux::G4PSCellFlux ( const G4String & name,
const G4String & unit,
G4int depth = 0 )

Definition at line 62 of file G4PSCellFlux.cc.

63 : G4VPrimitivePlotter(name, depth)
64 , HCID(-1)
65 , EvtMap(nullptr)
66 , weighted(true)
67{
69 SetUnit(unit);
70}
virtual void SetUnit(const G4String &unit)
virtual void DefineUnitAndCategory()

◆ ~G4PSCellFlux()

G4PSCellFlux::~G4PSCellFlux ( )
overridedefault

Member Function Documentation

◆ clear()

void G4PSCellFlux::clear ( )
overridevirtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 117 of file G4PSCellFlux.cc.

117{ EvtMap->clear(); }

◆ ComputeVolume()

G4double G4PSCellFlux::ComputeVolume ( G4Step * aStep,
G4int idx )
protectedvirtual

Reimplemented in G4PSCellFluxForCylinder3D.

Definition at line 147 of file G4PSCellFlux.cc.

148{
149 G4VSolid* solid = ComputeSolid(aStep, idx);
150 assert(solid);
151 return solid->GetCubicVolume();
152}
G4VSolid * ComputeSolid(G4Step *aStep, G4int replicaIdx)
virtual G4double GetCubicVolume()
Definition G4VSolid.cc:214

Referenced by ProcessHits().

◆ DefineUnitAndCategory()

void G4PSCellFlux::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 137 of file G4PSCellFlux.cc.

138{
139 // Per Unit Surface
140 new G4UnitDefinition("percentimeter2", "percm2", "Per Unit Surface",
141 (1. / cm2));
142 new G4UnitDefinition("permillimeter2", "permm2", "Per Unit Surface",
143 (1. / mm2));
144 new G4UnitDefinition("permeter2", "perm2", "Per Unit Surface", (1. / m2));
145}

Referenced by G4PSCellFlux().

◆ Initialize()

void G4PSCellFlux::Initialize ( G4HCofThisEvent * HCE)
overridevirtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 109 of file G4PSCellFlux.cc.

110{
111 EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName());
112 if(HCID < 0)
113 HCID = GetCollectionID(0);
114 HCE->AddHitsCollection(HCID, EvtMap);
115}
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)
const G4String & GetName() const
G4MultiFunctionalDetector * detector

◆ PrintAll()

void G4PSCellFlux::PrintAll ( )
overridevirtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 119 of file G4PSCellFlux.cc.

120{
121 G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
122 G4cout << " PrimitiveScorer " << GetName() << G4endl;
123 G4cout << " Number of entries " << EvtMap->entries() << G4endl;
124 for(const auto& [copy, flux]: *(EvtMap->GetMap()))
125 {
126 G4cout << " copy no.: " << copy
127 << " cell flux : " << *(flux) / GetUnitValue() << " ["
128 << GetUnit() << "]" << G4endl;
129 }
130}
#define G4endl
Definition G4ios.hh:67
G4GLOB_DLL std::ostream G4cout
const G4String & GetUnit() const
G4double GetUnitValue() const
void copy(G4double dst[], const G4double src[], std::size_t size=G4FieldTrack::ncompSVEC)

◆ ProcessHits()

G4bool G4PSCellFlux::ProcessHits ( G4Step * aStep,
G4TouchableHistory *  )
overrideprotectedvirtual

Implements G4VPrimitiveScorer.

Definition at line 72 of file G4PSCellFlux.cc.

73{
74 G4double stepLength = aStep->GetStepLength();
75 if(stepLength == 0.)
76 return false;
77
78 G4int idx = ((G4TouchableHistory*) (aStep->GetPreStepPoint()->GetTouchable()))
79 ->GetReplicaNumber(indexDepth);
80 G4double cubicVolume = ComputeVolume(aStep, idx);
81 G4double CellFlux = stepLength / cubicVolume;
82 if (scoreWeighted) {
83 CellFlux *= fScoreWeightCalculator(aStep);
84 }
85 if(weighted)
86 CellFlux *= aStep->GetPreStepPoint()->GetWeight();
87 G4int index = GetIndex(aStep);
88 EvtMap->add(index, CellFlux);
89
90 if(!hitIDMap.empty() && hitIDMap.find(index) != hitIDMap.end())
91 {
92 auto filler = G4VScoreHistFiller::Instance();
93 if(filler == nullptr)
94 {
96 "G4PSCellFlux::ProcessHits", "SCORER0123", JustWarning,
97 "G4TScoreHistFiller is not instantiated!! Histogram is not filled.");
98 }
99 else
100 {
101 filler->FillH1(hitIDMap[index],
102 aStep->GetPreStepPoint()->GetKineticEnergy(), CellFlux);
103 }
104 }
105
106 return true;
107}
@ JustWarning
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
double G4double
Definition G4Types.hh:83
int G4int
Definition G4Types.hh:85
virtual G4double ComputeVolume(G4Step *, G4int idx)
const G4VTouchable * GetTouchable() const
G4double GetKineticEnergy() const
G4double GetWeight() const
G4StepPoint * GetPreStepPoint() const
G4double GetStepLength() const
std::map< G4int, G4int > hitIDMap
virtual G4int GetIndex(G4Step *)
G4ScoreWeightCalculator fScoreWeightCalculator
static G4VScoreHistFiller * Instance()

◆ SetUnit()

void G4PSCellFlux::SetUnit ( const G4String & unit)
virtual

Definition at line 132 of file G4PSCellFlux.cc.

133{
134 CheckAndSetUnit(unit, "Per Unit Surface");
135}
void CheckAndSetUnit(const G4String &unit, const G4String &category)

Referenced by G4PSCellFlux(), G4PSCellFlux3D::G4PSCellFlux3D(), and G4ScoreQuantityMessenger::SetNewValue().

◆ Weighted()

void G4PSCellFlux::Weighted ( G4bool flg = true)
inline

Definition at line 64 of file G4PSCellFlux.hh.

64{ weighted = flg; }

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