Geant4-11
Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes
G4PSCylinderSurfaceFlux Class Reference

#include <G4PSCylinderSurfaceFlux.hh>

Inheritance diagram for G4PSCylinderSurfaceFlux:
G4VPrimitivePlotter G4VPrimitiveScorer G4PSCylinderSurfaceFlux3D

Public Member Functions

virtual void clear ()
 
void DivideByArea (G4bool flg=true)
 
virtual void DrawAll ()
 
virtual void EndOfEvent (G4HCofThisEvent *)
 
 G4PSCylinderSurfaceFlux (G4String name, G4int direction, const G4String &unit, G4int depth=0)
 
 G4PSCylinderSurfaceFlux (G4String name, G4int direction, G4int depth=0)
 
G4int GetCollectionID (G4int)
 
G4VSDFilterGetFilter () const
 
G4MultiFunctionalDetectorGetMultiFunctionalDetector () const
 
G4String GetName () const
 
G4int GetNumberOfHist () const
 
const G4StringGetUnit () const
 
G4double GetUnitValue () const
 
G4int GetVerboseLevel () const
 
virtual void Initialize (G4HCofThisEvent *)
 
void Plot (G4int copyNo, G4int histID)
 
virtual void PrintAll ()
 
void SetFilter (G4VSDFilter *f)
 
void SetMultiFunctionalDetector (G4MultiFunctionalDetector *d)
 
void SetNijk (G4int i, G4int j, G4int k)
 
virtual void SetUnit (const G4String &unit)
 
void SetVerboseLevel (G4int vl)
 
void Weighted (G4bool flg=true)
 
virtual ~G4PSCylinderSurfaceFlux ()
 

Protected Member Functions

void CheckAndSetUnit (const G4String &unit, const G4String &category)
 
G4VSolidComputeCurrentSolid (G4Step *aStep)
 
G4VSolidComputeSolid (G4Step *aStep, G4int replicaIdx)
 
virtual void DefineUnitAndCategory ()
 
virtual G4int GetIndex (G4Step *)
 
G4int IsSelectedSurface (G4Step *, G4Tubs *)
 
virtual G4bool ProcessHits (G4Step *, G4TouchableHistory *)
 

Protected Attributes

G4MultiFunctionalDetectordetector
 
G4VSDFilterfilter
 
G4int fNi
 
G4int fNj
 
G4int fNk
 
std::map< G4int, G4inthitIDMap
 
G4int indexDepth
 
G4String primitiveName
 
G4String unitName
 
G4double unitValue
 
G4int verboseLevel
 

Private Member Functions

G4bool HitPrimitive (G4Step *aStep, G4TouchableHistory *ROhis)
 

Private Attributes

G4bool divideByArea
 
G4THitsMap< G4double > * EvtMap
 
G4int fDirection
 
G4int HCID
 
G4bool weighted
 

Detailed Description

Definition at line 60 of file G4PSCylinderSurfaceFlux.hh.

Constructor & Destructor Documentation

◆ G4PSCylinderSurfaceFlux() [1/2]

G4PSCylinderSurfaceFlux::G4PSCylinderSurfaceFlux ( G4String  name,
G4int  direction,
G4int  depth = 0 
)

Definition at line 62 of file G4PSCylinderSurfaceFlux.cc.

65 , HCID(-1)
66 , fDirection(direction)
67 , EvtMap(0)
68 , weighted(true)
69 , divideByArea(true)
70{
72 SetUnit("percm2");
73}
G4THitsMap< G4double > * EvtMap
virtual void SetUnit(const G4String &unit)
G4VPrimitivePlotter(G4String name, G4int depth=0)
const char * name(G4int ptype)

References DefineUnitAndCategory(), and SetUnit().

◆ G4PSCylinderSurfaceFlux() [2/2]

G4PSCylinderSurfaceFlux::G4PSCylinderSurfaceFlux ( G4String  name,
G4int  direction,
const G4String unit,
G4int  depth = 0 
)

Definition at line 75 of file G4PSCylinderSurfaceFlux.cc.

79 , HCID(-1)
80 , fDirection(direction)
81 , EvtMap(0)
82 , weighted(true)
83 , divideByArea(true)
84{
86 SetUnit(unit);
87}

References DefineUnitAndCategory(), and SetUnit().

◆ ~G4PSCylinderSurfaceFlux()

G4PSCylinderSurfaceFlux::~G4PSCylinderSurfaceFlux ( )
virtual

Definition at line 89 of file G4PSCylinderSurfaceFlux.cc.

89{ ; }

Member Function Documentation

◆ CheckAndSetUnit()

void G4VPrimitiveScorer::CheckAndSetUnit ( const G4String unit,
const G4String category 
)
protectedinherited

Definition at line 81 of file G4VPrimitiveScorer.cc.

83{
84 if(G4UnitDefinition::GetCategory(unit) == category)
85 {
86 unitName = unit;
88 }
89 else
90 {
91 G4String msg = "Invalid unit [" + unit + "] (Current unit is [" +
92 GetUnit() + "] ) requested for " + GetName();
93 G4Exception("G4VPrimitiveScorer::CheckAndSetUnit", "Det0151", JustWarning,
94 msg);
95 }
96}
@ JustWarning
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:35
static G4double GetValueOf(const G4String &)
static G4String GetCategory(const G4String &)
G4String GetName() const
const G4String & GetUnit() const

