void AccumulatedAsymmetrieAlt(double xmin = 300., double xmax = 10000., const char *filename = "fp13.root") {

  double X1_statbox = 0.6;
  double X2_statbox = 0.95;
  double Y1_statbox = 0.65;
  double Y2_statbox = 0.95;

//---------------------------------------------------------------------//
// Laden der Histogramme                                               //
//---------------------------------------------------------------------//

    gStyle->SetOptFit(111);
    gStyle->SetOptStat("ni");
    gStyle->SetStatFormat("g");
    gStyle->SetMarkerStyle(20);
    gStyle->SetMarkerSize(0.5);
    gStyle->SetHistLineWidth(1);
    gROOT->ForceStyle();

    /* ---> Gruppe Bachmann/Haas/Mozer   */

    TFile *f1 = new TFile("ResultsAll.root");
    TH1D *hAsymSumAdd = new TH1D(*hAsymoSumAll);
    TH1D *hAsymDiffAdd = new TH1D(*hAsymoDiffAll);

    TFile *f2 = new TFile("Results.root");

    hAsymSumAdd->Add(ho);
    hAsymSumAdd->Add(hmo);
    hAsymDiffAdd->Add(ho,1.);
    hAsymDiffAdd->Add(hmo,-1.);
    hAsymSumAdd->Add(hu);
    hAsymSumAdd->Add(hmu);
    hAsymDiffAdd->Add(hu,-1.);
    hAsymDiffAdd->Add(hmu,1.);

    TH1D *hAsymAdd = new TH1D(*hAsymDiffAdd);
    hAsymAdd->Divide(hAsymSumAdd); 
    hAsymAdd->SetName("Alle Gruppen: Asymmetrie mit/ohne B-Feld");
    hAsymAdd->SetTitle("Alle Gruppen: Asymmetrie mit/ohneB-Feld");
    hAsymAdd->SetAxisRange(0.,8000.); 
    hAsymAdd->SetMinimum(-0.3);
    hAsymAdd->SetMaximum(0.3);
    hAsymAdd->Draw(); 

    TF1 *fitFuncAsym = new TF1("fitFuncAsym", fitFuncAsym, 0., 0., 4);  // ( Name, Titel, xminA, xmaxA, Anzahl der Parameter) 
    fitFuncAsym->SetParameters(0.04, 4.e-3, 3.14159, 0.0);           // Startwerte fuer Parameter festlegen
    fitFuncAsym->SetParName(0, "P*A");                             // Parameter benennen
    fitFuncAsym->SetParName(1, "#omega");
    fitFuncAsym->SetParName(2, "#varphi");
    fitFuncAsym->SetParName(3, "const");
    //    fit_mo->SetParLimits(0, 0., 1.);
    hAsymAdd->Fit(fitFuncAsym, "", "em", xmin, xmax);
    hAsymAdd->Draw(); 
}

//---------------------------------------------------------------------//
// definieren der Fitfunktion für Asymmetrie                          //
//---------------------------------------------------------------------//

double fitFuncAsym(double *x, double *par) {
    double value = 0.5*par[0] * TMath::Cos(par[1] * x[0] + par[2]) + par[3];
    return value;
}