rb::Track Class Reference

A rb::Prong with full reconstructed trajectory. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-07/RecoBase/Track.h"

Inheritance diagram for rb::Track:

Public Types
enum	ESaveGeVMode { kInitializeTotalGeV = 0, kResetTotalGeV = 1 }

Public Member Functions
virtual	~Track ()
void	AppendTrajectoryPoint (TVector3 pt)
void	AppendTrajectoryPoint (double pv, double pz)
void	PrependTrajectoryPoint (TVector3 pt)
	Support constructing tracks backwards. More...
void	PrependTrajectoryPoint (double pv, double pz)
size_t	NTrajectoryPoints () const
std::vector< TVector3 > const &	Trajectory () const
	return a constant reference to the track trajectory points More...
TVector3	TrajectoryPoint (unsigned int i) const
	The ith point on the trajectory, a 3-vector in cm. More...
virtual void	SetStart (TVector3 start)
virtual void	SetStart (double v0, double z0)
TVector3	Stop () const
	Position of the final trajectory point. More...
virtual TVector3	StopDir () const
	Get direction at the stop point. More...
virtual double	TotalLength () const
	Length (cm) of all the track segments. More...
virtual double	DistanceFromStart (double z) const
virtual double	DistanceFromEnd (double z) const
virtual double	W (const rb::CellHit *chit) const
	Estimate the unmeasured coordinate of chit. More...
virtual void	InterpolateXY (double z, double &x, double &y) const
virtual TVector3	InterpolateDir (double z) const
void	ClearTrajectoryPoints ()
	Forget about all trajectory points. More...
void	RemoveTrajectoryPoint (unsigned int i)
	Remove the ith trajectory point from the track. More...
rb::Track	ZipWith (const rb::Track &trk) const
	Combine with a 2D track from the other view to make a 3D track. More...
std::map< unsigned int, TVector3 >	PlaneDirMap () const
	map of the direction cosines at each plane (ie z) More...
virtual bool	Is3D () const
virtual TVector3	Start () const
virtual TVector3	Dir () const
	Unit vector describing prong direction. More...
virtual void	SetDir (TVector3 dir)
virtual void	SetDir (double dv, double dz)
bool	operator< (const Prong &other) const
bool	operator< (const Cluster &other) const
virtual void	Add (const art::Ptr< rb::CellHit > &cell, double weight=1)
virtual void	Add (const art::PtrVector< rb::CellHit > &cells, const std::vector< double > &weights=std::vector< double >())
void	SetID (int id)
const int	ID () const
bool	Is2D () const
virtual geo::View_t	View () const
	kXorY for 3D clusters. More...
rb::Cluster	Exclude (const rb::Cluster *excl) const
	Create a cluster from this one, but with the hits of excl removed. More...
rb::RecoHit	RecoHit (const art::Ptr< rb::CellHit > &chit) const
	Return calibrated hit based on assumed W coordinate. More...
rb::RecoHit	RecoHit (geo::View_t view, unsigned int viewIdx) const
rb::RecoHit	RecoHit (unsigned int globalIdx) const
virtual void	Clear ()
	Forget about all owned cell hits. More...
void	RemoveHit (const art::Ptr< rb::CellHit > hit)
	Remove hit from current cluster. More...
void	RemoveHit (unsigned int globalIdx)
	Remove the ith hit from current cluster. More...
void	SetNoise (bool noise)
	Declare the cluster to consist of noise hits or not. More...
bool	IsNoise () const
	Is the noise flag set? More...
void	SetWeight (unsigned int globalIdx, double weight)
	Set weight of the cell at this index. More...
void	SetWeight (geo::View_t view, unsigned int viewIdx, double weight)
	Set weight of the cell at this index in this view. More...
void	StandardSort ()
	Put the cells in the cluster into a standard order. More...
void	SavePrecalcTotalGeV (ESaveGeVMode savemode)
	Store the current result of TotalGeV / CalorimetricEnergy. More...
Constructors
See rb::Prong documentation
	Track ()
	Track (const art::PtrVector< rb::CellHit > &cells, TVector3 start, TVector3 dir, int id=0)
	Track (const art::PtrVector< rb::CellHit > &cells, geo::View_t view, double v0, double z0, double dv, double dz, int id=0)
	Track (const std::vector< art::Ptr< rb::CellHit > > &cells, geo::View_t view, double v0, double z0, double dv, double dz, int id=0)
	Track (const rb::Cluster &clust, TVector3 start, TVector3 dir, int id=0)
	Track (const rb::Cluster &clust, double v0, double z0, double dv, double dz, int id=0)
	Track (const rb::Prong &prong, int id=0)
	Track (art::PtrVector< rb::CellHit > &cells, std::vector< double > const &weights, std::vector< TVector3 > const &trajectory, TVector3 const &start, TVector3 const &dir, int const &id, geo::View_t view)
Cell counts
unsigned int	NCell (geo::View_t view) const
	Number of cells in view view. More...
unsigned int	NCell () const
	Number of cells in either view. More...
unsigned int	NXCell () const
	Number of cells in the x-view. More...
unsigned int	NYCell () const
	Number of cells in the y-view. More...
Cell Accessors
art::Ptr< rb::CellHit >	Cell (geo::View_t view, unsigned int viewIdx) const
	Get the ith cell from view view. More...
art::Ptr< rb::CellHit >	Cell (unsigned int globalIdx) const
	Get the ith cell from either view. More...
art::Ptr< rb::CellHit >	XCell (unsigned int xIdx) const
	Get the ith cell in the x-view. More...
art::Ptr< rb::CellHit >	YCell (unsigned int yIdx) const
	Get the ith cell in the y-view. More...
const art::PtrVector< rb::CellHit > &	XCells () const
	Get all cells from the x-view. More...
const art::PtrVector< rb::CellHit > &	YCells () const
	Get all cells from the x-view. More...
art::PtrVector< rb::CellHit >	AllCells () const
	Get all cells from both views. More...
std::vector< geo::OfflineChan >	OfflineChans () const
	Positions of all the CellHits. More...
std::vector< rb::WeightedHit >	WeightedHits () const
	Get all hits from both views, with weights attached. More...
double	Weight (unsigned int globalIdx) const
	Weight assigned to the cell. More...
double	Weight (geo::View_t view, unsigned int viewIdx) const
	Weight assigned to the cell. More...
Minima
The "lower" corner of a box containing all the hits
TVector3	MinXYZ () const
double	MinV (geo::View_t view) const
double	MinX () const
double	MinY () const
double	MinZ () const
double	MinTNS () const
unsigned int	MinPlane (geo::View_t view=geo::kXorY) const
unsigned int	MinCell (geo::View_t view) const
Maxima
The "upper" corner of a box containing all the hits
TVector3	MaxXYZ () const
double	MaxV (geo::View_t view) const
double	MaxX () const
double	MaxY () const
double	MaxZ () const
double	MaxTNS () const
unsigned int	MaxPlane (geo::View_t view=geo::kXorY) const
unsigned int	MaxCell (geo::View_t view) const
Means
The unweighted mean position and time of all the hits in the cluster
TVector3	MeanXYZ (rb::AveragingScheme=kDefaultScheme) const
double	MeanV (geo::View_t view, rb::AveragingScheme scheme=kDefaultScheme) const
double	MeanX (rb::AveragingScheme scheme=kDefaultScheme) const
double	MeanY (rb::AveragingScheme scheme=kDefaultScheme) const
double	MeanZ (rb::AveragingScheme scheme=kDefaultScheme) const
double	MeanTNS (rb::AveragingScheme scheme=kDefaultScheme) const
void	MinMaxMeanXYZ (TVector3 &lo, TVector3 &hi, TVector3 &mean, rb::AveragingScheme scheme=kDefaultScheme) const
	Gets the min/max/mean all at once, called by the functions above. More...
Extents
The size of a box containing all the hits
TVector3	ExtentXYZ () const
double	ExtentV (geo::View_t view) const
double	ExtentX () const
double	ExtentY () const
double	ExtentZ () const
unsigned int	ExtentPlane (geo::View_t view=geo::kXorY) const
unsigned int	ExtentCell (geo::View_t view) const
double	ExtentTNS () const
Connectedness measures
int	MostContiguousPlanes (geo::View_t view) const
	Longest run of adjacent planes with hits. More...
int	MostMissingPlanes (geo::View_t view) const
	Longest run of adjacent planes with no hits. More...
int	NMissingPlanes (geo::View_t view) const
	Total number of missing planes in cluster. More...

Protected Member Functions
void	FindNeighbouringPoints (double z, TVector3 &pt0, TVector3 &pt1) const
	Helper function for W. More...
void	FindNeighbouringPointIndices (double z, int &i0, int &i1) const
void	InterpolatePts (double z, const TVector3 &a, const TVector3 &b, double &x, double &y) const
	Helper function for InterpolateXY and ZipWith. More...
Internal helpers
void	EnsureWeightAlloc ()
	Helper. Resizes weights vectors to match cell vectors. More...

Protected Attributes
std::vector< TVector3 >	fTraj
TVector3	fStart
	Start location (xyz, cm) More...
TVector3	fDir
	Direction at starting point. More...
geo::View_t	fView
	view this cluster is in More...
art::PtrVector< rb::CellHit >	fXCell
	collection of x-view cells in cluster More...
art::PtrVector< rb::CellHit >	fYCell
	collection of y-view cells in cluster More...
std::vector< double >	fXWeights
	Weights, matching cell indexing. More...
std::vector< double >	fYWeights
	May be empty, means all weights are 1. More...
int	fID
	ID for cluster. More...
bool	fNoiseCluster
	flag for whether this is a noise cluster More...
double	fPrecalcTotalGeV
	-1 = uninitialized More...

Totals
double	TotalADC () const
	Sum of the ADC of all the contained hits. More...
double	TotalPE () const
	Sum of the PE value of all the contained hits. More...
double	TotalGeV (EEnergyCalcScheme escheme=kRecomputeEnergy) const
	Simple sum of the estimated GeV of all the hits. More...
double	TotalWeight () const
	Sum of all the weights. The effective number of hits. More...
double	CalorimetricEnergy (EEnergyCalcScheme escheme=kRecomputeEnergy) const
	Simple estimate of neutrino energy. More...
enum	EEnergyCalcScheme { kRecomputeEnergy = 0, kUsePrecalcEnergy = 1 }

Detailed Description

A rb::Prong with full reconstructed trajectory.

Definition at line 20 of file Track.h.

Member Enumeration Documentation

enum rb::Cluster::EEnergyCalcScheme

inherited

Enumerator
kRecomputeEnergy	Default. Ask Calibrator about each hit, and sum.
kUsePrecalcEnergy	Use a value computed by a previous call to SavePrecalcTotalGeV.

Definition at line 177 of file Cluster.h.

                          {
      kRecomputeEnergy = 0, ///< Default. Ask Calibrator about each hit, and sum
      kUsePrecalcEnergy = 1 ///< Use a value computed by a previous call to \ref SavePrecalcTotalGeV
    };

enum rb::Cluster::ESaveGeVMode

inherited

Enumerator
kInitializeTotalGeV	Default, this is the initial/only call.
kResetTotalGeV	Altering the already-set value is not an error.

Definition at line 281 of file Cluster.h.

                     {
      kInitializeTotalGeV = 0, ///< Default, this is the initial/only call
      kResetTotalGeV = 1 ///< Altering the already-set value is not an error
    };

Constructor & Destructor Documentation

rb::Track::Track ( )

inline

Definition at line 27 of file Track.h.

References AppendTrajectoryPoint(), dir, dz, PrependTrajectoryPoint(), gen_hdf5record::pt, febshutoff_auto::start, POTSpillRate::view, weights, and ~Track().

{}

rb::Track::Track	(	const art::PtrVector< rb::CellHit > &	cells,
		TVector3	start,
		TVector3	dir,
		int	id = 0
	)

Definition at line 22 of file Track.cxx.

References AppendTrajectoryPoint().

    : rb::Prong(cells, start, dir, id)
  {
    AppendTrajectoryPoint(start);
  }

rb::Track::Track	(	const art::PtrVector< rb::CellHit > &	cells,
		geo::View_t	view,
		double	v0,
		double	z0,
		double	dv,
		double	dz,
		int	id = 0
	)

Definition at line 32 of file Track.cxx.

References AppendTrajectoryPoint().

    : rb::Prong(cells, view, v0, z0, dv, dz, id)
  {
    AppendTrajectoryPoint(v0, z0);
  }

rb::Track::Track	(	const std::vector< art::Ptr< rb::CellHit > > &	cells,
		geo::View_t	view,
		double	v0,
		double	z0,
		double	dv,
		double	dz,
		int	id = 0
	)

Definition at line 45 of file Track.cxx.

References AppendTrajectoryPoint().

    : rb::Prong(cells, view, v0, z0, dv, dz, id)
  {
    AppendTrajectoryPoint(v0, z0);
  }

rb::Track::Track	(	const rb::Cluster &	clust,
		TVector3	start,
		TVector3	dir,
		int	id = 0
	)

Definition at line 58 of file Track.cxx.

References AppendTrajectoryPoint().

    : rb::Prong(clust, start, dir, id)
  {
    AppendTrajectoryPoint(start);
  }

rb::Track::Track	(	const rb::Cluster &	clust,
		double	v0,
		double	z0,
		double	dv,
		double	dz,
		int	id = 0
	)

Definition at line 68 of file Track.cxx.

References AppendTrajectoryPoint().

    : rb::Prong(clust, v0, z0, dv, dz, id)
  {
    AppendTrajectoryPoint(v0, z0);
  }

rb::Track::Track	(	const rb::Prong &	prong,
		int	id = 0
	)

Definition at line 80 of file Track.cxx.

References AppendTrajectoryPoint(), rb::Cluster::SetID(), and rb::Prong::Start().

    : rb::Prong(prong)
  {
    SetID(id);
    AppendTrajectoryPoint(Start());
  }

rb::Track::Track	(	art::PtrVector< rb::CellHit > &	cells,
		std::vector< double > const &	weights,
		std::vector< TVector3 > const &	trajectory,
		TVector3 const &	start,
		TVector3 const &	dir,
		int const &	id,
		geo::View_t	view
	)

Definition at line 89 of file Track.cxx.

References rb::Cluster::fView, rb::Cluster::SetWeight(), POTSpillRate::view, and w.

  : Prong(cells, start, dir, id)
  , fTraj(trajectory)
  {
    fView = view;
    // set the cell weights
    for(size_t w = 0; w < weights.size(); ++w) SetWeight(w, weights[w]);

    return;
  }

rb::Track::~Track ( )

virtual

Definition at line 107 of file Track.cxx.

Referenced by Track().

  {
  }

Member Function Documentation

void rb::Cluster::Add ( const art::Ptr< rb::CellHit > &  cell, double  weight = 1  )

virtualinherited

Definition at line 84 of file Cluster.cxx.

References ana::assert(), rb::Cluster::EnsureWeightAlloc(), rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fView, rb::Cluster::fXCell, rb::Cluster::fXWeights, rb::Cluster::fYCell, rb::Cluster::fYWeights, geo::kX, geo::kY, art::PtrVector< T >::push_back(), and rb::CellHit::View().

Referenced by rb::Cluster::Add(), fuzz::FuzzyKVertex::AddProng(), bpfit::BPFTmvaTrainer::analyze(), bpf::BPFCVNAna::analyze(), htk::HoughTrack::analyze(), rb::Cluster::Cluster(), ddthelpers::cluster_hits(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), tdslicer::TDSlicer::DBSCAN(), slicemergeviews::SliceMergeViews::DoMerge(), rb::Cluster::Exclude(), comi::NumiFiltering::filter(), mono::MonopoleCluster::filter(), bsf::BremShowerFilter::filter(), mono::MonopoleTrack::filter(), mono::MonopoleCluster::find_isolated_hits(), trk::KalmanTrack::FindTracks(), cosrej::CosRejFxs::getActivity(), numusand::NumuSandFxs::getActivityVtx(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackDedx(), trk::KalmanTrackMerge::JoinTracks(), trk::WindowTrackingAlg::LookForBremsstrahlungHits(), mono::Cluster::make(), me::MEFinder::MakeSlcME(), fuzz::ViewMatchAlg::MakeTestProng(), nerd::ViewMatchAlg::MakeTestProng(), trk::WindowTrackingAlg::MakeTrack(), airshower::AirSlicer::MakeTrackSlices(), me::MEFinder::MakeTrkME(), numue::NumuEAlg::MakeVertexCluster(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), trk::KalmanTrackMerge::MatchTracks(), cheat::BackTracker::MCTruthToCluster(), qeef::QeFinder::produce(), slicer::Slicer4D::produce(), nuesand::FillNueSandbox::produce(), slid::Recluster::produce(), numusand::FillSandbox::produce(), cvn::CVNCosmicMapper::produce(), zcl::SPCluster::produce(), airshower::AirSlicer::produce(), zcl::SMMCluster::produce(), sn::SNSlicer::produce(), sn::SupernovaMCCluster::produce(), zcl::FMMTracker::produce(), cvn::CVNMapper::produce(), dif::DiFShowerFinder::produce(), bpfit::BPFEnergyEstimator::produce(), bpfit::BPFEnergyEstimatorOnly::produce(), tdslicer::TDSlicer::produce(), trk::WindowTrackingAlg::ShortTrack(), trk::WindowTrackingAlg::ShortTrackExtraPlane(), airshower::AirSlicer::TemporalClusters(), dt::Cand::ToTrack(), htk::Track3D::Track3D(), mono::Track3D::Track3D(), zcl::FMMTracker::Tracking2D(), slicemergeviews::SliceMergeViews::TrimNoise(), tdslicer::TDSlicer::TrimNoise(), trk::KalmanTrackMerge::ViewMergeTracks(), slicer::MMSlicer::WindowSlice(), and slicer::Slicer::WindowSlice().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    assert(weight >= 0 && weight <= 1);

    // Are we optimizing (omitting) the weights storage?
    const bool optWeights = (weight == 1 && fXWeights.empty() && fYWeights.empty());
    if(!optWeights) EnsureWeightAlloc();

    if(cell->View() == geo::kX){
      assert(fView != geo::kY);
      fXCell.push_back(cell);
      if(!optWeights) fXWeights.push_back(weight);
    }
    else{
      assert(fView != geo::kX);
      fYCell.push_back(cell);
      if(!optWeights) fYWeights.push_back(weight);
    }
  }

void rb::Cluster::Add ( const art::PtrVector< rb::CellHit > &  cells, const std::vector< double > &  weights = std::vector<double>()  )

virtualinherited

Definition at line 107 of file Cluster.cxx.

References rb::Cluster::Add(), ana::assert(), art::PtrVector< T >::empty(), getGoodRuns4SAM::n, and art::PtrVector< T >::size().

  {
    const unsigned int N = cells.size();

    assert(weights.empty() || weights.size() == N);

    for(unsigned int n = 0; n < N; ++n){
      assert(weights.empty() || (weights[n] >= 0 && weights[n] <= 1));
      Add(cells[n], weights.empty() ? 1 : weights[n]);
    }
  }

art::PtrVector< rb::CellHit > rb::Cluster::AllCells ( ) const

inherited

Get all cells from both views.

Definition at line 180 of file Cluster.cxx.

