RunAction.cc

Go to the documentation of this file.

#include "RunAction.hh"

RunAction::RunAction()
{
  myMessenger = new RunActionMessenger(this);
  particleZ = 6;
  particleA = 17;
}

RunAction::~RunAction()
{
  delete myMessenger;
}

void RunAction::SetNewResultsFileName(G4String filename)
{
  resultsFileName = filename;
}

void RunAction::SetParticleZ(G4int newZ)
{
  particleZ=newZ;
}

void RunAction::SetParticleA(G4int newA)
{
  particleA=newA;
}

void RunAction::WriteResultsToFile(G4int eventID, AllAccumulatedStatistics * myEventStatistics)
{
  std::vector<EventHit> eventsToWrite = myEventStatistics->GetEventHits();
  for(std::vector<EventHit>::iterator it = eventsToWrite.begin(); it!=eventsToWrite.end(); it++)
    {
      fprintf(resultsFile, "%d %d %d ", eventID, particleZ, particleA);

      fprintf(resultsFile, "%.15e ",it->time/ns);
      if(it->upperFrontScintillatorPanelEnergy/MeV>1E-100)
        fprintf(resultsFile, "%.15e ",it->upperFrontScintillatorPanelEnergy/MeV);
      else
        fprintf(resultsFile, "%d ",0);
      
      if(it->lowerFrontScintillatorPanelEnergy/MeV>1E-100)
        fprintf(resultsFile, "%.15e ",it->lowerFrontScintillatorPanelEnergy/MeV);
      else
        fprintf(resultsFile, "%d ",0);

      if(it->backScintillatorPanelEnergy/MeV>1E-100)
        fprintf(resultsFile, "%.15e ",it->backScintillatorPanelEnergy/MeV);
      else
        fprintf(resultsFile, "%d ",0);

      if(it->geDetectorEnergy/MeV>1E-100)
        fprintf(resultsFile, "%.15e",it->geDetectorEnergy/MeV);
      else
        fprintf(resultsFile, "%d",0);

      for(std::vector<G4double>::iterator ig = it->scintillatorTubeEnergies.begin(); ig!=it->scintillatorTubeEnergies.end(); ig++)
        {
          if(*ig/MeV>1E-100)
            fprintf(resultsFile, " %.15e",*ig/MeV);
          else
            fprintf(resultsFile, " %d",0);
        }
      fprintf(resultsFile, "\n");
    }
}

void RunAction::BeginOfRunAction(const G4Run* aRun)
{ 
  G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
  resultsFile = fopen(resultsFileName.c_str(),"w");

  if(usingVerboseFile)
    {
      verboseFile = fopen(verboseResultsFileName, "w");  
    }
  
  //inform the runManager to save random number seed
  G4RunManager::GetRunManager()->SetRandomNumberStore(false);
    
}

void RunAction::EndOfRunAction(const G4Run* aRun)
{
  G4cout << "### End of run " << aRun->GetRunID() << "." << G4endl;
  fclose(resultsFile);
    if(usingVerboseFile)
    {
      fprintf(verboseFile, "\n");
      fclose(verboseFile);
    }
}


void RunAction::PrintToVerboseFile(G4String messageToPrint)
{
  if(usingVerboseFile)
    {
      fprintf(verboseFile,"%s",messageToPrint.c_str());
    }
}