My apologies for the list spam on this issue…
I am trying to use RooFit to code up the following pdf:
Clearly, for all i,j combinations, a new subdataset/pdf is required.
In previous HEP analyses, I coped with subdatasets/pdfs by using RooSimultaneous and manually defining subdatasets/pdfs:
RooCategory Category("Category","Category"); Category.defineType("PiPlusDataMagDown"); Category.defineType("PiMinusDataMagDown"); map<std::string,RooDataHist*> myMap; myMap.insert( pair<std::string,RooDataHist*>("PiPlusDataMagDown",Data_Pi_Offline_Positive_Hist_MagDown) ); myMap.insert( pair<std::string,RooDataHist*>("PiMinusDataMagDown",Data_Pi_Offline_Negative_Hist_MagDown) ); RooDataHist combData("combData","combined data", RooArgList(rd_Dplus_Dplus_DTF_M), Category, myMap); RooSimultaneous simPdf("simPdf","simultaneous pdf",Category) ; simPdf.addPdf(*sumPDF_Pi_MagDown,"PiPlusDataMagDown") ; // Pi plus data simPdf.addPdf(*sumPDF_PiM_MagDown,"PiMinusDataMagDown") ; // Pi minus data
With just a few subdatasets, this is clearly fine. However, in this case i and j run from 1-41 and so defining so many subdatasets/pdfs is clearly not viable. I kind of need a vector solution. I tried using RooSimPdfBuilder, but having chatted with Lorenzo (and others), it doesn’t seem to quite do what I need.
Could anyone offer any advice on how they would code up this pdf or point me to a similar tutorial/example?
Thank in advance,