05/14/2007
Make Default version 7.05
DHad.py
Get Yields
- Reprocess signal MC
dhad yield regular -t s —qsub
Problem in qsub:
failed opening input/output file because: 05/14/2007 10:08:21 [1821:31141]: error: can't open output file "/home/xs32/work/CLEO/analysis/DHad/ Use "qmod -c <jobid>" to clear job error state once the problem is fixed.
Uncomment lines in qsub() OK.
dhad yield regular -t s —qsub
dhad yield regular -t s —tag d —qsub
- As data and generic don't change, just link them
cd /home/xs32/work/CLEO/analysis/DHad/7.05/dat/evt/regular ln -sf /home/xs32/work/CLEO/analysis/DHad/7.03/dat/evt/regular/data_* . ln -sf /home/xs32/work/CLEO/analysis/DHad/7.03/dat/evt/regular/generic_* .
Post the chisq comparison on HN
chisq comparison for R=0 and R=7.5 Dear all, As David suggested on the last group meeting, I fitted the single tags for data with R=0 and R=7.5 respectively. The related plots can be found at: (Under the section "Single Tag") http://www.lepp.cornell.edu/~xs32/private/DHad/7.03/plots/ I also made tables comparing their chisq with the regular fits which has R=15: (Under section "Comparison" : "chisq") http://www.lepp.cornell.edu/~xs32/private/DHad/7.03/tables/ The trend looks like what we expected, R=7.5 has closer chisq with R=15. Xin
Fit for signal MC
dhad fit floatmass fitparams diag -t s —tag d —qsub
05/15/2007
Plot the params
Create table
dhad table line_shape_parameters
- For Mode 202, we used the version 7.03 value as initial value
s1a=new RooRealVar("s1a","s1a",2.25,1.5,4.0); f1a=new RooRealVar("f1a","f1a",0.166,0.0,1.0); s1b=new RooRealVar("s1b","s1b",4.00,1.5,4.0); f1b=new RooRealVar("f1b","f1b",0.049,0.0,0.4);
dhad fit floatmass fitparams diag -t s —tag d -m202,202 —qsub
Retriev RCS1.1 after done the fit
05/16/2007
Restore the double tag fit code to normal (1.1)
OK.
Message from BES
Thank you for your message. The "r" from the fit is 41+/-458 fm from the fit mentioned in the paper (hep-ex/0612056). Howevere, the "r" value from the different fits change a lot. From the fit to the data, we found that the other parameters are not sensitive to the "r" value. If we fixed the "r" at a value from several fm to several thousand fm in the fit, the quantities M_{psi(3770)}, Gamma_{psi(3770)}, Gamma^{ee}_{psi(3770)} and R_{uds} are almost unchanged. By the way, Did CLEO-c also make the cross section scan over the psi(3770) resonance? How may energy points at which CLEO-c collected the data? And what are the energies? best regards, Rong Gang
lineshpae para table
dhad table line_shape_parameters
Single tag fit
dhad fit regular -t s —qsub
dhad fit regular -t d —qsub
dhad fit regular -t s —tag d —qsub
dhad fit regular -t d —tag d —qsub
05/17/2007
Plots for regular fit
done.
Tables for regular fit
dhad table single_signal_eff
dhad table compare version 7.03 single_signal_eff 'Efficiency(%)' Indp
NoFSR yield
dhad yield noFSR -t s —qsub
dhad yield noFSR -t s —tag d —qsub
R= 0 and 7.5 for double tag fit
dhad fit R 0 -t d —tag d —qsub Log
dhad fit R 7.5 -t d —tag d —qsub Log
05/23/2007
Link the Old Dalitz file
signal/olddalitz
dhad yield olddalitz -t s —qsub Log
dhad yield olddalitz -t s —tag d —qsub Log
05/24/2007
Fit Old Dalitz MC
dhad fit olddalitz -t s —qsub
dhad fit olddalitz -t s —tag d —qsub
Plots done.
Tables
dhad table single_signal_eff olddalitz
dhad table compare column olddalitz_single_signal_eff,'Efficiency(%)' 7.03/single_signal_eff
Fit the file with the correct evt file
dhad fit olddalitz -t s —qsub
dhad fit olddalitz -t s —tag d —qsub
05/25/2007
Old Dalitz fit
Refresh plots done.
Table:
dhad table single_signal_eff olddalitz
dhad table compare column olddalitz_single_signal_eff,'Efficiency(%)' 7.03/single_signal_eff
dhad table compare column olddalitz_single_signal_eff,'Efficiency(%)' single_signal_eff
done.