38 #include <NovaTimingUtilities/TimingUtilities.h> 109 const HitList& hlist = *hls.at(i);
119 TVector3
end = track.
End();
120 if(start.Y() < end.Y())
std::swap(start, end);
124 (start.X() == 0 && end.X() == 0) ||
125 (start.Y() == 0 && end.Y() == 0) ||
126 (start.Z() == 0 && end.Z() == 0))
133 fOutputEvd <<
"(" << start.X() <<
"," << start.Y() <<
"," << start.Z() <<
") ";
134 fOutputEvd <<
"(" << end.X() <<
"," << end.Y() <<
"," << end.Z() <<
")";
142 for(
uint j = 0;
j < hitlist.size();
j++)
162 double xSpread = end.X() - start.X();
163 double zSpread = end.Z() - start.Z();
165 double zSpreadLoc = z - start.Z();
167 double x = start.X() + xSpread * zSpreadLoc/zSpread;
175 double ySpread = end.Y() - start.Y();
176 double zSpread = end.Z() - start.Z();
178 double zSpreadLoc = z - start.Z();
180 double y = start.Y() + ySpread * zSpreadLoc/zSpread;
DDTEvd(fhicl::ParameterSet const &p)
void WriteTrackInfo(const Track3D &track, long timeStart, const HitList &hitlist)
const uint32_t NOVA_EPOCH
int GetX(DAQHit &h, TVector3 start, TVector3 end)
novaddt::Plane const & Plane() const
novaddt::TDC const & TDC() const
std::vector< DAQHit > HitList
const uint64_t NOVA_TIME_FACTOR
DEFINE_ART_MODULE(TestTMapFile)
fvar< T > round(const fvar< T > &x)
TVector3 const & End() const
void swap(art::HLTGlobalStatus &lhs, art::HLTGlobalStatus &rhs)
void WriteEventInfo(const art::Event &event)
novaddt::View const & View() const
TVector3 const & Start() const
bool const & Is3D() const
EDAnalyzer(Table< Config > const &config)
std::string fInputHitInstance
std::string fInputHitLabel
void analyze(art::Event const &e) override
novaddt::Cell const & Cell() const
int GetY(DAQHit &h, TVector3 start, TVector3 end)
void WriteHitListInfo(const HitList &hitlist, long timeStart, TVector3 start, TVector3 end)
std::string fTrack3DLabel