errordist.cpp

Go to the documentation of this file.

#include "MTTensor.h"
#include <fstream>
#include <iostream>
#include <complex>
#include "MTSampleGenerator.h"
#include "Jacknife.h"

using namespace std;
int main()
{
        const complex<double> xx(1.405,2.23);
        const complex<double> xy(5.33,2.5);
        const complex<double> yx(-7.45,-4.23);
        const complex<double> yy(-1.45,-3.29);
        const MTTensor TestData(xx,xy,yx,yy);
        const size_t ntestcases = 100000;
        
        cerr << TestData.GetPhi11() << " " <<  TestData.GetPhi12() << " " <<  TestData.GetPhi21() << " " <<  TestData.GetPhi22() << endl;
        cerr << TestData.GetPhiMax() << " " << TestData.GetPhiMin() <<  " " << TestData.GetAlpha_phi() <<  " " << TestData.GetBeta_phi() << endl;
        cerr << TestData.GetPhixy() << endl;
        cerr << "Det Phi: " << TestData.GetdetPhi() << endl;
        
        const double maxnoise = 0.5;
        const double minnoise = 0.01;
        const double noisestep = 0.01;
        for (double noiselevel = minnoise; noiselevel <= maxnoise; noiselevel += noisestep)
        {
                MTTensor DistData(TestData);
                DistData.SetErrors(abs(TestData.GetZxx())*noiselevel,abs(TestData.GetZxy())*noiselevel,
                                        abs(TestData.GetZyx())*noiselevel,abs(TestData.GetZyy())*noiselevel);
                MTSampleGenerator Generator(&MTTensor::GetdetPhi,DistData);
                Jacknife<MTSampleGenerator> ErrEst(ntestcases,Generator);
                double currmean;
                double currvar;
                ErrEst.CalcErrors(currmean,currvar);
                //copy(ErrEst.GetSamples().begin(),ErrEst.GetSamples().end(),ostream_iterator<double>(cout,"\n"));
                cout << noiselevel << " " << currmean -TestData.GetdetPhi() << " " << sqrt(currvar)/TestData.GetdetPhi() << endl;
        }
}

---