OpHitAlg.cxx

Go to the documentation of this file.

// -*- mode: c++; c-basic-offset: 2; -*-
#include "OpHitAlg.h"

#include "larana/OpticalDetector/OpHitFinder/PulseRecoManager.h"
#include "larcorealg/Geometry/GeometryCore.h"
#include "lardataalg/DetectorInfo/DetectorClocks.h"
#include "lardataalg/DetectorInfo/ElecClock.h"
#include "lardataobj/RawData/OpDetWaveform.h"
#include "lardataobj/RecoBase/OpHit.h"
#include "larreco/Calibrator/IPhotonCalibrator.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

#include <vector>

namespace opdet {
  //----------------------------------------------------------------------------
  void RunHitFinder(std::vector<raw::OpDetWaveform> const& opDetWaveformVector,
                    std::vector<recob::OpHit>& hitVector,
                    pmtana::PulseRecoManager const& pulseRecoMgr,
                    pmtana::PMTPulseRecoBase const& threshAlg,
                    geo::GeometryCore const& geometry,
                    float hitThreshold,
                    detinfo::DetectorClocksData const& clocksData,
                    calib::IPhotonCalibrator const& calibrator,
                    bool use_start_time)
  {

    for (auto const& waveform : opDetWaveformVector) {

      const int channel = static_cast<int>(waveform.ChannelNumber());

      if (!geometry.IsValidOpChannel(channel)) {
        mf::LogError("OpHitFinder")
          << "Error! unrecognized channel number " << channel << ". Ignoring pulse";
        continue;
      }

      pulseRecoMgr.Reconstruct(waveform);

      // Get the result
      auto const& pulses = threshAlg.GetPulses();

      const double timeStamp = waveform.TimeStamp();

      for (auto const& pulse : pulses)
        ConstructHit(hitThreshold,
                     channel,
                     timeStamp,
                     pulse,
                     hitVector,
                     clocksData,
                     calibrator,
                     use_start_time);
    }
  }

  //----------------------------------------------------------------------------
  void ConstructHit(float hitThreshold,
                    int channel,
                    double timeStamp,
                    pmtana::pulse_param const& pulse,
                    std::vector<recob::OpHit>& hitVector,
                    detinfo::DetectorClocksData const& clocksData,
                    calib::IPhotonCalibrator const& calibrator,
                    bool use_start_time)
  {

    if (pulse.peak < hitThreshold) return;

    double absTime = timeStamp + clocksData.OpticalClock().TickPeriod() *
                                   (use_start_time ? pulse.t_start : pulse.t_max);

    double relTime = absTime - clocksData.TriggerTime();

    double startTime =
      timeStamp + clocksData.OpticalClock().TickPeriod() * pulse.t_start - clocksData.TriggerTime();

    double riseTime = clocksData.OpticalClock().TickPeriod() * pulse.t_rise;

    int frame = clocksData.OpticalClock().Frame(timeStamp);

    double PE = 0.0;
    if (calibrator.UseArea())
      PE = calibrator.PE(pulse.area, channel);
    else
      PE = calibrator.PE(pulse.peak, channel);

    double width = (pulse.t_end - pulse.t_start) * clocksData.OpticalClock().TickPeriod();

    hitVector.emplace_back(channel,
                           relTime,
                           absTime,
                           startTime,
                           riseTime,
                           frame,
                           width,
                           pulse.area,
                           pulse.peak,
                           PE,
                           0.0);
  }

} // End namespace opdet