Triplet1D.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2015 by Federico Marulli, Michele Moresco,        *
 *  and Alfonso Veropalumbo                                         *
 *                                                                  *
 *  federico.marulli3@unibo.it                                      *
 *                                                                  *
 *  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.,                                      *
 *  59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.       *
 *******************************************************************/
 
#ifndef __TRIPLET1D__
#define __TRIPLET1D__
 
 
#include "Triplet.h"
 
 
// ===================================================================================================
 
 
namespace cbl {
  
  namespace triplets {
    
    class Triplet1D : public Triplet {
 
    private:
      
  
      virtual void set_parameters () = 0;
  
 
      
    protected:
 
      std::vector<double> m_scale;
      
      std::vector<double> m_TT1D;
 
    
 
      double m_r12;
 
      double m_r12_binSize;
 
      double m_r13;
   
      double m_r13_binSize;   
 
      int m_nbins;
 
      double m_binSize;
      
 
  
    public:
    
 
      Triplet1D () 
    : m_r12(-1), m_r12_binSize(-1), m_r13(-1), m_r13_binSize(-1), m_nbins(-1), m_binSize(-1)
    { m_tripletDim = Dim::_1D_; }
                                                                
      Triplet1D (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int nbins) 
    : m_r12(r12), m_r12_binSize(r12_binSize), m_r13(r13), m_r13_binSize(r13_binSize) , m_nbins(nbins)
    { m_tripletDim = Dim::_1D_; }
 
      ~Triplet1D () = default;
      
 
    
 
      double scale (const int i) const override { return m_scale[i]; }
 
      std::vector<double> scale () const override { return m_scale; }
      
      double TT1D (const int i) const override { return m_TT1D[i]; }
 
      std::vector<double> TT1D () const override { return m_TT1D; }
      
      double r12 () const override { return m_r12; }
      
      double r12_binSize () const override { return m_r12_binSize; }
 
      double r13 () const override { return m_r13;  }
 
      double r13_binSize () const override { return m_r13_binSize; }
 
      int nbins () const override { return m_nbins; }
      
      double binSize () const override { return m_binSize; }
 
 
      
 
      void set_TT1D (const int i, const double tt) { checkDim(m_TT1D, i, "m_TT1D"); m_TT1D[i] = tt; }
 
      void add_TT1D (const int i, const double tt) { checkDim(m_TT1D, i, "m_TT1D"); m_TT1D[i] += tt; }
      
 
      
    
      void Sum (const std::shared_ptr<Triplet> tt, const double ww=1) override;
 
    
    };
 
    
    // ============================================================================================
    // ============================================================================================
 
    
    class Triplet1D_comoving : public Triplet1D {
 
    private:
      
  
      virtual void set_parameters () = 0;
  
 
      
    public:
    
 
      Triplet1D_comoving () = default;
                                                                
      Triplet1D_comoving (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int nbins) 
    : Triplet1D(r12, r12_binSize, r13, r13_binSize, nbins) {}
      
      virtual ~Triplet1D_comoving () = default;
 
    
    };
 
    
    // ============================================================================================
    // ============================================================================================
 
 
    class Triplet1D_comoving_theta : public Triplet1D_comoving {
      
    private:
 
  
      void set_parameters () override;
      
 
      
    public:
    
 
      Triplet1D_comoving_theta () 
    { m_tripletType = TripletType::_comoving_theta_; }
                                                                
      Triplet1D_comoving_theta (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int nbins) 
    : Triplet1D_comoving(r12, r12_binSize, r13, r13_binSize, nbins)
    { m_tripletType = TripletType::_comoving_theta_; set_parameters(); m_scale.resize(m_nbins+1, 0.); m_TT1D.resize(nbins+1, 0.); }
 
 
    
     
      void get_triplet (const double r12, const double r13, const double r23, int &klin) override;   
 
      void set_triplet (const int klin, const double ww=1.) override;
 
      void put (const double r12, const double r13, const double r23, const double ww=1.) override;
    
      void put (const std::shared_ptr<catalogue::Object> obj1, const std::shared_ptr<catalogue::Object> obj2, const std::shared_ptr<catalogue::Object> obj3) override;
 
    
    };
 