References G4Exception(), G4UnitDefinition::GetCategory(), G4VPrimitiveScorer::GetName(), G4VPrimitiveScorer::GetUnit(), G4UnitDefinition::GetValueOf(), JustWarning, G4VPrimitiveScorer::unitName, and G4VPrimitiveScorer::unitValue.

Referenced by G4PSCellCharge::SetUnit(), G4PSCellFlux::SetUnit(), G4PSCylinderSurfaceCurrent::SetUnit(), SetUnit(), G4PSDoseDeposit::SetUnit(), G4PSEnergyDeposit::SetUnit(), G4PSFlatSurfaceCurrent::SetUnit(), G4PSFlatSurfaceFlux::SetUnit(), G4PSMinKinEAtGeneration::SetUnit(), G4PSPassageCellFlux::SetUnit(), G4PSPassageTrackLength::SetUnit(), G4PSSphereSurfaceCurrent::SetUnit(), G4PSSphereSurfaceFlux::SetUnit(), and G4PSTrackLength::SetUnit().

◆ clear()

void G4PSCylinderSurfaceFlux::clear ( )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 247 of file G4PSCylinderSurfaceFlux.cc.

247{ EvtMap->clear(); }
void clear()
Definition: G4THitsMap.hh:524

References G4VTHitsMap< T, Map_t >::clear(), and EvtMap.

◆ ComputeCurrentSolid()

G4VSolid * G4VPrimitiveScorer::ComputeCurrentSolid ( G4Step aStep)
protectedinherited

Definition at line 128 of file G4VPrimitiveScorer.cc.

129{
130 G4StepPoint* preStep = aStep->GetPreStepPoint();
131 // The only difference: did not know the replica number
132 G4int replicaIdx =
133 (static_cast<const G4TouchableHistory*>(preStep->GetTouchable()))
134 ->GetReplicaNumber(indexDepth);
135
136 return ComputeSolid(aStep, replicaIdx);
137}
int G4int
Definition: G4Types.hh:85
const G4VTouchable * GetTouchable() const
G4StepPoint * GetPreStepPoint() const
G4VSolid * ComputeSolid(G4Step *aStep, G4int replicaIdx)

References G4VPrimitiveScorer::ComputeSolid(), G4Step::GetPreStepPoint(), G4StepPoint::GetTouchable(), and G4VPrimitiveScorer::indexDepth.

Referenced by G4PSCylinderSurfaceCurrent::ProcessHits(), ProcessHits(), and G4PSSphereSurfaceCurrent::ProcessHits().

◆ ComputeSolid()

G4VSolid * G4VPrimitiveScorer::ComputeSolid ( G4Step aStep,
G4int  replicaIdx 
)
protectedinherited

Definition at line 98 of file G4VPrimitiveScorer.cc.

99{
100 G4VSolid* solid = nullptr;
101 G4StepPoint* preStep = aStep->GetPreStepPoint();
102
103 auto physVol = preStep->GetPhysicalVolume();
104 G4VPVParameterisation* physParam = physVol->GetParameterisation();
105
106 if(physParam)
107 { // for parameterized volume
108 if(replicaIdx < 0)
109 {
111 desc << "Incorrect replica number --- GetReplicaNumber : " << replicaIdx
112 << G4endl;
113 G4Exception("G4VPrimitiveScorer::ComputeSolid", "DetPS0001", JustWarning,
114 desc);
115 // replicaIdx= 0; // You must ensure that it's in range !!!
116 }
117 solid = physParam->ComputeSolid(replicaIdx, physVol);
118 solid->ComputeDimensions(physParam, replicaIdx, physVol);
119 }
120 else
121 { // for ordinary volume
122 solid = physVol->GetLogicalVolume()->GetSolid();
123 }
124
125 return solid;
126}
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
#define G4endl
Definition: G4ios.hh:57
G4VPhysicalVolume * GetPhysicalVolume() const
virtual G4VSolid * ComputeSolid(const G4int, G4VPhysicalVolume *)
virtual void ComputeDimensions(G4VPVParameterisation *p, const G4int n, const G4VPhysicalVolume *pRep)
Definition: G4VSolid.cc:137

References G4VSolid::ComputeDimensions(), G4VPVParameterisation::ComputeSolid(), G4endl, G4Exception(), G4StepPoint::GetPhysicalVolume(), G4Step::GetPreStepPoint(), and JustWarning.

Referenced by G4VPrimitiveScorer::ComputeCurrentSolid(), G4PSCellFlux::ComputeVolume(), G4PSDoseDeposit::ComputeVolume(), and G4PSPassageCellFlux::ComputeVolume().

◆ DefineUnitAndCategory()

void G4PSCylinderSurfaceFlux::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 288 of file G4PSCylinderSurfaceFlux.cc.

289{
290 // Per Unit Surface
291 new G4UnitDefinition("percentimeter2", "percm2", "Per Unit Surface",
292 (1. / cm2));
293 new G4UnitDefinition("permillimeter2", "permm2", "Per Unit Surface",
294 (1. / mm2));
295 new G4UnitDefinition("permeter2", "perm2", "Per Unit Surface", (1. / m2));
296}
static constexpr double mm2
Definition: G4SIunits.hh:96
static constexpr double cm2
Definition: G4SIunits.hh:100
static constexpr double m2
Definition: G4SIunits.hh:110

