t0::IndirectHitParticleAssns Class Reference

Inheritance diagram for t0::IndirectHitParticleAssns:

Public Member Functions
	IndirectHitParticleAssns (fhicl::ParameterSet const &pset)
	Constructor. More...
void	reconfigure (fhicl::ParameterSet const &pset) override
void	CreateHitParticleAssociations (art::Event &, HitParticleAssociations *) override
	This rebuilds the internal maps. More...
Public Member Functions inherited from t0::IHitParticleAssociations
virtual	~IHitParticleAssociations () noexcept=default
	Virtual Destructor. More...

Private Attributes
std::vector< art::InputTag >	fHitModuleLabelVec
art::InputTag	fMCParticleModuleLabel
art::InputTag	fHitPartAssnsModuleLabel

Detailed Description

Definition at line 34 of file IndirectHitParticleAssns_tool.cc.

Constructor & Destructor Documentation

t0::IndirectHitParticleAssns::IndirectHitParticleAssns ( fhicl::ParameterSet const &  pset )

explicit

Constructor.

Parameters

pset	Constructor.

Arguments:

pset - Fcl parameters.

Definition at line 67 of file IndirectHitParticleAssns_tool.cc.

{
    reconfigure(pset);

    // Report.
    mf::LogInfo("IndirectHitParticleAssns") << "Configured\n";
}

Member Function Documentation

void t0::IndirectHitParticleAssns::CreateHitParticleAssociations ( art::Event &  evt, HitParticleAssociations *  hitPartAssns  )

overridevirtual

This rebuilds the internal maps.

Rebuild method -> rebuild the basic maps to get truth information

Arguments:

event - the art event used to extract all information

Implements t0::IHitParticleAssociations.

Definition at line 96 of file IndirectHitParticleAssns_tool.cc.

{
    // This function handles the "indirect" creation of hit<-->MCParticle associations by trying to
    // use the previously created Hit<-->MCParticle associations
    // Its pretty much a brute force approach here... but time is short!
    //
    // First step is to recover the preexisting associations

    // Get a handle for the associations...
    auto partHitAssnsHandle = evt.getValidHandle< HitParticleAssociations >(fHitPartAssnsModuleLabel);
    auto mcParticleHandle   = evt.getValidHandle< std::vector<simb::MCParticle> >(fMCParticleModuleLabel);

    if (!partHitAssnsHandle.isValid() || !mcParticleHandle.isValid())
    {
        throw cet::exception("IndirectHitParticleAssns") << "===>> NO MCParticle <--> Hit associations found for run/subrun/event: " << evt.run() << "/" << evt.subRun() << "/" << evt.id().event();
    }

    // Loop over input hit collections
    for(const auto& inputTag : fHitModuleLabelVec)
    {
        // Look up the hits we want to process as well, since if they are not there then no point in proceeding
        art::Handle< std::vector<recob::Hit> > hitListHandle;
        evt.getByLabel(inputTag,hitListHandle);

        if(!hitListHandle.isValid()){
            mf::LogInfo("IndirectHitParticleAssns") << "===>> NO Hit collection found to process for run/subrun/event: " << evt.run() << "/" << evt.subRun() << "/" << evt.id().event() << "\n";
            continue;
        }

        // Go through the associations and build out our (hopefully sparse) data structure
        using ParticleDataPair         = std::pair<size_t, const anab::BackTrackerHitMatchingData*>;
        using MCParticleDataSet        = std::set<ParticleDataPair>;
        using TickToPartDataMap        = std::unordered_map<raw::TDCtick_t,   MCParticleDataSet>;
        using ChannelToTickPartDataMap = std::unordered_map<raw::ChannelID_t, TickToPartDataMap>;
        
        ChannelToTickPartDataMap chanToTickPartDataMap;
        
        // Build out the maps between hits/particles
        for(HitParticleAssociations::const_iterator partHitItr = partHitAssnsHandle->begin(); partHitItr != partHitAssnsHandle->end(); ++partHitItr)
        {
            const art::Ptr<simb::MCParticle>&       mcParticle = partHitItr->first;
            const art::Ptr<recob::Hit>&             recoHit    = partHitItr->second;
            const anab::BackTrackerHitMatchingData* data       = &partHitAssnsHandle->data(partHitItr);
            
            TickToPartDataMap& tickToPartDataMap = chanToTickPartDataMap[recoHit->Channel()];
            
            for(raw::TDCtick_t tick = recoHit->PeakTimeMinusRMS(); tick <= recoHit->PeakTimePlusRMS(); tick++)
            {
                tickToPartDataMap[tick].insert(ParticleDataPair(mcParticle.key(),data));
            }
        }
        
        // Armed with the map, process the hit list
        for(size_t hitIdx = 0; hitIdx < hitListHandle->size(); hitIdx++)
        {
            art::Ptr<recob::Hit> hit(hitListHandle,hitIdx);
            
            TickToPartDataMap& tickToPartDataMap = chanToTickPartDataMap[hit->Channel()];
            
            if (tickToPartDataMap.empty())
            {
                mf::LogInfo("IndirectHitParticleAssns") << "No channel information found for hit " << hit << "\n";
                continue;
            }
            
            // Keep track of results
            MCParticleDataSet particleDataSet;
            
            // Go through the ticks in this hit and recover associations
            for(raw::TDCtick_t tick = hit->PeakTimeMinusRMS(); tick <= hit->PeakTimePlusRMS(); tick++)
            {
                TickToPartDataMap::iterator hitInfoItr = tickToPartDataMap.find(tick);
                
                if (hitInfoItr != tickToPartDataMap.end())
                {
                    for(const auto& partData : hitInfoItr->second) particleDataSet.insert(partData);
                }
            }
            
            // Now create new associations for the hit in question
            for(const auto& partData : particleDataSet)
                hitPartAssns->addSingle(art::Ptr<simb::MCParticle>(mcParticleHandle, partData.first), hit, *partData.second);
        }
    }

    return;
}

void t0::IndirectHitParticleAssns::reconfigure ( fhicl::ParameterSet const &  pset )

overridevirtual

Reconfigure method.

Arguments:

pset - Fcl parameter set.

Implements t0::IHitParticleAssociations.

Definition at line 82 of file IndirectHitParticleAssns_tool.cc.

{
    fHitPartAssnsModuleLabel = pset.get<art::InputTag             >("HitPartAssnsLabel");
    fMCParticleModuleLabel   = pset.get<art::InputTag             >("MCParticleModuleLabel");
    fHitModuleLabelVec       = pset.get<std::vector<art::InputTag>>("HitModuleLabelVec");
}

Member Data Documentation

std::vector<art::InputTag> t0::IndirectHitParticleAssns::fHitModuleLabelVec

private

Definition at line 54 of file IndirectHitParticleAssns_tool.cc.

art::InputTag t0::IndirectHitParticleAssns::fHitPartAssnsModuleLabel

private

Definition at line 56 of file IndirectHitParticleAssns_tool.cc.

art::InputTag t0::IndirectHitParticleAssns::fMCParticleModuleLabel

private

Definition at line 55 of file IndirectHitParticleAssns_tool.cc.

---

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

• IndirectHitParticleAssns_tool.cc