Public Member Functions | Private Member Functions | Private Attributes | List of all members
ems::EMShower3D Class Reference
Inheritance diagram for ems::EMShower3D:
art::EDProducer art::detail::Producer art::detail::LegacyModule art::Modifier art::ModuleBase art::ProductRegistryHelper

Public Member Functions

 EMShower3D (fhicl::ParameterSet const &p)
 
 EMShower3D (EMShower3D const &)=delete
 
 EMShower3D (EMShower3D &&)=delete
 
EMShower3D & operator= (EMShower3D const &)=delete
 
EMShower3D & operator= (EMShower3D &&)=delete
 
- Public Member Functions inherited from art::EDProducer
 EDProducer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 EDProducer (Table< Config > const &config)
 
std::string workerType () const
 
- Public Member Functions inherited from art::detail::Producer
virtual ~Producer () noexcept
 
 Producer (fhicl::ParameterSet const &)
 
 Producer (Producer const &)=delete
 
 Producer (Producer &&)=delete
 
Producer & operator= (Producer const &)=delete
 
Producer & operator= (Producer &&)=delete
 
void doBeginJob (SharedResources const &resources)
 
void doEndJob ()
 
void doRespondToOpenInputFile (FileBlock const &fb)
 
void doRespondToCloseInputFile (FileBlock const &fb)
 
void doRespondToOpenOutputFiles (FileBlock const &fb)
 
void doRespondToCloseOutputFiles (FileBlock const &fb)
 
bool doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doEndRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
 
- Public Member Functions inherited from art::Modifier
 ~Modifier () noexcept
 
 Modifier ()
 
 Modifier (Modifier const &)=delete
 
 Modifier (Modifier &&)=delete
 
Modifier & operator= (Modifier const &)=delete
 
Modifier & operator= (Modifier &&)=delete
 
- Public Member Functions inherited from art::ModuleBase
virtual ~ModuleBase () noexcept
 
 ModuleBase ()
 
ModuleDescription const & moduleDescription () const
 
void setModuleDescription (ModuleDescription const &)
 
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables () const
 
void sortConsumables (std::string const &current_process_name)
 
template<typename T , BranchType BT>
ViewToken< T > consumesView (InputTag const &tag)
 
template<typename T , BranchType BT>
ViewToken< T > mayConsumeView (InputTag const &tag)
 

Private Member Functions

void produce (art::Event &e) override
 
recob::Track ConvertFrom (detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, pma::Track3D &src)
 
recob::Track ConvertFrom2 (detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, pma::Track3D &src)
 
recob::Cluster ConvertFrom (const std::vector< art::Ptr< recob::Hit >> &src)
 
std::vector< ems::DirOfGamma * > CollectShower2D (detinfo::DetectorPropertiesData const &detProp, art::Event const &e)
 
void Link (art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > input)
 
void Reoptimize (detinfo::DetectorPropertiesData const &detProp)
 
void Make3DSeg (art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > pair)
 
bool Validate (detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > input, size_t id1, size_t id2, size_t c1, size_t c2, size_t plane3)
 
void FilterOutSmallParts (detinfo::DetectorPropertiesData const &detProp, double r2d, const std::vector< art::Ptr< recob::Hit >> &hits_in, std::vector< art::Ptr< recob::Hit >> &hits_out)
 
bool GetCloseHits (detinfo::DetectorPropertiesData const &detProp, double r2d, const std::vector< art::Ptr< recob::Hit >> &hits_in, std::vector< size_t > &used, std::vector< art::Ptr< recob::Hit >> &hits_out)
 
bool Has (const std::vector< size_t > &v, size_t idx)
 
size_t LinkCandidates (art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > input, size_t id)
 

Private Attributes

std::vector< IniSeg > fInisegs
 
std::vector< IniSeg > fSeltracks
 
std::vector< IniSeg > fPMA3D
 
std::vector< std::vector< art::Ptr< recob::Hit > > > fClusters
 
std::vector< size_t > fClustersNotUsed
 
std::vector< size_t > fTracksNotUsed
 
unsigned int fTrkIndex
 
unsigned int fClIndex
 
unsigned int fIniIndex
 
std::string fCluModuleLabel
 
std::string fTrk3DModuleLabel
 
pma::ProjectionMatchingAlg fProjectionMatchingAlg
 
calo::CalorimetryAlg fCalorimetryAlg
 
art::Handle< std::vector< recob::Cluster > > fCluListHandle
 

Additional Inherited Members

- Public Types inherited from art::EDProducer
using ModuleType = EDProducer
 
using WorkerType = WorkerT< EDProducer >
 
- Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using Table = Modifier::Table< UserConfig, KeysToIgnore >
 
- Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >
 
- Static Public Member Functions inherited from art::EDProducer
static void commitEvent (EventPrincipal &ep, Event &e)
 
- Protected Member Functions inherited from art::ModuleBase
ConsumesCollector & consumesCollector ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 

Detailed Description

Definition at line 52 of file EMShower3D_module.cc.

Constructor & Destructor Documentation

ems::EMShower3D::EMShower3D ( fhicl::ParameterSet const &  p)
explicit

Definition at line 134 of file EMShower3D_module.cc.

135  : EDProducer{p}
136  , fProjectionMatchingAlg(p.get<fhicl::ParameterSet>("ProjectionMatchingAlg"))
137  , fCalorimetryAlg(p.get<fhicl::ParameterSet>("CalorimetryAlg"))
138 {
139  fCluModuleLabel = p.get<std::string>("ClustersModuleLabel");
140  fTrk3DModuleLabel = p.get<std::string>("Trk3DModuleLabel");
141 
142  produces<std::vector<recob::Track>>();
143  produces<std::vector<recob::Vertex>>();
144  produces<std::vector<recob::Cluster>>();
145  produces<std::vector<recob::SpacePoint>>();
146  produces<art::Assns<recob::Track, recob::Hit>>();
147  produces<art::Assns<recob::Track, recob::Vertex>>();
148  produces<art::Assns<recob::Cluster, recob::Hit>>();
149  produces<art::Assns<recob::Track, recob::SpacePoint>>();
150  produces<art::Assns<recob::SpacePoint, recob::Hit>>();
151  produces<art::Assns<recob::Track, recob::Cluster>>();
152 }
string
std::string string
Definition: nybbler.cc:12
art::EDProducer::EDProducer
EDProducer(fhicl::ParameterSet const &pset)
Definition: EDProducer.h:20
ems::EMShower3D::fTrk3DModuleLabel
std::string fTrk3DModuleLabel
Definition: EMShower3D_module.cc:126
test.p
p
Definition: test.py:223
ems::EMShower3D::fCalorimetryAlg
calo::CalorimetryAlg fCalorimetryAlg
Definition: EMShower3D_module.cc:129
ems::EMShower3D::fProjectionMatchingAlg
pma::ProjectionMatchingAlg fProjectionMatchingAlg
Definition: EMShower3D_module.cc:128
ems::EMShower3D::fCluModuleLabel
std::string fCluModuleLabel
Definition: EMShower3D_module.cc:125
ems::EMShower3D::EMShower3D ( EMShower3D const &  )
delete
ems::EMShower3D::EMShower3D ( EMShower3D &&  )
delete

Member Function Documentation

std::vector< ems::DirOfGamma * > ems::EMShower3D::CollectShower2D ( detinfo::DetectorPropertiesData const &  detProp,
art::Event const &  e 
)
private

Definition at line 590 of file EMShower3D_module.cc.

592 {
593  std::vector<ems::DirOfGamma*> input;
594 
595  if (e.getByLabel(fCluModuleLabel, fCluListHandle)) {
596  art::FindManyP<recob::Hit> fb(fCluListHandle, e, fCluModuleLabel);
597  for (unsigned int c = 0; c < fCluListHandle->size(); c++) {
598  std::vector<art::Ptr<recob::Hit>> hitlist;
599  hitlist = fb.at(c);
600 
601  if (hitlist.size() > 5) {
602  std::vector<art::Ptr<recob::Hit>> hits_out;
603  FilterOutSmallParts(detProp, 2.0, hitlist, hits_out);
604 
605  if (hits_out.size() > 5) {
606  fClusters.push_back(hits_out);
607 
608  ems::DirOfGamma* sh = new ems::DirOfGamma(detProp, hits_out, 14, c);
609 
610  if (sh->GetHits2D().size()) input.push_back(sh);
611  }
612  }
613  }
614  }
615 
616  return input;
617 }
ems::EMShower3D::fClusters
std::vector< std::vector< art::Ptr< recob::Hit > > > fClusters
Definition: EMShower3D_module.cc:116
ems::EMShower3D::fCluListHandle
art::Handle< std::vector< recob::Cluster > > fCluListHandle
Definition: EMShower3D_module.cc:131
ems::EMShower3D::FilterOutSmallParts
void FilterOutSmallParts(detinfo::DetectorPropertiesData const &detProp, double r2d, const std::vector< art::Ptr< recob::Hit >> &hits_in, std::vector< art::Ptr< recob::Hit >> &hits_out)
Definition: EMShower3D_module.cc:957
ems::DirOfGamma::GetHits2D
std::vector< Hit2D * > const & GetHits2D() const
Definition: DirOfGamma.h:200
e
const double e
Definition: gUpMuFluxGen.cxx:165
input
static int input(void)
Definition: code.cpp:15695
ems::EMShower3D::fCluModuleLabel
std::string fCluModuleLabel
Definition: EMShower3D_module.cc:125
recob::Track ems::EMShower3D::ConvertFrom ( detinfo::DetectorClocksData const &  clock_data,
detinfo::DetectorPropertiesData const &  det_prop,
pma::Track3D &  src 
)
private

Definition at line 184 of file EMShower3D_module.cc.

187 {
188  auto const* geom = lar::providerFrom<geo::Geometry>();
189  double avdrift = (src.front()->Point3D().X() + src.back()->Point3D().X()) * 0.5;
190  unsigned int nplanes = geom->Nplanes(src.front()->TPC(), src.front()->Cryo());
191  size_t nusedhitsmax = 0;
192  int bestplane = -1;
193  for (unsigned int p = 0; p < nplanes; ++p) {
194  unsigned int nusedP = 0;
195  fProjectionMatchingAlg.selectInitialHits(src, p, &nusedP);
196 
197  if (nusedP > nusedhitsmax) {
198  nusedhitsmax = nusedP;
199  bestplane = int(p);
200  }
201  }
202 
203  std::vector<std::vector<double>> vdedx;
204  std::vector<double> dedx;
205 
206  for (unsigned int p = 0; p < nplanes; ++p) {
207  unsigned int nusedP = 0;
208  double dqdxplane = fProjectionMatchingAlg.selectInitialHits(src, p, &nusedP);
209  double timeP = detProp.ConvertXToTicks(avdrift, p, src.front()->TPC(), src.front()->Cryo());
210  double dEdxplane = fCalorimetryAlg.dEdx_AREA(clock_data, detProp, dqdxplane, timeP, p);
211  dedx.push_back(dEdxplane);
212  if (int(p) == bestplane)
213  dedx.push_back(1);
214  else
215  dedx.push_back(0);
216  vdedx.push_back(dedx);
217  }
218 
219  std::vector<TVector3> xyz, dircos;
220 
221  for (size_t i = 0; i < src.size(); i++) {
222  xyz.push_back(src[i]->Point3D());
223 
224  if (i < src.size() - 1) {
225  TVector3 dc(src[i + 1]->Point3D());
226  dc -= src[i]->Point3D();
227  dc *= 1.0 / dc.Mag();
228  dircos.push_back(dc);
229  }
230  else
231  dircos.push_back(dircos.back());
232  }
233 
234  return recob::Track(recob::TrackTrajectory(recob::tracking::convertCollToPoint(xyz),
235  recob::tracking::convertCollToVector(dircos),
236  recob::Track::Flags_t(xyz.size()),
237  false),
238  0,
239  -1.,
240  0,
241  recob::tracking::SMatrixSym55(),
242  recob::tracking::SMatrixSym55(),
243  fIniIndex);
244 }
recob::tracking::SMatrixSym55
ROOT::Math::SMatrix< Double32_t, 5, 5, ROOT::Math::MatRepSym< Double32_t, 5 > > SMatrixSym55
Definition: TrackingTypes.h:85
recob::tracking::convertCollToVector
std::vector< Vector_t > convertCollToVector(std::vector< Vector > const &coll)
Definition: TrackingTypes.h:58
pma::Track3D::front
pma::Hit3D const * front() const
Definition: PmaTrack3D.h:98
lar::fhicl::geo::Point3D
::fhicl::TupleAs< Point(::geo::Length_t,::geo::Length_t,::geo::Length_t)> Point3D
Atom object for reading a 3D point or vector (centimeters).
Definition: geo_vectors_fhicl.h:38
pma::Hit3D::Cryo
unsigned int Cryo() const noexcept
Definition: PmaHit3D.h:83
recob::Track::Flags_t
TrackTrajectory::Flags_t Flags_t
Definition: Track.h:68
pma::Hit3D::Point3D
TVector3 const & Point3D() const
Definition: PmaHit3D.h:55
ems::EMShower3D::fIniIndex
unsigned int fIniIndex
Definition: EMShower3D_module.cc:123
pma::ProjectionMatchingAlg::selectInitialHits
double selectInitialHits(pma::Track3D &trk, unsigned int view=geo::kZ, unsigned int *nused=0) const
Definition: ProjectionMatchingAlg.cxx:1376
recob::TrackTrajectory
A trajectory in space reconstructed from hits.
Definition: TrackTrajectory.h:62
test.p
p
Definition: test.py:223
recob::tracking::convertCollToPoint
std::vector< Point_t > convertCollToPoint(std::vector< Point > const &coll)
Definition: TrackingTypes.h:55
ems::EMShower3D::fCalorimetryAlg
calo::CalorimetryAlg fCalorimetryAlg
Definition: EMShower3D_module.cc:129
keras_to_tensorflow.int
int
Definition: keras_to_tensorflow.py:69
ems::EMShower3D::fProjectionMatchingAlg
pma::ProjectionMatchingAlg fProjectionMatchingAlg
Definition: EMShower3D_module.cc:128
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
pma::Track3D::back
pma::Hit3D const * back() const
Definition: PmaTrack3D.h:103
proxy::Track
TrackCollectionProxyElement< TrackCollProxy > Track
Proxy to an element of a proxy collection of recob::Track objects.
Definition: Track.h:1036
pma::Hit3D::TPC
unsigned int TPC() const noexcept
Definition: PmaHit3D.h:88
pma::Track3D::size
size_t size() const
Definition: PmaTrack3D.h:108
recob::Cluster ems::EMShower3D::ConvertFrom ( const std::vector< art::Ptr< recob::Hit >> &  src)
private

Definition at line 155 of file EMShower3D_module.cc.

156 {
157  return recob::Cluster(0.0F,
158  0.0F,
159  0.0F,
160  0.0F,
161  0.0F,
162  0.0F,
163  0.0F,
164  0.0F,
165  0.0F,
166  0.0F,
167  0.0F,
168  0.0F,
169  0.0F,
170  0.0F,
171  0.0F,
172  0.0F,
173  0.0F,
174  0.0F,
175  src.size(),
176  0.0F,
177  0.0F,
178  fClIndex,
179  src[0]->View(),
180  src[0]->WireID().planeID());
181 }
recob::Cluster
Set of hits with a 2D structure.
Definition: Cluster.h:71
ems::EMShower3D::fClIndex
unsigned int fClIndex
Definition: EMShower3D_module.cc:122
recob::Track ems::EMShower3D::ConvertFrom2 ( detinfo::DetectorClocksData const &  clock_data,
detinfo::DetectorPropertiesData const &  det_prop,
pma::Track3D &  src 
)
private

Definition at line 247 of file EMShower3D_module.cc.

250 {
251  auto const* geom = lar::providerFrom<geo::Geometry>();
252 
253  double avdrift = (src.front()->Point3D().X() + src.back()->Point3D().X()) * 0.5;
254  unsigned int nplanes = geom->Nplanes(src.front()->TPC(), src.front()->Cryo());
255  size_t nusedhitsmax = 0;
256  int bestplane = -1;
257  for (unsigned int p = 0; p < nplanes; ++p) {
258  unsigned int nusedP = 0;
259  fProjectionMatchingAlg.selectInitialHits(src, p, &nusedP);
260 
261  if (nusedP > nusedhitsmax) {
262  nusedhitsmax = nusedP;
263  bestplane = int(p);
264  }
265  }
266 
267  std::vector<std::vector<double>> vdedx;
268  std::vector<double> dedx;
269 
270  for (unsigned int p = 0; p < nplanes; ++p) {
271  unsigned int nusedP = 0;
272  double dqdxplane = fProjectionMatchingAlg.selectInitialHits(src, p, &nusedP);
273  double timeP = detProp.ConvertXToTicks(avdrift, p, src.front()->TPC(), src.front()->Cryo());
274  double dEdxplane = fCalorimetryAlg.dEdx_AREA(clockData, detProp, dqdxplane, timeP, p);
275  dedx.push_back(dEdxplane);
276  if (int(p) == bestplane)
277  dedx.push_back(1);
278  else
279  dedx.push_back(0);
280  vdedx.push_back(dedx);
281  }
282 
283  std::vector<TVector3> xyz, dircos;
284 
285  for (size_t i = 0; i < src.size(); i++) {
286  xyz.push_back(src[i]->Point3D());
287 
288  if (i < src.size() - 1) {
289  TVector3 dc(src[i + 1]->Point3D());
290  dc -= src[i]->Point3D();
291  dc *= 1.0 / dc.Mag();
292  dircos.push_back(dc);
293  }
294  else
295  dircos.push_back(dircos.back());
296  }
297 
298  return recob::Track(recob::TrackTrajectory(recob::tracking::convertCollToPoint(xyz),
299  recob::tracking::convertCollToVector(dircos),
300  recob::Track::Flags_t(xyz.size()),
301  false),
302  0,
303  -1.,
304  0,
305  recob::tracking::SMatrixSym55(),
306  recob::tracking::SMatrixSym55(),
307  fIniIndex);
308 }
recob::tracking::SMatrixSym55
ROOT::Math::SMatrix< Double32_t, 5, 5, ROOT::Math::MatRepSym< Double32_t, 5 > > SMatrixSym55
Definition: TrackingTypes.h:85
recob::tracking::convertCollToVector
std::vector< Vector_t > convertCollToVector(std::vector< Vector > const &coll)
Definition: TrackingTypes.h:58
pma::Track3D::front
pma::Hit3D const * front() const
Definition: PmaTrack3D.h:98
lar::fhicl::geo::Point3D
::fhicl::TupleAs< Point(::geo::Length_t,::geo::Length_t,::geo::Length_t)> Point3D
Atom object for reading a 3D point or vector (centimeters).
Definition: geo_vectors_fhicl.h:38
pma::Hit3D::Cryo
unsigned int Cryo() const noexcept
Definition: PmaHit3D.h:83
recob::Track::Flags_t
TrackTrajectory::Flags_t Flags_t
Definition: Track.h:68
pma::Hit3D::Point3D
TVector3 const & Point3D() const
Definition: PmaHit3D.h:55
ems::EMShower3D::fIniIndex
unsigned int fIniIndex
Definition: EMShower3D_module.cc:123
pma::ProjectionMatchingAlg::selectInitialHits
double selectInitialHits(pma::Track3D &trk, unsigned int view=geo::kZ, unsigned int *nused=0) const
Definition: ProjectionMatchingAlg.cxx:1376
recob::TrackTrajectory
A trajectory in space reconstructed from hits.
Definition: TrackTrajectory.h:62
test.p
p
Definition: test.py:223
recob::tracking::convertCollToPoint
std::vector< Point_t > convertCollToPoint(std::vector< Point > const &coll)
Definition: TrackingTypes.h:55
ems::EMShower3D::fCalorimetryAlg
calo::CalorimetryAlg fCalorimetryAlg
Definition: EMShower3D_module.cc:129
keras_to_tensorflow.int
int
Definition: keras_to_tensorflow.py:69
ems::EMShower3D::fProjectionMatchingAlg
pma::ProjectionMatchingAlg fProjectionMatchingAlg
Definition: EMShower3D_module.cc:128
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
pma::Track3D::back
pma::Hit3D const * back() const
Definition: PmaTrack3D.h:103
proxy::Track
TrackCollectionProxyElement< TrackCollProxy > Track
Proxy to an element of a proxy collection of recob::Track objects.
Definition: Track.h:1036
pma::Hit3D::TPC
unsigned int TPC() const noexcept
Definition: PmaHit3D.h:88
pma::Track3D::size
size_t size() const
Definition: PmaTrack3D.h:108
void ems::EMShower3D::FilterOutSmallParts ( detinfo::DetectorPropertiesData const &  detProp,
double  r2d,
const std::vector< art::Ptr< recob::Hit >> &  hits_in,
std::vector< art::Ptr< recob::Hit >> &  hits_out 
)
private

Definition at line 957 of file EMShower3D_module.cc.

961 {
962  size_t min_size = hits_in.size() / 5;
963  if (min_size < 3) min_size = 3;
964 
965  std::vector<size_t> used;
966  std::vector<art::Ptr<recob::Hit>> close_hits;
967 
968  while (GetCloseHits(detProp, r2d, hits_in, used, close_hits)) {
969  if (close_hits.size() > min_size)
970  for (auto h : close_hits)
971  hits_out.push_back(h);
972  }
973 }
reco_momentum_tuples.h
h
Definition: reco_momentum_tuples.py:66
ems::EMShower3D::GetCloseHits
bool GetCloseHits(detinfo::DetectorPropertiesData const &detProp, double r2d, const std::vector< art::Ptr< recob::Hit >> &hits_in, std::vector< size_t > &used, std::vector< art::Ptr< recob::Hit >> &hits_out)
Definition: EMShower3D_module.cc:881
bool ems::EMShower3D::GetCloseHits ( detinfo::DetectorPropertiesData const &  detProp,
double  r2d,
const std::vector< art::Ptr< recob::Hit >> &  hits_in,
std::vector< size_t > &  used,
std::vector< art::Ptr< recob::Hit >> &  hits_out 
)
private

Definition at line 881 of file EMShower3D_module.cc.

886 {
887 
888  hits_out.clear();
889 
890  const double gapMargin = 5.0; // can be changed to f(id_tpc1, id_tpc2)
891  size_t idx = 0;
892 
893  while ((idx < hits_in.size()) && Has(used, idx))
894  idx++;
895 
896  if (idx < hits_in.size()) {
897  hits_out.push_back(hits_in[idx]);
898  used.push_back(idx);
899 
900  double r2d2 = r2d * r2d;
901  double gapMargin2 = sqrt(2 * gapMargin * gapMargin);
902  gapMargin2 = (gapMargin2 + r2d) * (gapMargin2 + r2d);
903 
904  bool collect = true;
905  while (collect) {
906  collect = false;
907  for (size_t i = 0; i < hits_in.size(); i++)
908  if (!Has(used, i)) {
909  art::Ptr<recob::Hit> hi = hits_in[i];
910  TVector2 hi_cm = pma::WireDriftToCm(detProp,
911  hi->WireID().Wire,
912  hi->PeakTime(),
913  hi->WireID().Plane,
914  hi->WireID().TPC,
915  hi->WireID().Cryostat);
916 
917  bool accept = false;
918  // for (auto const& ho : hits_out)
919  for (size_t idx_o = 0; idx_o < hits_out.size(); idx_o++) {
920  art::Ptr<recob::Hit> ho = hits_out[idx_o];
921 
922  double d2 = pma::Dist2(hi_cm,
923  pma::WireDriftToCm(detProp,
924  ho->WireID().Wire,
925  ho->PeakTime(),
926  ho->WireID().Plane,
927  ho->WireID().TPC,
928  ho->WireID().Cryostat));
929 
930  if (hi->WireID().TPC == ho->WireID().TPC) {
931  if (d2 < r2d2) {
932  accept = true;
933  break;
934  }
935  }
936  else {
937  if (d2 < gapMargin2) {
938  accept = true;
939  break;
940  }
941  }
942  }
943  if (accept) {
944  collect = true;
945  hits_out.push_back(hi);
946  used.push_back(i);
947  }
948  }
949  }
950  return true;
951  }
952  else
953  return false;
954 }
pma::Dist2
double Dist2(const TVector2 &v1, const TVector2 &v2)
Definition: Utilities.cxx:37
recob::Hit::WireID
geo::WireID WireID() const
Definition: Hit.h:233
geo::CryostatID::Cryostat
CryostatID_t Cryostat
Index of cryostat.
Definition: geo_types.h:212
pma::WireDriftToCm
TVector2 WireDriftToCm(detinfo::DetectorPropertiesData const &detProp, unsigned int wire, float drift, unsigned int plane, unsigned int tpc, unsigned int cryo)
Definition: Utilities.cxx:294
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:580
geo::PlaneID::Plane
PlaneID_t Plane
Index of the plane within its TPC.
Definition: geo_types.h:493
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:218
ems::EMShower3D::Has
bool Has(const std::vector< size_t > &v, size_t idx)
Definition: EMShower3D_module.cc:873
geo::TPCID::TPC
TPCID_t TPC
Index of the TPC within its cryostat.
Definition: geo_types.h:406
bool ems::EMShower3D::Has ( const std::vector< size_t > &  v,
size_t  idx 
)
private

Definition at line 873 of file EMShower3D_module.cc.

874 {
875  for (auto c : v)
876  if (c == idx) return true;
877  return false;
878 }
void ems::EMShower3D::Link ( art::Event const &  e,
detinfo::DetectorPropertiesData const &  detProp,
std::vector< ems::DirOfGamma * >  input 
)
private

Definition at line 620 of file EMShower3D_module.cc.

623 {
624  std::vector<std::vector<size_t>> saveids;
625  std::vector<size_t> saveidsnotusedcls;
626  size_t i = 0;
627 
628  while (i < input.size()) {
629  if (!input[i]->GetCandidates().size()) {
630  i++;
631  continue;
632  }
633 
634  double mindist = 1.0; // cm
635  std::vector<ems::DirOfGamma*> pairs;
636 
637  size_t startview = input[i]->GetFirstHit()->WireID().Plane;
638  size_t tpc = input[i]->GetFirstHit()->WireID().TPC;
639  size_t cryo = input[i]->GetFirstHit()->WireID().Cryostat;
640 
641  float t1 = detProp.ConvertTicksToX(input[i]->GetFirstHit()->PeakTime(), startview, tpc, cryo);
642 
643  unsigned int idsave = 0;
644  for (unsigned int j = 0; j < input.size(); j++) {
645  if (!input[j]->GetCandidates().size()) continue;
646 
647  size_t secondview = input[j]->GetFirstHit()->WireID().Plane;
648  size_t tpc_j = input[j]->GetFirstHit()->WireID().TPC;
649  size_t cryo_j = input[j]->GetFirstHit()->WireID().Cryostat;
650 
651  if ((i != j) && (secondview != startview) && (tpc == tpc_j) && (cryo == cryo_j)) {
652  float t2 =
653  detProp.ConvertTicksToX(input[j]->GetFirstHit()->PeakTime(), secondview, tpc_j, cryo_j);
654  float dist = fabs(t2 - t1);
655 
656  if (dist < mindist) {
657  mindist = dist;
658  pairs.clear();
659  pairs.push_back(input[i]);
660  pairs.push_back(input[j]);
661  idsave = j;
662  }
663  }
664  }
665 
666  bool exist = false;
667  for (unsigned int v = 0; v < saveids.size(); v++)
668  if ((saveids[v][0] == i) || (saveids[v][0] == idsave))
669  if ((saveids[v][1] == i) || (saveids[v][1] == idsave)) exist = true;
670 
671  if (pairs.size()) {
672  if (!exist) Make3DSeg(e, detProp, pairs);
673 
674  std::vector<size_t> ids;
675  ids.push_back(i);
676  ids.push_back(idsave);
677  saveids.push_back(ids);
678  }
679  else {
680  saveidsnotusedcls.push_back(i);
681  }
682 
683  i++;
684  }
685 
686  i = 0;
687  while (i < saveidsnotusedcls.size()) {
688  LinkCandidates(e, detProp, input, i);
689  i++;
690  }
691 }
ems::EMShower3D::LinkCandidates
size_t LinkCandidates(art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > input, size_t id)
Definition: EMShower3D_module.cc:694
e
const double e
Definition: gUpMuFluxGen.cxx:165
ems::EMShower3D::Make3DSeg
void Make3DSeg(art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > pair)
Definition: EMShower3D_module.cc:774
test_gpu_visible.t1
t1
Definition: test_gpu_visible.py:16
larg4::dist
constexpr double dist(const TReal *x, const TReal *y, const unsigned int dimension)
Definition: OpFastScintillation.hh:563
size_t ems::EMShower3D::LinkCandidates ( art::Event const &  e,
detinfo::DetectorPropertiesData const &  detProp,
std::vector< ems::DirOfGamma * >  input,
size_t  id 
)
private

Definition at line 694 of file EMShower3D_module.cc.

698 {
699  art::ServiceHandle<geo::Geometry const> geom;
700 
701  size_t index = id;
702  bool found = false;
703 
704  if (input[id]->GetCandidates().size() < 2) { return index; }
705 
706  double mindist = 3.0; // cm
707  std::vector<ems::DirOfGamma*> pairs;
708 
709  size_t idcsave = 0;
710  size_t idcjsave = 0;
711  size_t c = 0;
712  size_t idsave = 0;
713  while (c < input[id]->GetCandidates().size()) {
714 
715  size_t startview = input[id]->GetCandidates()[c].GetPlane();
716  size_t tpc = input[id]->GetCandidates()[c].GetTPC();
717  size_t cryo = input[id]->GetCandidates()[c].GetCryo();
718 
719  float t1 = input[id]->GetCandidates()[c].GetPosition().Y(); // y --> drift in 2D space.
720 
721  // loop over 2D showers
722  for (size_t j = 0; j < input.size(); ++j) {
723  if (!input[j]->GetCandidates().size()) continue;
724  if (j == id) continue;
725 
726  // loop over candidates
727  for (size_t cj = 0; cj < input[j]->GetCandidates().size(); ++cj) {
728  size_t secondview = input[j]->GetCandidates()[cj].GetPlane();
729  size_t tpc_j = input[j]->GetCandidates()[cj].GetTPC();
730  size_t cryo_j = input[j]->GetCandidates()[cj].GetCryo();
731 
732  size_t thirdview = startview;
733 
734  const geo::CryostatGeo& cryostat = geom->Cryostat(cryo);
735  for (size_t p = 0; p < cryostat.MaxPlanes(); p++)
736  if ((p == startview) || (p == secondview)) { continue; }
737  else {
738  thirdview = p;
739  break;
740  }
741 
742  if ((startview != secondview) && (tpc == tpc_j) && (cryo == cryo_j)) {
743  float t2 = input[j]->GetCandidates()[cj].GetPosition().Y();
744  float dist = fabs(t2 - t1);
745 
746  if ((dist < mindist) && Validate(detProp, input, id, j, c, cj, thirdview)) {
747  mindist = dist;
748  pairs.clear();
749  pairs.push_back(input[id]);
750  pairs.push_back(input[j]);
751  idsave = j;
752  index = j;
753  idcsave = c;
754  idcjsave = cj;
755  found = true;
756  }
757  }
758  }
759  }
760 
761  c++;
762  }
763 
764  if (found && pairs.size()) {
765  input[id]->SetIdCandidate(idcsave);
766  input[idsave]->SetIdCandidate(idcjsave);
767  Make3DSeg(e, detProp, pairs);
768  }
769 
770  return index;
771 }
geo::CryostatGeo::MaxPlanes
unsigned int MaxPlanes() const
Returns the largest number of planes among the TPCs in this cryostat.
Definition: CryostatGeo.cxx:170
geo::CryostatGeo
Geometry information for a single cryostat.
Definition: CryostatGeo.h:43
ems::EMShower3D::Validate
bool Validate(detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > input, size_t id1, size_t id2, size_t c1, size_t c2, size_t plane3)
Definition: EMShower3D_module.cc:823
util::size
decltype(auto) constexpr size(T &&obj)
ADL-aware version of std::size.
Definition: StdUtils.h:92
e
const double e
Definition: gUpMuFluxGen.cxx:165
ems::EMShower3D::Make3DSeg
void Make3DSeg(art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > pair)
Definition: EMShower3D_module.cc:774
test_gpu_visible.t1
t1
Definition: test_gpu_visible.py:16
test.p
p
Definition: test.py:223
geo::GeometryCore::Cryostat
CryostatGeo const & Cryostat(geo::CryostatID const &cryoid) const
Returns the specified cryostat.
Definition: GeometryCore.cxx:321
larg4::dist
constexpr double dist(const TReal *x, const TReal *y, const unsigned int dimension)
Definition: OpFastScintillation.hh:563
void ems::EMShower3D::Make3DSeg ( art::Event const &  e,
detinfo::DetectorPropertiesData const &  detProp,
std::vector< ems::DirOfGamma * >  pair 
)
private

Definition at line 774 of file EMShower3D_module.cc.

777 {
778  if (pair.size() < 2) return;
779 
780  // to build a track correctly 2d hits must belong to the same tpc
781  size_t tpc1 = pair[0]->GetFirstHit()->WireID().TPC;
782  size_t tpc2 = pair[1]->GetFirstHit()->WireID().TPC;
783 
784  std::vector<art::Ptr<recob::Hit>> vec1 = pair[0]->GetIniHits();
785  std::vector<art::Ptr<recob::Hit>> vec2 = pair[1]->GetIniHits();
786 
787  if ((vec1.size() < 3) && (vec2.size() < 3)) return;
788 
789  std::vector<art::Ptr<recob::Hit>> hitscl1uniquetpc;
790  std::vector<art::Ptr<recob::Hit>> hitscl2uniquetpc;
791 
792  if (tpc1 == tpc2)
793  for (size_t i = 0; i < vec1.size(); ++i)
794  for (size_t j = 0; j < vec2.size(); ++j)
795  if ((vec1[i]->WireID().TPC == vec2[j]->WireID().TPC) && (tpc1 == vec2[j]->WireID().TPC)) {
796  hitscl1uniquetpc.push_back(vec1[i]);
797  hitscl2uniquetpc.push_back(vec2[j]);
798  }
799 
800  if ((hitscl1uniquetpc.size() > 2) && (hitscl2uniquetpc.size() > 2)) {
801  pma::Track3D* trk =
802  fProjectionMatchingAlg.buildSegment(detProp, hitscl1uniquetpc, hitscl2uniquetpc);
803 
804  // turn the track that front is at vertex - easier to handle associations.
805  if ((trk->back()->Hit2DPtr() == pair[0]->GetFirstHit()) ||
806  (trk->back()->Hit2DPtr() == pair[1]->GetFirstHit()))
807  trk->Flip();
808 
809  IniSeg initrack;
810  initrack.idcl1 = pair[0]->GetIdCl();
811  initrack.view1 = pair[0]->GetFirstHit()->WireID().Plane;
812  initrack.hits1 = hitscl1uniquetpc;
813  initrack.idcl2 = pair[1]->GetIdCl();
814  initrack.view2 = pair[1]->GetFirstHit()->WireID().Plane;
815  initrack.hits2 = hitscl2uniquetpc;
816  initrack.track = trk;
817 
818  fInisegs.push_back(initrack);
819  }
820 }
ems::EMShower3D::fInisegs
std::vector< IniSeg > fInisegs
Definition: EMShower3D_module.cc:112
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:407
pma::ProjectionMatchingAlg::buildSegment
pma::Track3D * buildSegment(const detinfo::DetectorPropertiesData &clockData, const std::vector< art::Ptr< recob::Hit >> &hits_1, const std::vector< art::Ptr< recob::Hit >> &hits_2={}) const
Definition: ProjectionMatchingAlg.cxx:899
ems::EMShower3D::fProjectionMatchingAlg
pma::ProjectionMatchingAlg fProjectionMatchingAlg
Definition: EMShower3D_module.cc:128
pma::Track3D::back
pma::Hit3D const * back() const
Definition: PmaTrack3D.h:103
pma::Hit3D::Hit2DPtr
art::Ptr< recob::Hit > const & Hit2DPtr() const
Definition: PmaHit3D.h:49
pma::Track3D::Flip
bool Flip(const detinfo::DetectorPropertiesData &detProp, std::vector< pma::Track3D * > &allTracks)
Definition: PmaTrack3D.cxx:533
EMShower3D& ems::EMShower3D::operator= ( EMShower3D const &  )
delete
EMShower3D& ems::EMShower3D::operator= ( EMShower3D &&  )
delete
void ems::EMShower3D::produce ( art::Event &  e)
overrideprivatevirtual

Implements art::EDProducer.

Definition at line 311 of file EMShower3D_module.cc.

312 {
313  art::ServiceHandle<geo::Geometry const> geom;
314  fSeltracks.clear();
315  fInisegs.clear();
316  fClusters.clear();
317  fPMA3D.clear();
318  fClustersNotUsed.clear();
319 
320  auto tracks = std::make_unique<std::vector<recob::Track>>();
321  auto vertices = std::make_unique<std::vector<recob::Vertex>>();
322  auto clusters = std::make_unique<std::vector<recob::Cluster>>();
323  auto allsp = std::make_unique<std::vector<recob::SpacePoint>>();
324 
325  auto trk2hit = std::make_unique<art::Assns<recob::Track, recob::Hit>>();
326  auto trk2vtx = std::make_unique<art::Assns<recob::Track, recob::Vertex>>();
327  auto cl2hit = std::make_unique<art::Assns<recob::Cluster, recob::Hit>>();
328  auto trk2cl = std::make_unique<art::Assns<recob::Track, recob::Cluster>>();
329  auto trk2sp = std::make_unique<art::Assns<recob::Track, recob::SpacePoint>>();
330  auto sp2hit = std::make_unique<art::Assns<recob::SpacePoint, recob::Hit>>();
331 
332  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
333  auto const detProp =
334  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(e, clockData);
335 
336  if (e.getByLabel(fCluModuleLabel, fCluListHandle)) {
337  fClustersNotUsed.clear();
338  fInisegs.clear();
339  art::FindManyP<recob::Hit> fb(fCluListHandle, e, fCluModuleLabel);
340 
341  for (size_t id = 0; id < fCluListHandle->size(); id++) {
342  std::vector<art::Ptr<recob::Hit>> hitlist;
343  hitlist = fb.at(id);
344 
345  if (hitlist.size() > 5) fClustersNotUsed.push_back(id);
346  }
347 
348  std::vector<ems::DirOfGamma*> showernviews = CollectShower2D(detProp, e);
349 
350  Link(e, detProp, showernviews);
351 
352  while (fInisegs.size()) {
353  fSeltracks.push_back(fInisegs[0]);
354  fInisegs.erase(fInisegs.begin() + 0);
355  }
356 
357  Reoptimize(detProp);
358 
359  // conversion from pma track to recob::track
360 
361  size_t spStart = 0, spEnd = 0;
362  double sp_pos[3], sp_err[6], vtx_pos[3];
363  for (size_t i = 0; i < 6; i++)
364  sp_err[i] = 1.0;
365 
366  fTrkIndex = 0;
367 
368  for (auto const trk : fSeltracks) {
369  tracks->push_back(ConvertFrom(clockData, detProp, *(trk.track)));
370 
371  vtx_pos[0] = trk.track->front()->Point3D().X();
372  vtx_pos[1] = trk.track->front()->Point3D().Y();
373  vtx_pos[2] = trk.track->front()->Point3D().Z();
374  vertices->emplace_back(vtx_pos, fTrkIndex);
375 
376  ++fTrkIndex;
377 
378  std::vector<art::Ptr<recob::Cluster>> cl2d;
379  cl2d.emplace_back(fCluListHandle, trk.idcl1);
380  cl2d.emplace_back(fCluListHandle, trk.idcl2);
381 
382  std::vector<art::Ptr<recob::Hit>> hits2d;
383  art::PtrVector<recob::Hit> sp_hits;
384 
385  spStart = allsp->size();
386  for (int h = trk.track->size() - 1; h >= 0; h--) {
387  pma::Hit3D* h3d = (*trk.track)[h];
388  hits2d.push_back(h3d->Hit2DPtr());
389 
390  if ((h == 0) || (sp_pos[0] != h3d->Point3D().X()) || (sp_pos[1] != h3d->Point3D().Y()) ||
391  (sp_pos[2] != h3d->Point3D().Z())) {
392  if (sp_hits.size()) // hits assigned to the previous sp
393  {
394  util::CreateAssn(*this, e, *allsp, sp_hits, *sp2hit);
395  sp_hits.clear();
396  }
397  sp_pos[0] = h3d->Point3D().X();
398  sp_pos[1] = h3d->Point3D().Y();
399  sp_pos[2] = h3d->Point3D().Z();
400  allsp->push_back(recob::SpacePoint(sp_pos, sp_err, 1.0));
401  }
402  sp_hits.push_back(h3d->Hit2DPtr());
403  }
404  if (sp_hits.size()) // hits assigned to the last sp
405  {
406  util::CreateAssn(*this, e, *allsp, sp_hits, *sp2hit);
407  }
408  spEnd = allsp->size();
409 
410  if (vertices->size()) {
411  size_t vtx_idx = (size_t)(vertices->size() - 1);
412  util::CreateAssn(*this, e, *tracks, *vertices, *trk2vtx, vtx_idx, vtx_idx + 1);
413  }
414 
415  if (cl2d.size()) { util::CreateAssn(*this, e, *tracks, cl2d, *trk2cl); }
416 
417  if (hits2d.size()) {
418  util::CreateAssn(*this, e, *tracks, *allsp, *trk2sp, spStart, spEnd);
419  util::CreateAssn(*this, e, *tracks, hits2d, *trk2hit);
420  }
421  }
422 
423  fIniIndex = fTrkIndex + 1;
424  for (auto const trk : fPMA3D) {
425  tracks->push_back(ConvertFrom2(clockData, detProp, *(trk.track)));
426 
427  fIniIndex++;
428 
429  std::vector<art::Ptr<recob::Cluster>> cl2d;
430  cl2d.push_back(art::Ptr<recob::Cluster>(fCluListHandle, trk.idcl1));
431  cl2d.push_back(art::Ptr<recob::Cluster>(fCluListHandle, trk.idcl2));
432 
433  std::vector<art::Ptr<recob::Hit>> hits2d;
434  art::PtrVector<recob::Hit> sp_hits;
435 
436  spStart = allsp->size();
437  for (int h = trk.track->size() - 1; h >= 0; h--) {
438  pma::Hit3D* h3d = (*trk.track)[h];
439  hits2d.push_back(h3d->Hit2DPtr());
440 
441  if ((h == 0) || (sp_pos[0] != h3d->Point3D().X()) || (sp_pos[1] != h3d->Point3D().Y()) ||
442  (sp_pos[2] != h3d->Point3D().Z())) {
443  if (sp_hits.size()) // hits assigned to the previous sp
444  {
445  util::CreateAssn(*this, e, *allsp, sp_hits, *sp2hit);
446  sp_hits.clear();
447  }
448  sp_pos[0] = h3d->Point3D().X();
449  sp_pos[1] = h3d->Point3D().Y();
450  sp_pos[2] = h3d->Point3D().Z();
451  allsp->push_back(recob::SpacePoint(sp_pos, sp_err, 1.0));
452  }
453  sp_hits.push_back(h3d->Hit2DPtr());
454  }
455  if (sp_hits.size()) // hits assigned to the last sp
456  {
457  util::CreateAssn(*this, e, *allsp, sp_hits, *sp2hit);
458  }
459  spEnd = allsp->size();
460 
461  if (cl2d.size()) { util::CreateAssn(*this, e, *tracks, cl2d, *trk2cl); }
462 
463  if (hits2d.size()) {
464  util::CreateAssn(*this, e, *tracks, *allsp, *trk2sp, spStart, spEnd);
465  util::CreateAssn(*this, e, *tracks, hits2d, *trk2hit);
466  }
467  }
468 
469  // create cluster from hits, which were an input to find initial part of the cascade.
470  fClIndex = 0;
471  for (auto const& cl : fClusters)
472  if (cl.size()) {
473  clusters->push_back(ConvertFrom(cl));
474  fClIndex++;
475 
476  util::CreateAssn(*this, e, *clusters, cl, *cl2hit);
477  }
478 
479  for (unsigned int i = 0; i < showernviews.size(); i++)
480  delete showernviews[i];
481 
482  for (unsigned int i = 0; i < fSeltracks.size(); i++)
483  delete fSeltracks[i].track;
484 
485  for (unsigned int i = 0; i < fInisegs.size(); i++)
486  delete fInisegs[i].track;
487 
488  for (unsigned int i = 0; i < fPMA3D.size(); i++)
489  delete fPMA3D[i].track;
490  }
491 
492  e.put(std::move(tracks));
493  e.put(std::move(vertices));
494  e.put(std::move(clusters));
495  e.put(std::move(allsp));
496 
497  e.put(std::move(trk2hit));
498  e.put(std::move(trk2vtx));
499  e.put(std::move(cl2hit));
500  e.put(std::move(trk2cl));
501  e.put(std::move(trk2sp));
502  e.put(std::move(sp2hit));
503 }
ems::EMShower3D::fClusters
std::vector< std::vector< art::Ptr< recob::Hit > > > fClusters
Definition: EMShower3D_module.cc:116
ems::EMShower3D::fCluListHandle
art::Handle< std::vector< recob::Cluster > > fCluListHandle
Definition: EMShower3D_module.cc:131
pma::Hit3D::Point3D
TVector3 const & Point3D() const
Definition: PmaHit3D.h:55
ems::EMShower3D::fIniIndex
unsigned int fIniIndex
Definition: EMShower3D_module.cc:123
cl
QAsciiDict< Entry > cl
Definition: tclscanner.cpp:1148
ems::EMShower3D::fClIndex
unsigned int fClIndex
Definition: EMShower3D_module.cc:122
ems::EMShower3D::fInisegs
std::vector< IniSeg > fInisegs
Definition: EMShower3D_module.cc:112
reco_momentum_tuples.h
h
Definition: reco_momentum_tuples.py:66
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
Definition: DataViewImpl.h:633
ems::EMShower3D::fTrkIndex
unsigned int fTrkIndex
Definition: EMShower3D_module.cc:121
art::PtrVector::push_back
void push_back(Ptr< U > const &p)
Definition: PtrVector.h:435
wirecell.gen.depos.move
def move(depos, offset)
Definition: depos.py:107
QAsciiDict::size
uint size() const
Definition: qasciidict.h:56
art::DataViewImpl::put
ProductID put(std::unique_ptr< PROD > &&edp, std::string const &instance={})
Definition: DataViewImpl.h:686
util::CreateAssn
bool CreateAssn(PRODUCER const &prod, art::Event &evt, std::vector< T > const &a, art::Ptr< U > const &b, art::Assns< U, T > &assn, std::string a_instance, size_t indx=UINT_MAX)
Creates a single one-to-one association.
Definition: AssociationUtil.h:601
ems::EMShower3D::CollectShower2D
std::vector< ems::DirOfGamma * > CollectShower2D(detinfo::DetectorPropertiesData const &detProp, art::Event const &e)
Definition: EMShower3D_module.cc:590
art::PtrVector::size
size_type size() const
Definition: PtrVector.h:302
tracks
Definition: tracks.py:1
ems::EMShower3D::fPMA3D
std::vector< IniSeg > fPMA3D
Definition: EMShower3D_module.cc:114
ems::EMShower3D::ConvertFrom
recob::Track ConvertFrom(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, pma::Track3D &src)
Definition: EMShower3D_module.cc:184
ems::EMShower3D::Link
void Link(art::Event const &e, detinfo::DetectorPropertiesData const &detProp, std::vector< ems::DirOfGamma * > input)
Definition: EMShower3D_module.cc:620
ems::EMShower3D::ConvertFrom2
recob::Track ConvertFrom2(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, pma::Track3D &src)
Definition: EMShower3D_module.cc:247
ems::EMShower3D::fClustersNotUsed
std::vector< size_t > fClustersNotUsed
Definition: EMShower3D_module.cc:118
pma::Hit3D::Hit2DPtr
art::Ptr< recob::Hit > const & Hit2DPtr() const
Definition: PmaHit3D.h:49
ems::EMShower3D::fSeltracks
std::vector< IniSeg > fSeltracks
Definition: EMShower3D_module.cc:113
ems::EMShower3D::Reoptimize
void Reoptimize(detinfo::DetectorPropertiesData const &detProp)
Definition: EMShower3D_module.cc:506
ems::EMShower3D::fCluModuleLabel
std::string fCluModuleLabel
Definition: EMShower3D_module.cc:125
art::PtrVector::clear
void clear()
Definition: PtrVector.h:533
void ems::EMShower3D::Reoptimize ( detinfo::DetectorPropertiesData const &  detProp)
private

Definition at line 506 of file EMShower3D_module.cc.

507 {
508  if (empty(fSeltracks)) return;
509  const float min_dist = 3.0F;
510  size_t ta = 0;
511  while (ta < (fSeltracks.size() - 1)) {
512  size_t tb = ta + 1;
513  bool found = false;
514  while (tb < fSeltracks.size()) {
515  if (ta == tb) {
516  tb++;
517  continue;
518  }
519 
520  TVector3 p1 = fSeltracks[ta].track->front()->Point3D();
521  TVector3 p2 = fSeltracks[tb].track->front()->Point3D();
522  float dist = std::sqrt(pma::Dist2(p1, p2));
523 
524  if (dist < min_dist)
525  if ((fSeltracks[ta].idcl1 == fSeltracks[tb].idcl1) ||
526  (fSeltracks[ta].idcl1 == fSeltracks[tb].idcl2) ||
527  (fSeltracks[ta].idcl2 == fSeltracks[tb].idcl1) ||
528  (fSeltracks[ta].idcl2 == fSeltracks[tb].idcl2)) {
529  found = true;
530  size_t view3 = fSeltracks[ta].view1;
531  size_t idcl3 = fSeltracks[ta].idcl1;
532  std::vector<art::Ptr<recob::Hit>> hits3 = fSeltracks[ta].hits1;
533  std::vector<art::Ptr<recob::Hit>> hits = fSeltracks[ta].hits1;
534  for (size_t h = 0; h < fSeltracks[ta].hits2.size(); ++h)
535  hits.push_back(fSeltracks[ta].hits2[h]);
536 
537  if ((fSeltracks[tb].view1 != fSeltracks[ta].view1) &&
538  (fSeltracks[tb].view1 != fSeltracks[ta].view2)) {
539  view3 = fSeltracks[tb].view1;
540  for (size_t h = 0; h < fSeltracks[tb].hits1.size(); ++h)
541  hits.push_back(fSeltracks[tb].hits1[h]);
542  }
543  if ((fSeltracks[tb].view2 != fSeltracks[ta].view1) &&
544  (fSeltracks[tb].view2 != fSeltracks[ta].view2)) {
545  view3 = fSeltracks[tb].view2;
546  for (size_t h = 0; h < fSeltracks[tb].hits2.size(); ++h)
547  hits.push_back(fSeltracks[tb].hits2[h]);
548  }
549 
550  if ((view3 == fSeltracks[ta].view1) || (view3 == fSeltracks[ta].view2)) {
551  delete fSeltracks[ta].track;
552  fSeltracks.erase(fSeltracks.begin() + ta);
553  }
554  else {
555  pma::Track3D* track = fProjectionMatchingAlg.buildSegment(detProp, hits);
556 
557  if (pma::Dist2(track->back()->Point3D(), fSeltracks[ta].track->front()->Point3D()) <
558  pma::Dist2(track->front()->Point3D(), fSeltracks[ta].track->front()->Point3D()))
559  track->Flip();
560 
561  IniSeg initrack;
562  initrack.idcl1 = fSeltracks[ta].idcl1;
563  initrack.idcl3 = idcl3;
564  initrack.view1 = fSeltracks[ta].view1;
565  initrack.view3 = view3;
566  initrack.hits1 = fSeltracks[ta].hits1;
567  initrack.hits3 = hits3;
568  initrack.idcl2 = fSeltracks[ta].idcl2;
569  initrack.view2 = fSeltracks[ta].view2;
570  initrack.hits2 = fSeltracks[ta].hits2;
571  initrack.track = track;
572 
573  delete fSeltracks[tb].track;
574  delete fSeltracks[ta].track;
575  fSeltracks.erase(fSeltracks.begin() + tb);
576  fSeltracks.erase(fSeltracks.begin() + ta);
577  fSeltracks.push_back(initrack);
578  }
579  }
580 
581  if (found) break;
582  tb++;
583  }
584 
585  if (!found) ta++;
586  }
587 }
pma::Dist2
double Dist2(const TVector2 &v1, const TVector2 &v2)
Definition: Utilities.cxx:37
pma::Track3D::front
pma::Hit3D const * front() const
Definition: PmaTrack3D.h:98
pma::Hit3D::Point3D
TVector3 const & Point3D() const
Definition: PmaHit3D.h:55
reco_momentum_tuples.h
h
Definition: reco_momentum_tuples.py:66
pma::ProjectionMatchingAlg::buildSegment
pma::Track3D * buildSegment(const detinfo::DetectorPropertiesData &clockData, const std::vector< art::Ptr< recob::Hit >> &hits_1, const std::vector< art::Ptr< recob::Hit >> &hits_2={}) const
Definition: ProjectionMatchingAlg.cxx:899
larg4::dist
constexpr double dist(const TReal *x, const TReal *y, const unsigned int dimension)
Definition: OpFastScintillation.hh:563
ems::EMShower3D::fProjectionMatchingAlg
pma::ProjectionMatchingAlg fProjectionMatchingAlg
Definition: EMShower3D_module.cc:128
pma::Track3D::back
pma::Hit3D const * back() const
Definition: PmaTrack3D.h:103
ems::EMShower3D::fSeltracks
std::vector< IniSeg > fSeltracks
Definition: EMShower3D_module.cc:113
pma::Track3D::Flip
bool Flip(const detinfo::DetectorPropertiesData &detProp, std::vector< pma::Track3D * > &allTracks)
Definition: PmaTrack3D.cxx:533
util::empty
decltype(auto) constexpr empty(T &&obj)
ADL-aware version of std::empty.
Definition: StdUtils.h:97
bool ems::EMShower3D::Validate ( detinfo::DetectorPropertiesData const &  detProp,
std::vector< ems::DirOfGamma * >  input,
size_t  id1,
size_t  id2,
size_t  c1,
size_t  c2,
size_t  plane3 
)
private

Definition at line 823 of file EMShower3D_module.cc.

830 {
831  bool result = false;
832  if (id1 == id2) return false;
833 
834  std::vector<art::Ptr<recob::Hit>> vec1 = input[id1]->GetCandidates()[c1].GetIniHits();
835  std::vector<art::Ptr<recob::Hit>> vec2 = input[id2]->GetCandidates()[c2].GetIniHits();
836 
837  if ((vec1.size() < 3) || (vec2.size() < 3)) return false;
838 
839  std::vector<art::Ptr<recob::Hit>> hitscl1uniquetpc;
840  std::vector<art::Ptr<recob::Hit>> hitscl2uniquetpc;
841 
842  size_t tpc = vec1[0]->WireID().TPC;
843  for (size_t i = 0; i < vec1.size(); ++i)
844  for (size_t j = 0; j < vec2.size(); ++j)
845  if ((vec1[i]->WireID().TPC == tpc) && (vec2[j]->WireID().TPC == tpc)) {
846  hitscl1uniquetpc.push_back(vec1[i]);
847  hitscl2uniquetpc.push_back(vec2[j]);
848  }
849 
850  if ((hitscl1uniquetpc.size() < 3) || (hitscl2uniquetpc.size() < 3)) return false;
851 
852  pma::Track3D* track =
853  fProjectionMatchingAlg.buildSegment(detProp, hitscl1uniquetpc, hitscl2uniquetpc);
854  for (size_t i = 0; i < input.size(); ++i) {
855  std::vector<Hit2D*> hits2dcl = input[i]->GetHits2D();
856  for (size_t h = 0; h < hits2dcl.size(); ++h) {
857  TVector2 pfront = pma::GetProjectionToPlane(
858  track->front()->Point3D(), plane3, track->FrontTPC(), track->FrontCryo());
859  TVector2 pback = pma::GetProjectionToPlane(
860  track->back()->Point3D(), plane3, track->BackTPC(), track->BackCryo());
861  if ((pma::Dist2(hits2dcl[h]->GetPointCm(), pfront) < 1.0F) &&
862  (pma::Dist2(hits2dcl[h]->GetPointCm(), pback) < 1.0F)) {
863  result = true;
864  break;
865  }
866  }
867  }
868  delete track;
869  return result;
870 }
result
static QCString result
Definition: fortranscanner.cpp:56614
pma::Dist2
double Dist2(const TVector2 &v1, const TVector2 &v2)
Definition: Utilities.cxx:37
pma::Track3D::front
pma::Hit3D const * front() const
Definition: PmaTrack3D.h:98
pma::Hit3D::Point3D
TVector3 const & Point3D() const
Definition: PmaHit3D.h:55
pma::Track3D::BackTPC
unsigned int BackTPC() const
Definition: PmaTrack3D.h:156
pma::Track3D::BackCryo
unsigned int BackCryo() const
Definition: PmaTrack3D.h:161
reco_momentum_tuples.h
h
Definition: reco_momentum_tuples.py:66
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:407
pma::Track3D::FrontTPC
unsigned int FrontTPC() const
Definition: PmaTrack3D.h:145
pma::Track3D::FrontCryo
unsigned int FrontCryo() const
Definition: PmaTrack3D.h:150
pma::ProjectionMatchingAlg::buildSegment
pma::Track3D * buildSegment(const detinfo::DetectorPropertiesData &clockData, const std::vector< art::Ptr< recob::Hit >> &hits_1, const std::vector< art::Ptr< recob::Hit >> &hits_2={}) const
Definition: ProjectionMatchingAlg.cxx:899
ems::EMShower3D::fProjectionMatchingAlg
pma::ProjectionMatchingAlg fProjectionMatchingAlg
Definition: EMShower3D_module.cc:128
pma::GetProjectionToPlane
TVector2 GetProjectionToPlane(const TVector3 &p, unsigned int plane, unsigned int tpc, unsigned int cryo)
Definition: Utilities.cxx:270
pma::Track3D::back
pma::Hit3D const * back() const
Definition: PmaTrack3D.h:103

Member Data Documentation

calo::CalorimetryAlg ems::EMShower3D::fCalorimetryAlg
private

Definition at line 129 of file EMShower3D_module.cc.

unsigned int ems::EMShower3D::fClIndex
private

Definition at line 122 of file EMShower3D_module.cc.

art::Handle<std::vector<recob::Cluster> > ems::EMShower3D::fCluListHandle
private

Definition at line 131 of file EMShower3D_module.cc.

std::string ems::EMShower3D::fCluModuleLabel
private

Definition at line 125 of file EMShower3D_module.cc.

std::vector<std::vector<art::Ptr<recob::Hit> > > ems::EMShower3D::fClusters
private

Definition at line 116 of file EMShower3D_module.cc.

std::vector<size_t> ems::EMShower3D::fClustersNotUsed
private

Definition at line 118 of file EMShower3D_module.cc.

unsigned int ems::EMShower3D::fIniIndex
private

Definition at line 123 of file EMShower3D_module.cc.

std::vector<IniSeg> ems::EMShower3D::fInisegs
private

Definition at line 112 of file EMShower3D_module.cc.

std::vector<IniSeg> ems::EMShower3D::fPMA3D
private

Definition at line 114 of file EMShower3D_module.cc.

pma::ProjectionMatchingAlg ems::EMShower3D::fProjectionMatchingAlg
private

Definition at line 128 of file EMShower3D_module.cc.

std::vector<IniSeg> ems::EMShower3D::fSeltracks
private

Definition at line 113 of file EMShower3D_module.cc.

std::vector<size_t> ems::EMShower3D::fTracksNotUsed
private

Definition at line 119 of file EMShower3D_module.cc.

std::string ems::EMShower3D::fTrk3DModuleLabel
private

Definition at line 126 of file EMShower3D_module.cc.

unsigned int ems::EMShower3D::fTrkIndex
private

Definition at line 121 of file EMShower3D_module.cc.

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