LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
ShowerRecoTools::ShowerTrackTrajToSpacePoint Class Reference
Inheritance diagram for ShowerRecoTools::ShowerTrackTrajToSpacePoint:
ShowerRecoTools::IShowerTool

Public Member Functions

 ShowerTrackTrajToSpacePoint (const fhicl::ParameterSet &pset)
 
int CalculateElement (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder) override
 
int RunShowerTool (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder, std::string evd_display_name_append="")
 
virtual void InitialiseProducers ()
 
void SetPtr (art::ProducesCollector *collector)
 
void InitaliseProducerPtr (reco::shower::ShowerProducedPtrsHolder &uniqueproducerPtrs)
 
virtual int AddAssociations (const art::Ptr< recob::PFParticle > &pfpPtr, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder)
 

Protected Member Functions

const shower::LArPandoraShowerAlgGetLArPandoraShowerAlg () 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

float fMaxDist
 
art::InputTag fPFParticleLabel
 
int fVerbose
 
std::string fInitialTrackSpacePointsOutputLabel
 
std::string fInitialTrackHitsOutputLabel
 
std::string fInitialTrackInputTag
 
std::string fShowerStartPositionInputTag
 
std::string fInitialTrackSpacePointsInputTag
 

Detailed Description

Definition at line 20 of file ShowerTrackTrajToSpacePoint_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerTrackTrajToSpacePoint::ShowerTrackTrajToSpacePoint ( const fhicl::ParameterSet pset)

Definition at line 43 of file ShowerTrackTrajToSpacePoint_tool.cc.

44  : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
45  , fMaxDist(pset.get<float>("MaxDist"))
46  , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
47  , fVerbose(pset.get<int>("Verbose"))
49  pset.get<std::string>("InitialTrackSpacePointsOutputLabel"))
50  , fInitialTrackHitsOutputLabel(pset.get<std::string>("InitialTrackHitsOutputLabel"))
51  , fInitialTrackInputTag(pset.get<std::string>("InitialTrackInputTag"))
52  , fShowerStartPositionInputTag(pset.get<std::string>("ShowerStartPositionInputTag"))
53  , fInitialTrackSpacePointsInputTag(pset.get<std::string>("InitialTrackSpacePointsInputTag"))
54  {}
T get(std::string const &key) const
Definition: ParameterSet.h:314
IShowerTool(const fhicl::ParameterSet &pset)
Definition: IShowerTool.h:33

Member Function Documentation

virtual int ShowerRecoTools::IShowerTool::AddAssociations ( const art::Ptr< recob::PFParticle > &  pfpPtr,
art::Event Event,
reco::shower::ShowerElementHolder ShowerEleHolder 
)
inlinevirtualinherited
template<class T , class A , class B >
void ShowerRecoTools::IShowerTool::AddSingle ( A &  a,
B &  b,
std::string  Name 
)
inlineprotectedinherited

Definition at line 152 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::AddSingle().

153  {
154  UniquePtrs->AddSingle<T>(a, b, Name);
155  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:85
void AddSingle(A &a, B &b, const std::string &Name)
int ShowerRecoTools::ShowerTrackTrajToSpacePoint::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle,
art::Event Event,
reco::shower::ShowerElementHolder ShowerEleHolder 
)
overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 56 of file ShowerTrackTrajToSpacePoint_tool.cc.

References reco::shower::ShowerElementHolder::CheckElement(), DEFINE_ART_CLASS_TOOL, larg4::dist(), fInitialTrackHitsOutputLabel, fInitialTrackInputTag, fInitialTrackSpacePointsInputTag, fInitialTrackSpacePointsOutputLabel, recob::Track::FlagsAtPoint(), fMaxDist, fPFParticleLabel, fShowerStartPositionInputTag, fVerbose, reco::shower::ShowerElementHolder::GetElement(), reco::shower::ShowerElementHolder::GetFindOneP(), art::ProductRetriever::getValidHandle(), recob::Track::LocationAtPoint(), recob::TrajectoryPointFlagTraits::NoPoint, recob::Track::NumberTrajectoryPoints(), recob::SpacePoint::position(), and reco::shower::ShowerElementHolder::SetElement().

