FlashHypothesisCalculator.cxx

Go to the documentation of this file.

#include "FlashHypothesisCalculator.h"

#include "larana/OpticalDetector/FlashHypothesis.h"
#include "lardataalg/Utilities/MappedContainer.h"

std::vector<double> opdet::FlashHypothesisCalculator::SegmentMidpoint(TVector3 const& pt1, TVector3 const& pt2, float XOffset)
{
  std::vector<double> xyz_segment(3);
  xyz_segment[0] = 0.5*(pt2.x()+pt1.x()) + XOffset;
  xyz_segment[1] = 0.5*(pt2.y()+pt1.y());
  xyz_segment[2] = 0.5*(pt2.z()+pt1.z());
  return xyz_segment;
}

void opdet::FlashHypothesisCalculator::FillFlashHypothesis(const float& yield,
                                                           const float& dEdx,
                                                           const TVector3& pt1,
                                                           const TVector3& pt2,
                                                           const std::vector<float>& qe_vector,
                                                           phot::MappedCounts_t const& vis_vector,
                                                           FlashHypothesis& hyp)
{

  if(qe_vector.size()!=hyp.GetVectorSize() || !vis_vector)
    throw std::runtime_error("ERROR in FlashHypothesisCalculator: vector sizes not equal!");

  const float total_yield = yield*dEdx*(pt2-pt1).Mag();

  for(size_t i_chan=0; i_chan<hyp.GetVectorSize(); i_chan++)
    hyp.SetHypothesisAndError(i_chan,total_yield*vis_vector[i_chan]*qe_vector[i_chan]);

}