Dear rooters
I am using vector to contain some pdfs and trying to add them afterwards using RooAddPdf. When I draw the RooAddPdf, the code crashes. If I do not use the vector to contain the pdfs, the code works fine. Do you know why the vector do not woks here? The number of pdfs is not constant when I have the RooAddPdf, so using vector here can make my life much easier. following are the code, the RooPdfs is a customer built pdf here. Do you have any ideas?
std::vector<RooAbsPdf*> single_pdfs;
for(int i=0; i<resonance_num; i++) {
RooRealVar opt_s(“opt_s”,“opt_s”,i);
RooAbsPdf* single_sig=new RooPdfs(pdfs_name[i], pdfs_title[i], pipi, kk, plane, f0mass, f0mean, gpi, factor, As500, A0, Apara, Aperp, D0, Dpara, Dperp, D1430, Dpara1430, Dperp1430, As1370, mean1430, width1430, UP, opt_s);
single_pdfs.push_back(single_sig);
}
std::vector<RooRealVar*> single_nums;
double sig_nums[7] = {num_980, num_1270, num_1430, num_500, num_770, num_1370, num_1370};
for(int i=0; i<resonance_num; i++) {
RooRealVar* Num_s = new RooRealVar(“Num_s”, “Num_s”, sig_nums[i]);
single_nums.push_back(Num_s);
}
RooArgList* pdfListss = new RooArgList();//“pdfLists”);
RooArgList* numListss = new RooArgList();//“numLists”);
for(int i=0; i<resonance_num; i++) {
pdfListss->add((single_pdfs[i]));
numListss->add((single_nums[i]));
}
RooAddPdf* Hel1 = new RooAddPdf(“All”, “All”, *pdfListss, *numListss);
Then I want to draw the Hel1, using
RooPlot* f0massplot = f0mass.frame(Range(f0dn, f0up),Title(""),Bins(20));
f0massplot -> GetXaxis()-> SetTitle(“m_{#pi^{+}#pi^{-}} MeV/c^{2}”);
f0massplot ->GetYaxis()->SetTitle(“Events / 20 MeV/c^{2}”);
data->plotOn(f0massplot);
Hel1 -> plotOn(f0massplot,LineColor(kBlue));
The code crashes
[#1] INFO:InputArguments – RooAbsData::plotOn(data) INFO: dataset has non-integer weights, auto-selecting SumW2 errors instead of Poisson errors
[#1] INFO:Plotting – RooAbsReal::plotOn(All) plot on f0_mass_GeV integrates over variables (pipi_hel,kk_hel,plane_ang)
*** Break *** segmentation violation
===========================================================
There was a crash.
This is the entire stack trace of all threads:
#0 0x00002b22a3b8c075 in waitpid () from /lib64/libc.so.6
#1 0x00002b22a3b2e741 in do_system () from /lib64/libc.so.6
#2 0x00002b22a1851a02 in TUnixSystem::StackTrace() () from /afs/cern.ch/sw/lcg/app/releases/ROOT/5.27.06/x86_64-slc5-gcc43-opt/root/lib/libCore.so
#3 0x00002b22a18524cc in TUnixSystem::DispatchSignals(ESignals) () from /afs/cern.ch/sw/lcg/app/releases/ROOT/5.27.06/x86_64-slc5-gcc43-opt/root/lib/libCore.so
#4