References plot_validation_datamc::c, rb::Cluster::fXCell, rb::Cluster::fYCell, art::PtrVector< T >::push_back(), and art::PtrVector< T >::size().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), qeef::QeFinderVal::analyze(), calib::TestParticleCorrections::analyze(), murem::MuonRemoveAna::analyze(), sn::SNSlicerAna::analyze(), validation::NoiseClusterAnalysis::analyze(), calib::HitEfficiency::analyze(), showere::ShowerEnergyAna::analyze(), hough::HoughValidate::analyze(), zcl::FastMMStudy::analyze(), upmuana::UpMuRecoAna::analyze(), upmuana::UpMuAnalysis::analyze(), htk::HoughTrack::analyze(), murem::MRCCAna::analyze(), trk::KalmanTrackAna::CheckAssociations(), air::AirKalmanAna::CheckAssociations(), cheat::RecoCheckAna::CheckRecoClusters(), trk::KalmanTrackAna::CheckRecoTracks(), trk::WindowTrackingAlg::CheckTrackDirectionInY(), bpfit::dEdxCalculator::computeDEDX(), skim::SkimmingUtils::CopySlice(), murem::TrackCleanUpAlg::DeDxInPlane(), nnbar::NNbarUtilities::energyBalancingVertex(), dt::ViewMerger::FakeThirdDimension(), caf::FillDiFShowerVars(), caf::FillDiFVars(), calib::StopperThreshold::FillHist(), calib::CosmicTrends::FillHistograms(), caf::FillProngVars(), caf::FillShowerVars(), CerenkovSelection::FillTree(), calib::StopperThreshold::FillTree(), trk::CosmicTrackAna::FillTrueInfo(), showere::ShowerEnergyFilterMC::filter(), sn::SNSlicer::FinalizeCluster(), cvn::Get2DProngEfficiencyByPDG(), cvn::Get2DProngEnergyByPDG(), cvn::Get2DProngPurityByPDG(), nnbar::NNbarUtilities::getCoordinateOfHits(), nuesand::FillNueSandbox::GetdEdx(), numusand::NumuSandFxs::getMissingE(), nnbar::NNbarUtilities::getMultipleCellFraction(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), cvn::GetProngEfficiencyByPDG(), cvn::GetProngEnergyByPDG(), cvn::GetProngPurityByPDG(), nnbar::NNbarUtilities::getSliceBoundary(), nnbar::NNbarUtilities::getSliceTimeWindow(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackDedx(), nnbar::NNbarUtilities::getXyAsymmetry(), nnbar::NNbarUtilities::isFebFlasher(), novaddt::HorizontalMuonAna::IsHorzMuonTrack(), caf::IsTrueDiF(), calib::RockMuonStopperSelection::IsTrulyContainedStopper(), calib::StopperSelection::IsTrulyContainedStopper(), calib::RockMuonStopperSelection::IsTrulyInFiducial(), calib::ThroughgoingSelection::IsTrulyUncontained(), trk::KalmanTrackMerge::JoinTracks(), dt::ViewMerger::JoinTracks(), trk::WindowTrackingAlg::MakeTrack(), airshower::AirSlicer::MakeTrackSlices(), murem::FindParent::MatchToOrigSlice(), numue::NumuEAlg::MCTruthEnergyVariables(), murem::MuonRemove::MuonInfo(), mono::Monopole::n_mc_hits(), dq::FlasherFinder::Nflash(), zcl::FastMonopoleTriggers::NumberOfCellsPerLength(), zcl::FastMMStudy::NumberOfCellsPerLength(), zcl::FastMonopoleTriggers::NumberOfHitsInOverlapPlanesCut(), zcl::FastMMStudy::NumberOfHitsInOverlapPlanesCut(), mono::operator<<(), skim::ParametersNumu::ParametersNumu(), calib::PCHitsList::ProcessTrackForBelowThresholdHits(), calib::PCHitsList::ProcessTrackForFLSHits(), trk::TrkAssn::produce(), slid::Recluster::produce(), numusand::FillSandbox::produce(), hough::MultiHoughT::produce(), cosrej::CosRej::produce(), numue::NumuEnergy::produce(), earms::ElasticArmsHS::produce(), ndreco::NDReco::produce(), vdt::VertexDT::produce(), upmuana::UpMuProb::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), caf::CAFMaker::produce(), cvn::ProngClassify(), sn::SliceRemover::remove(), murem::MuonRemove::RemoveByEfficiency(), trk::WindowTrackingAlg::SetTrackEndPoints(), trk::WindowTrackingAlg::ShortTrack(), cheat::BackTracker::SliceToMCTruth(), cheat::BackTracker::SliceToNeutrinoInteractions(), cosrej::MakeNueCosRej::SparsenessAsymmetry(), zcl::FastMonopoleTriggers::StdevCellsPerPlane(), zcl::FastMMStudy::StdevCellsPerPlane(), htk::Track3D::time_fit(), mono::Track3D::time_fit(), htk::Track3D::Track3D(), mono::Track3D::Track3D(), slid::NuEEnergyAlg::VertexEnergy(), zcl::FastMonopoleTriggers::WeightedCenterCut(), zcl::FastMMStudy::WeightedCenterCut(), moonshadowana::MoonShadowAna::WriteTrackInfo(), moonshadowana::MoonShadowAnaHough::WriteTrackInfo(), rb::Cluster::YCells(), dt::ViewMerger::ZipVerticalTracks(), and ZipWith().

  {
    art::PtrVector<rb::CellHit> all;
    for(size_t c = 0; c < fXCell.size(); ++c) all.push_back(fXCell[c]);
    for(size_t c = 0; c < fYCell.size(); ++c) all.push_back(fYCell[c]);

    return all;
  }

void rb::Track::AppendTrajectoryPoint

(

TVector3

pt

)

Definition at line 112 of file Track.cxx.

References rb::Cluster::fPrecalcTotalGeV, and fTraj.

Referenced by trk::WindowTrackingAlg::FindEndPoint(), trk::KalmanTrack::FindTracks(), bpfit::BreakPoint::FitTracks(), trk::CosmicTrackAlg::MakeTrack(), trk::WindowTrackingAlg::MakeViewTrack(), trk::KalmanTrackMerge::MatchTracks(), nuesand::FillNueSandbox::produce(), jmshower::JMTrackMerge::produce(), jmshower::JMClusterMerge::produce(), zcl::FMMTracker::produce(), cheat::MCCheater::produce(), trk::WindowTrackingAlg::ShortTrack(), trk::WindowTrackingAlg::ShortViewTrack(), dt::Chain::ToTrack(), dt::Cand::ToTrack(), Track(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    fTraj.push_back(pt);
  }

void rb::Track::AppendTrajectoryPoint	(	double	pv,
		double	pz
	)

Definition at line 120 of file Track.cxx.

References rb::Cluster::fPrecalcTotalGeV, fTraj, rb::Cluster::fView, geo::kX, and geo::kXorY.

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    if(fView == geo::kXorY)
      throw cet::exception("Track") << "view is geo::kXorY, cannot add a 2D trajectroy point";

    if(fView == geo::kX)
      fTraj.push_back(TVector3(pv, 0, pz));
    else
      fTraj.push_back(TVector3(0, pv, pz));
  }

double rb::Cluster::CalorimetricEnergy ( EEnergyCalcScheme  escheme = kRecomputeEnergy ) const

inherited

Simple estimate of neutrino energy.

Uses a simple scale factor to correct for dead material and thresholds, could be more advanced, but this is simple enough to use for preselection etc.

Definition at line 439 of file Cluster.cxx.

References calib::Calibrator::GetGeVToCalorimetricScale(), and rb::Cluster::TotalGeV().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), tbem::TBRecoAna::analyze(), ncid::CompareByCalEnergy(), cosrej::CompareByCalEnergy(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillHadClustVars(), caf::FillProngVars(), caf::FillShowerVars(), caf::FillSliceVars(), xsrec::MakeXSecCCPi0Inc::GetPhLL(), numue::NumuEAlg::GetUCMuonENonSingle(), numue::NumuEAlg::GetUCMuonESingle(), SliceLID::BasicProngVars::initBasicVars(), LSTME::BasicProngVars::initBasicVars(), me::MEFinder::MakeSlcME(), me::MEFinder::MakeTrkME(), skim::ParametersNue::ParametersNue(), skim::ParametersNumu::ParametersNumu(), me::MEFinder::PassesMEPresel(), presel::LEMPresel::produce(), numusand::FillSandbox::produce(), presel::NuePresel::produce(), cosrej::CosRej::produce(), cosrej::MakeNueCosRej::produce(), caf::CAFMaker::produce(), showerByEnergyNue(), showerByEnergyNumu(), lem::LEMSummarizer::SliceToLEMInput(), and ncid::NCNNKeras::TransMomFraction().

  {
    art::ServiceHandle<calib::Calibrator> cal;
    return cal->GetGeVToCalorimetricScale()*TotalGeV(escheme);
  }

art::Ptr< rb::CellHit > rb::Cluster::Cell ( geo::View_t  view, unsigned int  viewIdx  ) const

inherited

Get the ith cell from view view.

Definition at line 145 of file Cluster.cxx.

References ana::assert(), geo::kX, geo::kXorY, geo::kY, rb::Cluster::XCell(), and rb::Cluster::YCell().

Referenced by dqsf::DQSpillFlags::AccumulateSyncMetric(), rb::HitMap::Add(), fuzz::FuzzyKVertex::AddProng(), dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), tut::TutAnalyzer::analyze(), qeef::QeFinderVal::analyze(), rsim::RecordNoiseSpectrumFile::analyze(), chaninfo::CosmicEff::analyze(), remid::ReMIdValidate::analyze(), comi::CosmicMetrics::analyze(), calib::TestParticleCorrections::analyze(), dprf::ChannelPlots::analyze(), remid::ReMIdTrain::analyze(), mcchk::MonopoleAna::analyze(), dprf::TrackPlots::analyze(), align::Alignment::analyze(), mcchk::RockAna::analyze(), align::SplitTracks::analyze(), murem::MuonRemoveAna::analyze(), pa::ClusterEva::analyze(), vf::TrackEva::analyze(), comi::NumiFilteringAna::analyze(), zcl::FmmTriggerAna::analyze(), tbem::TBRecoAna::analyze(), zcl::FmmTrackerValidation::analyze(), zcl::FmmTrackerAna::analyze(), showere::ShowerEnergyAna::analyze(), bpfit::BPFTmvaTrainer::analyze(), bpf::BPFCVNAna::analyze(), calib::CosmicTrends::analyze(), calib::CalibAna::analyze(), ncs::Xeff::analyze(), ncs::ROCKMRE::analyze(), ncs::NCAna::analyze(), bpfit::DimuonFitter::AnaSlice(), trident::DimuonSliceAna::AnaSlice(), me::MEFinder::BestSlcMatch(), fuzz::ViewMatchAlg::CalcEnergyProfile(), nerd::ViewMatchAlg::CalcEnergyProfile(), trk::KalmanTrackMerge::CalcMatchScore(), slid::ParticleIDAlg::CalcPlaneHits(), trk::KalmanTrackMerge::CanJoinTracks(), slid::NuEEnergyAlg::CellEnergy(), murem::TrackCleanUpAlg::CleanUpTrack(), bpfit::dEdxCalculator::computeDEDX(), CerenkovSelection::Contained(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), cvn::RegPixelMapProducer::CreateMapGivenBoundary(), cvn::PixelMapProducer::CreateMapGivenBoundary(), cvn::RegPixelMapProducer::CreateMapGivenShowerVertex(), cvn::RegPixelMapProducer::CreateMapGivenVertex(), jmshower::RecoJMShower::DepositEnergy(), evd::RecoBaseDrawer::DrawTrack2D(), skim::ParametersNumu::EarliestLatestHitPos(), jmshower::RecoJMShower::Energy(), dif::DiFShowerFinder::eparm(), rb::Cluster::Exclude(), dt::ViewMerger::Extremes(), caf::FillHadClustVars(), calib::StopperThreshold::FillHist(), bpfit::BreakPoint::FillHitList(), sn::SupernovaAna::FillPlaneSeperations(), dqsf::DQSpillFlags::FillPopulatedDCMs(), caf::FillProngVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillSliceInfo(), caf::FillSliceVars(), trk::CosmicTrackAna::FillTrackHistograms(), CerenkovSelection::FillTree(), calib::StopperThreshold::FillTree(), comi::NumiFiltering::filter(), bsf::BremShowerFilter::filter(), showere::ShowerEnergyFilterMC::filter(), cvn::RegPixelMapProducer::FindCenterMedian(), cvn::PixelMapProducer::FindCenterMedian(), airshower::AirSlicer::FindNeighbors(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), trk::KalmanTrack::FindTracks(), trk::CosmicTrackUtilities::FindTriCells(), trident::DimuonSliceAna::FindVertex(), cvn::RegPixelMapProducer::FindVertexMaxGap(), cvn::PixelMapProducer::FindVertexMaxGap(), cvn::RegPixelMapProducer::FindVertexWindowThreshold(), cvn::PixelMapProducer::FindVertexWindowThreshold(), bpfit::DimuonFitter::FindVertexZ(), bpfit::BreakPoint::FitTracks(), bpfit::DimuonFitter::FitView(), nuesand::FillNueSandbox::FracAngChanges(), nuesand::FillNueSandbox::FracModalHits(), nuesand::FillNueSandbox::FracNonContiguous(), numusand::NumuSandFxs::getActivity(), numusand::NumuSandFxs::getAveTrackdEdx(), numusand::NumuSandFxs::getAveTrackdEdxLast4Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast6Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast8Cells(), trackinfo::TrackInfoFxs::getAveTrackdEdxLastCells(), cosrej::CosRejFxs::getBBC(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), evd::RecoBaseDrawer::GetClusterOutlines(), evd::RecoBaseDrawer::GetClusterOutlinesConvexHull(), nuesand::FillNueSandbox::GetdEdx(), cosrej::CosRejFxs::getFits(), jmshower::RecoJMShower::GetPlaneCentroidCell(), jmshower::RecoJMShower::GetPlaneE1Cell(), jmshower::JMShower::GetPlaneHits(), jmshower::RecoJMShower::GetRVPStats(), jmshower::JMShower::GetTransHits(), calib::PCHitsList::GetTrueEnergyPathAndLightForCell(), slid::DeconvolveAlg::GetWeights(), slid::NuEEnergyAlg::HadronicDepEnergy(), rb::HitList::HitList(), cvn::HitNuIndex(), tf::TimingFitAlg::HoughFit(), zcl::FastMonopoleTriggers::is_trigger_by_epoch1_fmmtrigger(), trk::CosmicTrackSelection::IsCalibTrack(), calib::RockMuonStopperSelection::IsContainedStopper(), calib::PCHitsList::IsGoodTrack(), trk::CosmicTrackAlg::IsTrackDownstreamFromTiming(), trk::KalmanTrackMerge::JoinTracks(), fuzz::ViewMatchAlg::MakeTestProng(), nerd::ViewMatchAlg::MakeTestProng(), trk::CosmicTrackAlg::MakeTrack(), me::MEFinder::MakeTrkME(), murem::FindParent::MatchToOrigSlice(), rb::Cluster::MaxCell(), nuesand::FillNueSandbox::MaxGap(), nuesand::FillNueSandbox::MaxHits(), rb::Cluster::MaxPlane(), rb::Cluster::MaxTNS(), rb::Cluster::MeanTNS(), dt::DiscreteTracker::MergeParasiteTracks(), rb::Cluster::MinCell(), nerd::NERDProng::MinHitClusterDist(), me::MEFinder::MinHitSlcDist(), rb::Cluster::MinMaxMeanXYZ(), rb::Cluster::MinPlane(), rb::Cluster::MinTNS(), rb::Cluster::MostContiguousPlanes(), rb::Cluster::MostMissingPlanes(), murem::MuonRemove::MuonInfo(), dif::DiFShowerFinder::muonstub(), rb::Cluster::NCell(), numue::NumuEAlg::NDEnergy(), numue::NumuEAlg::NDTrackEnergySplitting(), rb::Cluster::NMissingPlanes(), zcl::FmmTriggerEvd::NumberOfCellsPerLength(), zcl::FmmTriggerEvd::NumberOfHitsInOverlapPlanesCut(), zcl::FastMonopoleTriggers::NumberOfSurfaceHits(), zcl::FastMMStudy::NumberOfSurfaceHits(), rb::Cluster::OfflineChans(), skim::ParametersNumu::ParametersNumu(), zcl::FastMonopoleTriggers::passed_epoch2_fmmtrigger(), trident::DimuonSliceAna::PlaneCellCounts(), PlaneDirMap(), slid::ParticleIDAlg::PlaneRadius(), calib::PCHitsList::ProcessTrack(), calib::PCHitsList::ProcessTrackTrajectory(), trk::TrkAssn::produce(), vf::VertexFinder::produce(), qeef::QeFinder::produce(), hv::HoughVertex::produce(), nuesand::FillNueSandbox::produce(), slid::Recluster::produce(), lem::Preselection::produce(), presel::NuePresel::produce(), numusand::FillSandbox::produce(), slid::SliceLIDBuilder::produce(), zcl::SPCluster::produce(), earms::ElasticArmsHS::produce(), lem::MakeLibrary::produce(), dt::DiscreteTracker::produce(), zcl::SMMCluster::produce(), zcl::FMMTracker::produce(), zcl::FmmTriggerEvd::produce(), jmshower::NueSel::produce(), rvp::RecVarPID::produce(), dif::DiFShowerFinder::produce(), jmshower::NueSelLID::produce(), bpfit::BPFEnergyEstimator::produce(), bpfit::BPFEnergyEstimatorOnly::produce(), caf::CAFMaker::produce(), slid::ParticleIDAlg::Radius(), jmshower::RecoJMShower::Radius(), jmshower::RecoJMShower::RadiusV(), rb::Cluster::RecoHit(), jmshower::RecoJMShower::RecoShowers(), murem::MuonRemove::RemoveByTruth(), trk::KalmanTrack::RemoveHitsFromSignal(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), lem::MakeLibrary::SliceMeanPosEWeighted(), lem::MakeLibrary::SliceOrigin(), lem::LEMSummarizer::SliceToLEMInput(), zcl::FmmTriggerEvd::StdevCellsPerPlane(), calib::StopperThreshold::testPath(), rb::Cluster::TotalADC(), dt::ViewMerger::TotalChargePerView(), rb::Prong::TotalLength(), rb::Cluster::TotalPE(), trk::CosmicTrackUtilities::TrackResiduals(), trk::KalmanTrackMerge::ViewMergeTracks(), rb::Prong::W(), W(), and rb::Cluster::WeightedHits().

  {
    switch(view){
    case geo::kX: return XCell(viewIdx);
    case geo::kY: return YCell(viewIdx);
    case geo::kXorY: return Cell(viewIdx);
    default: assert(0 && "Unknown view");
    }
  }

art::Ptr< rb::CellHit > rb::Cluster::Cell ( unsigned int  globalIdx ) const

inherited

Get the ith cell from either view.

Definition at line 173 of file Cluster.cxx.

References rb::Cluster::NXCell(), rb::Cluster::XCell(), and rb::Cluster::YCell().

  {
    if(globalIdx < NXCell()) return XCell(globalIdx);
    return YCell(globalIdx-NXCell());
  }

void rb::Cluster::Clear ( )

virtualinherited

Forget about all owned cell hits.

Virtual function because derived classes may have more stuff to clear

Definition at line 279 of file Cluster.cxx.

References art::PtrVector< T >::clear(), rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fXCell, rb::Cluster::fXWeights, rb::Cluster::fYCell, and rb::Cluster::fYWeights.

Referenced by bsf::BremShowerFilter::filter(), mono::Cluster::make(), calib::PCHitsList::ProcessTrackForFLSHits(), nuesand::FillNueSandbox::produce(), cvn::CVNCosmicMapper::produce(), jmshower::RecoJMShower::RecoShowers(), and rb::Cluster::YCells().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    fXCell.clear();
    fYCell.clear();
    fXWeights.clear();
    fYWeights.clear();
  }

void rb::Track::ClearTrajectoryPoints

(

)

Forget about all trajectory points.

Definition at line 359 of file Track.cxx.

References rb::Cluster::fPrecalcTotalGeV, and fTraj.

Referenced by trk::KalmanTrack::FindTracks(), Trajectory(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    fTraj.clear();
    fTraj.resize(0);
  }

virtual TVector3 rb::Prong::Dir ( ) const

inlinevirtualinherited

Unit vector describing prong direction.

For 2D prongs the unmeasured component is zero

Definition at line 77 of file Prong.h.

References rb::Prong::DistanceFromStart(), rb::Prong::fDir, rb::Prong::operator<(), fhicl::other, rb::Prong::SetDir(), rb::Prong::SetStart(), rb::Prong::TotalLength(), rb::Prong::W(), and test::z.

