EventAction.cc

Go to the documentation of this file.

#include "EventAction.hh"

EventAction::EventAction(RunAction* run, DetectorConstruction* det, G4double accumulationtime)
  :runAct(run),printModulo(100),eventMessenger(0), myAccumulationTime(accumulationtime)
{
  detector=det;
  eventMessenger = new EventActionMessenger(this);
}

EventAction::~EventAction()
{
  delete eventMessenger;
}

// Just updates the eventnumber at the beginning of each event
void EventAction::BeginOfEventAction(const G4Event* evt)
{
  G4int evtNb = evt->GetEventID();
  std::stringstream ss;
  ss << "Event ID: " << evtNb << std::endl;
  runAct->PrintToVerboseFile(ss.str());

  if(evtNb%printModulo==0)
    G4cout << "Starting event " << evtNb << G4endl;
  
  myStatistics = new AllAccumulatedStatistics(((detector->GetNumberOfScintillatorBoxes())*(detector->GetNumberOfScintillatorsPerBox())*2),myAccumulationTime); 
}

void EventAction::PrintToVerboseFile(G4String messageToPrint)
{
  runAct->PrintToVerboseFile(messageToPrint);
}


void EventAction::EndOfEventAction(const G4Event* evt)
{
  G4int evtNb = evt->GetEventID();
 
  runAct->WriteResultsToFile(evtNb,myStatistics);
  runAct->PrintToVerboseFile("\n");
  
  delete myStatistics;
}  

void EventAction::AddScintillatorTubeEvent(G4int scintillatorID, G4double time, G4double energy)
{
  myStatistics->AddScintillatorTubeEvent(scintillatorID,time,energy); //a range check would probable be good here.
}

void EventAction::AddGeDetectorEvent(G4double time, G4double energy)
{
  myStatistics->AddGeDetectorEvent(time,energy);
}

void EventAction::AddScintillatorPlateEvent(AllAccumulatedStatistics::ScintillatorPlate plate, G4double time, G4double energy)
{
  myStatistics->AddScintillatorPlateEvent(plate,time,energy);
}