Inheritance diagram for t0::MCParticleTrackMatching:

Public Member Functions
	MCParticleTrackMatching (fhicl::ParameterSet const &p)

	MCParticleTrackMatching (MCParticleTrackMatching const &)=delete

	MCParticleTrackMatching (MCParticleTrackMatching &&)=delete

MCParticleTrackMatching &	operator= (MCParticleTrackMatching const &)=delete

MCParticleTrackMatching &	operator= (MCParticleTrackMatching &&)=delete

void	produce (art::Event &e) override

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 Attributes
art::InputTag	fTrackModuleLabel

art::InputTag	fTrackHitAssnLabel

art::InputTag	fHitModuleLabel

art::InputTag	fHitParticleAssnLabel

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 43 of file MCParticleTrackMatching_module.cc.

Constructor & Destructor Documentation

t0::MCParticleTrackMatching::MCParticleTrackMatching ( fhicl::ParameterSet const & p )

explicit

Definition at line 68 of file MCParticleTrackMatching_module.cc.

   : EDProducer{p}
 {
   fTrackModuleLabel = p.get<art::InputTag>("TrackModuleLabel");
   fTrackHitAssnLabel = p.get<art::InputTag>("TrackHitAssnLabel",fTrackModuleLabel);
   fHitModuleLabel = p.get<art::InputTag>("HitModuleLabel");
   fHitParticleAssnLabel = p.get<art::InputTag>("HitParticleAssnLabel");
 
   produces< art::Assns<recob::Track , simb::MCParticle, anab::BackTrackerMatchingData > > ();
 }

t0::MCParticleTrackMatching::MCParticleTrackMatching ( MCParticleTrackMatching const & )

delete

t0::MCParticleTrackMatching::MCParticleTrackMatching ( MCParticleTrackMatching && )

delete

Member Function Documentation

MCParticleTrackMatching& t0::MCParticleTrackMatching::operator= ( MCParticleTrackMatching const & )

delete

MCParticleTrackMatching& t0::MCParticleTrackMatching::operator= ( MCParticleTrackMatching && )

delete

void t0::MCParticleTrackMatching::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 79 of file MCParticleTrackMatching_module.cc.

 {
   if (evt.isRealData()) return;
 
   //auto mcpartHandle = evt.getValidHandle< std::vector<simb::MCParticle> >("largeant");
   std::unique_ptr< art::Assns<recob::Track, simb::MCParticle, anab::BackTrackerMatchingData > > MCPartTrackassn( new art::Assns<recob::Track, simb::MCParticle, anab::BackTrackerMatchingData >);
 
 
   double maxe = -1;
   double tote = 0;
   // int    trkid = -1;
   //int    maxtrkid = -1;
   //double maxn = -1;
   //double totn = 0;
   //int maxntrkid = -1;
 
   anab::BackTrackerMatchingData btdata;
   std::unordered_map<int,double> trkide;
 
   art::Handle< std::vector<recob::Track> > trackListHandle;
   evt.getByLabel(fTrackModuleLabel,trackListHandle);
 
   art::Handle< std::vector<recob::Hit> > hitListHandle;
   evt.getByLabel(fHitModuleLabel,hitListHandle);
 
   if(!trackListHandle.isValid()){
     std::cerr << "Track handle is not valid!" << std::endl;
     return;
   }
 
   if(!hitListHandle.isValid()){
     std::cerr << "Hit handle is not valid!" << std::endl;
     return;
   }
 
   auto const& trackList(*trackListHandle);
   art::FindManyP<recob::Hit> fmtht(trackListHandle, evt, fTrackHitAssnLabel);
   //auto const& mcpartList(*mcpartHandle);
 
   for(size_t i_t=0; i_t<trackList.size(); ++i_t){
     art::Ptr<recob::Track> trkPtr(trackListHandle,i_t);
     trkide.clear();
     tote = 0; maxe=-1; art::Ptr<simb::MCParticle> maxp;
 
     std::vector< art::Ptr<recob::Hit> > allHits = fmtht.at(i_t);
 
     std::vector<anab::BackTrackerHitMatchingData const*> bthmd_vec;
     std::vector< art::Ptr<simb::MCParticle> > matchedParticlePtrs;
 
     art::FindManyP<simb::MCParticle,anab::BackTrackerHitMatchingData>
       particles_per_hit(hitListHandle, evt, fHitParticleAssnLabel);
 
     for(size_t i_h=0; i_h<allHits.size(); ++i_h){
       bthmd_vec.clear(); matchedParticlePtrs.clear();
       particles_per_hit.get(allHits[i_h].key(),matchedParticlePtrs,bthmd_vec);
 
       for(size_t i_p=0; i_p<matchedParticlePtrs.size(); ++i_p){
         trkide[ matchedParticlePtrs[i_p]->TrackId() ] += bthmd_vec[i_p]->energy;
         tote += bthmd_vec[i_p]->energy;
         if( trkide[ matchedParticlePtrs[i_p]->TrackId() ] > maxe ){
           maxe = trkide[ matchedParticlePtrs[i_p]->TrackId() ];
           maxp = matchedParticlePtrs[i_p];
         }
       }//end loop over particles per hit
 
     }//end loop over hits
 
     btdata.cleanliness = maxe/tote;
     if(maxe>0)
       MCPartTrackassn->addSingle(trkPtr, maxp, btdata);
 
   }//end loop over tracks
 
   evt.put(std::move(MCPartTrackassn));
 } // Produce