quadtree.h

// -*- C++ -*-
// File: quadtree.h                                                          //
// Description: header file for quadtree management (Cquadtree class)        //
// This file is part of the SISCone project.                                 //
// For more details, see http://projects.hepforge.org/siscone                //
//                                                                           //
// Copyright (c) 2006 Gavin Salam and Gregory Soyez                          //
//                                                                           //
// 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 2 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, write to the Free Software               //
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA //
//                                                                           //
// $Revision::                                                              $//
// $Date::                                                                  $//
 
#ifndef __QUADTREE_H__
#define __QUADTREE_H__
 
#include "momentum.h"
#include <stdio.h>
 
namespace siscone{
 
class Cquadtree{
 public:
  Cquadtree();
 
  Cquadtree(double _x, double _y, double _half_size_x, double _half_size_y);
 
  ~Cquadtree();
 
  int init(double _x, double _y, double _half_size_x, double _half_size_y);
 
  int add(Cmomentum *v_add);
 
  Creference circle_intersect(double cx, double cy, double cR2);
 
  int save(FILE *flux);
 
  int save_leaves(FILE *flux);
 
  double centre_x;           
  double centre_y;           
  double half_size_x;        
  double half_size_y;        
 
  Cmomentum *v;              
 
  Cquadtree* children[2][2]; 
  bool has_child;            
};
 
}
#endif