AllAccumulatedStatistics.cc

Go to the documentation of this file.

#include "AllAccumulatedStatistics.hh"

AllAccumulatedStatistics::AllAccumulatedStatistics(G4double numberOfScintillatorTubes, G4double summationTime)
  :scintillatorTubes(numberOfScintillatorTubes),sumTime(summationTime)
{
  
}

AllAccumulatedStatistics::~AllAccumulatedStatistics()
{
  
}

bool AllAccumulatedStatistics::timeCheck(G4double time)
{
  if(eventHits.empty() || time > (eventHits.back().time + sumTime ))
    {
      eventHits.push_back(EventHit(scintillatorTubes,time));
      return true;
    }
  return false;
}

void AllAccumulatedStatistics::AddScintillatorTubeEvent(G4int ScintillatorID, G4double time, G4double energy)
{
  if(energy>0.1*eV) //don't add non-existent energies.
    {
      timeCheck(time);
      eventHits.back().scintillatorTubeEnergies.at(ScintillatorID)+=energy;
    }
}

void AllAccumulatedStatistics::AddGeDetectorEvent(G4double time, G4double energy)
{
  if(energy>0.1*eV) //don't add non-existent energies.
    {
      timeCheck(time);
      eventHits.back().geDetectorEnergy+=energy;
    }
}

void AllAccumulatedStatistics::AddScintillatorPlateEvent(ScintillatorPlate plate, G4double time, G4double energy)
{
    if(energy>0.1*eV) //don't add non-existent energies.
    {
      timeCheck(time);
      if(plate==UpperFront)
        eventHits.back().upperFrontScintillatorPanelEnergy+=energy;
      if(plate==LowerFront)
        eventHits.back().lowerFrontScintillatorPanelEnergy+=energy;
      if(plate==Back)
        eventHits.back().backScintillatorPanelEnergy+=energy;
    }
}

std::vector<EventHit> AllAccumulatedStatistics::GetEventHits()
{
  return eventHits;
}