HitHandScan Class Reference

Inheritance diagram for HitHandScan:

Public Member Functions
	HitHandScan (fhicl::ParameterSet const &p)
	HitHandScan (HitHandScan const &)=delete
	HitHandScan (HitHandScan &&)=delete
HitHandScan &	operator= (HitHandScan const &)=delete
HitHandScan &	operator= (HitHandScan &&)=delete
void	produce (art::Event &e) override
void	reconfigure (fhicl::ParameterSet const &p)
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
bool	insidePolygon (std::vector< Point > polygon, Point ptest)

Private Attributes
std::string	fPreviousHitModuleLabel
std::string	fHandScanFileName

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 49 of file HitHandScan_module.cc.

Constructor & Destructor Documentation

HitHandScan::HitHandScan ( fhicl::ParameterSet const &  p )

explicit

Definition at line 76 of file HitHandScan_module.cc.

                                                    : EDProducer{p}
{
  this->reconfigure(p);
  recob::HitCollectionCreator::declare_products(producesCollector());
}

HitHandScan::HitHandScan ( HitHandScan const &  )

delete

HitHandScan::HitHandScan ( HitHandScan &&  )

delete

Member Function Documentation

bool HitHandScan::insidePolygon ( std::vector< Point >  polygon, Point  ptest  )

private

Definition at line 191 of file HitHandScan_module.cc.

{
  size_t n = polygon.size();
  size_t i, j;
  bool c = false;
  for (i = 0, j = n - 1; i < n; j = i++)
    {
      if ( ((polygon[i].tick > ptest.tick) != (polygon[j].tick > ptest.tick)) &&
           (ptest.wire < (polygon[j].wire - polygon[i].wire) * 
            (ptest.tick - polygon[i].tick) / (polygon[j].tick - polygon[i].tick) + polygon[i].wire) )
        c = !c;
    }
  return c;
}

HitHandScan& HitHandScan::operator= ( HitHandScan const &  )

delete

HitHandScan& HitHandScan::operator= ( HitHandScan &&  )

delete

void HitHandScan::produce ( art::Event &  e )

overridevirtual

Implements art::EDProducer.

Definition at line 82 of file HitHandScan_module.cc.

{
  auto prevHitHandle = e.getHandle< std::vector< recob::Hit> >(fPreviousHitModuleLabel);

  art::FindOneP<recob::Wire> wires(prevHitHandle,e,fPreviousHitModuleLabel);
  art::FindOneP<raw::RawDigit> rawdigits(prevHitHandle,e,fPreviousHitModuleLabel);

  recob::HitCollectionCreator hcol(e, wires.isValid(), rawdigits.isValid());

  std::vector<Point> trackBox;
  
  std::string readrun,readevent,readtpc,readplane,readwire,readtick,readnumpoly;
  unsigned int run,event,tpc,plane,numpoly;
  double tick,wire;
  int linenumber = 0;
  std::string line;
  std::ifstream handscanfile(fHandScanFileName.c_str(),std::ifstream::in);
  if (handscanfile.is_open())
    {
      while (std::getline(handscanfile,line))
        {
          linenumber++;
          std::stringstream ss(line);
          ss >> std::ws >> readrun >> readevent >> readtpc >> readplane >> readwire >> readtick >> readnumpoly;
          try 
            {
              run = stoul(readrun);
              event = stoul(readevent);
              tick = stod(readtick);
              tpc = stoul(readtpc);
              plane = stoul(readplane);
              wire = stod(readwire);
              numpoly = stoul(readnumpoly);
              
              Point p;
              p.run = run;
              p.event = event;
              p.tpc = tpc;
              p.plane = plane;
              p.wire = wire;
              p.tick = tick;
              p.numpoly = numpoly;
              trackBox.push_back(p);
            } 
          catch (const std::exception& ex)
            {
              std::cout << "Line # " << linenumber << " is invalid, skipping..." << std::endl;
            }
        }
      handscanfile.close();
    }
  else
    {
      std::cout << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n"
                << "NO FILE FOUND WITH HAND SCANNED HIT DATA\n"
                << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!" << std::endl;
      return;
    }

  for (size_t i_hit = 0; i_hit < prevHitHandle->size(); i_hit++)
    {
      art::Ptr<recob::Hit> phit(prevHitHandle,i_hit);

      art::Ptr<recob::Wire> wire(wires.at(i_hit));
      art::Ptr<raw::RawDigit> rawdigit(rawdigits.at(i_hit));

      Point ptest;
      ptest.run = e.run();
      ptest.event = e.event();
      ptest.tpc = phit->WireID().TPC;
      ptest.plane = phit->WireID().Plane;
      ptest.wire = phit->WireID().Wire;
      ptest.tick = phit->PeakTime();
            
      std::map< unsigned int, std::vector<Point> > polygons;
      for (auto const &poly : trackBox)
        {
          if (poly.run == ptest.run 
              && poly.event == ptest.event
              && poly.tpc == ptest.tpc
              && poly.plane == ptest.plane)
            {
              polygons[poly.numpoly].push_back(poly);
            }
        }
      bool acceptPoint = false;
      for (auto const &p : polygons)
        {
          if (insidePolygon(p.second,ptest)) 
            {
              acceptPoint = true;
            }
        }
      if (acceptPoint) 
        {
          hcol.emplace_back(*phit,wire,rawdigit);
        }
    }
  
  hcol.put_into(e);

}

void HitHandScan::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 185 of file HitHandScan_module.cc.

{
  fHandScanFileName = p.get<std::string>("HandScanFileName");
  fPreviousHitModuleLabel = p.get<std::string>("PreviousHitModuleLabel");
}

Member Data Documentation

std::string HitHandScan::fHandScanFileName

private

Definition at line 71 of file HitHandScan_module.cc.

std::string HitHandScan::fPreviousHitModuleLabel

private

Definition at line 70 of file HitHandScan_module.cc.

---

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

• dunereco/dunereco/HitFinderDUNE/HitHandScan_module.cc