RunAction.cc

Go to the documentation of this file.

//! \file
//!
//! Source file for the RunAction class.
//!
//! Defines the RunAction class and the member routines which
//! control the actions taken during the run.
//!
//! \author D.K. Hasell
//! \version 1.0
//! \date 2010-10-14
//!
//! \ingroup control

// This must appear first ?

#include "G4Timer.hh"

// Include user header files referenced in this file.

#include "RunAction.h"

// Include the GEANT4 header files referenced in this file.

#include "G4Run.hh"

#include "Randomize.hh"

// Include the ROOT header file referenced here.

#include "TFile.h"
#include "TTree.h"

// Include the C++ header files referenced in this file.

#include <ctime>

// Use the std namespace.

using namespace std;

// Define the RunAction class.

// Initialise the static data members.

G4Timer * RunAction::Timer = new G4Timer();

TFile * RunAction::File = 0;

TTree * RunAction::Tree = 0;

// Constructor.

RunAction::RunAction() {}

// Destructor.

RunAction::~RunAction() {}

// Define what happens at the beginning of a run.

void RunAction::BeginOfRunAction( const G4Run * run ) {

   // Start run timer.

   Timer->Start();

   // Get the current local date and time.

   time_t DateTime = time( NULL );

   // Initialise random seed based on current local date and time.

   CLHEP::HepRandom::setTheSeed( DateTime );

   // Get the current GMT date and time.

   struct tm * GMTDateTime = gmtime( &DateTime );

   // Print the start of run information.

   G4cout << "Start of run " << run->GetRunID() << "\n"
          << "   Number of events requested "
          << run->GetNumberOfEventToBeProcessed() << " "
          << "   Number of events processed "
          << run->GetNumberOfEvent() << "\n"
          << "   GMT time " << asctime( GMTDateTime )
          << "   Random seed " << DateTime << "\n"
          << G4endl;

   // Open a TFile to store the MC data.

   File = new TFile( "MCData.root", "RECREATE",
                     "OLYMPUS MC data file" );

   // Create a TTree to hold the MC data.

   Tree = new TTree( "Tree", "OLYMPUS MC TTree" );

}

// Define what happens at the end of a run.

void RunAction::EndOfRunAction( const G4Run * run ) {

   // Get the current local date and time.

   time_t DateTime = time( NULL );

   // Get the current GMT date and time.

   struct tm * GMTDateTime = gmtime( &DateTime );

   // Get the current random number seed.

   long seed = CLHEP::HepRandom::getTheSeed();

   // Stop the timer.

   Timer->Stop();

   // Print the end of run information.

   G4cout << "\n\nEnd of run " << run->GetRunID() << "\n"
          << "   Number of events requested "
          << run->GetNumberOfEventToBeProcessed() << " "
          << "   Number of events processed "
          << run->GetNumberOfEvent() << "\n"
          << "   GMT time " << asctime( GMTDateTime )
          << "   Random seed " << seed << "\n\n"
          << *Timer
          << G4endl;

   // Print information about Tree.

   Tree->Print();

   // Write and close the MC data file.

   File->Write();

   File->Close();

}