Referenced by dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), comi::CosmicMetrics::analyze(), murem::MuonRemoveAna::analyze(), vf::TrackEva::analyze(), calib::HitEfficiency::analyze(), tbem::TBRecoAna::analyze(), trk::KalmanTrackAna::analyze(), bpf::BPFCVNAna::analyze(), calib::CalibAna::analyze(), upmuana::UpMuAnalysis::analyze(), red::EWCosmics::analyze(), slid::ParticleIDAlg::CalcAsymIneria(), slid::ParticleIDAlg::CalcCellPlaneTransverseDedx(), slid::ParticleIDAlg::CalcInterPlaneDedxProb(), slid::ParticleIDAlg::CalcPlaneHits(), slid::ParticleIDAlg::CalcPlaneTransverseDedx(), slid::ParticleIDAlg::CalcTrkHitPath(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxFD::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), murem::TrackCleanUpAlg::DeDxInPlane(), slid::ParticleIDAlg::DedxInverseLongLL(), slid::ParticleIDAlg::DedxLongLL(), dif::DiFShowerFinder::distancefromtrack(), slid::DeconvolveAlg::DistanceToCore(), vf::DOCAInfo::DOCAInfo(), dif::DiFShowerFinder::eparm(), dt::ViewMerger::FakeThirdDimension(), caf::FillDiFShowerVars(), calib::CosmicTrends::FillHistograms(), caf::FillProngVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillShowerVars(), filter::Filter::FillTrackVariables(), trk::CosmicTrackAna::FillTrueInfo(), comi::NumiFiltering::filter(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), bsf::BremShowerFilter::findShowerByTruth(), trk::KalmanTrack::FindTracks(), jmshower::RecoJMShower::GetCellTransDedx(), cosrej::CosRejFxs::getCVVars(), nuesand::FillNueSandbox::GetdEdx(), cosrej::CosRejFxs::getFits(), jmshower::RecoJMShower::GetInterPlaneDedxProb(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), ndreco::NDRecoFxs::getProngProximity(), jmshower::RecoJMShower::GetShwStop(), nuesand::FillNueSandbox::GetTrackAngle(), jmshower::JMShower::GetTransHits(), jmshower::RecoJMShower::GetTrkHitPath(), jmshower::RecoJMShower::GetTrkHitPos(), jmshower::RecoJMShower::GetTrkPlaneCell(), jmshower::RecoJMShower::GetTrkPlaneDistToEdge(), jmshower::RecoJMShower::GetTrkPlanePos(), trk::CosmicTrackSelection::GoodReco(), tf::TimingFitAlg::HoughFit(), LSTME::BasicProngVars::initBasicVars(), trk::CosmicTrackSelection::IsCalibTrack(), calib::PCHitsList::IsGoodTrack(), fuzz::ViewMatchAlg::MakeTestProng(), nerd::ViewMatchAlg::MakeTestProng(), trk::CosmicTrackAlg::MakeTrack(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), murem::MuonRemove::MuonInfo(), cosrej::MakeNueCosRej::MuonParentByDist(), cosrej::MakeNueCosRej::MuonParentByTime(), dif::DiFShowerFinder::muonstub(), skim::ParametersNumu::ParametersNumu(), slid::ParticleIDAlg::PlaneHitCell(), slid::ParticleIDAlg::PlaneHitXYZ(), slid::ParticleIDAlg::PlaneRadius(), wsnumu::WSNumu::produce(), slid::Recluster::produce(), nussand::FillNusSandbox::produce(), presel::CosmicVeto::produce(), trackinfo::TrackInfo::produce(), slid::SliceLIDBuilder::produce(), cosrej::CosRej::produce(), ndreco::NDReco::produce(), remid::RecoMuon::produce(), cosrej::MakeNueCosRej::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), slid::ParticleIDAlg::Radius(), jmshower::RecoJMShower::RecoShowers(), trk::WindowTrackingAlg::SetTrackEndPoints(), trk::WindowTrackingAlg::ShortTrack(), trk::WindowTrackingAlg::ShortTrackExtraPlane(), red::OverburdenCalculator::traceBack(), slid::NuEEnergyAlg::VertexEnergy(), and trk::KalmanTrackMerge::ViewMergeTracks().

{return fDir;}

double rb::Track::DistanceFromEnd ( double  z ) const

virtual

Definition at line 255 of file Track.cxx.

References std::abs(), dist, FindNeighbouringPointIndices(), rb::Prong::fStart, fTraj, MECModelEnuComparisons::i, Unit(), and test::z.

Referenced by trk::CosmicTrackAna::FillTrackHistograms(), cosrej::CosRej::produce(), and Trajectory().

  {
    if(fTraj.size() == 1) return std::abs(z-fStart.Z());

    int i0, i1;
    FindNeighbouringPointIndices(z, i0, i1);

    // Distance from the end the track to the first of the neighbouring points
    double dist = 0;
    for(int i = fTraj.size()-1; i > i1; --i)
      dist += (fTraj[i]-fTraj[i-1]).Mag();

    // The extra accumulated between the points. Also correct for the cases
    // where we're off the start or end of the track.

    // Skip this step if the two trajectory points have exactly the same
    // Z. This would happen if the track were exactly vertical, or more likely
    // if there's a repeated trajectory point somehow.
    const double denom = (fTraj[i1]-fTraj[i0]).Unit().Z();
    if(denom != 0)
      dist += (fTraj[i1].Z()-z) / denom;

    return dist;
  }

double rb::Track::DistanceFromStart ( double  z ) const

virtual

Parameters

z	z-coordinate to evaluate trajectory at

Returns

Distance from Start measured along the trajectory

Reimplemented from rb::Prong.

Definition at line 229 of file Track.cxx.

References std::abs(), dist, FindNeighbouringPointIndices(), rb::Prong::fStart, fTraj, MECModelEnuComparisons::i, and Unit().

Referenced by calib::PCHitsList::CellHitToPCHit(), cosrej::CosRej::produce(), and Trajectory().

  {
    if(fTraj.size() == 1) return std::abs(z-fStart.Z());

    int i0, i1;
    FindNeighbouringPointIndices(z, i0, i1);

    // Distance along the track to the first of the neighbouring points
    double dist = 0;
    for(int i = 0; i < i0; ++i)
      dist += (fTraj[i+1]-fTraj[i]).Mag();

    // The extra accumulated between the points. Also correct for the cases
    // where we're off the start or end of the track.

    // Skip this step if the two trajectory points have exactly the same
    // Z. This would happen if the track were exactly vertical, or more likely
    // if there's a repeated trajectory point somehow.
    const double denom = (fTraj[i1]-fTraj[i0]).Unit().Z();
    if(denom != 0)
      dist += (z-fTraj[i0].Z()) / denom;

    return dist;
  }

void rb::Cluster::EnsureWeightAlloc ( )

protectedinherited

Helper. Resizes weights vectors to match cell vectors.

Definition at line 712 of file Cluster.cxx.

References rb::Cluster::fXCell, rb::Cluster::fXWeights, rb::Cluster::fYCell, rb::Cluster::fYWeights, and art::PtrVector< T >::size().

Referenced by rb::Cluster::Add(), and rb::Cluster::SetWeight().

  {
    if(fXWeights.size() != fXCell.size()) fXWeights.resize(fXCell.size(), 1);
    if(fYWeights.size() != fYCell.size()) fYWeights.resize(fYCell.size(), 1);
  }

rb::Cluster rb::Cluster::Exclude ( const rb::Cluster *  excl ) const

inherited

Create a cluster from this one, but with the hits of excl removed.

Definition at line 233 of file Cluster.cxx.

References rb::Cluster::Add(), rb::Cluster::Cell(), rb::Cluster::fID, rb::Cluster::fNoiseCluster, rb::Cluster::fView, MECModelEnuComparisons::i, calib::j, genie::utils::mec::J(), rb::Cluster::NCell(), runNovaSAM::ret, rb::Cluster::SetNoise(), and rb::Cluster::Weight().

Referenced by mono::MonopoleCluster::filter(), slid::Recluster::produce(), slid::NuEEnergyAlg::VertexEnergy(), and rb::Cluster::YCells().

  {
    rb::Cluster ret(fView, fID);
    ret.SetNoise(fNoiseCluster);

    const unsigned int I = NCell();
    const unsigned int J = excl->NCell();
    for(unsigned int i = 0; i < I; ++i){
      const art::Ptr<rb::CellHit> chit = Cell(i);
      bool clash = false;
      for(unsigned int j = 0; j < J; ++j){
        if(excl->Cell(j) == chit){
          clash = true;
          break;
        }
      } // end for j
      if(!clash) ret.Add(chit, Weight(i));
    }

    if(ret.NCell() == 0)
      mf::LogWarning("Cluster") << "Exclude created empty cluster";

    return ret;
  }

unsigned int rb::Cluster::ExtentCell ( geo::View_t  view ) const

inherited

Definition at line 570 of file Cluster.cxx.

References rb::Cluster::MaxCell(), and rb::Cluster::MinCell().

Referenced by sn::SNSlicerAna::analyze(), sn::SNSlicer::ClusterIsGood(), and rb::Cluster::ExtentPlane().

  {
    return MaxCell(view)-MinCell(view)+1;
  }

unsigned int rb::Cluster::ExtentPlane ( geo::View_t  view = geo::kXorY ) const

inlineinherited

Definition at line 250 of file Cluster.h.

References rb::Cluster::ExtentCell(), rb::Cluster::MaxPlane(), and rb::Cluster::MinPlane().

Referenced by dt::ViewMerger::Adjacent(), align::AlignInspect::analyze(), align::Alignment::analyze(), align::SplitTracks::analyze(), sn::SNSlicerAna::analyze(), calib::PEResponse::analyze(), ncs::Xeff::analyze(), ncs::ROCKMRE::analyze(), ncs::NCAna::analyze(), trk::KalmanTrackMerge::CalcMatchScore(), sn::SNSlicer::ClusterIsGood(), dt::ViewMerger::FakeThirdDimension(), caf::FillDiFShowerVars(), caf::FillDiFVars(), calib::CosmicTrends::FillHistograms(), caf::FillProngVars(), caf::FillShowerVars(), jmshower::RecoJMShower::GetRVPStats(), LSTME::BasicProngVars::initBasicVars(), align::AlignInspect::IsGoodTrack(), align::Alignment::IsGoodTrack(), dt::ViewMerger::MatchScore(), dt::DiscreteTracker::MergeParasiteTracks(), skim::ParametersNue::ParametersNue(), presel::CosmicVeto::produce(), lem::Preselection::produce(), presel::NuePresel::produce(), rvp::RecVarPID::produce(), dprf::ChannelPlots::TrackOK(), dprf::TrackPlots::TrackOK(), and trk::KalmanTrackMerge::ViewMergeTracks().

{return MaxPlane(view)-MinPlane(view)+1;}

double rb::Cluster::ExtentTNS ( ) const

inlineinherited

Definition at line 252 of file Cluster.h.

References rb::Cluster::MaxTNS(), rb::Cluster::MinTNS(), rb::Cluster::MostContiguousPlanes(), rb::Cluster::MostMissingPlanes(), rb::Cluster::NMissingPlanes(), and rb::Cluster::StandardSort().

Referenced by sn::SNSlicerAna::analyze(), comi::DataCheck::analyze(), and sn::SNSlicer::ClusterIsGood().

{return MaxTNS()-MinTNS();}

double rb::Cluster::ExtentV ( geo::View_t  view ) const

inlineinherited

Definition at line 246 of file Cluster.h.

References rb::Cluster::MaxV(), and rb::Cluster::MinV().

{return MaxV(view)-MinV(view);}

double rb::Cluster::ExtentX ( ) const

inlineinherited

Definition at line 247 of file Cluster.h.

References rb::Cluster::ExtentXYZ().

{return ExtentXYZ().X();}

TVector3 rb::Cluster::ExtentXYZ ( ) const

inlineinherited

Definition at line 245 of file Cluster.h.

References rb::Cluster::MaxXYZ(), and rb::Cluster::MinXYZ().

Referenced by rb::Cluster::ExtentX(), rb::Cluster::ExtentY(), and rb::Cluster::ExtentZ().

{return MaxXYZ()-MinXYZ();}

double rb::Cluster::ExtentY ( ) const

inlineinherited

Definition at line 248 of file Cluster.h.

References rb::Cluster::ExtentXYZ().

{return ExtentXYZ().Y();}

double rb::Cluster::ExtentZ ( ) const

inlineinherited

Definition at line 249 of file Cluster.h.

References rb::Cluster::ExtentXYZ().

Referenced by air::AirFilter::parallel_muon_slice().

{return ExtentXYZ().Z();}

void rb::Track::FindNeighbouringPointIndices ( double  z, int &  i0, int &  i1  ) const

protected

Definition at line 380 of file Track.cxx.

References std::abs(), e, modifyFHiCL::found, fTraj, LOG_DEBUG, getGoodRuns4SAM::n, test::z, and Z.

Referenced by DistanceFromEnd(), DistanceFromStart(), FindNeighbouringPoints(), and Trajectory().

  {
    if(fTraj.size() < 2)
      throw cet::exception("Track") << "attempting to find neighboring points with "
                                    << "fewer than 2 trajectory points in track";

    // Track might be going upstream
    const int zDir = (fTraj.front().Z() > fTraj.back().Z()) ? -1 : 1;

    // If hit is off the end of the track extrapolate from the two closest points
    if( zDir*fTraj.front().Z() >= zDir*z ){
      i0 = 0;
      i1 = 1;
      return;
    }
    if( zDir*fTraj.back().Z() <= zDir*z ){
      i0 = fTraj.size()-2;
      i1 = fTraj.size()-1;
      return;
    }

    bool found = false;

    // We could do this more efficiently with std::[upper|lower]_bound, but the
    // presence of reversed tracks complicates that slightly.
    for(size_t n = 0; n < fTraj.size()-1; ++n){
      if( zDir*fTraj[n]  .Z() <= zDir*z && 
         zDir*fTraj[n+1].Z() >= zDir*z ){
        // If we've already found a pair of neighbours then this new pair is
        // ambiguous which is bad. Probably it means the track reverses
        // direction wrt Z.
        if(found){
          // But, it might also mean the value passed in was right on a
          // trajectory point. In which case there isn't really a problem.
          // Either choice will end up just giving the value at this point.
          if(std::abs(fTraj[n].Z()-z) > 1e-3 && std::abs(fTraj[n+1].Z()-z) > 1e-3){
            LOG_DEBUG("Track")
              << "Choice of neighbouring trajectory points is ambiguous";
          }
        }
        i0 = n;
        i1 = n+1;
        found = true;
        break;
      }
    }

    if(!found)
      throw cet::exception("Track") << "unable to find neighboring point index";
  }

void rb::Track::FindNeighbouringPoints ( double  z, TVector3 &  pt0, TVector3 &  pt1  ) const

protected

Helper function for W.

Definition at line 432 of file Track.cxx.

References FindNeighbouringPointIndices(), and fTraj.

Referenced by InterpolateDir(), InterpolateXY(), and Trajectory().

  {
    int i0, i1;
    FindNeighbouringPointIndices(z, i0, i1);
    pt0 = fTraj[i0];
    pt1 = fTraj[i1];
  }

const int rb::Cluster::ID ( ) const

inlineinherited

Definition at line 75 of file Cluster.h.

References rb::Cluster::fID, and rb::Cluster::W().

Referenced by align::SplitTracks::analyze(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), jmshower::RecoJMShower::ContStartPlane(), jmshower::RecoJMShower::DepositEnergy(), jmshower::RecoJMShower::Energy(), sn::SNSlicer::FinalizeCluster(), jmshower::RecoJMShower::GetCellTransDedx(), jmshower::RecoJMShower::GetDedxInvLongLL(), jmshower::RecoJMShower::GetDedxLongLL(), jmshower::RecoJMShower::GetPlaneDedx(), jmshower::RecoJMShower::GetTrkDistToEdge(), jmshower::RecoJMShower::GetTrkHitPath(), jmshower::RecoJMShower::IsFiducial(), trk::KalmanTrackMerge::JoinTracks(), jmshower::RecoJMShower::NMIPPlane(), rb::operator<<(), jmshower::RecoJMShower::Radius(), jmshower::RecoJMShower::RecoShowers(), slid::ParticleIDAlg::SetShower(), and trk::KalmanTrackMerge::ViewMergeTracks().

{ return fID; }

TVector3 rb::Track::InterpolateDir ( double  z ) const

virtual

Parameters

z	Coordinate at which trajectory is to be evaluated

Returns

The track direction estimated at that point

Definition at line 348 of file Track.cxx.

References a, b, and FindNeighbouringPoints().

Referenced by upmuana::UpMuRecoAna::analyze(), upmuana::UpMuAnalysis::analyze(), PlaneDirMap(), and Trajectory().

  {
    // The algorithm here is basically valid for z-directed tracks

    TVector3 a, b;
    FindNeighbouringPoints(z, a, b);

    return (b-a).Unit();
  }

void rb::Track::InterpolatePts ( double  z, const TVector3 &  a, const TVector3 &  b, double &  x, double &  y  ) const

protected

Helper function for InterpolateXY and ZipWith.

Definition at line 336 of file Track.cxx.

References makeHTMLView::dt.

Referenced by InterpolateXY(), and Trajectory().

  {
    const double denom = b.Z()-a.Z();
    const double dt = denom ? (z-a.Z())/denom : 0.5;

    x = (1-dt)*a.X()+dt*b.X();
    y = (1-dt)*a.Y()+dt*b.Y();
  }

void rb::Track::InterpolateXY ( double  z, double &  x, double &  y  ) const

virtual

Parameters

	z	Coordinate at which trajectory is to be evaluated
[out]	x	x-position of trajectory at this z
[out]	y	y-position of trajectory at this z

Definition at line 325 of file Track.cxx.

References a, b, FindNeighbouringPoints(), and InterpolatePts().

Referenced by align::SplitTracks::analyze(), murem::TrackCleanUpAlg::DeDxInPlane(), dt::ViewMerger::FakeThirdDimension(), cosrej::CosRej::produce(), murem::TrackCleanUpAlg::SortByDistFromTrack(), dt::ViewMerger::TotalChargePerView(), Trajectory(), trk::KalmanTrackMerge::ViewMergeTracks(), W(), and ZipWith().

  {
    // The algorithm here is basically valid for z-directed tracks

    TVector3 a, b;
    FindNeighbouringPoints(z, a, b);

    InterpolatePts(z, a, b, x, y);
  }

bool rb::Cluster::Is2D ( ) const

inlineinherited

Definition at line 96 of file Cluster.h.

References rb::Cluster::Is3D().

Referenced by dt::ViewMerger::Adjacent(), moonshadowana::MoonShadowAna::analyze(), slid::DeconvolveAlg::DistanceToCore(), evd::RecoBaseDrawer::DrawCluster2D(), bpfit::BreakPoint::GoodProng(), slid::Recluster::produce(), moonshadowana::MoonShadowAna::WriteEventInfo(), moonshadowana::MoonShadowAnaHough::WriteEventInfo(), and ZipWith().

{return !Is3D();}

virtual bool rb::Prong::Is3D ( ) const

inlinevirtualinherited

Reimplemented from rb::Cluster.

Definition at line 71 of file Prong.h.

References rb::Cluster::fView, and geo::kXorY.

Referenced by align::SplitTracks::analyze(), comi::NumiFilteringAna::analyze(), trk::KalmanTrackAna::analyze(), calib::StopperThreshold::analyze(), air::AirKalmanAna::analyze(), upmuana::UpMuRecoAna::analyze(), upmuana::UpMuAnalysis::analyze(), dt::ViewMerger::FakeThirdDimension(), comi::NumiFiltering::filter(), qeef::QeFinder::produce(), wsnumu::WSNumu::produce(), numusand::FillSandbox::produce(), and cosrej::CosRej::produce().

{return fView == geo::kXorY;}

bool rb::Cluster::IsNoise ( ) const

inlineinherited

Is the noise flag set?

Definition at line 163 of file Cluster.h.

References rb::Cluster::fNoiseCluster, rb::Cluster::SetWeight(), rb::Cluster::TotalADC(), rb::Cluster::TotalPE(), and ana::weight.

Referenced by qeef::QeFinderVal::analyze(), rsim::RecordNoiseSpectrumFile::analyze(), tut::RecoValidationTutorial::analyze(), calib::TestParticleCorrections::analyze(), mcchk::MonopoleAna::analyze(), sn::SupernovaAna::analyze(), comi::DataCheck::analyze(), validation::NoiseClusterAnalysis::analyze(), comi::NumiFilteringAna::analyze(), showere::ShowerEnergyAna::analyze(), hough::HoughValidate::analyze(), rockmuon::NDRockMuon::analyze(), ncid::NCNNKerasVal::analyze(), bpf::BPFCVNAna::analyze(), ncs::Xbeam::analyze(), ncs::Xeff::analyze(), ncs::ROCKMRE::analyze(), slid::LIDTraining::analyze(), ncs::NCAna::analyze(), me::MEFinder::BestSlcMatch(), trk::KalmanTrackAna::CheckAssociations(), air::AirKalmanAna::CheckAssociations(), evd::RecoBaseDrawer::DrawCluster2D(), LSTME::FillLSTME::fillEnergies(), SliceLID::FillSliceLID::fillSliceLIDPreds(), filter::Filter::FillSliceVariables(), comi::NumiFiltering::filter(), bsf::BremShowerFilter::filter(), showere::ShowerEnergyFilterMC::filter(), trident::DimuonSliceAna::filter(), me::MEFinder::GetPhysNoiseSlcs(), mono::HighEnergySliceRemoval::is_high_energy_slice(), novaddt::HorizontalMuonAna::IsHorzMuonTrack(), skim::EvaluatorNue::KeepSlice(), airshower::AirSlicer::MakeTrackSlices(), murem::FindParent::MatchToOrigSlice(), air::AirFilter::parallel_muon_slice(), lem::LEMAssociator::produce(), presel::LEMPresel::produce(), qeef::QeFinder::produce(), hough::HoughT::produce(), hv::HoughVertex::produce(), nuesand::FillNueSandbox::produce(), lem::UnTranspose::produce(), wsnumu::WSNumu::produce(), numue::TrackOverlapECalc::produce(), lem::LEMSummarizer::produce(), presel::RockPresel::produce(), slid::Recluster::produce(), nussand::FillNusSandbox::produce(), presel::CosmicVeto::produce(), trackinfo::TrackInfo::produce(), lem::Preselection::produce(), presel::NuePresel::produce(), jmshower::JMTrackMerge::produce(), presel::CosmicCVNVeto::produce(), hough::MultiHoughT::produce(), numusand::FillSandbox::produce(), cosrej::CosRej::produce(), numue::NumuEnergy::produce(), earms::ElasticArmsHS::produce(), dt::DiscreteTracker::produce(), ndreco::NDReco::produce(), jmshower::JMClusterMerge::produce(), vdt::VertexDT::produce(), ncpi0::MakeNCPi0BkgRej::produce(), jmshower::NueSel::produce(), rvp::RecVarPID::produce(), crvtx::CosmicRayVertex::produce(), dif::DiFShowerFinder::produce(), xnue::XnuePID::produce(), jmshower::NueSelLID::produce(), ncid::MakeNCCosRej::produce(), cosrej::MakeNueCosRej::produce(), ncid::NCNNKeras::produce(), murem::FindParent::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), caf::CAFMaker::produce(), bpfit::DimuonFitter::produce(), jmshower::RecoJMShower::RecoShowers(), and sn::SliceRemover::remove().

