opdet::OpHitFinder Class Reference

Inheritance diagram for opdet::OpHitFinder:

Public Member Functions
	OpHitFinder (const fhicl::ParameterSet &)
virtual	~OpHitFinder ()
void	produce (art::Event &)

Private Member Functions
std::map< int, int >	GetChannelMap ()
std::vector< double >	GetSPEScales ()
std::vector< double >	GetSPEShifts ()

Private Attributes
std::string	fInputModule
std::string	fGenModule
std::vector< std::string >	fInputLabels
std::set< unsigned int >	fChannelMasks
pmtana::PulseRecoManager	fPulseRecoMgr
pmtana::PMTPulseRecoBase *	fThreshAlg
pmtana::PMTPedestalBase *	fPedAlg
Float_t	fHitThreshold
unsigned int	fMaxOpChannel
bool	fUseStartTime
calib::IPhotonCalibrator const *	fCalib = nullptr

Detailed Description

Definition at line 63 of file OpHitFinder_module.cc.

Constructor & Destructor Documentation

opdet::OpHitFinder::OpHitFinder ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 104 of file OpHitFinder_module.cc.

                                                        : EDProducer{pset}, fPulseRecoMgr()
  {
    // Indicate that the Input Module comes from .fcl
    fInputModule = pset.get<std::string>("InputModule");
    fGenModule = pset.get<std::string>("GenModule");
    fInputLabels = pset.get<std::vector<std::string>>("InputLabels");
    fUseStartTime = pset.get<bool>("UseStartTime", false);

    for (auto const& ch :
         pset.get<std::vector<unsigned int>>("ChannelMasks", std::vector<unsigned int>()))
      fChannelMasks.insert(ch);

    fHitThreshold = pset.get<float>("HitThreshold");
    bool useCalibrator = pset.get<bool>("UseCalibrator", false);

    auto const& geometry(*lar::providerFrom<geo::Geometry>());
    fMaxOpChannel = geometry.MaxOpChannel();

    if (useCalibrator) {
      // If useCalibrator, get it from ART
      fCalib = lar::providerFrom<calib::IPhotonCalibratorService>();
    }
    else {
      // If not useCalibrator, make an internal one based
      // on fhicl settings to hit finder.
      bool areaToPE = pset.get<bool>("AreaToPE");
      float SPEArea = pset.get<float>("SPEArea");
      float SPEShift = pset.get<float>("SPEShift", 0.);

      // Reproduce behavior from GetSPEScales()
      if (!areaToPE) SPEArea = 20;

      // Delete and replace if we are reconfiguring
      if (fCalib) { delete fCalib; }

      fCalib = new calib::PhotonCalibratorStandard(SPEArea, SPEShift, areaToPE);
    }

    // Initialize the rise time calculator tool
    auto const rise_alg_pset = pset.get_if_present<fhicl::ParameterSet>("RiseTimeCalculator");

    // Initialize the hit finder algorithm
    auto const hit_alg_pset = pset.get<fhicl::ParameterSet>("HitAlgoPset");
    std::string threshAlgName = hit_alg_pset.get<std::string>("Name");
    if (threshAlgName == "Threshold")
      fThreshAlg = thresholdAlgorithm<pmtana::AlgoThreshold>(hit_alg_pset, rise_alg_pset);
    else if (threshAlgName == "SiPM")
      fThreshAlg = thresholdAlgorithm<pmtana::AlgoSiPM>(hit_alg_pset, rise_alg_pset);
    else if (threshAlgName == "SlidingWindow")
      fThreshAlg = thresholdAlgorithm<pmtana::AlgoSlidingWindow>(hit_alg_pset, rise_alg_pset);
    else if (threshAlgName == "FixedWindow")
      fThreshAlg = thresholdAlgorithm<pmtana::AlgoFixedWindow>(hit_alg_pset, rise_alg_pset);
    else if (threshAlgName == "CFD")
      fThreshAlg = thresholdAlgorithm<pmtana::AlgoCFD>(hit_alg_pset, rise_alg_pset);
    else
      throw art::Exception(art::errors::UnimplementedFeature)
        << "Cannot find implementation for " << threshAlgName << " algorithm.\n";

    // Initialize the pedestal estimation algorithm
    auto const ped_alg_pset = pset.get<fhicl::ParameterSet>("PedAlgoPset");
    std::string pedAlgName = ped_alg_pset.get<std::string>("Name");
    if (pedAlgName == "Edges")
      fPedAlg = new pmtana::PedAlgoEdges(ped_alg_pset);
    else if (pedAlgName == "RollingMean")
      fPedAlg = new pmtana::PedAlgoRollingMean(ped_alg_pset);
    else if (pedAlgName == "UB")
      fPedAlg = new pmtana::PedAlgoUB(ped_alg_pset);
    else
      throw art::Exception(art::errors::UnimplementedFeature)
        << "Cannot find implementation for " << pedAlgName << " algorithm.\n";

    produces<std::vector<recob::OpHit>>();

    fPulseRecoMgr.AddRecoAlgo(fThreshAlg);
    fPulseRecoMgr.SetDefaultPedAlgo(fPedAlg);
  }

