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

#include <PixelMapSimProducer.h>

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

	PixelMapSimProducer ()

void	SetUnwrapped (unsigned short unwrap)

void	SetProtoDUNE ()

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

unsigned int	NROI ()

void	GetDUNEGlobalWire (unsigned int localWire, unsigned int plane, unsigned int tpc, unsigned int &globalWire, unsigned int &globalPlane) const
	Function to convert to a global unwrapped wire number. More...

void	GetDUNEGlobalWireTDC (detinfo::DetectorPropertiesData const &detProp, unsigned int localWire, double localTDC, unsigned int plane, unsigned int tpc, unsigned int &globalWire, unsigned int &globalPlane, double &globalTDC) const

void	GetDUNE10ktGlobalWireTDC (detinfo::DetectorPropertiesData const &detProp, unsigned int localWire, double localTDC, unsigned int plane, unsigned int tpc, unsigned int &globalWire, unsigned int &globalPlane, double &globalTDC) const

void	GetProtoDUNEGlobalWire (unsigned int localWire, unsigned int plane, unsigned int tpc, unsigned int &globalWire, unsigned int &globalPlane) const

void	GetProtoDUNEGlobalWireTDC (unsigned int localWire, double localTDC, unsigned int plane, unsigned int tpc, unsigned int &globalWire, double &globalTDC, unsigned int &globalPlane) const

void	GetDUNEVertDrift3ViewGlobalWire (unsigned int localWire, unsigned int plane, unsigned int tpc, unsigned int &globalWire, unsigned int &globalPlane) const

unsigned int	NWire () const

unsigned int	NTdc () const

double	TRes () const

PixelMap	CreateMap (detinfo::DetectorPropertiesData const &detProp, const std::vector< art::Ptr< sim::SimChannel > > &slice)

PixelMap	CreateMap (detinfo::DetectorPropertiesData const &detProp, const std::vector< const sim::SimChannel * > &slice)

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

Private Member Functions
double	_getIntercept (geo::WireID wireid) const

void	_cacheIntercepts ()

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

double	fThreshold
	charge threshold for each time tick, below which isn't added to pixel map More...

unsigned short	fUnwrapped
	Use unwrapped pixel maps? More...

bool	fProtoDUNE
	Do we want to use this for particle extraction from protoDUNE? More...

unsigned int	fTotHits
	How many ROIs above threshold? More...

geo::GeometryCore const *	fGeometry

std::vector< double >	fVDPlane0

std::vector< double >	fVDPlane1

Detailed Description

Producer algorithm for PixelMap, input to CVN neural net.

Definition at line 29 of file PixelMapSimProducer.h.

Constructor & Destructor Documentation

cvn::PixelMapSimProducer::PixelMapSimProducer	(	unsigned int	nWire,
		unsigned int	nTdc,
		double	tRes,
		double	threshold = `0.`
	)

Definition at line 31 of file PixelMapSimProducer.cxx.

                                                                                                               :
     fNWire(nWire),
     fNTdc(nTdc),
     fTRes(tRes),
     fThreshold(threshold),
     fUnwrapped(2),
     fProtoDUNE(false),
     fTotHits(0)
   {
 
     fGeometry = &*(art::ServiceHandle<geo::Geometry>());  
     if (fGeometry->DetectorName().find("dunevd10kt_3view") != std::string::npos)
       _cacheIntercepts();
     
   }

cvn::PixelMapSimProducer::PixelMapSimProducer ( )

Definition at line 47 of file PixelMapSimProducer.cxx.

   {
     fGeometry = &*(art::ServiceHandle<geo::Geometry>());  
     if (fGeometry->DetectorName().find("dunevd10kt_3view") != std::string::npos)
       _cacheIntercepts();
   }

Member Function Documentation

void cvn::PixelMapSimProducer::_cacheIntercepts ( )

private