60  {
61 
62  //Check the Track has been defined
63  if (!ShowerEleHolder.CheckElement(fInitialTrackInputTag)) {
64  if (fVerbose)
65  mf::LogError("ShowerTrackTrajToSpacePoint") << "Initial track not set" << std::endl;
66  return 0;
67  }
68 
69  //Check the start position is set.
70  if (!ShowerEleHolder.CheckElement(fShowerStartPositionInputTag)) {
71  if (fVerbose)
72  mf::LogError("ShowerTrackTrajToSpacePoint")
73  << "Start position not set, returning " << std::endl;
74  return 0;
75  }
76 
77  //Check the Track Hits has been defined
78  if (!ShowerEleHolder.CheckElement(fInitialTrackSpacePointsInputTag)) {
79  if (fVerbose)
80  mf::LogError("ShowerTrackTrajToSpacePoint")
81  << "Initial track spacepoints not set" << std::endl;
82  return 0;
83  }
84 
85  //Get the start poistion
86  geo::Point_t ShowerStartPosition = {-999, -999, -999};
87  ShowerEleHolder.GetElement(fShowerStartPositionInputTag, ShowerStartPosition);
88 
89  //Get the initial track hits.
90  std::vector<art::Ptr<recob::SpacePoint>> intitaltrack_sp;
91  ShowerEleHolder.GetElement(fInitialTrackSpacePointsInputTag, intitaltrack_sp);
92 
93  //Get the track
94  recob::Track InitialTrack;
95  ShowerEleHolder.GetElement(fInitialTrackInputTag, InitialTrack);
96 
97  std::vector<art::Ptr<recob::SpacePoint>> new_intitaltrack_sp;
98  //Loop over the trajectory points
99  for (unsigned int traj = 0; traj < InitialTrack.NumberTrajectoryPoints(); ++traj) {
100 
101  //ignore bogus info.
102  auto flags = InitialTrack.FlagsAtPoint(traj);
103  if (flags.isSet(recob::TrajectoryPointFlagTraits::NoPoint)) { continue; }
104 
105  geo::Point_t TrajPosition = InitialTrack.LocationAtPoint(traj);
106  geo::Point_t TrajPositionStart = InitialTrack.LocationAtPoint(0);
107 
108  //Ignore values with 0 mag from the start position
109  if ((TrajPosition - TrajPositionStart).R() == 0) { continue; }
110  if ((TrajPosition - ShowerStartPosition).R() == 0) { continue; }
111 
112  float MinDist = 9999;
113  unsigned int index = 999;
114  for (unsigned int sp = 0; sp < intitaltrack_sp.size(); ++sp) {
115  //Find the spacepoint closest to the trajectory point.
116  art::Ptr<recob::SpacePoint> spacepoint = intitaltrack_sp[sp];
117  auto const dist = (spacepoint->position() - TrajPosition).R();
118  if (dist < MinDist && dist < fMaxDist) {
119  MinDist = dist;
120  index = sp;
121  }
122  }
123 
124  if (index == 999) { continue; }
125  //Add the spacepoint to the track spacepoints.
126  new_intitaltrack_sp.push_back(intitaltrack_sp[index]);
127 
128  //Delete the spacepoint so it can not be used again.
129  intitaltrack_sp.erase(intitaltrack_sp.begin() + index);
130  }
131 
132  // Get the spacepoints
133  auto const spHandle = Event.getValidHandle<std::vector<recob::SpacePoint>>(fPFParticleLabel);
134 
135  // Get the hits associated with the space points
136  const art::FindOneP<recob::Hit>& fohsp =
137  ShowerEleHolder.GetFindOneP<recob::Hit>(spHandle, Event, fPFParticleLabel);
138 
139  //Save the corresponding hits
140  std::vector<art::Ptr<recob::Hit>> trackHits;
141  for (auto const& spacePoint : new_intitaltrack_sp) {
142  //Get the hits
143  const art::Ptr<recob::Hit> hit = fohsp.at(spacePoint.key());
144  trackHits.push_back(hit);
145  }
146 
147  //Save the spacepoints.
148  ShowerEleHolder.SetElement(new_intitaltrack_sp, fInitialTrackSpacePointsOutputLabel);
149  ShowerEleHolder.SetElement(trackHits, fInitialTrackHitsOutputLabel);
150 
151  return 0;
152  }
static constexpr Flag_t NoPoint
The trajectory point is not defined.
Point_t const & LocationAtPoint(size_t i) const
Access to track position at different points.
Definition: Track.h:160
void SetElement(T &dataproduct, const std::string &Name, bool checktag=false)
geo::Point_t position() const
Returns the position of the point in world coordinates [cm].
Definition: SpacePoint.h:91
size_t NumberTrajectoryPoints() const
Various functions related to the presence and the number of (valid) points.
Definition: Track.h:136
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
bool CheckElement(const std::string &Name) const
int GetElement(const std::string &Name, T &Element) const
Detector simulation of raw signals on wires.
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space.
Definition: geo_vectors.h:180
ValidHandle< PROD > getValidHandle(InputTag const &tag) const
constexpr double dist(const TReal *x, const TReal *y, const unsigned int dimension)
PointFlags_t const & FlagsAtPoint(size_t i) const
Access to i-th TrajectoryPoint or its Flags.
Definition: Track.h:152
2D representation of charge deposited in the TDC/wire plane
Definition: Hit.h:46
const art::FindOneP< T1 > & GetFindOneP(const art::ValidHandle< std::vector< T2 >> &handle, const art::Event &evt, const art::InputTag &moduleTag)
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
const shower::LArPandoraShowerAlg& ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg ( ) const
inlineprotectedinherited

