#include <PhotonLibrary.h>

Inheritance diagram for phot::PhotonLibrary:

Public Member Functions
	PhotonLibrary (art::TFileDirectory *pDir=nullptr)
	If no valid `pDir` is provided, storage features will not be supported. More...

TTree *	ProduceTTree () const

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

void	SetCount (size_t Voxel, size_t OpChannel, float Count)

float	GetTimingPar (size_t Voxel, size_t OpChannel, size_t parnum) const

void	SetTimingPar (size_t Voxel, size_t OpChannel, float Count, size_t parnum)

void	SetTimingTF1 (size_t Voxel, size_t OpChannel, TF1 func)

virtual float	GetReflCount (size_t Voxel, size_t OpChannel) const override

void	SetReflCount (size_t Voxel, size_t OpChannel, float Count)

virtual float	GetReflT0 (size_t Voxel, size_t OpChannel) const override

void	SetReflT0 (size_t Voxel, size_t OpChannel, float reflT0)

virtual float const *	GetCounts (size_t Voxel) const override
	Returns a pointer to NOpChannels() visibility values, one per channel. More...

const std::vector< float > *	GetTimingPars (size_t Voxel) const

TF1 *	GetTimingTF1s (size_t Voxel) const

virtual float const *	GetReflCounts (size_t Voxel) const override

virtual float const *	GetReflT0s (size_t Voxel) const override

bool	hasTiming () const
	Returns whether the current library deals with time propagation distributions. More...

virtual bool	hasReflected () const override
	Returns whether the current library deals with reflected light count. More...

virtual bool	hasReflectedT0 () const override
	Returns whether the current library deals with reflected light timing. More...

void	StoreLibraryToFile (std::string LibraryFile, bool storeReflected=false, bool storeReflT0=false, size_t storeTiming=0) const

void	LoadLibraryFromFile (std::string LibraryFile, size_t NVoxels, bool storeReflected=false, bool storeReflT0=false, size_t storeTiming=0, int maxrange=200)

void	CreateEmptyLibrary (size_t NVoxels, size_t NChannels, bool storeReflected=false, bool storeReflT0=false, size_t storeTiming=0)

	PhotonLibrary ()

	~PhotonLibrary ()

float	GetCount (size_t Voxel, size_t OpChannel)

void	SetCount (size_t Voxel, size_t OpChannel, float Count)

const std::vector< float > *	GetCounts (size_t Voxel) const
	Returns a pointer to NOpChannels() visibility values, one per channel. More...

const std::vector< std::vector < float > > &	GetData () const

void	StoreLibraryToFile (std::string LibraryFile)

void	LoadLibraryFromFile (std::string LibraryFile, size_t NVoxels)

void	CreateEmptyLibrary (size_t NVoxels, size_t NChannels)

int	NOpChannels () const

int	NVoxels () const

Metadata: voxel information
bool	hasVoxelDef () const
	Returns whether voxel metadata is available. More...

sim::PhotonVoxelDef const &	GetVoxelDef () const

void	SetVoxelDef (sim::PhotonVoxelDef const &voxelDef)

Public Member Functions inherited from phot::IPhotonLibrary
virtual	~IPhotonLibrary ()=default

size_t	LibrarySize () const
	Returns the number of elements in the library. More...

Private Member Functions
virtual int	NOpChannels () const override

virtual int	NVoxels () const override

virtual bool	isVoxelValid (size_t Voxel) const override

bool	isVoxelValidImpl (size_t Voxel) const

size_t	uncheckedIndex (size_t Voxel, size_t OpChannel) const
	Returns the index of visibility of specified voxel and cell. More...

float	uncheckedAccess (size_t Voxel, size_t OpChannel) const
	Unchecked access to a visibility datum. More...

float &	uncheckedAccess (size_t Voxel, size_t OpChannel)
	Unchecked access to a visibility datum. More...

float	uncheckedAccessRefl (size_t Voxel, size_t OpChannel) const
	Unchecked access to a reflected visibility datum. More...

float &	uncheckedAccessRefl (size_t Voxel, size_t OpChannel)
	Unchecked access to a reflected visibility datum. More...

float	uncheckedAccessReflT (size_t Voxel, size_t OpChannel) const
	Unchecked access to a reflected T0 visibility datum. More...

float &	uncheckedAccessReflT (size_t Voxel, size_t OpChannel)
	Unchecked access to a reflected T0 visibility datum. More...

float	uncheckedAccessTimingPar (size_t Voxel, size_t OpChannel, size_t parnum) const
	Unchecked access to a parameter the time distribution. More...

float &	uncheckedAccessTimingPar (size_t Voxel, size_t OpChannel, size_t parnum)
	Unchecked access to a parameter of the time distribution. More...

TF1 &	uncheckedAccessTimingTF1 (size_t Voxel, size_t OpChannel)
	Unchecked access to a parameter of the time distribution. More...

const TF1 &	uncheckedAccessTimingTF1 (size_t Voxel, size_t OpChannel) const
	Unchecked access to a parameter of the time distribution. More...

void	LoadMetadata (TDirectory &srcDir)
	Reads the metadata from specified ROOT directory and sets it as current. More...

void	StoreMetadata () const
	Writes the current metadata (if any) into the ROOT output file. More...

Static Private Member Functions
static size_t	ExtractNOpChannels (TTree *tree)
	Returns the number of optical channels in the specified tree. More...

static int	size_t2int (size_t val)
	Converts size_t into integer. More...

