Hi,
I want to do a multidimenional extended ML fit.
I have here a small snippet of code:
[code]{
gSystem->Load(“libRooFit.so”);
RooRealVar mBC(“mBC”,“mBC”,1.8041,1.8865);
RooRealVar mBCtag(“mBCtag”,“mBCtag”,1.8401,1.8865);
RooRealVar var3(“var3”,“var3”,150,100,200);
RooRealVar var4(“var4”,“var4”,150,100,200);
RooRealVar var5(“var5”,“var5”,15,10,20);
RooRealVar var6(“var6”,“var6”,20,18,22);
RooRealVar var7(“var7”,“var7”,1,10,10);
RooRealVar CB2_m0(“CB2_m0”,“effSignal_mBC_CB2_m0”,1.865,1.86,1.87);
RooRealVar CB2_sigma(“CB2_sigma”,“effSignal_mBC_CB2_sigma”,0.001,0.0005,0.010);
RooGaussian gCB2(“gCB2”,“gCB2”,mBC,CB2_m0,CB2_sigma);
RooRealVar CB3_m0(“CB3_m0”,“effSignal_mBC_CB3_m0”,1.865,1.86,1.87);
RooRealVar CB3_sigma(“CB3_sigma”,“effSignal_mBC_CB3_sigma”,0.001,0.0005,0.010);
RooGaussian gCB3(“gCB3”,“gCB3”,mBCtag,CB3_m0,CB3_sigma);
RooPolynomial pol1(“pol1”,“pol1”,mBC,var5);
RooPolynomial pol2(“pol2”,“pol2”,mBCtag,var6);
RooProdPdf bkg(“bkg”,“bkg”,pol1,pol2);
RooProdPdf sig(“sig”,“sig”,gCB3,gCB2);
RooExtendPdf esig(“esig”,“esig”,sig,var3);
RooExtendPdf ebkg(“ebkg”,“ebkg”,bkg,var4);
RooRealVar f(“f”,“f”,0.5,0,1);
RooAddPdf tot(“tot”,“tot”,esig,ebkg);
//RooAddPdf tot(“tot”,“tot”,RooArgSet(sig,bkg),RooArgSet(var3,var4));
RooAddPdf tot2(“tot2”,“tot2”,sig,bkg,f);
//RooDataSet* data = tot.generate(RooArgSet(mBC,mBCtag),1000);//ERROR
RooDataSet* data = tot2.generate(RooArgSet(mBC,mBCtag),1000);
RooFitResult* res = tot.fitTo(data,Save());
res->Print(“v”);
TCanvas x;
RooPlot x1 = mBC.frame();
data->plotOn(x1);
tot.plotOn(x1); //ERROR
x1->Draw();[/code]
In case i built my extended PDF like this:
RooExtendPdf esig("esig","esig",sig,var3);
RooExtendPdf ebkg("ebkg","ebkg",bkg,var4);
RooAddPdf tot("tot","tot",esig,ebkg);I get an error when I try to plot this pdf: (generate gives also the same error)
[quote]
[#0] ERROR:InputArguments – RooAddPdf::checkObservables(tot): ERROR: coefficient ebkg and PDF esig have one or more dependents in common
[#0] ERROR:Plotting – RooAbsReal::plotOn(tot) error in checkObservables, abort[/quote]
When I just do:
Everything works fine.
Can someone explain this to me?
I’m using ROOT 5.34/03 and RooFit v3.55
Thanks in advance!
Peter