PhotonLibraryHybrid.h

Go to the documentation of this file.

#ifndef PHOTONLIBRARYHYBRID_H
#define PHOTONLIBRARYHYBRID_H

#include "larsim/PhotonPropagation/IPhotonLibrary.h"

#include <cmath>
#include <string>
#include <vector>

namespace sim {
  class PhotonVoxelDef;
}

namespace phot {
  class PhotonLibraryHybrid : public IPhotonLibrary {
  public:
    PhotonLibraryHybrid(const std::string& fname, const sim::PhotonVoxelDef& voxdef);
    virtual ~PhotonLibraryHybrid();

    virtual float GetCount(size_t Voxel, size_t OpChannel) const override;

    // This one is unimplemented
    virtual const float* GetCounts(size_t Voxel) const override;

    virtual bool hasReflected() const override { return false; }
    virtual const float* GetReflCounts(size_t Voxel) const override { return 0; }
    virtual float GetReflCount(size_t Voxel, size_t OpChannel) const override { return 0; }

    virtual bool hasReflectedT0() const override { return false; }
    virtual const float* GetReflT0s(size_t Voxel) const override { return 0; }
    virtual float GetReflT0(size_t Voxel, size_t OpChannel) const override { return 0; }

    virtual int NOpChannels() const override { return fRecords.size(); }
    virtual int NVoxels() const override;

  protected:
    const sim::PhotonVoxelDef& fVoxDef;

    struct FitFunc {
      FitFunc() {}
      FitFunc(float n, float d) : norm(n), decay(d) {}
      double Eval(double x) const { return exp(norm + decay * x) / (x * x); }

      float norm, decay;
    };

    struct Exception {
      Exception(size_t _vox, float _vis) : vox(_vox), vis(_vis) {}

      bool operator<(const Exception& e) const { return vox < e.vox; }
      bool operator<(size_t v) const { return vox < v; }

      size_t vox;
      float vis;
    };

    struct OpDetRecord {
      FitFunc fit;
      // std::map is particularly space-inefficient, and sorted vector can also
      // be searched in log(N) time.
      std::vector<Exception> exceptions;
    };

    std::vector<OpDetRecord> fRecords;
  };
} // namespace

#endif