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

Public Member Functions

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

Private Attributes

art::ServiceHandle< geo::Geometry > fGeom
 
calo::CalorimetryAlg fCalorimetryAlg
 
int fVerbose
 
double fdEdxTrackLength
 
double dEdxTrackLength
 
bool fMaxHitPlane
 
bool fMissFirstPoint
 
std::string fShowerStartPositionInputLabel
 
std::string fInitialTrackHitsInputLabel
 
std::string fShowerDirectionInputLabel
 
std::string fShowerdEdxOutputLabel
 
std::string fShowerBestPlaneOutputLabel
 

Additional Inherited Members

- 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="")
 
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::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)
 

Detailed Description

Definition at line 19 of file ShowerUnidirectiondEdx_tool.cc.

Constructor & Destructor Documentation

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

Definition at line 47 of file ShowerUnidirectiondEdx_tool.cc.

48  : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
49  , fCalorimetryAlg(pset.get<fhicl::ParameterSet>("CalorimetryAlg"))
50  , fVerbose(pset.get<int>("Verbose"))
51  , fdEdxTrackLength(pset.get<float>("dEdxTrackLength"))
52  , fMaxHitPlane(pset.get<bool>("MaxHitPlane"))
53  , fMissFirstPoint(pset.get<bool>("MissFirstPoint"))
54  , fShowerStartPositionInputLabel(pset.get<std::string>("ShowerStartPositionInputLabel"))
55  , fInitialTrackHitsInputLabel(pset.get<std::string>("InitialTrackHitsInputLabel"))
56  , fShowerDirectionInputLabel(pset.get<std::string>("ShowerDirectionInputLabel"))
57  , fShowerdEdxOutputLabel(pset.get<std::string>("ShowerdEdxOutputLabel"))
58  , fShowerBestPlaneOutputLabel(pset.get<std::string>("ShowerBestPlaneOutputLabel"))
59  {}
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::ShowerUnidirectiondEdx::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle,
art::Event &  Event,
reco::shower::ShowerElementHolder &  ShowerEleHolder 
)
overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 62 of file ShowerUnidirectiondEdx_tool.cc.

65  {
66 
67  dEdxTrackLength = fdEdxTrackLength;
68 
69  // Shower dEdx calculation
70  if (!ShowerEleHolder.CheckElement(fShowerStartPositionInputLabel)) {
71  if (fVerbose)
72  mf::LogError("ShowerUnidirectiondEdx") << "Start position not set, returning " << std::endl;
73  return 1;
74  }
75  if (!ShowerEleHolder.CheckElement(fInitialTrackHitsInputLabel)) {
76  if (fVerbose)
77  mf::LogError("ShowerUnidirectiondEdx")
78  << "Initial Track Hits not set returning" << std::endl;
79  return 1;
80  }
81  if (!ShowerEleHolder.CheckElement(fShowerDirectionInputLabel)) {
82  if (fVerbose)
83  mf::LogError("ShowerUnidirectiondEdx") << "Shower Direction not set" << std::endl;
84  return 1;
85  }
86 
87  //Get the initial track hits
88  std::vector<art::Ptr<recob::Hit>> trackhits;
89  ShowerEleHolder.GetElement(fInitialTrackHitsInputLabel, trackhits);
90 
91  if (trackhits.empty()) {
92  if (fVerbose)
93  mf::LogWarning("ShowerUnidirectiondEdx") << "Not Hits in the initial track" << std::endl;
94  return 0;
95  }
96 
97  TVector3 ShowerStartPosition = {-999, -999, -999};
98  ShowerEleHolder.GetElement(fShowerStartPositionInputLabel, ShowerStartPosition);
99 
100  TVector3 showerDir = {-999, -999, -999};
101  ShowerEleHolder.GetElement(fShowerDirectionInputLabel, showerDir);
102 
103  geo::TPCID vtxTPC = fGeom->FindTPCAtPosition(ShowerStartPosition);
104 
105  // Split the track hits per plane
106  std::vector<double> dEdxVec;
107  std::vector<std::vector<art::Ptr<recob::Hit>>> trackHits;
108  unsigned int numPlanes = fGeom->Nplanes();
109  trackHits.resize(numPlanes);
110 
111  // Loop over the track hits and split into planes
112  for (unsigned int hitIt = 0; hitIt < trackhits.size(); ++hitIt) {
113  art::Ptr<recob::Hit> hit = trackhits.at(hitIt);
114  geo::PlaneID hitWire = hit->WireID();
115  geo::TPCID TPC = hitWire.asTPCID();
116 
117  //only get hits from the same TPC as the vertex
118  if (TPC == vtxTPC) { (trackHits.at(hitWire.Plane)).push_back(hit); }
119  }
120 
121  int bestHitsPlane = 0;
122  int bestPlaneHits = 0;
123  int bestPlane = -999;
124  double minPitch = 999;
125 
126  auto const clockData =
127  art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(Event);
128  auto const detProp =
129  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(Event, clockData);
130 
131  for (unsigned int plane = 0; plane < numPlanes; ++plane) {
132  std::vector<art::Ptr<recob::Hit>> trackPlaneHits = trackHits.at(plane);
133 
134  if (trackPlaneHits.size()) {
135 
136  double dEdx = -999;
137  double totQ = 0;
138  double avgT = 0;
139  double pitch = 0;
140 
141  //Calculate the pitch
142  double wirepitch = fGeom->WirePitch(trackPlaneHits.at(0)->WireID().planeID());
143  double angleToVert = fGeom->WireAngleToVertical(fGeom->Plane(plane).View(),
144  trackPlaneHits[0]->WireID().planeID()) -
145  0.5 * TMath::Pi();
146  double cosgamma =
147  std::abs(std::sin(angleToVert) * showerDir.Y() + std::cos(angleToVert) * showerDir.Z());
148 
149  pitch = wirepitch / cosgamma;
150 
151  if (pitch) { // Check the pitch is calculated correctly
152  int nhits = 0;
153  std::vector<float> vQ;
154 
155  //Get the first wire
156  int w0 = trackPlaneHits.at(0)->WireID().Wire;
157 
158  for (auto const& hit : trackPlaneHits) {
159 
160  // Get the wire for each hit
161  int w1 = hit->WireID().Wire;
162  if (fMissFirstPoint && w0 == w1) { continue; }
163 
164  //Ignore hits that are too far away.
165  if (std::abs((w1 - w0) * pitch) < dEdxTrackLength) {
166  vQ.push_back(hit->Integral());
167  totQ += hit->Integral();
168  avgT += hit->PeakTime();
169  ++nhits;
170  }
171  }
172 
173  if (totQ) {
174  // Check if the pitch is the smallest yet (best plane)
175  if (pitch < minPitch) {
176  minPitch = pitch;
177  bestPlane = plane;
178  }
179 
180  //Get the median and calculate the dEdx using the algorithm.
181  double dQdx = TMath::Median(vQ.size(), &vQ[0]) / pitch;
182  dEdx = fCalorimetryAlg.dEdx_AREA(
183  clockData, detProp, dQdx, avgT / nhits, trackPlaneHits.at(0)->WireID().Plane);
184 
185  if (isinf(dEdx)) { dEdx = -999; };
186 
187  if (nhits > bestPlaneHits || ((nhits == bestPlaneHits) && (pitch < minPitch))) {
188  bestHitsPlane = plane;
189  bestPlaneHits = nhits;
190  }
191  }
192  dEdxVec.push_back(dEdx);
193  }
194  else {
195  throw cet::exception("ShowerUnidirectiondEdx")
196  << "pitch is 0. I can't think how it is 0? Stopping so I can tell you" << std::endl;
197  }
198  }
199  else { // if not (trackPlaneHits.size())
200  dEdxVec.push_back(-999);
201  }
202  trackPlaneHits.clear();
203  } //end loop over planes
204 
205  //TODO
206  std::vector<double> dEdxVecErr = {-999, -999, -999};
207 
208  ShowerEleHolder.SetElement(dEdxVec, dEdxVecErr, fShowerdEdxOutputLabel);
209 
210  //Set The best plane
211  if (fMaxHitPlane) { bestPlane = bestHitsPlane; }
212 
213  if (bestPlane == -999) {
214  throw cet::exception("ShowerUnidirectiondEdx") << "No best plane set";
215  }
216  else {
217  ShowerEleHolder.SetElement(bestPlane, fShowerBestPlaneOutputLabel);
218  }
219 
220  if (fVerbose > 1) {
221  std::cout << "Best Plane: " << bestPlane << std::endl;
222  for (unsigned int plane = 0; plane < dEdxVec.size(); plane++) {
223  std::cout << "Plane: " << plane << " with dEdx: " << dEdxVec[plane] << std::endl;
224  }
225  }
226 
227  return 0;
228  }
geo::GeometryCore::Plane
PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const
Returns the specified wire.
Definition: GeometryCore.h:2994
recob::Hit::WireID
geo::WireID WireID() const
Definition: Hit.h:233
geo::PlaneID
The data type to uniquely identify a Plane.
Definition: geo_types.h:472
reco::shower::ShowerElementHolder::SetElement
void SetElement(T &dataproduct, const std::string &Name, bool checktag=false)
Definition: ShowerElementHolder.hh:360
recob::Hit::Integral
float Integral() const
Integral under the calibrated signal waveform of the hit, in tick x ADC units.
Definition: Hit.h:224
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:580
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:211
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:824
geo::GeometryCore::FindTPCAtPosition
geo::TPCID FindTPCAtPosition(double const worldLoc[3]) const
Returns the ID of the TPC at specified location.
Definition: GeometryCore.h:2454
abs
T abs(T value)
Definition: sparse_vector_test.cc:30
geo::PlaneGeo::View
View_t View() const
Which coordinate does this plane measure.
Definition: PlaneGeo.h:184
dEdx
double dEdx(float dqdx, float Efield)
Definition: doAna.cpp:21
geo::GeometryCore::Nplanes
unsigned int Nplanes(unsigned int tpc=0, unsigned int cstat=0) const
Returns the total number of wire planes in the specified TPC.
Definition: GeometryCore.h:2848
push_back
fInnerVessel push_back(Point(-578.400000, 0.000000, 0.000000))
ShowerRecoTools::ShowerUnidirectiondEdx::fGeom
art::ServiceHandle< geo::Geometry > fGeom
Definition: ShowerUnidirectiondEdx_tool.cc:31
reco::shower::ShowerElementHolder::CheckElement
bool CheckElement(const std::string &Name) const
Definition: ShowerElementHolder.hh:415
geo::TPCID
The data type to uniquely identify a TPC.
Definition: geo_types.h:386
geo::PlaneID::Plane
PlaneID_t Plane
Index of the plane within its TPC.
Definition: geo_types.h:493
reco::shower::ShowerElementHolder::GetElement
int GetElement(const std::string &Name, T &Element) const
Definition: ShowerElementHolder.hh:219
hit
Detector simulation of raw signals on wires.
Definition: NumberOfHitsFilter_module.cc:25
geo::TPCID::asTPCID
constexpr TPCID const & asTPCID() const
Conversion to TPCID (for convenience of notation).
Definition: geo_types.h:446
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:218
mf::LogWarning
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
Definition: MessageLogger.h:210
calo::CalorimetryAlg::dEdx_AREA
double dEdx_AREA(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, recob::Hit const &hit, double pitch, double T0=0) const
Definition: CalorimetryAlg.cxx:131
geo::GeometryCore::WireAngleToVertical
double WireAngleToVertical(geo::View_t view, geo::TPCID const &tpcid) const
Returns the angle of the wires in the specified view from vertical.
Definition: GeometryCore.cxx:843
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030

Member Data Documentation

double ShowerRecoTools::ShowerUnidirectiondEdx::dEdxTrackLength
private

Definition at line 36 of file ShowerUnidirectiondEdx_tool.cc.

calo::CalorimetryAlg ShowerRecoTools::ShowerUnidirectiondEdx::fCalorimetryAlg
private

Definition at line 32 of file ShowerUnidirectiondEdx_tool.cc.

double ShowerRecoTools::ShowerUnidirectiondEdx::fdEdxTrackLength
private

Definition at line 36 of file ShowerUnidirectiondEdx_tool.cc.

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

Definition at line 31 of file ShowerUnidirectiondEdx_tool.cc.

std::string ShowerRecoTools::ShowerUnidirectiondEdx::fInitialTrackHitsInputLabel
private

Definition at line 41 of file ShowerUnidirectiondEdx_tool.cc.

bool ShowerRecoTools::ShowerUnidirectiondEdx::fMaxHitPlane
private

Definition at line 38 of file ShowerUnidirectiondEdx_tool.cc.

bool ShowerRecoTools::ShowerUnidirectiondEdx::fMissFirstPoint
private

Definition at line 39 of file ShowerUnidirectiondEdx_tool.cc.

std::string ShowerRecoTools::ShowerUnidirectiondEdx::fShowerBestPlaneOutputLabel
private

Definition at line 44 of file ShowerUnidirectiondEdx_tool.cc.

std::string ShowerRecoTools::ShowerUnidirectiondEdx::fShowerdEdxOutputLabel
private

Definition at line 43 of file ShowerUnidirectiondEdx_tool.cc.

std::string ShowerRecoTools::ShowerUnidirectiondEdx::fShowerDirectionInputLabel
private

Definition at line 42 of file ShowerUnidirectiondEdx_tool.cc.

std::string ShowerRecoTools::ShowerUnidirectiondEdx::fShowerStartPositionInputLabel
private

Definition at line 40 of file ShowerUnidirectiondEdx_tool.cc.

int ShowerRecoTools::ShowerUnidirectiondEdx::fVerbose
private

Definition at line 35 of file ShowerUnidirectiondEdx_tool.cc.

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