HitSelector Class Reference

Inheritance diagram for HitSelector:

Public Member Functions
	HitSelector (fhicl::ParameterSet const &pset)
virtual	~HitSelector ()
	Destructor. More...
virtual void	reconfigure (fhicl::ParameterSet const &pset)
virtual void	produce (art::Event &e)
virtual void	beginJob ()
	Begin job method. More...
virtual void	endJob ()
	End job method. More...

Private Types
using	HitPtrVector = std::vector< art::Ptr< recob::Hit >>

Private Attributes
art::InputTag	fHitProducerLabel
	The full collection of hits. More...
std::vector< float >	fMinMaxPulseHeighMulti
	Max pulse height of a pulse train must be this large. More...
std::vector< float >	fMinPulseHeightMulti
	Multi hit snippets, minimum pulse height per plane. More...
std::vector< float >	fMinPulseWidthMulti
	Multi hit snippets, minimum pulse width per plane. More...
std::vector< float >	fMinPulseHeightSingle
	Single hit snippets, minimum pulse height per plane. More...
std::vector< float >	fMinPulseWidthSingle
	Single hit snippets, minimum pulse width per plane. More...
int	fNumEvent
	Number of events seen. More...

Detailed Description

Definition at line 38 of file HitSelector_module.cc.

Member Typedef Documentation

using HitSelector::HitPtrVector = std::vector<art::Ptr<recob::Hit>>

private

Definition at line 54 of file HitSelector_module.cc.

Constructor & Destructor Documentation

HitSelector::HitSelector ( fhicl::ParameterSet const &  pset )

explicit

Constructor.

Arguments:

pset - Fcl parameters.

Definition at line 77 of file HitSelector_module.cc.

                                                       : EDProducer{pset},
fNumEvent(0)
{
    reconfigure(pset);
    
    // let HitCollectionCreator declare that we are going to produce
    // hits and associations with wires and raw digits
    // (with no particular product label)
    recob::HitCollectionCreator::declare_products(producesCollector());
    
    // Report.
    mf::LogInfo("HitSelector") << "HitSelector configured\n";
}

HitSelector::~HitSelector ( )

virtual

Destructor.

Definition at line 93 of file HitSelector_module.cc.

{}

Member Function Documentation

void HitSelector::beginJob ( )

virtual

Begin job method.

Definition at line 115 of file HitSelector_module.cc.

{
//    auto const* detp = lar::providerFrom<detinfo::DetectorPropertiesService>();
//    auto const* geo  = lar::providerFrom<geo::Geometry>();
//    auto const* ts   = lar::providerFrom<detinfo::DetectorClocksService>();
}

void HitSelector::endJob ( )

virtual

End job method.

Definition at line 240 of file HitSelector_module.cc.

{
    mf::LogInfo("HitSelector")
    << "HitSelector statistics:\n"
    << "  Number of events = " << fNumEvent;
}

void HitSelector::produce ( art::Event &  evt )

virtual

Produce method.

Arguments:

evt - Art event.

This is the primary method. The goal is to produce a list of recob::Hit objects which are a "clean" subset of all hits and which are believed to be due to a neutrino interaction. It does this by considering input CosmicTag objects, relating them to PFParticles/Tracks and removing the hits associated to those objects which are believed to be Cosmic Rays.

Definition at line 135 of file HitSelector_module.cc.

