// Hyperfeinstruktur Kalibration

void hyperfeinmess(){

    gROOT->Reset();
    gROOT->SetStyle("Plain");
    gStyle->SetOptFit();

    // Anzahl der Kanäle
    const int rows = 1024;


    // Kalibrationswerte

    double amp = 3607.7;
    double amperr = 1.53348;
    double phase = 0.0106183;
    double phaseerr = 0.000184971;

    double lambda = 632.8e-9;
    double enull = 14.39e3;



    // Daten öffnen
    ifstream in;

    in.open("data/moessbauer/hyperfein.ws5");

    //Daten einlesen
    double energydiff[rows];
    double data[rows];

    for (int i=0; i<rows; i++){
        in >> data[i];
        energydiff[i] = -1/TMath::C() * amp*cos( (6.135923152e-3) * i + phase) * lambda * rows/(2*180) * enull;
        //energydiff[i] = i;
        data[i] = (0-data[i])+380e3;

        //cout << energydiff[i] << "  " << data[i] << endl;
    }

    //Daten schließen
    in.close();


    //Graphen
    TCanvas *canvas = new TCanvas("canvas", "Canvas", 800, 400);
    canvas->cd();

    TGraph *messgraph = new TGraph( rows, energydiff, data);
    messgraph->SetMarkerStyle(6);

    TF1 *fit1 = new TF1("Fit1", "gaus", (2.77606e-008), (4.09103e-008));
    messgraph->Fit(fit1,"RM");






    messgraph->Draw("AP");
    canvas->Update();


}