03/01/2007
Alternative psi(3770) lineshape
- Mode 201: kpipipi0 Gamma = 31.8+/- 23.9
Try different initial value for gamma:
DHad.py fit BES2006lineshape Gamma float -t d -m 201 —init Gamma=31.8 —qsub Log
Gamma = 3.3 +/- 0.0
Try again: (Use the regular fit's value: 28.7 +/- 3.5 )
DHad.py fit BES2006lineshape Gamma float -t d -m 201 —init Gamma=28.7 —qsub
Still Gamma = 3.3 +/- 0.0.
- Raise this error on group meeting
Set R=0 for the fits and compare the width
update the process page ... done. psi(3770) width with different mass
Problem on mode 204
mode 204: chisq = 3424, Gamma = 7.6 +/- 0.4 (regular is 25.7 +/- 5.3)
Try again for this mode:
DHad.py fit R 0 Gamma float -t d -m 204 —init Gamma=25.7 —qsub Log
Gamma = 28.6 +/- 3.0 OK.update the plots page ... done.
Problem in mode 201
Gamma = mode 33.6 +/- 41.3
Try again: (Use the regular fit's value: 28.7 +/- 3.5 )
DHad.py fit R 0 Gamma float -t d -m 201 —init Gamma=28.7 —qsub Log
Gamma = 33.8 +/- 4.4 OK. update the plots page ... done.
Kspipi0 UPDDALITZ double tag fit
DHad.py fit Kspipi0/UPDDALITZ -t s —tag d —test
03/03/2007
Try different limit for the lineshape parameter Kspipi0
NO. NAME VALUE ERROR SIZE DERIVATIVE 2 Nbkgd1 2.31375e-06 1.11844e+00 1.61536e-04** at limit ** 3 Nbkgd2 1.47240e-05 1.53198e+00 1.88840e-04** at limit ** 4 Nbkgd3 2.53729e-06 2.11116e+00 2.21608e-04** at limit **
Do it again:
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub
BES lineshape with BES mass fit
DHad.py fit M 3.7724 BES2006lineshape Gamma float -t d —qsub
Create yields table
DHad.py table psiwidth_yields
03/05/2007
Check the fit for lineshape parameter Kspipi0
Background at limit.
2 Nbkgd1 2.31375e-06 1.11866e+00 3.23071e-05 -1.57080e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
3 Nbkgd2 1.47240e-05 1.53191e+00 3.77679e-05 -1.57079e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
4 Nbkgd3 2.53729e-06 2.12484e+00 4.43215e-05 -1.57080e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
Think about this ... Send out email ...
BES lineshape with BES mass fit
- Mass = 3.7724 GeV, BES 2006 lineshape, Gamma float
- psi(3770) width with different mass
- D Dbar yields of psi(3770) width float fit
Start the fitter job
Message from Peter
The scripts are in /nfs/cor/user/ponyisi/daf9/hadD/summerconf/werner_scripts . If I remember correctly, you want to execute the scripts newfit-cat_files.sh, newfit-cat_files_data.sh, or newfit-cat_files_data_statonly.sh (for generic MC, data with systematics, and data without systematics). These will produce files called newfit-generic.sh, newfit-data.sh, or newfit-data_statonly.sh. Then you can just run these files. The yield files are /cdat/daf9/ponyisi/hadD/summerconf/*_yields_for_werner. It includes single and double tags, in order K- pi+ yield K- pi+ stat error K- pi+ syst error K+ pi- yield K+ pi- stat error K+ pi- syst error .... In the double tags it goes (D0, D0bar) = (0, 0), (0, 1), (0, 3), (1, 0), ... The yield files are currently generated by the program yieldsprint.py in /nfs/cor/user/ponyisi/daf9/hadD/summerconf/. It is supposed to get the systematics right. It extracts the numbers from your postscript files.
Ref to 200609: 12
yieldsprint.py
DHad.py bf
done data.
03/06/2007
PRD draft
- Update table III in PRD
- Create the table III
DHad.py table dt_dz_eff_yield
DHad.py table dt_dp_eff_yield
DHad.py update prd table ...
lineshape parameter Kspipi0
Set the limit to -10 for the background and redo the fit
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub
1 N 1.18215e+03 3.42482e+01 -5.75263e-07 4.73690e-02 2 Nbkgd1 -1.00000e+01 1.42359e+00 1.88376e-05** at limit ** 3 Nbkgd2 -1.00000e+01 2.07156e+00 -9.49291e-06** at limit ** 4 Nbkgd3 -2.31247e+00 1.93438e+00 2.60321e-06 -7.91374e+00 5 NbkgdFlat 1.28760e+01 3.76811e+00 -3.02785e-06 -8.48819e-01 6 f1a 5.00000e-01 7.14758e-02 2.86767e-03** at limit ** 7 f1b 8.16379e-02 1.48507e-02 -2.50913e-04 3.60745e-03 8 md 1.86943e+00 4.70498e-05 -6.99567e-06 -2.32592e-03 9 s1a 1.91329e+00 3.34299e-01 3.61341e-03 7.97226e-03 10 s1b 3.42021e+00 2.14689e-01 -1.43682e-03 -1.30087e-02 11 sigmacommon1 5.48365e-03 7.06814e-04 -2.63099e-03 6.31133e-03
Set the initial value of the three bkgd as 0 and do the fit again:
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub Log
1 N 1.18208e+03 3.48177e+01 4.27241e-06 -2.86968e+00 2 Nbkgd1 -1.00000e+01 1.27635e+00 -2.42859e-05** at limit ** 3 Nbkgd2 -9.99992e+00 1.94156e+00 2.46771e-06** at limit ** 4 Nbkgd3 -2.30555e+00 1.83440e+00 -3.10544e-06 -6.36368e-02 5 NbkgdFlat 1.28705e+01 3.74892e+00 9.25884e-06 -1.22654e+00 6 f1a 5.00000e-01 6.93141e-02 -4.25560e-02 3.50341e-03 7 f1b 8.15784e-02 1.46673e-02 -6.58390e-04 -1.40336e-01 8 md 1.86943e+00 4.70193e-05 8.89342e-05 4.80355e-02 9 s1a 1.91280e+00 3.24422e-01 -3.34618e-03 1.01698e-01 10 s1b 3.42088e+00 2.14612e-01 3.47272e-03 -3.03541e-02 11 sigmacommon1 5.48438e-03 6.98716e-04 -3.32137e-03 -1.10351e-01
Set the limit from -100, do the fit again:
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub Log
2 Nbkgd1 -1.00000e+02 1.77144e+00 4.92357e-05 -1.57081e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
3 Nbkgd2 -1.00000e+02 1.75287e+00 4.86243e-05 -1.57079e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
4 Nbkgd3 -5.31103e+00 1.87503e+00 4.56183e-07 -1.55704e+00
5 NbkgdFlat 9.95888e+00 6.87317e+00 2.20609e-06 -1.56633e+00
6 f1a 2.62314e-01 4.37111e-02 2.38053e-04 4.92750e-02
7 f1b 1.13476e-01 8.39115e-03 4.77344e-05 -4.47395e-01
8 md 1.86951e+00 9.65254e-07 1.63273e-02 3.06308e-01
9 s1a 2.73372e+00 1.07873e-01 7.45788e-05 -6.29621e+00
10 s1b 4.00000e+00 2.68303e-02 5.24613e-03 1.57036e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
Still dosen't work.
PRD draft
DHad.py update prd table
For the table II and III done.
lineshape parameter Kspipi0
Set the background to zero and do the fit:
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub Log
1 N 1.17033e+03 3.42633e+01 3.06012e-06 -1.52241e+00
2 NbkgdFlat 8.67462e+00 3.51082e+00 3.57563e-06 -1.56663e+00
3 f1a 5.00000e-01 4.90684e-01 5.00000e-01 1.57030e+00
WARNING - - ABOVE PARAMETER IS AT LIMIT.
4 f1b 9.23809e-02 2.51834e-02 1.29450e-03 -5.68176e-01
5 md 1.86943e+00 4.73854e-05 2.13147e-04 2.89945e-01
6 s1a 1.62077e+00 1.64603e+00 4.49239e-04 -1.12760e+00
7 s1b 3.17367e+00 3.43020e-01 7.15538e-04 3.45788e-01
8 sigmacommon1 6.05996e-03 3.76913e-03 1.40403e-03 2.25900e-01
Make a larger limit for f1a:
f1a=new RooRealVar("f1a","f1a",0.05,0.0,0.8);
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub Log
1 N 1.17034e+03 3.42635e+01 3.06013e-06 -1.52241e+00 2 NbkgdFlat 8.65870e+00 3.50610e+00 3.57616e-06 -1.56663e+00 3 f1a 7.50831e-01 1.16256e-01 8.76236e-04 1.06974e+00 4 f1b 9.25594e-02 1.88453e-02 1.23191e-03 -5.67118e-01 5 md 1.86943e+00 4.80474e-05 2.13131e-04 2.89629e-01 6 s1a 2.73778e+00 1.79633e+00 1.26034e-03 -9.77625e-03 7 s1b 3.46047e+00 3.44394e-01 8.88855e-04 2.53706e+00 8 sigmacommon1 3.29580e-03 2.90926e-03 3.06706e-04 -1.07941e+00
BF fitter
DHad.py bf
- signal_single_efficiencies
Sample file:
/nfs/cor/user/ponyisi/daf9/hadD/summerconf/signal_single_efficiencies_for_werner
03/08/2007
BF fitter
- signal_single_efficiencies
Meanings of
/cdat/daf9/ponyisi/hadD/summerconf/yields_and_efficiencies
Format is: for single tags: efficiency denominator, signal MC yield, efficiency for crossfeeds: signal MC yield/error on yield, efficiency denominator, signal MC yield, efficiency, efficiency error
DHad.py bf
done.
- signal_single_efficiencies_errors
03/09/2007
Fit for the single tag, get the yields
DHad.py fit Kspipi0 UPDDALITZ -t s —qsub Log
Try to fit the lineshape paras using the old value as initial
s1a=new RooRealVar("s1a","s1a",2.27,1.5,4.0);
f1a=new RooRealVar("f1a","f1a",0.27,0.0,0.5);
RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",8.9 ,0.0,2000000.0);
RooRealVar Nbkgd2("Nbkgd2","Nbkgd2",10.3,0.0,2000000.0);
RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",1.7 ,0.0,2000000.0);
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub
1 N 1.17034e+03 3.38223e+01 -0.00000e+00 7.97759e-01 2 Nbkgd1 6.72389e-07 7.81295e-01 0.00000e+00** at limit ** 3 Nbkgd2 3.83525e-05 1.40937e+00 -0.00000e+00** at limit ** 4 Nbkgd3 1.01348e-04 1.30106e+00 0.00000e+00** at limit ** 5 NbkgdFlat 8.67849e+00 3.23682e+00 0.00000e+00 1.02681e+00 6 f1a 4.98356e-01 3.69272e-01 0.00000e+00 -5.65548e-02 7 f1b 9.30278e-02 1.93838e-02 0.00000e+00 -6.16653e-02 8 md 1.86943e+00 4.53323e-05 -0.00000e+00 1.03228e-01 9 s1a 1.57476e+00 3.97410e-01 -0.00000e+00 -3.02225e-02 10 s1b 3.19563e+00 3.03789e-01 0.00000e+00 2.35167e-02 11 sigmacommon1 6.18433e-03 7.71990e-04 0.00000e+00 8.21824e-02
Fit the lineshape paras using the old value as initial and background set to zero
f1a=new RooRealVar("f1a","f1a",0.27,0.0,0.8);
RooRealVar Nbkgd1("Nbkgd1","Nbkgd1",0.0);
RooRealVar Nbkgd2("Nbkgd2","Nbkgd2",0.0);
RooRealVar Nbkgd3("Nbkgd3","Nbkgd3",0.0);
DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub Log
Backup old regular
DHad.py cp regular 20070309
Fit mode 0, R=0 fit
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
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
dhad fit regular1 -t s —tag d —qsub Log
R= 2.5 fit
- BES line shape
- Gamma 28.6 GeV
- Mass 3771.8 MeV
- R = 2.5
/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
- Create the master sys table by mode done.
- reproduce the TABLE VI in prd draft
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
- Fix mode 0
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
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
- Fix mode 3
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
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
- Fix
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
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
- Fix
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
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
- Fix
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
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
- Fix mode 203
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
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
- Fix
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
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
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
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
....
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
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
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
- png file dosen't show up in mode 204 plots, re transform the two plots to png done.
- Noticed fitting problem in data on mode 203
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
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
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
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