void test(void){
	
	using namespace RooFit;
	using namespace RooStats;
	using namespace std;

	double XMIN = -5.;
	double XMAX = +5.;
	double XMIN2 = -0.;
	double XMAX2 = +5.;

	double par_mean  = 0.;
	double par_sigma = 1.;
	double par_norm  = 10.;

	RooRealVar x("x","x",XMIN,XMAX);
	x.setBins((int)XMAX-(int)XMIN);

	x.setRange("fullrange",XMIN,XMAX);
	x.setRange("testrange",XMIN2,XMAX2);

	RooRealVar mg("mg","mg",par_mean);
	RooRealVar sg("sg","sg",par_sigma);
	RooRealVar ng("ng","ng",par_norm);

	RooGaussian gaus("gaus","gaus",x,mg,sg);
	RooExtendPdf egaus("egaus","egaus",gaus,ng);

	TCanvas* c1 = new TCanvas("c1","c1");
	RooPlot* rp = x.frame(XMIN,XMAX);
	egaus.plotOn(rp,Range("fullrange"),Normalization(1.,RooAbsReal::RelativeExpected));
	rp->GetXaxis()->SetRangeUser(XMIN2,XMAX2);
	rp->Draw();

	double fraction = egaus.createIntegral(x,NormSet(x),Range("testrange"))->getVal();

	RooArgList pars(* egaus.getParameters(RooArgSet(x)));
	RooArgSet prodSet(egaus);
	RooProduct unNormPdf("gauss model","gauss model",prodSet);
	TF1* root_egaus = unNormPdf.asTF(RooArgList(x),pars);
	double integ = root_egaus->Integral(XMIN,XMAX);
	double integ_sub = root_egaus->Integral(XMIN2,XMAX2);

	cout<<"createIntegral "<<fraction<<endl;
	cout<<"TF1 conversion "<<integ_sub<<"  "<<integ<<" => "<<integ_sub/integ<<endl;

	//----------------------------------------------------------------------------

	RooRealVar x2("x2","x2",XMIN2,XMAX2);
	x2.setBins((int)XMAX2-(int)XMIN2);
	
	x2.setRange("testrange",XMIN2,XMAX2);

	RooRealVar mg2("mg2","mg2",par_mean);
	RooRealVar sg2("sg2","sg2",par_sigma);
	RooRealVar ng2("ng2","ng2",par_norm*fraction);

	RooGaussian gaus2("gaus2","gaus2",x2,mg2,sg2);
	RooExtendPdf egaus2("egaus2","egaus2",gaus2,ng2);

	TCanvas* c2 = new TCanvas("c2","c2");
	RooPlot* rp2 = x2.frame(XMIN2,XMAX2);
	egaus2.plotOn(rp2,Range("testrange"),LineColor(kRed),LineStyle(2),
				  Normalization(1.0,RooAbsReal::RelativeExpected));
	rp2->Draw();

}