18 #include "NovaDAQConventions/DAQConventions.h" 49 produces< std::vector<rb::Track> >();
53 for(
int i = 0;
i < 2; ++
i)
for(
int j = 0;
j < 2; ++
j)
fTable[
i][
j] = 0;
74 std::unique_ptr<std::vector<rb::Track> > outTracks(
new std::vector<rb::Track>);
81 for(
size_t trIdx = 0; trIdx < tracks->size(); ++trIdx){
83 if(!fmp.at(trIdx).empty())
continue;
89 if(thru) outTracks->push_back(track);
97 const bool ret = !outTracks->empty();
98 evt.
put(std::move(outTracks));
106 std::cout <<
"Throughgoing efficiency table. x-axis selection, y-axis truth:" <<
std::endl;
118 const TVector3 v1 = trk.
Start();
119 const TVector3 v2 = trk.
Stop();
120 const TVector3 stop = (v1.Y() < v2.Y()) ? v1 : v2;
122 if(v1.X() == v2.X() ||
123 v1.Y() == v2.Y())
return false;
125 const double x = stop.X();
126 const double y = stop.Y();
127 const double z = stop.Z();
133 if(z < 20 || z > geom->
DetLength()-20)
return true;
145 if(parts.empty())
return false;
156 if(stop.Z() < 0)
return true;
157 if(stop.Z() > geom->
DetLength())
return true;
back track the reconstruction to the simulation
bool IsTrulyUncontained(const rb::Track &trk) const
fvar< T > fabs(const fvar< T > &x)
const TLorentzVector & EndPosition() const
A rb::Prong with full reconstructed trajectory.
DEFINE_ART_MODULE(TestTMapFile)
virtual TVector3 Start() const
bool IsUncontained(const rb::Track &trk) const
art::PtrVector< rb::CellHit > AllCells() const
Get all cells from both views.
ProductID put(std::unique_ptr< PROD > &&product)
Track finder for cosmic rays.
Selects throughgoing cosmics.
double DetHalfHeight() const
ThroughgoingSelection(const fhicl::ParameterSet &pset)
double DetHalfWidth() const
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
void reconfigure(const fhicl::ParameterSet &pset)
virtual bool filter(art::Event &evt)
TVector3 Stop() const
Position of the final trajectory point.
Encapsulate the geometry of one entire detector (near, far, ndos)
std::vector< const sim::Particle * > HitsToParticle(const std::vector< const rb::CellHit * > &hits) const
Returns vector of sim::Particle objects contributing to the given collection of hits.