G4XTankPMT.cc

Go to the documentation of this file.

#include <utl/ParticleCases.h>

#include "G4XTankPMT.h"

using namespace G4XTankSimulatorAG;


void
G4XTankPMT::SetParticle(const utl::Particle& particle)
{
  fParticle = &particle;
  SetComponent();
  fExtraComponent = 0;
  if (particle.GetSource() != utl::Particle::eShowerFromMuonDecay)
    return;
  switch (fComponent) {
  case sevt::StationConstants::eElectron:
    fExtraComponent = sevt::StationConstants::eShowerMuonDecayElectron;
    break;
  case sevt::StationConstants::ePhoton:
    fExtraComponent = sevt::StationConstants::eShowerMuonDecayPhoton;
    break;
  default:
    break;
  }
}


void
G4XTankPMT::SetComponent()
{
  switch (fParticle->GetType()) {
  case OFFLINE_ELECTRONS:
    fComponent = sevt::StationConstants::eElectron;
    break;
  case OFFLINE_PHOTON:
    fComponent = sevt::StationConstants::ePhoton;
    break;
  case OFFLINE_MUONS:
    fComponent = sevt::StationConstants::eMuon;
    break;
  case OFFLINE_HADRONS:
    fComponent = sevt::StationConstants::eHadron;
    break;
  default:
    cerr << "Illegal Component in G4XTankPMT" << endl;
    exit(-1);
  }
}


void
G4XTankPMT::SetCurrentTank(const sevt::SEvent::StationIterator stationIt)
{
  fCurrentEventStationIt = stationIt;
  for (int i = 1; i <= 3; ++i)
    fPMTSimData[i] = &fCurrentEventStationIt->GetPMT(i).GetSimData();
}


void
G4XTankPMT::AddPhoton(const int nPMT, const double peTime)
{
  if ((peTime > 1e6) || (peTime < -1e6)) {
    cerr << "Dubious PE release time " << peTime << endl;
    exit(-1);
  }
  if (!fPMTSimData[nPMT]->HasPETimeDistribution())
    fPMTSimData[nPMT]->MakePETimeDistribution();
  fPMTSimData[nPMT]->GetPETimeDistribution().AddTime(peTime);
  if (!fPMTSimData[nPMT]->HasPETimeDistribution(fComponent))
    fPMTSimData[nPMT]->MakePETimeDistribution(fComponent);
  fPMTSimData[nPMT]->GetPETimeDistribution(fComponent).AddTime(peTime);

  if (fExtraComponent) {
    const sevt::StationConstants::SignalComponent comp =
      static_cast<sevt::StationConstants::SignalComponent>(fExtraComponent);
    if (!fPMTSimData[nPMT]->HasPETimeDistribution(comp))
      fPMTSimData[nPMT]->MakePETimeDistribution(comp);
    fPMTSimData[nPMT]->GetPETimeDistribution(comp).AddTime(peTime);
  }
}