03/14/2007

Update the single tags fit

dhad fit regular1 -t d —qsub Log

dhad fit regular1 -t s —qsub Log

dhad fit regular1 -t g —qsub Log

Update the double tags fit

dhad fit regular1 -t d —tag d —qsub Log

dhad fit regular1 -t s —tag d —qsub Log

dhad fit regular1 -t g —tag d —qsub Log

BF fitter

dhad bf

Result On HyperNews

03/15/2007

Sort out the BF procedure

need to parse the statonly

Change the name of regular1 to regular and make the page update

Resubmit the double tag fit for data and signal MC:

dhad fit regular1 -t d —tag d —qsub

Log

dhad fit regular1 -t s —tag d —qsub Log

R= 2.5 fit

/nfs/cor/user/ryd/RooFit/01-09-02/RooFitModels/RooDEnergyImp.cc

dhad fit BES2006lineshape Gamma 0.0286 M 3.7718 R 2.5 -t d —qsub Log

03/16/2007

Sort out the BF procedure

need to parse the statonly

done.

03/17/2007

Produce the Sys errors file : newfit-data_systerrors

03/19/2007

Update the BF fit

Back up the last result as: "bf_stat_sys_0"

dhad bf fit1 fit2

Copy the new one as bf_stat_sys_1

Parse the new result: bf_stat_sys

dhad bf parse

Need to transform to table ... dhad table bf ...

Submit the other fits

Double tag fits need backgrounds:

dhad fit regular1 -t d —tag d —qsub Log

dhad fit regular1 -t s —tag d —qsub Log

dhad fit regular1 -t g —tag d —qsub Log

03/20/2007

Create the data result table

dhad table data_results

Redo the lineshape para fit

backup the old fit and fig to fitparams0

dhad fit floatmass fitparams diag -t s —tag d —qsub Log

dhad fit Gamma 0.0225 -t d —tag d —qsub Log

dhad fit Gamma 0.0325 -t d —tag d —qsub Log

dhad fit P 0.5 -t d —tag d —qsub Log

Do the BF fit again

Compare the Signal eff change for the double tags

dhad table compare version 6.09 double_signal_eff 'Signal yield'

03/21/2007

Check the lineshape fits

