hash_shelf.h

Go to the documentation of this file.

//Copyright 2008,2009 Thomas A Caswell
//tcaswell@uchicago.edu
//http://jfi.uchicago.edu/~tcaswell
//
//This program is free software; you can redistribute it and/or modify
//it under the terms of the GNU General Public License as published by
//the Free Software Foundation; either version 3 of the License, or (at
//your option) any later version.
//
//This program is distributed in the hope that it will be useful, but
//WITHOUT ANY WARRANTY; without even the implied warranty of
//MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
//General Public License for more details.
//
//You should have received a copy of the GNU General Public License
//along with this program; if not, see <http://www.gnu.org/licenses>.
//
//Additional permission under GNU GPL version 3 section 7
//
//If you modify this Program, or any covered work, by linking or
//combining it with MATLAB (or a modified version of that library),
//containing parts covered by the terms of MATLAB User License, the
//licensors of this Program grant you additional permission to convey
//the resulting work.

#ifndef HASH_SHELF
#define HASH_SHELF


#include <vector>
#include <list>
#include <vector>




#include "part_def.h"

namespace utilities{

class Counted_vector;
class Array;
class Histogram2D;
class Wrapper_out;
class Md_store;



}

namespace tracking{
class hash_box;
class Master_box;

class Corr;
class Accumulator;

class Hash_shelf{
public:

  utilities::Tuplef get_img_dims()const
  {
    return utilities::Tuplef(img_dims_);
  }
  

  void push(particle * p);
  
  virtual unsigned int hash_function(const particle* p) const;


  Hash_shelf(utilities::Tuplef imgsz, float upb,int i_frame,bool own_part, int max_part_indx_=0, float z_offset = 0);
  


  virtual ~Hash_shelf();

  virtual void print() const;



  


  void get_region_px( int n, std::vector<const particle*> & out ,
                      float range_in_units) const;

  void get_region( int n,hash_box* box, int range) const ;
  void get_region( particle* n,hash_box* box, int range) const ;

  void get_region( int n,std::list<particle*>& out_list, int range) const ;

  void get_region( int n,std::vector<const particle*> & out_vector,int range) const ;

  int get_plane_num() const
  {
    return plane_number_;
  }
  
  
  void shelf_to_list(std::list<const particle*> &tmp) const;
  
  void shelf_to_list(std::list<particle*> &tmp) const;
  
  void compute_mean_forward_disp(utilities::Tuplef & cum_disp);
  const utilities::Tuplef& get_mean_forward_disp() const{
    return mean_forward_disp_;
  }
  const utilities::Tuplef&  get_cum_forward_disp() const{
    return cumulative_disp_;
  }

  
  void set_next_shelf(Hash_shelf* next);
  


  
  bool step_forwards(int n,const Hash_shelf* & next)const;
  
  void fill_in_neighborhood();

  void pass_fun_to_part(void(particle::*fun)());
  

  void pass_fun_to_part(void(particle::*fun)()const)const;
  
  virtual void output_to_wrapper(utilities::Wrapper_out & wrapper) const;
  
  
  void compute_corr(tracking::Corr &) const ;
  
  void compute_accum(Accumulator & in)const;
  

  void test() ;
  

  std::list<particle_track*> * shelf_to_list() const;

  const utilities::Md_store * md_store_;


private:
  std::vector<hash_box*> hash_;
  
  int plane_number_;
  
  utilities::Tuplef mean_forward_disp_;

  utilities::Tuplef cumulative_disp_;

  Hash_shelf * next_;

  unsigned int particle_count_;

  int max_part_indx_;
  
  utilities::Tuplei hash_dims_;
  utilities::Tuplei hash_cum_dim_;
  
  utilities::Tuplef img_dims_;

  
  
  float upb_;
  
  float z_offset_;
  
  
  bool own_part_;
  

  int tuple_to_indx(const utilities::Tuplei& cord)const;
  
  utilities::Tuplei indx_to_tuple(int indx) const;
  
  utilities::Tuplei range_indx_to_tuple(int indx,const utilities::Tuplei& side) const;
  
  void priv_init();

  void make_region_corners(const utilities::Tuplei & center, 
                           utilities::Tuplei & bottom_corner,
                           utilities::Tuplei & top_corner,
                           utilities::Tuplei & region_sides,
                           const int range)const;
  
  
  
};


}
#endif