Definition at line 82 of file IShowerTool.h.

Referenced by ShowerRecoTools::ShowerLengthPercentile::CalculateElement(), ShowerRecoTools::ShowerUnidirectiondEdx::CalculateElement(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPCAPropergationStartPosition::CalculateElement(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPFPVertexStartPosition::CalculateElement(), ShowerRecoTools::ShowerDirectionCheater::CalculateElement(), ShowerRecoTools::ShowerPCADirection::CalculateElement(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerTrajPointdEdx::CalculateElement(), ShowerRecoTools::ShowerPCADirection::CalculateShowerPCA(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::FindInitialTrackHits(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::FindTrackSpacePoints(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::FitSegmentAndCalculateResidual(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::PruneFrontOfSPSPool(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::PruneTrack(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::RunTestOfIncrementalSpacePointFinder(), ShowerRecoTools::ShowerTrackPCADirection::ShowerPCAVector(), and ShowerRecoTools::ShowerIncrementalTrackHitFinder::ShowerPCAVector().

83  {
84  return fLArPandoraShowerAlg;
85  };
shower::LArPandoraShowerAlg fLArPandoraShowerAlg
Definition: IShowerTool.h:92
template<class T >
art::Ptr<T> ShowerRecoTools::IShowerTool::GetProducedElementPtr ( std::string  Name,
reco::shower::ShowerElementHolder ShowerEleHolder,
int  iter = -1 
)
inlineprotectedinherited

Definition at line 102 of file IShowerTool.h.

References reco::shower::ShowerElementHolder::CheckElement(), reco::shower::ShowerProducedPtrsHolder::CheckUniqueProduerPtr(), reco::shower::ShowerProducedPtrsHolder::GetArtPtr(), and reco::shower::ShowerElementHolder::GetShowerNumber().

105  {
106 
107  //Check the element has been set
108  bool check_element = ShowerEleHolder.CheckElement(Name);
109  if (!check_element) {
110  throw cet::exception("IShowerTool") << "tried to get a element that does not exist. Failed "
111  "at making the art ptr for Element: "
112  << Name << std::endl;
113  }
114 
115  //Check the unique ptr has been set.
116  bool check_ptr = UniquePtrs->CheckUniqueProduerPtr(Name);
117  if (!check_ptr) {
118  throw cet::exception("IShowerTool")
119  << "tried to get a ptr that does not exist. Failed at making the art ptr for Element"
120  << Name;
121  }
122 
123  //Check if the user has defined an index if not just use the current shower index/
124  int index;
125  if (iter != -1) { index = iter; }
126  else {
127  index = ShowerEleHolder.GetShowerNumber();
128  }
129 
130  //Make the ptr
131  return UniquePtrs->GetArtPtr<T>(Name, index);
132  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:85
bool CheckUniqueProduerPtr(const std::string &Name) const
bool CheckElement(const std::string &Name) const
art::Ptr< T > GetArtPtr(const std::string &Name, const int &iter) const
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
int ShowerRecoTools::IShowerTool::GetVectorPtrSize ( std::string  Name)
inlineprotectedinherited
void ShowerRecoTools::IShowerTool::InitaliseProducerPtr ( reco::shower::ShowerProducedPtrsHolder uniqueproducerPtrs)
inlineinherited

Definition at line 68 of file IShowerTool.h.

69  {
70  UniquePtrs = &uniqueproducerPtrs;
71  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:85
virtual void ShowerRecoTools::IShowerTool::InitialiseProducers ( )
inlinevirtualinherited
template<class T >
void ShowerRecoTools::IShowerTool::InitialiseProduct ( std::string  Name,
std::string  InstanceName = "" 
)
inlineprotectedinherited

Definition at line 137 of file IShowerTool.h.

References art::ProducesCollector::produces(), and reco::shower::ShowerProducedPtrsHolder::SetShowerUniqueProduerPtr().

138  {
139 
140  if (collectorPtr == nullptr) {
141  mf::LogWarning("IShowerTool") << "The art::ProducesCollector ptr has not been set";
142  return;
143  }
144 
145  collectorPtr->produces<T>(InstanceName);
146  UniquePtrs->SetShowerUniqueProduerPtr(type<T>(), Name, InstanceName);
147  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:85
art::ProducesCollector * collectorPtr
Definition: IShowerTool.h:97
void produces(std::string const &instanceName={}, Persistable const persistable=Persistable::Yes)
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
int SetShowerUniqueProduerPtr(type< T >, const std::string &Name, const std::string &Instance="")
void ShowerRecoTools::IShowerTool::PrintPtr ( std::string  Name)
inlineprotectedinherited

Definition at line 162 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::PrintPtr().

Referenced by ShowerRecoTools::ShowerExampleTool::CalculateElement().

162 { UniquePtrs->PrintPtr(Name); }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:85
void PrintPtr(const std::string &Name) const
void ShowerRecoTools::IShowerTool::PrintPtrs ( )
inlineprotectedinherited
int ShowerRecoTools::IShowerTool::RunShowerTool ( const art::Ptr< recob::PFParticle > &  pfparticle,
art::Event Event,
reco::shower::ShowerElementHolder ShowerEleHolder,
std::string  evd_display_name_append = "" 
)
inlineinherited

Definition at line 46 of file IShowerTool.h.

50  {
51 
52  int calculation_status = CalculateElement(pfparticle, Event, ShowerEleHolder);
53  if (calculation_status != 0) return calculation_status;
54  if (fRunEventDisplay) {
56  pfparticle, Event, ShowerEleHolder, evd_display_name_append);
57  }
58  return calculation_status;
59  }
void DebugEVD(art::Ptr< recob::PFParticle > const &pfparticle, art::Event const &Event, const reco::shower::ShowerElementHolder &ShowerEleHolder, std::string const &evd_disp_name_append="") const
virtual int CalculateElement(const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder)=0
const shower::LArPandoraShowerAlg & GetLArPandoraShowerAlg() const
Definition: IShowerTool.h:82
void ShowerRecoTools::IShowerTool::SetPtr ( art::ProducesCollector collector)
inlineinherited

Definition at line 65 of file IShowerTool.h.

65 { collectorPtr = collector; }
art::ProducesCollector * collectorPtr
Definition: IShowerTool.h:97

Member Data Documentation

std::string ShowerRecoTools::ShowerTrackTrajToSpacePoint::fInitialTrackHitsOutputLabel
private

Definition at line 37 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackTrajToSpacePoint::fInitialTrackInputTag
private

Definition at line 38 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackTrajToSpacePoint::fInitialTrackSpacePointsInputTag
private

Definition at line 40 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackTrajToSpacePoint::fInitialTrackSpacePointsOutputLabel
private

Definition at line 36 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

float ShowerRecoTools::ShowerTrackTrajToSpacePoint::fMaxDist
private

Definition at line 31 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

art::InputTag ShowerRecoTools::ShowerTrackTrajToSpacePoint::fPFParticleLabel
private

Definition at line 33 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackTrajToSpacePoint::fShowerStartPositionInputTag
private

Definition at line 39 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().

int ShowerRecoTools::ShowerTrackTrajToSpacePoint::fVerbose
private

Definition at line 34 of file ShowerTrackTrajToSpacePoint_tool.cc.

Referenced by CalculateElement().


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