G4PSDoseDeposit Class Reference

#include <G4PSDoseDeposit.hh>

Inheritance diagram for G4PSDoseDeposit:

Public Member Functions
	G4PSDoseDeposit (G4String name, G4int depth=0)
	G4PSDoseDeposit (G4String name, const G4String &unit, G4int depth=0)
virtual	~G4PSDoseDeposit ()
virtual void	Initialize (G4HCofThisEvent *)
virtual void	EndOfEvent (G4HCofThisEvent *)
virtual void	clear ()
virtual void	DrawAll ()
virtual void	PrintAll ()
virtual void	SetUnit (const G4String &unit)
Protected Member Functions
virtual G4bool	ProcessHits (G4Step *, G4TouchableHistory *)
virtual G4double	ComputeVolume (G4Step *, G4int idx)

---

Detailed Description

Definition at line 46 of file G4PSDoseDeposit.hh.

---

Constructor & Destructor Documentation

G4PSDoseDeposit::G4PSDoseDeposit	(	G4String	name,
		G4int	depth = 0
	)

Definition at line 47 of file G4PSDoseDeposit.cc.

References SetUnit().

  :G4VPrimitiveScorer(name,depth),HCID(-1)
{
    SetUnit("Gy");
}

G4PSDoseDeposit::G4PSDoseDeposit	(	G4String	name,
		const G4String &	unit,
		G4int	depth = 0
	)

Definition at line 53 of file G4PSDoseDeposit.cc.

References SetUnit().

  :G4VPrimitiveScorer(name,depth),HCID(-1)
{
    SetUnit(unit);
}

G4PSDoseDeposit::~G4PSDoseDeposit

(

)

[virtual]

Definition at line 60 of file G4PSDoseDeposit.cc.

{;}

---

Member Function Documentation

void G4PSDoseDeposit::clear

(

)

[virtual]

Reimplemented from G4VPrimitiveScorer.

Definition at line 93 of file G4PSDoseDeposit.cc.

{
  EvtMap->clear();
}

G4double G4PSDoseDeposit::ComputeVolume	(	G4Step *	,
		G4int	idx
	)			[protected, virtual]

Reimplemented in G4PSDoseDepositForCylinder3D.

Definition at line 121 of file G4PSDoseDeposit.cc.

References G4VSolid::ComputeDimensions(), G4VPVParameterisation::ComputeSolid(), G4endl, G4Exception(), G4VSolid::GetCubicVolume(), G4VPhysicalVolume::GetLogicalVolume(), G4VPhysicalVolume::GetParameterisation(), G4StepPoint::GetPhysicalVolume(), G4Step::GetPreStepPoint(), G4LogicalVolume::GetSolid(), and JustWarning.

Referenced by ProcessHits().

                                                               {

  G4VPhysicalVolume* physVol = aStep->GetPreStepPoint()->GetPhysicalVolume();
  G4VPVParameterisation* physParam = physVol->GetParameterisation();
  G4VSolid* solid = 0;
  if(physParam)
  { // for parameterized volume
    if(idx<0)
    {
      G4ExceptionDescription ED;
      ED << "Incorrect replica number --- GetReplicaNumber : " << idx << G4endl;
      G4Exception("G4PSDoseDeposit::ComputeVolume","DetPS0004",JustWarning,ED);
    }
    solid = physParam->ComputeSolid(idx, physVol);
    solid->ComputeDimensions(physParam,idx,physVol);
  }
  else
  { // for ordinary volume
    solid = physVol->GetLogicalVolume()->GetSolid();
  }
  
  return solid->GetCubicVolume();
}

void G4PSDoseDeposit::DrawAll

(

)

[virtual]

Reimplemented from G4VPrimitiveScorer.

Definition at line 98 of file G4PSDoseDeposit.cc.

{;}

void G4PSDoseDeposit::EndOfEvent

(

G4HCofThisEvent *

)

[virtual]

Reimplemented from G4VPrimitiveScorer.