Mode 0:

   1  N            8.90090e+03   9.44855e+01   2.24301e-04  -7.71702e-01
   2  Nbkgd1       8.51209e+00   5.20022e+00   3.47375e-04  -7.52826e+00
   3  Nbkgd2       5.41287e-01   2.18171e+00   5.35051e-04   2.23562e+01
   4  Nbkgd3       4.17146e-04   1.17976e+00   4.86015e-04** at limit **
   5  NbkgdFlat    1.36554e-05   9.13395e-01   4.28352e-04** at limit **
   6  f1a          1.76019e-01   3.61810e-02   9.94313e-03  -2.13402e-01
   7  f1b          4.63210e-03   1.37274e-03   6.56477e-03  -4.40806e-01
   8  md           1.86461e+00   1.47423e-05   1.26713e-03  -7.16205e+00
   9  s1a          2.32596e+00   1.05684e-01   1.17948e-02  -2.16632e-01
  10  s1b          3.88259e+00   1.86121e+00   2.68001e-01   2.71088e-02
  11  sigmacommon1   4.06924e-03   1.48274e-04   5.81883e-03   2.98148e-01
  RooRealVar N("N","N",4000.0 , 0.0, 2000000.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m0,0 --qsub

Log

Mode 3:

   1  N            4.56047e+03   6.77933e+01   1.58447e-04  -1.33104e+00
   2  Nbkgd1       5.57433e+00   4.15637e+00   2.34727e-04  -4.42715e+00
   3  Nbkgd2       1.01957e+01   4.30013e+00   2.01120e-04  -6.24872e-02
   4  Nbkgd3       1.56948e+01   5.17565e+00   2.05646e-04   4.72069e-01
   5  NbkgdFlat    1.13796e-04   5.05703e+00   6.82654e-04** at limit **
   6  f1a          9.60734e-02   3.84683e-02   9.91784e-03  -9.98796e-03
   7  f1b          9.34468e-03   5.40078e-03   7.13956e-03   2.35385e-02
   8  md           1.86462e+00   2.08049e-05   1.26056e-03   9.20743e-02
   9  s1a          2.37371e+00   3.40921e-01   1.71874e-02  -4.75843e-03
  10  s1b          2.81025e+00   5.97724e-01   5.93599e-02  -2.08202e-02
  11  sigmacommon1   4.60213e-03   1.72918e-04   5.63813e-03  -1.67561e-03
  RooRealVar N("N","N",4000.0 , 0.0, 2000000.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m3,3 --qsub

Log

Mode 200:

   1  N            6.06787e+03   7.79536e+01   1.84131e-04  -7.38297e+00
   2  Nbkgd1       3.77657e-05   2.15891e+00   5.37943e-04** at limit **
   3  Nbkgd2       6.73142e+00   3.93965e+00   2.66803e-04   6.32671e+00
   4  Nbkgd3       2.99689e-06   8.40493e-01   3.37022e-04** at limit **
   5  NbkgdFlat    2.20972e-07   6.71032e-01   3.01107e-04** at limit **
   6  f1a          1.36354e-01   3.70993e-02   1.03410e-02  -7.71293e-02
   7  f1b          7.24351e-03   1.63045e-03   5.89933e-03   1.25620e-01
   8  md           1.86939e+00   1.73596e-05   9.38512e-04   7.73567e-01
   9  s1a          2.22630e+00   1.37866e-01   1.38090e-02  -2.32379e-02
  10  s1b          4.00000e+00   1.63323e-01   1.34123e-01** at limit **
  11  sigmacommon1   4.36350e-03   1.57928e-04   5.59139e-03  -9.17346e-02
  RooRealVar N("N","N",4000.0 , 0.0, 2000000.0);
  RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m200,200 --qsub

Log

Mode 201:

   1  N            1.60314e+03   4.15276e+01   9.18356e-05  -7.29652e-01
   2  Nbkgd1       2.39481e+01   8.21589e+00   1.36085e-04  -3.91226e+00
   3  Nbkgd2       1.28582e+01   7.11854e+00   1.40748e-04  -1.15451e-01
   4  Nbkgd3       3.91592e-04   4.25938e+00   3.76487e-04** at limit **
   5  NbkgdFlat    1.30136e+01   5.44388e+00   1.30150e-04  -3.36011e+00
   6  f1a          1.86121e-01   6.10091e-02   1.00793e-02  -1.72013e-02
   7  f1b          5.12582e-02   1.70833e-02   6.53120e-03  -1.90840e-01
   8  md           1.86943e+00   3.82175e-05   1.02525e-03   1.29232e-01
   9  s1a          2.47464e+00   3.36083e-01   1.31777e-02  -9.04256e-02
  10  s1b          2.91734e+00   5.51285e-01   2.99430e-02  -1.79573e-02
  11  sigmacommon1   6.21256e-03   4.22739e-04   8.34786e-03  -2.35217e-02
  RooRealVar N("N","N",4000.0 , 0.0, 2000000.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m201,201 --qsub

Log

Mode 202:

   1  N            4.11397e+03   1.25228e+06  -0.00000e+00   8.57089e-01
   2  Nbkgd1       3.47122e+00   1.75271e+05   0.00000e+00   5.39651e-01
   3  Nbkgd2       1.60248e+00   6.72875e+03   0.00000e+00   3.42025e+00
   4  Nbkgd3       1.13646e-08   1.35070e+06   0.00000e+00** at limit **
   5  NbkgdFlat    4.32597e-09   1.35842e+06  -0.00000e+00** at limit **
   6  f1a          4.62954e-01   2.59544e-01   0.00000e+00  -7.71131e-02
   7  f1b          8.59981e-03   2.40820e-01   0.00000e+00  -9.36453e-02
   8  md           1.86938e+00   7.79964e-03  -0.00000e+00   5.85993e-02
   9  s1a          2.28358e+00   1.97267e+00   0.00000e+00  -9.25623e-02
  10  s1b          3.99993e+00   1.85956e+00  -0.00000e+00   9.45448e-03
  11  sigmacommon1   3.01017e-03   5.96699e-04   0.00000e+00  -1.90118e-01
  RooRealVar N("N","N",4000.0 , 0.0, 2000000.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m202,202 --qsub

Log

Mode 203:

   1  N            1.17033e+03   3.42642e+01   7.65025e-05   1.42352e-01
   2  Nbkgd1       7.37854e-10   1.11899e+00   1.61634e-04** at limit **
   3  Nbkgd2       1.96959e-07   1.53285e+00   1.89280e-04** at limit **
   4  Nbkgd3       4.63112e-07   2.11239e+00   2.22310e-04** at limit **
   5  NbkgdFlat    8.67540e+00   3.51797e+00   8.93693e-05   9.52155e-01
   6  f1a          5.00000e-01   3.68429e-01   1.56056e-01** at limit **
   7  f1b          9.23497e-02   3.31511e-02   6.46907e-03  -3.13694e-03
   8  md           1.86943e+00   4.82263e-05   1.06560e-03   3.75127e-02
   9  s1a          1.62286e+00   1.90871e+00   1.09960e-02  -8.69856e-03
  10  s1b          3.17266e+00   5.41874e-01   1.78856e-02   7.14586e-03
  11  sigmacommon1   6.05537e-03   3.80065e-03   6.98973e-03  -1.73259e-02

We used the old value as initial value; increased f1a's upper limit from 0.5 to 1.0; set three background components to be zeros.

  RooRealVar N("N","N",1158.0 , 0.0, 2000000.0);
  s1a=new RooRealVar("s1a","s1a",2.27,1.5,4.0);
  f1a=new RooRealVar("f1a","f1a",0.27,0.0,1.0);
  s1b=new RooRealVar("s1b","s1b",3.66,1.5,4.0);
  f1b=new RooRealVar("f1b","f1b",0.029,0.0,0.4);
  RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
  RooRealVar Nbkgd2("Nbkgd2","Nbkgd2",0.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m203,203 --qsub

Log

Mode 205:

   1  N            4.46317e+03   6.68785e+01   1.56756e-04  -1.86592e-02
   2  Nbkgd1       5.04759e-06   1.67406e+00   3.96918e-04** at limit **
   3  Nbkgd2       2.32990e+00   2.44113e+00   2.38837e-04  -1.83276e-01
   4  Nbkgd3       8.50146e+00   3.75089e+00   2.01037e-04   1.26310e-01
   5  NbkgdFlat    2.57769e-10   8.83716e-01   2.94033e-04** at limit **
   6  f1a          1.20989e-01   6.98253e-02   1.05665e-02   3.51617e-03
   7  f1b          8.55749e-03   4.30960e-03   6.58077e-03  -1.25395e-02
   8  md           1.86943e+00   2.04441e-05   9.40261e-04  -8.95673e-02
   9  s1a          2.22056e+00   3.53702e-01   1.49025e-02  -4.58204e-03
  10  s1b          3.12904e+00   5.03372e-01   5.87580e-02  -1.20996e-03
  11  sigmacommon1   4.56459e-03   2.59030e-04   5.59804e-03  -4.49168e-03
  RooRealVar N("N","N",4000.0 , 0.0, 2000000.0);
  RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m205,205 --qsub

Log

Submit all the fits

dhad yield widede -t s —qsub Log

dhad yield widede -t d —qsub Log

dhad yield noFSR -t s —qsub Log

dhad yield noFSR -t s —tag d —qsub Log

dhad fit Gamma 0.0325 -t d —qsub Log

dhad fit Gamma 0.0225 -t d —qsub Log

dhad fit P 0.5 -t d —qsub Log

dhad fit P 0.5 -t s —qsub Log

dhad fit widede -t d —qsub Log

dhad fit widede -t s —qsub Log

dhad fit noFSR -t s —qsub Log

dhad fit noFSR -t s —tag d Log

Update the PRD tables

dhad update prd table

done.

03/22/2007

Check the line shape fits

Mode 0 OK

Mode 3 OK

Mode 200:

   1  N            6.06844e+03   7.79586e+01   1.84127e-04   2.77779e+00
   2  Nbkgd2       6.72038e+00   3.95819e+00   2.66658e-04   1.93575e+00
   3  f1a          1.36602e-01   4.22485e-02   1.03768e-02   2.32943e-02
   4  f1b          7.24691e-03   1.69866e-03   5.90514e-03   2.79040e-02
   5  md           1.86939e+00   1.73598e-05   9.38577e-04   1.59535e-01
   6  s1a          2.22542e+00   1.54279e-01   1.37751e-02   1.16161e-02
   7  s1b          4.00000e+00   1.62682e-01   1.34050e-01** at limit **
   8  sigmacommon1   4.36285e-03   1.71832e-04   5.59220e-03   5.13077e-02

Fix 200:

  s1b=new RooRealVar("s1b","s1b",3.66,1.5,8.0);
  RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m200,200 --qsub

Log

Mode 201 OK

Mode 202

   1  N            4.11744e+03   6.38574e+01   2.12309e-06   5.86916e-01
   2  Nbkgd1       3.84957e-03   1.82276e+00   8.68599e-06** at limit **
   3  Nbkgd2       1.59112e+00   1.61668e+00  -7.67877e-06  -9.73796e-01
   4  f1a          1.66337e-01   5.40961e-02   1.16943e-03  -3.53780e-02
   5  f1b          4.91775e-03   1.70534e-03   5.30546e-04   8.99437e-02
   6  md           1.86938e+00   2.02252e-05  -2.39446e-05  -5.24236e-02
   7  s1a          2.25266e+00   1.43034e-01  -4.88236e-04  -4.03447e-02
   8  s1b          3.99996e+00   2.69662e-01  -7.46819e-03   1.41601e-02
   9  sigmacommon1   4.10586e-03   2.19420e-04  -7.59094e-04  -8.37288e-02

Fix 202

  RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m202,202 --qsub

Log

   1  N            4.11741e+03   6.41746e+01   1.50402e-04  -3.87980e-03
   2  Nbkgd2       1.58996e+00   1.59084e+00   1.88864e-04   2.05693e+00
   3  f1a          1.65930e-01   6.26740e-02   1.01643e-02   3.05184e-02
   4  f1b          4.90664e-03   1.82894e-03   6.01316e-03   5.30215e-02
   5  md           1.86938e+00   2.10069e-05   9.27944e-04   1.10872e-01
   6  s1a          2.25374e+00   1.67505e-01   6.58399e-03   5.26580e-02
   7  s1b          4.00000e+00   3.30024e-01   1.57651e-01** at limit **
   8  sigmacommon1   4.10744e-03   2.44340e-04   5.64715e-03   4.67340e-02

Try again:

  s1a=new RooRealVar("s1a","s1a",2.0,1.5,8.0);
  RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
  RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);
  RooRealVar NbkgdFlat("NbkgdFlat","NbkgdFlat",0.0);

dhad fit floatmass fitparams diag -t s --tag d -m202,202 --qsub

Log

....

Mode 203 OK

Mode 205 OK

Update Fitting pages

dhad yield widede noFSR -t s —qsub Log

dhad fit widede noFSR -t s —qsub Log

dhad fit Gamma float -t d —qsub Log

No FSR Fits

BES line shape fits

dhad fit BES2006lineshape Gamma 0.0286 M 3.7718 -t d —qsub

Log

P= 0.5 for single tag

Resubmit the p=0.5 for single tag

dhad fit P 0.5 -t d —qsub Log

dhad fit P 0.5 -t s —qsub Log

03/23/2007

BES line shape table

dhad table compare yields BES2006lineshape/Gamma/0.0286/M/3.7718/R/2.5 -t d -v1

Fix problem in the mode 204

   1  N1           1.70528e+04   2.39573e+02  -0.00000e+00   9.30990e+00
   2  N2           1.71831e+04   2.34772e+02  -0.00000e+00   1.00764e+01
   3  Nbkgd1       3.24254e+04   2.73950e+02  -0.00000e+00   1.02751e+01
   4  Nbkgd2       3.18072e+04   2.69788e+02  -0.00000e+00   1.05710e+01
   5  md           1.87073e+00   5.28483e-05  -0.00000e+00   2.16850e+00
   6  p            1.50000e+00   1.11982e-04  -0.00000e+00** at limit **
   7  sigmap1      4.00000e-02   3.79182e-07  -0.00000e+00** at limit **
   8  xi          -2.86434e+01   6.16519e-01  -0.00000e+00   2.79986e-01

Change the initial value:

   RooRealVar N1("N1","N1",8163.0,0.0,2000000.0);
   RooRealVar N2("N2","N2",8168.0,0.0,2000000.0);
   p = new RooRealVar("p","p",0.5, 0.1, 2.0);
   sigmap1 = new RooRealVar("sigmap1","sigmap1",0.0045,0.002,0.080);

dhad fit BES2006lineshape Gamma 0.0286 M 3.7718 -t d -m204 --qsub

Log

Fit BES mass and width

dhad fit Gamma 0.0286 M 3.7718 -t d —qsub Log

Make tables done.

03/25/2007

BF fit for noFSR

dhad bf noFSR

03/26/2007

BF fit for noFSR

dhad bf noFSR

Compare the difference

dhad table compare bf noFSR

BF fit for widede

Redo the widede yield selection:

dhad yield widede -t s —qsub Log

dhad yield widede -t d —qsub Log

dhad yield widede -t s —tag d —qsub Log

dhad yield widede -t d —tag d —qsub Log

dhad fit widede -t s —qsub Log

dhad fit widede -t s —tag d —qsub Log

dhad fit widede -t d —qsub Log

dhad fit widede -t s —tag d —qsub Log

03/27/2007

Prepare result for the PRD table 1 and 2

dhad table data_results_prd

dhad table data_results_bf_ratio_prd

Check the noFSR BF fitter

Make sure the efficiency matrix is correct

Update all the noFSR related tables:

dhad table double_signal_eff dhad table noFSR_single_signal_eff done.

dhad bf noFSR fit1 fit2

dhad table compare bf noFSR

Widede plots update

dhad yield widede noFSR -t s —qsub Log

dhad fit widede noFSR -t s —qsub Log

Make sure the double tag parameters are standard values:

dhad fit widede -t s —tag d —qsub Log

dhad fit widede -t d —tag d —qsub Log

Redo the noFSR

Make sure all the parameters in the RooFit script had been restored back to the standard values.

dhad yield noFSR -t s —qsub Log

dhad yield noFSR -t s —tag d —qsub Log

dhad fit noFSR -t s —qsub Log

dhad fit noFSR -t s —tag d —qsub Log

03/28/2007

Check Fitting

dhad fit P 0.5 -t s —tag d —qsub Log

create plots page done.

No FSR Update

update plots done. update tables done. update bf fits done. update bf tables done.

BF fit for widede

dhad table single_widede

dhad bf widede

dhad table compare bf widede

   1  N1           4.04736e+04   3.19132e+02   8.67800e-04  -2.98980e+00
   2  N2           4.05257e+04   3.20152e+02   8.67697e-04  -2.49239e+00
   3  Nbkgd1       2.96811e+04   3.01742e+02   9.17166e-04   1.03330e+00
   4  Nbkgd2       2.94441e+04   3.02345e+02   9.18908e-04   1.76237e+00
   5  md           1.86968e+00   4.12122e-05   6.56069e-03  -3.82123e-01
   6  p            1.50000e+00   1.07032e-02   1.64041e-01  -3.65312e-02
   7  sigmap1      2.65865e-02   2.01352e-04   8.61981e-03  -9.21537e-01
   8  xi          -1.00004e-01   4.70102e-01   1.17571e-01** at limit **

Increase the upper limit for xi:

    xi= new RooRealVar("xi","xi",-20.0,-100.0, 0.);

dhad fit widede -t d -m203  --qsub

Log

   1  N1           4.05064e+04   2.80912e+02  -3.70950e-05  -2.41693e+00
   2  N2           4.05586e+04   2.84801e+02  -4.53726e-05   7.21370e-01
   3  Nbkgd1       2.96476e+04   2.59927e+02   6.84385e-05  -5.83369e+00
   4  Nbkgd2       2.94119e+04   2.67629e+02   4.54519e-05   2.90846e+00
   5  md           1.86967e+00   3.82788e-05   2.54253e-04  -3.96088e-01
   6  p            1.50000e+00   6.83137e-03   1.67194e-03** at limit **
   7  sigmap1      2.65986e-02   1.80666e-04   2.67232e-04   2.54061e-01
   8  xi          -4.99910e-06   2.00468e-01  -1.00705e-03** at limit **



    xi= new RooRealVar("xi","xi",-20.0,-100.0, 2.);
dhad fit widede -t d -m203  --qsub

Log

 FCN=-2.3778e+06 FROM MIGRAD    STATUS=FAILED        980 CALLS         981 TOTAL
                     EDM=80.3241    STRATEGY= 1      ERR MATRIX NOT POS-DEF
  EXT PARAMETER                APPROXIMATE        STEP         FIRST
  NO.   NAME      VALUE            ERROR          SIZE      DERIVATIVE
   1  N1           4.09943e+04   1.91486e+02   8.65057e-04   2.00879e+01
   2  N2           4.10600e+04   1.92334e+02   8.64935e-04   5.21147e+01
   3  Nbkgd1       2.91633e+04   1.77830e+02   9.18760e-04  -2.78834e+00
   4  Nbkgd2       2.89113e+04   1.78461e+02   9.20565e-04  -5.15998e+01
   5  md           1.86963e+00   4.86020e-05   6.56718e-03   1.15555e+00
   6  p            1.49587e+00   1.14336e+00   1.75842e-01   1.09795e+00
   7  sigmap1      2.67999e-02   1.13508e-04   8.60596e-03  -5.00747e+00
   8  xi           1.60273e+00   2.32104e+01   8.95584e-02  -6.80303e+00

Use all the standard initial values as input:

   RooRealVar N1("N1","N1",8164.0,0.0,2000000.0);
   RooRealVar N2("N2","N2",8168.0,0.0,2000000.0);

   RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",41312.0,0.0,2000000.0);
   RooRealVar Nbkgd2("Nbkgd2","Nbkgd2",40821.0,0.0,2000000.0);

   xi= new RooRealVar("xi","xi",-10.0,-100.0,-0.1);

dhad fit widede -t d -m203  --qsub

Log

Record the log file done.

Restore it to the standard stage done.

dhad bf widede

dhad table compare bf widede

BF fit for Gamma = 22.5 MeV

dhad bf Gamma 0.0225

dhad table compare bf Gamma/0.0225

BF fit for Gamma = 32.5 MeV

dhad bf Gamma 0.0325

dhad table compare bf Gamma/0.0325

BF fit for P = 0.5

dhad bf P 0.5

dhad table compare bf P/0.5

Create the fit Results Data Variations table

dhad table data_results_variations

Implement the get_sys_err func in DHadBF

dhad bf nofit1 nofit2