static std::vector< float > *	EmptyList ()

Private Attributes
bool	fHasReflected = false
	Whether the current library deals with reflected light counts. More...

bool	fHasReflectedT0
	Whether the current library deals with reflected light timing. More...

size_t	fHasTiming
	Whether the current library deals with time propagation distribution. More...

util::LazyVector< float >	fLookupTable

util::LazyVector< float >	fReflLookupTable

util::LazyVector< float >	fReflTLookupTable

util::LazyVector< std::vector < float > >	fTimingParLookupTable

util::LazyVector< TF1 >	fTimingParTF1LookupTable

std::string	fTimingParFormula

size_t	fTimingParNParameters

size_t	fNOpChannels

size_t	fNVoxels

std::optional < sim::PhotonVoxelDef >	fVoxelDef
	Voxel definition loaded from library metadata. More...

art::TFileDirectory *	fDir = nullptr
	ROOT directory where to write data. More...

std::vector< std::vector< float > >	fLookupTable

Static Private Attributes
static std::string const	OpChannelBranchName = "OpChannel"
	Name of the optical channel number in the input tree. More...

static std::vector< float >	EmptyChannelsList

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

Definition at line 25 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

Constructor & Destructor Documentation

phot::PhotonLibrary::PhotonLibrary ( art::TFileDirectory * pDir = nullptr )

If no valid pDir is provided, storage features will not be supported.

Definition at line 46 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

46 : fDir(pDir) {}

phot::PhotonLibrary::fDir

art::TFileDirectory * fDir

ROOT directory where to write data.

Definition: larsim/larsim/PhotonPropagation/PhotonLibrary.h:165

phot::PhotonLibrary::PhotonLibrary ( )

Definition at line 23 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.cxx.

   {
     fLookupTable.clear();
   }

phot::PhotonLibrary::~PhotonLibrary ( )

Definition at line 31 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.cxx.

   {
     fLookupTable.clear();
   }

Member Function Documentation

void phot::PhotonLibrary::CreateEmptyLibrary	(	size_t	NVoxels,
		size_t	NChannels
	)

Definition at line 77 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.cxx.

   {
     fLookupTable.clear();
 
     fNVoxels     = NVoxels;
     fNOpChannels = NOpChannels;
 
     fLookupTable.resize(NVoxels);
 
     for(size_t ivox=0; ivox!=NVoxels; ivox++)
       {
         fLookupTable[ivox].resize(NOpChannels,0);
       }
   }

void phot::PhotonLibrary::CreateEmptyLibrary	(	size_t	NVoxels,
		size_t	NChannels,
		bool	storeReflected = `false`,
		bool	storeReflT0 = `false`,
		size_t	storeTiming = `0`
	)

Definition at line 134 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     fLookupTable.clear();
     fReflLookupTable.clear();
     fReflTLookupTable.clear();
     fTimingParLookupTable.clear();
     fTimingParTF1LookupTable.clear();
 
     fNVoxels = NVoxels;
     fNOpChannels = NOpChannels;
 
     fLookupTable.resize(LibrarySize());
     fHasReflected = storeReflected;
     if (storeReflected) fReflLookupTable.resize(LibrarySize());
     fHasReflectedT0 = storeReflT0;
     if (storeReflT0) fReflTLookupTable.resize(LibrarySize());
     fHasTiming = storeTiming;
     if (storeTiming != 0) {
       fTimingParLookupTable.resize(LibrarySize());
       fTimingParTF1LookupTable.resize(LibrarySize());
     }
   }

static std::vector<float>* phot::PhotonLibrary::EmptyList ( )

inlinestaticprivate

Definition at line 45 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h.

46 { EmptyChannelsList.clear(); return &EmptyChannelsList; }

phot::PhotonLibrary::EmptyChannelsList

static std::vector< float > EmptyChannelsList

Definition: sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h:43

size_t phot::PhotonLibrary::ExtractNOpChannels ( TTree * tree )

staticprivate

Returns the number of optical channels in the specified tree.

Definition at line 577 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     TBranch* channelBranch = tree->GetBranch(OpChannelBranchName.c_str());
     if (!channelBranch) {
       throw art::Exception(art::errors::NotFound)
         << "Tree '" << tree->GetName() << "' has no branch 'OpChannel'";
     }
 
     // fix a new local address for the branch
     char* oldAddress = channelBranch->GetAddress();
     Int_t channel;
     channelBranch->SetAddress(&channel);
     Int_t maxChannel = -1;
 
     // read all the channel values and pick the largest one
     Long64_t iEntry = 0;
     while (channelBranch->GetEntry(iEntry++)) {
       if (channel > maxChannel) maxChannel = channel;
     } // while
 
     MF_LOG_DEBUG("PhotonLibrary") << "Detected highest channel to be " << maxChannel << " from "
                                   << iEntry << " tree entries";
 
     // restore the old branch address
     channelBranch->SetAddress(oldAddress);
 
     return size_t(maxChannel + 1);
 
   } // PhotonLibrary::ExtractNOpChannels()

float phot::PhotonLibrary::GetCount	(	size_t	Voxel,
		size_t	OpChannel
	)

Definition at line 186 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.cxx.

   {
     //if(/*(Voxel<0)||*/(Voxel>=fNVoxels)||/*(OpChannel<0)||*/(OpChannel>=fNOpChannels))
     //  return 0;
     //else
       return fLookupTable[Voxel][OpChannel]; 
   }

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

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 310 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       return 0;
     else
       return uncheckedAccess(Voxel, OpChannel);
   }

