Deprecated/UpgradeASCIITests/G4TankSimulatorASCII/G4TankStackingAction.cc

Go to the documentation of this file.

#include "G4TankStackingAction.h"
#include "G4TankTrackingAction.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 G4TankSimulatorASCII;
using CLHEP::RandFlat;


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


G4ClassificationOfNewTrack
G4TankStackingAction::ClassifyNewTrack(const G4Track* aTrack)
{
  //Kill Decay Products from the Primary particle 
  if ( aTrack->GetParentID() ==1 &&  aTrack->GetCreatorProcess()->GetProcessName() == "Decay")
    {
    
      if (aTrack->GetDefinition()->GetParticleName()=="e-" || aTrack->GetDefinition()->GetParticleName()=="e+" ) 
        {
          G4TankTrackingAction::MichelElectronKE=aTrack->GetKineticEnergy();
          G4TankTrackingAction::MichelElectronTheta=aTrack->GetMomentumDirection().theta();
          G4TankTrackingAction::MichelElectronZ=aTrack->GetPosition().z();
        
          
        }
      return fWaiting;
    }

  // 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();
}