BackTrackerService.h

Go to the documentation of this file.

// \file BackTrackerService.h
// \A service for backtracking reconstruction information to its truth information
//
// \author jason.stock@mines.sdsmt.edu
// Based on the original BackTracker by Brian Rebel (brebel@fnal.gov
//
#ifndef CHEAT_BACKTRACKERSERVICESERVICE_H
#define CHEAT_BACKTRACKERSERVICESERVICE_H

#include <vector>

#include "BackTracker.h"
#include "fhiclcpp/ParameterSet.h"
#include "art/Framework/Services/Registry/ActivityRegistry.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/View.h"

//Included Services
#include "larcore/Geometry/Geometry.h"
#include "larcorealg/Geometry/GeometryCore.h"
#include "lardata/DetectorInfoServices/DetectorClocksService.h"
#include "larsim/MCCheater/ParticleInventory.h"



namespace cheat{
  class BackTrackerService: private BackTracker
  {
    public:

      struct fhiclConfig{
        fhicl::Table < BackTracker::fhiclConfig > BackTrackerTable{
          fhicl::Name("BackTracker"),
            fhicl::Comment("This is the fhicl configuration of the BackTracker service provider.") };
      };


      using provider_type = BackTracker;
      const provider_type* provider() const 
      {return static_cast<provider_type const*>(this); }

      BackTrackerService( const fhicl::ParameterSet& pSet, art::ActivityRegistry& reg);
      BackTrackerService(const fhiclConfig& config, art::ActivityRegistry& reg);

      //Temporrary rebuild function (until I can make this lazy again.

      void Rebuild( const art::Event& evt );

      const std::vector < art::Ptr < sim::SimChannel > >& SimChannels() ;

      const std::vector < const sim::IDE* > TrackIdToSimIDEs_Ps(int const& id) ;
      const std::vector < const sim::IDE* > TrackIdToSimIDEs_Ps(int const& id, const geo::View_t view) ;

      const art::Ptr < sim::SimChannel > FindSimChannel( raw::ChannelID_t channel );

      const std::vector < sim::TrackIDE > ChannelToTrackIDEs(raw::ChannelID_t channel, const double hit_start_time, const double hit_end_time) ;

      const std::vector < sim::TrackIDE > HitToTrackIDEs(recob::Hit const& hit) ;
      const std::vector < sim::TrackIDE > HitToTrackIDEs(art::Ptr < recob::Hit > const& hit) ;

      const std::vector < int > HitToTrackIds(recob::Hit const& hit) ;

      const std::vector < sim::TrackIDE > HitToEveTrackIDEs(recob::Hit const& hit) ;
      const std::vector < sim::TrackIDE > HitToEveTrackIDEs(art::Ptr < recob::Hit > const& hit) ;

      const std::vector < art::Ptr < recob::Hit > > TrackIdToHits_Ps( const int& tkId, std::vector < art::Ptr < recob::Hit > > const& hitsIn ) ;
      const std::vector < art::Ptr < recob::Hit > > TrackIdToHits_Ps( const int& tkId ) ;

      const std::vector < std::vector < art::Ptr < recob::Hit > > > TrackIdsToHits_Ps( std::vector < int > const& tkIds, std:: vector < art::Ptr < recob::Hit > > const& hitsIn ) ;
      const std::vector < std::vector < art::Ptr < recob::Hit > > > TrackIdsToHits_Ps( std::vector < int > const& tkIds ) ;

      const std::vector< sim::IDE > HitToAvgSimIDEs ( recob::Hit const& hit) ;
      const std::vector< sim::IDE > HitToAvgSimIDEs ( art::Ptr<recob::Hit> hit) ;

      const std::vector< const sim::IDE* > HitToSimIDEs_Ps (recob::Hit const& hit) ;
      const std::vector< const sim::IDE* > HitToSimIDEs_Ps (art::Ptr< recob::Hit > const& hit) ;

      const std::vector<double> SimIDEsToXYZ( std::vector< sim::IDE > const& ides) ;
      const std::vector<double> SimIDEsToXYZ( std::vector< const sim::IDE* > const& ide_Ps) ;

      const std::vector<double> HitToXYZ(const recob::Hit& hit) ;
      const std::vector<double> HitToXYZ(art::Ptr<recob::Hit> const& hit) ;

      const double HitCollectionPurity( std::set<int> const& trackIds, std::vector< art::Ptr<recob::Hit> > const& hits) ;
      const double HitChargeCollectionPurity( std::set<int> const& trackIds, std::vector< art::Ptr<recob::Hit> > const&     hits) ;

      const double HitCollectionEfficiency( std::set<int> const& trackIds, std::vector< art::Ptr<recob::Hit> > const& hits, std::vector< art::Ptr<recob::Hit> > const& allhits, geo::View_t const& view) ;

      const double HitChargeCollectionEfficiency( std::set<int> trackIds, std::vector< art::Ptr<recob::Hit> > const&        hits,        std::vector< art::Ptr<recob::Hit> > const& allhits, geo::View_t const& view) ;

      const std::set<int> GetSetOfTrackIds() ;
      const std::set<int> GetSetOfEveIds() ;

      const std::set<int> GetSetOfTrackIds( std::vector< art::Ptr< recob::Hit > > const& hits ) ;
      const std::set<int> GetSetOfEveIds( std::vector< art::Ptr< recob::Hit > > const& hits ) ;

      const std::vector< double> SpacePointHitsToWeightedXYZ(std::vector<art::Ptr<recob::Hit>> const& hits);
      const std::vector< art::Ptr< recob::Hit > > SpacePointToHits_Ps(art::Ptr<recob::SpacePoint> const&      spt);
      const std::vector< double > SpacePointToXYZ( art::Ptr< recob::SpacePoint > const& spt);





    private:
      //FHICL parameters.
      //The BackTracker service has no parameters.

      //Configure services
//      cheat::BackTracker fBackTracker;

      const art::Event* fEvt=nullptr;

      //Prep functions go here.
      void priv_PrepEvent ( const art::Event& evt );
      void priv_PrepSimChannels (const art::Event& evt);
//      void priv_PrepAllHitList ();
      void priv_PrepFailed ();

      bool priv_CanRun (const art::Event& evt);

      bool priv_SimChannelsReady() { return BackTracker::SimChannelsReady();}
//      bool priv_AllHitListReady() { return BackTracker::AllHitListReady();}




  }; //class BackTrackerService

} //end namespace cheat
DECLARE_ART_SERVICE(cheat::BackTrackerService, LEGACY)



#endif //CHEAT_BACKTRACKERSERVICESERVICE_H