{
    ++fNumEvent;
    
    // Start by looking up the original hits
    art::Handle< std::vector<recob::Hit> > hitHandle;
    evt.getByLabel(fHitProducerLabel, hitHandle);
    
    // also get the associated wires and raw digits;
    // we assume they have been created by the same module as the hits
    art::FindOneP<raw::RawDigit> hitToRawDigitAssns(hitHandle, evt, fHitProducerLabel);
    art::FindOneP<recob::Wire>   hitToWireAssns(    hitHandle, evt, fHitProducerLabel);
    
    // this object contains the hit collection
    // and its associations to wires and raw digits:
    recob::HitCollectionCreator hcol(evt,
                                     hitToWireAssns.isValid(),
                                     hitToRawDigitAssns.isValid()
                                     );

    // If there are no hits then there should be no output
    if (hitHandle.isValid())
    {
        HitPtrVector hitPtrVec;
        HitPtrVector hitSnippetVec;
        
        art::fill_ptr_vector(hitPtrVec, hitHandle);

        int lastHitIndex = 0;
        
        // Loop the hits and make some plots
        for(const auto& curHitPtr : hitPtrVec)
        {
            // Have we collected all of the hits on the same snippet?
            if (!hitSnippetVec.empty() && lastHitIndex >= curHitPtr->LocalIndex())
            {
                // Recover plane
                int plane = hitSnippetVec.front()->WireID().Plane;

                // Only worried about multi hit snippets
                if (hitSnippetVec.size() > 1)
                {
                    // Sort in order of largest to smallest pulse height
                    std::sort(hitSnippetVec.begin(),hitSnippetVec.end(),[](const auto& left, const auto& right){return left->PeakAmplitude() > right->PeakAmplitude();});
                    
                    float maxPulseHeight = hitSnippetVec.front()->PeakAmplitude();
                    
                    // Require the largest pulse in the train to be above threshold
                    if (maxPulseHeight > fMinMaxPulseHeighMulti[plane])
                    {
                        // Loop over all hits...
                        for(size_t idx = 0; idx < hitSnippetVec.size(); idx++)
                        {
                            art::Ptr<recob::Hit> hitPtr = hitSnippetVec[idx];
                            
                            float pulseHeight = hitPtr->PeakAmplitude();
                            float pulseWid    = hitPtr->RMS();
                            int   numDOF      = hitPtr->DegreesOfFreedom();
                            
                            // Watch for case of long pulse trains (numDOF == 1) and keep all of these no matter
                            // Also always keep the largest (first) pulse
                            // Otherwise, select on minimum pulse height and width
                            if (idx == 0 || numDOF == 1 || (pulseHeight > fMinPulseHeightMulti[plane] && pulseWid > fMinPulseWidthMulti[plane]))
                            {
                                art::Ptr<recob::Wire>   wire      = hitToWireAssns.at(hitPtr.key());
                                art::Ptr<raw::RawDigit> rawdigits = hitToRawDigitAssns.at(hitPtr.key());
                                
                                hcol.emplace_back(*hitPtr, wire, rawdigits);
                            }
                        }
                    }
                }
                // Here we handle the case of single hits 
                else
                {
                    art::Ptr<recob::Hit> hitPtr = hitSnippetVec.front();
                    
                    float  pulseHeight = hitPtr->PeakAmplitude();
                    float  pulseWid    = hitPtr->RMS();
                    
                    if (pulseHeight > fMinPulseHeightSingle[plane] && pulseWid > fMinPulseWidthSingle[plane])
                    {
                        art::Ptr<recob::Wire>   wire      = hitToWireAssns.at(hitPtr.key());
                        art::Ptr<raw::RawDigit> rawdigits = hitToRawDigitAssns.at(hitPtr.key());
                        
                        hcol.emplace_back(*hitPtr, wire, rawdigits);
                   }
                }
                
                hitSnippetVec.clear();
            }
 
            lastHitIndex = curHitPtr->LocalIndex();
            hitSnippetVec.push_back(curHitPtr);
        }
    }
    
    // put the hit collection and associations into the event
    hcol.put_into(evt);
    
    return;
}

void HitSelector::reconfigure ( fhicl::ParameterSet const &  pset )

virtual

Reconfigure method.

Arguments:

pset - Fcl parameter set.

Definition at line 103 of file HitSelector_module.cc.

{
    fHitProducerLabel      = pset.get<art::InputTag>("HitProducerLabel");
    fMinMaxPulseHeighMulti = pset.get<std::vector<float>>("MinMaxPulseHeightMulti");
    fMinPulseHeightMulti   = pset.get<std::vector<float>>("MinPulseHeightMulti");
    fMinPulseWidthMulti    = pset.get<std::vector<float>>("MinPulseWidthMulti");
    fMinPulseHeightSingle  = pset.get<std::vector<float>>("MinPulseHeightSingle");
    fMinPulseWidthSingle   = pset.get<std::vector<float>>("MinPulseWidthSingle");
}

Member Data Documentation

art::InputTag HitSelector::fHitProducerLabel

private

The full collection of hits.

Definition at line 57 of file HitSelector_module.cc.

std::vector<float> HitSelector::fMinMaxPulseHeighMulti

private

Max pulse height of a pulse train must be this large.

Definition at line 58 of file HitSelector_module.cc.

std::vector<float> HitSelector::fMinPulseHeightMulti

private

Multi hit snippets, minimum pulse height per plane.

Definition at line 59 of file HitSelector_module.cc.

std::vector<float> HitSelector::fMinPulseHeightSingle

private

Single hit snippets, minimum pulse height per plane.

Definition at line 61 of file HitSelector_module.cc.

std::vector<float> HitSelector::fMinPulseWidthMulti

private

Multi hit snippets, minimum pulse width per plane.

Definition at line 60 of file HitSelector_module.cc.

std::vector<float> HitSelector::fMinPulseWidthSingle

private

Single hit snippets, minimum pulse width per plane.

Definition at line 62 of file HitSelector_module.cc.

int HitSelector::fNumEvent

private

Number of events seen.

Definition at line 65 of file HitSelector_module.cc.

---

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

• HitSelector_module.cc