References cm2, m2, and mm2.

Referenced by G4PSCylinderSurfaceFlux().

◆ DivideByArea()

void G4PSCylinderSurfaceFlux::DivideByArea ( G4bool  flg = true)
inline

Definition at line 71 of file G4PSCylinderSurfaceFlux.hh.

71{ divideByArea = flg; }

References divideByArea.

◆ DrawAll()

void G4PSCylinderSurfaceFlux::DrawAll ( )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 249 of file G4PSCylinderSurfaceFlux.cc.

249{ ; }

◆ EndOfEvent()

void G4PSCylinderSurfaceFlux::EndOfEvent ( G4HCofThisEvent )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 245 of file G4PSCylinderSurfaceFlux.cc.

245{ ; }

◆ GetCollectionID()

G4int G4VPrimitiveScorer::GetCollectionID ( G4int  )
inherited

Definition at line 55 of file G4VPrimitiveScorer.cc.

56{
57 if(detector)
59 "/" + primitiveName);
60 else
61 return -1;
62}
static G4SDManager * GetSDMpointer()
Definition: G4SDManager.cc:38
G4int GetCollectionID(G4String colName)
Definition: G4SDManager.cc:142
G4MultiFunctionalDetector * detector

References G4VPrimitiveScorer::detector, G4SDManager::GetCollectionID(), G4VSensitiveDetector::GetName(), G4SDManager::GetSDMpointer(), and G4VPrimitiveScorer::primitiveName.

Referenced by G4PSCellCharge::Initialize(), G4PSCellFlux::Initialize(), G4PSCylinderSurfaceCurrent::Initialize(), Initialize(), G4PSDoseDeposit::Initialize(), G4PSEnergyDeposit::Initialize(), G4PSFlatSurfaceCurrent::Initialize(), G4PSFlatSurfaceFlux::Initialize(), G4PSMinKinEAtGeneration::Initialize(), G4PSNofCollision::Initialize(), G4PSNofSecondary::Initialize(), G4PSNofStep::Initialize(), G4PSPassageCellCurrent::Initialize(), G4PSPassageCellFlux::Initialize(), G4PSPassageTrackLength::Initialize(), G4PSPopulation::Initialize(), G4PSSphereSurfaceCurrent::Initialize(), G4PSSphereSurfaceFlux::Initialize(), G4PSTermination::Initialize(), G4PSTrackCounter::Initialize(), G4PSTrackLength::Initialize(), and G4PSVolumeFlux::Initialize().

◆ GetFilter()

G4VSDFilter * G4VPrimitiveScorer::GetFilter ( ) const
inlineinherited

Definition at line 108 of file G4VPrimitiveScorer.hh.

108{ return filter; }

References G4VPrimitiveScorer::filter.

Referenced by G4VScoringMesh::List(), and G4VScoringMesh::SetFilter().

◆ GetIndex()

G4int G4VPrimitiveScorer::GetIndex ( G4Step aStep)
protectedvirtualinherited

Reimplemented in G4PSCellCharge3D, G4PSCellFlux3D, G4PSCylinderSurfaceCurrent3D, G4PSCylinderSurfaceFlux3D, G4PSDoseDeposit3D, G4PSEnergyDeposit3D, G4PSFlatSurfaceCurrent3D, G4PSFlatSurfaceFlux3D, G4PSMinKinEAtGeneration3D, G4PSNofCollision3D, G4PSNofSecondary3D, G4PSNofStep3D, G4PSPassageCellCurrent3D, G4PSPassageCellFlux3D, G4PSPassageTrackLength3D, G4PSPopulation3D, G4PSSphereSurfaceCurrent3D, G4PSSphereSurfaceFlux3D, G4PSStepChecker3D, G4PSTermination3D, G4PSTrackCounter3D, G4PSTrackLength3D, and G4PSVolumeFlux3D.

Definition at line 74 of file G4VPrimitiveScorer.cc.

75{
76 G4StepPoint* preStep = aStep->GetPreStepPoint();
78 return th->GetReplicaNumber(indexDepth);
79}
G4int GetReplicaNumber(G4int depth=0) const

References G4Step::GetPreStepPoint(), G4TouchableHistory::GetReplicaNumber(), G4StepPoint::GetTouchable(), and G4VPrimitiveScorer::indexDepth.

Referenced by G4PSCellCharge::ProcessHits(), G4PSCellFlux::ProcessHits(), G4PSCylinderSurfaceCurrent::ProcessHits(), ProcessHits(), G4PSDoseDeposit::ProcessHits(), G4PSEnergyDeposit::ProcessHits(), G4PSFlatSurfaceCurrent::ProcessHits(), G4PSFlatSurfaceFlux::ProcessHits(), G4PSMinKinEAtGeneration::ProcessHits(), G4PSNofCollision::ProcessHits(), G4PSNofSecondary::ProcessHits(), G4PSNofStep::ProcessHits(), G4PSPassageCellCurrent::ProcessHits(), G4PSPassageCellFlux::ProcessHits(), G4PSPassageTrackLength::ProcessHits(), G4PSPopulation::ProcessHits(), G4PSSphereSurfaceCurrent::ProcessHits(), G4PSSphereSurfaceFlux::ProcessHits(), G4PSStepChecker::ProcessHits(), G4PSTermination::ProcessHits(), G4PSTrackCounter::ProcessHits(), G4PSTrackLength::ProcessHits(), and G4PSVolumeFlux::ProcessHits().