{return fNoiseCluster;}

unsigned int rb::Cluster::MaxCell ( geo::View_t  view ) const

inherited

Definition at line 518 of file Cluster.cxx.

References ana::assert(), rb::CellHit::Cell(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), and runNovaSAM::ret.

Referenced by align::SplitTracks::analyze(), cvn::PixelMapProducer::DefineBoundary(), rb::Cluster::ExtentCell(), caf::FillHadClustVars(), caf::FillSliceVars(), hough::HoughT::is_good_slice(), rb::Cluster::MaxZ(), lem::MeanCell(), dt::View::View(), trk::KalmanTrackMerge::ViewMergeTracks(), and dt::ViewMerger::ZipVerticalTracks().

  {
    assert(NCell(view) > 0);
    unsigned int ret = 0;
    for(unsigned int i = 0; i < NCell(view); ++i)
      if(Cell(view, i)->Cell() > ret) ret = Cell(view, i)->Cell();
    return ret;
  }

unsigned int rb::Cluster::MaxPlane ( geo::View_t  view = geo::kXorY ) const

inherited

Definition at line 508 of file Cluster.cxx.

References ana::assert(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), rb::CellHit::Plane(), and runNovaSAM::ret.

Referenced by dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), align::SplitTracks::analyze(), calib::HitEfficiency::analyze(), tbem::TBRecoAna::analyze(), trk::KalmanTrackMerge::CalcMatchScore(), slid::ParticleIDAlg::CalcPlaneHits(), trk::KalmanTrackMerge::CanJoinTracks(), skim::Evaluator::CheckSliceQuality(), CerenkovSelection::Contained(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), murem::TrackCleanUpAlg::DeDxInPlane(), lem::DefaultVertex(), cvn::PixelMapProducer::DefineBoundary(), dif::DiFShowerFinder::eparm(), rb::Cluster::ExtentPlane(), caf::FillDiFVars(), caf::FillHadClustVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillSliceInfo(), caf::FillSliceVars(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), bsf::BremShowerFilter::findShowerByTruth(), trk::KalmanTrack::FindTracks(), cvn::RegPixelMapProducer::FindVertexWindowThreshold(), cvn::PixelMapProducer::FindVertexWindowThreshold(), hough::HoughT::is_good_slice(), trk::CosmicTrackSelection::IsCalibTrack(), calib::PCHitsList::IsGoodTrack(), trk::KalmanTrackMerge::JoinTracks(), dt::ViewMerger::JoinTracks(), dt::ViewMerger::MatchScore(), murem::FindParent::MatchToOrigSlice(), rb::Cluster::MaxZ(), dif::DiFShowerFinder::muonstub(), skim::ParametersNumu::ParametersNumu(), slid::Recluster::produce(), remid::RecoMuon::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), dt::ViewMerger::ScoreForExtension(), dt::ViewMerger::ScoreForJoinPlusExtension(), trk::WindowTrackingAlg::SetTrackEndPoints(), cosrej::MakeNueCosRej::SparsenessAsymmetry(), slid::NuEEnergyAlg::VertexEnergy(), dt::View::View(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    assert(NCell(view) > 0);
    unsigned int ret = 0;
    for(unsigned int i = 0; i < NCell(view); ++i)
      if(Cell(view, i)->Plane() > ret) ret = Cell(view, i)->Plane();
    return ret;
  }

double rb::Cluster::MaxTNS ( ) const

inherited

Definition at line 528 of file Cluster.cxx.

References ana::assert(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), runNovaSAM::ret, and rb::CellHit::TNS().

Referenced by align::AlignInspect::analyze(), nnbar::ShowerMuonCoincidence::analyze(), ncs::NCAna::analyze(), evd::RecoBaseDrawer::DrawCluster2D(), evd::RecoBaseDrawer::DrawTrack2D(), rb::Cluster::ExtentTNS(), caf::FillProngVars(), caf::FillSliceVars(), bsf::BremShowerFilter::findShowerByTruth(), rb::Cluster::MaxZ(), skim::ParametersNumu::ParametersNumu(), trk::TrkAssn::produce(), jmshower::JMTrackMerge::produce(), jmshower::JMClusterMerge::produce(), caf::CAFMaker::produce(), and sn::TrackActivityRemover::remove().

  {
    assert(NCell() > 0);
    double ret = -DBL_MAX;
    for(unsigned int i = 0; i < NCell(); ++i)
      if(Cell(i)->TNS() > ret) ret = Cell(i)->TNS();
    return ret;
  }

double rb::Cluster::MaxV ( geo::View_t  view ) const

inherited

Definition at line 500 of file Cluster.cxx.

References ana::assert(), geo::kX, geo::kY, and rb::Cluster::MaxXYZ().

Referenced by rb::Cluster::ExtentV(), rb::Cluster::MinZ(), fuzz::FuzzyKVertex::produce(), and nerd::NERDProng::produce().

  {
    assert(view == geo::kX || view == geo::kY);
    if(view == geo::kX) return MaxXYZ().X();
    return MaxXYZ().Y();
  }

double rb::Cluster::MaxX ( ) const

inlineinherited

Definition at line 217 of file Cluster.h.

References rb::Cluster::MaxXYZ().

Referenced by caf::FillProngVars(), and presel::CosmicVeto::produce().

{return MaxXYZ().X();}

TVector3 rb::Cluster::MaxXYZ ( ) const

inherited

Definition at line 492 of file Cluster.cxx.

References hi(), and rb::Cluster::MinMaxMeanXYZ().

Referenced by comi::NumiFilteringAna::analyze(), rb::Cluster::ExtentXYZ(), comi::NumiFiltering::filter(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), skim::EvaluatorNue::KeepSlice(), rb::Cluster::MaxV(), rb::Cluster::MaxX(), rb::Cluster::MaxY(), rb::Cluster::MaxZ(), rb::Cluster::MinZ(), skim::ParametersNumu::ParametersNumu(), nuesand::FillNueSandbox::produce(), presel::CosmicVeto::produce(), cosrej::CosRej::produce(), and caf::CAFMaker::produce().

  {
    TVector3 hi, junk;
    MinMaxMeanXYZ(junk, hi, junk);
    return hi;
  }

double rb::Cluster::MaxY ( ) const

inlineinherited

Definition at line 218 of file Cluster.h.

References rb::Cluster::MaxXYZ().

Referenced by caf::FillProngVars(), and presel::CosmicVeto::produce().

{return MaxXYZ().Y();}

double rb::Cluster::MaxZ ( ) const

inlineinherited

Definition at line 219 of file Cluster.h.

References geo::kXorY, rb::Cluster::MaxCell(), rb::Cluster::MaxPlane(), rb::Cluster::MaxTNS(), rb::Cluster::MaxXYZ(), rb::Cluster::MeanV(), and rb::Cluster::MeanXYZ().

Referenced by fuzz::FuzzyKVertex::produce(), and nerd::NERDProng::produce().

{return MaxXYZ().Z();}

double rb::Cluster::MeanTNS ( rb::AveragingScheme  scheme = kDefaultScheme ) const

inherited

Definition at line 554 of file Cluster.cxx.

References ana::assert(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::kByEnergy, rb::Cluster::NCell(), rb::CellHit::PE(), runNovaSAM::ret, rb::CellHit::TNS(), w, and rb::Cluster::Weight().

Referenced by me::MEFinder::AddMEHit(), sn::SNSlicerAna::analyze(), tbem::TBRecoAna::analyze(), ncs::NCAna::analyze(), sn::SNMichelAnalyzer::analyze(), me::MEFinder::BestSlcMatch(), skim::Evaluator::CheckTiming(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), CerenkovSelection::Contained(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), caf::FillProngVars(), caf::FillSliceVars(), comi::NumiFiltering::filter(), bsf::BremShowerFilter::filter(), tf::TimingFitAlg::HoughFit(), skim::EvaluatorNue::KeepSlice(), me::MEFinder::MakeSlcME(), me::MEFinder::MakeTrkME(), rb::Cluster::MeanZ(), cosrej::MakeNueCosRej::MuonParentByDist(), cosrej::MakeNueCosRej::MuonParentByTime(), calib::PCHitsList::ProcessTrackForFLSHits(), recodemo::RecoDemo::produce(), hv::HoughVertex::produce(), presel::CosmicVeto::produce(), presel::CosmicCVNVeto::produce(), sn::SupernovaMCCluster::produce(), vdt::VertexDT::produce(), jmshower::NueSel::produce(), dif::DiFShowerFinder::produce(), jmshower::NueSelLID::produce(), and caf::CAFMaker::produce().

  {
    assert(NCell() > 0);
    double ret = 0;
    double totW = 0;
    for(unsigned int i = 0; i < NCell(); ++i){
      double w = Weight(i);
      if(scheme == kByEnergy) w *= Cell(i)->PE();

      ret += Cell(i)->TNS()*w;
      totW += w;
    }
    return ret/totW;
  }

double rb::Cluster::MeanV ( geo::View_t  view, rb::AveragingScheme  scheme = kDefaultScheme  ) const

inherited

Definition at line 546 of file Cluster.cxx.

References ana::assert(), geo::kX, geo::kY, and rb::Cluster::MeanXYZ().

Referenced by nerd::NERDProng::AddProngs(), rb::Cluster::MaxZ(), and rb::Cluster::W().

  {
    assert(view == geo::kX || view == geo::kY);
    if(view == geo::kX) return MeanXYZ(scheme).X();
    return MeanXYZ(scheme).Y();
  }

double rb::Cluster::MeanX ( rb::AveragingScheme  scheme = kDefaultScheme ) const

inlineinherited

Definition at line 231 of file Cluster.h.

References rb::Cluster::MeanXYZ().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), calib::TestParticleCorrections::analyze(), sn::SNSlicerAna::analyze(), me::TrkME::DistToTrk(), slid::NuEEnergyAlg::HadronicDepEnergy(), recodemo::RecoDemo::produce(), sn::SupernovaMCCluster::produce(), and me::MEFinder::TrkMEDist().

{return MeanXYZ(scheme).X();}

TVector3 rb::Cluster::MeanXYZ ( rb::AveragingScheme  scheme = kDefaultScheme ) const

inherited

Definition at line 538 of file Cluster.cxx.

References extractScale::mean, and rb::Cluster::MinMaxMeanXYZ().

Referenced by sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), jmshower::RecoJMShower::GetRVPStats(), tf::TimingFitAlg::HoughFit(), rb::Cluster::MaxZ(), rb::Cluster::MeanV(), rb::Cluster::MeanX(), rb::Cluster::MeanY(), rb::Cluster::MeanZ(), nuesand::FillNueSandbox::produce(), presel::NuePresel::produce(), lem::Preselection::produce(), and rb::Cluster::TotalGeVFastClusterOnly().

  {
    TVector3 mean, junk;
    MinMaxMeanXYZ(junk, junk, mean, scheme);
    return mean;
  }

double rb::Cluster::MeanY ( rb::AveragingScheme  scheme = kDefaultScheme ) const

inlineinherited

Definition at line 232 of file Cluster.h.

References rb::Cluster::MeanXYZ().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), sn::SNSlicerAna::analyze(), me::TrkME::DistToTrk(), slid::NuEEnergyAlg::HadronicDepEnergy(), recodemo::RecoDemo::produce(), sn::SupernovaMCCluster::produce(), and me::MEFinder::TrkMEDist().

{return MeanXYZ(scheme).Y();}

double rb::Cluster::MeanZ ( rb::AveragingScheme  scheme = kDefaultScheme ) const

inlineinherited

Definition at line 233 of file Cluster.h.

References hi(), lo(), extractScale::mean, rb::Cluster::MeanTNS(), rb::Cluster::MeanXYZ(), and rb::Cluster::MinMaxMeanXYZ().

Referenced by nerd::NERDProng::AddProngs(), caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), sn::SNSlicerAna::analyze(), trk::KalmanTrackMerge::CanJoinTracks(), me::TrkME::DistToTrk(), rb::Cluster::operator<(), recodemo::RecoDemo::produce(), sn::SupernovaMCCluster::produce(), and me::MEFinder::TrkMEDist().

{return MeanXYZ(scheme).Z();}

unsigned int rb::Cluster::MinCell ( geo::View_t  view ) const

inherited

Definition at line 472 of file Cluster.cxx.

References ana::assert(), rb::CellHit::Cell(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), and runNovaSAM::ret.

Referenced by calib::TestParticleCorrections::analyze(), align::SplitTracks::analyze(), cvn::PixelMapProducer::DefineBoundary(), rb::Cluster::ExtentCell(), caf::FillHadClustVars(), caf::FillSliceVars(), hough::HoughT::is_good_slice(), lem::MeanCell(), rb::Cluster::MinZ(), dt::View::View(), trk::KalmanTrackMerge::ViewMergeTracks(), and dt::ViewMerger::ZipVerticalTracks().

  {
    assert(NCell(view) > 0);
    unsigned int ret = UINT_MAX;
    for(unsigned int i = 0; i < NCell(view); ++i)
      if(Cell(view, i)->Cell() < ret) ret = Cell(view, i)->Cell();
    return ret;
  }

void rb::Cluster::MinMaxMeanXYZ ( TVector3 &  lo, TVector3 &  hi, TVector3 &  mean, rb::AveragingScheme  scheme = kDefaultScheme  ) const

inherited

Gets the min/max/mean all at once, called by the functions above.

Todo:

Won't be quite right in the FD when cells aren't at right angles?

Definition at line 582 of file Cluster.cxx.

References ana::assert(), rb::CellHit::Cell(), rb::Cluster::Cell(), geom(), make_syst_table_plots::h, hi(), MECModelEnuComparisons::i, rb::kByEnergy, geo::kX, geo::kY, lo(), extractScale::mean, getGoodRuns4SAM::n, rb::Cluster::NCell(), rb::CellHit::PE(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), POTSpillRate::view, rb::Cluster::View(), ana::weight, rb::Cluster::Weight(), submit_syst::x, submit_syst::y, and test::z.

Referenced by caf::FillHadClustVars(), caf::FillSliceVars(), rb::Cluster::MaxXYZ(), rb::Cluster::MeanXYZ(), rb::Cluster::MeanZ(), and rb::Cluster::MinXYZ().

  {
    // If this becomes a bottleneck we can cache the results

    /// \todo Won't be quite right in the FD when cells aren't at right angles?

    assert(NCell() > 0);

    art::ServiceHandle<geo::Geometry> geom;

    double Wx = 0;
    double Wy = 0;
    for(int n = 0; n < 3; ++n){
      lo(n) = DBL_MAX;
      hi(n) = -DBL_MAX;
      mean(n) = 0;
    }

    for(unsigned int i = 0; i < NCell(); ++i){
      const art::Ptr<rb::CellHit>& h = Cell(i);
      double xyz[3];
      geom->Plane(h->Plane())->Cell(h->Cell())->GetCenter(xyz);
      const double x = xyz[0];
      const double y = xyz[1];
      const double z = xyz[2];
      const geo::View_t view = geom->Plane(h->Plane())->View();
      double weight = Weight(i);
      if(scheme == kByEnergy) weight *= h->PE();

      if(view == geo::kX){
        if(x < lo.X()) lo.SetX(x);
        if(x > hi.X()) hi.SetX(x);
        mean.SetX(mean.X()+x*weight);
        Wx += weight;
      }
      if(view == geo::kY){
        if(y < lo.Y()) lo.SetY(y);
        if(y > hi.Y()) hi.SetY(y);
        mean.SetY(mean.Y()+y*weight);
        Wy += weight;
      }
      if(z < lo.Z()) lo.SetZ(z);
      if(z > hi.Z()) hi.SetZ(z);
      mean.SetZ(mean.Z()+z*weight);
    }

    if(Wx) mean.SetX(mean.X()/Wx);
    if(Wy) mean.SetY(mean.Y()/Wy);
    if(Wx || Wy) mean.SetZ(mean.Z()/(Wx+Wy));
  }

unsigned int rb::Cluster::MinPlane ( geo::View_t  view = geo::kXorY ) const

inherited

Definition at line 462 of file Cluster.cxx.

References ana::assert(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), rb::CellHit::Plane(), and runNovaSAM::ret.

Referenced by dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), align::AlignInspect::analyze(), calib::TestParticleCorrections::analyze(), align::SplitTracks::analyze(), calib::HitEfficiency::analyze(), ncs::ROCKMRE::analyze(), ncs::NCAna::analyze(), trk::KalmanTrackMerge::CalcMatchScore(), slid::ParticleIDAlg::CalcPlaneHits(), trk::KalmanTrackMerge::CanJoinTracks(), skim::Evaluator::CheckSliceQuality(), CerenkovSelection::Contained(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), murem::TrackCleanUpAlg::DeDxInPlane(), lem::DefaultVertex(), cvn::PixelMapProducer::DefineBoundary(), dif::DiFShowerFinder::eparm(), rb::Cluster::ExtentPlane(), caf::FillDiFVars(), caf::FillHadClustVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillSliceInfo(), caf::FillSliceVars(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), bsf::BremShowerFilter::findShowerByTruth(), trk::KalmanTrack::FindTracks(), cvn::RegPixelMapProducer::FindVertexMaxGap(), cvn::PixelMapProducer::FindVertexMaxGap(), cvn::RegPixelMapProducer::FindVertexWindowThreshold(), cvn::PixelMapProducer::FindVertexWindowThreshold(), jmshower::RecoJMShower::GetRVPStats(), hough::HoughT::is_good_slice(), trk::CosmicTrackSelection::IsCalibTrack(), calib::PCHitsList::IsGoodTrack(), trk::KalmanTrackMerge::JoinTracks(), dt::ViewMerger::JoinTracks(), dt::ViewMerger::MatchScore(), murem::FindParent::MatchToOrigSlice(), rb::Cluster::MinZ(), skim::ParametersNue::ParametersNue(), skim::ParametersNumu::ParametersNumu(), slid::Recluster::produce(), rvp::RecVarPID::produce(), remid::RecoMuon::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), jmshower::RecoJMShower::RecoShowers(), dt::ViewMerger::ScoreForExtension(), dt::ViewMerger::ScoreForJoinPlusExtension(), trk::WindowTrackingAlg::SetTrackEndPoints(), cosrej::MakeNueCosRej::SparsenessAsymmetry(), slid::NuEEnergyAlg::VertexEnergy(), dt::View::View(), trk::KalmanTrackMerge::ViewMergeTracks(), and dt::ViewMerger::ZipVerticalTracks().

  {
    assert(NCell(view) > 0);
    unsigned int ret = UINT_MAX;
    for(unsigned int i = 0; i < NCell(view); ++i)
      if(Cell(view, i)->Plane() < ret) ret = Cell(view, i)->Plane();
    return ret;
  }

double rb::Cluster::MinTNS ( ) const

inherited

Definition at line 482 of file Cluster.cxx.

References ana::assert(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), runNovaSAM::ret, and rb::CellHit::TNS().

Referenced by align::AlignInspect::analyze(), sn::SNSlicerAna::analyze(), tbem::TBRecoAna::analyze(), nnbar::ShowerMuonCoincidence::analyze(), ncs::NCAna::analyze(), rb::CompareTracksByTime(), evd::RecoBaseDrawer::DrawCluster2D(), evd::RecoBaseDrawer::DrawTrack2D(), rb::Cluster::ExtentTNS(), caf::FillProngVars(), caf::FillSliceVars(), bsf::BremShowerFilter::findShowerByTruth(), rb::Cluster::MinZ(), murem::MuonRemove::MuonInfo(), skim::ParametersNumu::ParametersNumu(), trk::TrkAssn::produce(), jmshower::JMTrackMerge::produce(), earms::ElasticArmsHS::produce(), jmshower::JMClusterMerge::produce(), caf::CAFMaker::produce(), sn::TrackActivityRemover::remove(), sn::SortClustersByTime(), and sn::SortTracksByTime().

  {
    assert(NCell() > 0);
    double ret = DBL_MAX;
    for(unsigned int i = 0; i < NCell(); ++i)
      if(Cell(i)->TNS() < ret) ret = Cell(i)->TNS();
    return ret;
  }

double rb::Cluster::MinV ( geo::View_t  view ) const

inherited

Definition at line 454 of file Cluster.cxx.

References ana::assert(), geo::kX, geo::kY, and rb::Cluster::MinXYZ().

Referenced by rb::Cluster::ExtentV(), fuzz::FuzzyKVertex::produce(), and nerd::NERDProng::produce().

  {
    assert(view == geo::kX || view == geo::kY);
    if(view == geo::kX) return MinXYZ().X();
    return MinXYZ().Y();
  }

