Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PSSphereSurfaceFlux Class Reference

#include <G4PSSphereSurfaceFlux.hh>

Inheritance diagram for G4PSSphereSurfaceFlux:
Collaboration diagram for G4PSSphereSurfaceFlux:

Public Member Functions

 G4PSSphereSurfaceFlux (G4String name, G4int direction, G4int depth=0)
 
 G4PSSphereSurfaceFlux (G4String name, G4int direction, const G4String &unit, G4int depth=0)
 
virtual ~G4PSSphereSurfaceFlux ()
 
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 *, G4Sphere *)
 
virtual void DefineUnitAndCategory ()
 
- Protected Member Functions inherited from G4VPrimitiveScorer
virtual G4int GetIndex (G4Step *)
 
void CheckAndSetUnit (const G4String &unit, const G4String &category)
 

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 60 of file G4PSSphereSurfaceFlux.hh.

Constructor & Destructor Documentation

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

Definition at line 60 of file G4PSSphereSurfaceFlux.cc.

62  : G4VPrimitiveScorer(name,depth),HCID(-1),fDirection(direction),
63  EvtMap(0),weighted(true),divideByArea(true)
64 {
66  SetUnit("percm2");
67 }
virtual void DefineUnitAndCategory()
virtual void SetUnit(const G4String &unit)
G4VPrimitiveScorer(G4String name, G4int depth=0)

Here is the call graph for this function:

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

Definition at line 69 of file G4PSSphereSurfaceFlux.cc.

73  : G4VPrimitiveScorer(name,depth),HCID(-1),fDirection(direction),
74  EvtMap(0),weighted(true),divideByArea(true)
75 {
77  SetUnit(unit);
78 }
virtual void DefineUnitAndCategory()
virtual void SetUnit(const G4String &unit)
G4VPrimitiveScorer(G4String name, G4int depth=0)

Here is the call graph for this function:

G4PSSphereSurfaceFlux::~G4PSSphereSurfaceFlux ( )
virtual

Definition at line 80 of file G4PSSphereSurfaceFlux.cc.

81 {;}

Member Function Documentation

void G4PSSphereSurfaceFlux::clear ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 211 of file G4PSSphereSurfaceFlux.cc.

211  {
212  EvtMap->clear();
213 }
void clear()
Definition: G4THitsMap.hh:267

Here is the call graph for this function:

void G4PSSphereSurfaceFlux::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 247 of file G4PSSphereSurfaceFlux.cc.

247  {
248  // Per Unit Surface
249  new G4UnitDefinition("percentimeter2","percm2","Per Unit Surface",(1./cm2));
250  new G4UnitDefinition("permillimeter2","permm2","Per Unit Surface",(1./mm2));
251  new G4UnitDefinition("permeter2","perm2","Per Unit Surface",(1./m2));
252 }
static constexpr double cm2
Definition: G4SIunits.hh:120
static constexpr double mm2
Definition: G4SIunits.hh:116
static constexpr double m2
Definition: G4SIunits.hh:130

Here is the caller graph for this function:

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

Definition at line 72 of file G4PSSphereSurfaceFlux.hh.

72 { divideByArea = flg; }
void G4PSSphereSurfaceFlux::DrawAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 215 of file G4PSSphereSurfaceFlux.cc.

216 {;}
void G4PSSphereSurfaceFlux::EndOfEvent ( G4HCofThisEvent )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 208 of file G4PSSphereSurfaceFlux.cc.

209 {;}
void G4PSSphereSurfaceFlux::Initialize ( G4HCofThisEvent HCE)
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 201 of file G4PSSphereSurfaceFlux.cc.

202 {
203  EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName());
204  if ( HCID < 0 ) HCID = GetCollectionID(0);
205  HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
206 }
G4String GetName() const
G4int GetCollectionID(G4int)
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)
G4MultiFunctionalDetector * detector

Here is the call graph for this function:

G4int G4PSSphereSurfaceFlux::IsSelectedSurface ( G4Step aStep,
G4Sphere sphereSolid 
)
protected

Definition at line 156 of file G4PSSphereSurfaceFlux.cc.

