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

Public Member Functions

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

Private Member Functions

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)
 
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

int fVerbose
 
bool fUsePandoraVertex
 
art::InputTag fHitModuleLabel
 
art::InputTag fPFParticleLabel
 
std::string fInitialTrackHitsInputLabel
 
std::string fShowerStartPositionInputLabel
 
std::string fInitialTrackInputLabel
 
std::string fShowerDirectionOutputLabel
 

Detailed Description

Definition at line 21 of file ShowerTrackHitDirection_tool.cc.

Constructor & Destructor Documentation

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

Definition at line 45 of file ShowerTrackHitDirection_tool.cc.

46  : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
47  , fVerbose(pset.get<int>("Verbose"))
48  , fUsePandoraVertex(pset.get<bool>("UsePandoraVertex"))
49  , fHitModuleLabel(pset.get<art::InputTag>("HitModuleLabel"))
50  , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
51  , fInitialTrackHitsInputLabel(pset.get<std::string>("InitialTrackHitsInputLabel"))
52  , fShowerStartPositionInputLabel(pset.get<std::string>("ShowerStartPositionInputLabel"))
53  , fInitialTrackInputLabel(pset.get<std::string>("InitialTrackInputLabel"))
54  , fShowerDirectionOutputLabel(pset.get<std::string>("ShowerDirectionOutputLabel"))
55  {}
T get(std::string const &key) const
Definition: ParameterSet.h:314
IShowerTool(const fhicl::ParameterSet &pset)
Definition: IShowerTool.h:31

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 150 of file IShowerTool.h.

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

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

Implements ShowerRecoTools::IShowerTool.

Definition at line 57 of file ShowerTrackHitDirection_tool.cc.

References util::abs(), reco::shower::ShowerElementHolder::CheckElement(), DEFINE_ART_CLASS_TOOL, fHitModuleLabel, fInitialTrackHitsInputLabel, fInitialTrackInputLabel, fPFParticleLabel, fShowerDirectionOutputLabel, fShowerStartPositionInputLabel, fUsePandoraVertex, fVerbose, reco::shower::ShowerElementHolder::GetElement(), reco::shower::ShowerElementHolder::GetFindManyP(), art::ProductRetriever::getValidHandle(), reco::shower::ShowerElementHolder::SetElement(), recob::Track::Start(), X, Y, and Z.