double rb::Cluster::MinX ( ) const

inlineinherited

Definition at line 204 of file Cluster.h.

References rb::Cluster::MinXYZ().

Referenced by caf::FillProngVars(), and presel::CosmicVeto::produce().

{return MinXYZ().X();}

TVector3 rb::Cluster::MinXYZ ( ) const

inherited

Definition at line 446 of file Cluster.cxx.

References lo(), and rb::Cluster::MinMaxMeanXYZ().

Referenced by comi::NumiFilteringAna::analyze(), rb::Cluster::ExtentXYZ(), comi::NumiFiltering::filter(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), rb::Cluster::MinV(), rb::Cluster::MinX(), rb::Cluster::MinY(), rb::Cluster::MinZ(), skim::ParametersNumu::ParametersNumu(), nuesand::FillNueSandbox::produce(), presel::CosmicVeto::produce(), cosrej::CosRej::produce(), and caf::CAFMaker::produce().

  {
    TVector3 lo, junk;
    MinMaxMeanXYZ(lo, junk, junk);
    return lo;
  }

double rb::Cluster::MinY ( ) const

inlineinherited

Definition at line 205 of file Cluster.h.

References rb::Cluster::MinXYZ().

Referenced by caf::FillProngVars(), and presel::CosmicVeto::produce().

{return MinXYZ().Y();}

double rb::Cluster::MinZ ( ) const

inlineinherited

Definition at line 206 of file Cluster.h.

References geo::kXorY, rb::Cluster::MaxV(), rb::Cluster::MaxXYZ(), rb::Cluster::MinCell(), rb::Cluster::MinPlane(), rb::Cluster::MinTNS(), and rb::Cluster::MinXYZ().

Referenced by fuzz::FuzzyKVertex::produce(), and nerd::NERDProng::produce().

{return MinXYZ().Z();}

int rb::Cluster::MostContiguousPlanes ( geo::View_t  view ) const

inherited

Longest run of adjacent planes with hits.

view may be kXorY to require hits on every plane no matter the view

Definition at line 635 of file Cluster.cxx.

References rb::Cluster::Cell(), MECModelEnuComparisons::i, geo::kXorY, std::max(), rb::Cluster::NCell(), and NDAPDHVSetting::plane.

Referenced by ncid::NCNNKerasVal::analyze(), rb::Cluster::ExtentTNS(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillHadClustVars(), caf::FillProngVars(), caf::FillShowerVars(), caf::FillSliceInfo(), caf::FillSliceVars(), and skim::ParametersNumu::ParametersNumu().

  {
    std::set<int> planes;

    for(unsigned int i = 0; i < NCell(view); ++i)
      planes.insert(Cell(view, i)->Plane());

    if(planes.empty()) return 0;

    // In 3D expect a hit every plane, in 2D only every other plane.
    //
    // This isn't right if the planes aren't alternating in views, but
    // thankfully we stopped building detectors like that.
    const int expGap = (view == geo::kXorY) ? 1 : 2;

    int prevPlane = -999;
    int curRun = 1;
    int longestRun = 0;
    for(auto plane: planes){
      if(plane - prevPlane == expGap){
        ++curRun;
        longestRun = std::max(longestRun, curRun);
      }
      else{
        curRun = 1;
      }
      prevPlane = plane;
    } // end for plane

    return longestRun;
  }

int rb::Cluster::MostMissingPlanes ( geo::View_t  view ) const

inherited

Longest run of adjacent planes with no hits.

view may be kXorY to allow a hit in either view to end the run

Definition at line 668 of file Cluster.cxx.

References rb::Cluster::Cell(), MECModelEnuComparisons::i, geo::kXorY, std::max(), rb::Cluster::NCell(), and NDAPDHVSetting::plane.

Referenced by rb::Cluster::ExtentTNS(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillProngVars(), caf::FillShowerVars(), and xsrec::MakeXSecCCPi0Inc::GetPhLL().

  {
    std::set<int> planes;

    for(unsigned int i = 0; i < NCell(view); ++i)
      planes.insert(Cell(view, i)->Plane());

    if(planes.empty()) return 0;

    int maxDiff = 1;
    int prevPlane = *planes.begin();
    for(int plane: planes){
      maxDiff = std::max(maxDiff, plane-prevPlane);
      prevPlane = plane;
    }

    // One plane difference as expected in 3D isn't a gap
    if(view == geo::kXorY) return maxDiff-1;

    // Two plane difference in 2D isn't a gap. Beyond that, every two planes
    // difference only indicates one missing plane of hits.
    return (maxDiff-2)/2;
  }

unsigned int rb::Cluster::NCell ( geo::View_t  view ) const

inherited

Number of cells in view view.

Definition at line 134 of file Cluster.cxx.

References ana::assert(), geo::kX, geo::kXorY, geo::kY, rb::Cluster::NCell(), rb::Cluster::NXCell(), and rb::Cluster::NYCell().

Referenced by dqsf::DQSpillFlags::AccumulateSyncMetric(), rb::HitMap::Add(), fuzz::FuzzyKVertex::AddProng(), dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), tut::TutAnalyzer::analyze(), qeef::QeFinderVal::analyze(), rsim::RecordNoiseSpectrumFile::analyze(), chaninfo::CosmicEff::analyze(), remid::ReMIdValidate::analyze(), comi::CosmicMetrics::analyze(), align::AlignInspect::analyze(), tut::RecoValidationTutorial::analyze(), dprf::ChannelPlots::analyze(), calib::TestParticleCorrections::analyze(), remid::ReMIdTrain::analyze(), dprf::TrackPlots::analyze(), mcchk::MonopoleAna::analyze(), align::Alignment::analyze(), mcchk::RockAna::analyze(), sn::SNSlicerAna::analyze(), sn::SupernovaAna::analyze(), murem::MuonRemoveAna::analyze(), align::SplitTracks::analyze(), pa::ClusterEva::analyze(), comi::DataCheck::analyze(), vf::TrackEva::analyze(), comi::NumiFilteringAna::analyze(), zcl::FmmTriggerAna::analyze(), tbem::TBRecoAna::analyze(), zcl::FmmTrackerValidation::analyze(), zcl::FmmTrackerAna::analyze(), trk::KalmanTrackAna::analyze(), showere::ShowerEnergyAna::analyze(), hough::HoughValidate::analyze(), rockmuon::NDRockMuon::analyze(), zcl::FastMMStudy::analyze(), bpfit::BPFTmvaTrainer::analyze(), bpf::BPFCVNAna::analyze(), calib::CosmicTrends::analyze(), ncs::Xeff::analyze(), calib::CalibAna::analyze(), ncs::ROCKMRE::analyze(), moonshadowana::MoonShadowAna::analyze(), ncs::NCAna::analyze(), htk::HoughTrack::analyze(), bpfit::DimuonFitter::AnaSlice(), trident::DimuonSliceAna::AnaSlice(), me::MEFinder::BestSlcMatch(), fuzz::ViewMatchAlg::CalcEnergyProfile(), nerd::ViewMatchAlg::CalcEnergyProfile(), trk::KalmanTrackMerge::CalcMatchScore(), slid::ParticleIDAlg::CalcPlaneHits(), trk::KalmanTrackMerge::CanJoinTracks(), trk::KalmanTrackAna::CheckRecoTracks(), murem::TrackCleanUpAlg::CleanUpTrack(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), sn::SNSlicer::ClusterIsGood(), murem::TrackCleanUpAlg::ComputeVertexRegion(), CerenkovSelection::Contained(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxFD::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), nnbar::NNbarUtilities::correlationTCell(), nnbar::NNbarUtilities::correlationTPlane(), cvn::RegPixelMapProducer::CreateMapGivenBoundary(), cvn::PixelMapProducer::CreateMapGivenBoundary(), cvn::RegPixelMapProducer::CreateMapGivenShowerVertex(), cvn::RegPixelMapProducer::CreateMapGivenVertex(), jmshower::RecoJMShower::DepositEnergy(), evd::RecoBaseDrawer::DrawTrack2D(), skim::ParametersNumu::EarliestLatestHitPos(), jmshower::RecoJMShower::Energy(), dif::DiFShowerFinder::eparm(), rb::Cluster::Exclude(), dt::ViewMerger::Extremes(), numue::NumuEAlg::FDEnergy(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillHadClustVars(), calib::StopperThreshold::FillHist(), bpfit::BreakPoint::FillHitList(), sn::SupernovaAna::FillPlaneSeperations(), dqsf::DQSpillFlags::FillPopulatedDCMs(), caf::FillProngVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillShowerVars(), caf::FillSliceInfo(), filter::Filter::FillSliceVariables(), caf::FillSliceVars(), trk::CosmicTrackAna::FillTrackHistograms(), CerenkovSelection::FillTree(), calib::StopperThreshold::FillTree(), comi::NumiFiltering::filter(), bsf::BremShowerFilter::filter(), showere::ShowerEnergyFilterMC::filter(), mono::MonopoleTrack::filter(), cvn::RegPixelMapProducer::FindCenterMedian(), cvn::PixelMapProducer::FindCenterMedian(), airshower::AirSlicer::FindNeighbors(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), trk::KalmanTrack::FindTracks(), trk::CosmicTrackUtilities::FindTriCells(), trident::DimuonSliceAna::FindVertex(), cvn::RegPixelMapProducer::FindVertexMaxGap(), cvn::PixelMapProducer::FindVertexMaxGap(), cvn::RegPixelMapProducer::FindVertexWindowThreshold(), cvn::PixelMapProducer::FindVertexWindowThreshold(), bpfit::DimuonFitter::FindVertexZ(), bpfit::BreakPoint::FitTracks(), bpfit::DimuonFitter::FitView(), nuesand::FillNueSandbox::FracAngChanges(), nuesand::FillNueSandbox::FracModalHits(), nuesand::FillNueSandbox::FracNonContiguous(), cosrej::CosRejFxs::getActivity(), numusand::NumuSandFxs::getActivity(), numusand::NumuSandFxs::getActivityVtx(), numusand::NumuSandFxs::getAveTrackdEdx(), numusand::NumuSandFxs::getAveTrackdEdxLast4Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast6Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast8Cells(), trackinfo::TrackInfoFxs::getAveTrackdEdxLastCells(), cosrej::CosRejFxs::getBBC(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), xsrec::MakeXSecCCPi0Inc::GetBPI(), evd::RecoBaseDrawer::GetClusterOutlines(), evd::RecoBaseDrawer::GetClusterOutlinesConvexHull(), cosrej::CosRejFxs::getCVVars(), cosrej::CosRejFxs::getFits(), xsrec::MakeXSecCCPi0Inc::GetPhLL(), me::MEFinder::GetPhysNoiseSlcs(), jmshower::RecoJMShower::GetPlaneCentroidCell(), jmshower::RecoJMShower::GetPlaneE1Cell(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), jmshower::RecoJMShower::GetRVPStats(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackDedx(), numue::NumuEAlg::GetUCMuonENonSingle(), numue::NumuEAlg::GetUCMuonESingle(), fuzz::FuzzyKValidate::GeVCalc(), slid::NuEEnergyAlg::HadronicDepEnergy(), dt::DiscreteTracker::HighActivity(), nnbar::NNbarUtilities::hitCountInView(), rb::HitList::HitList(), cvn::HitNuIndex(), jmshower::hitSort(), tf::TimingFitAlg::HoughFit(), LSTME::BasicProngVars::initBasicVars(), hough::HoughT::is_good_slice(), zcl::FastMonopoleTriggers::is_trigger_by_epoch1_fmmtrigger(), trk::CosmicTrackSelection::IsCalibTrack(), calib::RockMuonStopperSelection::IsContainedStopper(), calib::PCHitsList::IsGoodTrack(), trk::CosmicTrackAlg::IsTrackDownstreamFromTiming(), trk::KalmanTrackMerge::JoinTracks(), me::MEFinder::MakeSlcME(), fuzz::ViewMatchAlg::MakeTestProng(), nerd::ViewMatchAlg::MakeTestProng(), trk::CosmicTrackAlg::MakeTrack(), trk::WindowTrackingAlg::MakeTrack(), airshower::AirSlicer::MakeTrackSlices(), me::MEFinder::MakeTrkME(), murem::FindParent::MatchToOrigSlice(), nuesand::FillNueSandbox::MaxGap(), nuesand::FillNueSandbox::MaxHits(), dt::DiscreteTracker::MergeParasiteTracks(), nerd::NERDProng::MinHitClusterDist(), me::MEFinder::MinHitSlcDist(), dt::ViewMerger::MissingChargeByExtension(), murem::MuonRemove::MuonInfo(), dif::DiFShowerFinder::muonstub(), numue::NumuEAlg::NDEnergy(), numue::NumuEAlg::NDTrackEnergySplitting(), zcl::FmmTriggerEvd::NumberOfCellsPerLength(), zcl::FastMonopoleTriggers::NumberOfHitsInOverlapPlanesCut(), zcl::FmmTriggerEvd::NumberOfHitsInOverlapPlanesCut(), zcl::FastMMStudy::NumberOfHitsInOverlapPlanesCut(), zcl::FastMonopoleTriggers::NumberOfSurfaceHits(), zcl::FastMMStudy::NumberOfSurfaceHits(), rb::operator<<(), air::AirFilter::parallel_muon_slice(), skim::ParametersNumu::ParametersNumu(), zcl::FastMonopoleTriggers::passed_epoch2_fmmtrigger(), me::MEFinder::PassesMEPresel(), trident::DimuonSliceAna::PlaneCellCounts(), calib::PCHitsList::ProcessTrack(), calib::PCHitsList::ProcessTrackTrajectory(), trk::TrkAssn::produce(), vf::VertexFinder::produce(), qeef::QeFinder::produce(), hv::HoughVertex::produce(), nuesand::FillNueSandbox::produce(), wsnumu::WSNumu::produce(), slid::Recluster::produce(), presel::CosmicVeto::produce(), presel::CosmicCVNVeto::produce(), presel::NuePresel::produce(), lem::Preselection::produce(), numusand::FillSandbox::produce(), slid::SliceLIDBuilder::produce(), zcl::SPCluster::produce(), cosrej::CosRej::produce(), numue::NumuEnergy::produce(), lem::MakeLibrary::produce(), earms::ElasticArmsHS::produce(), dt::DiscreteTracker::produce(), ndreco::NDReco::produce(), calib::PCHitsList::produce(), sn::SNSlicer::produce(), sn::SupernovaMCCluster::produce(), vdt::VertexDT::produce(), zcl::FmmTriggerEvd::produce(), jmshower::NueSel::produce(), dif::DiFShowerFinder::produce(), rvp::RecVarPID::produce(), jmshower::NueSelLID::produce(), bpfit::BPFEnergyEstimator::produce(), bpfit::BPFEnergyEstimatorOnly::produce(), cosrej::MakeNueCosRej::produce(), ncid::NCNNKeras::produce(), murem::MuonRemove::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), caf::CAFMaker::produce(), slid::ParticleIDAlg::Radius(), jmshower::RecoJMShower::Radius(), jmshower::RecoJMShower::RadiusV(), jmshower::RecoJMShower::RecoShowers(), murem::MuonRemove::RemoveByTruth(), trk::KalmanTrack::RemoveHitsFromSignal(), dt::ViewMerger::ScoreForJoinPlusExtension(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), slid::NuEEnergyAlg::ShowerDepEnergy(), lem::MakeLibrary::SliceMeanPosEWeighted(), lem::MakeLibrary::SliceOrigin(), lem::LEMSummarizer::SliceToLEMInput(), remid::ReMIdValidate::SortByHits(), remid::ReMIdTrain::SortByHits(), htk::sortByNCell(), zcl::FmmTriggerEvd::StdevCellsPerPlane(), nnbar::NNbarUtilities::temporalExpand(), calib::StopperThreshold::testPath(), htk::Track3D::time_fit(), mono::Track3D::time_fit(), dt::ViewMerger::TotalChargePerView(), nnbar::NNbarUtilities::totalVisibleEnergyInView(), murem::TrackCleanUpAlg::TrackEinMIP(), zcl::FMMTracker::Tracking2D(), dprf::ChannelPlots::TrackOK(), dprf::TrackPlots::TrackOK(), trk::CosmicTrackUtilities::TrackResiduals(), slid::NuEEnergyAlg::VertexEnergy(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    switch(view){
    case geo::kX: return NXCell();
    case geo::kY: return NYCell();
    case geo::kXorY: return NCell();
    default: assert(0 && "Unknown view");
    }
  }

unsigned int rb::Cluster::NCell ( ) const

inlineinherited

Number of cells in either view.

Definition at line 110 of file Cluster.h.

References rb::Cluster::Cell(), rb::Cluster::NXCell(), rb::Cluster::NYCell(), rb::Cluster::XCell(), and rb::Cluster::YCell().

Referenced by rb::Cluster::Exclude(), jmshower::JMShower::GetPlaneHits(), jmshower::JMShower::GetTransHits(), rb::Cluster::MaxCell(), rb::Cluster::MaxPlane(), rb::Cluster::MaxTNS(), rb::Cluster::MeanTNS(), rb::Cluster::MinCell(), rb::Cluster::MinMaxMeanXYZ(), rb::Cluster::MinPlane(), rb::Cluster::MinTNS(), rb::Cluster::MostContiguousPlanes(), rb::Cluster::MostMissingPlanes(), rb::Cluster::NCell(), rb::Cluster::NMissingPlanes(), rb::Cluster::OfflineChans(), PlaneDirMap(), rb::Cluster::RemoveHit(), rb::Cluster::TotalADC(), rb::Cluster::TotalGeV(), rb::Prong::TotalLength(), rb::Cluster::TotalPE(), rb::Cluster::TotalWeight(), rb::Cluster::View(), and rb::Cluster::WeightedHits().

{return NXCell()+NYCell();}

int rb::Cluster::NMissingPlanes ( geo::View_t  view ) const

inherited

Total number of missing planes in cluster.

Put simply, number of planes spanned minus number of hit planes. If track is 2D, only planes in that view are counted.

Definition at line 693 of file Cluster.cxx.

References rb::Cluster::Cell(), MECModelEnuComparisons::i, geo::kXorY, rb::Cluster::NCell(), update_sam_good_runs_metadata::stride, and rb::Cluster::View().

Referenced by bpf::BPFCVNAna::analyze(), rb::Cluster::ExtentTNS(), caf::FillTrackVars(), skim::ParametersNue::ParametersNue(), wsnumu::WSNumu::produce(), and slid::SliceLIDBuilder::produce().

  {
    if(view == geo::kXorY) view = View();

    const int stride = (view == geo::kXorY) ? 1 : 2;
    
    std::set<int> planes;

    for(unsigned int i = 0; i < NCell(view); ++i)
      planes.insert(Cell(view, i)->Plane());

    if(planes.empty()) return 0;
    
    return (*planes.rbegin() - *planes.begin())/stride - (int)planes.size() + 1; 

  }

size_t rb::Track::NTrajectoryPoints ( ) const

inline

Definition at line 83 of file Track.h.

References fTraj.

Referenced by trk::KalmanTrackMerge::CalcMatchScore(), trk::KalmanTrackMerge::CanJoinTracks(), trk::WindowTrackingAlg::CheckTrackDirectionInY(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxFD::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), trk::WindowTrackingAlg::DetermineInitialDirection(), evd::RecoBaseDrawer::DrawTrack2D(), dt::ViewMerger::FakeThirdDimension(), calib::StopperThreshold::FillHist(), calib::StopperThreshold::FillTree(), trk::KalmanTrack::FindTracks(), calib::StopperThreshold::FindTrajectoryPoints(), cosrej::CosRejFxs::FScatterEstim(), numusand::NumuSandFxs::getScatt(), cosrej::CosRejFxs::getScatt(), ndreco::NDRecoFxs::getTrackDedx(), ndreco::NDRecoFxs::getTrackScatt(), trk::CosmicTrackSelection::GoodSteps(), calib::PCHitsList::GoodSteps(), remid::RecoMuon::IsTrackContained(), trk::KalmanTrackMerge::JoinTracks(), dt::ViewMerger::JoinTracks(), numue::NumuEAlg::NDTrackLength(), cosrej::CosRej::produce(), remid::RecoMuon::produce(), trk::WindowTrackingAlg::SetTrackEndPoints(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), StopDir(), bsf::BremShowerFilter::trackEndContained(), dif::DiFShowerFinder::trackEndContained(), trk::KalmanTrackMerge::ViewMergeTracks(), and ZipWith().

{return fTraj.size();}

unsigned int rb::Cluster::NXCell ( ) const

inlineinherited

Number of cells in the x-view.

Definition at line 106 of file Cluster.h.

References rb::Cluster::fXCell, and art::PtrVector< T >::size().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), chaninfo::CosmicEff::analyze(), comi::CosmicMetrics::analyze(), align::AlignInspect::analyze(), align::Alignment::analyze(), align::SplitTracks::analyze(), sn::SNSlicerAna::analyze(), calib::HitEfficiency::analyze(), calib::PEResponse::analyze(), trk::KalmanTrackAna::analyze(), hough::HoughValidate::analyze(), slid::LIDTraining::analyze(), red::EWCosmics::analyze(), bpfit::DimuonFitter::AnaSlice(), trident::DimuonSliceAna::AnaSlice(), rb::Cluster::Cell(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), sn::SNSlicer::ClusterIsGood(), cvn::PixelMapProducer::CreateMapGivenBoundary(), me::TrkME::DistToTrk(), evd::RecoBaseDrawer::DrawCluster2D(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillHadClustVars(), caf::FillProngVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillShowerVars(), filter::Filter::FillSliceVariables(), caf::FillSliceVars(), sn::SNSlicer::FinalizeCluster(), cvn::Get2DProngEfficiencyByPDG(), cvn::Get2DProngEnergyByPDG(), cvn::Get2DProngPurityByPDG(), jmshower::RecoJMShower::GetCentroid(), cvn::GetProngEfficiencyByPDG(), cvn::GetProngEnergyByPDG(), cvn::GetProngPurityByPDG(), jmshower::RecoJMShower::GetRVPStats(), bpfit::BreakPoint::GoodProng(), LSTME::BasicProngVars::initBasicVars(), trk::CosmicTrackSelection::IsCalibTrack(), calib::PCHitsList::IsGoodTrack(), skim::EvaluatorNue::KeepSlice(), nerd::NERDProng::MakeCluster(), trk::CosmicTrackAlg::MakeTrack(), rb::Cluster::NCell(), skim::ParametersNue::ParametersNue(), skim::ParametersNumu::ParametersNumu(), hough::HoughT::produce(), nuesand::FillNueSandbox::produce(), nussand::FillNusSandbox::produce(), lem::Preselection::produce(), presel::NuePresel::produce(), hough::MultiHoughT::produce(), earms::ElasticArmsHS::produce(), sn::SNSlicer::produce(), jmshower::JMClusterMerge::produce(), sn::SupernovaMCCluster::produce(), vdt::VertexDT::produce(), jmshower::NueSel::produce(), cvntf::CVNProngEvaluatorTF::produce(), rvp::RecVarPID::produce(), jmshower::NueSelLID::produce(), remid::RecoMuon::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), caf::CAFMaker::produce(), cvn::ProngClassify(), jmshower::RecoJMShower::RecoShowers(), rb::Cluster::RemoveHit(), rb::Cluster::SetWeight(), lem::LEMSummarizer::SliceToLEMInput(), rb::Cluster::TotalGeVFastClusterOnly(), dprf::ChannelPlots::TrackOK(), dprf::TrackPlots::TrackOK(), me::MEFinder::TrkMEDist(), rb::Cluster::Weight(), and slicer::MMSlicer::WindowSlice().

{return fXCell.size();}

unsigned int rb::Cluster::NYCell ( ) const

inlineinherited

Number of cells in the y-view.

Definition at line 108 of file Cluster.h.

References rb::Cluster::fYCell, and art::PtrVector< T >::size().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), chaninfo::CosmicEff::analyze(), comi::CosmicMetrics::analyze(), align::AlignInspect::analyze(), align::Alignment::analyze(), align::SplitTracks::analyze(), sn::SNSlicerAna::analyze(), calib::HitEfficiency::analyze(), calib::PEResponse::analyze(), trk::KalmanTrackAna::analyze(), hough::HoughValidate::analyze(), slid::LIDTraining::analyze(), red::EWCosmics::analyze(), bpfit::DimuonFitter::AnaSlice(), trident::DimuonSliceAna::AnaSlice(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), sn::SNSlicer::ClusterIsGood(), cvn::PixelMapProducer::CreateMapGivenBoundary(), me::TrkME::DistToTrk(), evd::RecoBaseDrawer::DrawCluster2D(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillHadClustVars(), caf::FillProngVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillShowerVars(), filter::Filter::FillSliceVariables(), caf::FillSliceVars(), sn::SNSlicer::FinalizeCluster(), cvn::Get2DProngEfficiencyByPDG(), cvn::Get2DProngEnergyByPDG(), cvn::Get2DProngPurityByPDG(), jmshower::RecoJMShower::GetCentroid(), cvn::GetProngEfficiencyByPDG(), cvn::GetProngEnergyByPDG(), cvn::GetProngPurityByPDG(), jmshower::RecoJMShower::GetRVPStats(), bpfit::BreakPoint::GoodProng(), LSTME::BasicProngVars::initBasicVars(), trk::CosmicTrackSelection::IsCalibTrack(), calib::PCHitsList::IsGoodTrack(), skim::EvaluatorNue::KeepSlice(), nerd::NERDProng::MakeCluster(), trk::CosmicTrackAlg::MakeTrack(), rb::Cluster::NCell(), skim::ParametersNue::ParametersNue(), skim::ParametersNumu::ParametersNumu(), hough::HoughT::produce(), nuesand::FillNueSandbox::produce(), nussand::FillNusSandbox::produce(), presel::NuePresel::produce(), hough::MultiHoughT::produce(), lem::Preselection::produce(), earms::ElasticArmsHS::produce(), jmshower::JMClusterMerge::produce(), sn::SNSlicer::produce(), sn::SupernovaMCCluster::produce(), vdt::VertexDT::produce(), jmshower::NueSel::produce(), cvntf::CVNProngEvaluatorTF::produce(), rvp::RecVarPID::produce(), jmshower::NueSelLID::produce(), remid::RecoMuon::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), caf::CAFMaker::produce(), cvn::ProngClassify(), jmshower::RecoJMShower::RecoShowers(), rb::Cluster::RemoveHit(), lem::LEMSummarizer::SliceToLEMInput(), rb::Cluster::TotalGeVFastClusterOnly(), dprf::ChannelPlots::TrackOK(), dprf::TrackPlots::TrackOK(), me::MEFinder::TrkMEDist(), and slicer::MMSlicer::WindowSlice().