const std::vector<float>* phot::PhotonLibrary::GetCounts ( size_t Voxel ) const

virtual

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

Implements phot::IPhotonLibrary.

const std::vector< float > * phot::PhotonLibrary::GetCounts ( size_t Voxel ) const

overridevirtual

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

Implements phot::IPhotonLibrary.

Definition at line 406 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if (Voxel >= fNVoxels)
       return nullptr;
     else
       return fLookupTable.data_address(uncheckedIndex(Voxel, 0));
   }

const std::vector<std::vector<float> >& phot::PhotonLibrary::GetData ( ) const

inline

Definition at line 26 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h.

27 { return fLookupTable; }

phot::PhotonLibrary::fLookupTable

util::LazyVector< float > fLookupTable

Definition: larsim/larsim/PhotonPropagation/PhotonLibrary.h:150

float phot::PhotonLibrary::GetReflCount	(	size_t	Voxel,
		size_t	OpChannel
	)		const

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 329 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       return 0;
     else
       return uncheckedAccessRefl(Voxel, OpChannel);
   }

float const * phot::PhotonLibrary::GetReflCounts ( size_t Voxel ) const

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 457 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if (Voxel >= fNVoxels)
       return nullptr;
     else
       return fReflLookupTable.data_address(uncheckedIndex(Voxel, 0));
   }

float phot::PhotonLibrary::GetReflT0	(	size_t	Voxel,
		size_t	OpChannel
	)		const

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 339 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       return 0;
     else
       return uncheckedAccessReflT(Voxel, OpChannel);
   }

float const * phot::PhotonLibrary::GetReflT0s ( size_t Voxel ) const

overridevirtual

Implements phot::IPhotonLibrary.

Definition at line 468 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if (Voxel >= fNVoxels)
       return nullptr;
     else
       return fReflTLookupTable.data_address(uncheckedIndex(Voxel, 0));
   }

float phot::PhotonLibrary::GetTimingPar	(	size_t	Voxel,
		size_t	OpChannel,
		size_t	parnum
	)		const

Definition at line 320 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       return 0;
     else
       return uncheckedAccessTimingPar(Voxel, OpChannel, parnum);
   } //----------------------------------------------------

const std::vector< float > * phot::PhotonLibrary::GetTimingPars ( size_t Voxel ) const

Definition at line 417 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if (Voxel >= fNVoxels)
       return nullptr;
     else
       return fTimingParLookupTable.data_address(uncheckedIndex(Voxel, 0));
   }

TF1 * phot::PhotonLibrary::GetTimingTF1s ( size_t Voxel ) const

Definition at line 428 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if (Voxel >= fNVoxels) return nullptr;
     /*
      * Sorry, Universe, but we can't undergo ROOT's bad design hell.
      * TF1::GetRandom() is non-const member, because it uses some internal
      * integral information which can be produced on the spot instead than
      * always be present. That's called caching, it's Good, but it must not
      * interfere with constantness of the interface (in fact, this is one of
      * the few acceptable uses of `mutable` members).
      * Because of this, this method can't return a constant `TF1`, therefore
      * it can't be constant, therefore the underlying access returning a
      * constant object is not acceptable.
      * So I do the Bad thing.
      * Plus I opened JIRA ROOT-9549
      * (https://sft.its.cern.ch/jira/browse/ROOT-9549).
      * After that is solved, this method should become:
      *
      * TF1 const* PhotonLibrary::GetTimingTF1s(size_t Voxel) const
      *
      * and the users should update their code accordingly.
      */
     else
       return const_cast<TF1*>(fTimingParTF1LookupTable.data_address(uncheckedIndex(Voxel, 0)));
   }

sim::PhotonVoxelDef const& phot::PhotonLibrary::GetVoxelDef ( ) const

inline

Returns the current voxel metadata (undefined behaviour if none).

See Also: hasVoxelDef()

Definition at line 106 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       assert(fVoxelDef);
       return *fVoxelDef;
     }

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

inlineoverridevirtual

Returns whether the current library deals with reflected light count.

Implements phot::IPhotonLibrary.

Definition at line 64 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fHasReflected;
     }

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

inlineoverridevirtual

Returns whether the current library deals with reflected light timing.

Implements phot::IPhotonLibrary.

Definition at line 71 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fHasReflectedT0;
     }

bool phot::PhotonLibrary::hasTiming ( ) const

inline

Returns whether the current library deals with time propagation distributions.

Definition at line 57 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fHasTiming;
     }

bool phot::PhotonLibrary::hasVoxelDef ( ) const

inline

Returns whether voxel metadata is available.

Definition at line 98 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fVoxelDef.has_value();
     }

virtual bool phot::PhotonLibrary::isVoxelValid ( size_t Voxel ) const

inlineoverrideprivatevirtual

Reimplemented from phot::IPhotonLibrary.

Definition at line 136 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return isVoxelValidImpl(Voxel);
     }

bool phot::PhotonLibrary::isVoxelValidImpl ( size_t Voxel ) const

inlineprivate

Definition at line 168 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return Voxel < fNVoxels;
     }

void phot::PhotonLibrary::LoadLibraryFromFile	(	std::string	LibraryFile,
		size_t	NVoxels
	)

Definition at line 95 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.cxx.

   {
     fLookupTable.clear();
 
     std::cout<< "Reading photon library from input file: " << LibraryFile.c_str()<<std::endl;
 
     TFile *f = nullptr;
     TTree *tt = nullptr;
 
     try
       {
         f  =  TFile::Open(LibraryFile.c_str());
         if(!f) {
           std::cerr<<"\033[95m<<"<<__FUNCTION__<<">>\033[00m " << "Failed to open a ROOT file: " << LibraryFile.c_str()<<std::endl;
           std::cerr<<"If you don't have photon library data file, download from below URL..."<<std::endl;
           std::cerr<<"/cvmfs/icarus.opensciencegrid.org/products/icarus/icarus_data/v08_28_00/icarus_data/PhotonLibrary/PhotonLibrary-20180801.root"<<std::endl<<std::endl;
           throw std::exception();
         }
         tt =  (TTree*)f->Get("PhotonLibraryData");
         if (!tt) { // Library not in the top directory
             TKey *key = f->FindKeyAny("PhotonLibraryData");
             if (key)
                 tt = (TTree*)key->ReadObj();
             else {
               std::cerr << "PhotonLibraryData not found in file" <<LibraryFile<<std::endl;
             }
         }
       }
     catch(...)
       {
         std::cerr << "Error in ttree load, reading photon library: " << LibraryFile.c_str()<<std::endl;
       }
 
     Int_t     Voxel;
     Int_t     OpChannel;
     Float_t   Visibility;
 
     tt->SetBranchAddress("Voxel",      &Voxel);
     tt->SetBranchAddress("OpChannel",  &OpChannel);
     tt->SetBranchAddress("Visibility", &Visibility);
 
 
 
 
     fNVoxels     = NVoxels;
     fNOpChannels = 1;      // Minimum default, overwritten by library reading
 
 
     fLookupTable.resize(NVoxels);
 
 
     size_t NEntries = tt->GetEntries();
 
     for(size_t i=0; i!=NEntries; ++i) {
       tt->GetEntry(i);
 
       // Set # of optical channels to 1 more than largest one seen
       if (OpChannel >= (int)fNOpChannels)
         fNOpChannels = OpChannel+1;
 
       // Expand this voxel's vector if needed
       if (fLookupTable[Voxel].size() < fNOpChannels)
         fLookupTable[Voxel].resize(fNOpChannels, 0);
 
       // Set the visibility at this optical channel
       fLookupTable[Voxel][OpChannel] = Visibility;
     }
 
     // Go through the table and fill in any missing 0's
     for(size_t ivox=0; ivox!=NVoxels; ivox++)
     {
       if (fLookupTable[ivox].size() < fNOpChannels)
         fLookupTable[ivox].resize(fNOpChannels,0);
     }
 
 
     std::cout <<  NVoxels << " voxels,  " << fNOpChannels<<" channels" <<std::endl;
 
 
     try
       {
         f->Close();
       }
     catch(...)
       {
         std::cerr << "Error in closing file : " << LibraryFile.c_str()<<std::endl;
       }
   }

void phot::PhotonLibrary::LoadLibraryFromFile	(	std::string	LibraryFile,
		size_t	NVoxels,
		bool	storeReflected = `false`,
		bool	storeReflT0 = `false`,
		size_t	storeTiming = `0`,
		int	maxrange = `200`
	)

Definition at line 165 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     fLookupTable.clear();
     fReflLookupTable.clear();
     fReflTLookupTable.clear();
     fTimingParLookupTable.clear();
     fTimingParTF1LookupTable.clear();
 
     mf::LogInfo("PhotonLibrary") << "Reading photon library from input file: "
                                  << LibraryFile.c_str() << std::endl;
 
     TFile* f = nullptr;
     TTree* tt = nullptr;
     TDirectory* pSrcDir = nullptr;
 
     try {
       f = TFile::Open(LibraryFile.c_str());
       tt = (TTree*)f->Get("PhotonLibraryData");
       if (tt) { pSrcDir = f; }
       else { // Library not in the top directory
         TKey* key = f->FindKeyAny("PhotonLibraryData");
         if (key) {
           tt = (TTree*)key->ReadObj();
           pSrcDir = key->GetMotherDir();
         }
         else {
           mf::LogError("PhotonLibrary") << "PhotonLibraryData not found in file" << LibraryFile;
         }
       }
     }
     catch (...) {
       throw cet::exception("PhotonLibrary")
         << "Error in ttree load, reading photon library: " << LibraryFile << "\n";
     }
 
     Int_t Voxel;
     Int_t OpChannel;
     Float_t Visibility;
     Float_t ReflVisibility;
     Float_t ReflTfirst;
     std::vector<Float_t> timing_par;
 
     tt->SetBranchAddress("Voxel", &Voxel);
     tt->SetBranchAddress("OpChannel", &OpChannel);
     tt->SetBranchAddress("Visibility", &Visibility);
 
     fHasTiming = getTiming;
 
     fHasReflected = getReflected;
     if (getReflected) tt->SetBranchAddress("ReflVisibility", &ReflVisibility);
     fHasReflectedT0 = getReflT0;
     if (getReflT0) tt->SetBranchAddress("ReflTfirst", &ReflTfirst);
 
     fNVoxels = NVoxels;
     fNOpChannels = PhotonLibrary::ExtractNOpChannels(tt); // EXPENSIVE!!!
 
     // with STL vectors, where `resize()` directly controls the allocation of
     // memory, reserving the space is redundant; not so with `util::LazyVector`,
     // where `resize()` never increases the memory; `data_init()` allocates
     // all the storage we need at once, effectively suppressing the laziness
     // of the vector (by design, that was only relevant in `CreateEmptyLibrary()`)
     fLookupTable.resize(LibrarySize());
     fLookupTable.data_init(LibrarySize());
 
     if (fHasTiming != 0) {
       timing_par.resize(getTiming);
       tt->SetBranchAddress("timing_par", timing_par.data());
       fTimingParNParameters = fHasTiming;
       // should be pSrcDir->Get()? kept as is for backward compatibility
       TNamed* n = (TNamed*)f->Get("fTimingParFormula");
       if (!n)
         mf::LogError("PhotonLibrary")
           << "Error reading the photon propagation formula. Please check the photon library."
           << std::endl;
       fTimingParFormula = n->GetTitle();
       fTimingParTF1LookupTable.resize(LibrarySize());
       fTimingParTF1LookupTable.data_init(LibrarySize());
       mf::LogInfo("PhotonLibrary")
         << "Time parametrization is activated. Using the formula: " << fTimingParFormula << " with "
         << fTimingParNParameters << " parameters." << std::endl;
     }
     if (fHasReflected) {
       fReflLookupTable.resize(LibrarySize());
       fReflLookupTable.data_init(LibrarySize());
     }
     if (fHasReflectedT0) {
       fReflTLookupTable.resize(LibrarySize());
       fReflTLookupTable.data_init(LibrarySize());
     }
 
     size_t NEntries = tt->GetEntries();
 
     for (size_t i = 0; i != NEntries; ++i) {
 
       tt->GetEntry(i);
 
       // Set the visibility at this optical channel
       uncheckedAccess(Voxel, OpChannel) = Visibility;
 
       if (fHasReflected) uncheckedAccessRefl(Voxel, OpChannel) = ReflVisibility;
       if (fHasReflectedT0) uncheckedAccessReflT(Voxel, OpChannel) = ReflTfirst;
       if (fHasTiming != 0) {
         // TODO: use TF1::Copy
         TF1 timingfunction(Form("timing_%i_%i", Voxel, OpChannel),
                            fTimingParFormula.c_str(),
                            timing_par[0],
                            fTimingMaxRange);
         timingfunction.SetParameter(
           0,
           0.0); //first parameter is now in the range. Let's do this to keep compatible with old libraries.
         for (size_t k = 1; k < fTimingParNParameters; k++) {
           timingfunction.SetParameter(k, timing_par[k]);
         }
 
         uncheckedAccessTimingTF1(Voxel, OpChannel) = timingfunction;
       }
     } // for entries
 
     LoadMetadata(*pSrcDir);
     {
       mf::LogInfo log("PhotonLibrary");
       log << "Photon lookup table size : " << NVoxels << " voxels,  " << fNOpChannels
           << " channels";
       if (hasVoxelDef())
         log << "; " << GetVoxelDef();
       else
         log << " (no voxel geometry included)";
     }
 
     try {
       f->Close();
     }
     catch (...) {
       mf::LogError("PhotonLibrary") << "Error in closing file : " << LibraryFile;
     }
   }

void phot::PhotonLibrary::LoadMetadata ( TDirectory & srcDir )

private

Reads the metadata from specified ROOT directory and sets it as current.

Definition at line 478 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
 
     constexpr std::size_t NExpectedKeys = 9U;
 
     std::vector<std::string> missingKeys;
 
     RooReader<RooInt, Int_t> readInt{srcDir, missingKeys};
     RooReader<RooDouble, Double_t> readDouble{srcDir, missingKeys};
 
     double xMin;
     double xMax;
     int xN;
     double yMin;
     double yMax;
     int yN;
     double zMin;
     double zMax;
     int zN;
     if (auto metaValue = readDouble("MinX")) xMin = *metaValue;
     if (auto metaValue = readDouble("MaxX")) xMax = *metaValue;
     if (auto metaValue = readInt("NDivX")) xN = *metaValue;
     if (auto metaValue = readDouble("MinY")) yMin = *metaValue;
     if (auto metaValue = readDouble("MaxY")) yMax = *metaValue;
     if (auto metaValue = readInt("NDivY")) yN = *metaValue;
     if (auto metaValue = readDouble("MinZ")) zMin = *metaValue;
     if (auto metaValue = readDouble("MaxZ")) zMax = *metaValue;
     if (auto metaValue = readInt("NDivZ")) zN = *metaValue;
 
     if (!missingKeys.empty()) {
       if (missingKeys.size() != NExpectedKeys) {
         mf::LogError log("PhotonLibrary");
         log << "Photon library at '" << srcDir.GetPath() << "' is missing " << missingKeys.size()
             << " metadata elements:";
         for (auto const& key : missingKeys)
           log << " '" << key << "'";
       }
       else {
         mf::LogTrace("PhotonLibrary") << "No voxel metadata found in '" << srcDir.GetPath() << "'";
       }
       return;
     } // if missing keys
 
     fVoxelDef.emplace(xMin, xMax, xN, yMin, yMax, yN, zMin, zMax, zN);
 
   } // PhotonLibrary::LoadMetadata()

int phot::PhotonLibrary::NOpChannels ( ) const

inlinevirtual

Implements phot::IPhotonLibrary.

Definition at line 34 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h.

34 { return fNOpChannels; }

phot::PhotonLibrary::fNOpChannels

size_t fNOpChannels

Definition: larsim/larsim/PhotonPropagation/PhotonLibrary.h:158

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

inlineoverrideprivatevirtual

Implements phot::IPhotonLibrary.

Definition at line 125 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fNOpChannels;
     }

int phot::PhotonLibrary::NVoxels ( ) const

inlinevirtual

Implements phot::IPhotonLibrary.

Definition at line 35 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h.

35 { return fNVoxels; }

phot::PhotonLibrary::fNVoxels

size_t fNVoxels

Definition: larsim/larsim/PhotonPropagation/PhotonLibrary.h:159

virtual int phot::PhotonLibrary::NVoxels ( ) const

inlineoverrideprivatevirtual

Implements phot::IPhotonLibrary.

Definition at line 130 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fNVoxels;
     }

TTree* phot::PhotonLibrary::ProduceTTree ( ) const

void phot::PhotonLibrary::SetCount	(	size_t	Voxel,
		size_t	OpChannel,
		float	Count
	)

void phot::PhotonLibrary::SetCount	(	size_t	Voxel,
		size_t	OpChannel,
		float	Count
	)

Definition at line 350 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       mf::LogError("PhotonLibrary")
         << "Error - attempting to set count in voxel " << Voxel << " which is out of range";
     else
       uncheckedAccess(Voxel, OpChannel) = Count;
   }

void phot::PhotonLibrary::SetReflCount	(	size_t	Voxel,
		size_t	OpChannel,
		float	Count
	)

Definition at line 383 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       mf::LogError("PhotonLibrary")
         << "Error - attempting to set count in voxel " << Voxel << " which is out of range";
     else
       uncheckedAccessRefl(Voxel, OpChannel) = Count;
   }

void phot::PhotonLibrary::SetReflT0	(	size_t	Voxel,
		size_t	OpChannel,
		float	reflT0
	)

Definition at line 394 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       mf::LogError("PhotonLibrary")
         << "Error - attempting to set count in voxel " << Voxel << " which is out of range";
     else
       uncheckedAccessReflT(Voxel, OpChannel) = Count;
   }

void phot::PhotonLibrary::SetTimingPar	(	size_t	Voxel,
		size_t	OpChannel,
		float	Count,
		size_t	parnum
	)

Definition at line 361 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       mf::LogError("PhotonLibrary") << "Error - attempting to set timing t0 count in voxel "
                                     << Voxel << " which is out of range";
     else
       uncheckedAccessTimingPar(Voxel, OpChannel, parnum) = Count;
   }

void phot::PhotonLibrary::SetTimingTF1	(	size_t	Voxel,
		size_t	OpChannel,
		TF1	func
	)

Definition at line 372 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     if ((Voxel >= fNVoxels) || (OpChannel >= fNOpChannels))
       mf::LogError("PhotonLibrary") << "Error - attempting to set a propagation function in voxel "
                                     << Voxel << " which is out of range";
     else
       uncheckedAccessTimingTF1(Voxel, OpChannel) = func;
   }

void phot::PhotonLibrary::SetVoxelDef ( sim::PhotonVoxelDef const & voxelDef )

inline

Copies the specified voxel definition into our own (overwrites the existing metadata if any).

Definition at line 115 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       fVoxelDef = voxelDef;
     }

static int phot::PhotonLibrary::size_t2int ( size_t val )

inlinestaticprivate

Converts size_t into integer.

Definition at line 265 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       constexpr size_t int_max_as_size_t{std::numeric_limits<int>::max()};
       return val <= int_max_as_size_t ? static_cast<int>(val) : -1;
     }

void phot::PhotonLibrary::StoreLibraryToFile ( std::string LibraryFile )

Definition at line 38 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.cxx.

   {
     std::cout << "Writing photon library to input file: " << LibraryFile.c_str()<<std::endl;
 
     TFile fout(LibraryFile.c_str(),"RECREATE");
 
     TTree *tt = new TTree("PhotonLibraryData","PhotonLibraryData");
 
     Int_t     Voxel;
     Int_t     OpChannel;
     Float_t   Visibility;
 
 
     tt->Branch("Voxel",      &Voxel,      "Voxel/I");
     tt->Branch("OpChannel",  &OpChannel,  "OpChannel/I");
     tt->Branch("Visibility", &Visibility, "Visibility/F");
 
 
     for(size_t ivox=0; ivox!=fLookupTable.size(); ++ivox)
       {
         for(size_t ichan=0; ichan!=fLookupTable.at(ivox).size(); ++ichan)
           {
             if(fLookupTable[ivox].at(ichan) > 0)
               {
                 Voxel      = ivox;
                 OpChannel  = ichan;
                 Visibility = fLookupTable[ivox][ichan];
                 tt->Fill();
               }
           }
       }
     fout.cd();
     tt->Write();
     fout.Close();
   }

void phot::PhotonLibrary::StoreLibraryToFile	(	std::string	LibraryFile,
		bool	storeReflected = `false`,
		bool	storeReflT0 = `false`,
		size_t	storeTiming = `0`
	)		const

Definition at line 51 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
     mf::LogInfo("PhotonLibrary") << "Writing photon library to file";
 
     if (!fDir) {
       throw cet::exception("PhotonLibrary")
         << "StoreLibraryToFile(): no ROOT file provided, can't store anything.\n";
     }
 
     TTree* tt = fDir->make<TTree>("PhotonLibraryData", "PhotonLibraryData");
 
     Int_t Voxel = 0;
     Int_t OpChannel = 0;
     Float_t Visibility = 0;
     Float_t ReflVisibility = 0;
     Float_t ReflTfirst = 0;
     Float_t* timing_par = nullptr;
 
     tt->Branch("Voxel", &Voxel, "Voxel/I");
     tt->Branch(OpChannelBranchName.c_str(), &OpChannel, (OpChannelBranchName + "/I").c_str());
     tt->Branch("Visibility", &Visibility, "Visibility/F");
 
     if (storeTiming != 0) {
       if (!hasTiming()) {
         // if this happens, you need to call CreateEmptyLibrary() with storeReflected set true
         throw cet::exception("PhotonLibrary")
           << "StoreLibraryToFile() requested to store the time propagation distribution "
              "parameters, which was not simulated.";
       }
       tt->Branch("timing_par", timing_par, Form("timing_par[%i]/F", size_t2int(storeTiming)));
       if (fLookupTable.size() != fTimingParLookupTable.size())
         throw cet::exception(" Photon Library ")
           << "Time propagation lookup table is different size than Direct table \n"
           << "this should not be happening. ";
     }
 
     if (storeReflected) {
       if (!hasReflected()) {
         // if this happens, you need to call CreateEmptyLibrary() with storeReflected set true
         throw cet::exception("PhotonLibrary")
           << "StoreLibraryToFile() requested to store reflected light, which was not simulated.";
       }
       tt->Branch("ReflVisibility", &ReflVisibility, "ReflVisibility/F");
       if (fLookupTable.size() != fReflLookupTable.size())
         throw cet::exception(" Photon Library ")
           << "Reflected light lookup table is different size than Direct table \n"
           << "this should not be happening. ";
     }
     if (storeReflT0) {
       if (!hasReflectedT0()) {
         // if this happens, you need to call CreateEmptyLibrary() with storeReflectedT0 set true
         throw cet::exception("PhotonLibrary") << "StoreLibraryToFile() requested to store "
                                                  "reflected light timing, which was not simulated.";
       }
       tt->Branch("ReflTfirst", &ReflTfirst, "ReflTfirst/F");
     }
     for (size_t ivox = 0; ivox != fNVoxels; ++ivox) {
       for (size_t ichan = 0; ichan != fNOpChannels; ++ichan) {
         Visibility = uncheckedAccess(ivox, ichan);
         if (storeReflected) ReflVisibility = uncheckedAccessRefl(ivox, ichan);
         if (storeReflT0) ReflTfirst = uncheckedAccessReflT(ivox, ichan);
         if (storeTiming != 0) {
           for (size_t i = 0; i < storeTiming; i++)
             timing_par[i] = uncheckedAccessTimingPar(ivox, ichan, i);
         }
         if (Visibility > 0 || ReflVisibility > 0) {
           Voxel = ivox;
           OpChannel = ichan;
           // visibility(ies) is(are) already set
           tt->Fill();
         }
       }
     }
 
     StoreMetadata();
   }

void phot::PhotonLibrary::StoreMetadata ( ) const

private

Writes the current metadata (if any) into the ROOT output file.

Definition at line 527 of file larsim/larsim/PhotonPropagation/PhotonLibrary.cxx.

   {
 
     assert(fDir);
 
     // NVoxels
     fDir->makeAndRegister<RooInt>("NVoxels", "Total number of voxels in the library", fNVoxels);
 
     // NChannels
     fDir->makeAndRegister<RooInt>(
       "NChannels", "Total number of optical detector channels in the library", fNOpChannels);
 
     if (!hasVoxelDef()) return;
     sim::PhotonVoxelDef const& voxelDef = GetVoxelDef();
 
     // lower point
     geo::Point_t const& lower = voxelDef.GetRegionLowerCorner();
     fDir->makeAndRegister<RooDouble>(
       "MinX", "Lower x coordinate covered by the library (world coordinates, cm)", lower.X());
     fDir->makeAndRegister<RooDouble>(
       "MinY", "Lower y coordinate covered by the library (world coordinates, cm)", lower.Y());
     fDir->makeAndRegister<RooDouble>(
       "MinZ", "Lower z coordinate covered by the library (world coordinates, cm)", lower.Z());
 
     // upper point
     geo::Point_t const& upper = voxelDef.GetRegionUpperCorner();
     fDir->makeAndRegister<RooDouble>(
       "MaxX", "Upper x coordinate covered by the library (world coordinates, cm)", upper.X());
     fDir->makeAndRegister<RooDouble>(
       "MaxY", "Upper y coordinate covered by the library (world coordinates, cm)", upper.Y());
     fDir->makeAndRegister<RooDouble>(
       "MaxZ", "Upper z coordinate covered by the library (world coordinates, cm)", upper.Z());
 
     // steps
     geo::Vector_t const& stepSizes = voxelDef.GetVoxelSize();
     fDir->makeAndRegister<RooDouble>("StepX", "Size on x direction of a voxel (cm)", stepSizes.X());
     fDir->makeAndRegister<RooDouble>("StepY", "Size on y direction of a voxel (cm)", stepSizes.Y());
     fDir->makeAndRegister<RooDouble>("StepZ", "Size on z direction of a voxel (cm)", stepSizes.Z());
 
     // divisions
     auto const& steps = voxelDef.GetSteps();
     fDir->makeAndRegister<RooInt>("NDivX", "Steps on the x direction", steps[0]);
     fDir->makeAndRegister<RooInt>("NDivY", "Steps on the y direction", steps[1]);
     fDir->makeAndRegister<RooInt>("NDivZ", "Steps on the z direction", steps[2]);
 
   } // PhotonLibrary::StoreMetadata()

float phot::PhotonLibrary::uncheckedAccess	(	size_t	Voxel,
		size_t	OpChannel
	)		const

inlineprivate

Unchecked access to a visibility datum.

Definition at line 182 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fLookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

float& phot::PhotonLibrary::uncheckedAccess	(	size_t	Voxel,
		size_t	OpChannel
	)

inlineprivate

Unchecked access to a visibility datum.

Definition at line 189 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fLookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

float phot::PhotonLibrary::uncheckedAccessRefl	(	size_t	Voxel,
		size_t	OpChannel
	)		const

inlineprivate

Unchecked access to a reflected visibility datum.

Definition at line 196 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fReflLookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

float& phot::PhotonLibrary::uncheckedAccessRefl	(	size_t	Voxel,
		size_t	OpChannel
	)

inlineprivate

Unchecked access to a reflected visibility datum.

Definition at line 203 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fReflLookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

float phot::PhotonLibrary::uncheckedAccessReflT	(	size_t	Voxel,
		size_t	OpChannel
	)		const

inlineprivate

Unchecked access to a reflected T0 visibility datum.

Definition at line 210 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fReflTLookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

float& phot::PhotonLibrary::uncheckedAccessReflT	(	size_t	Voxel,
		size_t	OpChannel
	)

inlineprivate

Unchecked access to a reflected T0 visibility datum.

Definition at line 217 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fReflTLookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

float phot::PhotonLibrary::uncheckedAccessTimingPar	(	size_t	Voxel,
		size_t	OpChannel,
		size_t	parnum
	)		const

inlineprivate

Unchecked access to a parameter the time distribution.

Definition at line 224 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fTimingParLookupTable[uncheckedIndex(Voxel, OpChannel)][parnum];
     }

float& phot::PhotonLibrary::uncheckedAccessTimingPar	(	size_t	Voxel,
		size_t	OpChannel,
		size_t	parnum
	)

inlineprivate

Unchecked access to a parameter of the time distribution.

Definition at line 231 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fTimingParLookupTable[uncheckedIndex(Voxel, OpChannel)][parnum];
     }

TF1& phot::PhotonLibrary::uncheckedAccessTimingTF1	(	size_t	Voxel,
		size_t	OpChannel
	)

inlineprivate

Unchecked access to a parameter of the time distribution.

Definition at line 238 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return fTimingParTF1LookupTable[uncheckedIndex(Voxel, OpChannel)];
     }

const TF1& phot::PhotonLibrary::uncheckedAccessTimingTF1	(	size_t	Voxel,
		size_t	OpChannel
	)		const

inlineprivate

Unchecked access to a parameter of the time distribution.

Definition at line 245 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       // note that this will produce a segmentation fault if the formula is not there
       return *(fTimingParTF1LookupTable.data_address(uncheckedIndex(Voxel, OpChannel)));
     }

size_t phot::PhotonLibrary::uncheckedIndex	(	size_t	Voxel,
		size_t	OpChannel
	)		const

inlineprivate

Returns the index of visibility of specified voxel and cell.

Definition at line 175 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

     {
       return Voxel * fNOpChannels + OpChannel;
     }

Member Data Documentation

std::vector< float > phot::PhotonLibrary::EmptyChannelsList

staticprivate

Definition at line 43 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h.

art::TFileDirectory* phot::PhotonLibrary::fDir = nullptr

private

ROOT directory where to write data.

Definition at line 165 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

bool phot::PhotonLibrary::fHasReflected = false

private

Whether the current library deals with reflected light counts.

Definition at line 141 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

bool phot::PhotonLibrary::fHasReflectedT0

private

Initial value:

=

false

Whether the current library deals with reflected light timing.

Definition at line 142 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

size_t phot::PhotonLibrary::fHasTiming

private

Initial value:

=

0

Whether the current library deals with time propagation distribution.

Definition at line 145 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

std::vector<std::vector<float> > phot::PhotonLibrary::fLookupTable

private

Definition at line 39 of file sbncode/sbncode/OpT0Finder/flashmatch/Base/FMWKTools/PhotonLibrary.h.

util::LazyVector<float> phot::PhotonLibrary::fLookupTable

private

Definition at line 150 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

size_t phot::PhotonLibrary::fNOpChannels

private

Definition at line 158 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

size_t phot::PhotonLibrary::fNVoxels

private

Definition at line 159 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

util::LazyVector<float> phot::PhotonLibrary::fReflLookupTable

private

Definition at line 151 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

util::LazyVector<float> phot::PhotonLibrary::fReflTLookupTable

private

Definition at line 152 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

std::string phot::PhotonLibrary::fTimingParFormula

private

Definition at line 155 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

util::LazyVector<std::vector<float> > phot::PhotonLibrary::fTimingParLookupTable

private

Definition at line 153 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

size_t phot::PhotonLibrary::fTimingParNParameters

private

Definition at line 156 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

util::LazyVector<TF1> phot::PhotonLibrary::fTimingParTF1LookupTable

private

Definition at line 154 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

std::optional<sim::PhotonVoxelDef> phot::PhotonLibrary::fVoxelDef

private

Voxel definition loaded from library metadata.

Definition at line 162 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

std::string const phot::PhotonLibrary::OpChannelBranchName = "OpChannel"

staticprivate

Name of the optical channel number in the input tree.

Definition at line 258 of file larsim/larsim/PhotonPropagation/PhotonLibrary.h.

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

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation