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

#include <G4PSCylinderSurfaceCurrent.hh>

Inheritance diagram for G4PSCylinderSurfaceCurrent:
Collaboration diagram for G4PSCylinderSurfaceCurrent:

Public Member Functions

 G4PSCylinderSurfaceCurrent (G4String name, G4int direction, G4int depth=0)
 
 G4PSCylinderSurfaceCurrent (G4String name, G4int direction, const G4String &unit, G4int depth=0)
 
virtual ~G4PSCylinderSurfaceCurrent ()
 
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 58 of file G4PSCylinderSurfaceCurrent.hh.

Constructor & Destructor Documentation

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

Definition at line 57 of file G4PSCylinderSurfaceCurrent.cc.

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

Here is the call graph for this function:

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

Definition at line 66 of file G4PSCylinderSurfaceCurrent.cc.

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

Here is the call graph for this function:

G4PSCylinderSurfaceCurrent::~G4PSCylinderSurfaceCurrent ( )
virtual

Definition at line 77 of file G4PSCylinderSurfaceCurrent.cc.

78 {;}

Member Function Documentation

void G4PSCylinderSurfaceCurrent::clear ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 171 of file G4PSCylinderSurfaceCurrent.cc.

171  {
172  EvtMap->clear();
173 }
void clear()
Definition: G4THitsMap.hh:267

Here is the call graph for this function:

void G4PSCylinderSurfaceCurrent::DefineUnitAndCategory ( )
protectedvirtual

Definition at line 213 of file G4PSCylinderSurfaceCurrent.cc.

213  {
214  // Per Unit Surface
215  new G4UnitDefinition("percentimeter2","percm2","Per Unit Surface",(1./cm2));
216  new G4UnitDefinition("permillimeter2","permm2","Per Unit Surface",(1./mm2));
217  new G4UnitDefinition("permeter2","perm2","Per Unit Surface",(1./m2));
218 }
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 G4PSCylinderSurfaceCurrent::DivideByArea ( G4bool  flg = true)
inline

Definition at line 71 of file G4PSCylinderSurfaceCurrent.hh.

71 { divideByArea = flg; }
void G4PSCylinderSurfaceCurrent::DrawAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 175 of file G4PSCylinderSurfaceCurrent.cc.

176 {;}
void G4PSCylinderSurfaceCurrent::EndOfEvent ( G4HCofThisEvent )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 168 of file G4PSCylinderSurfaceCurrent.cc.

169 {;}
void G4PSCylinderSurfaceCurrent::Initialize ( G4HCofThisEvent HCE)
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 161 of file G4PSCylinderSurfaceCurrent.cc.

162 {
163  EvtMap = new G4THitsMap<G4double>(detector->GetName(), GetName());
164  if ( HCID < 0 ) HCID = GetCollectionID(0);
165  HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
166 }
G4String GetName() const
G4int GetCollectionID(G4int)
void AddHitsCollection(G4int HCID, G4VHitsCollection *aHC)
G4MultiFunctionalDetector * detector

Here is the call graph for this function:

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

Definition at line 124 of file G4PSCylinderSurfaceCurrent.cc.

124  {
125 
126  G4TouchableHandle theTouchable =
129 
130  if (aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary ){
131  // Entering Geometry
132  G4ThreeVector stppos1= aStep->GetPreStepPoint()->GetPosition();
133  G4ThreeVector localpos1 =
134  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos1);
135  if ( std::fabs(localpos1.z()) > tubsSolid->GetZHalfLength() ) return -1;
136  G4double localR2 = localpos1.x()*localpos1.x()+localpos1.y()*localpos1.y();
137  G4double InsideRadius = tubsSolid->GetInnerRadius();
138  if (localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
139  &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
140  return fCurrent_In;
141  }
142  }
143 
144  if (aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary ){
145  // Exiting Geometry
146  G4ThreeVector stppos2= aStep->GetPostStepPoint()->GetPosition();
147  G4ThreeVector localpos2 =
148  theTouchable->GetHistory()->GetTopTransform().TransformPoint(stppos2);
149  if ( std::fabs(localpos2.z()) > tubsSolid->GetZHalfLength() ) return -1;
150  G4double localR2 = localpos2.x()*localpos2.x()+localpos2.y()*localpos2.y();
151  G4double InsideRadius = tubsSolid->GetInnerRadius();
152  if (localR2 > (InsideRadius-kCarTolerance)*(InsideRadius-kCarTolerance)
153  &&localR2 < (InsideRadius+kCarTolerance)*(InsideRadius+kCarTolerance)){
154  return fCurrent_Out;
155  }
156  }
157 
158  return -1;
159 }
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 G4PSCylinderSurfaceCurrent::PrintAll ( void  )
virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 178 of file G4PSCylinderSurfaceCurrent.cc.

179 {
180  G4cout << " MultiFunctionalDet " << detector->GetName() << G4endl;
181  G4cout << " PrimitiveScorer " << GetName() <<G4endl;
182  G4cout << " Number of entries " << EvtMap->entries() << G4endl;
183  std::map<G4int,G4double*>::iterator itr = EvtMap->GetMap()->begin();
184  for(; itr != EvtMap->GetMap()->end(); itr++) {
185  G4cout << " copy no.: " << itr->first
186  << " current : " ;
187  if ( divideByArea ) {
188  G4cout << *(itr->second)/GetUnitValue()
189  << " ["<<GetUnit()<<"]";
190  } else {
191  G4cout << *(itr->second) << " [tracks]";
192  }
193  G4cout << G4endl;
194  }
195 }
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 G4PSCylinderSurfaceCurrent::ProcessHits ( G4Step aStep,
G4TouchableHistory  
)
protectedvirtual

Implements G4VPrimitiveScorer.

Definition at line 80 of file G4PSCylinderSurfaceCurrent.cc.

81 {
82  G4StepPoint* preStep = aStep->GetPreStepPoint();
83  G4VPhysicalVolume* physVol = preStep->GetPhysicalVolume();
84  G4VPVParameterisation* physParam = physVol->GetParameterisation();
85  G4VSolid * solid = 0;
86  if(physParam)
87  { // for parameterized volume
89  ->GetReplicaNumber(indexDepth);
90  solid = physParam->ComputeSolid(idx, physVol);
91  solid->ComputeDimensions(physParam,idx,physVol);
92  }
93  else
94  { // for ordinary volume
95  solid = physVol->GetLogicalVolume()->GetSolid();
96  }
97 
98  G4Tubs* tubsSolid = (G4Tubs*)(solid);
99 
100  G4int dirFlag =IsSelectedSurface(aStep,tubsSolid);
101  // G4cout << " pos " << preStep->GetPosition() <<" dirFlag " << G4endl;
102  if ( dirFlag > 0 ) {
103  if ( fDirection == fCurrent_InOut || fDirection == dirFlag ){
104  G4TouchableHandle theTouchable = preStep->GetTouchableHandle();
105  //
106  G4double current = 1.0;
107  if ( weighted ) current = preStep->GetWeight(); // Current (Particle Weight)
108  //
109  if ( divideByArea ){
110  G4double square = 2.*tubsSolid->GetZHalfLength()
111  *tubsSolid->GetInnerRadius()* tubsSolid->GetDeltaPhiAngle()/radian;
112  current = current/square; // Current normalized by Area
113  }
114 
115  G4int index = GetIndex(aStep);
116  EvtMap->add(index,current);
117  }
118 
119  }
120 
121  return TRUE;
122 }
G4double GetWeight() const
virtual void ComputeDimensions(G4VPVParameterisation *p, const G4int n, const G4VPhysicalVolume *pRep)
Definition: G4VSolid.cc:138
G4int IsSelectedSurface(G4Step *, G4Tubs *)
virtual G4VSolid * ComputeSolid(const G4int, G4VPhysicalVolume *)
Definition: G4Tubs.hh:85
G4VSolid * GetSolid() const
const G4VTouchable * GetTouchable() const
int G4int
Definition: G4Types.hh:78
G4StepPoint * GetPreStepPoint() const
G4double GetDeltaPhiAngle() const
G4VPhysicalVolume * GetPhysicalVolume() const
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
G4double GetZHalfLength() const
double G4double
Definition: G4Types.hh:76
const G4TouchableHandle & GetTouchableHandle() const

Here is the call graph for this function:

void G4PSCylinderSurfaceCurrent::SetUnit ( const G4String unit)
virtual

Definition at line 197 of file G4PSCylinderSurfaceCurrent.cc.

198 {
199  if ( divideByArea ) {
200  CheckAndSetUnit(unit,"Per Unit Surface");
201  } else {
202  if (unit == "" ){
203  unitName = unit;
204  unitValue = 1.0;
205  }else{
206  G4String msg = "Invalid unit ["+unit+"] (Current unit is [" +GetUnit()+"] ) for " + GetName();
207  G4Exception("G4PSCylinderSurfaceCurrent::SetUnit","DetPS0002",
208  JustWarning,msg);
209  }
210  }
211 }
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 G4PSCylinderSurfaceCurrent::Weighted ( G4bool  flg = true)
inline

Definition at line 68 of file G4PSCylinderSurfaceCurrent.hh.

68 { weighted = flg; }

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