Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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)
 

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::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 }
virtual void SetUnit(const G4String &unit)
G4VPrimitiveScorer(G4String name, G4int depth=0)

Here is the call graph for this function:

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 }
virtual void SetUnit(const G4String &unit)
G4VPrimitiveScorer(G4String name, G4int depth=0)

Here is the call graph for this function:

G4PSCylinderSurfaceFlux::~G4PSCylinderSurfaceFlux ( )
virtual

Definition at line 78 of file G4PSCylinderSurfaceFlux.cc.

79 {;}

Member Function Documentation

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

Here is the call graph for this function:

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 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 G4PSCylinderSurfaceFlux::DivideByArea ( G4bool  flg = true)
inline

Definition at line 70 of file G4PSCylinderSurfaceFlux.hh.

70 { divideByArea = flg; }
void G4PSCylinderSurfaceFlux::DrawAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 208 of file G4PSCylinderSurfaceFlux.cc.

209 {;}
void G4PSCylinderSurfaceFlux::EndOfEvent ( G4HCofThisEvent )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 201 of file G4PSCylinderSurfaceFlux.cc.

202 {;}
void G4PSCylinderSurfaceFlux::Initialize ( G4HCofThisEvent HCE)
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 193 of file G4PSCylinderSurfaceFlux.cc.

194 {
196  GetName());
197  if ( HCID < 0 ) HCID = GetCollectionID(0);
198  HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
199 }
G4String GetName() const
G4MultiFunctionalDetector * GetMultiFunctionalDetector() const
G4int GetCollectionID(G4int)
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)

Here is the call graph for this function:

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

Definition at line 150 of file G4PSCylinderSurfaceFlux.cc.

150  {
151 
152  G4TouchableHandle theTouchable =
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 }
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
G4double GetZHalfLength() 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 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 }
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 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
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 }
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 *)
Definition: G4Tubs.hh:85
G4VSolid * GetSolid() const
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:86
const G4VTouchable * GetTouchable() const
int G4int
Definition: G4Types.hh:78
double z() const
G4StepPoint * GetPreStepPoint() const
const G4ThreeVector & GetMomentumDirection() const
G4double GetDeltaPhiAngle() const
G4VPhysicalVolume * GetPhysicalVolume() const
const G4ThreeVector & GetPosition() const
#define FALSE
Definition: globals.hh:52
virtual G4VPVParameterisation * GetParameterisation() const =0
virtual G4int GetIndex(G4Step *)
#define TRUE
Definition: globals.hh:55
G4double GetInnerRadius() const
static constexpr double radian
Definition: G4SIunits.hh:142
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
G4double GetZHalfLength() const
const G4AffineTransform & GetTopTransform() const
G4int IsSelectedSurface(G4Step *, G4Tubs *)
double G4double
Definition: G4Types.hh:76
const G4TouchableHandle & GetTouchableHandle() const

Here is the call graph for this function:

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 }
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 G4PSCylinderSurfaceFlux::Weighted ( G4bool  flg = true)
inline

Definition at line 67 of file G4PSCylinderSurfaceFlux.hh.

67 { weighted = flg; }

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