Geant4  10.02.p03
G4PSCylinderSurfaceFlux Class Reference

#include <G4PSCylinderSurfaceFlux.hh>

Inheritance diagram for G4PSCylinderSurfaceFlux:
Collaboration diagram for G4PSCylinderSurfaceFlux:

Public Member Functions

 G4PSCylinderSurfaceFlux (G4String name, G4int direction, G4int depth=0)
 
 G4PSCylinderSurfaceFlux (G4String name, G4int direction, const G4String &unit, G4int depth=0)
 
virtual ~G4PSCylinderSurfaceFlux ()
 
void Weighted (G4bool flg=true)
 
void DivideByArea (G4bool flg=true)
 
virtual void Initialize (G4HCofThisEvent *)
 
virtual void EndOfEvent (G4HCofThisEvent *)
 
virtual void clear ()
 
virtual void DrawAll ()
 
virtual void PrintAll ()
 
virtual void SetUnit (const G4String &unit)
 
- Public Member Functions inherited from G4VPrimitiveScorer
 G4VPrimitiveScorer (G4String name, G4int depth=0)
 
virtual ~G4VPrimitiveScorer ()
 
G4int GetCollectionID (G4int)
 
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

virtual G4bool ProcessHits (G4Step *, G4TouchableHistory *)
 
G4int IsSelectedSurface (G4Step *, G4Tubs *)
 
virtual void DefineUnitAndCategory ()
 
- Protected Member Functions inherited from G4VPrimitiveScorer
virtual G4int GetIndex (G4Step *)
 
void CheckAndSetUnit (const G4String &unit, const G4String &category)
 

Private Attributes

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

Additional Inherited Members

- 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 59 of file G4PSCylinderSurfaceFlux.hh.

Constructor & Destructor Documentation

◆ G4PSCylinderSurfaceFlux() [1/2]

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

Definition at line 58 of file G4PSCylinderSurfaceFlux.cc.

60  : G4VPrimitiveScorer(name,depth),HCID(-1),fDirection(direction),EvtMap(0),
61  weighted(true),divideByArea(true)
62 {
64  SetUnit("percm2");
65 }
G4THitsMap< G4double > * EvtMap
virtual void SetUnit(const G4String &unit)
G4VPrimitiveScorer(G4String name, G4int depth=0)
Here is the call graph for this function:

◆ G4PSCylinderSurfaceFlux() [2/2]

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

Definition at line 67 of file G4PSCylinderSurfaceFlux.cc.

71  : G4VPrimitiveScorer(name,depth),HCID(-1),fDirection(direction),EvtMap(0),
72  weighted(true),divideByArea(true)
73 {
75  SetUnit(unit);
76 }
G4THitsMap< G4double > * EvtMap
virtual void SetUnit(const G4String &unit)
G4VPrimitiveScorer(G4String name, G4int depth=0)
Here is the call graph for this function:

◆ ~G4PSCylinderSurfaceFlux()

G4PSCylinderSurfaceFlux::~G4PSCylinderSurfaceFlux ( )
virtual

Definition at line 78 of file G4PSCylinderSurfaceFlux.cc.

79 {;}

Member Function Documentation

◆ clear()

void G4PSCylinderSurfaceFlux::clear ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 204 of file G4PSCylinderSurfaceFlux.cc.

204  {
205  EvtMap->clear();
206 }
void clear()
Definition: G4THitsMap.hh:209
G4THitsMap< G4double > * EvtMap
Here is the call graph for this function:
Here is the caller graph for this function:

◆ DefineUnitAndCategory()

void G4PSCylinderSurfaceFlux::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 240 of file G4PSCylinderSurfaceFlux.cc.

240  {
241  // Per Unit Surface
242  new G4UnitDefinition("percentimeter2","percm2","Per Unit Surface",(1./cm2));
243  new G4UnitDefinition("permillimeter2","permm2","Per Unit Surface",(1./mm2));
244  new G4UnitDefinition("permeter2","perm2","Per Unit Surface",(1./m2));
245 }
static const double cm2
Definition: G4SIunits.hh:119
static const double m2
Definition: G4SIunits.hh:129
static const double mm2
Definition: G4SIunits.hh:115
Here is the caller graph for this function:

◆ DivideByArea()

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

Definition at line 70 of file G4PSCylinderSurfaceFlux.hh.

Here is the call graph for this function:

◆ DrawAll()

void G4PSCylinderSurfaceFlux::DrawAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 208 of file G4PSCylinderSurfaceFlux.cc.

209 {;}
Here is the caller graph for this function:

◆ EndOfEvent()

void G4PSCylinderSurfaceFlux::EndOfEvent ( G4HCofThisEvent )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 201 of file G4PSCylinderSurfaceFlux.cc.

202 {;}
Here is the caller graph for this function:

◆ Initialize()