156  {
157 
158  G4TouchableHandle theTouchable =
161 
162  if (aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary ){
163  // Entering Geometry
164  G4ThreeVector stppos1= aStep->GetPreStepPoint()->GetPosition();
165  G4ThreeVector localpos1 =
166  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos1);
167  G4double localR2 = localpos1.x()*localpos1.x()
168  +localpos1.y()*localpos1.y()
169  +localpos1.z()*localpos1.z();
170  //G4double InsideRadius2 =
171  // sphereSolid->GetInsideRadius()*sphereSolid->GetInsideRadius();
172  //if(std::fabs( localR2 - InsideRadius2 ) < kCarTolerance ){
173  G4double InsideRadius = sphereSolid->GetInnerRadius();
174  if ( localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
175  &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
176  return fFlux_In;
177  }
178  }
179 
180  if (aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary ){
181  // Exiting Geometry
182  G4ThreeVector stppos2= aStep->GetPostStepPoint()->GetPosition();
183  G4ThreeVector localpos2 =
184  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos2);
185  G4double localR2 = localpos2.x()*localpos2.x()
186  +localpos2.y()*localpos2.y()
187  +localpos2.z()*localpos2.z();
188  //G4double InsideRadius2 =
189  // sphereSolid->GetInsideRadius()*sphereSolid->GetInsideRadius();
190  //if(std::facb(localR2 - InsideRadius2) ) < kCarTolerance ){
191  G4double InsideRadius = sphereSolid->GetInnerRadius();
192  if ( localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
193  &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
194  return fFlux_Out;
195  }
196  }
197 
198  return -1;
199 }
double x() const
G4StepStatus GetStepStatus() const
G4double GetSurfaceTolerance() const
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:86
double z() const
G4StepPoint * GetPreStepPoint() const
const G4ThreeVector & GetPosition() const
const G4double kCarTolerance
G4double GetInnerRadius() const
G4ThreeVector TransformPoint(const G4ThreeVector &vec) const
G4StepPoint * GetPostStepPoint() const
double y() const
const G4AffineTransform & GetTopTransform() const
double G4double
Definition: G4Types.hh:76
const G4TouchableHandle & GetTouchableHandle() const
static G4GeometryTolerance * GetInstance()

Here is the call graph for this function:

Here is the caller graph for this function:

void G4PSSphereSurfaceFlux::PrintAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 218 of file G4PSSphereSurfaceFlux.cc.

219 {
220  G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
221  G4cout << " PrimitiveScorer " << GetName() <<G4endl;
222  G4cout << " Number of entries " << EvtMap->entries() << G4endl;
223  std::map<G4int,G4double*>::iterator itr = EvtMap->GetMap()->begin();
224  for(; itr != EvtMap->GetMap()->end(); itr++) {
225  G4cout << " copy no.: " << itr->first
226  << " Flux : " << *(itr->second)/GetUnitValue()
227  << " ["<<GetUnit()<<"]"
228  << G4endl;
229  }
230 }
G4String GetName() const
G4GLOB_DLL std::ostream G4cout
G4double GetUnitValue() const
G4int entries() const
Definition: G4THitsMap.hh:200
std::map< G4int, T * > * GetMap() const
Definition: G4THitsMap.hh:99
#define G4endl
Definition: G4ios.hh:61
G4MultiFunctionalDetector * detector
const G4String & GetUnit() const

Here is the call graph for this function:

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

Implements G4VPrimitiveScorer.

Definition at line 83 of file G4PSSphereSurfaceFlux.cc.

