#include <ApaCrossCosmicIdAlg.h>

Classes
struct	BeamTime

struct	Config

Public Member Functions
	ApaCrossCosmicIdAlg (const core::ProviderManager &manager, const Config &config)

	ApaCrossCosmicIdAlg (const core::ProviderManager &manager, const fhicl::ParameterSet &pset)

	~ApaCrossCosmicIdAlg ()

void	reconfigure (const core::ProviderManager &manager, const Config &config)

std::pair< double, double >	MinApaDistance (const recob::Track &track, std::vector< double > &t0List, geo::TPCID &tpcid)

double	T0FromApaCross (const recob::Track &track, std::vector< art::Ptr< recob::Hit >> hits, std::map< geo::CryostatID, std::vector< double >> &t_zeros)

double	ApaDistance (recob::Track track, double t0, std::vector< art::Ptr< recob::Hit >> hits)

bool	ApaCrossCosmicId (const recob::Track &track, std::vector< art::Ptr< recob::Hit >> &hits, std::map< geo::CryostatID, std::vector< double >> &t_zeros)

Private Attributes
double	fDistanceLimit

double	fMaxApaDistance

double	fBeamTimeMin

double	fBeamTimeMax

geo::GeometryCore const *	fGeometry

detinfo::DetectorPropertiesStandard const *	fDetectorProperties

Detailed Description

Definition at line 35 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.h.

Constructor & Destructor Documentation

ana::ApaCrossCosmicIdAlg::ApaCrossCosmicIdAlg	(	const core::ProviderManager &	manager,
		const Config &	config
	)

Definition at line 5 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

                                                                                                 {
   this->reconfigure(manager, config);
 }

ana::ApaCrossCosmicIdAlg::ApaCrossCosmicIdAlg	(	const core::ProviderManager &	manager,
		const fhicl::ParameterSet &	pset
	)

inline

Definition at line 77 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.h.

77 :

78 ApaCrossCosmicIdAlg(manager, fhicl::Table<Config>(pset, {})()) {}

ana::ApaCrossCosmicIdAlg::ApaCrossCosmicIdAlg

ApaCrossCosmicIdAlg(const core::ProviderManager &manager, const Config &config)

Definition: sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc:5

ana::ApaCrossCosmicIdAlg::~ApaCrossCosmicIdAlg ( )

Definition at line 14 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

14 {}

Member Function Documentation

bool ana::ApaCrossCosmicIdAlg::ApaCrossCosmicId	(	const recob::Track &	track,
		std::vector< art::Ptr< recob::Hit >> &	hits,
		std::map< geo::CryostatID, std::vector< double >> &	t_zeros
	)

Definition at line 112 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

                                                                                                                                                         {
   double crossTimeThrough = T0FromApaCross(track, hits, t_zeros);
   // If the matched time is outside of the beam time then tag as a cosmic
   if(crossTimeThrough != -99999 && (crossTimeThrough < fBeamTimeMin || crossTimeThrough > fBeamTimeMax)) return true;
   return false;
 }

double ana::ApaCrossCosmicIdAlg::ApaDistance	(	recob::Track	track,
		double	t0,
		std::vector< art::Ptr< recob::Hit >>	hits
	)

Definition at line 100 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

                                                                                                     {
 
   std::vector<double> t0List {t0};
   // Determine the TPC from hit collection
   geo::TPCID tpc = CosmicIdUtils::DetectedInTPC(hits);
   // Get the distance to the APA at the given time
   std::pair<double, double> min = MinApaDistance(track, t0List, tpc);
   return min.first;
 
 }

std::pair< double, double > ana::ApaCrossCosmicIdAlg::MinApaDistance	(	const recob::Track &	track,
		std::vector< double > &	t0List,
		geo::TPCID &	tpcid
	)

Definition at line 29 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

                                                                                                                               {
   double crossTime = -99999;
 
   double minDist = 99999;
   double startX = track.Vertex().X();
   double endX = track.End().X();
   geo::Point_t point = track.Vertex();
 
 
   // if TPCID is invalid, return null values
   if (!tpcid) {
     return std::make_pair(minDist, crossTime);
   }
 
   // get the TPC drift direction
   const geo::TPCGeo &tpc = fGeometry->GetElement(tpcid);
   int drift_direction_x = tpc.DriftDirection() == geo::driftdir::kPos ? 1 : -1;
 
   // get the location of the APA in the TPC -- should be the x-position downstream of the drift direction
   double apa_X = (drift_direction_x == 1) ? tpc.MaxX() : tpc.MinX();
 
   // Don't try to shift tracks near the APA (may give artificially small distances)
   if (std::abs(startX - apa_X) < fMaxApaDistance || std::abs(endX - apa_X) < fMaxApaDistance)
     return std::make_pair(minDist, crossTime);
 
   // if drift direction positive use start/end with highest X
   if (drift_direction_x == 1) {
     point = endX < startX ? track.End() : track.Vertex();
   }
   else {
     point = endX > startX ? track.End() : track.Vertex();
   }
 
   // Shift track by all t0's
   for(auto const& t0 : t0List){
     // If particle crosses the APA before t = 0 the crossing point won't be reconstructed
     if(t0 < 0) continue;
     double shift = t0 * fDetectorProperties->DriftVelocity() * drift_direction_x;
     double shiftedX = point.X() + shift;
 
     //Check track still in TPC
     if (shiftedX - apa_X > fDistanceLimit) continue;
     // Calculate distance between start/end and APA
     double dist = std::abs(shiftedX - apa_X);
     if(dist < minDist) {
       minDist = dist;
       crossTime = t0;
     }
   }
 
   return std::make_pair(minDist, crossTime);
 
 }

void ana::ApaCrossCosmicIdAlg::reconfigure	(	const core::ProviderManager &	manager,
		const Config &	config
	)

Definition at line 16 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

                                                                                              {
   fDetectorProperties = manager.GetDetectorPropertiesProvider();
   fGeometry = manager.GetGeometryProvider();
   fDistanceLimit = config.DistanceLimit(); 
   fMaxApaDistance = config.MaxApaDistance();
   fBeamTimeMin = config.BeamTimeLimits().BeamTimeMin();
   fBeamTimeMax = config.BeamTimeLimits().BeamTimeMax();
 
   return;
 }

double ana::ApaCrossCosmicIdAlg::T0FromApaCross	(	const recob::Track &	track,
		std::vector< art::Ptr< recob::Hit >>	hits,
		std::map< geo::CryostatID, std::vector< double >> &	t_zeros
	)

Definition at line 85 of file sbnana/sbnanalysis/ana/SBNOscReco/CosmicIDAlgs/ApaCrossCosmicIdAlg.cc.

                                                                                                                                                        {
   // Determine the TPC from hit collection
   geo::TPCID tpc = CosmicIdUtils::DetectedInTPC(hits);
   if (tpc) {
     // Get the minimum distance to the APA and corresponding time
     std::pair<double, double> min = MinApaDistance(track, t_zeros.at(tpc), tpc);
     // Check the distance is within allowed limit
     if(min.first < fDistanceLimit) return min.second;
   }
 
   return -99999;
 
 }