mmgData.SetName('mmgData')
        # w.Import(mmgData)

        ## Define the translated mass
        mmgMassFunc = w.factory('''FormulaVar::mmgMassFunc(
            "((mmgMass - mmgMode) / mmgScale + mmgMode)  - fPho*phoScale/100.",
            {mmgMass, mmgMode[91.2,70,110], mmgScale[1.,0.1,5], fPho, phoScale}
            )''')
        mmgMode = w.var('mmgMode')
        mmgScale = w.var('mmgScale')

        ## Get the nominal model for translation, use rho=2
        model = RooKeysPdf('model', 'model', mmgMass, mmgData, RooKeysPdf.NoMirror, 2)

        ## Sample nominal model to binned data, no Poisson fluctuation
        modelHist = model.createHistogram('mmgMass', 10000)
        modelDH = RooDataHist('modelDH', 'modelDH', RooArgList(mmgMass), modelHist)

        ## Use the sampled model to build a transformed model, interpolation order=2
        tmodel = RooHistPdf('tmodel', 'tmodel', RooArgList(mmgMassFunc),
                            RooArgList(mmgMass), modelDH, 2)


        ## Find the mode of the tmodel.
        mmgMode.setVal(modelHist.GetBinCenter(modelHist.GetMaximumBin()))
        mmgMode.setConstant(True)
        fPho.setConstant(True)

        ## Fit the photon scale
        phoScale.setConstant(False)
        mmgScale.setConstant(True)
Ejemplo n.º 2
0
data1 = g.generate(ArgSet(x), 10000)

w.var('gm').setVal(-3.)
w.var('gs').setVal(2)
data2 = g.generate(ArgSet(x), 10000)

w.var('gm').setVal(0.)
w.var('gs').setVal(1.)

## 3. Use data 1 to build Keys PDF k1(x|s1,m1)
k = KeysPdf('k', 'k', x, data1, KeysPdf.NoMirror, 2)

## 4. Use k to build HistPDF with _transformed_ x variable t = (x-m)/s
##    h1((x-hm)/hs|m1,s1)
t = w.factory('FormulaVar::t("(x - hm) / hs", {x, hm[0,-5,5], hs[1,0.1,3]})')
hist = k.createHistogram('x', 10000)
dh = DataHist('dh', 'dh', ArgList(x), hist)
h = HistPdf('h', 'h', ArgList(t), ArgList(x), dh, 2)

## 5. Plot results
f1 = x.frame()
f1.SetTitle('Training Data')
data1.plotOn(f1)
g.plotOn(f1)
k.plotOn(f1, LineColor(kRed), LineStyle(kDashed))
h.plotOn(f1, LineColor(kBlack), LineStyle(kDashed))
canvases.next('Training_Data')
f1.Draw()

## 6. Fit h1 to data2 to simulate a "measurement" of m2, s2.
h.fitTo(data2)