#ifndef __CINT__
#include "RooGlobalFunc.h"
#endif
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooGaussian.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "RooPlot.h"
#include "RooFitResult.h"
#include "RooGenericPdf.h"
#include "RooConstVar.h"
#include "RooDataHist.h"
#include "RooHistPdf.h"
#include "RooFitResult.h"
#include "RooArgSet.h"
#include "RooAddPdf.h"
#include "RooCategory.h"
#include "RooSimultaneous.h"
#include "RooChebychev.h"
#include "RooAbsReal.h"
#include "RooMCStudy.h"
#include "RooMinuit.h"
#include "RooNLLVar.h"
#include "RooTrace.h"
#include "RooAddition.h"
#include "RooProdPdf.h"
#include "RooGaussian.h"
#include "RooFFTConvPdf.h"

#include <string>
#include <iostream>
#include <sstream>

void test(){
  
  // Define observable of interest
  RooRealVar signal ("signal", "signal", -10.0, 10.0);

  // Define/Create first Gaussian
  RooRealVar mean1 ("mean1", "mean1", 1.0, -10.0, 10.0);
  RooRealVar sig1 ("sig1", "sig1", 2.0, -10.0, 10.0);
  RooGaussian g1 ("g1", "g1", signal, mean1, sig1);

  // Define/Create second Gaussian
  RooRealVar mean2 ("mean2", "mean2", 5.0, -10.0, 10.0);
  RooRealVar sig2 ("sig2", "sig2", 2.0, -10.0, 10.0);
  RooGaussian g2 ("g2", "g2", signal, mean2, sig2);

  // Define/Create Gaussian to be use in convolution
  RooRealVar meanc ("meanc", "meanc", 0.0, -10.0, 10.0);
  RooRealVar sigc ("sigc", "sigc", 2.0, -10.0, 10.0);
  RooGaussian gconv ("gconv", "gconv", signal, meanc, sigc);


  // Vector to store g1 and g2 before convolution.
  // The size of the vector is used in the 'for' loop below
  std::vector<RooGaussian> vec_gauss;
  vec_gauss.push_back(g1);
  vec_gauss.push_back(g2);

  // Create RooArgList for the g1 and g2 pdfs after convolution
  RooArgList convlist;

  // Loops through vec_gauss, applying convolution to each Gaussian.
  // Adds them to convlist after convolution.
  for(int i=1;i<vec_gauss.size()+1;i++)
    {
      std::ostringstream ss; //convert for loop counter to string
      ss << i;
      
      RooFFTConvPdf temp (("t"+ss.str()).c_str(), ("t"+ss.str()).c_str(), signal, vec_gauss[i-1], gconv);
      
      convlist.add(temp);
      convlist.Print("v");      
    }

  convlist.Print("s");

  // Define RooRealVars to be associated with the Gaussians. (n1 with g1 ; n2 with g2)
  RooRealVar n1 ("n1", "n1", 1.0, 0.0, 10.0);
  RooRealVar n2 ("n2", "n2", 1.0, 0.0, 10.0);
  RooArgList fitvars(n1,n2);

  // Add pdfs together to form a combined model
  RooAddPdf model("model", "model", convlist, fitvars);

}