WireCell::SigProc::PMTNoiseROI Class Reference

#include <PMTNoiseROI.h>

Public Member Functions
	PMTNoiseROI (int start_bin, int end_bin, int peak, int channel, float peak_height)
	~PMTNoiseROI ()
int	get_start_bin ()
int	get_end_bin ()
void	insert_peak (int peak)
void	insert_uwires (int wire_no, float peak_height)
void	insert_vwires (int wire_no, float peak_height)
std::vector< int > &	get_peaks ()
std::vector< int > &	get_uwires ()
std::vector< int > &	get_vwires ()
std::vector< int > &	get_wwires ()
std::map< int, float > &	get_uwires_peak_heights ()
std::map< int, float > &	get_vwires_peak_heights ()
std::map< int, float > &	get_wwires_peak_heights ()
float	get_average_uwires_peak_height (int group)
float	get_average_vwires_peak_height (int group)
float	get_max_uwires_peak_height (int group)
float	get_max_vwires_peak_height (int group)
float	get_average_wwires_peak_height ()
float	get_max_wwires_peak_height ()
std::vector< std::vector< int > > &	get_sorted_uwires ()
std::vector< std::vector< int > > &	get_sorted_vwires ()
std::vector< int > &	get_sorted_wwires ()
void	sort_wires (int nwire=1)
bool	merge_ROI (PMTNoiseROI &ROI)

Private Attributes
std::vector< int >	peaks
int	start_bin
int	end_bin
std::map< int, float >	uwires_peak_heights
std::map< int, float >	vwires_peak_heights
std::map< int, float >	wwires_peak_heights
std::vector< int >	induction_uwires
std::vector< int >	induction_vwires
std::vector< int >	collection_wwires
std::vector< std::vector< int > >	sorted_ind_uwires
std::vector< std::vector< int > >	sorted_ind_vwires
std::vector< int >	sorted_col_wwires

Detailed Description

Definition at line 9 of file PMTNoiseROI.h.

Constructor & Destructor Documentation

PMTNoiseROI::PMTNoiseROI	(	int	start_bin,
		int	end_bin,
		int	peak,
		int	channel,
		float	peak_height
	)

Definition at line 10 of file PMTNoiseROI.cxx.

  : start_bin(start_bin)
  , end_bin(end_bin)
{
  peaks.push_back(peak);
  collection_wwires.push_back(channel);
  wwires_peak_heights[channel] = fabs(peak_height);
}

PMTNoiseROI::~PMTNoiseROI

(

)

Definition at line 20 of file PMTNoiseROI.cxx.

                         {
  
}

Member Function Documentation

float PMTNoiseROI::get_average_uwires_peak_height

(

int

group

)

Definition at line 90 of file PMTNoiseROI.cxx.

                                                          {
  float ave = 0 ;
  float ave1 = 0 ;
  if (group < int(sorted_ind_uwires.size())){
    for (int i=0;i!=int(sorted_ind_uwires.at(group).size());i++){
      ave += uwires_peak_heights[sorted_ind_uwires.at(group).at(i)];
      ave1 ++;
    }
  }
      
  if (ave1 >0){
    return ave/ave1;
  }else{
    return 0;
  }
}

float PMTNoiseROI::get_average_vwires_peak_height

(

int

group

)

Definition at line 107 of file PMTNoiseROI.cxx.

                                                          {
  float ave = 0 ;
  float ave1 = 0 ;
  if (group < int(sorted_ind_vwires.size())){
    for (int i=0;i!=int(sorted_ind_vwires.at(group).size());i++){
      ave += vwires_peak_heights[sorted_ind_vwires.at(group).at(i)];
      ave1 ++;
    }
  }
  if (ave1 >0){
    return ave/ave1;
  }else{
    return 0;
  }
}

float PMTNoiseROI::get_average_wwires_peak_height

(

)

