geom_2d.h

// -*- C++ -*-
// File: geom_2d.h                                                           //
// Description: header file for two-dimensional geometry tools               //
// This file is part of the SISCone project.                                 //
// WARNING: this is not the main SISCone trunk but                           //
//          an adaptation to spherical coordinates                           //
// For more details, see http://projects.hepforge.org/siscone                //
//                                                                           //
// Copyright (c) 2006-2008 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 __SPH_GEOM_2D_H__
#define __SPH_GEOM_2D_H__
 
#include <iostream>
#include <math.h>
#include <siscone/defines.h>
#include <siscone/geom_2d.h>
 
#ifndef M_PI
#define M_PI 3.141592653589793238462643383279502884197
#endif
 
namespace siscone_spherical{
 
class CSphtheta_phi_range{
public:
  CSphtheta_phi_range();
 
  CSphtheta_phi_range(double c_theta, double c_phi, double R);
 
  CSphtheta_phi_range& operator = (const CSphtheta_phi_range &r);
 
  int add_particle(const double theta, const double phi);
 
  unsigned int theta_range;     
 
  unsigned int phi_range;     
 
  static double theta_min;  
  static double theta_max;  
 
private:
  inline unsigned int get_theta_cell(double theta){
    if (theta>=theta_max) return 1u<<31;
    return (unsigned int) (1u << ((int) (32*((theta-theta_min)/(theta_max-theta_min)))));
  }
 
  inline unsigned int get_phi_cell(double phi){
    return (unsigned int) (1u << ((int) (32*phi/twopi+16)%32));
  }
};
 
bool is_range_overlap(const CSphtheta_phi_range &r1, const CSphtheta_phi_range &r2);
 
const CSphtheta_phi_range range_union(const CSphtheta_phi_range &r1, const CSphtheta_phi_range &r2);
 
}
 
#endif