Double_t myfunc(Double_t *x, Double_t *par)
{
	Float_t phi =x[0]/180*TMath::Pi();
	Float_t dh =par[0];
	Float_t din =par[1];
	Float_t nc =par[2];
	Float_t nt =par[3];
	Float_t L =par[4];
	Float_t index =par[5];
	Float_t Nteff;Float_t Nceff;Float_t dD;Float_t dphi;
	Float_t z=4;Float_t Cd =209;
	Nteff=nt*L/cos(phi)*index;Nceff=nc*L/cos(phi)*index;
	dD=Cd*sqrt(z)/cos(phi)/sqrt(Nteff);
	dphi=L*1e4/sqrt(12)*tan(phi)/sqrt(Nceff);
	Double_t f = sqrt(pow(dh,2)+pow(dD,2)+pow(dphi,2));
	return f;
}
Double_t myfunc1(Double_t *x, Double_t *par)
{
	Float_t phi =x[0]/180*TMath::Pi();
	Float_t s =par[0];
	Float_t L =par[1];
	Float_t ph =par[2]/180*TMath::Pi();
	Float_t nc =par[3];
	Float_t ifc =par[4];
	Double_t f = sqrt(pow(s,2)+pow(L,2)/12*pow(tan(phi-ph),2)/pow(nc,ifc));
	return f;
}


void posrebs()
{

	TF1 *f1 = new TF1("f1",myfunc,-30,30,6);
	f1->SetParameter(0,50);
	f1->SetParameter(1,218);
	f1->SetParameter(2,24.3);
	f1->SetParameter(3,94);
	f1->SetParameter(4,0.6);
	f1->SetParameter(5,0.4);
	f1->SetTitle("");
	f1->SetLineColor(1);
	f1->GetXaxis()->SetTitle("Angle [deg]");
	f1->GetYaxis()->SetTitle("#sigma [#mum]");
	f1->Draw("");
	TF1 *f2 = new TF1("f2",myfunc1,-30,30,5);
	f2->SetParameter(0,0.131);
	f2->SetParameter(1,6);
	f2->SetParameter(2,0.1);
	f2->SetParameter(3,24.3);
	f2->SetParameter(4,0.64);
	f2->SetLineColor(2);
	f2->Draw("same");


}