void G4PSCylinderSurfaceFlux::Initialize ( G4HCofThisEvent HCE)
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 193 of file G4PSCylinderSurfaceFlux.cc.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ IsSelectedSurface()

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

Definition at line 150 of file G4PSCylinderSurfaceFlux.cc.

150  {
151 
152  G4TouchableHandle theTouchable =
153  aStep->GetPreStepPoint()->GetTouchableHandle();
155 
156  if (aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary ){
157  // Entering Geometry
158  G4ThreeVector stppos1= aStep->GetPreStepPoint()->GetPosition();
159  G4ThreeVector localpos1 =
160  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos1);
161  if ( std::fabs(localpos1.z()) > tubsSolid->GetZHalfLength() ) return -1;
162  //if(std::fabs( localpos1.x()*localpos1.x()+localpos1.y()*localpos1.y()
163  // - (tubsSolid->GetInnerRadius()*tubsSolid->GetInnerRadius()))
164  // <kCarTolerance ){
165  G4double localR2 = localpos1.x()*localpos1.x()+localpos1.y()*localpos1.y();
166  G4double InsideRadius = tubsSolid->GetInnerRadius();
167  if (localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
168  &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
169  return fFlux_In;
170  }
171  }
172 
173  if (aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary ){
174  // Exiting Geometry
175  G4ThreeVector stppos2= aStep->GetPostStepPoint()->GetPosition();
176  G4ThreeVector localpos2 =
177  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos2);
178  if ( std::fabs(localpos2.z()) > tubsSolid->GetZHalfLength() ) return -1;
179  //if(std::fabs( localpos2.x()*localpos2.x()+localpos2.y()*localpos2.y()
180  // - (tubsSolid->GetInnerRadius()*tubsSolid->GetInnerRadius()))
181  // <kCarTolerance ){
182  G4double localR2 = localpos2.x()*localpos2.x()+localpos2.y()*localpos2.y();
183  G4double InsideRadius = tubsSolid->GetInnerRadius();
184  if (localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
185  &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
186  return fFlux_Out;
187  }
188  }
189 
190  return -1;
191 }
G4ThreeVector TransformPoint(const G4ThreeVector &vec) const
G4double GetSurfaceTolerance() const
G4double GetZHalfLength() const
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:86
double x() const
double y() const
double z() const
const G4AffineTransform & GetTopTransform() const
G4double GetInnerRadius() const
double G4double
Definition: G4Types.hh:76
static G4GeometryTolerance * GetInstance()
Here is the call graph for this function:
Here is the caller graph for this function:

◆ PrintAll()

void G4PSCylinderSurfaceFlux::PrintAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 211 of file G4PSCylinderSurfaceFlux.cc.

212 {
213  G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
214  G4cout << " PrimitiveScorer" << GetName() <<G4endl;
215  G4cout << " Number of entries " << EvtMap->entries() << G4endl;
216  std::map<G4int,G4double*>::iterator itr = EvtMap->GetMap()->begin();
217  for(; itr != EvtMap->GetMap()->end(); itr++) {
218  G4cout << " copy no.: " << itr->first
219  << " flux : " << *(itr->second)/GetUnitValue()
220  << " ["<<GetUnit()<<"]"
221  << G4endl;
222  }
223 }
G4double GetUnitValue() const
std::map< G4int, T * > * GetMap() const
Definition: G4THitsMap.hh:68
G4THitsMap< G4double > * EvtMap
G4String GetName() const
G4GLOB_DLL std::ostream G4cout
G4int entries() const
Definition: G4THitsMap.hh:79
#define G4endl
Definition: G4ios.hh:61
G4MultiFunctionalDetector * detector
const G4String & GetUnit() const
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ProcessHits()

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

Implements G4VPrimitiveScorer.

Definition at line 81 of file G4PSCylinderSurfaceFlux.cc.

82 {
83  G4StepPoint* preStep = aStep->GetPreStepPoint();
84 
85  G4VPhysicalVolume* physVol = preStep->GetPhysicalVolume();
86  G4VPVParameterisation* physParam = physVol->GetParameterisation();
87  G4VSolid * solid = 0;
88  if(physParam)
89  { // for parameterized volume
90  G4int idx = ((G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable()))
91  ->GetReplicaNumber(indexDepth);
92  solid = physParam->ComputeSolid(idx, physVol);
93  solid->ComputeDimensions(physParam,idx,physVol);
94  }
95  else
96  { // for ordinary volume
97  solid = physVol->GetLogicalVolume()->GetSolid();
98  }
99 
100  G4Tubs* tubsSolid = (G4Tubs*)(solid);
101 
102  G4int dirFlag =IsSelectedSurface(aStep,tubsSolid);
103 
104  if ( dirFlag > 0 ){
105  if (fDirection == fFlux_InOut || dirFlag == fDirection ){
106 
107  G4StepPoint* thisStep=0;
108  if ( dirFlag == fFlux_In ){
109  thisStep = preStep;
110  }else if ( dirFlag == fFlux_Out ){
111  thisStep = aStep->GetPostStepPoint();
112  }else{
113  return FALSE;
114  }
115 
116  G4TouchableHandle theTouchable = thisStep->GetTouchableHandle();
117  G4ThreeVector pdirection = thisStep->GetMomentumDirection();
118  G4ThreeVector localdir =
119  theTouchable->GetHistory()->GetTopTransform().TransformAxis(pdirection);
120  G4ThreeVector position = thisStep->GetPosition();
121  G4ThreeVector localpos =
122  theTouchable->GetHistory()->GetTopTransform().TransformAxis(position);
123  G4double angleFactor = (localdir.x()*localpos.x()+localdir.y()*localpos.y())
124  /std::sqrt(localdir.x()*localdir.x()
125  +localdir.y()*localdir.y()+localdir.z()*localdir.z())
126  /std::sqrt(localpos.x()*localpos.x()+localpos.y()*localpos.y());
127 
128  if ( angleFactor < 0 ) angleFactor *= -1.;
129  G4double square = 2.*tubsSolid->GetZHalfLength()
130  *tubsSolid->GetInnerRadius()* tubsSolid->GetDeltaPhiAngle()/radian;
131 
132  G4double flux = 1.0;
133  if ( weighted ) flux *=preStep->GetWeight();
134  // Current (Particle Weight)
135 
136  flux = flux/angleFactor;
137  if ( divideByArea ) flux /= square;
138  //Flux with angle.
139  G4int index = GetIndex(aStep);
140  EvtMap->add(index,flux);
141  return TRUE;
142  }else{
143  return FALSE;
144  }
145  }else{
146  return FALSE;
147  }
148 }
G4ThreeVector TransformAxis(const G4ThreeVector &axis) const
Int_t index
G4int add(const G4int &key, T *&aHit) const
Definition: G4THitsMap.hh:138
virtual void ComputeDimensions(G4VPVParameterisation *p, const G4int n, const G4VPhysicalVolume *pRep)
Definition: G4VSolid.cc:137
virtual G4VSolid * ComputeSolid(const G4int, G4VPhysicalVolume *)
Definition: G4Tubs.hh:85
G4THitsMap< G4double > * EvtMap
int G4int
Definition: G4Types.hh:78
G4double GetZHalfLength() const
#define FALSE
Definition: globals.hh:52
virtual G4VPVParameterisation * GetParameterisation() const =0
virtual G4int GetIndex(G4Step *)
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:86
#define TRUE
Definition: globals.hh:55
double x() const
double y() const
double z() const
const G4AffineTransform & GetTopTransform() const
static const double radian
Definition: G4SIunits.hh:141
G4double GetInnerRadius() const
G4int IsSelectedSurface(G4Step *, G4Tubs *)
double G4double
Definition: G4Types.hh:76
G4LogicalVolume * GetLogicalVolume() const
G4VSolid * GetSolid() const
G4double GetDeltaPhiAngle() const
Here is the call graph for this function:
Here is the caller graph for this function:

◆ SetUnit()

void G4PSCylinderSurfaceFlux::SetUnit ( const G4String unit)
virtual

Definition at line 225 of file G4PSCylinderSurfaceFlux.cc.

226 {
227  if ( divideByArea ) {
228  CheckAndSetUnit(unit,"Per Unit Surface");
229  } else {
230  if (unit == "" ){
231  unitName = unit;
232  unitValue = 1.0;
233  }else{
234  G4String msg = "Invalid unit ["+unit+"] (Current unit is [" +GetUnit()+"] ) for " + GetName();
235  G4Exception("G4PSCylinderSurfaceFlux::SetUnit","DetPS0003",JustWarning,msg);
236  }
237  }
238 }
void CheckAndSetUnit(const G4String &unit, const G4String &category)
G4String GetName() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
const G4String & GetUnit() const
Here is the call graph for this function:
Here is the caller graph for this function:

◆ Weighted()

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

Definition at line 67 of file G4PSCylinderSurfaceFlux.hh.

Member Data Documentation

◆ divideByArea

G4bool G4PSCylinderSurfaceFlux::divideByArea
private

Definition at line 95 of file G4PSCylinderSurfaceFlux.hh.

◆ EvtMap

G4THitsMap<G4double>* G4PSCylinderSurfaceFlux::EvtMap
private

Definition at line 93 of file G4PSCylinderSurfaceFlux.hh.

◆ fDirection

G4int G4PSCylinderSurfaceFlux::fDirection
private

Definition at line 92 of file G4PSCylinderSurfaceFlux.hh.

◆ HCID

G4int G4PSCylinderSurfaceFlux::HCID
private

Definition at line 91 of file G4PSCylinderSurfaceFlux.hh.

◆ weighted

G4bool G4PSCylinderSurfaceFlux::weighted
private

Definition at line 94 of file G4PSCylinderSurfaceFlux.hh.


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