opdet::FlashHypothesisCreator Class Reference

#include "FlashHypothesisCreator.h"

Public Types
using	Providers_t = lar::ProviderPack< geo::GeometryCore, detinfo::LArProperties >
	Set of service providers used in the common(est) interface. More...

Public Member Functions
	FlashHypothesisCreator ()
FlashHypothesisCollection	GetFlashHypothesisCollection (recob::Track const &track, std::vector< float > const &dEdxVector, Providers_t providers, phot::PhotonVisibilityService const &pvs, opdet::OpDigiProperties const &opdigip, float XOffset=0)
FlashHypothesisCollection	GetFlashHypothesisCollection (sim::MCTrack const &mctrack, std::vector< float > const &dEdxVector, Providers_t providers, phot::PhotonVisibilityService const &pvs, opdet::OpDigiProperties const &opdigip, float XOffset=0)
FlashHypothesisCollection	GetFlashHypothesisCollection (std::vector< TVector3 > const &trajVector, std::vector< float > const &dEdxVector, Providers_t providers, phot::PhotonVisibilityService const &pvs, opdet::OpDigiProperties const &opdigip, float XOffset=0)
FlashHypothesisCollection	GetFlashHypothesisCollection (TVector3 const &pt1, TVector3 const &pt2, float const &dEdx, Providers_t providers, phot::PhotonVisibilityService const &pvs, opdet::OpDigiProperties const &opdigip, float XOffset=0)

Private Member Functions
FlashHypothesisCollection	CreateFlashHypothesesFromSegment (TVector3 const &pt1, TVector3 const &pt2, float const &dEdx, Providers_t providers, phot::PhotonVisibilityService const &pvs, opdet::OpDigiProperties const &opdigip, float XOffset)

Private Attributes
FlashHypothesisCalculator	_calc

Detailed Description

Definition at line 34 of file FlashHypothesisCreator.h.

Member Typedef Documentation

using opdet::FlashHypothesisCreator::Providers_t = lar::ProviderPack<geo::GeometryCore, detinfo::LArProperties>

Set of service providers used in the common(est) interface.

Definition at line 39 of file FlashHypothesisCreator.h.

Constructor & Destructor Documentation

opdet::FlashHypothesisCreator::FlashHypothesisCreator ( )

inline

Definition at line 41 of file FlashHypothesisCreator.h.

References CreateFlashHypothesesFromSegment(), GetFlashHypothesisCollection(), and track.

{}

Member Function Documentation

opdet::FlashHypothesisCollection opdet::FlashHypothesisCreator::CreateFlashHypothesesFromSegment ( TVector3 const &  pt1, TVector3 const &  pt2, float const &  dEdx, Providers_t  providers, phot::PhotonVisibilityService const &  pvs, opdet::OpDigiProperties const &  opdigip, float  XOffset  )

private

Definition at line 123 of file FlashHypothesisCreator.cxx.

References _calc, opdet::FlashHypothesisCalculator::FillFlashHypothesis(), lar::ProviderPack< Providers >::get(), phot::PhotonVisibilityService::GetAllVisibilities(), opdet::OpDigiProperties::QE(), opdet::FlashHypothesisCalculator::SegmentMidpoint(), and opdet::FlashHypothesisCollection::SetPromptHypAndPromptFraction().

Referenced by FlashHypothesisCreator(), and GetFlashHypothesisCollection().

{
  auto const* geom = providers.get<geo::GeometryCore>();
  auto const* larp = providers.get<detinfo::LArProperties>();
  auto const nOpDets = geom->NOpDets();
  FlashHypothesisCollection fhc(nOpDets);

  FlashHypothesis prompt_hyp = FlashHypothesis(nOpDets);
  
  std::vector<double> xyz_segment(_calc.SegmentMidpoint(pt1,pt2,XOffset));
  
  //get the visibility vector
  const float* PointVisibility = pvs.GetAllVisibilities(&xyz_segment[0]);
  
  //check visibility pointer, as it may be null if given a y/z outside some range
  if (!PointVisibility) return fhc;
  
  //klugey ... right now, set a qe_vector that gives constant qe across all opdets
  std::vector<float> qe_vector(nOpDets,opdigip.QE());
  _calc.FillFlashHypothesis(larp->ScintYield()*larp->ScintYieldRatio(),
                            dEdx,
                            pt1,pt2,
                            qe_vector,
                            PointVisibility,
                            prompt_hyp);
  
  fhc.SetPromptHypAndPromptFraction(prompt_hyp,larp->ScintYieldRatio());
  return fhc;
}

opdet::FlashHypothesisCollection opdet::FlashHypothesisCreator::GetFlashHypothesisCollection	(	recob::Track const &	track,
		std::vector< float > const &	dEdxVector,
		Providers_t	providers,
		phot::PhotonVisibilityService const &	pvs,
		opdet::OpDigiProperties const &	opdigip,
		float	XOffset = 0
	)

Title: FlashHypothesis Creator Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov)

Description: Class that produces a flash hypothesis for a trajectory. Input: Trajectory (std::vector<TVector3> objects) Output: FlashHypotheses

Definition at line 13 of file FlashHypothesisCreator.cxx.

References CreateFlashHypothesesFromSegment(), lar::ProviderPack< Providers >::get(), recob::Track::LocationAtPoint(), recob::Track::NumberTrajectoryPoints(), and pt.

Referenced by FlashHypothesisCreator(), and opdet::FlashHypothesisAnaAlg::RunComparison().

{
  bool interpolate_dEdx=false;
  if(track.NumberTrajectoryPoints() == dEdxVector.size())
    interpolate_dEdx=true;
  else if(track.NumberTrajectoryPoints() == dEdxVector.size()+1)
    interpolate_dEdx=false;
  else
    throw "ERROR in FlashHypothesisCreator: dEdx vector size not compatible with track size.";

  auto const* geom = providers.get<geo::GeometryCore>();
  FlashHypothesisCollection fhc(geom->NOpDets());
  for(size_t pt=1; pt<track.NumberTrajectoryPoints(); pt++){
    if(interpolate_dEdx)
      fhc = fhc + CreateFlashHypothesesFromSegment(track.LocationAtPoint(pt-1),
                                                   track.LocationAtPoint(pt),
                                                   0.5*(dEdxVector[pt]+dEdxVector[pt-1]),
                                                   providers,pvs,opdigip,XOffset);
    else
      fhc = fhc + CreateFlashHypothesesFromSegment(track.LocationAtPoint(pt-1),
                                                   track.LocationAtPoint(pt),
                                                   dEdxVector[pt-1],
                                                   providers,pvs,opdigip,XOffset);
  }
  return fhc;
}

opdet::FlashHypothesisCollection opdet::FlashHypothesisCreator::GetFlashHypothesisCollection	(	sim::MCTrack const &	mctrack,
		std::vector< float > const &	dEdxVector,
		Providers_t	providers,
		phot::PhotonVisibilityService const &	pvs,
		opdet::OpDigiProperties const &	opdigip,
		float	XOffset = 0
	)

Definition at line 46 of file FlashHypothesisCreator.cxx.

References CreateFlashHypothesesFromSegment(), lar::ProviderPack< Providers >::get(), and pt.

{
  bool interpolate_dEdx=false;
  if(mctrack.size() == dEdxVector.size())
    interpolate_dEdx=true;
  else if(mctrack.size() == dEdxVector.size()+1)
    interpolate_dEdx=false;
  else
    throw "ERROR in FlashHypothesisCreator: dEdx vector size not compatible with mctrack size.";

  auto const* geom = providers.get<geo::GeometryCore>();
  FlashHypothesisCollection fhc(geom->NOpDets());
  for(size_t pt=1; pt<mctrack.size(); pt++){
    if(interpolate_dEdx)
      fhc = fhc + CreateFlashHypothesesFromSegment(mctrack[pt-1].Position().Vect(),
                                                   mctrack[pt].Position().Vect(),
                                                   0.5*(dEdxVector[pt]+dEdxVector[pt-1]),
                                                   providers,pvs,opdigip,XOffset);
    else
      fhc = fhc + CreateFlashHypothesesFromSegment(mctrack[pt-1].Position().Vect(),
                                                   mctrack[pt].Position().Vect(),
                                                   dEdxVector[pt-1],
                                                   providers,pvs,opdigip,XOffset);
  }
  return fhc;
}

opdet::FlashHypothesisCollection opdet::FlashHypothesisCreator::GetFlashHypothesisCollection	(	std::vector< TVector3 > const &	trajVector,
		std::vector< float > const &	dEdxVector,
		Providers_t	providers,
		phot::PhotonVisibilityService const &	pvs,
		opdet::OpDigiProperties const &	opdigip,
		float	XOffset = 0
	)

Definition at line 79 of file FlashHypothesisCreator.cxx.

References CreateFlashHypothesesFromSegment(), lar::ProviderPack< Providers >::get(), and pt.

{
  bool interpolate_dEdx=false;
  if(trajVector.size() == dEdxVector.size())
    interpolate_dEdx=true;
  else if(trajVector.size() == dEdxVector.size()+1)
    interpolate_dEdx=false;
  else
    throw "ERROR in FlashHypothesisCreator: dEdx vector size not compatible with trajVector size.";

  auto const* geom = providers.get<geo::GeometryCore>();
  FlashHypothesisCollection fhc(geom->NOpDets());
  for(size_t pt=1; pt<trajVector.size(); pt++){
    if(interpolate_dEdx)
      fhc = fhc + CreateFlashHypothesesFromSegment(trajVector[pt-1],
                                                   trajVector[pt],
                                                   0.5*(dEdxVector[pt]+dEdxVector[pt-1]),
                                                   providers,pvs,opdigip,XOffset);
    else
      fhc = fhc + CreateFlashHypothesesFromSegment(trajVector[pt-1],
                                                   trajVector[pt],
                                                   dEdxVector[pt-1],
                                                   providers,pvs,opdigip,XOffset);
  }
  return fhc;
}

opdet::FlashHypothesisCollection opdet::FlashHypothesisCreator::GetFlashHypothesisCollection	(	TVector3 const &	pt1,
		TVector3 const &	pt2,
		float const &	dEdx,
		Providers_t	providers,
		phot::PhotonVisibilityService const &	pvs,
		opdet::OpDigiProperties const &	opdigip,
		float	XOffset = 0
	)

Definition at line 112 of file FlashHypothesisCreator.cxx.

References CreateFlashHypothesesFromSegment().

{
  return CreateFlashHypothesesFromSegment(pt1,pt2,dEdx,providers,pvs,opdigip,XOffset);
}

Member Data Documentation

FlashHypothesisCalculator opdet::FlashHypothesisCreator::_calc

private

Definition at line 79 of file FlashHypothesisCreator.h.

Referenced by CreateFlashHypothesesFromSegment().

---

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

• larana/v06_21_01/source/larana/OpticalDetector/FlashHypothesisCreator.h
• larana/v06_21_01/source/larana/OpticalDetector/FlashHypothesisCreator.cxx