Producer algorithm for PixelMap, input to CVN neural net. More...

#include "PixelMapProducer.h"

Public Member Functions
	PixelMapProducer (unsigned int nWire, unsigned int nTdc, double tRes, double threshold=0.)

	PixelMapProducer ()=default

	PixelMapProducer (const fhicl::ParameterSet &pset)

unsigned int	NWire () const

unsigned int	NTdc () const

double	TRes () const

unsigned int	TotHits () const

double	Threshold () const

bool	MultipleDrifts () const

void	SetTotHits (unsigned int tothits)

void	SetMultipleDrifts (bool multipledrifts=true)

virtual Boundary	DefineBoundary (detinfo::DetectorPropertiesData const &detProp, const std::vector< const T * > &cluster)
	Get boundaries for pixel map representation of cluster. More...

virtual void	ConvertLocaltoGlobal (geo::WireID wireid, unsigned int &globalWire, unsigned int &globalPlane) const

virtual void	ConvertLocaltoGlobalTDC (geo::WireID wireid, double localTDC, unsigned int &globalWire, unsigned int &globalPlane, double &globalTDC) const

virtual PixelMap	CreateMap (detinfo::DetectorPropertiesData const &detProp, const std::vector< art::Ptr< T >> &cluster)

virtual PixelMap	CreateMap (detinfo::DetectorPropertiesData const &detProp, const std::vector< const T * > &cluster)

virtual PixelMap	CreateMapGivenBoundary (detinfo::DetectorPropertiesData const &detProp, const std::vector< const T * > &cluster, const Boundary &bound)

Private Attributes
unsigned int	fNWire
	Number of wires, length for pixel maps. More...

unsigned int	fNTdc
	Number of tdcs, width of pixel map. More...

double	fTRes
	Timing resolution for pixel map. More...

unsigned int	fTotHits
	Total hits in the pixel map. More...

double	fThreshold
	Charge threshold to consider for hits/waveforms etc. More...

bool	fMultipleDrifts
	True if making the pixel map requires handling for multiple drift regions. More...

Detailed Description

template<class T, class U>
class lcvn::PixelMapProducer< T, U >

Producer algorithm for PixelMap, input to CVN neural net.

Definition at line 81 of file PixelMapProducer.h.

Constructor & Destructor Documentation

template<class T , class U >

lcvn::PixelMapProducer< T, U >::PixelMapProducer	(	unsigned int	nWire,
		unsigned int	nTdc,
		double	tRes,
		double	threshold = `0.`
	)

Definition at line 102 of file PixelMapProducer.cxx.

106 : fNWire(nWire), fNTdc(nTdc), fTRes(tRes), fThreshold(threshold), fMultipleDrifts(false)

107 {}

lcvn::PixelMapProducer::fThreshold

double fThreshold

Charge threshold to consider for hits/waveforms etc.

Definition: PixelMapProducer.h:127

lcvn::PixelMapProducer::fNTdc

unsigned int fNTdc

Number of tdcs, width of pixel map.

Definition: PixelMapProducer.h:124

lcvn::PixelMapProducer::fNWire

unsigned int fNWire

Number of wires, length for pixel maps.

Definition: PixelMapProducer.h:123

lcvn::PixelMapProducer::fTRes

double fTRes

Timing resolution for pixel map.

Definition: PixelMapProducer.h:125

lcvn::PixelMapProducer::fMultipleDrifts

bool fMultipleDrifts

True if making the pixel map requires handling for multiple drift regions.

Definition: PixelMapProducer.h:129

template<class T , class U >

lcvn::PixelMapProducer< T, U >::PixelMapProducer ( )

default

template<class T , class U >

lcvn::PixelMapProducer< T, U >::PixelMapProducer ( const fhicl::ParameterSet & pset )

Definition at line 110 of file PixelMapProducer.cxx.

     : fNWire(pset.get<unsigned int>("WireLength"))
     , fNTdc(pset.get<unsigned int>("TdcWidth"))
     , fTRes(pset.get<double>("TimeResolution"))
     , fThreshold(pset.get<double>("Threshold"))
     , fMultipleDrifts(pset.get<bool>("MultipleDrifts"))
   {}

Member Function Documentation

template<class T , class U >

void lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobal	(	geo::WireID	wireid,
		unsigned int &	globalWire,
		unsigned int &	globalPlane
	)		const

virtual

Definition at line 282 of file PixelMapProducer.cxx.

References geo::PlaneID::Plane, and geo::WireID::Wire.

Referenced by lcvn::PixelMapProducer< T, U >::CreateMapGivenBoundary(), and lcvn::PixelMapProducer< T, U >::DefineBoundary().

   {
     globalWire = wireid.Wire;
     globalPlane = wireid.Plane;
   }

template<class T , class U >

void lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobalTDC	(	geo::WireID	wireid,
		double	localTDC,
		unsigned int &	globalWire,
		unsigned int &	globalPlane,
		double &	globalTDC
	)		const

virtual

Definition at line 291 of file PixelMapProducer.cxx.

References geo::PlaneID::Plane, and geo::WireID::Wire.

Referenced by lcvn::PixelMapProducer< T, U >::CreateMapGivenBoundary(), and lcvn::PixelMapProducer< T, U >::DefineBoundary().

   {
     globalWire = wireid.Wire;
     globalPlane = wireid.Plane;
     globalTDC = localTDC;
   }

template<class T , class U >

PixelMap lcvn::PixelMapProducer< T, U >::CreateMap	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< art::Ptr< T >> &	cluster
	)

virtual

Definition at line 119 of file PixelMapProducer.cxx.

   {
     std::vector<const T*> newCluster;
     for (const art::Ptr<T> hit : cluster) {
       newCluster.push_back(hit.get());
     }
     return CreateMap(detProp, newCluster);
   }

template<class T , class U >

PixelMap lcvn::PixelMapProducer< T, U >::CreateMap	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< const T * > &	cluster
	)

virtual

Definition at line 130 of file PixelMapProducer.cxx.

References lcvn::PixelMapProducer< T, U >::CreateMapGivenBoundary(), and lcvn::PixelMapProducer< T, U >::DefineBoundary().

   {
     Boundary bound = DefineBoundary(detProp, cluster);
     return CreateMapGivenBoundary(detProp, cluster, bound);
   }

template<class T , class U >

PixelMap lcvn::PixelMapProducer< T, U >::CreateMapGivenBoundary	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< const T * > &	cluster,
		const Boundary &	bound
	)

virtual

Definition at line 138 of file PixelMapProducer.cxx.

References lcvn::PixelMap::Add(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobal(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobalTDC(), geo::WireID::first(), lcvn::PixelMapProducer< T, U >::fMultipleDrifts, lcvn::PixelMapProducer< T, U >::fNTdc, lcvn::PixelMapProducer< T, U >::fNWire, lcvn::PixelMapProducer< T, U >::fThreshold, lcvn::PixelMapProducer< T, U >::fTotHits, geo::PlaneID::Plane, lcvn::PixelMap::SetTotHits(), and geo::WireID::Wire.

Referenced by lcvn::PixelMapProducer< T, U >::CreateMap().

   {
     PixelMap pm(fNWire, fNTdc, bound);
 
     for (size_t iHit = 0; iHit < cluster.size(); ++iHit) {
 
       U wraphit(*(cluster[iHit]), fThreshold);
       Waveform wf = wraphit.GetWaveform();
       geo::WireID wireid = wraphit.GetID();
 
       unsigned int tempWire = wireid.Wire;
       unsigned int tempPlane = wireid.Plane;
 
       if (!fMultipleDrifts) ConvertLocaltoGlobal(wireid, tempWire, tempPlane);
 
       for (auto& pulse : wf) {
         // Leigh: Simple modification to unwrap the collection view wire plane
         for (auto& i : pulse) {
           const double pe = i.second;
           double temptdc = i.first;
           if (fMultipleDrifts)
             ConvertLocaltoGlobalTDC(wireid, i.first, tempWire, tempPlane, temptdc);
 
           const unsigned int wire = tempWire;
           const unsigned int wirePlane = tempPlane;
           const double tdc = temptdc;
 
           pm.Add(wire, tdc, wirePlane, pe);
         }
       }
     }
     pm.SetTotHits(fTotHits);
     return pm;
   }

template<class T , class U >

Boundary lcvn::PixelMapProducer< T, U >::DefineBoundary	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< const T * > &	cluster
	)

virtual

Get boundaries for pixel map representation of cluster.

Definition at line 184 of file PixelMapProducer.cxx.

References util::abs(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobal(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobalTDC(), geo::WireID::first(), lcvn::PixelMapProducer< T, U >::fMultipleDrifts, lcvn::PixelMapProducer< T, U >::fNWire, lcvn::PixelMapProducer< T, U >::fThreshold, lcvn::PixelMapProducer< T, U >::fTotHits, lcvn::PixelMapProducer< T, U >::fTRes, geo::PlaneID::Plane, and geo::WireID::Wire.

Referenced by lcvn::PixelMapProducer< T, U >::CreateMap().

   {
     std::vector<double> tmin_0;
     std::vector<double> tmin_1;
     std::vector<double> tmin_2;
 
     std::vector<int> wire_0, bwire_0;
     std::vector<int> wire_1, bwire_1;
     std::vector<int> wire_2, bwire_2;
 
     std::vector<double> tsum = {0., 0., 0.};
     std::vector<double> tsize = {0., 0., 0.};
 
     for (size_t iHit = 0; iHit < cluster.size(); ++iHit) {
       U wraphit(*(cluster[iHit]), fThreshold);
       Waveform wf = wraphit.GetWaveform();
       geo::WireID wireid = wraphit.GetID();
 
       unsigned int tempWire = wireid.Wire;
       unsigned int tempPlane = wireid.Plane;
 
       if (!fMultipleDrifts) ConvertLocaltoGlobal(wireid, tempWire, tempPlane);
 
       for (auto& pulse : wf) {
         double min_tick = (double)INT_MAX;
         for (auto& i : pulse) {
           double temptdc = i.first;
 
           if (fMultipleDrifts)
             ConvertLocaltoGlobalTDC(wireid, i.first, tempWire, tempPlane, temptdc);
 
           if (temptdc < min_tick) min_tick = temptdc;
 
           tsum[tempPlane] += temptdc;
           tsize[tempPlane] += 1.;
         }
 
         if (!(pulse.empty())) {
           if (tempPlane == 0) {
             tmin_0.push_back(min_tick);
             wire_0.push_back(tempWire);
           }
           if (tempPlane == 1) {
             tmin_1.push_back(min_tick);
             wire_0.push_back(tempWire);
           }
           if (tempPlane == 2) {
             tmin_2.push_back(min_tick);
             wire_0.push_back(tempWire);
           }
         }
       } // end loop over pulses on single wire
     }   // end loop over struck wires
 
     double tmean_0 = tsum[0] / tsize[0];
     double tmean_1 = tsum[1] / tsize[1];
     double tmean_2 = tsum[2] / tsize[2];
 
     for (int i = 0; i < (int)wire_0.size(); i++) {
       if (std::abs(tmin_0[i] - tmean_0) < (double)fTRes) bwire_0.push_back(wire_0[i]);
     }
     for (int i = 0; i < (int)wire_1.size(); i++) {
       if (std::abs(tmin_1[i] - tmean_1) < (double)fTRes) bwire_1.push_back(wire_1[i]);
     }
     for (int i = 0; i < (int)wire_2.size(); i++) {
       if (std::abs(tmin_2[i] - tmean_2) < (double)fTRes) bwire_2.push_back(wire_2[i]);
     }
 
     std::cout << "Boundary wire vector sizes: " << bwire_0.size() << ", " << bwire_1.size() << ", "
               << bwire_2.size() << std::endl;
 
     int minwire_0 = 0;
     int minwire_1 = 0;
     int minwire_2 = 0;
     auto minwireelement_0 = std::min_element(bwire_0.begin(), bwire_0.end());
     auto minwireelement_1 = std::min_element(bwire_1.begin(), bwire_1.end());
     auto minwireelement_2 = std::min_element(bwire_2.begin(), bwire_2.end());
 
     if (bwire_0.size() > 0) {
       minwire_0 = *minwireelement_0 - 1;
       std::cout << "minwire 0: " << (*minwireelement_0 - 1) << std::endl;
     }
     if (bwire_1.size() > 0) {
       minwire_1 = *minwireelement_1 - 1;
       std::cout << "minwire 1: " << (*minwireelement_1 - 1) << std::endl;
     }
     if (bwire_2.size() > 0) {
       minwire_2 = *minwireelement_2 - 1;
       std::cout << "minwire 2: " << (*minwireelement_2 - 1) << std::endl;
     }
 
     fTotHits = bwire_0.size() + bwire_1.size() + bwire_2.size();
 
     return Boundary(fNWire, fTRes, minwire_0, minwire_1, minwire_2, tmean_0, tmean_1, tmean_2);
   }