RegCNNBoundary.cxx

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \file    RegCNNBoundary.cxx
/// \brief   RegCNNBoundary for RegCNN PixelMap modified from CVNBoundary.cxx
/// \author  Ilsoo Seong - iseong@uci.edu
//           Wenjie Wu   - wenjieww@uci.edu
////////////////////////////////////////////////////////////////////////

#include  <iostream>
#include  <ostream>
#include  <utility>

#include "messagefacility/MessageLogger/MessageLogger.h"

#include "dunereco/RegCNN/func/RegCNNBoundary.h"

namespace cnn
{

  RegCNNBoundary::RegCNNBoundary(const int& nWire, // # of wires of pixel map
                     const int& nTDC,              // # of TDC of pixel map
                     const int& wRes,              // # of Wire merging
                             const int& tRes,              // # of TDC merging
                     const int& WireMeanX,
                     const int& WireMeanY,
                     const int& WireMeanZ,
                     const int& TDCMeanX,
                     const int& TDCMeanY,
                     const int& TDCMeanZ):
    fFirstWire{WireMeanX-(nWire*wRes)/2,
      WireMeanY-(nWire*wRes)/2,
      WireMeanZ-(nWire*wRes)/2},
    fLastWire{WireMeanX+(nWire*wRes)/2-wRes/2,
      WireMeanY+(nWire*wRes)/2-wRes/2,
      WireMeanZ+(nWire*wRes)/2-wRes/2},
    fFirstTDC{TDCMeanX - (nTDC*tRes)/2,
        TDCMeanY  - (nTDC*tRes)/2,
        TDCMeanZ  - (nTDC*tRes)/2},  
    fLastTDC{TDCMeanX + (nTDC*tRes)/2-tRes/2,
        TDCMeanY + (nTDC*tRes)/2-tRes/2,
        TDCMeanZ + (nTDC*tRes)/2-tRes/2}
  
  {
      if (fLastWire[0] - fFirstWire[0] != nWire*wRes - wRes/2)
          mf::LogError("RegCNNBoundary::RegCNNBoundary") << "RegCNN 1st Map Boundary Wrong"<<std::endl;
      if (fLastWire[1] - fFirstWire[1] != nWire*wRes - wRes/2)
          mf::LogError("RegCNNBoundary::RegCNNBoundary") << "RegCNN 2nd Map Boundary Wrong"<<std::endl;
      if (fLastWire[2] - fFirstWire[2] != nWire*wRes - wRes/2)
          mf::LogError("RegCNNBoundary::RegCNNBoundary") << "RegCNN 3rd Map Boundary Wrong"<<std::endl;
  }

  bool RegCNNBoundary::IsWithin(const int& wire, const int& tdc, const unsigned int& view)
  {
    bool inWireRcnne = (int) wire >= fFirstWire[view] && (int) wire < fLastWire[view];
    bool inTDCRcnne = (int) tdc >= fFirstTDC[view] &&
                      (int) tdc < fLastTDC[view];
    return inWireRcnne && inTDCRcnne;
  }


  std::ostream& operator<<(std::ostream& os, const RegCNNBoundary& b)
  {
    os<<"RegCNNBoundary with "
      <<"(first,last) wire X: ("    <<b.FirstWire(0)<<", "<<b.LastWire(0)
      <<"), (first,last) wire Y: (" <<b.FirstWire(1)<<", "<<b.LastWire(1)
      <<"), (first,last) wire Z: (" <<b.FirstWire(2)<<", "<<b.LastWire(2)
      <<"), (first,last) tdc X: ("  <<b.FirstTDC(0) <<", "<<b.LastTDC(0)
      <<"), (first,last) tdc Y: ("  <<b.FirstTDC(1) <<", "<<b.LastTDC(1)
      <<"), (first,last) tdc Z: ("  <<b.FirstTDC(2) <<", "<<b.LastTDC(2)
      <<")";

    return os;
  }
}