60  {
61 
62  //Check the Track Hits has been defined
63  if (!ShowerEleHolder.CheckElement(fInitialTrackHitsInputLabel)) {
64  if (fVerbose)
65  mf::LogError("ShowerTrackHitDirection") << "Initial track hits not set" << std::endl;
66  return 0;
67  }
68 
69  //Check the start position is set.
71  if (fVerbose)
72  mf::LogError("ShowerTrackHitDirection")
73  << "Start position not set, returning " << std::endl;
74  return 0;
75  }
76 
77  //Get the start poistion
78  geo::Point_t StartPosition = {-999, -999, -99};
79  if (fUsePandoraVertex) {
80  ShowerEleHolder.GetElement(fShowerStartPositionInputLabel, StartPosition);
81  }
82  else {
83  //Check the Tracks has been defined
84  if (!ShowerEleHolder.CheckElement(fInitialTrackInputLabel)) {
85  if (fVerbose)
86  mf::LogError("ShowerTrackHitDirection") << "Initial track not set" << std::endl;
87  return 0;
88  }
89  recob::Track InitialTrack;
90  ShowerEleHolder.GetElement(fInitialTrackInputLabel, InitialTrack);
91  geo::Point_t Start_point = InitialTrack.Start();
92  StartPosition = {Start_point.X(), Start_point.Y(), Start_point.Z()};
93  }
94 
95  //Get the spacepoints associated to hits
96  auto const hitHandle = Event.getValidHandle<std::vector<recob::Hit>>(fHitModuleLabel);
97 
98  //Get the spacepoint handle. We need to do this in 3D.
100  ShowerEleHolder.GetFindManyP<recob::SpacePoint>(hitHandle, Event, fPFParticleLabel);
101 
102  //Get the initial track hits.
103  std::vector<art::Ptr<recob::Hit>> InitialTrackHits;
104  ShowerEleHolder.GetElement(fInitialTrackHitsInputLabel, InitialTrackHits);
105 
106  //Calculate the mean direction and the the standard deviation
107  float sumX = 0, sumX2 = 0;
108  float sumY = 0, sumY2 = 0;
109  float sumZ = 0, sumZ2 = 0;
110 
111  //Get the spacepoints associated to the track hit
112  std::vector<art::Ptr<recob::SpacePoint>> intitaltrack_sp;
113  for (auto const hit : InitialTrackHits) {
114  std::vector<art::Ptr<recob::SpacePoint>> sps = fmsp.at(hit.key());
115  for (auto const sp : sps) {
116  intitaltrack_sp.push_back(sp);
117 
118  //Get the direction relative to the start positon
119  auto const pos = sp->position() - StartPosition;
120  if (pos.R() == 0) { continue; }
121 
122  sumX = pos.X();
123  sumX2 += pos.X() * pos.X();
124  sumY = pos.Y();
125  sumY2 += pos.Y() * pos.Y();
126  sumZ = pos.Z();
127  sumZ2 += pos.Z() * pos.Z();
128  }
129  }
130 
131  float NumSps = intitaltrack_sp.size();
132  auto const Mean = geo::Vector_t{sumX / NumSps, sumY / NumSps, sumZ / NumSps}.Unit();
133 
134  float RMSX = 999;
135  float RMSY = 999;
136  float RMSZ = 999;
137  if (sumX2 / NumSps - ((sumX / NumSps) * ((sumX / NumSps))) > 0) {
138  RMSX = std::sqrt(sumX2 / NumSps - ((sumX / NumSps) * ((sumX / NumSps))));
139  }
140  if (sumY2 / NumSps - ((sumY / NumSps) * ((sumY / NumSps))) > 0) {
141  RMSY = std::sqrt(sumY2 / NumSps - ((sumY / NumSps) * ((sumY / NumSps))));
142  }
143  if (sumZ2 / NumSps - ((sumZ / NumSps) * ((sumZ / NumSps))) > 0) {
144  RMSZ = std::sqrt(sumZ2 / NumSps - ((sumZ / NumSps) * ((sumZ / NumSps))));
145  }
146 
147  //Loop over the spacepoints and remove ones the relative direction is not within one sigma.
148  geo::Vector_t Direction_Mean{0, 0, 0};
149  int N = 0;
150  for (auto const sp : intitaltrack_sp) {
151  auto const Direction = sp->position() - StartPosition;
152  if ((std::abs((Direction - Mean).X()) < 1 * RMSX) &&
153  (std::abs((Direction - Mean).Y()) < 1 * RMSY) &&
154  (std::abs((Direction - Mean).Z()) < 1 * RMSZ)) {
155  if (Direction.R() == 0) { continue; }
156  ++N;
157  Direction_Mean += Direction;
158  }
159  }
160 
161  if (N > 0) {
162  //Take the mean value
163  Direction_Mean = Direction_Mean.Unit();
164  ShowerEleHolder.SetElement(Direction_Mean, fShowerDirectionOutputLabel);
165  }
166  else {
167  if (fVerbose)
168  mf::LogError("ShowerTrackHitDirection")
169  << "None of the points are within 1 sigma" << std::endl;
170  return 1;
171  }
172  return 0;
173  }
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Vector_t
Type for representation of momenta in 3D space.
Definition: geo_vectors.h:160
Float_t Y
Definition: plot.C:37
void SetElement(T &dataproduct, const std::string &Name, bool checktag=false)
constexpr auto abs(T v)
Returns the absolute value of the argument.
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Float_t Z
Definition: plot.C:37
Point_t const & Start() const
Access to track position at different points.
Definition: Track.h:157
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
Direction
Definition: types.h:12
Float_t X
Definition: plot.C:37
const art::FindManyP< T1 > & GetFindManyP(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 80 of file IShowerTool.h.

Referenced by ShowerRecoTools::ShowerLengthPercentile::CalculateElement(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPCAPropergationStartPosition::CalculateElement(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPFPVertexStartPosition::CalculateElement(), ShowerRecoTools::ShowerUnidirectiondEdx::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().

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

Definition at line 100 of file IShowerTool.h.

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

103  {
104 
105  //Check the element has been set
106  bool check_element = ShowerEleHolder.CheckElement(Name);
107  if (!check_element) {
108  throw cet::exception("IShowerTool") << "tried to get a element that does not exist. Failed "
109  "at making the art ptr for Element: "
110  << Name << std::endl;
111  }
112 
113  //Check the unique ptr has been set.
114  bool check_ptr = UniquePtrs->CheckUniqueProduerPtr(Name);
115  if (!check_ptr) {
116  throw cet::exception("IShowerTool")
117  << "tried to get a ptr that does not exist. Failed at making the art ptr for Element"
118  << Name;
119  }
120 
121  //Check if the user has defined an index if not just use the current shower index/
122  int index;
123  if (iter != -1) { index = iter; }
124  else {
125  index = ShowerEleHolder.GetShowerNumber();
126  }
127 
128  //Make the ptr
129  return UniquePtrs->GetArtPtr<T>(Name, index);
130  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
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 66 of file IShowerTool.h.

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

Definition at line 135 of file IShowerTool.h.

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

136  {
137 
138  if (collectorPtr == nullptr) {
139  mf::LogWarning("IShowerTool") << "The art::ProducesCollector ptr has not been set";
140  return;
141  }
142 
143  collectorPtr->produces<T>(InstanceName);
144  UniquePtrs->SetShowerUniqueProduerPtr(type<T>(), Name, InstanceName);
145  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
art::ProducesCollector * collectorPtr
Definition: IShowerTool.h:95
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 160 of file IShowerTool.h.

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

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

160 { UniquePtrs->PrintPtr(Name); }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
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 44 of file IShowerTool.h.

48  {
49 
50  int calculation_status = CalculateElement(pfparticle, Event, ShowerEleHolder);
51  if (calculation_status != 0) return calculation_status;
52  if (fRunEventDisplay) {
54  pfparticle, Event, ShowerEleHolder, evd_display_name_append);
55  }
56  return calculation_status;
57  }
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:80
void ShowerRecoTools::IShowerTool::SetPtr ( art::ProducesCollector collector)
inlineinherited

Definition at line 63 of file IShowerTool.h.

63 { collectorPtr = collector; }
art::ProducesCollector * collectorPtr
Definition: IShowerTool.h:95

Member Data Documentation

art::InputTag ShowerRecoTools::ShowerTrackHitDirection::fHitModuleLabel
private

Definition at line 36 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackHitDirection::fInitialTrackHitsInputLabel
private

Definition at line 39 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackHitDirection::fInitialTrackInputLabel
private

Definition at line 41 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

art::InputTag ShowerRecoTools::ShowerTrackHitDirection::fPFParticleLabel
private

Definition at line 37 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackHitDirection::fShowerDirectionOutputLabel
private

Definition at line 42 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerTrackHitDirection::fShowerStartPositionInputLabel
private

Definition at line 40 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

bool ShowerRecoTools::ShowerTrackHitDirection::fUsePandoraVertex
private

Definition at line 34 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().

int ShowerRecoTools::ShowerTrackHitDirection::fVerbose
private

Definition at line 33 of file ShowerTrackHitDirection_tool.cc.

Referenced by CalculateElement().


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