Inheritance diagram for t0::IndirectHitParticleAssns:

Public Member Functions
	IndirectHitParticleAssns (fhicl::ParameterSet const &pset)
	Constructor. More...

	~IndirectHitParticleAssns ()
	Destructor. More...

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

void	CreateHitParticleAssociations (art::Event &, HitParticleAssociations *) override
	This rebuilds the internal maps. More...

Private Attributes
art::InputTag	fHitModuleLabel

art::InputTag	fMCParticleModuleLabel

art::InputTag	fHitPartAssnsModuleLabel

Detailed Description

Definition at line 36 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 74 of file IndirectHitParticleAssns_tool.cc.

References reconfigure().

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

t0::IndirectHitParticleAssns::~IndirectHitParticleAssns ( )

Destructor.

Definition at line 84 of file IndirectHitParticleAssns_tool.cc.

85 {}

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 108 of file IndirectHitParticleAssns_tool.cc.

References art::Assns< L, R, D >::addSingle(), recob::Hit::Channel(), DEFINE_ART_CLASS_TOOL, art::EventID::event(), fHitModuleLabel, fHitPartAssnsModuleLabel, fMCParticleModuleLabel, art::DataViewImpl::getByLabel(), art::DataViewImpl::getValidHandle(), art::Event::id(), art::Handle< T >::isValid(), art::Ptr< T >::key(), recob::Hit::PeakTimeMinusRMS(), recob::Hit::PeakTimePlusRMS(), art::Event::run(), and art::Event::subRun().

 {
     // 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();
     }
     
     // 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(fHitModuleLabel,hitListHandle);
     
     if(!hitListHandle.isValid()){
         throw cet::exception("IndirectHitParticleAssns") << "===>> NO Hit collection found to process for run/subrun/event: " << evt.run() << "/" << evt.subRun() << "/" << evt.id().event();
     }
     
     // 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 94 of file IndirectHitParticleAssns_tool.cc.

References fHitModuleLabel, fHitPartAssnsModuleLabel, fMCParticleModuleLabel, and fhicl::ParameterSet::get().

Referenced by IndirectHitParticleAssns().

 {
     fHitPartAssnsModuleLabel = pset.get<art::InputTag>("HitPartAssnsLabel");
     fMCParticleModuleLabel   = pset.get<art::InputTag>("MCParticleModuleLabel");
     fHitModuleLabel          = pset.get<art::InputTag>("HitModuleLabel");
 }