Fitting Progress

 -  Templates

     -  Msig(x) : InvMass-signal from MC reconstructed (CrystalBall function) 

     -  MBkg(x) : InvMass-Bkg from Data (3rd order Polynomial) 

     -  Fbkg(x) : PsProper Decay length distribution (x) for Bkg candidates 

     -  Fsig(x)** :

       - R(x) : Prompt Jpsi Template from MC-truth reconstructed

       - CsiB(x) : Non-Prompt Jpsi Template from MC-truth 

 - Model for UnBinned two dimensional-fitting

     - Signal Part (Msig(m), Fsig(x))

     - Bkg Part (Mbkg(m), Fbkg(x))

     --> F(m,x) = f.Msig(m).Fsig(x) + (1-f).Mbkg(m).Fbkg(x)   <-- Simplistic representation of Fitting Function

     -> Results - NonPrompt Fraction (from x-dimension)

                - S/(S+B) (from mass-dimension)

​

%jsroot on

#ifndef __CINT__

#include "RooGlobalFunc.h"

#endif

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooDataHist.h"

#include "RooGaussian.h"

#include "RooConstVar.h"

#include "RooCategory.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "RooPlot.h"

#include "TFile.h"

using namespace RooFit ;

TFile *fMsig = TFile::Open("/eos/user/h/hsharma/SWAN_projects/Distributions/MC_signal/InvTrueMass_CBfit_pT_3_5GeV.root");

TFile *fMBkg = TFile::Open("MbkgFunc_3_5GeV.root");

TFile *fFB = TFile::Open("FB_3_5GeV.root");

​

​

RooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby 
                Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University
                All rights reserved, please read http://roofit.sourceforge.net/license.txt

## Making Msig(m) and Mbkg(m) PDFs

RooRealVar M("M","M",2.4,4.0);

​

double SigPars[5];

double BkgPars[4];

​

TF1 *MSigFunc = (TF1*)fMsig->Get("invMassSignal");

TF1 *MBkgFunc = (TF1*)fMBkg->Get("fitPol3");

​

MSigFunc->GetParameters(SigPars);

MBkgFunc->GetParameters(BkgPars);

​

cout << SigPars[0] <<  " : " <<  BkgPars[0] << endl;

cout << SigPars[1] <<  " : " <<  BkgPars[1] << endl;

cout << SigPars[2] <<  " : " <<  BkgPars[2] << endl;

cout << SigPars[3] <<  " : " <<  BkgPars[3] << endl;

cout << SigPars[4]  << endl;

​

​

​

RooRealVar s1("s1","s1",SigPars[0]);

RooRealVar s2("s2","s2",SigPars[1]);

RooRealVar s3("s3","s3",SigPars[2]);

RooRealVar s4("s4","s4",SigPars[3]);

RooRealVar s5("s5","s5",SigPars[4]);

​

RooRealVar b1("b1","b1",BkgPars[0]);

RooRealVar b2("b2","b2",BkgPars[1]);

RooRealVar b3("b3","b3",BkgPars[2]);

RooRealVar b4("b4","b4",BkgPars[3]);

​

​

​

​

​

​

RooTFnPdfBinding Msig("Msig","Msig",MSigFunc,RooArgList(M,s1,s2,s3,s4,s5));   // Signal : Blue

RooTFnPdfBinding Mbkg("Mbkg","Mbkg",MBkgFunc,RooArgList(M,b1,b2,b3,b4));     // Bkg : Reg

​

RooPlot *ff1 =M.frame((Title("InvMass PDF templates 3 < pT < 5")));

Msig.plotOn(ff1);

Mbkg.plotOn(ff1);

ff1->Draw();

c1->Draw();

​

​

​