Definition at line 165 of file PixelMapSimProducer.cxx.

                                             {
    
     // double spacing = 0.847;
     for(int plane = 0; plane < 2; plane++){
       
       int nCRM_row = 6;
       for(int diag_tpc = 0; diag_tpc < nCRM_row; diag_tpc++){
         
         unsigned int nWiresTPC = fGeometry->Nwires(plane, 0, 0);
         int tpc_id = plane == 0 ? (nCRM_row+1)*diag_tpc : (nCRM_row-1)*(nCRM_row-diag_tpc);
         geo::WireID start = geo::WireID(0, tpc_id, plane, 0); 
         geo::WireID end = geo::WireID(0, tpc_id, plane, nWiresTPC-1);
 
         double start_intercept = _getIntercept(start);
         double end_intercept = _getIntercept(end);
         if(plane == 0){
           fVDPlane0.push_back(start_intercept);
           fVDPlane0.push_back(end_intercept);
         }
         else{
           fVDPlane1.push_back(_getIntercept(end));
           fVDPlane1.push_back(_getIntercept(start));
         }
       }
 
     }
   }

double cvn::PixelMapSimProducer::_getIntercept ( geo::WireID wireid ) const

private

Definition at line 152 of file PixelMapSimProducer.cxx.

   {
     const geo::WireGeo* pwire = fGeometry->WirePtr(wireid);
     geo::Point_t center = pwire->GetCenter<geo::Point_t>();
     double slope = 0.;
     if(!pwire->isVertical()) slope = pwire->TanThetaZ();
     
     double intercept = center.Y() - slope*center.Z();
     if(wireid.Plane == 2) intercept = 0.;
     
     return intercept;
   }

PixelMap cvn::PixelMapSimProducer::CreateMap	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< art::Ptr< sim::SimChannel > > &	slice
	)

Definition at line 54 of file PixelMapSimProducer.cxx.

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

PixelMap cvn::PixelMapSimProducer::CreateMap	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< const sim::SimChannel * > &	slice
	)

Definition at line 64 of file PixelMapSimProducer.cxx.

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

PixelMap cvn::PixelMapSimProducer::CreateMapGivenBoundary	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< const sim::SimChannel * > &	cluster,
		const Boundary &	bound
	)

