phot::PhotonLibraryHybrid Class Reference

TODO doc. More...

#include <PhotonLibraryHybrid.h>

Inheritance diagram for phot::PhotonLibraryHybrid:

Classes
struct	Exception
struct	FitFunc
struct	OpDetRecord

Public Member Functions
	PhotonLibraryHybrid (const std::string &fname, const sim::PhotonVoxelDef &voxdef)
virtual	~PhotonLibraryHybrid ()
virtual float	GetCount (size_t Voxel, size_t OpChannel) const override
virtual const float *	GetCounts (size_t Voxel) const override
	Returns a pointer to NOpChannels() visibility values, one per channel. More...
virtual bool	hasReflected () const override
	Don't implement reflected light. More...
virtual const float *	GetReflCounts (size_t Voxel) const override
virtual float	GetReflCount (size_t Voxel, size_t OpChannel) const override
virtual bool	hasReflectedT0 () const override
	Don't implement reflected light timing. More...
virtual const float *	GetReflT0s (size_t Voxel) const override
virtual float	GetReflT0 (size_t Voxel, size_t OpChannel) const override
virtual int	NOpChannels () const override
virtual int	NVoxels () const override
Public Member Functions inherited from phot::IPhotonLibrary
virtual	~IPhotonLibrary ()=default
virtual bool	isVoxelValid (size_t Voxel) const
size_t	LibrarySize () const
	Returns the number of elements in the library. More...

Protected Attributes
const sim::PhotonVoxelDef &	fVoxDef
std::vector< OpDetRecord >	fRecords

Additional Inherited Members
Public Types inherited from phot::IPhotonLibrary
using	Counts_t = const float *
	Type for visibility count per optical channel. More...
using	T0s_t = const float *
	Type for time of arrival per optical channel. More...
using	Params_t = std::vector< float > const *
using	Functions_t = TF1 *

Detailed Description

TODO doc.

Definition at line 18 of file PhotonLibraryHybrid.h.

Constructor & Destructor Documentation

phot::PhotonLibraryHybrid::PhotonLibraryHybrid	(	const std::string &	fname,
		const sim::PhotonVoxelDef &	voxdef
	)

Definition at line 27 of file PhotonLibraryHybrid.cxx.

    : fVoxDef(voxdef)
  {
    TFile f(fname.c_str());
    !f.IsZombie() || fatal("Could not open PhotonLibrary "+fname);

    for(int opdetIdx = 0; true; ++opdetIdx){
      const std::string dirname = TString::Format("opdet_%d", opdetIdx).Data();
      TDirectory* dir = (TDirectory*)f.Get(dirname.c_str());
      if(!dir) break; // Ran out of opdets

      OpDetRecord rec;

      TVectorD* fit = (TVectorD*)dir->Get("fit");
      fit || fatal("Didn't find "+dirname+"/fit in "+fname);
      rec.fit = FitFunc((*fit)[0], (*fit)[1]);

      TTree* tr = (TTree*)dir->Get("tr");
      tr || fatal("Didn't find "+dirname+"/tr in "+fname);
      int vox;
      float vis;
      tr->SetBranchAddress("vox", &vox);
      tr->SetBranchAddress("vis", &vis);

      rec.exceptions.reserve(tr->GetEntries());
      for(int i = 0; i < tr->GetEntries(); ++i){
        tr->GetEntry(i);
        vox < NVoxels() || fatal("Voxel out of range");
        rec.exceptions.emplace_back(vox, vis);
      }
      rec.exceptions.shrink_to_fit(); // In case we don't trust reserve()

      // In case they weren't sorted in the file
      std::sort(rec.exceptions.begin(), rec.exceptions.end());

      fRecords.push_back(rec);
    } // end for opdetIdx

    !fRecords.empty() || fatal("No opdet_*/ directories in "+fname);

    art::ServiceHandle<geo::Geometry const> geom;
    geom->NOpDets() == fRecords.size() || fatal("Number of opdets mismatch");

    fRecords.shrink_to_fit(); // save memory
  }

phot::PhotonLibraryHybrid::~PhotonLibraryHybrid ( )

virtual

Definition at line 75 of file PhotonLibraryHybrid.cxx.

  {
  }

Member Function Documentation

float phot::PhotonLibraryHybrid::GetCount ( size_t  Voxel, size_t  OpChannel  ) const

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 98 of file PhotonLibraryHybrid.cxx.

  {
    int(vox) < NVoxels() || fatal("GetCount(): Voxel out of range");
    int(opchan) < NOpChannels() || fatal("GetCount(): OpChan out of range");

    const OpDetRecord& rec = fRecords[opchan];

    auto it2 = std::lower_bound(rec.exceptions.begin(),
                                rec.exceptions.end(),
                                vox);
    if(it2->vox == vox) return it2->vis;

    // Otherwise, we use the interpolation, which requires a distance

    static art::ServiceHandle<geo::Geometry const> geom;
    const geo::OpDetGeo& opdet = geom->OpDetGeoFromOpDet(opchan);

    const auto voxvec = fVoxDef.GetPhotonVoxel(vox).GetCenter();

    const double dist = opdet.DistanceToPoint(voxvec);

    return rec.fit.Eval(dist);
  }

const float * phot::PhotonLibraryHybrid::GetCounts ( size_t  Voxel ) const

overridevirtual

Returns a pointer to NOpChannels() visibility values, one per channel.

Implements phot::IPhotonLibrary.

Definition at line 86 of file PhotonLibraryHybrid.cxx.

  {
    // The concept of GetCounts() conflicts fairly badly with how this library
    // works. This is probably the best we can do.

    static float* counts = 0;
    if(!counts) counts = new float[NOpChannels()];
    for(int od = 0; od < NOpChannels(); ++od) counts[od] = GetCount(vox, od);
    return counts;
  }

virtual float phot::PhotonLibraryHybrid::GetReflCount ( size_t  Voxel, size_t  OpChannel  ) const

inlineoverridevirtual

Implements phot::IPhotonLibrary.

Definition at line 33 of file PhotonLibraryHybrid.h.

{return 0;}

virtual const float* phot::PhotonLibraryHybrid::GetReflCounts ( size_t  Voxel ) const

inlineoverridevirtual

Implements phot::IPhotonLibrary.

Definition at line 32 of file PhotonLibraryHybrid.h.

{return 0;}

virtual float phot::PhotonLibraryHybrid::GetReflT0 ( size_t  Voxel, size_t  OpChannel  ) const

inlineoverridevirtual

Implements phot::IPhotonLibrary.

Definition at line 38 of file PhotonLibraryHybrid.h.

{return 0;}

virtual const float* phot::PhotonLibraryHybrid::GetReflT0s ( size_t  Voxel ) const

inlineoverridevirtual

Implements phot::IPhotonLibrary.

Definition at line 37 of file PhotonLibraryHybrid.h.

{return 0;}

virtual bool phot::PhotonLibraryHybrid::hasReflected ( ) const

inlineoverridevirtual

Don't implement reflected light.

Implements phot::IPhotonLibrary.

Definition at line 31 of file PhotonLibraryHybrid.h.

{return false;}

virtual bool phot::PhotonLibraryHybrid::hasReflectedT0 ( ) const

inlineoverridevirtual

Don't implement reflected light timing.

Implements phot::IPhotonLibrary.

Definition at line 36 of file PhotonLibraryHybrid.h.

{return false;}

virtual int phot::PhotonLibraryHybrid::NOpChannels ( ) const

inlineoverridevirtual

Implements phot::IPhotonLibrary.

Definition at line 40 of file PhotonLibraryHybrid.h.

{return fRecords.size();}

int phot::PhotonLibraryHybrid::NVoxels ( ) const

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 80 of file PhotonLibraryHybrid.cxx.

  {
    return fVoxDef.GetNVoxels();
  }

Member Data Documentation

std::vector<OpDetRecord> phot::PhotonLibraryHybrid::fRecords

protected

Definition at line 74 of file PhotonLibraryHybrid.h.

const sim::PhotonVoxelDef& phot::PhotonLibraryHybrid::fVoxDef

protected

Definition at line 44 of file PhotonLibraryHybrid.h.

---

The documentation for this class was generated from the following files:

• PhotonLibraryHybrid.h
• PhotonLibraryHybrid.cxx