#include "PhotonVisibilityService.h"

Public Member Functions
	PhotonVisibilityService (fhicl::ParameterSet const &pset, art::ActivityRegistry &reg)

void	reconfigure (fhicl::ParameterSet const &p)

double	GetQuenchingFactor (double dQdx) const

float	GetVisibility (double const *xyz, unsigned int OpChannel, bool wantReflected=false) const

float const *	GetAllVisibilities (double const *xyz, bool wantReflected=false) const

void	LoadLibrary () const

void	StoreLibrary ()

void	StoreLightProd (int VoxID, double N)

void	RetrieveLightProd (int &VoxID, double &N) const

void	SetLibraryEntry (int VoxID, int OpChannel, float N, bool wantReflected=false)

float	GetLibraryEntry (int VoxID, int OpChannel, bool wantReflected=false) const

float const *	GetLibraryEntries (int VoxID, bool wantReflected=false) const

float const *	GetReflT0s (double const *xyz) const

void	SetLibraryReflT0Entry (int VoxID, int OpChannel, float value)

float const *	GetLibraryReflT0Entries (int VoxID) const

float	GetLibraryReflT0Entry (int VoxID, int Channel) const

const std::vector< float > *	GetTimingPar (double const *xyz) const

void	SetLibraryTimingParEntry (int VoxID, int OpChannel, float value, size_t parnum)

const std::vector< float > *	GetLibraryTimingParEntries (int VoxID) const

float	GetLibraryTimingParEntry (int VoxID, int Channel, size_t npar) const

TF1 *	GetTimingTF1 (double const *xyz) const

void	SetLibraryTimingTF1Entry (int VoxID, int OpChannel, TF1 func)

TF1 *	GetLibraryTimingTF1Entries (int VoxID) const

void	SetDirectLightPropFunctions (TF1 const *functions[8], double &d_break, double &d_max, double &tf1_sampling_factor) const

void	SetReflectedCOLightPropFunctions (TF1 const *functions[5], double &t0_max, double &t0_break_point) const

bool	IsBuildJob () const

bool	UseParameterization () const

bool	StoreReflected () const

bool	StoreReflT0 () const

bool	IncludeParPropTime () const

size_t	ParPropTimeNpar () const

std::string	ParPropTimeFormula () const

bool	IncludePropTime () const

const sim::PhotonVoxelDef &	GetVoxelDef () const

size_t	NOpChannels () const

Static Public Member Functions
static double	DistanceToOpDet (double const *xyz, unsigned int OpDet)

static double	SolidAngleFactor (double const *xyz, unsigned int OpDet)

Private Attributes
int	fCurrentVoxel

double	fCurrentValue

float	fXmin

float	fXmax

float	fYmin

float	fYmax

float	fZmin

float	fZmax

int	fNx

int	fNy

int	fNz

bool	fUseCryoBoundary

bool	fLibraryBuildJob

bool	fDoNotLoadLibrary

bool	fParameterization

bool	fHybrid

bool	fStoreReflected

bool	fStoreReflT0

bool	fIncludePropTime

bool	fParPropTime

size_t	fParPropTime_npar

std::string	fParPropTime_formula

bool	fInterpolate

TF1 *	fparslogNorm

TF1 *	fparslogNorm_far

TF1 *	fparsMPV

TF1 *	fparsMPV_far

TF1 *	fparsWidth

TF1 *	fparsCte

TF1 *	fparsCte_far

TF1 *	fparsSlope

double	fD_break

double	fD_max

double	fTF1_sampling_factor

TF1 *	fparslogNorm_refl

TF1 *	fparsMPV_refl

TF1 *	fparsWidth_refl

TF1 *	fparsCte_refl

TF1 *	fparsSlope_refl

double	fT0_max

double	fT0_break_point

std::string	fLibraryFile

IPhotonLibrary *	fTheLibrary

sim::PhotonVoxelDef	fVoxelDef

Detailed Description

Definition at line 25 of file PhotonVisibilityService.h.

Constructor & Destructor Documentation

phot::PhotonVisibilityService::PhotonVisibilityService	(	fhicl::ParameterSet const &	pset,
		art::ActivityRegistry &	reg
	)

Definition at line 43 of file PhotonVisibilityService_service.cc.

References reconfigure().

                                                                                                         :
     
     fCurrentVoxel(0),
     fCurrentValue(0.),
     fXmin(0.),
     fXmax(0.),
     fYmin(0.),
     fYmax(0.),
     fZmin(0.),
     fZmax(0.),
     fNx(0),
     fNy(0),
     fNz(0),
     fUseCryoBoundary(false),
     fLibraryBuildJob(false),
     fDoNotLoadLibrary(false),
     fParameterization(false),
     fHybrid(false),
     fStoreReflected(false),
     fStoreReflT0(false),
     fIncludePropTime(false),
     fParPropTime(false),
     fTheLibrary(0)
     
   {
     this->reconfigure(pset);
     mf::LogInfo("PhotonVisibilityService")<<"PhotonVisbilityService initializing"<<std::endl;
   }

Member Function Documentation

double phot::PhotonVisibilityService::DistanceToOpDet	(	double const *	xyz,
		unsigned int	OpDet
	)

static

Definition at line 327 of file PhotonVisibilityService_service.cc.

References geo::OpDetGeo::DistanceToPoint(), and geo::GeometryCore::OpDetGeoFromOpDet().

   {
     art::ServiceHandle<geo::Geometry> geom;
     return geom->OpDetGeoFromOpDet(OpDet).DistanceToPoint(xyz);
       
   }