Definition at line 71 of file PixelMapSimProducer.cxx.

   {
 
     PixelMap pm(fNWire, fNTdc, bound);
     
     for(size_t iHit = 0; iHit < cluster.size(); ++iHit)
     {
       
       const sim::SimChannel* reco_wire = cluster[iHit];
       auto& ROIs = reco_wire->TDCIDEMap();
       if(!(ROIs.size())) continue;
       // auto ROIs = reco_wire->SignalROI();
 
       std::vector<geo::WireID> wireids = fGeometry->ChannelToWire(reco_wire->Channel());
       if(!wireids.size()) continue;
       geo::WireID wireid = wireids[0];
       
       // if(wireids.size() > 1){
       //   for(auto iwire : wireids)
       //     if(iwire.Plane == reco_wire->View()) wireid = iwire;
       // }
       unsigned int tempWire  = wireid.Wire;
       unsigned int tempPlane = wireid.Plane;
 
       if(!fProtoDUNE){
         if(fUnwrapped == 1){
           if (fGeometry->DetectorName().find("dunevd10kt_3view") != std::string::npos){
             GetDUNEVertDrift3ViewGlobalWire(wireid.Wire, wireid.Plane,wireid.TPC,tempWire,tempPlane);
           }
         }
         else if(fUnwrapped == 2){
           // Old method that has problems with the APA crossers, kept for old times' sake
           GetDUNEGlobalWire(wireid.Wire,wireid.Plane,wireid.TPC,tempWire,tempPlane);
         }
       }
       else{
         GetProtoDUNEGlobalWire(wireid.Wire,wireid.Plane,wireid.TPC,tempWire,tempPlane);
       }
 
       for(auto iROI = ROIs.begin(); iROI != ROIs.end(); ++iROI){
         auto& ROI = *iROI;
         auto tick = ROI.first;
         double temptdc  = (double)tick;
         double charge =  0.005*reco_wire->Charge(tick); 
         if(!(charge > fThreshold)) continue;   
         // Leigh: Simple modification to unwrap the collection view wire plane
         if(!fProtoDUNE){
           if(fUnwrapped == 1){
             // Jeremy: Autodetect geometry for DUNE 10kt module. Is this a bad idea??
             if (fGeometry->DetectorName() == "dune10kt_v1") {
               if (wireid.TPC%6 == 0 or wireid.TPC%6 == 5) continue; // Skip dummy TPCs in 10kt module
               GetDUNE10ktGlobalWireTDC(detProp, wireid.Wire,(double)tick,
                 wireid.Plane,wireid.TPC,tempWire,tempPlane,temptdc);
             }
             else if (fGeometry->DetectorName().find("dunevd10kt_3view") == std::string::npos){
               GetDUNEGlobalWireTDC(detProp, wireid.Wire,(double)tick,
                 wireid.Plane,wireid.TPC,tempWire,tempPlane,temptdc);
             }
           }
         }
         const double pe  = charge;
         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;
   }

Boundary cvn::PixelMapSimProducer::DefineBoundary	(	detinfo::DetectorPropertiesData const &	detProp,
		const std::vector< const sim::SimChannel * > &	cluster
	)

Get boundaries for pixel map representation of cluster.

Definition at line 193 of file PixelMapSimProducer.cxx.

   {
 
     std::vector<int> wire_0;
     std::vector<int> wire_1;
     std::vector<int> wire_2;
     
     std::vector<double> twire_0;
     std::vector<double> twire_1;
     std::vector<double> twire_2;
    
     double tsum_0 = 0., tsum_1 = 0., tsum_2 = 0.;
     int total_t0 = 0, total_t1 = 0, total_t2 = 0;
 
     for(size_t iHit = 0; iHit < cluster.size(); ++iHit)
     {
       const sim::SimChannel* reco_wire = cluster[iHit];
       auto& ROIs = reco_wire->TDCIDEMap();
       if(!(ROIs.size())) continue;
 
       std::vector<geo::WireID> wireids = fGeometry->ChannelToWire(reco_wire->Channel());
       if(!wireids.size()) continue;
       geo::WireID wireid = wireids[0];
       
       // if(wireids.size() > 1){
       //   for(auto iwire : wireids)
       //     if(iwire.Plane == reco_wire->View()) wireid = iwire;
       // }
       unsigned int globalWire  = wireid.Wire;
       unsigned int globalPlane = wireid.Plane;
      
       if(!fProtoDUNE){
         if(fUnwrapped == 1){
           if (fGeometry->DetectorName().find("dunevd10kt_3view") != std::string::npos){
             GetDUNEVertDrift3ViewGlobalWire(wireid.Wire, wireid.Plane,wireid.TPC,globalWire,globalPlane);
           }
         }
         else if(fUnwrapped == 2){
           // Old method that has problems with the APA crossers, kept for old times' sake
           GetDUNEGlobalWire(wireid.Wire,wireid.Plane,wireid.TPC,globalWire,globalPlane);
         }
       }
       else{
         GetProtoDUNEGlobalWire(wireid.Wire,wireid.Plane,wireid.TPC,globalWire,globalPlane);
       }
 
       for(auto iROI = ROIs.begin(); iROI != ROIs.end(); ++iROI){
         auto& ROI = *iROI;
         auto tick = ROI.first;
         // for(int tick = ROI.begin_index(); tick < (int)ROI.end_index(); tick++){
           
         double globalTime  = (double)tick;
         double charge = 0.005*reco_wire->Charge(tick);
         if(!(charge > fThreshold)) continue;  
         // Leigh: Simple modification to unwrap the collection view wire plane
         if(!fProtoDUNE){
           if(fUnwrapped == 1){
             // Jeremy: Autodetect geometry for DUNE 10kt module. Is this a bad idea??
             if (fGeometry->DetectorName() == "dune10kt_v1") {
               if (wireid.TPC%6 == 0 or wireid.TPC%6 == 5) continue; // Skip dummy TPCs in 10kt module
               GetDUNE10ktGlobalWireTDC(detProp, wireid.Wire,(double)tick,
                 wireid.Plane,wireid.TPC,globalWire,globalPlane,globalTime);
             }
             else if (fGeometry->DetectorName().find("dunevd10kt_3view") == std::string::npos){
               GetDUNEGlobalWireTDC(detProp, wireid.Wire,(double)tick,
                 wireid.Plane,wireid.TPC,globalWire,globalPlane,globalTime);
             }
           }
         }
 
         if(globalPlane==0){
           tsum_0 += globalTime;
           total_t0++;
           wire_0.push_back(globalWire);
           twire_0.push_back((double)tick);
         }
         if(globalPlane==1){
           tsum_1 += globalTime;
           total_t1++;
           wire_1.push_back(globalWire);
           twire_1.push_back((double)tick);
         }
         if(globalPlane==2){
           tsum_2 += globalTime;
           total_t2++;
           wire_2.push_back(globalWire);
           twire_2.push_back((double)tick);
         }
         
         // }
       }
    
       // if(!none_threshold){
       //   if(globalPlane==0){
       //     wire_0.push_back(globalWire);
       //     twire_0.push_back((double)min_tick);
       //   }
       //   if(globalPlane==1){
       //     wire_1.push_back(globalWire);
       //     twire_1.push_back((double)min_tick);
       //   }
       //   if(globalPlane==2){
       //     wire_2.push_back(globalWire);
       //     twire_2.push_back((double)min_tick);
       //   }
       // }
 
     }
     double tmean_0 = tsum_0/total_t0; 
     double tmean_1 = tsum_1/total_t1; 
     double tmean_2 = tsum_2/total_t2; 
     
 
     //auto maxwireelement_0= std::max_element(wire_0.begin(), wire_0.end());
     //std::cout<<"maxwire 0: "<<*maxwireelement_0<<std::endl;
     //auto maxwireelement_1= std::max_element(wire_1.begin(), wire_1.end());
     //std::cout<<"maxwire 1: "<<*maxwireelement_1<<std::endl;
     //auto maxwireelement_2= std::max_element(wire_2.begin(), wire_2.end());
     //std::cout<<"maxwire 2: "<<*maxwireelement_2<<std::endl;
 
 
     std::vector<int> bwire_0;
     std::vector<int> bwire_1;
     std::vector<int> bwire_2;
     for(int i = 0; i < (int)wire_0.size(); i++){
       double t = twire_0[i];
       if(std::abs(t-tmean_0) < (double)fTRes)
         bwire_0.push_back(wire_0[i]);
     }
     for(int i = 0; i < (int)wire_1.size(); i++){
       double t = twire_1[i];
       if(std::abs(t-tmean_1) < (double)fTRes)
         bwire_1.push_back(wire_1[i]);
     }
     for(int i = 0; i < (int)wire_2.size(); i++){
       double t = twire_2[i];
       if(std::abs(t-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) <<std::endl;}
     if(bwire_1.size() > 0) { minwire_1 = *minwireelement_1-1; std::cout<<"minwire 1: "<<(*minwireelement_1) <<std::endl;}
     if(bwire_2.size() > 0) { minwire_2 = *minwireelement_2-1; std::cout<<"minwire 2: "<<(*minwireelement_2) <<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;
   }

void cvn::PixelMapSimProducer::GetDUNE10ktGlobalWireTDC	(	detinfo::DetectorPropertiesData const &	detProp,
		unsigned int	localWire,
		double	localTDC,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int &	globalWire,
		unsigned int &	globalPlane,
		double &	globalTDC
	)		const

Definition at line 456 of file PixelMapSimProducer.cxx.

   {
     unsigned int nWiresTPC = 400;
     unsigned int wireGap = 4;
     double driftLen = fGeometry->TPC(tpc).DriftDistance();
     double apaLen = fGeometry->TPC(tpc).Width() - fGeometry->TPC(tpc).ActiveWidth();
     double driftVel = detProp.DriftVelocity();
     unsigned int drift_size = (driftLen / driftVel) * 2; // Time in ticks to cross a TPC 
     unsigned int apa_size   = 4*(apaLen / driftVel) * 2; // Width of the whole APA in TDC
 
 
     globalWire = 0;
     globalPlane = 0;
 
     // Collection plane has more wires
     if(plane == 2){
       nWiresTPC = 480;
       wireGap = 5;
       globalPlane = 2;
     }
 
     bool includeZGap = true;
     if(includeZGap) nWiresTPC += wireGap;
 
     // 10kt has four real TPCs and two dummies in each slice
     //
     //                 |--|-----|-----|-----|-----|--| /  /
     //      y ^        |11| 10  |  9  |  8  |  7  | 6|/  /
     //        | -| z   |--|-----|-----|-----|-----|--|  /
     //        | /      | 5|  4  |  3  |  2  |  1  | 0| /
     //  x <---|/       |--|-----|-----|-----|-----|--|/
     //                     ^  wires  ^ ^  wires  ^
     //
     // We already filtered out the dummies, so we can assume 0->3 as follows:
     //
     //                 |-----|-----|-----|-----| /  /
     //      y ^        |  7  |  6  |  5  |  4  |/  /
     //        | -| z   |-----|-----|-----|-----|  /
     //        | /      |  3  |  2  |  1  |  0  | /
     //  x <---|/       |-----|-----|-----|-----|/
     //                  ^  wires  ^ ^  wires  ^
     //
 
     size_t tpc_x = (tpc%6) - 1;   // x coordinate in 0->4 range
     size_t tpc_xy = (tpc%12) - 1; // xy coordinate as 0->3 & 6->9 (converted from 1->4, 7->10)
     if (tpc_xy > 3) tpc_xy -= 2;  // now subtract 2 so it's in the 0->7 range
 
     // Induction views depend on the drift direction
     if (plane < 2 and tpc%2 == 1) globalPlane = !plane;
     else globalPlane = plane;
 
     int offset = 752; // Offset between upper and lower modules in induction views, from Robert & Dorota's code
     // Second induction plane gets offset from the back of the TPC
     if (globalPlane != 1) globalWire += (tpc/12)*nWiresTPC;
     else globalWire += ((300-tpc)/12)*nWiresTPC;
     // Reverse wires and add offset for upper modules in induction views
     if (tpc_xy > 3 and globalPlane < 2) globalWire += fGeometry->Nwires(globalPlane, tpc, 0) + offset - localWire;
     else globalWire += localWire;
 
     if (tpc_x % 2 == 0) globalTDC = localTDC;
     else globalTDC = (2*drift_size) - localTDC;
     if (tpc_x > 1) globalTDC += 2 * (drift_size + apa_size);
 
   } // function PixelMapSimProducer::GetDUNE10ktGlobalWireTDC

void cvn::PixelMapSimProducer::GetDUNEGlobalWire	(	unsigned int	localWire,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int &	globalWire,
		unsigned int &	globalPlane
	)		const

Function to convert to a global unwrapped wire number.

Definition at line 354 of file PixelMapSimProducer.cxx.

   {
     unsigned int nWiresTPC = 400;
 
     globalWire = localWire;
     globalPlane = 0;
 
     // Collection plane has more wires
     if(plane == 2){
       nWiresTPC=480;
       globalPlane = 2;
     }
 
     // Workspace geometry has two drift regions
     //                  |-----|-----| /  /
     //      y ^         |  3  |  2  |/  /
     //        | -| z    |-----|-----|  /
     //        | /       |  1  |  0  | /
     //  x <---|/        |-----|-----|/
     //
 
     int tpcMod4 = tpc%4;
 
     // Induction views depend on the drift direction
     if(plane < 2){
       // For drift in negative x direction keep U and V as defined.
       if(tpcMod4 == 0 || tpcMod4 == 3){
         globalPlane = plane;
       }
       // For drift in positive x direction, swap U and V.
       else{
         if(plane == 0) globalPlane = 1;
         else globalPlane = 0;
       }
     }
 
     if(globalPlane != 1){
       globalWire += (tpc/4)*nWiresTPC;
     }
     else{
       globalWire += ((23-tpc)/4)*nWiresTPC;
     }
 
   }

void cvn::PixelMapSimProducer::GetDUNEGlobalWireTDC	(	detinfo::DetectorPropertiesData const &	detProp,
		unsigned int	localWire,
		double	localTDC,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int &	globalWire,
		unsigned int &	globalPlane,
		double &	globalTDC
	)		const

Definition at line 400 of file PixelMapSimProducer.cxx.

   {
 
     unsigned int nWiresTPC = 400;
     unsigned int wireGap = 4;
     double driftLen = fGeometry->TPC(tpc,0).DriftDistance();
     double apaLen = fGeometry->TPC(tpc,0).Width() - fGeometry->TPC(tpc,0).ActiveWidth();
     double driftVel = detProp.DriftVelocity();
     unsigned int drift_size = (driftLen / driftVel) * 2; // Time in ticks to cross a TPC 
     unsigned int apa_size   = 4*(apaLen / driftVel) * 2; // Width of the whole APA in TDC
 
     globalWire = 0;
     globalPlane = 0;
 //    int dir = fGeometry->TPC(tpc,0).DetectDriftDirection();
 
     // Collection plane has more wires
     if(plane == 2){
       nWiresTPC = 480;
       wireGap = 5;
       globalPlane = 2;
     }
 
     bool includeZGap = true;
     if(includeZGap) nWiresTPC += wireGap;
 
     // Workspace geometry has two drift regions
     //                  |-----|-----| /  /
     //      y ^         |  3  |  2  |/  /
     //        | -| z    |-----|-----|  /
     //        | /       |  1  |  0  | /
     //  x <---|/        |-----|-----|/
     //
     int tpcMod4 = tpc%4;
     // Induction views depend on the drift direction
     if (plane < 2 and tpc%2 == 1) globalPlane = !plane;
     else globalPlane = plane;
 
     int offset = 752; // Offset between upper and lower modules in induction views, from Robert & Dorota's code
     // Second induction plane gets offset from the back of the TPC
     if (globalPlane != 1) globalWire += (tpc/4)*nWiresTPC;
     else globalWire += ((23-tpc)/4)*nWiresTPC;
     // Reverse wires and add offset for upper modules in induction views
     // Nitish : what's the difference between Nwires here and nWiresTPC?
     if (tpcMod4 > 1 and globalPlane < 2) globalWire += fGeometry->Nwires(globalPlane, tpc, 0) + offset - localWire;
     else globalWire += localWire;
 
     if(tpcMod4 == 0 || tpcMod4 == 2){
       globalTDC = drift_size - localTDC;
     }
     else{
       globalTDC = localTDC + drift_size + apa_size;
     }
   }

void cvn::PixelMapSimProducer::GetDUNEVertDrift3ViewGlobalWire	(	unsigned int	localWire,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int &	globalWire,
		unsigned int &	globalPlane
	)		const

Definition at line 583 of file PixelMapSimProducer.cxx.

   {
     // Preliminary function for VD Geometries (3 View for now)
     // 1x6x6 -- single drift volume should make things significantly simpler
     
     int nCRM_row = 6;
     // spacing between y-intercepts of parallel wires in a given plane. 
     double spacing = 0.847; 
     
     globalPlane = plane;
     unsigned int nWiresTPC = fGeometry->Nwires(globalPlane, tpc, 0);
     
     if(globalPlane < 2){
       
       geo::WireID wire_id = geo::WireID(0, tpc, globalPlane, localWire);
       double wire_intercept = _getIntercept(wire_id);
       double low_bound, upper_bound; 
       int start, end, diag_tpc;
       // get wires on diagonal CRMs and their intercepts which bound the current wire's intercept 
       if(globalPlane == 0){
         start = std::lower_bound(fVDPlane0.begin(), fVDPlane0.end(), wire_intercept) - fVDPlane0.begin() - 1;
         end = std::upper_bound(fVDPlane0.begin(), fVDPlane0.end(), wire_intercept) - fVDPlane0.begin();
         low_bound = fVDPlane0[start];
         upper_bound = fVDPlane0[end];
         diag_tpc = (start/2);
       }
       else{
         end = std::lower_bound(fVDPlane1.begin(), fVDPlane1.end(), wire_intercept) - fVDPlane1.begin() - 1;
         start = std::upper_bound(fVDPlane1.begin(), fVDPlane1.end(), wire_intercept) - fVDPlane1.begin();
         low_bound = fVDPlane1[end];
         upper_bound = fVDPlane1[start];
         diag_tpc = (nCRM_row-(end/2) - 1);
       }
       // if the intercept of the wire is in between two diagonal CRMs, assign it to the diagonal CRM its closest to 
       if((start % 2)^globalPlane){
         
         int diag_idx = diag_tpc + !globalPlane;
         globalWire = (wire_intercept > (low_bound+upper_bound)*0.5) ? (nWiresTPC-1)*diag_idx + !globalPlane : (nWiresTPC-1)*diag_idx + globalPlane;
       }
       // otherwise assign it to the closest wire within the same CRM
       else{
         int diag_idx = diag_tpc;
         int offset = globalPlane ? std::round((upper_bound - wire_intercept)/spacing) : std::round((wire_intercept-low_bound)/spacing);
         globalWire = (nWiresTPC-1)*diag_idx + offset + 1;
         
       }
     }
     else{
       int tpc_z = tpc/6;
       globalWire = localWire + tpc_z*nWiresTPC;
     }
  
   }

void cvn::PixelMapSimProducer::GetProtoDUNEGlobalWire	(	unsigned int	localWire,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int &	globalWire,
		unsigned int &	globalPlane
	)		const

Definition at line 525 of file PixelMapSimProducer.cxx.

   { 
     unsigned int nWiresTPC = 400;
     
     globalWire = localWire;
     globalPlane = 0;
     
     // Collection plane has more wires
     if(plane == 2){
       nWiresTPC=480;
       globalPlane = 2;
     }
     
     // ProtoDUNE has a central CPA so times are fine
     // It (annoyingly) has two dummy TPCs on the sides
     //                  
     //      y ^       |-|-----|-----|-|   /
     //        | -| z  | |     |     | |  /
     //        | /     |3|  2  |  1  |0| /
     //  x <---|/      |-|-----|-----|-|/
     //
     
     int tpcMod4 = tpc%4;
     // tpcMod4: 1 for -ve drift, 2 for +ve drift
     // Induction views depend on the drift direction
     if(plane < 2){
       // For drift in negative x direction keep U and V as defined.
       if(tpcMod4 == 1){
         globalPlane = plane;
       }
       // For drift in positive x direction, swap U and V.
       else{
         if(plane == 0) globalPlane = 1;
         else globalPlane = 0;
       }
     }
     
     if(globalPlane != 1){
       globalWire += (tpc/4)*nWiresTPC;
     }
     else{
       globalWire += ((12-tpc)/4)*nWiresTPC;
     }
   
   } // function PixelMapSimProducer::GetProtoDUNEGlobalWire

void cvn::PixelMapSimProducer::GetProtoDUNEGlobalWireTDC	(	unsigned int	localWire,
		double	localTDC,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int &	globalWire,
		double &	globalTDC,
		unsigned int &	globalPlane
	)		const

Definition at line 573 of file PixelMapSimProducer.cxx.

   {
     // We can just use the existing function to get the global wire & plane
  //   GetProtoDUNEGlobalWire(localWire, plane, tpc, globalWire, globalPlane);
     GetDUNEGlobalWire(localWire, plane, tpc, globalWire, globalPlane);
     // Implement additional mirroring here?
 
   } // function GetProtoDUNEGlobalWireTDC

unsigned int cvn::PixelMapSimProducer::NROI ( )

inline

Definition at line 42 of file PixelMapSimProducer.h.

42 {return fTotHits;};

cvn::PixelMapSimProducer::fTotHits

unsigned int fTotHits

How many ROIs above threshold?

Definition: PixelMapSimProducer.h:82

unsigned int cvn::PixelMapSimProducer::NTdc ( ) const

inline

Definition at line 62 of file PixelMapSimProducer.h.

62 {return fNTdc;};

cvn::PixelMapSimProducer::fNTdc

unsigned int fNTdc

Number of tdcs, width of pixel map.

Definition: PixelMapSimProducer.h:76

unsigned int cvn::PixelMapSimProducer::NWire ( ) const

inline

Definition at line 61 of file PixelMapSimProducer.h.

61 {return fNWire;};

cvn::PixelMapSimProducer::fNWire

unsigned int fNWire

Number of wires, length for pixel maps.

Definition: PixelMapSimProducer.h:75

void cvn::PixelMapSimProducer::SetProtoDUNE ( )

inline

Definition at line 36 of file PixelMapSimProducer.h.

36 {fProtoDUNE = true;};

cvn::PixelMapSimProducer::fProtoDUNE

bool fProtoDUNE

Do we want to use this for particle extraction from protoDUNE?

Definition: PixelMapSimProducer.h:80

void cvn::PixelMapSimProducer::SetUnwrapped ( unsigned short unwrap )

inline

Definition at line 35 of file PixelMapSimProducer.h.

35 {fUnwrapped = unwrap;};

cvn::PixelMapSimProducer::fUnwrapped

unsigned short fUnwrapped

Use unwrapped pixel maps?

Definition: PixelMapSimProducer.h:79

double cvn::PixelMapSimProducer::TRes ( ) const

inline

Definition at line 63 of file PixelMapSimProducer.h.

63 {return fTRes;};

cvn::PixelMapSimProducer::fTRes

double fTRes

Timing resolution for pixel map.

Definition: PixelMapSimProducer.h:77

Member Data Documentation

geo::GeometryCore const* cvn::PixelMapSimProducer::fGeometry

private

Definition at line 84 of file PixelMapSimProducer.h.

unsigned int cvn::PixelMapSimProducer::fNTdc

private

Number of tdcs, width of pixel map.

Definition at line 76 of file PixelMapSimProducer.h.

unsigned int cvn::PixelMapSimProducer::fNWire

private

Number of wires, length for pixel maps.

Definition at line 75 of file PixelMapSimProducer.h.

bool cvn::PixelMapSimProducer::fProtoDUNE

private

Do we want to use this for particle extraction from protoDUNE?

Definition at line 80 of file PixelMapSimProducer.h.

double cvn::PixelMapSimProducer::fThreshold

private

charge threshold for each time tick, below which isn't added to pixel map

Definition at line 78 of file PixelMapSimProducer.h.

unsigned int cvn::PixelMapSimProducer::fTotHits

private

How many ROIs above threshold?

Definition at line 82 of file PixelMapSimProducer.h.

double cvn::PixelMapSimProducer::fTRes

private

Timing resolution for pixel map.

Definition at line 77 of file PixelMapSimProducer.h.

unsigned short cvn::PixelMapSimProducer::fUnwrapped

private

Use unwrapped pixel maps?

Definition at line 79 of file PixelMapSimProducer.h.

std::vector<double> cvn::PixelMapSimProducer::fVDPlane0

private

Definition at line 85 of file PixelMapSimProducer.h.

std::vector<double> cvn::PixelMapSimProducer::fVDPlane1

private

Definition at line 86 of file PixelMapSimProducer.h.

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

dunereco/dunereco/CVN/art/PixelMapSimProducer.h
dunereco/dunereco/CVN/art/PixelMapSimProducer.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation