Inheritance diagram for opdet::OpHitFinder:

Public Member Functions
	OpHitFinder (const fhicl::ParameterSet &)

virtual	~OpHitFinder ()

void	produce (art::Event &)

Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDProducer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept

	Producer (fhicl::ParameterSet const &)

	Producer (Producer const &)=delete

	Producer (Producer &&)=delete

Producer &	operator= (Producer const &)=delete

Producer &	operator= (Producer &&)=delete

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

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

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 50 of file OpHitFinder_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 91 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 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 = new pmtana::AlgoThreshold(hit_alg_pset);
     else if (threshAlgName == "SiPM")
       fThreshAlg = new pmtana::AlgoSiPM(hit_alg_pset);
     else if (threshAlgName == "SlidingWindow")
       fThreshAlg = new pmtana::AlgoSlidingWindow(hit_alg_pset);
     else if (threshAlgName == "FixedWindow")
       fThreshAlg = new pmtana::AlgoFixedWindow(hit_alg_pset);
     else if (threshAlgName == "CFD")
       fThreshAlg = new pmtana::AlgoCFD(hit_alg_pset);
     else
       throw art::Exception(art::errors::UnimplementedFeature)
         << "Cannot find implementation for " << threshAlgName << " algorithm.\n";
 
     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 166 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 )

virtual

Implements art::EDProducer.

Definition at line 175 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 78 of file OpHitFinder_module.cc.

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

private

Definition at line 68 of file OpHitFinder_module.cc.

std::string opdet::OpHitFinder::fGenModule

private

Definition at line 66 of file OpHitFinder_module.cc.

Float_t opdet::OpHitFinder::fHitThreshold

private

Definition at line 74 of file OpHitFinder_module.cc.

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

private

Definition at line 67 of file OpHitFinder_module.cc.

std::string opdet::OpHitFinder::fInputModule

private

Definition at line 65 of file OpHitFinder_module.cc.

unsigned int opdet::OpHitFinder::fMaxOpChannel

private

Definition at line 75 of file OpHitFinder_module.cc.

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

private

Definition at line 72 of file OpHitFinder_module.cc.

pmtana::PulseRecoManager opdet::OpHitFinder::fPulseRecoMgr

private

Definition at line 70 of file OpHitFinder_module.cc.

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

private

Definition at line 71 of file OpHitFinder_module.cc.

bool opdet::OpHitFinder::fUseStartTime

private

Definition at line 76 of file OpHitFinder_module.cc.

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

larana/larana/OpticalDetector/OpHitFinder_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation