Dear all,
I am continuing to use the RooSimPdfBuilder class and the tutorial here:
roofit.sourceforge.net/docs/clas … t2.cc.html
…and have everything working happily
When I build the pdf:
// Build the PDF
RooSimultaneous* simSum = mgr.buildPdf(*config,&data) ;
simSum->Print("v") ;
simSum->getParameters(&data)->Print("v") ;
and print, I see output like the following:
Plot label is "simPdf"
--- RooAbsPdf ---
Cached value = 0
1) 0x46e4aa0 RooRealVar:: a_a_param_a1 = 0.5 L(0 - 1) "a_a_param (a1)"
2) 0x46e6f20 RooRealVar:: a_a_param_a2 = 0.5 L(0 - 1) "a_a_param (a2)"
3) 0x46e3d80 RooRealVar:: a_d_param_a1 = 0.5 L(0 - 1) "a_d_param (a1)"
4) 0x46fb020 RooRealVar:: a_d_param_a2 = 0.5 L(0 - 1) "a_d_param (a2)"
5) 0x7fff9465a020 RooRealVar:: gamma = 0.5 L(0 - 1) "gamma"
6) 0x46ef570 RooRealVar:: h_a_param_h1 = 0.5 L(0 - 1) "h_a_param (h1)"
7) 0x4704bc0 RooRealVar:: h_a_param_h2 = 0.5 L(0 - 1) "h_a_param (h2)"
8) 0x46fc120 RooRealVar:: h_d_param_h1 = 0.5 L(0 - 1) "h_d_param (h1)"
9) 0x47176d0 RooRealVar:: h_d_param_h2 = 0.5 L(0 - 1) "h_d_param (h2)"
But, I’m failing to work out how to access/edit the parameters that have been automatically generated e.g. a_a_param_a1
e.g:
a_a_param_a1.getVal();
doesn’t work…
Could anyone help me with the correct syntax?
Follow up question: I’d like to get at these values so I can e.g.
- constrain them in the fit ( sum over a_a_param_a(i) = 1) - is this possible?
- print fitted parameters etc etc
Thanks!
S