FastTankSimulator.h

Go to the documentation of this file.

#ifndef _FastTankSimulatorOG_FastTankSimulator_h_
#define _FastTankSimulatorOG_FastTankSimulator_h_

#include <fwk/VModule.h>

#include <utl/CoordinateSystem.h>
#include <utl/TabulatedFunction.h>
#include <utl/Particle.h>

#include <vector>


namespace utl {
  class RandomEngine;
}

namespace sdet {
  class Station;
}

namespace sevt {
  class Station;
}


namespace FastTankSimulatorOG {

  class FastTankSimulator : public fwk::VModule {

  public:
    FastTankSimulator();

    VModule::ResultFlag Init();
    VModule::ResultFlag Run(evt::Event& event);
    VModule::ResultFlag Finish();

  private:
    double fdEdXMuon = 0;
    double fdEdXElectron = 0;
    double fCerenkovMinMuon = 0;
    double fCerenkovMinElectron = 0;
    double fEMinPhoton = 0;
    double fPhotoElectronRate = 0;
    double fWaterRIndex = 0;

    std::vector<double> fPECumulativeProbability;

    utl::TabulatedFunction fPhotonInteractionLength;
    utl::TabulatedFunction fPairProductionProbability;

    double fTankRadiusSquared = 0;
    double fTankHeight = 0;

    double CalculateDistanceInTank(const utl::Particle& p);
    double CalculateIntegratedCerenkovRate(const utl::Particle& p, const double);
    void CalculatePhotoElectrons(const utl::Particle& p, const double);

    void SimulateElectrons(const utl::Particle& p);
    void SimulateMuons(const utl::Particle& p);
    void SimulatePhotons(const utl::Particle& p);

    void SimulatePairProduction(const utl::Particle& p);
    void SimulateComptonScattering(const utl::Particle& p);

    utl::Particle fInitialParticle;
    const sdet::Station* fCurrentDetectorStation = nullptr;
    sevt::Station* fCurrentEventStation = nullptr;

    utl::RandomEngine* fRandomEngine = nullptr;

    REGISTER_MODULE("FastTankSimulatorOG", FastTankSimulator);

  };

}

#endif