opdet::OpHitFinder::~OpHitFinder ( )

virtual

Definition at line 183 of file OpHitFinder_module.cc.

  {

    delete fThreshAlg;
    delete fPedAlg;
  }

Member Function Documentation

std::map<int, int> opdet::OpHitFinder::GetChannelMap ( )

private

std::vector<double> opdet::OpHitFinder::GetSPEScales ( )

private

std::vector<double> opdet::OpHitFinder::GetSPEShifts ( )

private

void opdet::OpHitFinder::produce

(

art::Event &

evt

)

Definition at line 192 of file OpHitFinder_module.cc.

  {

    // These is the storage pointer we will put in the event
    std::unique_ptr<std::vector<recob::OpHit>> HitPtr(new std::vector<recob::OpHit>);

    std::vector<const sim::BeamGateInfo*> beamGateArray;
    try {
      evt.getView(fGenModule, beamGateArray);
    }
    catch (art::Exception const& err) {
      if (err.categoryCode() != art::errors::ProductNotFound) throw;
    }

    auto const& geometry(*lar::providerFrom<geo::Geometry>());
    auto const clock_data =
      art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
    auto const& calibrator(*fCalib);
    //
    // Get the pulses from the event
    //

    // Load pulses into WaveformVector
    if (fChannelMasks.empty() && fInputLabels.size() < 2) {
      art::Handle<std::vector<raw::OpDetWaveform>> wfHandle;
      if (fInputLabels.empty())
        evt.getByLabel(fInputModule, wfHandle);
      else
        evt.getByLabel(fInputModule, fInputLabels.front(), wfHandle);
      assert(wfHandle.isValid());
      RunHitFinder(*wfHandle,
                   *HitPtr,
                   fPulseRecoMgr,
                   *fThreshAlg,
                   geometry,
                   fHitThreshold,
                   clock_data,
                   calibrator,
                   fUseStartTime);
    }
    else {

      // Reserve a large enough array
      int totalsize = 0;
      for (auto label : fInputLabels) {
        art::Handle<std::vector<raw::OpDetWaveform>> wfHandle;
        evt.getByLabel(fInputModule, label, wfHandle);
        if (!wfHandle.isValid()) continue; // Skip non-existent collections
        totalsize += wfHandle->size();
      }

      std::vector<raw::OpDetWaveform> WaveformVector;
      WaveformVector.reserve(totalsize);

      for (auto label : fInputLabels) {
        art::Handle<std::vector<raw::OpDetWaveform>> wfHandle;
        evt.getByLabel(fInputModule, label, wfHandle);
        if (!wfHandle.isValid()) continue; // Skip non-existent collections

        //WaveformVector.insert(WaveformVector.end(),
        //                      wfHandle->begin(), wfHandle->end());
        for (auto const& wf : *wfHandle) {
          if (fChannelMasks.find(wf.ChannelNumber()) != fChannelMasks.end()) continue;
          WaveformVector.push_back(wf);
        }
      }

      RunHitFinder(WaveformVector,
                   *HitPtr,
                   fPulseRecoMgr,
                   *fThreshAlg,
                   geometry,
                   fHitThreshold,
                   clock_data,
                   calibrator,
                   fUseStartTime);
    }
    // Store results into the event
    evt.put(std::move(HitPtr));
  }

Member Data Documentation

calib::IPhotonCalibrator const* opdet::OpHitFinder::fCalib = nullptr

private

Definition at line 91 of file OpHitFinder_module.cc.

std::set<unsigned int> opdet::OpHitFinder::fChannelMasks

private

Definition at line 81 of file OpHitFinder_module.cc.

std::string opdet::OpHitFinder::fGenModule

private

Definition at line 79 of file OpHitFinder_module.cc.

Float_t opdet::OpHitFinder::fHitThreshold

private

Definition at line 87 of file OpHitFinder_module.cc.

std::vector<std::string> opdet::OpHitFinder::fInputLabels

private

Definition at line 80 of file OpHitFinder_module.cc.

std::string opdet::OpHitFinder::fInputModule

private

Definition at line 78 of file OpHitFinder_module.cc.

unsigned int opdet::OpHitFinder::fMaxOpChannel

private

Definition at line 88 of file OpHitFinder_module.cc.

pmtana::PMTPedestalBase* opdet::OpHitFinder::fPedAlg

private

Definition at line 85 of file OpHitFinder_module.cc.

pmtana::PulseRecoManager opdet::OpHitFinder::fPulseRecoMgr

private

Definition at line 83 of file OpHitFinder_module.cc.

pmtana::PMTPulseRecoBase* opdet::OpHitFinder::fThreshAlg

private

Definition at line 84 of file OpHitFinder_module.cc.

bool opdet::OpHitFinder::fUseStartTime

private

Definition at line 89 of file OpHitFinder_module.cc.

---

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

• OpHitFinder_module.cc