◆ GetMultiFunctionalDetector()

G4MultiFunctionalDetector * G4VPrimitiveScorer::GetMultiFunctionalDetector ( ) const
inlineinherited

◆ GetName()

G4String G4VPrimitiveScorer::GetName ( ) const
inlineinherited

Definition at line 106 of file G4VPrimitiveScorer.hh.

106{ return primitiveName; }

References G4VPrimitiveScorer::primitiveName.

Referenced by G4VPrimitiveScorer::CheckAndSetUnit(), G4VScoringMesh::GetPrimitiveScorer(), G4PSCellCharge::Initialize(), G4PSCellFlux::Initialize(), G4PSCylinderSurfaceCurrent::Initialize(), Initialize(), G4PSDoseDeposit::Initialize(), G4PSEnergyDeposit::Initialize(), G4PSFlatSurfaceCurrent::Initialize(), G4PSFlatSurfaceFlux::Initialize(), G4PSMinKinEAtGeneration::Initialize(), G4PSNofCollision::Initialize(), G4PSNofSecondary::Initialize(), G4PSNofStep::Initialize(), G4PSPassageCellCurrent::Initialize(), G4PSPassageCellFlux::Initialize(), G4PSPassageTrackLength::Initialize(), G4PSPopulation::Initialize(), G4PSSphereSurfaceCurrent::Initialize(), G4PSSphereSurfaceFlux::Initialize(), G4PSTermination::Initialize(), G4PSTrackCounter::Initialize(), G4PSTrackLength::Initialize(), G4PSVolumeFlux::Initialize(), G4VScoringMesh::List(), G4PSCellCharge::PrintAll(), G4PSCellFlux::PrintAll(), G4PSCylinderSurfaceCurrent::PrintAll(), PrintAll(), G4PSDoseDeposit::PrintAll(), G4PSEnergyDeposit::PrintAll(), G4PSFlatSurfaceCurrent::PrintAll(), G4PSFlatSurfaceFlux::PrintAll(), G4PSMinKinEAtGeneration::PrintAll(), G4PSNofCollision::PrintAll(), G4PSNofSecondary::PrintAll(), G4PSNofStep::PrintAll(), G4PSPassageCellCurrent::PrintAll(), G4PSPassageCellFlux::PrintAll(), G4PSPassageTrackLength::PrintAll(), G4PSPopulation::PrintAll(), G4PSSphereSurfaceCurrent::PrintAll(), G4PSSphereSurfaceFlux::PrintAll(), G4PSTermination::PrintAll(), G4PSTrackCounter::PrintAll(), G4PSTrackLength::PrintAll(), G4PSVolumeFlux::PrintAll(), G4MultiFunctionalDetector::RegisterPrimitive(), G4MultiFunctionalDetector::RemovePrimitive(), G4VScoringMesh::SetFilter(), G4VScoringMesh::SetPrimitiveScorer(), G4PSCylinderSurfaceCurrent::SetUnit(), SetUnit(), G4PSFlatSurfaceCurrent::SetUnit(), G4PSFlatSurfaceFlux::SetUnit(), G4PSNofCollision::SetUnit(), G4PSNofSecondary::SetUnit(), G4PSNofStep::SetUnit(), G4PSPassageCellCurrent::SetUnit(), G4PSPopulation::SetUnit(), G4PSSphereSurfaceCurrent::SetUnit(), G4PSSphereSurfaceFlux::SetUnit(), G4PSTermination::SetUnit(), and G4PSTrackCounter::SetUnit().

◆ GetNumberOfHist()

G4int G4VPrimitivePlotter::GetNumberOfHist ( ) const
inlineinherited

Definition at line 51 of file G4VPrimitivePlotter.hh.

51{ return hitIDMap.size(); }
std::map< G4int, G4int > hitIDMap

References G4VPrimitivePlotter::hitIDMap.

◆ GetUnit()

const G4String & G4VPrimitiveScorer::GetUnit ( ) const
inlineinherited

◆ GetUnitValue()

G4double G4VPrimitiveScorer::GetUnitValue ( ) const
inlineinherited

◆ GetVerboseLevel()

G4int G4VPrimitiveScorer::GetVerboseLevel ( ) const
inlineinherited

Definition at line 110 of file G4VPrimitiveScorer.hh.

110{ return verboseLevel; }

References G4VPrimitiveScorer::verboseLevel.

◆ HitPrimitive()

G4bool G4VPrimitiveScorer::HitPrimitive ( G4Step aStep,
G4TouchableHistory ROhis 
)
inlineprivateinherited

Definition at line 113 of file G4VPrimitiveScorer.hh.

114 {
115 if(filter)
116 {
117 if(!(filter->Accept(aStep)))
118 return false;
119 }
120 return ProcessHits(aStep, ROhis);
121 }
virtual G4bool ProcessHits(G4Step *, G4TouchableHistory *)=0
virtual G4bool Accept(const G4Step *) const =0

References G4VSDFilter::Accept(), G4VPrimitiveScorer::filter, and G4VPrimitiveScorer::ProcessHits().

