void costheta_resonance(TString input=""){

	TLorentzVector * PLamData = new TLorentzVector();
	TLorentzVector * PKData = new TLorentzVector();

	TLorentzVector * PLamFit = new TLorentzVector();
	TLorentzVector * PKFit = new TLorentzVector();

	TLorentzVector PLamKFit, PLamKData;


	TFile * infile = new TFile(input, "READ");

	TH1D * h_data = (TH1D*) infile->Get("DataThetaKm_HeliLambdaKm");
	TH1D * h_fit = (TH1D*) infile->Get("FitThetaKm_HeliLambdaKm");


	TCanvas * c = new TCanvas("c", "costheta", 0,0,800,500);

	h_data->SetLineColor(kBlue);
	h_fit->SetLineColor(kRed);

	TLegend * l = new TLegend(0.745,0.727,0.927,0.886);

	l->AddEntry(h_data ,"Data", "lep");
	l->AddEntry(h_fit ,"Fit");

	h_data->Draw();
	h_fit->Draw("SAME HIST");
	l->Draw("SAME");



}