#include <math.h>

#include "cor_const.h"

void kaiser(){
   gROOT->SetStyle("Plain");
   gStyle->SetOptStat(0);
   gStyle->SetOptFit(1111);
   gStyle->SetPalette(1);

   double cpibin, csigmabin;
   int value;

   TCanvas *wake = new TCanvas("c","kaiser plot",0,0,800,800); 
   TH2F *h = new TH2F("h","kaiser",kbins+1,-10,10,(kbins)/2. + 1,0,10);

   int n = 0;
   ifstream widthfile("kaiser.txt");
   while((widthfile >> cpibin >> csigmabin >> value)){
      //cout<<cpibin<<" "<<csigmabin<<" "<<value<< endl;
      h->Fill(cpibin,csigmabin,value);
      n++;
   }  
   widthfile.close();

   h->SetTitle(""); 
   h->GetXaxis()->SetTitle("s_{#parallel} [h^{-1}Mpc]");
   h->GetYaxis()->SetTitle("s_{#perp} [h^{-1}Mpc]");
   h->GetYaxis()->SetTitleOffset(1.2);

   h->Draw("COLZ");
}