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

#include <G4PSFlatSurfaceCurrent.hh>

Inheritance diagram for G4PSFlatSurfaceCurrent:
Collaboration diagram for G4PSFlatSurfaceCurrent:

Public Member Functions

 G4PSFlatSurfaceCurrent (G4String name, G4int direction, G4int depth=0)
 
 G4PSFlatSurfaceCurrent (G4String name, G4int direction, const G4String &unit, G4int depth=0)
 
virtual ~G4PSFlatSurfaceCurrent ()
 
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 *, G4Box *)
 
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 G4PSFlatSurfaceCurrent.hh.

Constructor & Destructor Documentation

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

Definition at line 59 of file G4PSFlatSurfaceCurrent.cc.

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

Here is the call graph for this function:

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

Definition at line 68 of file G4PSFlatSurfaceCurrent.cc.

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

Here is the call graph for this function:

G4PSFlatSurfaceCurrent::~G4PSFlatSurfaceCurrent ( )
virtual

Definition at line 79 of file G4PSFlatSurfaceCurrent.cc.

80 {;}

Member Function Documentation

void G4PSFlatSurfaceCurrent::clear ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 159 of file G4PSFlatSurfaceCurrent.cc.

159  {
160  EvtMap->clear();
161 }
void clear()
Definition: G4THitsMap.hh:267

Here is the call graph for this function:

void G4PSFlatSurfaceCurrent::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 199 of file G4PSFlatSurfaceCurrent.cc.

199  {
200  // Per Unit Surface
201  new G4UnitDefinition("percentimeter2","percm2","Per Unit Surface",(1./cm2));
202  new G4UnitDefinition("permillimeter2","permm2","Per Unit Surface",(1./mm2));
203  new G4UnitDefinition("permeter2","perm2","Per Unit Surface",(1./m2));
204 }
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 G4PSFlatSurfaceCurrent::DivideByArea ( G4bool  flg = true)
inline

Definition at line 72 of file G4PSFlatSurfaceCurrent.hh.

72 { divideByArea = flg; }

Here is the caller graph for this function:

void G4PSFlatSurfaceCurrent::DrawAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 163 of file G4PSFlatSurfaceCurrent.cc.

164 {;}
void G4PSFlatSurfaceCurrent::EndOfEvent ( G4HCofThisEvent )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 156 of file G4PSFlatSurfaceCurrent.cc.

157 {;}
void G4PSFlatSurfaceCurrent::Initialize ( G4HCofThisEvent HCE)
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 149 of file G4PSFlatSurfaceCurrent.cc.

150 {
151  EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName());
152  if ( HCID < 0 ) HCID = GetCollectionID(0);
153  HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
154 }
G4String GetName() const
G4int GetCollectionID(G4int)
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)
G4MultiFunctionalDetector * detector

Here is the call graph for this function:

G4int G4PSFlatSurfaceCurrent::IsSelectedSurface ( G4Step aStep,
G4Box boxSolid 
)
protected

Definition at line 120 of file G4PSFlatSurfaceCurrent.cc.

120  {
121 
122  G4TouchableHandle theTouchable =
125 
126  if (aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary ){
127  // Entering Geometry
128  G4ThreeVector stppos1= aStep->GetPreStepPoint()->GetPosition();
129  G4ThreeVector localpos1 =
130  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos1);
131  if(std::fabs( localpos1.z() + boxSolid->GetZHalfLength())<kCarTolerance ){
132  return fCurrent_In;
133  }
134  }
135 
136  if (aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary ){
137  // Exiting Geometry
138  G4ThreeVector stppos2= aStep->GetPostStepPoint()->GetPosition();
139  G4ThreeVector localpos2 =
140  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos2);
141  if(std::fabs( localpos2.z() + boxSolid->GetZHalfLength())<kCarTolerance ){
142  return fCurrent_Out;
143  }
144  }
145 
146  return -1;
147 }
G4StepStatus GetStepStatus() const
G4double GetSurfaceTolerance() const
virtual const G4NavigationHistory * GetHistory() const
Definition: G4VTouchable.cc:86
G4double GetZHalfLength() const
double z() const
G4StepPoint * GetPreStepPoint() const
const G4ThreeVector & GetPosition() const
const G4double kCarTolerance
G4ThreeVector TransformPoint(const G4ThreeVector &vec) const
G4StepPoint * GetPostStepPoint() 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 G4PSFlatSurfaceCurrent::PrintAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 166 of file G4PSFlatSurfaceCurrent.cc.

167 {
168  G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
169  G4cout << " PrimitiveScorer " << GetName() <<G4endl;
170  G4cout << " Number of entries " << EvtMap->entries() << G4endl;
171  std::map<G4int,G4double*>::iterator itr = EvtMap->GetMap()->begin();
172  for(; itr != EvtMap->GetMap()->end(); itr++) {
173  G4cout << " copy no.: " << itr->first << " current : " ;
174  if ( divideByArea ) {
175  G4cout << *(itr->second)/GetUnitValue()
176  << " ["<<GetUnit()<<"]";
177  }else {
178  G4cout << *(itr->second)/GetUnitValue() << " [tracks]";
179  }
180  G4cout << G4endl;
181  }
182 }
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 G4PSFlatSurfaceCurrent::ProcessHits ( G4Step aStep,
G4TouchableHistory  
)
protectedvirtual

Implements G4VPrimitiveScorer.

Definition at line 82 of file G4PSFlatSurfaceCurrent.cc.

83 {
84  G4StepPoint* preStep = aStep->GetPreStepPoint();
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  G4Box* boxSolid = (G4Box*)(solid);
101 
102  G4int dirFlag =IsSelectedSurface(aStep,boxSolid);
103  if ( dirFlag > 0 ) {
104  if ( fDirection == fCurrent_InOut || fDirection == dirFlag ){
105  G4int index = GetIndex(aStep);
106  G4TouchableHandle theTouchable = preStep->GetTouchableHandle();
107  G4double current = 1.0;
108  if ( weighted ) current=preStep->GetWeight(); // Current (Particle Weight)
109  if ( divideByArea ){
110  G4double square = 4.*boxSolid->GetXHalfLength()*boxSolid->GetYHalfLength();
111  current = current/square; // Normalized by Area
112  }
113  EvtMap->add(index,current);
114  }
115  }
116 
117  return TRUE;
118 }
G4double GetXHalfLength() const
G4double GetWeight() const
virtual void ComputeDimensions(G4VPVParameterisation *p, const G4int n, const G4VPhysicalVolume *pRep)
Definition: G4VSolid.cc:138
Definition: G4Box.hh:64
virtual G4VSolid * ComputeSolid(const G4int, G4VPhysicalVolume *)
G4VSolid * GetSolid() const
const G4VTouchable * GetTouchable() const
int G4int
Definition: G4Types.hh:78
G4StepPoint * GetPreStepPoint() const
G4int IsSelectedSurface(G4Step *, G4Box *)
G4VPhysicalVolume * GetPhysicalVolume() const
virtual G4VPVParameterisation * GetParameterisation() const =0
virtual G4int GetIndex(G4Step *)
G4double GetYHalfLength() const
#define TRUE
Definition: globals.hh:55
G4int add(const G4int &key, T *&aHit) const
Definition: G4THitsMap.hh:106
G4LogicalVolume * GetLogicalVolume() const
double G4double
Definition: G4Types.hh:76
const G4TouchableHandle & GetTouchableHandle() const

Here is the call graph for this function:

void G4PSFlatSurfaceCurrent::SetUnit ( const G4String unit)
virtual

Definition at line 184 of file G4PSFlatSurfaceCurrent.cc.

185 {
186  if ( divideByArea ) {
187  CheckAndSetUnit(unit,"Per Unit Surface");
188  } else {
189  if (unit == "" ){
190  unitName = unit;
191  unitValue = 1.0;
192  }else{
193  G4String msg = "Invalid unit ["+unit+"] (Current unit is [" +GetUnit()+"] ) for " + GetName();
194  G4Exception("G4PSFlatSurfaceCurrent::SetUnit","DetPS0007",JustWarning,msg);
195  }
196  }
197 }
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 G4PSFlatSurfaceCurrent::Weighted ( G4bool  flg = true)
inline

Definition at line 70 of file G4PSFlatSurfaceCurrent.hh.

70 { weighted = flg; }

Here is the caller graph for this function:


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