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

Public Member Functions

 ShowerPFPVertexStartPosition (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

art::InputTag fPFParticleLabel
 
int fVerbose
 
std::string fShowerStartPositionOutputLabel
 
std::string fShowerDirectionInputLabel
 

Detailed Description

Definition at line 23 of file ShowerPFPVertexStartPosition_tool.cc.

Constructor & Destructor Documentation

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

Definition at line 41 of file ShowerPFPVertexStartPosition_tool.cc.

42  : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
43  , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
44  , fVerbose(pset.get<int>("Verbose"))
45  , fShowerStartPositionOutputLabel(pset.get<std::string>("ShowerStartPositionOutputLabel"))
46  , fShowerDirectionInputLabel(pset.get<std::string>("ShowerDirectionInputLabel"))
47  {}
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::ShowerPFPVertexStartPosition::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle,
art::Event Event,
reco::shower::ShowerElementHolder ShowerEleHolder 
)
overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 49 of file ShowerPFPVertexStartPosition_tool.cc.

References reco::shower::ShowerElementHolder::CheckElement(), DEFINE_ART_CLASS_TOOL, fPFParticleLabel, fShowerDirectionInputLabel, fShowerStartPositionOutputLabel, fVerbose, reco::shower::ShowerElementHolder::GetElement(), reco::shower::ShowerElementHolder::GetFindManyP(), ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg(), art::ProductRetriever::getValidHandle(), art::Ptr< T >::key(), shower::LArPandoraShowerAlg::OrderShowerSpacePoints(), recob::Vertex::position(), reco::shower::ShowerElementHolder::SetElement(), and shower::LArPandoraShowerAlg::ShowerCentre().

53  {
54 
55  // Get the assocated pfParicle vertex PFParticles
56  auto const pfpHandle = Event.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);
57 
59  ShowerEleHolder.GetFindManyP<recob::Vertex>(pfpHandle, Event, fPFParticleLabel);
60  // art::FindManyP<recob::Vertex> fmv(pfpHandle, Event, fPFParticleLabel);
61  if (!fmv.isValid()) {
62  throw cet::exception("ShowerPFPVertexStartPosition")
63  << "Vertex and PF particle association is somehow not valid. Stopping";
64  }
65 
66  std::vector<art::Ptr<recob::Vertex>> vtx_cand;
67  try {
68  vtx_cand = fmv.at(pfparticle.key());
69  }
70  catch (...) {
71  if (fVerbose)
72  mf::LogError("ShowerPFPVertexStartPosition") << "PFP-Vertex assan not set, returning";
73  return 1;
74  }
75  //If there is more than one then fail becuase I don't think that this can be the case
76  if (vtx_cand.size() != 1) {
77  if (fVerbose)
78  mf::LogError("ShowerPFPVertexStartPosition")
79  << "Wrong number of vertices: " << vtx_cand.size() << ", returning";
80  return 1;
81  }
82 
83  //If there is only one vertex good news we just say that is the start of the shower.
84  if (vtx_cand.size() == 1) {
85  art::Ptr<recob::Vertex> StartPositionVertex = vtx_cand[0];
86  auto ShowerStartPosition(StartPositionVertex->position());
87  geo::Point_t ShowerStartPositionErr = {-999, -999, -999};
88  ShowerEleHolder.SetElement(
89  ShowerStartPosition, ShowerStartPositionErr, fShowerStartPositionOutputLabel);
90  return 0;
91  }
92 
93  //If we there have none then use the direction to find the neutrino vertex
94  if (ShowerEleHolder.CheckElement(fShowerDirectionInputLabel)) {
95 
96  geo::Vector_t ShowerDirection = {-999, -999, -999};
97  ShowerEleHolder.GetElement(fShowerDirectionInputLabel, ShowerDirection);
98 
100  ShowerEleHolder.GetFindManyP<recob::SpacePoint>(pfpHandle, Event, fPFParticleLabel);
101 
102  //Get the spacepoints handle and the hit assoication
103  auto const spHandle = Event.getValidHandle<std::vector<recob::SpacePoint>>(fPFParticleLabel);
104  const art::FindManyP<recob::Hit>& fmh =
105  ShowerEleHolder.GetFindManyP<recob::Hit>(spHandle, Event, fPFParticleLabel);
106 
107  //Get the spacepoints
108  std::vector<art::Ptr<recob::SpacePoint>> spacePoints_pfp = fmspp.at(pfparticle.key());
109 
110  //Cannot continue if we have no spacepoints
111  if (spacePoints_pfp.empty()) { return 0; }
112 
113  //Get the Shower Center
114  auto const clockData =
116  auto const detProp =
118 
120  clockData, detProp, spacePoints_pfp, fmh);
121 
122  //Order the Hits from the shower centre. The most negative will be the start position.
124  spacePoints_pfp, ShowerCentre, ShowerDirection);
125 
126  //Set the start position.
127  auto ShowerStartPosition = spacePoints_pfp[0]->position();
128 
129  geo::Point_t ShowerStartPositionErr = {-999, -999, -999};
130  ShowerEleHolder.SetElement(
131  ShowerStartPosition, ShowerStartPositionErr, fShowerStartPositionOutputLabel);
132 
133  return 0;
134  }
135 
136  if (fVerbose)
137  mf::LogWarning("ShowerPFPVertexStartPosition")
138  << "Start Position has not been set yet. If you are not calculating the start position "
139  "again then maybe you should stop";
140  return 0;
141  }
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
void SetElement(T &dataproduct, const std::string &Name, bool checktag=false)
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition of vertex object for LArSoft.
Definition: Vertex.h:35
key_type key() const noexcept
Definition: Ptr.h:166
bool CheckElement(const std::string &Name) const
int GetElement(const std::string &Name, T &Element) const
const shower::LArPandoraShowerAlg & GetLArPandoraShowerAlg() const
Definition: IShowerTool.h:82
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
geo::Point_t ShowerCentre(std::vector< art::Ptr< recob::SpacePoint >> const &showersps) const
void OrderShowerSpacePoints(std::vector< art::Ptr< recob::SpacePoint >> &showersps, geo::Point_t const &vertex, geo::Vector_t const &direction) const
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
2D representation of charge deposited in the TDC/wire plane
Definition: Hit.h:46
const art::FindManyP< T1 > & GetFindManyP(const art::ValidHandle< std::vector< T2 >> &handle, const art::Event &evt, const art::InputTag &moduleTag)
Definition: fwd.h:26
const Point_t & position() const
Return vertex 3D position.
Definition: Vertex.h:64
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
const shower::LArPandoraShowerAlg& ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg ( ) const
inlineprotectedinherited
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

art::InputTag ShowerRecoTools::ShowerPFPVertexStartPosition::fPFParticleLabel
private

Definition at line 35 of file ShowerPFPVertexStartPosition_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerPFPVertexStartPosition::fShowerDirectionInputLabel
private

Definition at line 38 of file ShowerPFPVertexStartPosition_tool.cc.

Referenced by CalculateElement().

std::string ShowerRecoTools::ShowerPFPVertexStartPosition::fShowerStartPositionOutputLabel
private

Definition at line 37 of file ShowerPFPVertexStartPosition_tool.cc.

Referenced by CalculateElement().

int ShowerRecoTools::ShowerPFPVertexStartPosition::fVerbose
private

Definition at line 36 of file ShowerPFPVertexStartPosition_tool.cc.

Referenced by CalculateElement().


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