{return fYCell.size();}

std::vector< geo::OfflineChan > rb::Cluster::OfflineChans ( ) const

inherited

Positions of all the CellHits.

Definition at line 190 of file Cluster.cxx.

References rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), and runNovaSAM::ret.

Referenced by calib::PCHitsList::ProcessTrack(), calib::PCHitsList::ProcessTrackForBelowThresholdHits(), calib::PCHitsList::ProcessTrackForFLSHits(), and rb::Cluster::YCells().

  {
    std::vector<geo::OfflineChan> ret;
    ret.reserve(NCell());
    for(unsigned int i = 0; i < NCell(); ++i)
      ret.push_back(Cell(i)->OfflineChan());
    return ret;
  }

bool rb::Prong::operator< ( const Prong &  other ) const

inherited

Definition at line 226 of file Prong.cxx.

References rb::Prong::Start().

Referenced by rb::Prong::Dir().

  {
    return this->Start().Z() < other.Start().Z();
  }

bool rb::Cluster::operator< ( const Cluster &  other ) const

inherited

Definition at line 576 of file Cluster.cxx.

References rb::Cluster::MeanZ().

  {
    return this->MeanZ() < other.MeanZ();
  }

std::map< unsigned int, TVector3 > rb::Track::PlaneDirMap

(

)

const

map of the direction cosines at each plane (ie z)

Definition at line 496 of file Track.cxx.

References plot_validation_datamc::c, geo::PlaneGeo::Cell(), rb::Cluster::Cell(), geom(), geo::CellGeo::GetCenter(), InterpolateDir(), rb::Cluster::NCell(), rb::CellHit::Plane(), geo::GeometryBase::Plane(), and NDAPDHVSetting::plane.

Referenced by calib::PCHitsList::ProcessTrack(), calib::PCHitsList::ProcessTrackForBelowThresholdHits(), calib::PCHitsList::ProcessTrackForFLSHits(), Trajectory(), and trk::CosmicTrackUtilities::TriCellPathLength().

  {
    std::map<unsigned int, TVector3> planeDirMap;

    // loop over the planes in the track and get the points on either side
    // of the plane z
    art::ServiceHandle<geo::Geometry> geom;

    // looping over the cells this way does the X and Y
    // view separately, but the map orders by plane
    // number using the < operator, so the entries
    // in the map are linearly increasing, with X and Y
    // views mixed together.
    for(size_t c = 0; c < NCell(); ++c){
      unsigned int plane = Cell(c)->Plane();

      // check if we already have the direction at this plane
      // no need to overwrite it if we do.
      if(planeDirMap.count(plane) > 0) continue;

      // get the z position for this cell
      double xyz[3] = {0.};
      geom->Plane(plane)->Cell(Cell(c)->Cell())->GetCenter(xyz);

      // find the direction based on the points on either
      // side of this z location.
      // set the direction of travel as the direction
      // at this plane.
      planeDirMap[plane] = InterpolateDir(xyz[2]);

    }// end loop over cells

    return planeDirMap;
  }

void rb::Track::PrependTrajectoryPoint

(

TVector3

pt

)

Support constructing tracks backwards.

Note: This function also sets the track start to pt

Definition at line 134 of file Track.cxx.

References rb::Cluster::fPrecalcTotalGeV, rb::Prong::fStart, fTraj, and gen_hdf5record::pt.

Referenced by trk::WindowTrackingAlg::FindEndPoint(), PrependTrajectoryPoint(), and Track().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    fTraj.insert(fTraj.begin(), pt);

    fStart = pt;
  }

void rb::Track::PrependTrajectoryPoint	(	double	pv,
		double	pz
	)

Definition at line 144 of file Track.cxx.

References rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fView, geo::kX, geo::kXorY, and PrependTrajectoryPoint().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    if(fView == geo::kXorY)
      throw cet::exception("Track") << "view is geo::kXorY, cannot add a 2D trajectroy point";

    if(fView == geo::kX)
      PrependTrajectoryPoint(TVector3(pv, 0, pz));
    else
      PrependTrajectoryPoint(TVector3(0, pv, pz));
  }

rb::RecoHit rb::Cluster::RecoHit ( const art::Ptr< rb::CellHit > &  chit ) const

inherited

Return calibrated hit based on assumed W coordinate.

To calibrate a hit you must know its distance from the readout. This is a convenience function which uses the W function to estimate the position of a CellHit and then uses calib::Calibrator to produce a calibrated hit from that information.

Definition at line 259 of file Cluster.cxx.

References art::Ptr< T >::get(), and calib::Calibrator::MakeRecoHit().

Referenced by dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), tut::TutAnalyzer::analyze(), calib::TestParticleCorrections::analyze(), tbem::TBRecoAna::analyze(), trk::KalmanTrackAna::analyze(), showere::ShowerEnergyAna::analyze(), upmuana::UpMuRecoAna::analyze(), upmuana::UpMuAnalysis::analyze(), slid::NuEEnergyAlg::CellEnergy(), murem::TrackCleanUpAlg::CleanUpTrack(), nnbar::NNbarUtilities::correlationTCell(), nnbar::NNbarUtilities::correlationTPlane(), cvn::RegPixelMapProducer::CreateMapGivenBoundary(), cvn::PixelMapProducer::CreateMapGivenBoundary(), cvn::RegPixelMapProducer::CreateMapGivenShowerVertex(), cvn::RegPixelMapProducer::CreateMapGivenVertex(), murem::TrackCleanUpAlg::DeDxInPlane(), nnbar::NNbarUtilities::energyBalancingVertex(), dif::DiFShowerFinder::eparm(), calib::StopperThreshold::FillHist(), trk::CosmicTrackAna::FillTrackHistograms(), calib::StopperThreshold::FillTree(), showere::ShowerEnergyFilterMC::filter(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), numusand::NumuSandFxs::getActivity(), numusand::NumuSandFxs::getAveTrackdEdx(), numusand::NumuSandFxs::getAveTrackdEdxLast4Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast6Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast8Cells(), trackinfo::TrackInfoFxs::getAveTrackdEdxLastCells(), xsrec::MakeXSecCCPi0Inc::GetBPI(), nnbar::NNbarUtilities::getCoordinateOfHits(), cosrej::CosRejFxs::getFits(), nuesand::FillNueSandbox::GetPlaneEnergy(), fuzz::FuzzyKValidate::GeVCalc(), tf::TimingFitAlg::HoughFit(), dif::DiFShowerFinder::muonstub(), numue::NumuEAlg::NDEnergy(), numue::NumuEAlg::NDTrackEnergySplitting(), nussand::FillNusSandbox::produce(), numusand::FillSandbox::produce(), slid::SliceLIDBuilder::produce(), lem::MakeLibrary::produce(), ndreco::NDReco::produce(), remid::RecoMuon::produce(), upmuana::UpMuProb::produce(), caf::CAFMaker::produce(), rb::Cluster::RecoHit(), lem::MakeLibrary::SliceMeanPosEWeighted(), lem::LEMSummarizer::SliceToLEMInput(), nnbar::NNbarUtilities::temporalExpand(), calib::StopperThreshold::testPath(), dt::ViewMerger::TotalChargePerView(), rb::Cluster::TotalGeV(), nnbar::NNbarUtilities::totalVisibleEnergyInView(), murem::TrackCleanUpAlg::TrackEinMIP(), trk::CosmicTrackUtilities::TrackResiduals(), and rb::Cluster::YCells().

  {
    art::ServiceHandle<calib::Calibrator> cal;
    rb::RecoHit rb(cal->MakeRecoHit(*chit, this->W(chit.get())));
    return rb;
  }

rb::RecoHit rb::Cluster::RecoHit ( geo::View_t  view, unsigned int  viewIdx  ) const

inherited

Definition at line 267 of file Cluster.cxx.

References rb::Cluster::Cell(), and rb::Cluster::RecoHit().

  {
    return RecoHit(Cell(view, viewIdx));
  }

rb::RecoHit rb::Cluster::RecoHit ( unsigned int  globalIdx ) const

inherited

Definition at line 273 of file Cluster.cxx.

References rb::Cluster::Cell(), and rb::Cluster::RecoHit().

  {
    return RecoHit(Cell(globalIdx));
  }

void rb::Cluster::RemoveHit ( const art::Ptr< rb::CellHit >  hit )

inherited

Remove hit from current cluster.

Definition at line 290 of file Cluster.cxx.

References art::PtrVector< T >::begin(), art::PtrVector< T >::erase(), rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fXCell, rb::Cluster::fXWeights, rb::Cluster::fYCell, rb::Cluster::fYWeights, MECModelEnuComparisons::i, geo::kX, geo::kY, rb::Cluster::NXCell(), rb::Cluster::NYCell(), and rb::CellHit::View().

Referenced by slid::NuEEnergyAlg::HadronicDepEnergy(), murem::FindParent::MatchToOrigSlice(), slid::Recluster::produce(), sn::SNSlicer::produce(), slid::NuEEnergyAlg::VertexEnergy(), and rb::Cluster::YCells().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    if(hit->View() == geo::kX){
      // try to find the hit in the xcells
      for(unsigned int i = 0;i<NXCell();++i){
        if(fXCell[i] == hit){
          fXCell.erase(fXCell.begin()+i);
          if(!fXWeights.empty()) fXWeights.erase(fXWeights.begin()+i);
          break;
        }
      }
    }
    else if(hit->View() == geo::kY){
     // try to find the hit in the ycells
      for(unsigned int i = 0;i<NYCell();++i){
        if(fYCell[i] == hit){
          fYCell.erase(fYCell.begin()+i);
          if(!fYWeights.empty()) fYWeights.erase(fYWeights.begin()+i);
          break;
        }
      }
    }
  }

void rb::Cluster::RemoveHit ( unsigned int  globalIdx )

inherited

Remove the ith hit from current cluster.

Definition at line 317 of file Cluster.cxx.

References ana::assert(), art::PtrVector< T >::begin(), art::PtrVector< T >::erase(), rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fXCell, rb::Cluster::fYCell, rb::Cluster::NCell(), and rb::Cluster::NXCell().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    assert(globalIdx < NCell());
    if(globalIdx < NXCell()) fXCell.erase(fXCell.begin()+globalIdx);
    else fYCell.erase(fYCell.begin()+globalIdx-NXCell());
  }

void rb::Track::RemoveTrajectoryPoint

(

unsigned int

i

)

Remove the ith trajectory point from the track.

Definition at line 367 of file Track.cxx.

References rb::Cluster::fPrecalcTotalGeV, fTraj, and MECModelEnuComparisons::i.

Referenced by trk::KalmanTrack::FindTracks(), and Trajectory().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    if(i > fTraj.size())
      throw cet::exception("Track") << "attepmt to remove trajectory point " 
                                    << i << " while there are only " << fTraj.size()
                                    << " trajectory points";

    fTraj.erase(fTraj.begin()+i);
  }

void rb::Cluster::SavePrecalcTotalGeV ( ESaveGeVMode  savemode )

inherited

Store the current result of TotalGeV / CalorimetricEnergy.

If you intend to lock in the current calibrations, call this method immediately before storing the product in the event. The result can be retrieved by passing kUsePrecalcEnergy to TotalGeV or CalorimetricEnergy.

Multiple calls are considered to be an error. If you really do need to update the energy, pass kResetTotalGeV.

Definition at line 751 of file Cluster.cxx.

References allTimeWatchdog::endl, stan::math::fabs(), rb::Cluster::fPrecalcTotalGeV, rb::Cluster::kRecomputeEnergy, rb::Cluster::kResetTotalGeV, and rb::Cluster::TotalGeV().

  {
    const double gev = TotalGeV(kRecomputeEnergy);
    if(fPrecalcTotalGeV < 0 || savemode == kResetTotalGeV){
      fPrecalcTotalGeV = gev;
    }
    else{
      if(fabs(gev-fPrecalcTotalGeV) > .01){
        mf::LogError("rb::Cluster") << "SaveTotalGeV() already called. Old value was " << fPrecalcTotalGeV << ", now " << gev
                                    << ". If you intended to reset the stored value, you must pass kResetTotalGeV. Aborting." << std::endl;
        abort();
      }
      else{
        mf::LogWarning("rb::Cluster") << "Called SaveTotalGeV() twice for the same cluster, which is inefficient." << std::endl;
      }
    }
  }

void rb::Prong::SetDir ( TVector3  dir )

virtualinherited

Definition at line 104 of file Prong.cxx.

References rb::Prong::fDir, and rb::Cluster::fPrecalcTotalGeV.

Referenced by slid::ParticleIDAlg::DedxInverseLongLL(), trk::WindowTrackingAlg::DetermineInitialDirection(), rb::Prong::Dir(), bpfit::BreakPoint::FitTracks(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), slid::Recluster::produce(), cheat::MCCheater::produce(), rb::Prong::Prong(), jmshower::RecoJMShower::RecoShowers(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    // I've removed this assert to allow for 2D prongs to set a 3D start
    // direction if they so choose.
    // assert(fView == geo::kXorY);

    fDir = dir.Unit(); // Ensure normalized
  }

void rb::Prong::SetDir ( double  dv, double  dz  )

virtualinherited

Parameters

dv	Component of prong direction in View coordinate
dz	z-component of prong direction

Definition at line 117 of file Prong.cxx.

References ana::assert(), rb::Prong::fDir, rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fView, geo::kX, and geo::kY.

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    assert(fView == geo::kX || fView == geo::kY);

    if(fView == geo::kX){
      fDir = TVector3(dv, 0, dz);
    }
    else{
      fDir = TVector3(0, dv, dz);
    }

    fDir = fDir.Unit(); // Ensure normalized
  }

void rb::Cluster::SetID ( int  id )

inlineinherited

Definition at line 74 of file Cluster.h.

References rb::Cluster::fID.

Referenced by sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), trk::KalmanTrackMerge::MatchTracks(), slid::Recluster::produce(), sn::SNSlicer::produce(), rb::Prong::Prong(), jmshower::RecoJMShower::RecoShowers(), rb::Shower::Shower(), Track(), trk::KalmanTrackMerge::ViewMergeTracks(), and rb::WeightedProng::WeightedProng().

{ fID = id;   }

void rb::Cluster::SetNoise ( bool  noise )

inlineinherited

Declare the cluster to consist of noise hits or not.

Definition at line 161 of file Cluster.h.

References rb::Cluster::fNoiseCluster.

Referenced by rb::Cluster::Exclude(), slicer::Slicer4D::produce(), airshower::AirSlicer::produce(), jmshower::JMClusterMerge::produce(), tdslicer::TDSlicer::produce(), and slicer::Slicer::WindowSlice().

{fNoiseCluster = noise;}

void rb::Track::SetStart ( TVector3  start )

virtual

Reimplemented from rb::Prong.

Definition at line 168 of file Track.cxx.

References rb::Prong::fStart, fTraj, and rb::Prong::SetStart().

Referenced by trk::KalmanTrack::FindTracks(), bpfit::BreakPoint::FitTracks(), trk::KalmanTrackMerge::JoinTracks(), cheat::MCCheater::produce(), Trajectory(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    Prong::SetStart(start);

    if(fTraj.empty()) fTraj.resize(1);
    fTraj[0] = fStart;
  }

void rb::Track::SetStart ( double  v0, double  z0  )

