detdist.cpp

Go to the documentation of this file.

#include "MTTensor.h"
#include <iostream>
#include <complex>

#include "Jacknife.h"
#include "Bootstrap.h"
#include "MTSampleGenerator.h"
#include "ArraySampleGenerator.h"
using namespace std;
int main()
{
        const complex<double> xx(0.0,0.0);
        const complex<double> xy(0.0555,0.0402);
        const complex<double> yx(-0.0736,-0.0695);
        const complex<double> yy(-0.0,-0.0);
        MTTensor TestData(xx,xy,yx,yy);
        
        const double noiselevel = 0.1;
        TestData.SetErrors(abs(xx)*noiselevel,abs(xy)*noiselevel,abs(yx)*noiselevel,abs(yy)*noiselevel);
        const int nsamples = 1000;
        MTSampleGenerator Generator(&MTTensor::GetRhoxy,TestData);
        gplib::rvec array(nsamples);
        
        std::generate_n(array.begin(),nsamples,Generator);
        //copy(array.begin(),array.end(),ostream_iterator<double>(cout,"\n"));
        ArraySampleGenerator ArrGen(array);
        Jacknife<MTSampleGenerator> JackEst(nsamples,Generator);
        Bootstrap<ArraySampleGenerator> ArrayEst(nsamples,ArrGen);
        Bootstrap<MTSampleGenerator> BootEst(nsamples,Generator);
        
        double jackmean, jackvar, bootmean, bootvar, arrmean, arrvar;
        JackEst.CalcErrors(jackmean,jackvar);
        BootEst.CalcErrors(bootmean,bootvar);
        ArrayEst.CalcErrors(arrmean,arrvar);
        cout << "Jm: " << jackmean << " Jv: " << jackvar << endl;
        cout <<" Bm: " << bootmean << " Bv: " << bootvar << std::endl;
        cout <<" Am: " << arrmean << " Av: " << arrvar << std::endl;
        //copy(ArrayEst.GetSamples().begin(),ArrayEst.GetSamples().end(),ostream_iterator<double>(cout,"\n"));
}

---