Geant4 11.1.1
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 (G4String name, G4int depth=0)
 
virtual ~G4VPrimitivePlotter ()
 
void Plot (G4int copyNo, G4int histID)
 
G4int GetNumberOfHist () const
 
- Public Member Functions inherited from G4VPrimitiveScorer
 G4VPrimitiveScorer (G4String name, G4int depth=0)
 
virtual ~G4VPrimitiveScorer ()
 
G4int GetCollectionID (G4int)
 
virtual void Initialize (G4HCofThisEvent *)
 
virtual void EndOfEvent (G4HCofThisEvent *)
 
virtual void clear ()
 
virtual void DrawAll ()
 
virtual void PrintAll ()
 
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
virtual G4bool ProcessHits (G4Step *, G4TouchableHistory *)=0
 
virtual G4int GetIndex (G4Step *)
 
void CheckAndSetUnit (const G4String &unit, const G4String &category)
 
G4VSolidComputeSolid (G4Step *aStep, G4int replicaIdx)
 
G4VSolidComputeCurrentSolid (G4Step *aStep)
 

Additional Inherited Members

- Protected Attributes inherited from G4VPrimitivePlotter
std::map< G4int, G4inthitIDMap
 
- Protected Attributes inherited from G4VPrimitiveScorer
G4String primitiveName
 
G4MultiFunctionalDetectordetector
 
G4VSDFilterfilter
 
G4int verboseLevel
 
G4int indexDepth
 
G4String unitName
 
G4double unitValue
 
G4int fNi
 
G4int fNj
 
G4int fNk
 

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() [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(); }
void clear()
Definition: G4THitsMap.hh:524

◆ 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
G4int GetCollectionID(G4int)

◆ 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:57
G4GLOB_DLL std::ostream G4cout
const G4String & GetUnit() const
G4double GetUnitValue() const
Map_t * GetMap() const
Definition: G4THitsMap.hh:155
size_t entries() const
Definition: G4THitsMap.hh:158
void copy(G4double dst[], const G4double src[], std::size_t size=G4FieldTrack::ncompSVEC)
Definition: G4FieldUtils.cc:98

◆ 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)
Definition: G4Exception.cc:59
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
Definition: G4THitsMap.hh:177

◆ 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: