shower::ShowerEnergyAlg Class Reference

#include "ShowerEnergyAlg.h"

Classes
struct	LinearFunction

Public Member Functions
	ShowerEnergyAlg (fhicl::ParameterSet const &pset)
double	ShowerEnergy (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit >> const &hits, geo::PlaneID::PlaneID_t plane) const

Private Attributes
std::array< LinearFunction, 3 > const	fLinearFunctions

Detailed Description

Definition at line 32 of file ShowerEnergyAlg.h.

Constructor & Destructor Documentation

shower::ShowerEnergyAlg::ShowerEnergyAlg ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 26 of file ShowerEnergyAlg.cxx.

References fhicl::ParameterSet::get().

  : fLinearFunctions{{{pset.get<double>("UGradient"), pset.get<double>("UIntercept")},
                      {pset.get<double>("VGradient"), pset.get<double>("VIntercept")},
                      {pset.get<double>("ZGradient"), pset.get<double>("ZIntercept")}}}
{}

Member Function Documentation

double shower::ShowerEnergyAlg::ShowerEnergy	(	detinfo::DetectorClocksData const &	clockData,
		detinfo::DetectorPropertiesData const &	detProp,
		std::vector< art::Ptr< recob::Hit >> const &	hits,
		geo::PlaneID::PlaneID_t	plane
	)		const

This overload is preferred as it does not rely on the cached DetectorProperties data member.

Definition at line 32 of file ShowerEnergyAlg.cxx.

References detinfo::DetectorPropertiesData::ElectronLifetime(), fLinearFunctions, hits(), detinfo::sampling_rate(), and lar::to_element.

Referenced by shower::EMShowerAlg::MakeShower().

{
  // Should we throw instead if the plane is not in the range [0,3)?
  if (plane >= fLinearFunctions.size()) { return 0.; }

  auto const coeff = sampling_rate(clockData) / (detprop.ElectronLifetime() * 1e3);

  auto in_plane = [plane](auto const& hit) { return hit.WireID().Plane == plane; };
  auto charge = [coeff](auto const& hit) {
    return hit.Integral() * std::exp(coeff * hit.PeakTime());
  };

  double const totalCharge =
    ranges::accumulate(hits | ranges::views::transform(to_element) |
                         ranges::views::filter(in_plane) | ranges::views::transform(charge),
                       0.);

  return fLinearFunctions[plane].energy_from(totalCharge);
}

Member Data Documentation

std::array<LinearFunction, 3> const shower::ShowerEnergyAlg::fLinearFunctions

private

Definition at line 49 of file ShowerEnergyAlg.h.

Referenced by ShowerEnergy().

---

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

• larreco/v09_25_00/source/larreco/RecoAlg/ShowerEnergyAlg.h
• larreco/v09_25_00/source/larreco/RecoAlg/ShowerEnergyAlg.cxx