G4PSVolumeFlux Class Reference

#include <G4PSVolumeFlux.hh>

Inheritance diagram for G4PSVolumeFlux:

Public Member Functions
	G4PSVolumeFlux (G4String name, G4int direction=1, G4int depth=0)
virtual	~G4PSVolumeFlux ()
virtual void	Initialize (G4HCofThisEvent *)
virtual void	EndOfEvent (G4HCofThisEvent *)
virtual void	clear ()
virtual void	DrawAll ()
virtual void	PrintAll ()
void	SetDivAre (G4bool val)
void	SetDivCos (G4bool val)
Public Member Functions inherited from G4VPrimitivePlotter
	G4VPrimitivePlotter (G4String name, G4int depth=0)
virtual	~G4VPrimitivePlotter ()
void	Plot (G4int copyNo, G4int histID)
G4int	GetNumberOfHist () const
Public Member Functions inherited from G4VPrimitiveScorer
	G4VPrimitiveScorer (G4String name, G4int depth=0)
virtual	~G4VPrimitiveScorer ()
G4int	GetCollectionID (G4int)
virtual void	Initialize (G4HCofThisEvent *)
virtual void	EndOfEvent (G4HCofThisEvent *)
virtual void	clear ()
virtual void	DrawAll ()
virtual void	PrintAll ()
void	SetUnit (const G4String &unit)
const G4String &	GetUnit () const
G4double	GetUnitValue () const
void	SetMultiFunctionalDetector (G4MultiFunctionalDetector *d)
G4MultiFunctionalDetector *	GetMultiFunctionalDetector () const
G4String	GetName () const
void	SetFilter (G4VSDFilter *f)
G4VSDFilter *	GetFilter () const
void	SetVerboseLevel (G4int vl)
G4int	GetVerboseLevel () const
void	SetNijk (G4int i, G4int j, G4int k)

Protected Member Functions
virtual G4bool	ProcessHits (G4Step *, G4TouchableHistory *)
Protected Member Functions inherited from G4VPrimitiveScorer
virtual G4bool	ProcessHits (G4Step *, G4TouchableHistory *)=0
virtual G4int	GetIndex (G4Step *)
void	CheckAndSetUnit (const G4String &unit, const G4String &category)
G4VSolid *	ComputeSolid (G4Step *aStep, G4int replicaIdx)
G4VSolid *	ComputeCurrentSolid (G4Step *aStep)

Additional Inherited Members
Protected Attributes inherited from G4VPrimitivePlotter
std::map< G4int, G4int >	hitIDMap
Protected Attributes inherited from G4VPrimitiveScorer
G4String	primitiveName
G4MultiFunctionalDetector *	detector
G4VSDFilter *	filter
G4int	verboseLevel
G4int	indexDepth
G4String	unitName
G4double	unitValue
G4int	fNi
G4int	fNj
G4int	fNk

Detailed Description

Definition at line 43 of file G4PSVolumeFlux.hh.

Constructor & Destructor Documentation

G4PSVolumeFlux()

G4PSVolumeFlux::G4PSVolumeFlux	(	G4String	name,
		G4int	direction = 1,
		G4int	depth = 0
	)

Definition at line 49 of file G4PSVolumeFlux.cc.

: G4VPrimitivePlotter(name,depth),HCID(-1),fDirection(direction),EvtMap(nullptr)
{;}

~G4PSVolumeFlux()

G4PSVolumeFlux::~G4PSVolumeFlux ( )

virtual

Definition at line 53 of file G4PSVolumeFlux.cc.

{;}

Member Function Documentation

clear()

void G4PSVolumeFlux::clear ( )

virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 140 of file G4PSVolumeFlux.cc.

                          {
  EvtMap->clear();
}

DrawAll()

void G4PSVolumeFlux::DrawAll ( )

virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 144 of file G4PSVolumeFlux.cc.

{;}

EndOfEvent()

void G4PSVolumeFlux::EndOfEvent ( G4HCofThisEvent *  )

virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 137 of file G4PSVolumeFlux.cc.

{;}

Initialize()

void G4PSVolumeFlux::Initialize ( G4HCofThisEvent *  HCE )

virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 130 of file G4PSVolumeFlux.cc.

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

PrintAll()

void G4PSVolumeFlux::PrintAll ( )

virtual

Reimplemented from G4VPrimitiveScorer.

Definition at line 147 of file G4PSVolumeFlux.cc.

{
  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
       << "  flux  : " << *(itr->second)
       << G4endl;
  }
}

ProcessHits()

G4bool G4PSVolumeFlux::ProcessHits ( G4Step *  aStep, G4TouchableHistory *    )

protectedvirtual

Implements G4VPrimitiveScorer.

Definition at line 56 of file G4PSVolumeFlux.cc.

{
  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
  G4StepPoint* thisStepPoint = nullptr;
  if(fDirection==1) // Score only the inward particle
  {
    if(preStepPoint->GetStepStatus()==fGeomBoundary)
    { thisStepPoint = preStepPoint; }
    else
    { return false; }
  }
  else if(fDirection==2) // Score only the outward particle
  {
    if(postStepPoint->GetStepStatus()==fGeomBoundary)
    { thisStepPoint = postStepPoint; }
    else
    { return false; }
  }
 
  G4double flux = preStepPoint->GetWeight();
 
  if(divare||divcos)
  {
    G4VPhysicalVolume* physVol = preStepPoint->GetPhysicalVolume();
    G4VPVParameterisation* physParam = physVol->GetParameterisation();
    G4VSolid * solid = nullptr;
    if(physParam)
    { // for parameterized volume
      auto idx = ((G4TouchableHistory*)(preStepPoint->GetTouchable()))->GetReplicaNumber(indexDepth);
      solid = physParam->ComputeSolid(idx, physVol);
      solid->ComputeDimensions(physParam,idx,physVol);
    }
    else
    { // for ordinary volume
      solid = physVol->GetLogicalVolume()->GetSolid();
    }
 
    if(divare)
    { flux /= solid->GetSurfaceArea(); }
 
    if(divcos)
    {
      G4TouchableHandle theTouchable = thisStepPoint->GetTouchableHandle();
      G4ThreeVector pdirection = thisStepPoint->GetMomentumDirection();
      G4ThreeVector localdir = theTouchable->GetHistory()->GetTopTransform().TransformAxis(pdirection);
      G4ThreeVector globalPos = thisStepPoint->GetPosition();
      G4ThreeVector localPos = theTouchable->GetHistory()->GetTopTransform().TransformPoint(globalPos);
      G4ThreeVector surfNormal = solid->SurfaceNormal(localPos);
      G4double cosT = surfNormal.cosTheta(localdir);
      if(cosT!=0.) flux /= std::abs(cosT);
    }
  }
 
  G4int index = GetIndex(aStep);
  EvtMap->add(index,flux);
 
  if(hitIDMap.size()>0 && hitIDMap.find(index)!=hitIDMap.end())
  {
    auto filler = G4VScoreHistFiller::Instance();
    if(!filler)
    {
      G4Exception("G4PSVolumeFlux::ProcessHits","SCORER0123",JustWarning,
             "G4TScoreHistFiller is not instantiated!! Histogram is not filled.");
    }
    else
    {
      filler->FillH1(hitIDMap[index],thisStepPoint->GetKineticEnergy(),flux);
    }
  }
 
  return true;
}

SetDivAre()

void G4PSVolumeFlux::SetDivAre ( G4bool  val )

inline

Definition at line 67 of file G4PSVolumeFlux.hh.

      { divare = val; }

SetDivCos()

void G4PSVolumeFlux::SetDivCos ( G4bool  val )

inline

Definition at line 69 of file G4PSVolumeFlux.hh.

      { divcos = val; }

Referenced by G4ScoreQuantityMessenger::SetNewValue().

---

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

• G4PSVolumeFlux.hh
• G4PSVolumeFlux.cc