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