MultiVariateGaussian

nsahoo · January 18, 2017, 11:52am

Hi,

I’m trying to generate a 10 efficiency functions where the parameters are varied in multi-Gaussian phase space (the parameters are varied by their errors). The snipped code as well as the .txt file containing all the input parameters are attached below.

void MultGausT(int iBin, const char outfile[] = “accXrecoEff”){

//------cosThetaL-----
int dimL = 7;
RooArgList yVecL;
RooArgList muVecL;
int i,j;
RooRealVar* yL;
RooRealVar* mu_yL;
for (i = 0; i < dimL; i++) {
char* name = Form(“yL%d”, i);
double mmean;
mmean = readParam(iBin,“accXrecoEff2_ord0_L”,i);
double emean;
emean = readParam(iBin,“accXrecoEff2Err_ord0_L”,i);
yL = new RooRealVar(name, name, mmean-emean, mmean+emean);
yVecL.add(yL);
char mu_name = Form(“mu_yL%d”,i);
mu_yL = new RooRealVar(mu_name, mu_name, readParam(iBin,“accXrecoEff2_ord0_L”,i));
muVecL.add(*mu_yL);
}

TMatrixDSym covL(dimL);
for (i = 0; i < dimL; i++) {
for (j = 0; j < dimL; j++) {
double errParL = readParam(iBin,TString::Format(“covMat_L%d”,i+1), j);
covL(i,j) = errParL;
}
}

RooMultiVarGaussian mvgL(“mvgL”, “mvgL”, yVecL, muVecL, covL);
RooDataSet* meeeL = mvgL.generate(yVecL, 10);
meeeL->RooDataSet::write(TString::Format(“eff_parL_%d.txt”,iBin));

//return;
}

Actually I’m running over two bins, the .txt file attached is for bin4. It’s reading correctly all the input parameters but is getting stuck at the generation step(see the last 3 lines of the above code).

The printout is also attached here:

[nsahoo@lxplus005 plugins]$ ./fit MultGausT x

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

INFO : ScaleFactor for input data is 1.000
INFO : Plots will be stored to ./plots
INFO : Datacards will be stored to .
INFO : Workspaces will be stored to .
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[0]=6.986006e-05
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[0]=8.186301e-07
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[0]=6.986006e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[1]=8.831249e-06
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[1]=3.192351e-06
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[1]=8.831249e-06
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[2]=4.726557e-05
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[2]=1.106756e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[2]=4.726557e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[3]=-2.107568e-07
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[3]=1.412946e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[3]=-2.107568e-07
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[4]=-4.997589e-05
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[4]=3.287010e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[4]=-4.997589e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[5]=-2.357662e-06
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[5]=1.352825e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[5]=-2.357662e-06
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[6]=-2.020628e-05
INFO: readParam, matched accXrecoEff2Err_ord0_L!
INFO: readParam, get accXrecoEff2Err_ord0_L[6]=2.546049e-05
INFO: readParam, matched accXrecoEff2_ord0_L!
INFO: readParam, get accXrecoEff2_ord0_L[6]=-2.020628e-05
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[0]=7.822817e-04
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[1]=7.987476e-04
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[2]=-3.377330e-04
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[3]=-5.313775e-04
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[4]=-4.882346e-04
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[5]=3.013436e-04
INFO: readParam, matched covMat_L1!
INFO: readParam, get covMat_L1[6]=1.222489e-04
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[0]=7.987476e-04
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[1]=1.657146e-03
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[2]=5.730430e-04
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[3]=-5.745665e-04
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[4]=-5.542780e-05
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[5]=-3.402064e-04
INFO: readParam, matched covMat_L2!
INFO: readParam, get covMat_L2[6]=-2.117853e-04
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[0]=-3.377330e-04
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[1]=5.730430e-04
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[2]=2.987204e-03
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[3]=2.456434e-03
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[4]=3.179746e-06
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[5]=7.535262e-05
INFO: readParam, matched covMat_L3!
INFO: readParam, get covMat_L3[6]=-1.413877e-03
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[0]=-5.313775e-04
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[1]=-5.745665e-04
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[2]=2.456434e-03
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[3]=4.486167e-03
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[4]=-3.126766e-04
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[5]=8.195578e-05
INFO: readParam, matched covMat_L4!
INFO: readParam, get covMat_L4[6]=-2.407801e-04
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[0]=-4.882346e-04
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[1]=-5.542780e-05
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[2]=3.179746e-06
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[3]=-3.126766e-04
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[4]=5.428042e-03
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[5]=1.244709e-03
INFO: readParam, matched covMat_L5!
INFO: readParam, get covMat_L5[6]=-1.490179e-03
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[0]=3.013436e-04
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[1]=-3.402064e-04
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[2]=7.535262e-05
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[3]=8.195578e-05
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[4]=1.244709e-03
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[5]=1.242629e-02
INFO: readParam, matched covMat_L6!
INFO: readParam, get covMat_L6[6]=2.499442e-03
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[0]=1.222489e-04
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[1]=-2.117853e-04
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[2]=-1.413877e-03
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[3]=-2.407801e-04
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[4]=-1.490179e-03
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[5]=2.499442e-03
INFO: readParam, matched covMat_L7!
INFO: readParam, get covMat_L7[6]=1.765099e-02

Please let me know how to sort this out.

cheers, Niladri
fitParameters4.txt (7.5 KB)

moneta · January 18, 2017, 3:02pm

HI,

We cannot sort out your problem without a minimal runnable program, the text output is not enough. Please attach it

Best Regards

Lorenzo

nsahoo · January 18, 2017, 7:56pm

Hi Lorenzo,

Here are the files attached. Please do “make” for compilation and then type

./fit MultivariateGaussianTest 4

to run over bin#4 and the .txt file for bin#4 is “fitParameters4.txt”. I tried to attach the Makefile, but it showed me “wrong input file format”, so I put it as .txt file. Please remove .txt extension before doing compilation.

cheers - Niladri
Makefile.txt (774 Bytes)
SingleBuToKstarMuMuFitters.cc (24.1 KB)
fitParameters4.txt (6.07 KB)