    // ============================================================================================
    // ============================================================================================
 
 
    class Triplet1D_multipoles_direct : public Triplet1D_comoving {
      
    private:
 
  
      void set_parameters () override;
      
 
      
    public:
    
 
      Triplet1D_multipoles_direct () 
    { m_tripletType = TripletType::_multipoles_direct_; }
                                                                
      Triplet1D_multipoles_direct (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int norders) 
    : Triplet1D_comoving(r12, r12_binSize, r13, r13_binSize, norders)
    { m_tripletType = TripletType::_multipoles_direct_; set_parameters(); m_scale.resize(m_nbins+1, 0.); m_TT1D.resize(m_nbins+1, 0.); }
 
 
    
 
      void get_triplet (const double r12, const double r13, const double r23, int &klin)    
      { (void)r12; (void)r13; (void)r23; (void)klin; ErrorCBL("", "get_triplet", "Triplet1D.h!"); }
 
      void set_triplet (const int klin, const double ww=1.)   
      { (void)klin; (void)ww; cbl::ErrorCBL("", "get_triplet", "Triplet1D.h!"); }
   
      void put (const double r12, const double r13, const double r23, const double ww=1.) override;
    
      void put (const std::shared_ptr<catalogue::Object> obj1, const std::shared_ptr<catalogue::Object> obj2, const std::shared_ptr<catalogue::Object> obj3) override;
 
    
    };
 
 
    // ============================================================================================
    // ============================================================================================
 
 
    class Triplet1D_comoving_side : public Triplet1D_comoving {
 
    private:
 
      double m_min;
 
      double m_max;
 
  
      void set_parameters () override;
 
 
      
    public:
      
 
      Triplet1D_comoving_side () 
      { m_tripletType = TripletType::_comoving_side_; }
 
      Triplet1D_comoving_side (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int nbins) 
    : Triplet1D_comoving(r12, r12_binSize, r13, r13_binSize, nbins)
      { m_tripletType = TripletType::_comoving_side_; set_parameters(); m_scale.resize(m_nbins+1, 0.); m_TT1D.resize(nbins+1, 0.); }
 
 
      
     
      void get_triplet (const double r12, const double r13, const double r23, int &klin) override;   
 
      void set_triplet (const int klin, const double ww=1.) override;
 
      void put (const double r12, const double r13, const double r23, const double ww=1.) override;
    
      void put (const std::shared_ptr<catalogue::Object> obj1, const std::shared_ptr<catalogue::Object> obj2, const std::shared_ptr<catalogue::Object> obj3) override;
 
      
    };
 
    
    // ============================================================================================
    // ============================================================================================
 
 
    class Triplet1D_comoving_costheta : public Triplet1D_comoving {
      
    private:
 
  
      void set_parameters () override;
      
 
      
    public:
    
 
      Triplet1D_comoving_costheta () 
    { m_tripletType = TripletType::_comoving_costheta_; }
                                                                
      Triplet1D_comoving_costheta (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int nbins) 
    : Triplet1D_comoving(r12, r12_binSize, r13, r13_binSize, nbins)
    { m_tripletType = TripletType::_comoving_costheta_; set_parameters(); m_scale.resize(m_nbins+1, 0.); m_TT1D.resize(m_nbins+1, 0.); }
 
 
    
     
      void get_triplet (const double r12, const double r13, const double r23, int &klin) override;   
 
      void set_triplet (const int klin, const double ww=1.) override;
 
      void put (const double r12, const double r13, const double r23, const double ww=1.) override;
    
      void put (const std::shared_ptr<catalogue::Object> obj1, const std::shared_ptr<catalogue::Object> obj2, const std::shared_ptr<catalogue::Object> obj3) override;
 
    
    };
 
    
    // ============================================================================================
    // ============================================================================================
 
 
    class Triplet1D_angular : public Triplet1D {
 
    private:
 
  
      virtual void set_parameters () = 0;
 
 
      
    public:
    
 
      Triplet1D_angular () = default;
           
      Triplet1D_angular (const double r12, const double r12_binSize, const double r13, const double r13_binSize, const int nbins) 
    : Triplet1D(r12, r12_binSize, r13, r13_binSize, nbins) {}
 
      virtual ~Triplet1D_angular () = default;
 
    
    };
 
  }
}
 
#endif