Geant4 11.2.2
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 (G4String name, G4int depth=0)
 
 G4PSCellFlux (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 (G4String name, G4int depth=0)
 
- Public Member Functions inherited from G4VPrimitiveScorer
 G4VPrimitiveScorer (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 SetMultiFunctionalDetector (G4MultiFunctionalDetector *d)
 
G4MultiFunctionalDetectorGetMultiFunctionalDetector () const
 
G4String GetName () const
 
void SetFilter (G4VSDFilter *f)
 
G4VSDFilterGetFilter () const
 
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}
 

Detailed Description

Definition at line 57 of file G4PSCellFlux.hh.

Constructor & Destructor Documentation

◆ G4PSCellFlux() [1/2]

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

Definition at line 58 of file G4PSCellFlux.cc.

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

◆ G4PSCellFlux() [2/2]

G4PSCellFlux::G4PSCellFlux ( 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 115 of file G4PSCellFlux.cc.

115{ EvtMap->clear(); }

◆ ComputeVolume()

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

Reimplemented in G4PSCellFluxForCylinder3D.

Definition at line 145 of file G4PSCellFlux.cc.

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

Referenced by ProcessHits().

◆ DefineUnitAndCategory()

void G4PSCellFlux::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 135 of file G4PSCellFlux.cc.

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

Referenced by G4PSCellFlux().

◆ Initialize()

void G4PSCellFlux::Initialize ( G4HCofThisEvent * HCE)
overridevirtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 107 of file G4PSCellFlux.cc.

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

◆ PrintAll()

void G4PSCellFlux::PrintAll ( )
overridevirtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 117 of file G4PSCellFlux.cc.

118{
119 G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
120 G4cout << " PrimitiveScorer " << GetName() << G4endl;
121 G4cout << " Number of entries " << EvtMap->entries() << G4endl;
122 for(const auto& [copy, flux]: *(EvtMap->GetMap()))
123 {
124 G4cout << " copy no.: " << copy
125 << " cell flux : " << *(flux) / GetUnitValue() << " ["
126 << GetUnit() << "]" << G4endl;
127 }
128}
#define G4endl
Definition G4ios.hh:67
G4GLOB_DLL std::ostream G4cout
const G4String & GetUnit() const
G4double GetUnitValue() const
Map_t * GetMap() const
size_t entries() 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
79 ->GetReplicaNumber(indexDepth);
80 G4double cubicVolume = ComputeVolume(aStep, idx);
81
82 G4double CellFlux = stepLength / cubicVolume;
83 if(weighted)
84 CellFlux *= aStep->GetPreStepPoint()->GetWeight();
85 G4int index = GetIndex(aStep);
86 EvtMap->add(index, CellFlux);
87
88 if(!hitIDMap.empty() && hitIDMap.find(index) != hitIDMap.end())
89 {
90 auto filler = G4VScoreHistFiller::Instance();
91 if(filler == nullptr)
92 {
94 "G4PSCellFlux::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(), CellFlux);
101 }
102 }
103
104 return true;
105}
@ 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 *)
static G4VScoreHistFiller * Instance()
size_t add(const G4int &key, U *&aHit) const

◆ SetUnit()

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

Definition at line 130 of file G4PSCellFlux.cc.

131{
132 CheckAndSetUnit(unit, "Per Unit Surface");
133}
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: