G4TankSimulatorOG/G4TankStackingAction.cc

Go to the documentation of this file.

#include "G4TankStackingAction.h"
#include "G4TankSimulator.h"

#include <fwk/RandomEngineRegistry.h>

#include <sdet/Station.h>
#include <sdet/PMT.h>

#include <utl/TabulatedFunction.h>

#include <CLHEP/Random/RandFlat.h>

#include <G4Track.hh>
#include <G4TrackingManager.hh>
#include <G4TrackStatus.hh>
#include <G4OpticalPhoton.hh>

using namespace sdet;
using namespace fwk;
using namespace G4TankSimulatorOG;
using CLHEP::RandFlat;


G4TankStackingAction::G4TankStackingAction() :
  fCurrentDetectorStation(0)
{
  fRandomEngine =
    &RandomEngineRegistry::GetInstance().Get(RandomEngineRegistry::eDetector);
}


G4ClassificationOfNewTrack
G4TankStackingAction::ClassifyNewTrack(const G4Track* aTrack)
{
  // NOTE: for now, we assume that each pmt has the same collection
  // and quantum efficiencies

  if (aTrack->GetDefinition() != G4OpticalPhoton::OpticalPhoton())
    return fWaiting;

  // optical photon

  const double energy = aTrack->GetKineticEnergy() *
    (utl::MeV / CLHEP::MeV);

  // Randomly kill photons according to quantum and collection
  // efficiency to speed up simulation

  if (RandFlat::shoot(&fRandomEngine->GetEngine(), 0, 1) >
      fCollectionEfficiency * fQuantumEff.InterpolateY(energy, 1))
    return fKill;

  return fWaiting;
}


void
G4TankStackingAction::SetStation(const sdet::Station* const station)
{
  fCurrentDetectorStation = station;

  fCollectionEfficiency = fCurrentDetectorStation->GetPMT(1).GetCollectionEfficiency();

  fQuantumEff = fCurrentDetectorStation->GetPMT(1).GetQuantumEfficiency();
}