3.08686 : 79395.4
1.76049 : -27163.3
0.0285696 : 416.651
0.373257 : -464178
44.3565
[#1] INFO:NumericIntegration -- RooRealIntegral::init(Msig_Int[M]) using numeric integrator RooIntegrator1D to calculate Int(M)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(Mbkg_Int[M]) using numeric integrator RooIntegrator1D to calculate Int(M)

xxxxxxxxxx

## R(x) PDF

​

double m1=-0.000398836;

double m2=-6.24013e-05;

double s1=0.00952525;

double s2=0.00174459;

​

double alp=0.00500007;

double lam=2.71545 ;

​

double fG1 = 0.272099;

double fG2 = 0.131904;

​

RooRealVar lxy("lxy","lxy",-0.6,0.6);

​

RooRealVar sigma1("sigma1","sigma1",s1);

RooRealVar mean1("mean1","mean1",m1);

RooGaussian gaus1("gaus1","gaus1",lxy,mean1,sigma1);

​

RooRealVar sigma2("sigma2","sigma2",s2);

RooRealVar mean2("mean2","mean2",m2);

RooGaussian gaus2("gaus2","gaus2",lxy,mean2,sigma2);

​

RooRealVar alpha("alpha", "alpha", alp);

RooRealVar lambda("lambda", "lambda",lam);

RooGenericPdf powerlaw("powerlaw", "powerlaw", "(TMath::Abs(lxy) > alpha ? TMath::Power(TMath::Abs(lxy),-lambda) : TMath::Power(alpha,-lambda))", RooArgSet(lxy, alpha, lambda));

​

RooRealVar *fracResolution1 = new RooRealVar("fracResolution1", "fracResolution1",fG1);

RooRealVar *fracResolution2 = new RooRealVar("fracResolution2", "fracResolution2",fG2);   

RooAddPdf *Rx = new RooAddPdf("Rx", "Rx", RooArgSet(gaus1, gaus2, powerlaw), RooArgList(*fracResolution1, *fracResolution2), kTRUE);

​

RooPlot * f5 = lxy.frame(Title("R(x) Template 3 < pT < 5"));

f5->SetMinimum(0.00001);

TCanvas cc;

​

Rx->plotOn(f5);

f5->Draw();

​

cc.Draw();

//cc.SetLogy();

​

[#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)

xxxxxxxxxx

## CsiB(x) : PsProperDecay length distribution for NonPrompt Jpsi  3.0 < pT < 5.0

TFile *fIn = TFile::Open("/eos/user/h/hsharma/SWAN_projects/Distributions/MC_signal/AnalysisHistograms_jpsi2ee_NP_PsProper.root");

TH2F * x_pt=(TH2F*)fIn->Get("jpsi2eeHistos")->FindObject("mcTruthJpsiNonPrompt_PureMCTruth_AfterSelection")->FindObject("kPseudoProperDecayTimeMC_pTMC");

​

int pT1 = 3;

int pT2 = 5;

​

​

x_pt->GetXaxis()->SetRange(x_pt->GetXaxis()->FindBin(3),x_pt->GetXaxis()->FindBin(pT2-1));

TCanvas c2;

c2.cd();

x_pt->ProjectionY()->Draw();

c2.Draw();

xxxxxxxxxx

## NonPrompt x-PDF template CsiB(x) * R(x) 

RooFFTConvPdf *FB = (RooFFTConvPdf*)fFB->Get("FB");

RooPlot *ff2 = .frame();

TCanvas c3;

FB->plotOn(ff2);

ff2->Draw();

c3.Draw();

[#1] INFO:Eval -- RooRealVar::setRange(lxy) new range named 'refrange_fft_FB' created with bounds [-0.6,0.6]
[#0] WARNING:Eval -- The FFT convolution 'FB' will run with 100 bins. A decent accuracy for difficult convolutions is typically only reached with n >= 1000. Suggest to increase the number of bins of the observable 'lxy'.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(CsiB_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 110 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 111 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 112 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 113 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 114 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 115 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 116 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 117 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 118 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 119 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 120 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 121 is out of range (0,109)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy|refrange_fft_FB]_Norm[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 110 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 111 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 112 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 113 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 114 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 115 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 116 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 117 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 118 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 119 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 120 is out of range (0,109)
[#0] ERROR:InputArguments -- RooAbsRealLValue::setBin(lxy) ERROR: bin index 121 is out of range (0,109)

[#1] INFO:Caching -- RooAbsCachedPdf::getCache(FB) creating new cache 0x7f693dd24980 with pdf CsiB_CONV_Rx_CACHE_Obs[lxy] for nset (lxy) with code 0

x

​

RooRealVar fB("fB","fB",0.1,0.0000001,0.999999);

RooAddPdf Fsig("Fsig","Fsig",RooArgSet(*Rx,*FB),RooArgList(fB));

lxy.setRange("ConRange",-0.6,0.6);

lxy.setBins(5);

lxy.setBin(0,"ConRange");

​

TCanvas c2;

​

RooPlot *ff0 = lxy.frame();

Fsig.plotOn(ff0);

ff0->Draw();

c2.Draw();

​

[#1] INFO:Eval -- RooRealVar::setRange(lxy) new range named 'ConRange' created with bounds [-0.6,0.6]
[#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)
[#1] INFO:Eval -- RooRealVar::setRange(lxy) new range named 'refrange_fft_FB' created with bounds [-0.6,0.6]
[#0] WARNING:Eval -- The FFT convolution 'FB' will run with 5 bins. A decent accuracy for difficult convolutions is typically only reached with n >= 1000. Suggest to increase the number of bins of the observable 'lxy'.
[#1] INFO:NumericIntegration -- RooRealIntegral::init(CsiB_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy)

[#1] INFO:Caching -- RooAbsCachedPdf::getCache(FB) creating new cache 0x7fa54d121a60 with pdf CsiB_CONV_Rx_CACHE_Obs[lxy] for nset (lxy) with code 0

​

RooGenericPdf *testexpC;

 RooFFTConvPdf *lxybkgC;

​

 RooGenericPdf *testexpL;

 RooFFTConvPdf *lxybkgL;

​

 RooGenericPdf *testexpR;

 RooFFTConvPdf *lxybkgR;

​

​

RooRealVar *fracS;

 RooRealVar *fracL;

 RooAddPdf *convpdf_lxybkg;

​

//RooRealVar lxy("lxy","lxy",-0.4,0.4);

​

 PsProperDL.setRange("Right", 0.01, 0.2);

 PsProperDL.setRange("Central", -0.005, 0.005);

 PsProperDL.setRange("Left", -0.2, -0.01);

​

RooRealVar muC("muC", "muC", 0.2, 0.002, 5.);

 testexpC = new RooGenericPdf("testexpC", "testexpC", "TMath::Exp(-1. * TMath::Abs(PsProperDL) / muC)", RooArgSet(PsProperDL, muC));

 lxybkgC = new RooFFTConvPdf("lxybkgC", "lxybkgC", PsProperDL, *testexpC, *model_resolution);

// lxybkgC->fitTo(*ds, Range("Central"), Timer(kTRUE));

​

 RooRealVar muR("muR", "muR", 0.02, 0.002, 2.);

 testexpR = new RooGenericPdf("testexpR", "testexpR", "(PsProperDL > 0 ? TMath::Exp(-1. * TMath::Abs(PsProperDL) / muR) : 0 )", RooArgSet(PsProperDL, muR));

 lxybkgR = new RooFFTConvPdf("lxybkgR", "lxybkgR", PsProperDL, *testexpR, *model_resolution);

// lxybkgR->fitTo(*ds, Range("Right"), Timer(kTRUE));

​

 RooRealVar muL("muL", "slope for leftside", 0.02, 0.002, 2.);

 testexpL = new RooGenericPdf("testexpL", "testexpL", "(PsProperDL < 0 ? TMath::Exp(-1. * TMath::Abs(PsProperDL) / muL) : 0 )", RooArgSet(PsProperDL, muL));

 lxybkgL = new RooFFTConvPdf("lxybkgL", "lxybkgL", PsProperDL, *testexpL, *model_resolution);

// lxybkgL->fitTo(*ds, Range("Left"), Timer(kTRUE));

​

 fracS = new RooRealVar("fracS","fracS", 0.5, 0, 1.);

 fracL = new RooRealVar("fracL","fracL", 0.5, 0, 1.);

​

xxxxxxxxxx

## 2D-Fitting model in (m,x) dimensions

RooProdPdf SignalPart("SignalPart","SignalPart",RooArgList(Msig,Fsig)) ;

RooProdPdf BkgPart("BkgPart","BkgPart",RooArgList(Mbkg,Fbkg)) ;

​

RooRealVar fsig("fsig","fsig",0.1,0,1);

​

RooAddPdf Total_Func("Total_Func","Total_Func",RooArgSet(SignalPart,BkgPart),RooArgList(fsig));

xxxxxxxxxx

# Next Step towards loading unbinned data and fitting the model (to be contd...)

 Fitting Progress - Templates     - Msig(x) : InvMass-signal from MC reconstructed (CrystalBall function)     - MBkg(x) : InvMass-Bkg from Data (3rd order Polynomial)     - Fbkg(x) : PsProper Decay length distribution (x) for Bkg candidates     - Fsig(x)** :       - R(x) : Prompt Jpsi Template from MC-truth reconstructed       - CsiB(x) : Non-Prompt Jpsi Template from MC-truth           - Model for UnBinned two dimensional-fitting     - Signal Part (Msig(m), Fsig(x))     - Bkg Part (Mbkg(m), Fbkg(x))         --> F(m,x) = f.Msig(m).Fsig(x) + (1-f).Mbkg(m).Fbkg(x)   <-- Simplistic representation of Fitting Function         -> Results - NonPrompt Fraction (from x-dimension)               - S/(S+B) (from mass-dimension)                               In [ ]: %jsroot on#ifndef __CINT__#include "RooGlobalFunc.h"#endif#include "RooRealVar.h"#include "RooDataSet.h"#include "RooDataHist.h"#include "RooGaussian.h"#include "RooConstVar.h"#include "RooCategory.h"#include "TCanvas.h"#include "TAxis.h"#include "RooPlot.h"#include "TFile.h"using namespace RooFit ;In [2]: TFile *fMsig = TFile::Open("/eos/user/h/hsharma/SWAN_projects/Distributions/MC_signal/InvTrueMass_CBfit_pT_3_5GeV.root");TFile *fMBkg = TFile::Open("MbkgFunc_3_5GeV.root");TFile *fFB = TFile::Open("FB_3_5GeV.root");​​ RooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University All rights reserved, please read http://roofit.sourceforge.net/license.txt Making Msig(m) and Mbkg(m) PDFs In [19]: RooRealVar M("M","M",2.4,4.0);​double SigPars[5];double BkgPars[4];​TF1 *MSigFunc = (TF1*)fMsig->Get("invMassSignal");TF1 *MBkgFunc = (TF1*)fMBkg->Get("fitPol3");​MSigFunc->GetParameters(SigPars);MBkgFunc->GetParameters(BkgPars);​cout << SigPars[0] << " : " << BkgPars[0] << endl;cout << SigPars[1] << " : " << BkgPars[1] << endl;cout << SigPars[2] << " : " << BkgPars[2] << endl;cout << SigPars[3] << " : " << BkgPars[3] << endl;cout << SigPars[4] << endl;​​​RooRealVar s1("s1","s1",SigPars[0]);RooRealVar s2("s2","s2",SigPars[1]);RooRealVar s3("s3","s3",SigPars[2]);RooRealVar s4("s4","s4",SigPars[3]);RooRealVar s5("s5","s5",SigPars[4]);​RooRealVar b1("b1","b1",BkgPars[0]);RooRealVar b2("b2","b2",BkgPars[1]);RooRealVar b3("b3","b3",BkgPars[2]);RooRealVar b4("b4","b4",BkgPars[3]);​​​​​​RooTFnPdfBinding Msig("Msig","Msig",MSigFunc,RooArgList(M,s1,s2,s3,s4,s5));RooTFnPdfBinding Mbkg("Mbkg","Mbkg",MBkgFunc,RooArgList(M,b1,b2,b3,b4));​RooPlot *ff1 =M.frame((Title("InvMass PDF templates")));Msig.plotOn(ff1);Mbkg.plotOn(ff1);ff1->Draw();c1->Draw(); 3.08686 : 79395.4 1.76049 : -27163.3 0.0285696 : 416.651 0.373257 : -464178 44.3565 [#1] INFO:NumericIntegration -- RooRealIntegral::init(Msig_Int[M]) using numeric integrator RooIntegrator1D to calculate Int(M) [#1] INFO:NumericIntegration -- RooRealIntegral::init(Mbkg_Int[M]) using numeric integrator RooIntegrator1D to calculate Int(M) Fsig(x) PDF In [3]: ​double m1=-0.000398836;double m2=-6.24013e-05;double s1=0.00952525;double s2=0.00174459;​double alp=0.00500007;double lam=2.71545 ;​double fG1 = 0.272099;double fG2 = 0.131904;​RooRealVar lxy("lxy","lxy",-0.6,0.6);​RooRealVar sigma1("sigma1","sigma1",s1);RooRealVar mean1("mean1","mean1",m1);RooGaussian gaus1("gaus1","gaus1",lxy,mean1,sigma1);​RooRealVar sigma2("sigma2","sigma2",s2);RooRealVar mean2("mean2","mean2",m2);RooGaussian gaus2("gaus2","gaus2",lxy,mean2,sigma2);​RooRealVar alpha("alpha", "alpha", alp);RooRealVar lambda("lambda", "lambda",lam);RooGenericPdf powerlaw("powerlaw", "powerlaw", "(TMath::Abs(lxy) > alpha ? TMath::Power(TMath::Abs(lxy),-lambda) : TMath::Power(alpha,-lambda))", RooArgSet(lxy, alpha, lambda));​RooRealVar *fracResolution1 = new RooRealVar("fracResolution1", "fracResolution1",fG1);RooRealVar *fracResolution2 = new RooRealVar("fracResolution2", "fracResolution2",fG2); RooAddPdf *Rx = new RooAddPdf("Rx", "Rx", RooArgSet(gaus1, gaus2, powerlaw), RooArgList(*fracResolution1, *fracResolution2), kTRUE);​RooPlot * f5 = lxy.frame();f5->SetMinimum(0.00001);TCanvas cc;​Rx->plotOn(f5);f5->Draw();​cc.Draw();//cc.SetLogy();​ ROOT canvas-0.6-0.4-0.200.20.40.6lxy00.10.20.30.40.50.60.70.80.9Projection of RxA RooPlot of "lxy"Toggle tool buttons [#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy) In [10]:x ​RooFFTConvPdf *FB = (RooFFTConvPdf*)fFB->Get("FB");​RooRealVar fB("fB","fB",0.1,0.0000001,0.999999);​RooAddPdf Fsig("Fsig","Fsig",RooArgSet(*Rx,*FB),RooArgList(fB));​In [5]: lxy.setRange("ConRange",-0.6,0.6);lxy.setBins(5);lxy.setBin(0,"ConRange");​TCanvas c2;​RooPlot *ff0 = lxy.frame();Fsig.plotOn(ff0);ff0->Draw();c2.Draw();​[#1] INFO:Eval -- RooRealVar::setRange(lxy) new range named 'ConRange' created with bounds [-0.6,0.6] [#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy) [#1] INFO:Eval -- RooRealVar::setRange(lxy) new range named 'refrange_fft_FB' created with bounds [-0.6,0.6] [#0] WARNING:Eval -- The FFT convolution 'FB' will run with 5 bins. A decent accuracy for difficult convolutions is typically only reached with n >= 1000. Suggest to increase the number of bins of the observable 'lxy'. [#1] INFO:NumericIntegration -- RooRealIntegral::init(CsiB_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy) [#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[lxy]) using numeric integrator RooIntegrator1D to calculate Int(lxy) ROOT canvas-0.6-0.4-0.200.20.40.6lxy10-310-210-1100Projection of FsigA RooPlot of "lxy"Toggle tool buttons [#1] INFO:Caching -- RooAbsCachedPdf::getCache(FB) creating new cache 0x7fa54d121a60 with pdf CsiB_CONV_Rx_CACHE_Obs[lxy] for nset (lxy) with code 0 In [ ]: ​In [ ]: RooGenericPdf *testexpC; RooFFTConvPdf *lxybkgC;​ RooGenericPdf *testexpL; RooFFTConvPdf *lxybkgL;​ RooGenericPdf *testexpR; RooFFTConvPdf *lxybkgR;​​RooRealVar *fracS; RooRealVar *fracL; RooAddPdf *convpdf_lxybkg;​//RooRealVar lxy("lxy","lxy",-0.4,0.4);​ PsProperDL.setRange("Right", 0.01, 0.2); PsProperDL.setRange("Central", -0.005, 0.005); PsProperDL.setRange("Left", -0.2, -0.01);​RooRealVar muC("muC", "muC", 0.2, 0.002, 5.); testexpC = new RooGenericPdf("testexpC", "testexpC", "TMath::Exp(-1. * TMath::Abs(PsProperDL) / muC)", RooArgSet(PsProperDL, muC)); lxybkgC = new RooFFTConvPdf("lxybkgC", "lxybkgC", PsProperDL, *testexpC, *model_resolution);// lxybkgC->fitTo(*ds, Range("Central"), Timer(kTRUE));​ RooRealVar muR("muR", "muR", 0.02, 0.002, 2.); testexpR = new RooGenericPdf("testexpR", "testexpR", "(PsProperDL > 0 ? TMath::Exp(-1. * TMath::Abs(PsProperDL) / muR) : 0 )", RooArgSet(PsProperDL, muR)); lxybkgR = new RooFFTConvPdf("lxybkgR", "lxybkgR", PsProperDL, *testexpR, *model_resolution);// lxybkgR->fitTo(*ds, Range("Right"), Timer(kTRUE));​ RooRealVar muL("muL", "slope for leftside", 0.02, 0.002, 2.); testexpL = new RooGenericPdf("testexpL", "testexpL", "(PsProperDL < 0 ? TMath::Exp(-1. * TMath::Abs(PsProperDL) / muL) : 0 )", RooArgSet(PsProperDL, muL)); lxybkgL = new RooFFTConvPdf("lxybkgL", "lxybkgL", PsProperDL, *testexpL, *model_resolution);// lxybkgL->fitTo(*ds, Range("Left"), Timer(kTRUE));​ fracS = new RooRealVar("fracS","fracS", 0.5, 0, 1.); fracL = new RooRealVar("fracL","fracL", 0.5, 0, 1.);​2D-Fitting model in (m,x) dimensions In [ ]: RooProdPdf SignalPart("SignalPart","SignalPart",RooArgList(Msig,Fsig)) ;RooProdPdf BkgPart("BkgPart","BkgPart",RooArgList(Mbkg,Fbkg)) ;​RooRealVar fsig("fsig","fsig",0.1,0,1);​RooAddPdf Total_Func("Total_Func","Total_Func",RooArgSet(SignalPart,BkgPart),RooArgList(fsig));