vicinity.h

// -*- C++ -*-
// File: vicinity.h                                                          //
// Description: header file for particle vicinity (Cvicinity 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 __VICINITY_H__
#define __VICINITY_H__
 
#include <vector>
#include <list>
#include "momentum.h"
#include "defines.h"
#include "quadtree.h"
 
namespace siscone{
 
  
 
class Cvicinity_inclusion {
public:
  Cvicinity_inclusion() : cone(false), cocirc(false) {}
 
  bool cone;    
  bool cocirc;  
};
 
 
class Cvicinity_elm{
 public:
  Cmomentum *v;
 
  Cvicinity_inclusion *is_inside;   
 
  // centre variables
  double eta;              
  double phi;              
  double angle;            
  bool side;               
  double cocircular_range; 
 
  std::list<Cvicinity_elm * > cocircular;                                          
};
 
bool ve_less(Cvicinity_elm *ve1, Cvicinity_elm *ve2);
 
 
class Cvicinity{
 public:
  Cvicinity();
 
  Cvicinity(std::vector<Cmomentum> &_particle_list);
 
  ~Cvicinity();
 
  void set_particle_list(std::vector<Cmomentum> &_particle_list);
 
  void build(Cmomentum *_parent, double _VR);
 
  // cone kinematical information
  Cmomentum *parent;         
  double VR;                 
  double VR2;                
  double R;                  
  double R2;                 
  double inv_R_EPS_COCIRC;   
  double inv_R_2EPS_COCIRC;  
 
  // particle list information
  int n_part;                                 
  std::vector<Cmomentum> plist;               
  std::vector<Cvicinity_inclusion> pincluded; 
  Cvicinity_elm *ve_list;                     
#ifdef USE_QUADTREE_FOR_STABILITY_TEST
  Cquadtree *quadtree;                        
#endif
 
  // vicinity information
  std::vector<Cvicinity_elm*> vicinity;       
  unsigned int vicinity_size;                 
 
 protected:
  void append_to_vicinity(Cmomentum *v);
 
  // internal variables
  double pcx;    
  double pcy;    
};
 
}
 
#endif