Definition at line 90 of file G4PSDoseDeposit.cc.

{;}

void G4PSDoseDeposit::Initialize

(

G4HCofThisEvent *

)

[virtual]

Reimplemented from G4VPrimitiveScorer.

Definition at line 82 of file G4PSDoseDeposit.cc.

References G4HCofThisEvent::AddHitsCollection(), G4VPrimitiveScorer::GetCollectionID(), G4VPrimitiveScorer::GetMultiFunctionalDetector(), G4VPrimitiveScorer::GetName(), and G4VSensitiveDetector::GetName().

{
  EvtMap = new G4THitsMap<G4double>(GetMultiFunctionalDetector()->GetName(),
                                    GetName());
  if(HCID < 0) {HCID = GetCollectionID(0);}
  HCE->AddHitsCollection(HCID, (G4VHitsCollection*)EvtMap);
}

void G4PSDoseDeposit::PrintAll

(

)

[virtual]

Reimplemented from G4VPrimitiveScorer.

Definition at line 101 of file G4PSDoseDeposit.cc.

References G4VPrimitiveScorer::detector, G4cout, G4endl, G4VPrimitiveScorer::GetName(), G4VSensitiveDetector::GetName(), G4VPrimitiveScorer::GetUnit(), and G4VPrimitiveScorer::GetUnitValue().

{
  G4cout << " MultiFunctionalDet  " << detector->GetName() << G4endl;
  G4cout << " PrimitiveScorer " << GetName() << G4endl;
  G4cout << " Number of entries " << EvtMap->entries() << G4endl;
  std::map<G4int,G4double*>::iterator itr = EvtMap->GetMap()->begin();
  for(; itr != EvtMap->GetMap()->end(); itr++) {
    G4cout << "  copy no.: " << itr->first
           << "  dose deposit: " 
           << *(itr->second)/GetUnitValue()
           << " ["<<GetUnit() <<"]"
           << G4endl;
  }
}

G4bool G4PSDoseDeposit::ProcessHits	(	G4Step *	,
		G4TouchableHistory *
	)			[protected, virtual]

Implements G4VPrimitiveScorer.

Definition at line 63 of file G4PSDoseDeposit.cc.

References ComputeVolume(), FALSE, G4Material::GetDensity(), G4VPrimitiveScorer::GetIndex(), G4StepPoint::GetMaterial(), G4Step::GetPreStepPoint(), G4Track::GetStep(), G4Step::GetTotalEnergyDeposit(), G4StepPoint::GetTouchable(), G4Step::GetTrack(), G4StepPoint::GetWeight(), G4VPrimitiveScorer::indexDepth, and TRUE.

{
  G4double edep = aStep->GetTotalEnergyDeposit();
  if ( edep == 0. ) return FALSE;

  G4int idx = ((G4TouchableHistory*)
               (aStep->GetPreStepPoint()->GetTouchable()))
               ->GetReplicaNumber(indexDepth);
  G4double cubicVolume = ComputeVolume(aStep, idx);


  G4double density = aStep->GetTrack()->GetStep()->GetPreStepPoint()->GetMaterial()->GetDensity();
  G4double dose    = edep / ( density * cubicVolume );
  dose *= aStep->GetPreStepPoint()->GetWeight(); 
  G4int  index = GetIndex(aStep);
  EvtMap->add(index,dose);  
  return TRUE;
}

void G4PSDoseDeposit::SetUnit

(

const G4String &

unit

)

[virtual]

Reimplemented from G4VPrimitiveScorer.

Definition at line 116 of file G4PSDoseDeposit.cc.

References G4VPrimitiveScorer::CheckAndSetUnit().

Referenced by G4PSDoseDeposit(), and G4PSDoseDeposit3D::G4PSDoseDeposit3D().

{
        CheckAndSetUnit(unit,"Dose");
}

---

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

• G4PSDoseDeposit.hh
• G4PSDoseDeposit.cc

---