G4TankPMT.cc

Go to the documentation of this file.

#include "G4TankPMT.h"
#include "G4TankTrackingAction.h"
#include "ParticleCases.h"

using namespace G4TankSimulatorASCII;


void
G4TankPMT::SetParticle(const utl::Particle& particle)
{
  fParticle = &particle;
  SetComponent();
  fExtraComponent = 0;
 
}


void
G4TankPMT::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;
  default:
    fComponent = sevt::StationConstants::eHadron;
    break; 
  }
}


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


void
G4TankPMT::AddPhoton(const int nPMT, const double peTime, const bool IsFromMuonDecay )
{
  if ((peTime > 1.0e6) || (peTime < -1.0e6)) {
    cerr << "Dubious PE release time " << peTime << endl;
    exit(-1);
  }

  if ( IsFromMuonDecay ) 
    {
      G4TankTrackingAction::MichelElectronNPE[1]+=1;
      G4TankTrackingAction::MichelElectronNPE_time[1]+=(peTime);
    }

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