track_list.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.
#include <list>
#include <set>
#include <vector>
#include "particle_track.h"
//#include "track_shelf.h"

#ifndef TRACK_LIST
#define TRACK_LIST
namespace tracking{

class Track_shelf;
class track_list{

public:

  void link_next(std::list<particle_track*>* new_next);

  track_list(std::list<particle_track*>* first_list,
             float i_max_disp,Track_shelf* in_tracks);

  ~track_list();
  
protected:
  std::list<particle_track*> * p_list;
  std::list<particle_track*> * n_list;
  
  std::list<particle_track*>* store_list;
  
  //std::set<particle_track*,bool(*)(particle_track*,particle_track*)> p_sub_net;
  std::set<particle_track*> p_sub_net;
  //std::set<particle_track*,bool(*)(particle_track*,particle_track*)> n_sub_net;
  std::set<particle_track*> n_sub_net;
  Track_shelf* tracks;

  float max_disp;
  float max_disp_sq;

  
  void link_sub_ntwrk();
  
  
  void trivial_bonds();

  void set_up_sub_ntwrk();

  void clean_lists();
  
  std::list<particle_track*>::iterator erase(std::list<particle_track*>* in,
                                        std::list<particle_track*>::iterator it);
  
  static bool part_lt(particle_track* a, particle_track* b);
  
  void recur_fnc_np(std::vector<std::pair<particle_track*, particle_track*> >* min,
                    std::vector<std::pair<particle_track*, particle_track*> >* cur,
                    float disp,float& min_disp,
                    std::set<particle_track*>::iterator it,
                    const std::set<particle_track*>::iterator& itend);
  void recur_fnc_pn(std::vector<std::pair<particle_track*, particle_track*> >* min,
                    std::vector<std::pair<particle_track*, particle_track*> >* cur,
                    float disp,float & min_disp,
                    std::set<particle_track*>::iterator it,
                    const std::set<particle_track*>::iterator& itend);
  void  link_pairs(std::vector<std::pair<particle_track*, particle_track*> >& in);
  
};

// tac 2009-10-19
// took this out in September because I got rid of the unq_id
// inline bool track_list::part_lt(particle_track* a, particle_track* b){
//   return (a->unq_id_)<(b->unq_id_);
//}
}

#endif