Definition at line 53 of file PMTNoiseROI.cxx.

                                                 {
  float ave = 0 ;
  float ave1 = 0 ;
  for (int i=0;i!=int(sorted_col_wwires.size());i++){
    ave += wwires_peak_heights[sorted_col_wwires.at(i)];
    ave1 ++;
  }
  if (ave1 >0){
    return ave/ave1;
  }else{
    return 0;
  }
}

int WireCell::SigProc::PMTNoiseROI::get_end_bin ( )

inline

Definition at line 15 of file PMTNoiseROI.h.

{return end_bin;};

float PMTNoiseROI::get_max_uwires_peak_height

(

int

group

)

Definition at line 68 of file PMTNoiseROI.cxx.

                                                      {
  float max = 0;
  if (group < int(sorted_ind_uwires.size())){
    for (int i=0;i!=int(sorted_ind_uwires.at(group).size());i++){
      if (uwires_peak_heights[sorted_ind_uwires.at(group).at(i)] > max)
        max= uwires_peak_heights[sorted_ind_uwires.at(group).at(i)];
    }
  }
  return max;
}

float PMTNoiseROI::get_max_vwires_peak_height

(

int

group

)

Definition at line 79 of file PMTNoiseROI.cxx.

                                                      {
  float max = 0;
  if (group < int(sorted_ind_vwires.size())){
    for (int i=0;i!=int(sorted_ind_vwires.at(group).size());i++){
      if (vwires_peak_heights[sorted_ind_vwires.at(group).at(i)] > max)
        max= vwires_peak_heights[sorted_ind_vwires.at(group).at(i)];
    }
  }
  return max;
}

float PMTNoiseROI::get_max_wwires_peak_height

(

)

Definition at line 44 of file PMTNoiseROI.cxx.

                                             {
  float max = 0;
  for (int i=0;i!=int(sorted_col_wwires.size());i++){
    if (wwires_peak_heights[sorted_col_wwires.at(i)] > max)
      max= wwires_peak_heights[sorted_col_wwires.at(i)];
  }
  return max;
}

std::vector<int>& WireCell::SigProc::PMTNoiseROI::get_peaks ( )

inline

Definition at line 23 of file PMTNoiseROI.h.

{return peaks;}

std::vector<std::vector<int> >& WireCell::SigProc::PMTNoiseROI::get_sorted_uwires ( )

inline

Definition at line 44 of file PMTNoiseROI.h.

{return sorted_ind_uwires;}     

std::vector<std::vector<int> >& WireCell::SigProc::PMTNoiseROI::get_sorted_vwires ( )

inline

Definition at line 45 of file PMTNoiseROI.h.

{return sorted_ind_vwires;}     

std::vector<int>& WireCell::SigProc::PMTNoiseROI::get_sorted_wwires ( )

inline

Definition at line 46 of file PMTNoiseROI.h.

{return sorted_col_wwires;}    

int WireCell::SigProc::PMTNoiseROI::get_start_bin ( )

inline

Definition at line 14 of file PMTNoiseROI.h.

{return start_bin;};

std::vector<int>& WireCell::SigProc::PMTNoiseROI::get_uwires ( )

inline

Definition at line 25 of file PMTNoiseROI.h.

{return induction_uwires;}

std::map<int,float>& WireCell::SigProc::PMTNoiseROI::get_uwires_peak_heights ( )

inline

Definition at line 29 of file PMTNoiseROI.h.

{return uwires_peak_heights;} 

std::vector<int>& WireCell::SigProc::PMTNoiseROI::get_vwires ( )

inline

Definition at line 26 of file PMTNoiseROI.h.

{return induction_vwires;}

std::map<int,float>& WireCell::SigProc::PMTNoiseROI::get_vwires_peak_heights ( )

inline

Definition at line 30 of file PMTNoiseROI.h.

{return vwires_peak_heights;} 

std::vector<int>& WireCell::SigProc::PMTNoiseROI::get_wwires ( )

inline

Definition at line 27 of file PMTNoiseROI.h.

{return collection_wwires;}

std::map<int,float>& WireCell::SigProc::PMTNoiseROI::get_wwires_peak_heights ( )

inline

Definition at line 31 of file PMTNoiseROI.h.

{return wwires_peak_heights;} 

void PMTNoiseROI::insert_peak

(

int

peak

)

Definition at line 24 of file PMTNoiseROI.cxx.

                                     {
  if (find(peaks.begin(),peaks.end(),peak)==peaks.end())
    peaks.push_back(peak);
}

void PMTNoiseROI::insert_uwires	(	int	wire_no,
		float	peak_height
	)

Definition at line 30 of file PMTNoiseROI.cxx.

                                                             {
  if (find(induction_uwires.begin(),induction_uwires.end(),wire_no)==induction_uwires.end()){
    induction_uwires.push_back(wire_no);
    uwires_peak_heights[wire_no] = peak_height;
  }
}

void PMTNoiseROI::insert_vwires	(	int	wire_no,
		float	peak_height
	)

Definition at line 37 of file PMTNoiseROI.cxx.

                                                             {
  if (find(induction_vwires.begin(),induction_vwires.end(),wire_no)==induction_vwires.end()){
    induction_vwires.push_back(wire_no);
    vwires_peak_heights[wire_no] = peak_height;
  }
}

bool PMTNoiseROI::merge_ROI

(

PMTNoiseROI &

ROI

)

Definition at line 124 of file PMTNoiseROI.cxx.

                                           {
  // decide how to merge two ROIs? ... 

  // one peak is contained in the other's range then merge
  if ( (peaks.at(0)>=ROI.get_start_bin() && peaks.at(0)<=ROI.get_end_bin())
       || (start_bin <= ROI.get_peaks().at(0) && end_bin >= ROI.get_peaks().at(0))){
    if (ROI.get_start_bin() < start_bin)
      start_bin = ROI.get_start_bin();
    if (ROI.get_end_bin() > end_bin)
      end_bin = ROI.get_end_bin();
    
    for (int i=0;i!=int(ROI.get_peaks().size());i++){
      if (find(peaks.begin(),peaks.end(),ROI.get_peaks().at(i)) == peaks.end())
        peaks.push_back(ROI.get_peaks().at(i));
    }

    for (int i=0;i!=int(ROI.get_wwires().size());i++){
      if (find(collection_wwires.begin(),collection_wwires.end(),ROI.get_wwires().at(i)) == collection_wwires.end()){
        collection_wwires.push_back(ROI.get_wwires().at(i));
      }
      
      // calculate wire height 
      if (wwires_peak_heights.find(ROI.get_wwires().at(i))==wwires_peak_heights.end()){
        wwires_peak_heights[ROI.get_wwires().at(i)] = ROI.get_wwires_peak_heights()[ROI.get_wwires().at(i)];
      }else{
        if (ROI.get_wwires_peak_heights()[ROI.get_wwires().at(i)] > wwires_peak_heights[ROI.get_wwires().at(i)])
          wwires_peak_heights[ROI.get_wwires().at(i)] = ROI.get_wwires_peak_heights()[ROI.get_wwires().at(i)];
      }
      
    }
    
    //    for (int i=0;i!=ROI.get_wwires_peak_heights().size();i++){
    // wwires_peak_heights.push_back(ROI.get_wwires_peak_heights().at(i));
    //}
    

    
    return true;
  }
 
  return false;
}

void PMTNoiseROI::sort_wires

(

int

nwire = 1

)

Definition at line 168 of file PMTNoiseROI.cxx.

                                     {

  std::vector<int> temp_wires;
  if (int(induction_uwires.size()) >=nwire){
    std::sort(induction_uwires.begin(),induction_uwires.end());
    temp_wires.push_back(induction_uwires.at(0));
    // do u wires;
    for (int i=1;i<int(induction_uwires.size());i++){
      if (induction_uwires.at(i) -  temp_wires.back() == 1){
        temp_wires.push_back(induction_uwires.at(i));
      }else{
        if (int(temp_wires.size()) >= nwire){
          sorted_ind_uwires.push_back(temp_wires);
        }
        temp_wires.clear();
        temp_wires.push_back(induction_uwires.at(i));
      }
    }
    if (int(temp_wires.size()) >= nwire){
      sorted_ind_uwires.push_back(temp_wires);
    }
    temp_wires.clear();
  }

  if (int(induction_vwires.size()) >=nwire){
    std::sort(induction_vwires.begin(),induction_vwires.end());
    temp_wires.push_back(induction_vwires.at(0));
    // do u wires;
    for (int i=1;i<int(induction_vwires.size());i++){
      if (induction_vwires.at(i) -  temp_wires.back() == 1){
        temp_wires.push_back(induction_vwires.at(i));
      }else{
        if (int(temp_wires.size()) >= nwire){
          sorted_ind_vwires.push_back(temp_wires);
        }
        temp_wires.clear();
        temp_wires.push_back(induction_vwires.at(i));
      }
    }
    if (int(temp_wires.size()) >= nwire){
      sorted_ind_vwires.push_back(temp_wires);
    }
    temp_wires.clear();
  }


  if (int(collection_wwires.size()) >=nwire){
    std::sort(collection_wwires.begin(),collection_wwires.end());
    temp_wires.push_back(collection_wwires.at(0));
    // do u wires;
    for (int i=1;i<int(collection_wwires.size());i++){
      if (collection_wwires.at(i) -  temp_wires.back() == 1){
        temp_wires.push_back(collection_wwires.at(i));
      }else{
        if (int(temp_wires.size()) >= nwire){
          for (int j=0;j!=int(temp_wires.size());j++){
            sorted_col_wwires.push_back(temp_wires.at(j));
          }
        }
        temp_wires.clear();
        temp_wires.push_back(collection_wwires.at(i));
      }
    }
    if (int(temp_wires.size()) >= nwire){
      for (int j=0;j!=int(temp_wires.size());j++){
        sorted_col_wwires.push_back(temp_wires.at(j));
      }
    }
    temp_wires.clear();
  }


}

Member Data Documentation

std::vector<int> WireCell::SigProc::PMTNoiseROI::collection_wwires

private

Definition at line 63 of file PMTNoiseROI.h.

int WireCell::SigProc::PMTNoiseROI::end_bin

private

Definition at line 54 of file PMTNoiseROI.h.

std::vector<int> WireCell::SigProc::PMTNoiseROI::induction_uwires

private

Definition at line 61 of file PMTNoiseROI.h.

std::vector<int> WireCell::SigProc::PMTNoiseROI::induction_vwires

private

Definition at line 62 of file PMTNoiseROI.h.

std::vector<int> WireCell::SigProc::PMTNoiseROI::peaks

private

Definition at line 52 of file PMTNoiseROI.h.

std::vector<int> WireCell::SigProc::PMTNoiseROI::sorted_col_wwires

private

Definition at line 67 of file PMTNoiseROI.h.

std::vector<std::vector<int> > WireCell::SigProc::PMTNoiseROI::sorted_ind_uwires

private

Definition at line 65 of file PMTNoiseROI.h.

std::vector<std::vector<int> > WireCell::SigProc::PMTNoiseROI::sorted_ind_vwires

private

Definition at line 66 of file PMTNoiseROI.h.

int WireCell::SigProc::PMTNoiseROI::start_bin

private

Definition at line 53 of file PMTNoiseROI.h.

std::map<int,float> WireCell::SigProc::PMTNoiseROI::uwires_peak_heights

private

Definition at line 56 of file PMTNoiseROI.h.

std::map<int,float> WireCell::SigProc::PMTNoiseROI::vwires_peak_heights

private

Definition at line 57 of file PMTNoiseROI.h.

std::map<int,float> WireCell::SigProc::PMTNoiseROI::wwires_peak_heights

private

Definition at line 58 of file PMTNoiseROI.h.

---

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

• wire-cell-build/sigproc/src/PMTNoiseROI.h
• wire-cell-build/sigproc/src/PMTNoiseROI.cxx