◆ Initialize()

void G4PSCylinderSurfaceFlux::Initialize ( G4HCofThisEvent HCE)
virtual

◆ IsSelectedSurface()

G4int G4PSCylinderSurfaceFlux::IsSelectedSurface ( G4Step aStep,
G4Tubs tubsSolid 
)
protected

Definition at line 179 of file G4PSCylinderSurfaceFlux.cc.

181{
182 G4TouchableHandle theTouchable =
186
188 {
189 // Entering Geometry
190 G4ThreeVector stppos1 = aStep->GetPreStepPoint()->GetPosition();
191 G4ThreeVector localpos1 =
192 theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos1);
193 if(std::fabs(localpos1.z()) > tubsSolid->GetZHalfLength())
194 return -1;
195 // if(std::fabs( localpos1.x()*localpos1.x()+localpos1.y()*localpos1.y()
196 // - (tubsSolid->GetInnerRadius()*tubsSolid->GetInnerRadius()))
197 // <kCarTolerance ){
198 G4double localR2 =
199 localpos1.x() * localpos1.x() + localpos1.y() * localpos1.y();
200 G4double InsideRadius = tubsSolid->GetInnerRadius();
201 if(localR2 >
202 (InsideRadius - kCarTolerance) * (InsideRadius - kCarTolerance) &&
203 localR2 <
204 (InsideRadius + kCarTolerance) * (InsideRadius + kCarTolerance))
205 {
206 return fFlux_In;
207 }
208 }
209
211 {
212 // Exiting Geometry
213 G4ThreeVector stppos2 = aStep->GetPostStepPoint()->GetPosition();
214 G4ThreeVector localpos2 =
215 theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos2);
216 if(std::fabs(localpos2.z()) > tubsSolid->GetZHalfLength())
217 return -1;
218 // if(std::fabs( localpos2.x()*localpos2.x()+localpos2.y()*localpos2.y()
219 // - (tubsSolid->GetInnerRadius()*tubsSolid->GetInnerRadius()))
220 // <kCarTolerance ){
221 G4double localR2 =
222 localpos2.x() * localpos2.x() + localpos2.y() * localpos2.y();
223 G4double InsideRadius = tubsSolid->GetInnerRadius();
224 if(localR2 >
225 (InsideRadius - kCarTolerance) * (InsideRadius - kCarTolerance) &&
226 localR2 <
227 (InsideRadius + kCarTolerance) * (InsideRadius + kCarTolerance))
228 {
229 return fFlux_Out;
230 }
231 }
232
233 return -1;
234}
const G4double kCarTolerance
@ fFlux_Out
@ fFlux_In
@ fGeomBoundary
Definition: G4StepStatus.hh:43
double G4double
Definition: G4Types.hh:83
double z() const
double x() const
double y() const
G4ThreeVector TransformPoint(const G4ThreeVector &vec) const
G4double GetSurfaceTolerance() const
static G4GeometryTolerance * GetInstance()
const G4AffineTransform & GetTopTransform() const
G4StepStatus GetStepStatus() const
const G4ThreeVector & GetPosition() const
const G4TouchableHandle & GetTouchableHandle() const
G4StepPoint * GetPostStepPoint() const
G4double GetZHalfLength() const
G4double GetInnerRadius() const
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:82

References fFlux_In, fFlux_Out, fGeomBoundary, G4VTouchable::GetHistory(), G4Tubs::GetInnerRadius(), G4GeometryTolerance::GetInstance(), G4StepPoint::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4StepPoint::GetStepStatus(), G4GeometryTolerance::GetSurfaceTolerance(), G4NavigationHistory::GetTopTransform(), G4StepPoint::GetTouchableHandle(), G4Tubs::GetZHalfLength(), kCarTolerance, G4AffineTransform::TransformPoint(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by ProcessHits().

◆ Plot()

void G4VPrimitivePlotter::Plot ( G4int  copyNo,
G4int  histID 
)
inlineinherited

Definition at line 45 of file G4VPrimitivePlotter.hh.

45{ hitIDMap[copyNo] = histID; }

References G4VPrimitivePlotter::hitIDMap.

◆ PrintAll()

void G4PSCylinderSurfaceFlux::PrintAll ( )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 251 of file G4PSCylinderSurfaceFlux.cc.

252{
253 G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
254 G4cout << " PrimitiveScorer" << GetName() << G4endl;
255 G4cout << " Number of entries " << EvtMap->entries() << G4endl;
256 std::map<G4int, G4double*>::iterator itr = EvtMap->GetMap()->begin();
257 for(; itr != EvtMap->GetMap()->end(); itr++)
258 {
259 G4cout << " copy no.: " << itr->first
260 << " flux : " << *(itr->second) / GetUnitValue() << " ["
261 << GetUnit() << "]" << G4endl;
262 }
263}
G4GLOB_DLL std::ostream G4cout
G4double GetUnitValue() const
Map_t * GetMap() const
Definition: G4THitsMap.hh:155
size_t entries() const
Definition: G4THitsMap.hh:158

References G4VPrimitiveScorer::detector, G4VTHitsMap< T, Map_t >::entries(), EvtMap, G4cout, G4endl, G4VTHitsMap< T, Map_t >::GetMap(), G4VPrimitiveScorer::GetName(), G4VSensitiveDetector::GetName(), G4VPrimitiveScorer::GetUnit(), and G4VPrimitiveScorer::GetUnitValue().

◆ ProcessHits()

G4bool G4PSCylinderSurfaceFlux::ProcessHits ( G4Step aStep,
G4TouchableHistory  
)
protectedvirtual

Implements G4VPrimitiveScorer.

Definition at line 91 of file G4PSCylinderSurfaceFlux.cc.

92{
93 G4StepPoint* preStep = aStep->GetPreStepPoint();
94 G4VSolid* solid = ComputeCurrentSolid(aStep);
95 assert(dynamic_cast<G4Tubs*>(solid));
96
97 G4Tubs* tubsSolid = static_cast<G4Tubs*>(solid);
98
99 G4int dirFlag = IsSelectedSurface(aStep, tubsSolid);
100
101 if(dirFlag > 0)
102 {
103 if(fDirection == fFlux_InOut || dirFlag == fDirection)
104 {
105 G4StepPoint* thisStep = 0;
106 if(dirFlag == fFlux_In)
107 {
108 thisStep = preStep;
109 }
110 else if(dirFlag == fFlux_Out)
111 {
112 thisStep = aStep->GetPostStepPoint();
113 }
114 else
115 {
116 return FALSE;
117 }
118
119 G4TouchableHandle theTouchable = thisStep->GetTouchableHandle();
120 G4ThreeVector pdirection = thisStep->GetMomentumDirection();
121 G4ThreeVector localdir =
122 theTouchable->GetHistory()->GetTopTransform().TransformAxis(pdirection);
123 G4ThreeVector position = thisStep->GetPosition();
124 G4ThreeVector localpos =
126 G4double angleFactor =
127 (localdir.x() * localpos.x() + localdir.y() * localpos.y()) /
128 std::sqrt(localdir.x() * localdir.x() + localdir.y() * localdir.y() +
129 localdir.z() * localdir.z()) /
130 std::sqrt(localpos.x() * localpos.x() + localpos.y() * localpos.y());
131
132 if(angleFactor < 0)
133 angleFactor *= -1.;
134 G4double square = 2. * tubsSolid->GetZHalfLength() *
135 tubsSolid->GetInnerRadius() *
136 tubsSolid->GetDeltaPhiAngle() / radian;
137
138 G4double flux = 1.0;
139 if(weighted)
140 flux *= preStep->GetWeight();
141 // Current (Particle Weight)
142
143 flux = flux / angleFactor;
144 if(divideByArea)
145 flux /= square;
146 // Flux with angle.
147 G4int index = GetIndex(aStep);
148 EvtMap->add(index, flux);
149
150 if(hitIDMap.size() > 0 && hitIDMap.find(index) != hitIDMap.end())
151 {
152 auto filler = G4VScoreHistFiller::Instance();
153 if(!filler)
154 {
155 G4Exception("G4PSCylinderSurfaceFlux::ProcessHits", "SCORER0123",
157 "G4TScoreHistFiller is not instantiated!! Histogram is "
158 "not filled.");
159 }
160 else
161 {
162 filler->FillH1(hitIDMap[index], thisStep->GetKineticEnergy(), flux);
163 }
164 }
165
166 return TRUE;
167 }
168 else
169 {
170 return FALSE;
171 }
172 }
173 else
174 {
175 return FALSE;
176 }
177}
@ fFlux_InOut
static constexpr double radian
Definition: G4SIunits.hh:122
#define TRUE
Definition: Globals.hh:27
#define FALSE
Definition: Globals.hh:23
G4ThreeVector TransformAxis(const G4ThreeVector &axis) const
G4int IsSelectedSurface(G4Step *, G4Tubs *)
const G4ThreeVector & GetMomentumDirection() const
G4double GetKineticEnergy() const
G4double GetWeight() const
Definition: G4Tubs.hh:75
G4double GetDeltaPhiAngle() const
virtual G4int GetIndex(G4Step *)
G4VSolid * ComputeCurrentSolid(G4Step *aStep)
static G4VScoreHistFiller * Instance()
size_t add(const G4int &key, U *&aHit) const
Definition: G4THitsMap.hh:177

References G4VTHitsMap< T, Map_t >::add(), G4VPrimitiveScorer::ComputeCurrentSolid(), divideByArea, EvtMap, FALSE, fDirection, fFlux_In, fFlux_InOut, fFlux_Out, G4Exception(), G4Tubs::GetDeltaPhiAngle(), G4VTouchable::GetHistory(), G4VPrimitiveScorer::GetIndex(), G4Tubs::GetInnerRadius(), G4StepPoint::GetKineticEnergy(), G4StepPoint::GetMomentumDirection(), G4StepPoint::GetPosition(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4NavigationHistory::GetTopTransform(), G4StepPoint::GetTouchableHandle(), G4StepPoint::GetWeight(), G4Tubs::GetZHalfLength(), G4VPrimitivePlotter::hitIDMap, G4VScoreHistFiller::Instance(), IsSelectedSurface(), JustWarning, radian, G4AffineTransform::TransformAxis(), TRUE, weighted, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

◆ SetFilter()

void G4VPrimitiveScorer::SetFilter ( G4VSDFilter f)
inlineinherited

Definition at line 107 of file G4VPrimitiveScorer.hh.

107{ filter = f; }

References G4VPrimitiveScorer::filter.

Referenced by G4VScoringMesh::SetFilter().

◆ SetMultiFunctionalDetector()

void G4VPrimitiveScorer::SetMultiFunctionalDetector ( G4MultiFunctionalDetector d)
inlineinherited

◆ SetNijk()

void G4VPrimitiveScorer::SetNijk ( G4int  i,
G4int  j,
G4int  k 
)
inlineinherited

◆ SetUnit()

void G4PSCylinderSurfaceFlux::SetUnit ( const G4String unit)
virtual

Definition at line 265 of file G4PSCylinderSurfaceFlux.cc.

266{
267 if(divideByArea)
268 {
269 CheckAndSetUnit(unit, "Per Unit Surface");
270 }
271 else
272 {
273 if(unit == "")
274 {
275 unitName = unit;
276 unitValue = 1.0;
277 }
278 else
279 {
280 G4String msg = "Invalid unit [" + unit + "] (Current unit is [" +
281 GetUnit() + "] ) for " + GetName();
282 G4Exception("G4PSCylinderSurfaceFlux::SetUnit", "DetPS0003", JustWarning,
283 msg);
284 }
285 }
286}
void CheckAndSetUnit(const G4String &unit, const G4String &category)

References G4VPrimitiveScorer::CheckAndSetUnit(), divideByArea, G4Exception(), G4VPrimitiveScorer::GetName(), G4VPrimitiveScorer::GetUnit(), JustWarning, G4VPrimitiveScorer::unitName, and G4VPrimitiveScorer::unitValue.

Referenced by G4PSCylinderSurfaceFlux(), and G4PSCylinderSurfaceFlux3D::G4PSCylinderSurfaceFlux3D().

◆ SetVerboseLevel()

void G4VPrimitiveScorer::SetVerboseLevel ( G4int  vl)
inlineinherited

Definition at line 109 of file G4VPrimitiveScorer.hh.

109{ verboseLevel = vl; }

References G4VPrimitiveScorer::verboseLevel.

◆ Weighted()

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

Definition at line 68 of file G4PSCylinderSurfaceFlux.hh.

68{ weighted = flg; }

References weighted.

Field Documentation

◆ detector

G4MultiFunctionalDetector* G4VPrimitiveScorer::detector
protectedinherited

◆ divideByArea

G4bool G4PSCylinderSurfaceFlux::divideByArea
private

Definition at line 95 of file G4PSCylinderSurfaceFlux.hh.

Referenced by DivideByArea(), ProcessHits(), and SetUnit().

◆ EvtMap

G4THitsMap<G4double>* G4PSCylinderSurfaceFlux::EvtMap
private

Definition at line 93 of file G4PSCylinderSurfaceFlux.hh.

Referenced by clear(), Initialize(), PrintAll(), and ProcessHits().

◆ fDirection

G4int G4PSCylinderSurfaceFlux::fDirection
private

Definition at line 92 of file G4PSCylinderSurfaceFlux.hh.

Referenced by ProcessHits().

◆ filter

G4VSDFilter* G4VPrimitiveScorer::filter
protectedinherited

◆ fNi

G4int G4VPrimitiveScorer::fNi
protectedinherited

◆ fNj

G4int G4VPrimitiveScorer::fNj
protectedinherited

Definition at line 131 of file G4VPrimitiveScorer.hh.

Referenced by G4PSCellCharge3D::G4PSCellCharge3D(), G4PSCellFlux3D::G4PSCellFlux3D(), G4PSCylinderSurfaceCurrent3D::G4PSCylinderSurfaceCurrent3D(), G4PSCylinderSurfaceFlux3D::G4PSCylinderSurfaceFlux3D(), G4PSDoseDeposit3D::G4PSDoseDeposit3D(), G4PSEnergyDeposit3D::G4PSEnergyDeposit3D(), G4PSFlatSurfaceCurrent3D::G4PSFlatSurfaceCurrent3D(), G4PSFlatSurfaceFlux3D::G4PSFlatSurfaceFlux3D(), G4PSMinKinEAtGeneration3D::G4PSMinKinEAtGeneration3D(), G4PSNofCollision3D::G4PSNofCollision3D(), G4PSNofSecondary3D::G4PSNofSecondary3D(), G4PSNofStep3D::G4PSNofStep3D(), G4PSPassageCellCurrent3D::G4PSPassageCellCurrent3D(), G4PSPassageCellFlux3D::G4PSPassageCellFlux3D(), G4PSPassageTrackLength3D::G4PSPassageTrackLength3D(), G4PSPopulation3D::G4PSPopulation3D(), G4PSSphereSurfaceCurrent3D::G4PSSphereSurfaceCurrent3D(), G4PSSphereSurfaceFlux3D::G4PSSphereSurfaceFlux3D(), G4PSStepChecker3D::G4PSStepChecker3D(), G4PSTermination3D::G4PSTermination3D(), G4PSTrackCounter3D::G4PSTrackCounter3D(), G4PSTrackLength3D::G4PSTrackLength3D(), G4PSVolumeFlux3D::G4PSVolumeFlux3D(), G4PSCellCharge3D::GetIndex(), G4PSCellFlux3D::GetIndex(), G4PSCylinderSurfaceCurrent3D::GetIndex(), G4PSCylinderSurfaceFlux3D::GetIndex(), G4PSDoseDeposit3D::GetIndex(), G4PSEnergyDeposit3D::GetIndex(), G4PSFlatSurfaceCurrent3D::GetIndex(), G4PSFlatSurfaceFlux3D::GetIndex(), G4PSMinKinEAtGeneration3D::GetIndex(), G4PSNofCollision3D::GetIndex(), G4PSNofSecondary3D::GetIndex(), G4PSNofStep3D::GetIndex(), G4PSPassageCellCurrent3D::GetIndex(), G4PSPassageCellFlux3D::GetIndex(), G4PSPassageTrackLength3D::GetIndex(), G4PSPopulation3D::GetIndex(), G4PSSphereSurfaceCurrent3D::GetIndex(), G4PSSphereSurfaceFlux3D::GetIndex(), G4PSStepChecker3D::GetIndex(), G4PSTermination3D::GetIndex(), G4PSTrackCounter3D::GetIndex(), G4PSTrackLength3D::GetIndex(), G4PSVolumeFlux3D::GetIndex(), and G4VPrimitiveScorer::SetNijk().

◆ fNk

G4int G4VPrimitiveScorer::fNk
protectedinherited

Definition at line 131 of file G4VPrimitiveScorer.hh.

Referenced by G4PSCellCharge3D::G4PSCellCharge3D(), G4PSCellFlux3D::G4PSCellFlux3D(), G4PSCylinderSurfaceCurrent3D::G4PSCylinderSurfaceCurrent3D(), G4PSCylinderSurfaceFlux3D::G4PSCylinderSurfaceFlux3D(), G4PSDoseDeposit3D::G4PSDoseDeposit3D(), G4PSEnergyDeposit3D::G4PSEnergyDeposit3D(), G4PSFlatSurfaceCurrent3D::G4PSFlatSurfaceCurrent3D(), G4PSFlatSurfaceFlux3D::G4PSFlatSurfaceFlux3D(), G4PSMinKinEAtGeneration3D::G4PSMinKinEAtGeneration3D(), G4PSNofCollision3D::G4PSNofCollision3D(), G4PSNofSecondary3D::G4PSNofSecondary3D(), G4PSNofStep3D::G4PSNofStep3D(), G4PSPassageCellCurrent3D::G4PSPassageCellCurrent3D(), G4PSPassageCellFlux3D::G4PSPassageCellFlux3D(), G4PSPassageTrackLength3D::G4PSPassageTrackLength3D(), G4PSPopulation3D::G4PSPopulation3D(), G4PSSphereSurfaceCurrent3D::G4PSSphereSurfaceCurrent3D(), G4PSSphereSurfaceFlux3D::G4PSSphereSurfaceFlux3D(), G4PSStepChecker3D::G4PSStepChecker3D(), G4PSTermination3D::G4PSTermination3D(), G4PSTrackCounter3D::G4PSTrackCounter3D(), G4PSTrackLength3D::G4PSTrackLength3D(), G4PSVolumeFlux3D::G4PSVolumeFlux3D(), G4PSCellCharge3D::GetIndex(), G4PSCellFlux3D::GetIndex(), G4PSCylinderSurfaceCurrent3D::GetIndex(), G4PSCylinderSurfaceFlux3D::GetIndex(), G4PSDoseDeposit3D::GetIndex(), G4PSEnergyDeposit3D::GetIndex(), G4PSFlatSurfaceCurrent3D::GetIndex(), G4PSFlatSurfaceFlux3D::GetIndex(), G4PSMinKinEAtGeneration3D::GetIndex(), G4PSNofCollision3D::GetIndex(), G4PSNofSecondary3D::GetIndex(), G4PSNofStep3D::GetIndex(), G4PSPassageCellCurrent3D::GetIndex(), G4PSPassageCellFlux3D::GetIndex(), G4PSPassageTrackLength3D::GetIndex(), G4PSPopulation3D::GetIndex(), G4PSSphereSurfaceCurrent3D::GetIndex(), G4PSSphereSurfaceFlux3D::GetIndex(), G4PSStepChecker3D::GetIndex(), G4PSTermination3D::GetIndex(), G4PSTrackCounter3D::GetIndex(), G4PSTrackLength3D::GetIndex(), G4PSVolumeFlux3D::GetIndex(), and G4VPrimitiveScorer::SetNijk().

◆ HCID

G4int G4PSCylinderSurfaceFlux::HCID
private

Definition at line 91 of file G4PSCylinderSurfaceFlux.hh.

Referenced by Initialize().

◆ hitIDMap

std::map<G4int, G4int> G4VPrimitivePlotter::hitIDMap
protectedinherited

◆ indexDepth

G4int G4VPrimitiveScorer::indexDepth
protectedinherited

◆ primitiveName

G4String G4VPrimitiveScorer::primitiveName
protectedinherited

◆ unitName

G4String G4VPrimitiveScorer::unitName
protectedinherited

◆ unitValue

G4double G4VPrimitiveScorer::unitValue
protectedinherited

◆ verboseLevel

G4int G4VPrimitiveScorer::verboseLevel
protectedinherited

◆ weighted

G4bool G4PSCylinderSurfaceFlux::weighted
private

Definition at line 94 of file G4PSCylinderSurfaceFlux.hh.

Referenced by ProcessHits(), and Weighted().


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