virtual

Reimplemented from rb::Prong.

Definition at line 177 of file Track.cxx.

References rb::Prong::fStart, fTraj, and rb::Prong::SetStart().

  {
    Prong::SetStart(v0, z0);

    if(fTraj.empty()) fTraj.resize(1);
    fTraj[0] = fStart;
  }

void rb::Cluster::SetWeight ( unsigned int  globalIdx, double  weight  )

inherited

Set weight of the cell at this index.

Definition at line 327 of file Cluster.cxx.

References rb::Cluster::fPrecalcTotalGeV, geo::kX, geo::kY, and rb::Cluster::NXCell().

Referenced by rb::Cluster::IsNoise(), rb::Shower::Shower(), and Track().

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    if(globalIdx < NXCell())
      SetWeight(geo::kX, globalIdx, weight);
    else
      SetWeight(geo::kY, globalIdx-NXCell(), weight);
  }

void rb::Cluster::SetWeight ( geo::View_t  view, unsigned int  viewIdx, double  weight  )

inherited

Set weight of the cell at this index in this view.

Definition at line 338 of file Cluster.cxx.

References ana::assert(), rb::Cluster::EnsureWeightAlloc(), rb::Cluster::fPrecalcTotalGeV, rb::Cluster::fXCell, rb::Cluster::fXWeights, rb::Cluster::fYCell, rb::Cluster::fYWeights, geo::kX, art::PtrVector< T >::size(), and ana::weight.

  {
    fPrecalcTotalGeV = -1; // Invalidate any cached value

    assert(weight >= 0 && weight <= 1);

    // Are we optimizing (omitting) the weights storage?
    const bool optWeights = (weight == 1 && fXWeights.empty() && fYWeights.empty());
    if(!optWeights) EnsureWeightAlloc();

    if(view == geo::kX){
      assert(viewIdx < fXCell.size());
      if(!optWeights) fXWeights[viewIdx] = weight;
    }
    else{
      assert(viewIdx < fYCell.size());
      if(!optWeights) fYWeights[viewIdx] = weight;
    }
  }

void rb::Cluster::StandardSort ( )

inherited

Put the cells in the cluster into a standard order.

Helps ensure that two clusters composed of the same cell hits present in exactly the same way.

Definition at line 731 of file Cluster.cxx.

References art::PtrVector< T >::begin(), om::cerr, art::PtrVector< T >::end(), allTimeWatchdog::endl, rb::Cluster::fXCell, rb::Cluster::fXWeights, rb::Cluster::fYCell, rb::Cluster::fYWeights, and rb::standard_compare().

Referenced by rb::Cluster::ExtentTNS().

  {
    //
    // Unweighted clusters are an easy sorting job
    //
    if (fXWeights.size()==0 && fYWeights.size()==0) {
      sort(fXCell.begin(), fXCell.end(), standard_compare);
      sort(fYCell.begin(), fYCell.end(), standard_compare);
      return;
    }
    //
    // \todo Sorting weighted hits requires more complicated code
    //
    std::cerr << __FILE__ << ":" << __LINE__ 
              << " Sort for weighted hits not implemented"
              << std::endl;
    abort();
  }

virtual TVector3 rb::Prong::Start ( ) const

inlinevirtualinherited

Definition at line 73 of file Prong.h.

References rb::Prong::fStart.

Referenced by chaninfo::CosmicEff::analyze(), comi::CosmicMetrics::analyze(), align::AlignInspect::analyze(), align::Alignment::analyze(), murem::MuonRemoveAna::analyze(), align::SplitTracks::analyze(), calib::HitEfficiency::analyze(), tbem::TBRecoAna::analyze(), trk::KalmanTrackAna::analyze(), air::AirKalmanAna::analyze(), calib::CalibAna::analyze(), moonshadowana::MoonShadowAna::analyze(), upmuana::UpMuAnalysis::analyze(), red::EWCosmics::analyze(), slid::ParticleIDAlg::CalcAsymIneria(), slid::ParticleIDAlg::CalcDetectorXYRegionIndex(), trk::KalmanTrackMerge::CalcMatchScore(), slid::ParticleIDAlg::CalcTrkHitPath(), trk::KalmanTrackMerge::CanJoinTracks(), CerenkovSelection::Contained(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), upmuana::UpMuAnalysis::containmentType(), upmuana::UpMuRecoAna::containmentType(), dif::DiFShowerFinder::distancefromtrack(), slid::DeconvolveAlg::DistanceToCore(), vf::DOCAInfo::DOCAInfo(), jmshower::RecoJMShower::Energy(), dt::ViewMerger::FakeThirdDimension(), caf::FillDiFShowerVars(), calib::CosmicTrends::FillHistograms(), caf::FillNumuInfo(), caf::FillProngVars(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillShowerVars(), trk::CosmicTrackAna::FillTrueInfo(), comi::NumiFiltering::filter(), trk::KalmanTrack::FindTracks(), jmshower::RecoJMShower::Gap(), numusand::NumuSandFxs::getActivity(), xsrec::MakeXSecCCPi0Inc::GetBPI(), cosrej::CosRejFxs::getCVVars(), jmshower::RecoJMShower::GetDedxInvLongLL(), jmshower::RecoJMShower::GetDedxLongLL(), jmshower::RecoJMShower::GetDedxTransLL(), cosrej::CosRejFxs::getFits(), slid::ParticleIDAlg::GetGapVertexToShowerStart(), nuesand::FillNueSandbox::GetMincell(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), ndreco::NDRecoFxs::getProngProximity(), jmshower::RecoJMShower::GetShwStop(), air::AirKalmanAna::GetThetaXYZ(), ndreco::NDRecoFxs::getTrackProximity(), sn::TrackActivityRemover::GetTrackType(), sn::SNMichelAnalyzer::GetTrackType(), jmshower::JMShower::GetTransHits(), jmshower::RecoJMShower::GetTrkHitPath(), jmshower::RecoJMShower::GetTrkHitPos(), jmshower::RecoJMShower::GetTrkPlaneCell(), jmshower::RecoJMShower::GetTrkPlaneDistToEdge(), jmshower::RecoJMShower::GetTrkPlanePos(), CerenkovSelection::GoodCosmic(), trk::CosmicTrackSelection::GoodReco(), LSTME::BasicProngVars::initBasicVars(), trk::CosmicTrackSelection::IsCalibTrack(), calib::EnergyStandardCandles::IsContained(), remid::ReMIdValidate::IsContained(), calib::RockMuonStopperSelection::IsContainedStopper(), calib::StopperSelection::IsContainedStopper(), calib::PCHitsList::IsGoodTrack(), calib::ThroughgoingSelection::IsUncontained(), filter::Filter::isVertexInFiducial(), trk::KalmanTrackMerge::JoinTracks(), trk::WindowTrackingAlg::Make3DTrack(), fuzz::ViewMatchAlg::MakeTestProng(), nerd::ViewMatchAlg::MakeTestProng(), dt::ViewMerger::MissingChargeByExtension(), murem::MuonRemove::MuonInfo(), cosrej::MakeNueCosRej::MuonParentByDist(), cosrej::MakeNueCosRej::MuonParentByTime(), rb::Prong::operator<(), air::AirKalmanAna::OppSignTrkCount(), skim::ParametersNumu::ParametersNumu(), slid::ParticleIDAlg::PlaneHitCell(), slid::ParticleIDAlg::PlaneHitXYZ(), slid::ParticleIDAlg::PlaneRadius(), vf::VertexFinder::produce(), slid::Recluster::produce(), presel::CosmicVeto::produce(), calib::FiberCalibration::produce(), cosrej::CosRej::produce(), calib::TimingCalibration::produce(), jmshower::JMClusterMerge::produce(), remid::RecoMuon::produce(), fuzz::FuzzyKVertex::produce(), nerd::NERDProng::produce(), slid::ParticleIDAlg::Radius(), jmshower::RecoJMShower::Radius(), jmshower::RecoJMShower::RecoShowers(), dt::ViewMerger::ScoreForJoinPlusExtension(), trk::WindowTrackingAlg::SetTrackEndPoints(), slid::NuEEnergyAlg::ShowerEnergy(), Track(), bsf::BremShowerFilter::trackStartContained(), dif::DiFShowerFinder::trackStartContained(), air::AirKalmanAna::Unitize(), progbar.ProgressBar::Update(), slid::NuEEnergyAlg::VertexEnergy(), moonshadowana::MoonShadowAna::WriteTrackInfo(), and moonshadowana::MoonShadowAnaHough::WriteTrackInfo().

{return fStart;}

TVector3 rb::Track::Stop

(

)

const

Position of the final trajectory point.

Definition at line 186 of file Track.cxx.

References fTraj.

Referenced by chaninfo::CosmicEff::analyze(), comi::CosmicMetrics::analyze(), align::AlignInspect::analyze(), align::Alignment::analyze(), align::SplitTracks::analyze(), calib::HitEfficiency::analyze(), tbem::TBRecoAna::analyze(), trk::KalmanTrackAna::analyze(), bpfit::BreakPointProtonAna::analyze(), air::AirKalmanAna::analyze(), calib::CalibAna::analyze(), moonshadowana::MoonShadowAna::analyze(), upmuana::UpMuAnalysis::analyze(), sn::SNMichelAnalyzer::analyze(), red::EWCosmics::analyze(), trk::KalmanTrackMerge::CalcMatchScore(), trk::KalmanTrackMerge::CanJoinTracks(), CerenkovSelection::Contained(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxFD::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), upmuana::UpMuAnalysis::containmentType(), upmuana::UpMuRecoAna::containmentType(), me::TrkME::DistToTrk(), calib::CosmicTrends::FillHistograms(), caf::FillNumuInfo(), trk::CosmicTrackAna::FillRecoInfo(), caf::FillTrackVars(), trk::CosmicTrackAna::FillTrueInfo(), comi::NumiFiltering::filter(), trk::KalmanTrack::FindTracks(), calib::StopperThreshold::FindTrajectoryPoints(), numusand::NumuSandFxs::getActivity(), cosrej::CosRejFxs::getCVVars(), cosrej::CosRejFxs::getFits(), nuesand::FillNueSandbox::GetMincell(), air::AirKalmanAna::GetThetaXYZ(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackProximity(), sn::TrackActivityRemover::GetTrackType(), sn::SNMichelAnalyzer::GetTrackType(), CerenkovSelection::GoodCosmic(), trk::CosmicTrackSelection::GoodReco(), trk::CosmicTrackSelection::IsCalibTrack(), calib::EnergyStandardCandles::IsContained(), remid::ReMIdValidate::IsContained(), calib::RockMuonStopperSelection::IsContainedStopper(), calib::StopperSelection::IsContainedStopper(), calib::PCHitsList::IsGoodTrack(), remid::RecoMuon::IsTrackContained(), calib::ThroughgoingSelection::IsUncontained(), filter::Filter::isVertexInFiducial(), trk::KalmanTrackMerge::JoinTracks(), trk::WindowTrackingAlg::Make3DTrack(), dt::ViewMerger::MissingChargeByExtension(), murem::MuonRemove::MuonInfo(), cosrej::MakeNueCosRej::MuonParentByDist(), cosrej::MakeNueCosRej::MuonParentByTime(), air::AirKalmanAna::OppSignTrkCount(), skim::ParametersNumu::ParametersNumu(), presel::CosmicVeto::produce(), calib::FiberCalibration::produce(), cosrej::CosRej::produce(), calib::TimingCalibration::produce(), jmshower::JMClusterMerge::produce(), remid::RecoMuon::produce(), sn::TrackActivityRemover::remove(), dt::ViewMerger::ScoreForJoinPlusExtension(), trk::WindowTrackingAlg::SetTrackEndPoints(), red::OverburdenCalculator::traceBack(), Trajectory(), air::AirKalmanAna::Unitize(), moonshadowana::MoonShadowAna::WriteTrackInfo(), and moonshadowana::MoonShadowAnaHough::WriteTrackInfo().

  {
    return fTraj.back();
  }

TVector3 rb::Track::StopDir ( ) const

virtual

Get direction at the stop point.

If number of trajectory points is less than 2, then can't figure out directions

Definition at line 193 of file Track.cxx.

References release_diff::diff, NTrajectoryPoints(), and TrajectoryPoint().

Referenced by caf::FillTrackVars(), trk::WindowTrackingAlg::SetTrackEndPoints(), and Trajectory().

  {
    const int ntrajectory_points = NTrajectoryPoints();

    /// If number of trajectory points is less than 2,
    /// then can't figure out directions
    if(ntrajectory_points < 2){
        const TVector3 nullvector3;
        return nullvector3;
    }

    const TVector3 stop_point    = TrajectoryPoint(ntrajectory_points - 1);
    const TVector3 stop_point_m1 = TrajectoryPoint(ntrajectory_points - 2);

    const TVector3 diff          = stop_point - stop_point_m1;

    return diff.Unit();
  }

double rb::Cluster::TotalADC ( ) const

inherited

Sum of the ADC of all the contained hits.

Definition at line 360 of file Cluster.cxx.

References rawdata::RawDigit::ADC(), rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), runNovaSAM::ret, and rb::Cluster::Weight().

Referenced by caf::AddSlcMEToVec(), caf::AddTrkMEToVec(), sn::SupernovaAna::analyze(), sn::SNSlicerAna::analyze(), tbem::TBRecoAna::analyze(), air::AirKalmanAna::analyze(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), sn::SNSlicer::ClusterIsGood(), mono::HighEnergySliceRemoval::is_high_energy_slice(), rb::Cluster::IsNoise(), rb::operator<<(), air::AirFilter::parallel_muon_slice(), and skim::ParametersNumu::ParametersNumu().

  {
    double ret = 0;
    for(unsigned int i = 0; i < NCell(); ++i)
      ret += Cell(i)->ADC()*Weight(i);
    return ret;
  }

double rb::Cluster::TotalGeV ( EEnergyCalcScheme  escheme = kRecomputeEnergy ) const

inherited

Simple sum of the estimated GeV of all the hits.

WARNING: This is not a sensible energy estimator by itself. If you just want a simple calorimetric estimate look at CalorimetricEnergy. Otherwise, a range of more sophisticated/elaborate estimators are available elsewhere.

Definition at line 378 of file Cluster.cxx.

References allTimeWatchdog::endl, rb::Cluster::fPrecalcTotalGeV, rb::RecoHit::GeV(), MECModelEnuComparisons::i, rb::RecoHit::IsCalibrated(), rb::Cluster::kUsePrecalcEnergy, rb::Cluster::NCell(), rb::Cluster::RecoHit(), runNovaSAM::ret, rb::Cluster::TotalGeVFastClusterOnly(), and rb::Cluster::Weight().

Referenced by calib::TestParticleCorrections::analyze(), sn::SNSlicerAna::analyze(), vf::TrackEva::analyze(), tbem::TBRecoAna::analyze(), bpfit::BPFTmvaTrainer::analyze(), bpf::BPFCVNAna::analyze(), htk::HoughTrack::analyze(), rb::Cluster::CalorimetricEnergy(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), bpf::BPFCVNAna::CompareE(), nnbar::NNbarUtilities::energyBalancingVertex(), numue::NumuEAlg::FDEnergy(), caf::FillNueEnergyVars(), trk::CosmicTrackAna::FillRecoInfo(), cosrej::CosRejFxs::getActivity(), numusand::NumuSandFxs::getActivityVtx(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackDedx(), mono::MonopoleCluster::is_cosmic_track(), numue::NumuEAlg::NDEnergy(), skim::ParametersNue::ParametersNue(), skim::ParametersNumu::ParametersNumu(), qeef::QeFinder::produce(), numusand::FillSandbox::produce(), slid::SliceLIDBuilder::produce(), cosrej::CosRej::produce(), numue::NumuEnergy::produce(), sn::SupernovaMCCluster::produce(), bpfit::BPFEnergyEstimator::produce(), bpfit::BPFEnergyEstimatorOnly::produce(), cvn::ProngClassify(), and rb::Cluster::SavePrecalcTotalGeV().

  {
    if(escheme == kUsePrecalcEnergy){
      if(fPrecalcTotalGeV >= 0) return fPrecalcTotalGeV;

      mf::LogWarning("rb::Cluster::TotalGeV()") << "Passed kUsePrecalcEnergy but no energy was cached in the Cluster. Falling back to recomputing." << std::endl;
    }

    // Optimization that's only safe for literal Clusters
    if(typeid(*this) == typeid(rb::Cluster)) return TotalGeVFastClusterOnly();

    // For derived classes, fall through to the generic implementation

    double ret = 0;
    for(unsigned int i = 0; i < NCell(); ++i){
      const rb::RecoHit rhit = RecoHit(i);
      if(rhit.IsCalibrated()) ret += rhit.GeV()*Weight(i);
    }
    return ret;
  }

double rb::Track::TotalLength ( ) const

virtual

Length (cm) of all the track segments.

ie This measure follows any curvature in the track

Reimplemented from rb::Prong.

Definition at line 213 of file Track.cxx.

References dist, fTraj, Mag(), and getGoodRuns4SAM::n.

Referenced by murem::MuonRemoveAna::analyze(), vf::TrackEva::analyze(), trk::KalmanTrackAna::analyze(), bpf::BPFCVNAna::analyze(), bpfit::BreakPointProtonAna::analyze(), air::AirKalmanAna::analyze(), calib::CalibAna::analyze(), moonshadowana::MoonShadowAna::analyze(), upmuana::UpMuRecoAna::analyze(), upmuana::UpMuAnalysis::analyze(), trk::KalmanTrackAna::CheckRecoTracks(), dt::CompareByLength(), dt::CompareByLengthVM(), murem::TrackCleanUpAlg::DeDxInPlane(), trk::CosmicTrackAna::FillRecoInfo(), CerenkovSelection::FillTree(), numusand::NumuSandFxs::getAveTrackdEdx(), mono::MonopoleCluster::is_cosmic_track(), calib::RockMuonStopperSelection::IsContainedStopper(), novaddt::HorizontalMuonAna::IsHorzMuonTrack(), nuesand::FillNueSandbox::produce(), trackinfo::TrackInfo::produce(), numusand::FillSandbox::produce(), cosrej::CosRej::produce(), rvp::RecVarPID::produce(), remid::RecoMuon::produce(), trk::CosmicTrackUtilities::RangeMomentum(), caf::CAFMaker::sortRBTrackLength(), and Trajectory().

  {
    if(fTraj.size() == 1) return 0;

    // We take the sum of the straight lines between the control points
    double dist = 0;

    const int N = fTraj.size();
    for(int n = 0; n < N-1; ++n){
      dist += (fTraj[n+1]-fTraj[n]).Mag();
    }

    return dist;
  }

double rb::Cluster::TotalPE ( ) const

inherited

Sum of the PE value of all the contained hits.

Definition at line 369 of file Cluster.cxx.

References rb::Cluster::Cell(), MECModelEnuComparisons::i, rb::Cluster::NCell(), rb::CellHit::PE(), runNovaSAM::ret, and rb::Cluster::Weight().

Referenced by rb::Cluster::IsNoise().

  {
    double ret = 0;
    for(unsigned int i = 0; i < NCell(); ++i)
      ret += Cell(i)->PE()*Weight(i);
    return ret;
  }

double rb::Cluster::TotalWeight ( ) const

inherited

Sum of all the weights. The effective number of hits.

Definition at line 430 of file Cluster.cxx.

References rb::Cluster::NCell(), runNovaSAM::ret, and rb::Cluster::Weight().

  {
    double ret = 0;
    for(unsigned int cellIdx = 0; cellIdx < NCell(); ++cellIdx)
      ret += Weight(cellIdx);
    return ret;
  }

std::vector<TVector3> const& rb::Track::Trajectory ( ) const

inline

return a constant reference to the track trajectory points

Definition at line 86 of file Track.h.

References a, b, ClearTrajectoryPoints(), DistanceFromEnd(), DistanceFromStart(), FindNeighbouringPointIndices(), FindNeighbouringPoints(), fTraj, MECModelEnuComparisons::i, InterpolateDir(), InterpolatePts(), InterpolateXY(), PlaneDirMap(), RemoveTrajectoryPoint(), SetStart(), Stop(), StopDir(), TotalLength(), TrajectoryPoint(), W(), submit_syst::x, submit_syst::y, test::z, and ZipWith().

Referenced by bpfit::dEdxCalculator::computeDEDX(), sn::TrackActivityRemover::DistanceHitToTrackBody(), calib::StopperThreshold::FillHist(), trk::CosmicTrackAna::FillRecoInfo(), trk::CosmicTrackAna::FillTrackHistograms(), calib::StopperThreshold::FillTree(), numusand::NumuSandFxs::getAveTrackdEdxLast4Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast6Cells(), numusand::NumuSandFxs::getAveTrackdEdxLast8Cells(), trackinfo::TrackInfoFxs::getAveTrackdEdxLastCells(), trk::WindowTrackingAlg::LookForBremsstrahlungHits(), calib::PCHitsList::ProcessTrackTrajectory(), trk::CosmicTrackUtilities::RangeMomentum(), trk::WindowTrackingAlg::SetTrackEndPoints(), calib::StopperThreshold::testPath(), and trk::CosmicTrackUtilities::TrackResiduals().

{ return fTraj; }

TVector3 rb::Track::TrajectoryPoint

(

unsigned int

i

)

const

The ith point on the trajectory, a 3-vector in cm.

Definition at line 158 of file Track.cxx.

References fTraj, and MECModelEnuComparisons::i.

Referenced by trk::KalmanTrackAna::analyze(), trk::KalmanTrackMerge::CalcMatchScore(), trk::KalmanTrackMerge::CanJoinTracks(), trk::WindowTrackingAlg::CheckTrackDirectionInY(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxFD::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), trk::WindowTrackingAlg::DetermineInitialDirection(), evd::RecoBaseDrawer::DrawTrack2D(), dt::ViewMerger::FakeThirdDimension(), calib::StopperThreshold::FillHist(), calib::StopperThreshold::FillTree(), trk::KalmanTrack::FindTracks(), calib::StopperThreshold::FindTrajectoryPoints(), cosrej::CosRejFxs::FScatterEstim(), numusand::NumuSandFxs::getScatt(), cosrej::CosRejFxs::getScatt(), ndreco::NDRecoFxs::getTrackDedx(), ndreco::NDRecoFxs::getTrackScatt(), trk::CosmicTrackSelection::GoodSteps(), calib::PCHitsList::GoodSteps(), remid::RecoMuon::IsTrackContained(), trk::KalmanTrackMerge::JoinTracks(), dt::ViewMerger::JoinTracks(), numue::NumuEAlg::NDTrackLength(), cosrej::CosRej::produce(), remid::RecoMuon::produce(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), StopDir(), bsf::BremShowerFilter::trackEndContained(), dif::DiFShowerFinder::trackEndContained(), Trajectory(), and trk::KalmanTrackMerge::ViewMergeTracks().

  {
    if(i > fTraj.size())
      throw cet::exception("Track") << "attempt to get trajectory point " 
                                    << i << " is out of range: " << fTraj.size();

    return fTraj[i];
  }

virtual geo::View_t rb::Cluster::View ( ) const

inlinevirtualinherited

kXorY for 3D clusters.

Definition at line 99 of file Cluster.h.

References rb::Cluster::fView, and rb::Cluster::NCell().

Referenced by nerd::NERDProng::AddProngs(), dt::ViewMerger::Adjacent(), calib::HitEfficiency::analyze(), trk::KalmanTrackAna::analyze(), trk::KalmanTrackMerge::CalcMatchScore(), trk::KalmanTrackMerge::CanJoinTracks(), cheat::RecoCheckAna::CheckRecoClusters(), trk::KalmanTrackAna::CheckRecoTracks(), skim::SkimmingUtils::CopySlice(), slid::DeconvolveAlg::DistanceToCore(), evd::RecoBaseDrawer::DrawTrack2D(), dt::ViewMerger::FakeThirdDimension(), trk::KalmanTrackMerge::JoinTracks(), dt::ViewMerger::JoinTracks(), trk::WindowTrackingAlg::LookForBremsstrahlungHits(), trk::WindowTrackingAlg::Make3DTrack(), rb::Cluster::MinMaxMeanXYZ(), dt::ViewMerger::MissingChargeByExtension(), rb::Cluster::NMissingPlanes(), rb::operator<<(), nussand::FillNusSandbox::produce(), nerd::NERDProng::produce(), dt::ViewMerger::ScoreForJoinPlusExtension(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), trk::WindowTrackingAlg::ShortTrackExtraPlane(), and ZipWith().

{return fView;}

double rb::Track::W ( const rb::CellHit *  chit ) const

virtual

Estimate the unmeasured coordinate of chit.

For Track, this is calculated by taking the position of the straight line defined by the two neighbouring trajectory points at the z-coordinate of the cell. In the case of a hit beyond the end of the track, the line defined by the two hits at the end is extrapolated.

Parameters

chit	Any CellHit, not necessarily owned by this container

Returns

Estimated unmeasured coordinate of this hit (cm) assuming it to be associated with this reconstructed object.

Todo:

CJB: do we want to allow the case with only one hit?

Reimplemented from rb::Prong.

Definition at line 281 of file Track.cxx.

References rb::CellHit::Cell(), rb::Cluster::Cell(), geo::GeometryBase::DetId(), fTraj, geom(), InterpolateXY(), novadaq::cnv::kTESTBEAM, geo::kX, rb::CellHit::Plane(), geo::GeometryBase::Plane(), rb::CellHit::View(), and rb::Prong::W().

Referenced by calib::PCHitsList::CellHitToPCHit(), bpfit::dEdxCalculator::computeDEDX(), evd::RecoBaseDrawer::DrawTrack2D(), calib::PCHitsList::ProcessTrack(), calib::PCHitsList::ProcessTrackTrajectory(), and Trajectory().

  {
    if(fTraj.empty()) return 0;
    /// \todo CJB: do we want to allow the case with only one hit?
    if(fTraj.size() == 1) return Prong::W(chit);

    art::ServiceHandle<geo::Geometry> geom;
    double xyz[3];
    geom->Plane(chit->Plane())->Cell(chit->Cell())->GetCenter(xyz);

    InterpolateXY(xyz[2], xyz[0], xyz[1]);

    if(chit->View() == geo::kX) return xyz[1];
    if(geom->DetId() == novadaq::cnv::kTESTBEAM) return -xyz[0]; // TB horizontal modules (Y-view) have readout on opposite end as ND/FD
    return xyz[0];
    /*
    // Find the points before and after. This part does still assume that
    // searching along Z is the best way to find such points.
    TVector3 prev, next;
    FindNeighbouringPoints(xyz[2], prev, next);

    // Solve the problem of the closest approach of the trajectory to the line
    // made by the cell by projecting down to a 2D trajectory and 2D point for
    // the cell.
    TVector3 start2D = prev;
    TVector3 dir2D = next-prev;
    TVector3 pt2D(xyz);
    const int otherView = 1-chit->View();
    start2D[otherView] = dir2D[otherView] = pt2D[otherView] = 0;

    // Actually find the closest point
    TVector3 closest;
    geo::ClosestApproach(pt2D, start2D, dir2D, closest);

    // Need to solve for the W position of the line, the exact part that we
    // projected away above. So just figure out Z and evaluate the original
    // curve there.
    const double dt = (closest.Z()-prev.Z())/(next-prev).Z();

    return ((1-dt)*prev+dt*next)[otherView];
    */
  }

double rb::Cluster::Weight ( unsigned int  globalIdx ) const

inherited

Weight assigned to the cell.

Definition at line 209 of file Cluster.cxx.

References rb::Cluster::fXWeights, rb::Cluster::fYWeights, geo::kX, geo::kY, and rb::Cluster::NXCell().

Referenced by fuzz::ViewMatchAlg::CalcEnergyProfile(), nerd::ViewMatchAlg::CalcEnergyProfile(), slid::NuEEnergyAlg::CellEnergy(), skim::SkimmingUtils::CopySlice(), rb::Cluster::Exclude(), trk::CosmicTrackAna::FillRecoInfo(), fuzz::FuzzyKValidate::GeVCalc(), fuzz::ViewMatchAlg::MakeTestProng(), nerd::ViewMatchAlg::MakeTestProng(), rb::Cluster::MeanTNS(), rb::Cluster::MinMaxMeanXYZ(), nussand::FillNusSandbox::produce(), rb::Cluster::TotalADC(), rb::Cluster::TotalGeV(), rb::Cluster::TotalGeVFastClusterOnly(), rb::Cluster::TotalPE(), rb::Cluster::TotalWeight(), murem::TrackCleanUpAlg::TrackEinMIP(), rb::Cluster::WeightedHits(), and rb::Cluster::YCells().

  {
    if(fXWeights.empty() && fYWeights.empty()) return 1;

    if(globalIdx < NXCell()) return Weight(geo::kX, globalIdx);
    return Weight(geo::kY, globalIdx-NXCell());
  }

double rb::Cluster::Weight ( geo::View_t  view, unsigned int  viewIdx  ) const

inherited

Weight assigned to the cell.

Definition at line 218 of file Cluster.cxx.

References ana::assert(), rb::Cluster::fXWeights, rb::Cluster::fYWeights, and geo::kX.

  {
    if(fXWeights.empty() && fYWeights.empty()) return 1;

    if(view == geo::kX){
      assert(viewIdx < fXWeights.size());
      return fXWeights[viewIdx];
    }
    else{
      assert(viewIdx < fYWeights.size());
      return fYWeights[viewIdx];
    }
  }

std::vector< rb::WeightedHit > rb::Cluster::WeightedHits ( ) const

inherited

Get all hits from both views, with weights attached.

Definition at line 200 of file Cluster.cxx.

References rb::Cluster::Cell(), rb::Cluster::NCell(), and rb::Cluster::Weight().

Referenced by rb::Cluster::YCells().

  {
    std::vector<rb::WeightedHit> all;
    for(unsigned int cellIdx = 0; cellIdx < NCell(); ++cellIdx)
      all.push_back(rb::WeightedHit(Cell(cellIdx), Weight(cellIdx)));
    return all;
  }

art::Ptr< rb::CellHit > rb::Cluster::XCell ( unsigned int  xIdx ) const

inherited

Get the ith cell in the x-view.

Definition at line 157 of file Cluster.cxx.

References ana::assert(), rb::Cluster::fXCell, and art::PtrVector< T >::size().

Referenced by align::SplitTracks::analyze(), trk::KalmanTrackAna::analyze(), zcl::FastMMStudy::analyze(), slid::LIDTraining::analyze(), rb::Cluster::Cell(), evd::RecoBaseDrawer::DrawCluster2D(), caf::FillHadClustVars(), caf::FillSliceVars(), jmshower::RecoJMShower::GetCentroid(), jmshower::RecoJMShower::GetRVPStats(), skim::EvaluatorNue::KeepSlice(), nerd::NERDProng::MakeCluster(), trk::CosmicTrackAlg::MakeTrack(), rb::Cluster::NCell(), nuesand::FillNueSandbox::produce(), nussand::FillNusSandbox::produce(), lem::Preselection::produce(), presel::NuePresel::produce(), vdt::VertexDT::produce(), jmshower::NueSel::produce(), rvp::RecVarPID::produce(), jmshower::NueSelLID::produce(), remid::RecoMuon::produce(), caf::CAFMaker::produce(), jmshower::RecoJMShower::RecoShowers(), lem::LEMSummarizer::SliceToLEMInput(), and rb::Cluster::TotalGeVFastClusterOnly().

  {
    assert(xIdx < fXCell.size());

    return fXCell[xIdx];
  }

const art::PtrVector<rb::CellHit>& rb::Cluster::XCells ( ) const

inlineinherited

Get all cells from the x-view.

Definition at line 124 of file Cluster.h.

References rb::Cluster::fXCell.

Referenced by upmuana::UpMuAnalysis::analyze(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillProngVars(), caf::FillShowerVars(), mono::MonopoleCluster::find_isolated_hits(), zcl::FMMTracker::produce(), cvn::ProngClassify(), and mono::Track3D::Track3D().

{ return fXCell; }

art::Ptr< rb::CellHit > rb::Cluster::YCell ( unsigned int  yIdx ) const

inherited

Get the ith cell in the y-view.

Definition at line 165 of file Cluster.cxx.

References ana::assert(), rb::Cluster::fYCell, and art::PtrVector< T >::size().

Referenced by align::SplitTracks::analyze(), trk::KalmanTrackAna::analyze(), slid::LIDTraining::analyze(), rb::Cluster::Cell(), evd::RecoBaseDrawer::DrawCluster2D(), caf::FillHadClustVars(), caf::FillSliceVars(), jmshower::RecoJMShower::GetCentroid(), jmshower::RecoJMShower::GetRVPStats(), skim::EvaluatorNue::KeepSlice(), nerd::NERDProng::MakeCluster(), trk::CosmicTrackAlg::MakeTrack(), rb::Cluster::NCell(), nuesand::FillNueSandbox::produce(), nussand::FillNusSandbox::produce(), lem::Preselection::produce(), presel::NuePresel::produce(), vdt::VertexDT::produce(), jmshower::NueSel::produce(), rvp::RecVarPID::produce(), jmshower::NueSelLID::produce(), remid::RecoMuon::produce(), caf::CAFMaker::produce(), jmshower::RecoJMShower::RecoShowers(), lem::LEMSummarizer::SliceToLEMInput(), and rb::Cluster::TotalGeVFastClusterOnly().

  {
    assert(yIdx < fYCell.size());

    return fYCell[yIdx];
  }

const art::PtrVector<rb::CellHit>& rb::Cluster::YCells ( ) const

inlineinherited

Get all cells from the x-view.

Definition at line 126 of file Cluster.h.

References rb::Cluster::AllCells(), rb::Cluster::Clear(), rb::Cluster::Exclude(), rb::Cluster::fYCell, rb::Cluster::OfflineChans(), rb::Cluster::RecoHit(), rb::Cluster::RemoveHit(), rb::Cluster::Weight(), and rb::Cluster::WeightedHits().

Referenced by upmuana::UpMuAnalysis::analyze(), caf::FillDiFShowerVars(), caf::FillDiFVars(), caf::FillProngVars(), caf::FillShowerVars(), mono::MonopoleCluster::find_isolated_hits(), zcl::FMMTracker::produce(), cvn::ProngClassify(), and mono::Track3D::Track3D().

{ return fYCell; }

rb::Track rb::Track::ZipWith

(

const rb::Track &

trk

)

const

Combine with a 2D track from the other view to make a 3D track.

Definition at line 442 of file Track.cxx.

References a, rb::Cluster::AllCells(), b, fTraj, InterpolateXY(), rb::Cluster::Is2D(), geo::kX, geo::kY, NTrajectoryPoints(), gen_hdf5record::pt, runNovaSAM::ret, std::swap(), rb::Cluster::View(), submit_syst::x, and submit_syst::y.

Referenced by trk::WindowTrackingAlg::Make3DTrack(), Trajectory(), and dt::ViewMerger::ZipTracks().

  {
    if(!Is2D() || !trk.Is2D())
      throw cet::exception("Track") << "attempting to zip a tracks with unknown dimensions; this track: "
      << this->View() << " zipping track: " << trk.View();
    
    if(this->NTrajectoryPoints() < 2 || trk.NTrajectoryPoints() < 2)
      throw cet::exception("Track") << "too few points to attempt ZipWith: "
      << this->NTrajectoryPoints() << " "
      << trk.NTrajectoryPoints();
    
    const rb::Track* xtrk = this;
    const rb::Track* ytrk = &trk;
    
    if(this->View() != geo::kX) std::swap(xtrk, ytrk);
    
    if(xtrk->View() != geo::kX || ytrk->View() != geo::kY)
      throw cet::exception("Track") << "views of 2D tracks are wrong";
    
    
    std::vector<TVector3> zipTraj;
    zipTraj.reserve(xtrk->NTrajectoryPoints() + ytrk->NTrajectoryPoints());
    // Give each point a 3rd coordinate interpolated from the other track
    for(const TVector3& pt: xtrk->fTraj){
      double junk, y;
      ytrk->InterpolateXY(pt.Z(), junk, y);
      zipTraj.emplace_back(pt.X(), y, pt.Z());
    }
    for(const TVector3& pt: ytrk->fTraj){
      double junk, x;
      xtrk->InterpolateXY(pt.Z(), x, junk);
      zipTraj.emplace_back(x, pt.Y(), pt.Z());
    }
    
    // Sort output points upstream to downstream. User could reverse result if
    // they wanted.
    std::sort(zipTraj.begin(), zipTraj.end(),
              [](const TVector3& a, const TVector3& b)
              {
                return a.Z() < b.Z();
              });
    
    Track ret(xtrk->AllCells(), zipTraj[0], zipTraj[1]-zipTraj[0]);
    ret.Add(ytrk->AllCells());
    
      // First point is already stored as the vertex
    for(size_t p = 1; p < zipTraj.size(); ++p){
      ret.AppendTrajectoryPoint(zipTraj[p]);
    }
    
    return ret;
  }

Member Data Documentation

TVector3 rb::Prong::fDir

protectedinherited

Direction at starting point.

Definition at line 112 of file Prong.h.

Referenced by rb::Prong::Dir(), rb::Prong::DistanceFromStart(), rb::Prong::SetDir(), rb::Shower::Stop(), rb::Prong::TotalLength(), and rb::Prong::W().

int rb::Cluster::fID

protectedinherited

ID for cluster.

Definition at line 321 of file Cluster.h.

Referenced by rb::Cluster::Exclude(), rb::Cluster::ID(), and rb::Cluster::SetID().

bool rb::Cluster::fNoiseCluster

protectedinherited

flag for whether this is a noise cluster

Definition at line 322 of file Cluster.h.

Referenced by rb::Cluster::Exclude(), rb::Cluster::IsNoise(), and rb::Cluster::SetNoise().

double rb::Cluster::fPrecalcTotalGeV

protectedinherited

-1 = uninitialized

Definition at line 324 of file Cluster.h.

Referenced by rb::Cluster::Add(), AppendTrajectoryPoint(), rb::Cluster::Clear(), ClearTrajectoryPoints(), PrependTrajectoryPoint(), rb::Cluster::RemoveHit(), RemoveTrajectoryPoint(), rb::Cluster::SavePrecalcTotalGeV(), rb::Prong::SetDir(), rb::Prong::SetStart(), rb::Cluster::SetWeight(), and rb::Cluster::TotalGeV().

TVector3 rb::Prong::fStart

protectedinherited

Start location (xyz, cm)

Definition at line 111 of file Prong.h.

Referenced by DistanceFromEnd(), rb::Prong::DistanceFromStart(), DistanceFromStart(), PrependTrajectoryPoint(), rb::Prong::SetStart(), SetStart(), rb::Prong::Start(), rb::Shower::Stop(), rb::Prong::TotalLength(), and rb::Prong::W().

std::vector<TVector3> rb::Track::fTraj

protected

Definition at line 160 of file Track.h.

Referenced by AppendTrajectoryPoint(), ClearTrajectoryPoints(), DistanceFromEnd(), DistanceFromStart(), FindNeighbouringPointIndices(), FindNeighbouringPoints(), NTrajectoryPoints(), PrependTrajectoryPoint(), RemoveTrajectoryPoint(), SetStart(), Stop(), TotalLength(), Trajectory(), TrajectoryPoint(), W(), and ZipWith().

geo::View_t rb::Cluster::fView

protectedinherited

view this cluster is in

Definition at line 316 of file Cluster.h.

Referenced by rb::Cluster::Add(), AppendTrajectoryPoint(), rb::Cluster::Cluster(), rb::Cluster::Exclude(), rb::Prong::Is3D(), PrependTrajectoryPoint(), rb::Prong::SetDir(), rb::Prong::SetStart(), Track(), and rb::Cluster::View().

art::PtrVector<rb::CellHit> rb::Cluster::fXCell

protectedinherited

collection of x-view cells in cluster

Definition at line 317 of file Cluster.h.

Referenced by rb::Cluster::Add(), rb::Cluster::AllCells(), rb::Cluster::Clear(), rb::Cluster::Cluster(), rb::Cluster::EnsureWeightAlloc(), rb::Cluster::Is3D(), rb::Cluster::NXCell(), rb::Cluster::RemoveHit(), rb::Cluster::SetWeight(), rb::Cluster::StandardSort(), rb::Cluster::XCell(), and rb::Cluster::XCells().

std::vector<double> rb::Cluster::fXWeights

protectedinherited

Weights, matching cell indexing.

Definition at line 319 of file Cluster.h.

Referenced by rb::Cluster::Add(), rb::Cluster::Clear(), rb::Cluster::EnsureWeightAlloc(), rb::Cluster::RemoveHit(), rb::Cluster::SetWeight(), rb::Cluster::StandardSort(), and rb::Cluster::Weight().

art::PtrVector<rb::CellHit> rb::Cluster::fYCell

protectedinherited

collection of y-view cells in cluster

Definition at line 318 of file Cluster.h.

Referenced by rb::Cluster::Add(), rb::Cluster::AllCells(), rb::Cluster::Clear(), rb::Cluster::Cluster(), rb::Cluster::EnsureWeightAlloc(), rb::Cluster::Is3D(), rb::Cluster::NYCell(), rb::Cluster::RemoveHit(), rb::Cluster::SetWeight(), rb::Cluster::StandardSort(), rb::Cluster::YCell(), and rb::Cluster::YCells().

std::vector<double> rb::Cluster::fYWeights

protectedinherited

May be empty, means all weights are 1.

Definition at line 320 of file Cluster.h.

Referenced by rb::Cluster::Add(), rb::Cluster::Clear(), rb::Cluster::EnsureWeightAlloc(), rb::Cluster::RemoveHit(), rb::Cluster::SetWeight(), rb::Cluster::StandardSort(), and rb::Cluster::Weight().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-07/RecoBase/Track.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-07/RecoBase/Track.cxx