Public Member Functions | Private Member Functions | Private Attributes | List of all members
ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder Class Reference
Inheritance diagram for ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder:
ShowerRecoTools::IShowerTool

Public Member Functions

 ShowerPandoraSlidingFitTrackFinder (const fhicl::ParameterSet &pset)
 
int CalculateElement (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder) override
 

Private Member Functions

void InitialiseProducers () override
 
int AddAssociations (const art::Ptr< recob::PFParticle > &pfpPtr, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder) override
 
- Private Member Functions inherited from ShowerRecoTools::IShowerTool
 IShowerTool (const fhicl::ParameterSet &pset)
 
virtual ~IShowerTool () noexcept=default
 
int RunShowerTool (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder, std::string evd_display_name_append="")
 
void SetPtr (art::ProducesCollector *collector)
 
void InitaliseProducerPtr (reco::shower::ShowerProducedPtrsHolder &uniqueproducerPtrs)
 
const shower::LArPandoraShowerAlg & GetLArPandoraShowerAlg () const
 
template<class T >
art::Ptr< T > GetProducedElementPtr (std::string Name, reco::shower::ShowerElementHolder &ShowerEleHolder, int iter=-1)
 
template<class T >
void InitialiseProduct (std::string Name, std::string InstanceName="")
 
template<class T , class A , class B >
void AddSingle (A &a, B &b, std::string Name)
 
int GetVectorPtrSize (std::string Name)
 
void PrintPtrs ()
 
void PrintPtr (std::string Name)
 

Private Attributes

art::ServiceHandle< geo::Geometry > fGeom
 
int fVerbose
 
float fSlidingFitHalfWindow
 
float fMinTrajectoryPoints
 
std::string fInitialTrackOutputLabel
 
std::string fInitialTrackLengthOutputLabel
 
std::string fShowerStartPositionInputLabel
 
std::string fShowerDirectionInputLabel
 
std::string fInitialTrackSpacePointsInputLabel
 
std::string fInitialTrackHitsInputLabel
 

Detailed Description

Definition at line 22 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::ShowerPandoraSlidingFitTrackFinder ( const fhicl::ParameterSet &  pset)

Definition at line 54 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

56  : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
57  , fVerbose(pset.get<int>("Verbose"))
58  , fSlidingFitHalfWindow(pset.get<float>("SlidingFitHalfWindow"))
59  , fMinTrajectoryPoints(pset.get<float>("MinTrajectoryPoints"))
60  , fInitialTrackOutputLabel(pset.get<std::string>("InitialTrackOutputLabel"))
61  , fInitialTrackLengthOutputLabel(pset.get<std::string>("InitialTrackLengthOutputLabel"))
62  , fShowerStartPositionInputLabel(pset.get<std::string>("ShowerStartPositionInputLabel"))
63  , fShowerDirectionInputLabel(pset.get<std::string>("ShowerDirectionInputLabel"))
64  , fInitialTrackSpacePointsInputLabel(pset.get<std::string>("InitialTrackSpacePointsInputLabel"))
65  , fInitialTrackHitsInputLabel(pset.get<std::string>("InitialTrackHitsInputLabel"))
66  {}
string
std::string string
Definition: nybbler.cc:12
fhicl::ParameterSet::get
T get(std::string const &key) const
Definition: ParameterSet.h:271
ShowerRecoTools::IShowerTool::IShowerTool
IShowerTool(const fhicl::ParameterSet &pset)
Definition: IShowerTool.h:28

Member Function Documentation

int ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::AddAssociations ( const art::Ptr< recob::PFParticle > &  pfpPtr,
art::Event &  Event,
reco::shower::ShowerElementHolder &  ShowerEleHolder 
)
overrideprivatevirtual

Reimplemented from ShowerRecoTools::IShowerTool.

Definition at line 209 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

