G4TankSimulatorOG/G4TankPMTAction.cc

Go to the documentation of this file.

#include "G4TankPMTAction.h"
#include "G4TankEventAction.h"
#include "G4TankSimulator.h"
#include "G4TankSteppingAction.h"

#include <sevt/PMT.h>
#include <sevt/PMTSimData.h>
#include <sevt/SEvent.h>
#include <sevt/Station.h>
#include <sevt/StationSimData.h>

#include <utl/ErrorLogger.h>
#include <utl/Particle.h>
#include <utl/ParticleCases.h>
#include <utl/Trace.h>

#include <fwk/RunController.h>

#include <G4Event.hh>
#include <G4HCofThisEvent.hh>
#include <G4SDManager.hh>
#include <G4Step.hh>
#include <G4TouchableHistory.hh>
#include <G4Track.hh>
#include <G4Trajectory.hh>
#include <G4TrajectoryContainer.hh>
#include <G4VHitsCollection.hh>

#include <iostream>

using namespace std;
using namespace G4TankSimulatorOG;
using namespace sevt;
using namespace utl;
using namespace fwk;


G4TankPMTAction::G4TankPMTAction(const G4String name, const int pmtIndex) :
  G4VSensitiveDetector(name),
  fPMTIndex(pmtIndex)
{
  if (fPMTIndex < 1 || fPMTIndex > 3) {
    FATAL("G4PMTAction() : invalid PMT index.");
    exit(EXIT_FAILURE);
  }

  fG4TankSimulator =
    dynamic_cast<G4TankSimulator*>(&RunController::GetInstance().GetModule("G4TankSimulatorOG"));
}


G4bool
G4TankPMTAction::ProcessHits(G4Step* theStep, G4TouchableHistory* /*R0Hist*/)
{
  const double time = theStep->GetPreStepPoint()->GetGlobalTime();

  // get current particle from G4TankStepping action.
  // use to reject hits registered in sensitive volume
  // in case the particle is not a photon
  const int currentParticle = G4TankSteppingAction::GetCurrentParticleId();

  // This is related to bug #156
  if (time >= 1*utl::second || currentParticle)
    return true;

  fG4TankSimulator->AddPhoton(fPMTIndex, time);
  
  return true;
}