G4AtlasFluxRecorder Class Reference

#include <G4AtlasFluxRecorder.h>

Inheritance diagram for G4AtlasFluxRecorder:

Collaboration diagram for G4AtlasFluxRecorder:

Public Member Functions
Methods related to flux recording
Initialize flux recording
void	InitializeFluxRecording () override final
void	RecordFlux (const G4Event *) override final
	Record fluxes from current event. More...
void	WriteFluxInformation () override final
	Dump flux information to text files. More...

Detailed Description

Definition at line 10 of file G4AtlasFluxRecorder.h.

Member Function Documentation

InitializeFluxRecording()

void G4AtlasFluxRecorder::InitializeFluxRecording ( )

finaloverridevirtual

Implements IFluxRecorder.

Definition at line 23 of file G4AtlasFluxRecorder.cxx.

{
  G4UImanager *ui = G4UImanager::GetUIpointer();
  ui->ApplyCommand("/run/setCutForAGivenParticle proton 0 mm");
 
  G4ScoringManager* ScM = G4ScoringManager::GetScoringManagerIfExist();
 
  if(!ScM) { return; }
 
  ui->ApplyCommand("/score/create/cylinderMesh cylMesh_1");
  //                        R  Z(-24 to 24)
  ui->ApplyCommand("/score/mesh/cylinderSize 12. 24. m");
  //                iR iZ
  //ui->ApplyCommand("/score/mesh/nBin 1 1 1");
  ui->ApplyCommand("/score/mesh/nBin 120 480 1");
 
  ui->ApplyCommand("/score/quantity/energyDeposit eDep");
 
  ui->ApplyCommand("/score/quantity/cellFlux CF_photon");
  ui->ApplyCommand("/score/filter/particle photonFilter gamma");
  // above 2 line crete tally for cell flux for gamma
 
  ui->ApplyCommand("/score/quantity/cellFlux CF_neutron");
  ui->ApplyCommand("/score/filter/particle neutronFilter neutron");
 
  ui->ApplyCommand("/score/quantity/cellFlux CF_HEneutron");
  ui->ApplyCommand("/score/filter/particleWithKineticEnergy HEneutronFilter 20 7000000 MeV neutron");
 
  ui->ApplyCommand("/score/quantity/doseDeposit dose");
 
  ui->ApplyCommand("/score/close");
  ui->ApplyCommand("/score/list");
 
  G4int nPar = ScM->GetNumberOfMesh();
 
  if(nPar<1) { return; }
 
  G4ParticleTable::G4PTblDicIterator* particleIterator
    = G4ParticleTable::GetParticleTable()->GetIterator();
 
  for(G4int iw=0;iw<nPar;iw++) {
    G4VScoringMesh* mesh = ScM->GetMesh(iw);
    G4VPhysicalVolume* pWorld
      = G4TransportationManager::GetTransportationManager()
      ->IsWorldExisting(ScM->GetWorldName(iw));
    if(!pWorld) {
      pWorld = G4TransportationManager::GetTransportationManager()
        ->GetParallelWorld(ScM->GetWorldName(iw));
      pWorld->SetName(ScM->GetWorldName(iw));
 
      G4ParallelWorldScoringProcess* theParallelWorldScoringProcess
        = new G4ParallelWorldScoringProcess(ScM->GetWorldName(iw));
      theParallelWorldScoringProcess->SetParallelWorld(ScM->GetWorldName(iw));
 
      particleIterator->reset();
      while( (*particleIterator)() ) {
        G4ParticleDefinition* particle = particleIterator->value();
        G4ProcessManager* pmanager = particle->GetProcessManager();
        if(pmanager) {
          pmanager->AddProcess(theParallelWorldScoringProcess);
          pmanager->SetProcessOrderingToLast(theParallelWorldScoringProcess, idxAtRest);
          pmanager->SetProcessOrderingToSecond(theParallelWorldScoringProcess, idxAlongStep);
          pmanager->SetProcessOrderingToLast(theParallelWorldScoringProcess, idxPostStep);
        }
      }
    }
 
    mesh->Construct(pWorld);
  }
  return;
}

RecordFlux()

void G4AtlasFluxRecorder::RecordFlux ( const G4Event *  currentEvent )

finaloverridevirtual

Record fluxes from current event.

Implements IFluxRecorder.

Definition at line 95 of file G4AtlasFluxRecorder.cxx.

{
  G4ScoringManager* ScM = G4ScoringManager::GetScoringManagerIfExist();
  if(ScM) {
    G4int nPar = ScM->GetNumberOfMesh();
    G4HCofThisEvent* HCE = currentEvent->GetHCofThisEvent();
    if(HCE && nPar>0) {
      G4int nColl = HCE->GetCapacity();
      for(G4int i=0;i<nColl;i++) {
        G4VHitsCollection* HC = HCE->GetHC(i);
        if(HC) { ScM->Accumulate(HC); }
      }
    }
  }
  return;
}

WriteFluxInformation()

void G4AtlasFluxRecorder::WriteFluxInformation ( )

finaloverridevirtual

Dump flux information to text files.

Implements IFluxRecorder.

Definition at line 112 of file G4AtlasFluxRecorder.cxx.

{
  G4UImanager *ui=G4UImanager::GetUIpointer();
  ui->ApplyCommand("/score/dumpQuantityToFile cylMesh_1 eDep edep.txt");
  ui->ApplyCommand("/score/dumpQuantityToFile cylMesh_1 CF_neutron neutron.txt");
  ui->ApplyCommand("/score/dumpQuantityToFile cylMesh_1 CF_HEneutron HEneutron.txt");
  ui->ApplyCommand("/score/dumpQuantityToFile cylMesh_1 CF_photon photon.txt");
  ui->ApplyCommand("/score/dumpQuantityToFile cylMesh_1 dose dose.txt");
  return;
}

---

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

• G4AtlasFluxRecorder.h
• G4AtlasFluxRecorder.cxx