int farad_1(){ gROOT->SetStyle("Plain"); gStyle->SetOptFit(00001); //FARADAY john 1 TGraphErrors fj1("messdaten/faraday_j_1.txt","%lg %lg"); fj1.SetTitle("Faradayeffekt-john 1;I/A;Angel/#circ"); fj1.SetMarkerStyle(3); //fit TF1 fitfj1("Linear law","[0]+x*[1]",-5,5); fitfj1.SetLineColor(kRed); fitfj1.SetLineStyle(2); fitfj1.SetParName(0,"Y-Intersect") fitfj1.SetParName(1,"Slope") fj1.Fit(&fitfj1); double par1 = fitfj1.GetParameter(0); double err1 = fitfj1.GetParError(0); double par2 = fitfj1.GetParameter(1); double err2 = fitfj1.GetParError(1); std::cout << endl; std::cout << endl; std::cout << "Y-intersect: " << par1 << endl; std::cout << "Y-intersect error: " << err1 << endl; std::cout << "Slope " << par2 << endl; std::cout << "Slope error " << err2 << endl; std::cout << endl; std::cout << endl; //FARADAY john 2 TGraphErrors fj2("messdaten/faraday_j_2.txt","%lg %lg"); fj2.SetTitle("Faradayeffekt-john 2;I/A;Angel/#circ"); fj2.SetMarkerStyle(3); //fit TF1 fitfj2("Linear law","[0]+x*[1]",-5,5); fitfj2.SetLineColor(kRed); fitfj2.SetLineStyle(2); fj2.Fit(&fitfj2); //FARADAY smith 1 TGraphErrors fs1("messdaten/faraday_s_1.txt","%lg %lg"); fs1.SetTitle("Faradayeffekt-smith 1;I/A;Angel/#circ"); fs1.SetMarkerStyle(3); //fit TF1 fitfs1("Linear law","[0]+x*[1]",-5,5); fitfs1.SetLineColor(kRed); fitfs1.SetLineStyle(2); fs1.Fit(&fitfs1); //FARADAY smith 2 TGraphErrors fs2("messdaten/faraday_s_2.txt","%lg %lg"); fs2.SetTitle("Faradayeffekt-smith 2;I/A;Angel/#circ"); fs2.SetMarkerStyle(3); //fit TF1 fitfs2("Linear law","[0]+x*[1]",-5,5); fitfs2.SetLineColor(kRed); fitfs2.SetLineStyle(2); fs2.Fit(&fitfs2); TCanvas* c=new TCanvas("Canvas","Canvas",800,800); c->Divide(2,2); c->cd(1);fj1.DrawClone("AP"); // c->cd(1);legfj1.DrawClone("Same"); c->cd(2);fj2.DrawClone("AP"); // c->cd(2);fitfj2.DrawClone("Same"); c->cd(3);fs1.DrawClone("AP"); // c->cd(3);fitfs1.DrawClone("Same"); c->cd(4);fs2.DrawClone("AP"); // c->cd(4);fitfs2.DrawClone("Same"); c->cd(1)->SetGridx(); c->cd(1)->SetGridy(); c->cd(2)->SetGridx(); c->cd(2)->SetGridy(); c->cd(3)->SetGridx(); c->cd(3)->SetGridy(); c->cd(4)->SetGridx(); c->cd(4)->SetGridy(); c->Print("Faradayeffekt.pdf"); }