ModelFunction_ThreePointCorrelation_comoving_connected.cpp

Go to the documentation of this file.

/*******************************************************************
 *  Copyright (C) 2017 by Federico Marulli                         *
 *  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.      *
 *******************************************************************/
 
#include "ModelFunction_ThreePointCorrelation.h"
 
#include "ModelFunction_ThreePointCorrelation_comoving_connected.h"
 
using namespace std;
 
using namespace cbl;
 
 
// ============================================================================================
 
 
std::vector<double> cbl::modelling::threept::zeta_RSD (const std::vector<double> theta, const std::shared_ptr<void> inputs, std::vector<double> &parameter)
{
  // structure contaning the required input data
  shared_ptr<STR_data_model_threept> pp = static_pointer_cast<STR_data_model_threept>(inputs);
 
  // input likelihood parameters
 
  // bias
  double b1 = parameter[0];
  double b2 = parameter[1];
  double bt = parameter[2];
 
  //beta
  double beta = parameter[3];
 
  vector<double> model = pp->cosmology->zeta_RSD (pp->r1, pp->r2, theta.size(), b1, b2, bt, beta, pp->rr, pp->kk, pp->Pk_matter, false, pp->max_ll, pp->use_k);
  return model;
  /*
  double theta_binSize = 1./theta.size();
 
  vector<double> xx = cosmobl::linear_bin_vector(model.size(), 0., 1.);
  cosmobl::glob::FuncGrid interp_zeta(xx, model, "Spline");
 
  vector<double> zeta(theta.size());
 
  for (size_t i=0; i<theta.size(); i++)
    zeta[i] = interp_zeta.integrate_qag(double(i)*theta_binSize, double(i+1)*theta_binSize, 1.e-4)/theta_binSize;
 
  return zeta;
  */
  //return pp->cosmology->zeta_RSD (pp->r1, pp->r2, theta.size(), b1, b2, bt, beta, pp->rr, pp->kk, pp->Pk_matter, false, pp->max_ll, pp->use_k);
}