213  {
214 
215  //Check the track has been set
216  if (!ShowerEleHolder.CheckElement(fInitialTrackOutputLabel)) {
217  if (fVerbose)
218  mf::LogError("ShowerPandoraSlidingFitTrackFinderAddAssn")
219  << "Track not set so the assocation can not be made " << std::endl;
220  return 1;
221  }
222 
223  //Get the size of the ptr as it is.
224  int trackptrsize = GetVectorPtrSize(fInitialTrackOutputLabel);
225 
226  const art::Ptr<recob::Track> trackptr = GetProducedElementPtr<recob::Track>(
227  fInitialTrackOutputLabel, ShowerEleHolder, trackptrsize - 1);
228  const art::Ptr<recob::Shower> showerptr =
229  GetProducedElementPtr<recob::Shower>("shower", ShowerEleHolder);
230 
231  AddSingle<art::Assns<recob::Shower, recob::Track>>(showerptr, trackptr, "ShowerTrackAssn");
232 
233  std::vector<art::Ptr<recob::Hit>> TrackHits;
234  ShowerEleHolder.GetElement(fInitialTrackHitsInputLabel, TrackHits);
235 
236  for (auto const& TrackHit : TrackHits) {
237  AddSingle<art::Assns<recob::Track, recob::Hit>>(trackptr, TrackHit, "ShowerTrackHitAssn");
238  }
239 
240  return 0;
241  }
ShowerRecoTools::IShowerTool::GetVectorPtrSize
int GetVectorPtrSize(std::string Name)
Definition: IShowerTool.h:167
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:211
reco::shower::ShowerElementHolder::CheckElement
bool CheckElement(const std::string &Name) const
Definition: ShowerElementHolder.hh:415
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030
int ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle,
art::Event &  Event,
reco::shower::ShowerElementHolder &  ShowerEleHolder 
)
overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 79 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

83  {
84  //This is all based on the shower vertex being known. If it is not lets not do the track
85  if (!ShowerEleHolder.CheckElement(fShowerStartPositionInputLabel)) {
86  if (fVerbose)
87  mf::LogError("ShowerPandoraSlidingFitTrackFinder")
88  << "Start position not set, returning " << std::endl;
89  return 1;
90  }
91  if (!ShowerEleHolder.CheckElement(fShowerDirectionInputLabel)) {
92  if (fVerbose)
93  mf::LogError("ShowerPandoraSlidingFitTrackFinder")
94  << "Direction not set, returning " << std::endl;
95  return 1;
96  }
97  if (!ShowerEleHolder.CheckElement(fInitialTrackSpacePointsInputLabel)) {
98  if (fVerbose)
99  mf::LogError("ShowerPandoraSlidingFitTrackFinder")
100  << "Initial Spacepoints not set, returning " << std::endl;
101  return 1;
102  }
103 
104  TVector3 ShowerStartPosition = {-999, -999, -999};
105  ShowerEleHolder.GetElement(fShowerStartPositionInputLabel, ShowerStartPosition);
106 
107  TVector3 ShowerDirection = {-999, -999, -999};
108  ShowerEleHolder.GetElement(fShowerDirectionInputLabel, ShowerDirection);
109 
110  std::vector<art::Ptr<recob::SpacePoint>> spacepoints;
111  ShowerEleHolder.GetElement(fInitialTrackSpacePointsInputLabel, spacepoints);
112 
113  // The track fitter tries to create a traj point from each spacepoint so if we don't have enough
114  // spacepoints we will not get enough traj points, so let's not even try
115  if (spacepoints.size() < fMinTrajectoryPoints) {
116  if (fVerbose)
117  mf::LogWarning("ShowerPandoraSlidingFitTrackFinder")
118  << "Insufficient space points points to build track: " << spacepoints.size();
119  return 1;
120  }
121  lar_pandora::LArPandoraDetectorType* detType(lar_pandora::detector_functions::GetDetectorType());
122  // 'wirePitchW` is here used only to provide length scale for binning hits and performing sliding/local linear fits.
123  const float wirePitchW(detType->WirePitchW());
124 
125  const pandora::CartesianVector vertexPosition(
126  ShowerStartPosition.X(), ShowerStartPosition.Y(), ShowerStartPosition.Z());
127 
128  pandora::CartesianPointVector cartesianPointVector;
129  for (const art::Ptr<recob::SpacePoint> spacePoint : spacepoints)
130  cartesianPointVector.emplace_back(
131  pandora::CartesianVector(spacePoint->XYZ()[0], spacePoint->XYZ()[1], spacePoint->XYZ()[2]));
132 
133  lar_content::LArTrackStateVector trackStateVector;
134  pandora::IntVector indexVector;
135  try {
136  lar_content::LArPfoHelper::GetSlidingFitTrajectory(cartesianPointVector,
137  vertexPosition,
138  fSlidingFitHalfWindow,
139  wirePitchW,
140  trackStateVector,
141  &indexVector);
142  }
143  catch (const pandora::StatusCodeException&) {
144  if (fVerbose)
145  mf::LogWarning("ShowerPandoraSlidingFitTrackFinder")
146  << "Unable to extract sliding fit trajectory" << std::endl;
147  return 1;
148  }
149  if (trackStateVector.size() < fMinTrajectoryPoints) {
150  if (fVerbose)
151  mf::LogWarning("ShowerPandoraSlidingFitTrackFinder")
152  << "Insufficient input trajectory points to build track: " << trackStateVector.size();
153  return 1;
154  }
155 
156  if (trackStateVector.empty())
157  throw cet::exception("ShowerPandoraSlidingFitTrackFinder")
158  << "BuildTrack - No input trajectory points provided " << std::endl;
159 
160  recob::tracking::Positions_t xyz;
161  recob::tracking::Momenta_t pxpypz;
162  recob::TrackTrajectory::Flags_t flags;
163 
164  for (const lar_content::LArTrackState& trackState : trackStateVector) {
165  xyz.emplace_back(recob::tracking::Point_t(trackState.GetPosition().GetX(),
166  trackState.GetPosition().GetY(),
167  trackState.GetPosition().GetZ()));
168  pxpypz.emplace_back(recob::tracking::Vector_t(trackState.GetDirection().GetX(),
169  trackState.GetDirection().GetY(),
170  trackState.GetDirection().GetZ()));
171 
172  // Set flag NoPoint if point has bogus coordinates, otherwise use clean flag set
173  if (std::fabs(trackState.GetPosition().GetX() - util::kBogusF) <
174  std::numeric_limits<float>::epsilon() &&
175  std::fabs(trackState.GetPosition().GetY() - util::kBogusF) <
176  std::numeric_limits<float>::epsilon() &&
177  std::fabs(trackState.GetPosition().GetZ() - util::kBogusF) <
178  std::numeric_limits<float>::epsilon()) {
179  flags.emplace_back(recob::TrajectoryPointFlags(recob::TrajectoryPointFlags::InvalidHitIndex,
180  recob::TrajectoryPointFlagTraits::NoPoint));
181  }
182  else {
183  flags.emplace_back(recob::TrajectoryPointFlags());
184  }
185  }
186 
187  // note from gc: eventually we should produce a TrackTrajectory, not a Track with empty covariance matrix and bogus chi2, etc.
188  recob::Track InitialTrack = recob::Track(
189  recob::TrackTrajectory(std::move(xyz), std::move(pxpypz), std::move(flags), false),
190  util::kBogusI,
191  util::kBogusF,
192  util::kBogusI,
193  recob::tracking::SMatrixSym55(),
194  recob::tracking::SMatrixSym55(),
195  pfparticle.key());
196 
197  ShowerEleHolder.SetElement(InitialTrack, fInitialTrackOutputLabel);
198 
199  TVector3 Start = {InitialTrack.Start().X(), InitialTrack.Start().Y(), InitialTrack.Start().Z()};
200  TVector3 End = {InitialTrack.End().X(), InitialTrack.End().Y(), InitialTrack.End().Z()};
201  float tracklength = (Start - End).Mag();
202 
203  ShowerEleHolder.SetElement(tracklength, fInitialTrackLengthOutputLabel);
204 
205  return 0;
206  }
recob::TrajectoryPointFlagTraits::NoPoint
static constexpr Flag_t NoPoint
The trajectory point is not defined.
Definition: TrajectoryPointFlags.h:119
recob::tracking::SMatrixSym55
ROOT::Math::SMatrix< Double32_t, 5, 5, ROOT::Math::MatRepSym< Double32_t, 5 > > SMatrixSym55
Definition: TrackingTypes.h:85
lar_content::LArTrackState
LArTrackState class.
Definition: LArPfoObjects.h:29
lar_pandora::LArPandoraDetectorType
Empty interface to map pandora to specifics in the LArSoft geometry.
Definition: LArPandoraDetectorType.h:25
reco::shower::ShowerElementHolder::SetElement
void SetElement(T &dataproduct, const std::string &Name, bool checktag=false)
Definition: ShowerElementHolder.hh:360
lar_content::LArPfoHelper::GetSlidingFitTrajectory
static void GetSlidingFitTrajectory(const pandora::CartesianPointVector &pointVector, const pandora::CartesianVector &vertexPosition, const unsigned int layerWindow, const float layerPitch, LArTrackStateVector &trackStateVector, pandora::IntVector *const pIndexVector=nullptr)
Apply 3D sliding fit to a set of 3D points and return track trajectory.
util::kBogusI
constexpr int kBogusI
obviously bogus integer value
Definition: PhysicalConstants.h:72
pandora::IntVector
std::vector< int > IntVector
Definition: LArPandoraHelper.h:24
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:211
recob::tracking::Vector_t
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< Coord_t >, ROOT::Math::GlobalCoordinateSystemTag > Vector_t
Type for representation of momenta in 3D space. See recob::tracking::Coord_t for more details on the ...
Definition: TrackingTypes.h:29
recob::Track::Start
Point_t const & Start() const
Access to track position at different points.
Definition: Track.h:123
recob::TrackTrajectory
A trajectory in space reconstructed from hits.
Definition: TrackTrajectory.h:62
wirecell.gen.depos.move
def move(depos, offset)
Definition: depos.py:107
art::Ptr::key
key_type key() const noexcept
Definition: Ptr.h:216
reco::shower::ShowerElementHolder::CheckElement
bool CheckElement(const std::string &Name) const
Definition: ShowerElementHolder.hh:415
recob::tracking::Momenta_t
std::vector< Vector_t > Momenta_t
Type of momentum list.
Definition: TrackingTypes.h:35
recob::TrackTrajectory::Flags_t
std::vector< PointFlags_t > Flags_t
Type of point flag list.
Definition: TrackTrajectory.h:88
reco::shower::ShowerElementHolder::GetElement
int GetElement(const std::string &Name, T &Element) const
Definition: ShowerElementHolder.hh:219
recob::TrajectoryPointFlags::InvalidHitIndex
static constexpr HitIndex_t InvalidHitIndex
Value marking an invalid hit index.
Definition: TrajectoryPointFlags.h:336
recob::tracking::Positions_t
std::vector< Point_t > Positions_t
Type of trajectory point list.
Definition: TrackingTypes.h:32
lar_pandora::detector_functions::GetDetectorType
LArPandoraDetectorType * GetDetectorType()
Factory class that returns the correct detector type interface.
Definition: LArPandoraDetectorType.cxx:22
util::kBogusF
constexpr float kBogusF
obviously bogus float value
Definition: PhysicalConstants.h:73
End
void End(void)
Definition: gXSecComp.cxx:210
mf::LogWarning
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
Definition: MessageLogger.h:210
recob::Track::End
Point_t const & End() const
Definition: Track.h:125
proxy::Track
TrackCollectionProxyElement< TrackCollProxy > Track
Proxy to an element of a proxy collection of recob::Track objects.
Definition: Track.h:1036
lar_content::LArTrackStateVector
std::vector< LArTrackState > LArTrackStateVector
Definition: LArPfoObjects.h:67
recob::tracking::Point_t
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< Coord_t >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space. See recob::tracking::Coord_t for more detai...
Definition: TrackingTypes.h:26
recob::TrajectoryPointFlags
Set of flags pertaining a point of the track.
Definition: TrajectoryPointFlags.h:315
Start
#define Start
Definition: config.cpp:1229
recob::Track
Track from a non-cascading particle.A recob::Track consists of a recob::TrackTrajectory, plus additional members relevant for a "fitted" track:
Definition: Track.h:49
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030
void ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::InitialiseProducers ( )
overrideprivatevirtual

Reimplemented from ShowerRecoTools::IShowerTool.

Definition at line 69 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

70  {
71 
72  InitialiseProduct<std::vector<recob::Track>>(fInitialTrackOutputLabel);
73  InitialiseProduct<art::Assns<recob::Shower, recob::Track>>("ShowerTrackAssn");
74  InitialiseProduct<art::Assns<recob::Track, recob::Hit>>("ShowerTrackHitAssn");
75  }

Member Data Documentation

art::ServiceHandle<geo::Geometry> ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fGeom
private

Definition at line 40 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

std::string ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fInitialTrackHitsInputLabel
private

Definition at line 51 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

std::string ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fInitialTrackLengthOutputLabel
private

Definition at line 47 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

std::string ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fInitialTrackOutputLabel
private

Definition at line 46 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

std::string ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fInitialTrackSpacePointsInputLabel
private

Definition at line 50 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

float ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fMinTrajectoryPoints
private

Definition at line 45 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

std::string ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fShowerDirectionInputLabel
private

Definition at line 49 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

std::string ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fShowerStartPositionInputLabel
private

Definition at line 48 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

float ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fSlidingFitHalfWindow
private

Definition at line 44 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

int ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::fVerbose
private

Definition at line 43 of file ShowerPandoraSlidingFitTrackFinder_tool.cc.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.11