84 {
85  G4StepPoint* preStep = aStep->GetPreStepPoint();
86 
87  G4VPhysicalVolume* physVol = preStep->GetPhysicalVolume();
88  G4VPVParameterisation* physParam = physVol->GetParameterisation();
89  G4VSolid * solid = 0;
90  if(physParam)
91  { // for parameterized volume
93  ->GetReplicaNumber(indexDepth);
94  solid = physParam->ComputeSolid(idx, physVol);
95  solid->ComputeDimensions(physParam,idx,physVol);
96  }
97  else
98  { // for ordinary volume
99  solid = physVol->GetLogicalVolume()->GetSolid();
100  }
101 
102  G4Sphere* sphereSolid = (G4Sphere*)(solid);
103 
104  G4int dirFlag =IsSelectedSurface(aStep,sphereSolid);
105  if ( dirFlag > 0 ) {
106  if ( fDirection == fFlux_InOut || fDirection == dirFlag ){
107 
108  G4StepPoint* thisStep=0;
109  if ( dirFlag == fFlux_In ){
110  thisStep = preStep;
111  }else if ( dirFlag == fFlux_Out ){
112  thisStep = aStep->GetPostStepPoint();
113  }else{
114  return FALSE;
115  }
116 
117  G4TouchableHandle theTouchable = thisStep->GetTouchableHandle();
118  G4ThreeVector pdirection = thisStep->GetMomentumDirection();
119  G4ThreeVector localdir =
120  theTouchable->GetHistory()->GetTopTransform().TransformAxis(pdirection);
121  G4double localdirL2 = localdir.x()*localdir.x()
122  +localdir.y()*localdir.y()
123  +localdir.z()*localdir.z();
124  G4ThreeVector stppos1= aStep->GetPreStepPoint()->GetPosition();
125  G4ThreeVector localpos1 =
126  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos1);
127  G4double localR2 = localpos1.x()*localpos1.x()
128  +localpos1.y()*localpos1.y()
129  +localpos1.z()*localpos1.z();
130  G4double anglefactor = (localdir.x()*localpos1.x()
131  +localdir.y()*localpos1.y()
132  +localdir.z()*localpos1.z())
133  /std::sqrt(localdirL2)/std::sqrt(localR2);
134  if ( anglefactor < 0.0 ) anglefactor *= -1.0;
135 
136  G4double radi = sphereSolid->GetInnerRadius();
137  G4double dph = sphereSolid->GetDeltaPhiAngle()/radian;
138  G4double stth = sphereSolid->GetStartThetaAngle()/radian;
139  G4double enth = stth+sphereSolid->GetDeltaThetaAngle()/radian;
140  G4double square = radi*radi*dph*( -std::cos(enth) + std::cos(stth) );
141 
142  G4double current = 1.0;
143  if ( weighted ) thisStep->GetWeight(); // Flux (Particle Weight)
144  if ( divideByArea ) current = current/square; // Flux with angle.
145 
146  current /= anglefactor;
147 
148  G4int index = GetIndex(aStep);
149  EvtMap->add(index,current);
150  }
151  }
152 
153  return TRUE;
154 }
G4double GetWeight() const
double x() const
virtual void ComputeDimensions(G4VPVParameterisation *p, const G4int n, const G4VPhysicalVolume *pRep)
Definition: G4VSolid.cc:138
virtual G4VSolid * ComputeSolid(const G4int, G4VPhysicalVolume *)
G4VSolid * GetSolid() const
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:86
const G4VTouchable * GetTouchable() const
G4double GetDeltaPhiAngle() const
int G4int
Definition: G4Types.hh:78
double z() const
G4StepPoint * GetPreStepPoint() const
const G4ThreeVector & GetMomentumDirection() const
G4VPhysicalVolume * GetPhysicalVolume() const
G4double GetStartThetaAngle() const
const G4ThreeVector & GetPosition() const
#define FALSE
Definition: globals.hh:52
virtual G4VPVParameterisation * GetParameterisation() const =0
virtual G4int GetIndex(G4Step *)
G4int IsSelectedSurface(G4Step *, G4Sphere *)
#define TRUE
Definition: globals.hh:55
G4double GetInnerRadius() const
static constexpr double radian
Definition: G4SIunits.hh:142
G4ThreeVector TransformPoint(const G4ThreeVector &vec) const
G4int add(const G4int &key, T *&aHit) const
Definition: G4THitsMap.hh:106
G4LogicalVolume * GetLogicalVolume() const
G4StepPoint * GetPostStepPoint() const
G4ThreeVector TransformAxis(const G4ThreeVector &axis) const
double y() const
const G4AffineTransform & GetTopTransform() const
double G4double
Definition: G4Types.hh:76
G4double GetDeltaThetaAngle() const
const G4TouchableHandle & GetTouchableHandle() const

Here is the call graph for this function:

void G4PSSphereSurfaceFlux::SetUnit ( const G4String unit)
virtual

Definition at line 232 of file G4PSSphereSurfaceFlux.cc.

233 {
234  if ( divideByArea ) {
235  CheckAndSetUnit(unit,"Per Unit Surface");
236  } else {
237  if (unit == "" ){
238  unitName = unit;
239  unitValue = 1.0;
240  }else{
241  G4String msg = "Invalid unit ["+unit+"] (Current unit is [" +GetUnit()+"] ) for " + GetName();
242  G4Exception("G4PSSphereSurfaceFlux::SetUnit","DetPS0016",JustWarning,msg);
243  }
244  }
245 }
G4String GetName() const
void CheckAndSetUnit(const G4String &unit, const G4String &category)
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:

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

Definition at line 69 of file G4PSSphereSurfaceFlux.hh.

69 { weighted = flg; }

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