float const * phot::PhotonVisibilityService::GetAllVisibilities	(	double const *	xyz,
		bool	wantReflected = `false`
	)		const

Definition at line 309 of file PhotonVisibilityService_service.cc.

References fInterpolate, fVoxelDef, GetLibraryEntries(), GetVisibility(), sim::PhotonVoxelDef::GetVoxelID(), and NOpChannels().

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::AddLightFromSegment(), opdet::FlashHypothesisCreator::CreateFlashHypothesesFromSegment(), larg4::OpFastScintillation::RecordPhotonsProduced(), and opdet::TrackTimeAssoc::ScanMIPHypotheses().

   {
     if(fInterpolate){
       static std::vector<float> ret;
       if(ret.size() != NOpChannels()) ret.resize(NOpChannels());
       for(unsigned int i = 0; i < NOpChannels(); ++i) ret[i] = GetVisibility(xyz, i, wantReflected);
       return &ret.front();
     }
     else{
       size_t VoxID = fVoxelDef.GetVoxelID(xyz);
       return GetLibraryEntries(VoxID, wantReflected);
     }
   }

float const * phot::PhotonVisibilityService::GetLibraryEntries	(	int	VoxID,
		bool	wantReflected = `false`
	)		const

Definition at line 413 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::IPhotonLibrary::GetCounts(), phot::IPhotonLibrary::GetReflCounts(), and LoadLibrary().

Referenced by GetAllVisibilities().

   {
     if(fTheLibrary == 0)
       LoadLibrary();
 
     if(!wantReflected)
       return fTheLibrary->GetCounts(VoxID);
     else
       return fTheLibrary->GetReflCounts(VoxID);
   }

float phot::PhotonVisibilityService::GetLibraryEntry	(	int	VoxID,
		int	OpChannel,
		bool	wantReflected = `false`
	)		const

Definition at line 426 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::IPhotonLibrary::GetCount(), phot::IPhotonLibrary::GetReflCount(), and LoadLibrary().

Referenced by GetVisibility().

   {
     if(fTheLibrary == 0)
       LoadLibrary();
 
     if(!wantReflected)
       return fTheLibrary->GetCount(VoxID, Channel);
     else
       return fTheLibrary->GetReflCount(VoxID, Channel);
   }

float const * phot::PhotonVisibilityService::GetLibraryReflT0Entries ( int VoxID ) const

Definition at line 450 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::IPhotonLibrary::GetReflT0s(), and LoadLibrary().

Referenced by GetReflT0s().

   {
     if(fTheLibrary == 0)
       LoadLibrary();
 
     return fTheLibrary->GetReflT0s(VoxID);
   }

float phot::PhotonVisibilityService::GetLibraryReflT0Entry	(	int	VoxID,
		int	Channel
	)		const

Definition at line 473 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::IPhotonLibrary::GetReflT0(), and LoadLibrary().

   {
     if(fTheLibrary == 0)
       LoadLibrary();
 
     return fTheLibrary->GetReflT0(VoxID, Channel);
   }

const std::vector< float > * phot::PhotonVisibilityService::GetLibraryTimingParEntries ( int VoxID ) const

Definition at line 501 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::PhotonLibrary::GetTimingPars(), and LoadLibrary().

Referenced by GetTimingPar().

   {
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
     if(fTheLibrary == 0)
       LoadLibrary();
 
     return lib->GetTimingPars(VoxID);
   }

float phot::PhotonVisibilityService::GetLibraryTimingParEntry	(	int	VoxID,
		int	Channel,
		size_t	npar
	)		const

Definition at line 550 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::PhotonLibrary::GetTimingPar(), and LoadLibrary().

   {
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
     if(fTheLibrary == 0)
       LoadLibrary();
 
     return lib->GetTimingPar(VoxID, Channel,npar);
   }

TF1 * phot::PhotonVisibilityService::GetLibraryTimingTF1Entries ( int VoxID ) const

Definition at line 512 of file PhotonVisibilityService_service.cc.

References fTheLibrary, phot::PhotonLibrary::GetTimingTF1s(), and LoadLibrary().

Referenced by GetTimingTF1().

   {
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
     if(fTheLibrary == 0)
       LoadLibrary();
 
     return lib->GetTimingTF1s(VoxID);
   }

double phot::PhotonVisibilityService::GetQuenchingFactor ( double dQdx ) const

Definition at line 296 of file PhotonVisibilityService_service.cc.

   {
     // for now, no quenching
     return 1.0;
 
   }

float const * phot::PhotonVisibilityService::GetReflT0s ( double const * xyz ) const

Definition at line 442 of file PhotonVisibilityService_service.cc.

References fVoxelDef, GetLibraryReflT0Entries(), and sim::PhotonVoxelDef::GetVoxelID().

Referenced by larg4::OpFastScintillation::RecordPhotonsProduced().

   {
     int VoxID = fVoxelDef.GetVoxelID(xyz);
     return GetLibraryReflT0Entries(VoxID);
   }

const std::vector< float > * phot::PhotonVisibilityService::GetTimingPar ( double const * xyz ) const

Definition at line 486 of file PhotonVisibilityService_service.cc.

References fVoxelDef, GetLibraryTimingParEntries(), and sim::PhotonVoxelDef::GetVoxelID().

   {
     int VoxID = fVoxelDef.GetVoxelID(xyz);
     return GetLibraryTimingParEntries(VoxID);
   }

TF1 * phot::PhotonVisibilityService::GetTimingTF1 ( double const * xyz ) const

Definition at line 492 of file PhotonVisibilityService_service.cc.

References fVoxelDef, GetLibraryTimingTF1Entries(), and sim::PhotonVoxelDef::GetVoxelID().

Referenced by larg4::OpFastScintillation::RecordPhotonsProduced().

   {
     int VoxID = fVoxelDef.GetVoxelID(xyz);
     return GetLibraryTimingTF1Entries(VoxID);
   }

float phot::PhotonVisibilityService::GetVisibility	(	double const *	xyz,
		unsigned int	OpChannel,
		bool	wantReflected = `false`
	)		const

Definition at line 347 of file PhotonVisibilityService_service.cc.

References fInterpolate, fVoxelDef, GetLibraryEntry(), sim::PhotonVoxelDef::GetNeighboringVoxelIDs(), sim::PhotonVoxelDef::GetVoxelID(), and n.

Referenced by phot::CreateHybridLibrary::CreateHybridLibrary(), and GetAllVisibilities().

   {
     // Static to avoid reallocating this buffer between calls
     // (GetNeighbouringVoxelIDs makes sure to clear it).
     static std::vector<sim::PhotonVoxelDef::NeiInfo> neis;
 
     if(fInterpolate){
       // In case we're outside the bounding box we'll get an empty vector here
       // and return visibility 0, which seems OK.
       fVoxelDef.GetNeighboringVoxelIDs(xyz, neis);
     }
     else{
       // For no interpolation, use a single entry with weight 1
       neis.clear();
       neis.emplace_back(fVoxelDef.GetVoxelID(xyz), 1);
     }
 
     // Sum up all the weighted neighbours to get interpolation behaviour
     float vis = 0;
     for(const sim::PhotonVoxelDef::NeiInfo& n: neis){
       vis += n.weight * GetLibraryEntry(n.id, OpChannel, wantReflected);
     }
     return vis;
   }

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

inline

Definition at line 78 of file PhotonVisibilityService.h.

References fVoxelDef, and NOpChannels().

Referenced by phot::CreateHybridLibrary::CreateHybridLibrary(), evgen::LightSource::LightSource(), and LoadLibrary().

78 {return fVoxelDef; }

phot::PhotonVisibilityService::fVoxelDef

sim::PhotonVoxelDef fVoxelDef

Definition: PhotonVisibilityService.h:127

bool phot::PhotonVisibilityService::IncludeParPropTime ( ) const

inline

Definition at line 72 of file PhotonVisibilityService.h.

References fParPropTime.

Referenced by larg4::OpFastScintillation::RecordPhotonsProduced().

72 { return fParPropTime; }

phot::PhotonVisibilityService::fParPropTime

bool fParPropTime

Definition: PhotonVisibilityService.h:103

bool phot::PhotonVisibilityService::IncludePropTime ( ) const

inline

Definition at line 76 of file PhotonVisibilityService.h.

References fIncludePropTime.

Referenced by larg4::OpFastScintillation::OpFastScintillation(), and larg4::OpFastScintillation::RecordPhotonsProduced().

76 { return fIncludePropTime; }

phot::PhotonVisibilityService::fIncludePropTime

bool fIncludePropTime

Definition: PhotonVisibilityService.h:101

bool phot::PhotonVisibilityService::IsBuildJob ( ) const

inline

Definition at line 68 of file PhotonVisibilityService.h.

References fLibraryBuildJob.

Referenced by opdet::SimPhotonCounter::analyze(), opdet::SimPhotonCounter::endJob(), and evgen::LightSource::produce().

68 { return fLibraryBuildJob; }

phot::PhotonVisibilityService::fLibraryBuildJob

bool fLibraryBuildJob

Definition: PhotonVisibilityService.h:95

void phot::PhotonVisibilityService::LoadLibrary ( ) const

Definition at line 73 of file PhotonVisibilityService_service.cc.

References phot::PhotonLibrary::CreateEmptyLibrary(), fDoNotLoadLibrary, fHybrid, fLibraryBuildJob, fLibraryFile, fParameterization, fParPropTime_npar, fStoreReflected, fStoreReflT0, fTheLibrary, sim::PhotonVoxelDef::GetNVoxels(), GetVoxelDef(), phot::PhotonLibrary::LoadLibraryFromFile(), and geo::GeometryCore::NOpDets().

Referenced by GetLibraryEntries(), GetLibraryEntry(), GetLibraryReflT0Entries(), GetLibraryReflT0Entry(), GetLibraryTimingParEntries(), GetLibraryTimingParEntry(), GetLibraryTimingTF1Entries(), NOpChannels(), SetLibraryEntry(), SetLibraryReflT0Entry(), SetLibraryTimingParEntry(), SetLibraryTimingTF1Entry(), and StoreLibrary().

   {
     // Don't do anything if the library has already been loaded.
 
     if(fTheLibrary == 0) {
 
       if((!fLibraryBuildJob)&&(!fDoNotLoadLibrary)) {
         std::string LibraryFileWithPath;
         cet::search_path sp("FW_SEARCH_PATH");
 
         if( !sp.find_file(fLibraryFile, LibraryFileWithPath) )
           throw cet::exception("PhotonVisibilityService") << "Unable to find photon library in "  << sp.to_string() << "\n";
 
         if(!fParameterization) {
           art::ServiceHandle<geo::Geometry> geom;
 
           mf::LogInfo("PhotonVisibilityService") << "PhotonVisibilityService Loading photon library from file "
                                                  << LibraryFileWithPath
                                                  << " for "
                                                  << GetVoxelDef().GetNVoxels()
                                                  << " voxels and "
                                                  << geom->NOpDets()
                                                  << " optical detectors."
                                                  << std::endl;
 
           if(fHybrid){
             fTheLibrary = new PhotonLibraryHybrid(LibraryFileWithPath,
                                                   GetVoxelDef());
           }
           else{
             PhotonLibrary* lib = new PhotonLibrary;
             fTheLibrary = lib;
 
             size_t NVoxels = GetVoxelDef().GetNVoxels();
             lib->LoadLibraryFromFile(LibraryFileWithPath, NVoxels, fStoreReflected, fStoreReflT0, fParPropTime_npar);
           }
         }
       }
       else {
         art::ServiceHandle<geo::Geometry> geom;
 
         size_t NOpDets = geom->NOpDets();
         size_t NVoxels = GetVoxelDef().GetNVoxels();
         mf::LogInfo("PhotonVisibilityService") << " Vis service running library build job.  Please ensure " 
                                                << " job contains LightSource, LArG4, SimPhotonCounter"<<std::endl;
         PhotonLibrary* lib = new PhotonLibrary;
         fTheLibrary = lib;
 
         lib->CreateEmptyLibrary(NVoxels, NOpDets, fStoreReflected, fStoreReflT0, fParPropTime_npar);
       }
 
     }
   }

size_t phot::PhotonVisibilityService::NOpChannels ( ) const

Definition at line 561 of file PhotonVisibilityService_service.cc.

References fTheLibrary, LoadLibrary(), and phot::IPhotonLibrary::NOpChannels().

Referenced by cosmic::BeamFlashTrackMatchTaggerAlg::AddLightFromSegment(), GetAllVisibilities(), GetVoxelDef(), larg4::OpFastScintillation::RecordPhotonsProduced(), and opdet::TrackTimeAssoc::ScanMIPHypotheses().

   {
     if(fTheLibrary == 0)
       LoadLibrary();
     
     return fTheLibrary->NOpChannels();
   }

std::string phot::PhotonVisibilityService::ParPropTimeFormula ( ) const

inline

Definition at line 74 of file PhotonVisibilityService.h.

References fParPropTime_formula.

74 { return fParPropTime_formula; }

phot::PhotonVisibilityService::fParPropTime_formula

std::string fParPropTime_formula

Definition: PhotonVisibilityService.h:105

size_t phot::PhotonVisibilityService::ParPropTimeNpar ( ) const

inline

Definition at line 73 of file PhotonVisibilityService.h.

References fParPropTime_npar.

73 { return fParPropTime_npar; }

phot::PhotonVisibilityService::fParPropTime_npar

size_t fParPropTime_npar

Definition: PhotonVisibilityService.h:104

void phot::PhotonVisibilityService::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 148 of file PhotonVisibilityService_service.cc.

References geo::GeometryCore::CryostatBoundaries(), fD_break, fD_max, fDoNotLoadLibrary, fHybrid, fIncludePropTime, fInterpolate, fLibraryBuildJob, fLibraryFile, fNx, fNy, fNz, fParameterization, fParPropTime, fParPropTime_formula, fParPropTime_npar, fparsCte, fparsCte_far, fparsCte_refl, fparslogNorm, fparslogNorm_far, fparslogNorm_refl, fparsMPV, fparsMPV_far, fparsMPV_refl, fparsSlope, fparsSlope_refl, fparsWidth, fparsWidth_refl, fStoreReflected, fStoreReflT0, fT0_break_point, fT0_max, fTF1_sampling_factor, fUseCryoBoundary, fVoxelDef, fXmax, fXmin, fYmax, fYmin, fZmax, fZmin, and fhicl::ParameterSet::get().

Referenced by PhotonVisibilityService().

   {
 
     art::ServiceHandle<geo::Geometry> geom;
     
     // Library details
     fLibraryBuildJob      = p.get< bool        >("LibraryBuildJob"     );
     fParameterization     = p.get< bool        >("DUNE10ktParameterization", false);
     fHybrid               = p.get< bool        >("HybridLibrary", false);
     fLibraryFile          = p.get< std::string >("LibraryFile"         );
     fDoNotLoadLibrary     = p.get< bool        >("DoNotLoadLibrary"    );
     fStoreReflected       = p.get< bool        >("StoreReflected", false);
     fStoreReflT0          = p.get< bool        >("StoreReflT0",  false);
     fIncludePropTime      = p.get< bool        >("IncludePropTime", false);
     // Voxel parameters
     fUseCryoBoundary      = p.get< bool        >("UseCryoBoundary"     );
     fInterpolate          = p.get< bool        >("Interpolate", false);
 
     fParPropTime          = p.get< bool        >("ParametrisedTimePropagation", false);
     fParPropTime_npar     = p.get< size_t      >("ParametrisedTimePropagationNParameters", 0);
     fParPropTime_formula  = p.get< std::string >("ParametrisedTimePropagationFittedFormula","");
     
     if (!fParPropTime)
     {
       fParPropTime_npar=0;
     }
 
     if(fUseCryoBoundary)
       {
         double CryoBounds[6];
         geom->CryostatBoundaries(CryoBounds);
         fXmin = CryoBounds[0];
         fXmax = CryoBounds[1];
         fYmin = CryoBounds[2];
         fYmax = CryoBounds[3];
         fZmin = CryoBounds[4];
         fZmax = CryoBounds[5];
       }
     else
       {
         fXmin      = p.get< double       >("XMin"     );
         fXmax      = p.get< double       >("XMax"     );
         fYmin      = p.get< double       >("YMin"     );
         fYmax      = p.get< double       >("YMax"     );
         fZmin      = p.get< double       >("ZMin"     );
         fZmax      = p.get< double       >("ZMax"     );
       }
 
     fNx        = p.get< int          >("NX"       );
     fNy        = p.get< int          >("NY"       );
     fNz        = p.get< int          >("NZ"       );
     
     fVoxelDef = sim::PhotonVoxelDef(fXmin, fXmax, fNx, fYmin, fYmax, fNy, fZmin, fZmax, fNz);
 
     if(fIncludePropTime)
       {
         // Construct parameterized model parameter functions.         
         std::cout<< "Getting direct light parameters from .fcl file"<<std::endl;
         std::vector<std::string> direct_functions = p.get<std::vector<std::string> >("Direct_functions");
         //range of distances where the parametrization is valid                                            
         fD_break = p.get<double>("D_break");
         fD_max = p.get<double>("D_max");
 
         fTF1_sampling_factor = p.get<double>("TF1_sampling_factor");    
 
         std::vector<double> direct_landauNormpars = p.get<std::vector<double> >("Direct_landauNormpars");
         fparslogNorm = new TF1("fparslogNorm", direct_functions[0].c_str(), 0., fD_break);
         for(unsigned int i=0; i<direct_landauNormpars.size(); ++i)
           fparslogNorm->SetParameter(i, direct_landauNormpars[i]);
         
         std::vector<double> direct_landauMPVpars = p.get<std::vector<double> >("Direct_landauMPVpars");
         fparsMPV = new TF1("fparsMPV", direct_functions[1].c_str(), 0., fD_break);
         for(unsigned int i=0; i<direct_landauMPVpars.size(); ++i)       
           fparsMPV->SetParameter(i, direct_landauMPVpars[i]);
 
         std::vector<double> direct_landauWidthpars = p.get<std::vector<double> >("Direct_landauWidthpars");
         fparsWidth = new TF1("fparsWidth", direct_functions[2].c_str(), 0., fD_break);
         for(unsigned int i=0; i<direct_landauWidthpars.size(); ++i)                
           fparsWidth->SetParameter(i, direct_landauWidthpars[i]);
 
         std::vector<double> direct_expoCtepars = p.get<std::vector<double> >("Direct_expoCtepars");
         fparsCte = new TF1("fparsCte", direct_functions[3].c_str(), 0., fD_break);
         for(unsigned int i=0; i<direct_expoCtepars.size(); ++i)  
           fparsCte->SetParameter(i, direct_expoCtepars[i]);
 
         std::vector<double> direct_expoSlopepars = p.get<std::vector<double> >("Direct_expoSlopepars");
         fparsSlope = new TF1("fparsSlope", direct_functions[4].c_str(), 0., fD_break);
         for(unsigned int i=0; i<direct_expoSlopepars.size(); ++i)
           fparsSlope->SetParameter(i, direct_expoSlopepars[i]);
 
         std::vector<double> direct_landauNormpars_far = p.get<std::vector<double> >("Direct_landauNormpars_far");
         fparslogNorm_far = new TF1("fparslogNorm_far", direct_functions[5].c_str(), fD_break, fD_max);
         for(unsigned int i=0; i<direct_landauNormpars_far.size(); ++i)
           fparslogNorm_far->SetParameter(i, direct_landauNormpars_far[i]);
 
         std::vector<double> direct_landauMPVpars_far = p.get<std::vector<double> >("Direct_landauMPVpars_far");
         fparsMPV_far = new TF1("fparsMPV_far", direct_functions[6].c_str(), fD_break, fD_max);
         for(unsigned int i=0; i<direct_landauMPVpars_far.size(); ++i)
           fparsMPV_far->SetParameter(i, direct_landauMPVpars_far[i]);
 
         std::vector<double> direct_expoCtepars_far = p.get<std::vector<double> >("Direct_expoCtepars_far");
         fparsCte_far = new TF1("fparsCte_far", direct_functions[7].c_str(), fD_break - 50., fD_max);
         for(unsigned int i=0; i<direct_expoCtepars_far.size(); ++i)
           fparsCte_far->SetParameter(i, direct_expoCtepars_far[i]);
 
         std::vector<std::string> reflected_functions = p.get<std::vector<std::string> >("Reflected_functions");
         //times where the parametrizations are valid or change
         fT0_max = p.get<double>("T0_max");
         fT0_break_point = p.get<double>("T0_break_point");
 
         std::vector<double> reflected_landauNormpars = p.get<std::vector<double> >("Reflected_landauNormpars");
         fparslogNorm_refl = new TF1("fparslogNorm_refl", reflected_functions[0].c_str(), 0., fT0_max);
         for(unsigned int i=0; i<reflected_landauNormpars.size(); ++i)
           fparslogNorm_refl->SetParameter(i, reflected_landauNormpars[i]);
 
         std::vector<double> reflected_landauMPVpars = p.get<std::vector<double> >("Reflected_landauMPVpars");
         fparsMPV_refl = new TF1("fparsMPV_refl", reflected_functions[1].c_str(), 0., fT0_max);
         for(unsigned int i=0; i<reflected_landauMPVpars.size(); ++i)
           fparsMPV_refl->SetParameter(i, reflected_landauMPVpars[i]);
 
         std::vector<double> reflected_landauWidthpars = p.get<std::vector<double> >("Reflected_landauWidthpars");
         fparsWidth_refl = new TF1("fparsWidth_refl", reflected_functions[2].c_str(), 0., fT0_max);
         for(unsigned int i=0; i<reflected_landauWidthpars.size(); ++i)
           fparsWidth_refl->SetParameter(i, reflected_landauWidthpars[i]);
 
         std::vector<double> reflected_expoCtepars = p.get<std::vector<double> >("Reflected_expoCtepars");
         fparsCte_refl = new TF1("fparsCte_refl", reflected_functions[3].c_str(), 0., fT0_max);
         for(unsigned int i=0; i<reflected_expoCtepars.size(); ++i)
           fparsCte_refl->SetParameter(i, reflected_expoCtepars[i]);
 
         std::vector<double> reflected_expoSlopepars = p.get<std::vector<double> >("Reflected_expoSlopepars");
         fparsSlope_refl = new TF1("fparsSlope_refl", reflected_functions[4].c_str(), 0., fT0_max);
         for(unsigned int i=0; i<reflected_expoSlopepars.size(); ++i)
           fparsSlope_refl->SetParameter(i, reflected_expoSlopepars[i]);
 
 
       }
 
 
     return;
         
   }

void phot::PhotonVisibilityService::RetrieveLightProd	(	int &	VoxID,
		double &	N
	)		const

Definition at line 386 of file PhotonVisibilityService_service.cc.

References fCurrentValue, and fCurrentVoxel.

Referenced by opdet::SimPhotonCounter::analyze().

   {
     N     = fCurrentValue;
     VoxID = fCurrentVoxel;
   }

void phot::PhotonVisibilityService::SetDirectLightPropFunctions	(	TF1 const *	functions[8],
		double &	d_break,
		double &	d_max,
		double &	tf1_sampling_factor
	)		const

Definition at line 570 of file PhotonVisibilityService_service.cc.

References fD_break, fD_max, fparsCte, fparsCte_far, fparslogNorm, fparslogNorm_far, fparsMPV, fparsMPV_far, fparsSlope, fparsWidth, and fTF1_sampling_factor.

Referenced by larg4::OpFastScintillation::OpFastScintillation().

   {
     functions[0] = fparslogNorm;
     functions[1] = fparsMPV;
     functions[2] = fparsWidth;
     functions[3] = fparsCte;
     functions[4] = fparsSlope; 
     functions[5] = fparslogNorm_far;
     functions[6] = fparsMPV_far;
     functions[7] = fparsCte_far;
 
     d_break = fD_break;
     d_max = fD_max;
     tf1_sampling_factor = fTF1_sampling_factor;
   }

void phot::PhotonVisibilityService::SetLibraryEntry	(	int	VoxID,
		int	OpChannel,
		float	N,
		bool	wantReflected = `false`
	)

Definition at line 394 of file PhotonVisibilityService_service.cc.

References fTheLibrary, LoadLibrary(), LOG_DEBUG, phot::PhotonLibrary::SetCount(), and phot::PhotonLibrary::SetReflCount().

Referenced by opdet::SimPhotonCounter::analyze().

   {   
     if(fTheLibrary == 0)
       LoadLibrary();
     
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
 
     if(!wantReflected)
     lib->SetCount(VoxID,OpChannel, N);
     
     else
     lib->SetReflCount(VoxID,OpChannel, N);
     
     //std::cout<< " PVS logging " << VoxID << " " << OpChannel<<std::endl;
     LOG_DEBUG("PhotonVisibilityService") << " PVS logging " << VoxID << " " << OpChannel<<std::endl;
   }

void phot::PhotonVisibilityService::SetLibraryReflT0Entry	(	int	VoxID,
		int	OpChannel,
		float	value
	)

Definition at line 460 of file PhotonVisibilityService_service.cc.

References fTheLibrary, LoadLibrary(), LOG_DEBUG, and phot::PhotonLibrary::SetReflT0().

Referenced by opdet::SimPhotonCounter::analyze().

   {
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
     if(fTheLibrary == 0)
       LoadLibrary();
 
     lib->SetReflT0(VoxID,OpChannel,T0);
 
     LOG_DEBUG("PhotonVisibilityService") << " PVS logging " << VoxID << " " << OpChannel<<std::endl;
   }

void phot::PhotonVisibilityService::SetLibraryTimingParEntry	(	int	VoxID,
		int	OpChannel,
		float	value,
		size_t	parnum
	)

Definition at line 523 of file PhotonVisibilityService_service.cc.

References fTheLibrary, LoadLibrary(), LOG_DEBUG, and phot::PhotonLibrary::SetTimingPar().

   {
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
     if(fTheLibrary == 0)
       LoadLibrary();
 
     lib->SetTimingPar(VoxID,OpChannel,par, parnum);
 
     LOG_DEBUG("PhotonVisibilityService") << " PVS logging " << VoxID << " " << OpChannel<<std::endl;
   }

void phot::PhotonVisibilityService::SetLibraryTimingTF1Entry	(	int	VoxID,
		int	OpChannel,
		TF1	func
	)

Definition at line 536 of file PhotonVisibilityService_service.cc.

References fTheLibrary, LoadLibrary(), LOG_DEBUG, and phot::PhotonLibrary::SetTimingTF1().

   {
     PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
     if(fTheLibrary == 0)
       LoadLibrary();
 
     lib->SetTimingTF1(VoxID,OpChannel,func);
 
     LOG_DEBUG("PhotonVisibilityService") << " PVS logging " << VoxID << " " << OpChannel<<std::endl;
   }

void phot::PhotonVisibilityService::SetReflectedCOLightPropFunctions	(	TF1 const *	functions[5],
		double &	t0_max,
		double &	t0_break_point
	)		const

Definition at line 587 of file PhotonVisibilityService_service.cc.

References DEFINE_ART_SERVICE, fparsCte_refl, fparslogNorm_refl, fparsMPV_refl, fparsSlope_refl, fparsWidth_refl, fT0_break_point, and fT0_max.

Referenced by larg4::OpFastScintillation::OpFastScintillation().

   {
     functions[0] = fparslogNorm_refl;
     functions[1] = fparsMPV_refl;
     functions[2] = fparsWidth_refl;
     functions[3] = fparsCte_refl;
     functions[4] = fparsSlope_refl;
 
     t0_max = fT0_max;
     t0_break_point = fT0_break_point;
   }

double phot::PhotonVisibilityService::SolidAngleFactor	(	double const *	xyz,
		unsigned int	OpDet
	)

static

Definition at line 339 of file PhotonVisibilityService_service.cc.

References geo::OpDetGeo::CosThetaFromNormal(), and geo::GeometryCore::OpDetGeoFromOpDet().

   {
     art::ServiceHandle<geo::Geometry> geom;
     return geom->OpDetGeoFromOpDet(OpDet).CosThetaFromNormal(xyz);
   }

void phot::PhotonVisibilityService::StoreLibrary ( )

Definition at line 128 of file PhotonVisibilityService_service.cc.

References fHybrid, fLibraryBuildJob, fLibraryFile, fParPropTime_npar, fStoreReflected, fStoreReflT0, fTheLibrary, LoadLibrary(), and phot::PhotonLibrary::StoreLibraryToFile().

Referenced by opdet::SimPhotonCounter::endJob().

   {
     if(fTheLibrary == 0)
       LoadLibrary();
 
     if(fLibraryBuildJob )
       {
           
           if(fHybrid){
               std::cout<< "This is would be building a Hybrid Library. Not defined. "<<std::endl;
           }
           mf::LogInfo("PhotonVisibilityService") << " Vis service "
                                                << " Storing Library entries to file..." <<std::endl;
           PhotonLibrary* lib = dynamic_cast<PhotonLibrary*>(fTheLibrary);
           lib->StoreLibraryToFile(fLibraryFile, fStoreReflected, fStoreReflT0, fParPropTime_npar);
       }
   }

void phot::PhotonVisibilityService::StoreLightProd	(	int	VoxID,
		double	N
	)

Definition at line 375 of file PhotonVisibilityService_service.cc.

References fCurrentValue, and fCurrentVoxel.

Referenced by evgen::LightSource::produce().

   {
     fCurrentVoxel = VoxID;
     fCurrentValue = N;
     mf::LogInfo("PhotonVisibilityService") << " PVS notes production of " << N << " photons at Vox " << VoxID<<std::endl; 
   }

bool phot::PhotonVisibilityService::StoreReflected ( ) const

inline

Definition at line 70 of file PhotonVisibilityService.h.

References fStoreReflected.

Referenced by opdet::SimPhotonCounter::analyze(), opdet::SimPhotonCounter::beginJob(), and larg4::OpFastScintillation::RecordPhotonsProduced().

70 { return fStoreReflected; }

phot::PhotonVisibilityService::fStoreReflected

bool fStoreReflected

Definition: PhotonVisibilityService.h:99

bool phot::PhotonVisibilityService::StoreReflT0 ( ) const

inline

Definition at line 71 of file PhotonVisibilityService.h.

References fStoreReflT0.

Referenced by opdet::SimPhotonCounter::analyze(), and larg4::OpFastScintillation::RecordPhotonsProduced().

71 { return fStoreReflT0; }

phot::PhotonVisibilityService::fStoreReflT0

bool fStoreReflT0

Definition: PhotonVisibilityService.h:100

bool phot::PhotonVisibilityService::UseParameterization ( ) const

inline

Definition at line 69 of file PhotonVisibilityService.h.

References fParameterization.

69 {return fParameterization;}

phot::PhotonVisibilityService::fParameterization

bool fParameterization

Definition: PhotonVisibilityService.h:97

Member Data Documentation

double phot::PhotonVisibilityService::fCurrentValue

private

Definition at line 84 of file PhotonVisibilityService.h.

Referenced by RetrieveLightProd(), and StoreLightProd().

int phot::PhotonVisibilityService::fCurrentVoxel

private

Definition at line 83 of file PhotonVisibilityService.h.

Referenced by RetrieveLightProd(), and StoreLightProd().

double phot::PhotonVisibilityService::fD_break

private

Definition at line 117 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

double phot::PhotonVisibilityService::fD_max

private

Definition at line 117 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

bool phot::PhotonVisibilityService::fDoNotLoadLibrary

private

Definition at line 96 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), and reconfigure().

bool phot::PhotonVisibilityService::fHybrid

private

Definition at line 98 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), reconfigure(), and StoreLibrary().

bool phot::PhotonVisibilityService::fIncludePropTime

private

Definition at line 101 of file PhotonVisibilityService.h.

Referenced by IncludePropTime(), and reconfigure().

bool phot::PhotonVisibilityService::fInterpolate

private

Definition at line 107 of file PhotonVisibilityService.h.

Referenced by GetAllVisibilities(), GetVisibility(), and reconfigure().

bool phot::PhotonVisibilityService::fLibraryBuildJob

private

Definition at line 95 of file PhotonVisibilityService.h.

Referenced by IsBuildJob(), LoadLibrary(), reconfigure(), and StoreLibrary().

std::string phot::PhotonVisibilityService::fLibraryFile

private

Definition at line 125 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), reconfigure(), and StoreLibrary().

int phot::PhotonVisibilityService::fNx

private

Definition at line 91 of file PhotonVisibilityService.h.

Referenced by reconfigure().

int phot::PhotonVisibilityService::fNy

private

Definition at line 91 of file PhotonVisibilityService.h.

Referenced by reconfigure().

int phot::PhotonVisibilityService::fNz

private

Definition at line 91 of file PhotonVisibilityService.h.

Referenced by reconfigure().

bool phot::PhotonVisibilityService::fParameterization

private

Definition at line 97 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), reconfigure(), and UseParameterization().

bool phot::PhotonVisibilityService::fParPropTime

private

Definition at line 103 of file PhotonVisibilityService.h.

Referenced by IncludeParPropTime(), and reconfigure().

std::string phot::PhotonVisibilityService::fParPropTime_formula

private

Definition at line 105 of file PhotonVisibilityService.h.

Referenced by ParPropTimeFormula(), and reconfigure().

size_t phot::PhotonVisibilityService::fParPropTime_npar

private

Definition at line 104 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), ParPropTimeNpar(), reconfigure(), and StoreLibrary().

TF1* phot::PhotonVisibilityService::fparsCte

private

Definition at line 114 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsCte_far

private

Definition at line 115 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsCte_refl

private

Definition at line 121 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparslogNorm

private

Definition at line 109 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparslogNorm_far

private

Definition at line 110 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparslogNorm_refl

private

Definition at line 118 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsMPV

private

Definition at line 111 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsMPV_far

private

Definition at line 112 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsMPV_refl

private

Definition at line 119 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsSlope

private

Definition at line 116 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsSlope_refl

private

Definition at line 122 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsWidth

private

Definition at line 113 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

TF1* phot::PhotonVisibilityService::fparsWidth_refl

private

Definition at line 120 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

bool phot::PhotonVisibilityService::fStoreReflected

private

Definition at line 99 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), reconfigure(), StoreLibrary(), and StoreReflected().

bool phot::PhotonVisibilityService::fStoreReflT0

private

Definition at line 100 of file PhotonVisibilityService.h.

Referenced by LoadLibrary(), reconfigure(), StoreLibrary(), and StoreReflT0().

double phot::PhotonVisibilityService::fT0_break_point

private

Definition at line 123 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

double phot::PhotonVisibilityService::fT0_max

private

Definition at line 123 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetReflectedCOLightPropFunctions().

double phot::PhotonVisibilityService::fTF1_sampling_factor

private

Definition at line 117 of file PhotonVisibilityService.h.

Referenced by reconfigure(), and SetDirectLightPropFunctions().

IPhotonLibrary* phot::PhotonVisibilityService::fTheLibrary

mutableprivate

Definition at line 126 of file PhotonVisibilityService.h.

Referenced by GetLibraryEntries(), GetLibraryEntry(), GetLibraryReflT0Entries(), GetLibraryReflT0Entry(), GetLibraryTimingParEntries(), GetLibraryTimingParEntry(), GetLibraryTimingTF1Entries(), LoadLibrary(), NOpChannels(), SetLibraryEntry(), SetLibraryReflT0Entry(), SetLibraryTimingParEntry(), SetLibraryTimingTF1Entry(), and StoreLibrary().

bool phot::PhotonVisibilityService::fUseCryoBoundary

private

Definition at line 93 of file PhotonVisibilityService.h.

Referenced by reconfigure().

sim::PhotonVoxelDef phot::PhotonVisibilityService::fVoxelDef

private

Definition at line 127 of file PhotonVisibilityService.h.

Referenced by GetAllVisibilities(), GetReflT0s(), GetTimingPar(), GetTimingTF1(), GetVisibility(), GetVoxelDef(), and reconfigure().

float phot::PhotonVisibilityService::fXmax

private

Definition at line 88 of file PhotonVisibilityService.h.

Referenced by reconfigure().

float phot::PhotonVisibilityService::fXmin

private

Definition at line 88 of file PhotonVisibilityService.h.

Referenced by reconfigure().

float phot::PhotonVisibilityService::fYmax

private

Definition at line 89 of file PhotonVisibilityService.h.

Referenced by reconfigure().

float phot::PhotonVisibilityService::fYmin

private

Definition at line 89 of file PhotonVisibilityService.h.

Referenced by reconfigure().

float phot::PhotonVisibilityService::fZmax

private

Definition at line 90 of file PhotonVisibilityService.h.

Referenced by reconfigure().

float phot::PhotonVisibilityService::fZmin

private

Definition at line 90 of file PhotonVisibilityService.h.

Referenced by reconfigure().

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

larsim/v06_53_00/source/larsim/PhotonPropagation/PhotonVisibilityService.h
larsim/v06_53_00/source/larsim/PhotonPropagation/PhotonVisibilityService_service.cc

Public Member Functions

Static Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation