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 ()

	PixelMapProducer (const fhicl::ParameterSet &pset)

void	SetMultipleDrifts ()

unsigned int	NROI ()

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

unsigned int	NWire () const

unsigned int	NTdc () const

double	TRes () 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)

Protected 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...

geo::GeometryCore const *	fGeometry

Detailed Description

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

Producer algorithm for PixelMap, input to CVN neural net.

Definition at line 86 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 107 of file PixelMapProducer.cxx.

References lcvn::PixelMapProducer< T, U >::fGeometry.

     : fNWire(nWire), fNTdc(nTdc), fTRes(tRes), fThreshold(threshold), fMultipleDrifts(false)
   {
 
     fGeometry = &*(art::ServiceHandle<geo::Geometry>());
   }

template<class T , class U >

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

Definition at line 118 of file PixelMapProducer.cxx.

References lcvn::PixelMapProducer< T, U >::fGeometry.

   {
     fGeometry = &*(art::ServiceHandle<geo::Geometry>());
   }

template<class T , class U >

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

Definition at line 124 of file PixelMapProducer.cxx.

References lcvn::PixelMapProducer< T, U >::fGeometry.

     : 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"))
   {
     fGeometry = &*(art::ServiceHandle<geo::Geometry>());
   }

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 302 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 311 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 135 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 146 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 154 of file PixelMapProducer.cxx.

References lcvn::PixelMap::Add(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobal(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobalTDC(), 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 201 of file PixelMapProducer.cxx.

References util::abs(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobal(), lcvn::PixelMapProducer< T, U >::ConvertLocaltoGlobalTDC(), 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();
 
     Boundary bound(fNWire, fTRes, minwire_0, minwire_1, minwire_2, tmean_0, tmean_1, tmean_2);
 
     return bound;
   }