D Had Log - February 2007

02/01/2007

Put the current contents into CBX

02/02/2007

Refine table

02/05/2007

Refine table

DHad.py table compare output mbc/calc/pi0/0.997 -t s -m 1

Update CBX table for this

DHad.py table compare output mbc/calc/pi0/0.997 -t s -m 1 —cbx

DHad.py update cbx table /home/xs32/work/CLEO/analysis/DHad/dat/tab/compare_output_mbc_calc_pi0_0.997.tex

02/06/2007

Update table in CBX

test

DHad.py update cbx table /home/xs32/work/CLEO/analysis/DHad/dat/tab/compare_output_mbc_calc_pi0_0.997.tex

02/07/2007

update the table

DHad.py table compare output mbc/calc/pi0/0.997 -t s -m 1 —label tab:mbc_pi0

issue:

• location not correct
• line/line/ need to be more flexiable

02/08/2007

update the table

DHad.py table compare output mbc/calc/pi0/0.997 -t s -m 1 —label tab:mbc_pi0

save to tmp file and do minor edit.

Input into CBX

02/09/2007

Update tables into CBX

DHad.py update cbx table

02/12/2007

Update tables into CBX

DHad.py update cbx table -v1

done.

02/13/2007

Update plots into CBX

DHad.py update cbx plot -v1

02/14/2007

Fit the new BES mass of psi(3770)

Message from Anders:

 As a crosscheck it was asked that we repeated the yield fits
 for alternative values of Gamma and the mass of the psi(3770). We
 should use the new BES numbers

 Gamma=28.6 MeV
 M=3772.2 MeV

 Could you run these fits?

• Double tag

DHad.py fit M 3.7722 Gamma 0.0286 -t d —tag d —qsub Log

DHad.py fit M 3.7722 Gamma 0.0286 -t s —tag d —qsub Log

DHad.py fit M 3.7722 Gamma 0.0286 -t g —tag d —qsub Log

• Single tag

DHad.py fit M 3.7722 Gamma 0.0286 -t d —qsub Log

DHad.py fit M 3.7722 Gamma 0.0286 -t s —qsub Log

DHad.py fit M 3.7722 Gamma 0.0286 -t g —qsub Log

02/15/2007

Check the fits

Done.

Make table of the new BES mass of psi(3770)

DHad.py table compare yields M/3.7722/Gamma/0.0286

Fit the BES mass of psi(3770) 3.7724

DHad.py fit M 3.7724 Gamma 0.0286 -t d —qsub Log

Fit the data with mass = 3.770

DHad.py fit M 3.77 -t d —qsub Log

Fit the psi(3770) width with different mass

DHad.py fit M 3.77 Gamma float -t d —qsub Log

DHad.py fit M 3.7724 Gamma float -t d —qsub Log

DHad.py fit Gamma float -t d —qsub Log

02/19/2007

Plot the fit of psi(3770) with 3.7724

done.

Make table of the mass of psi(3770) 3.7724

DHad.py table compare yields M/3.7724/Gamma/0.0286 -t d

done.

Fit the data with mass = 3.770

Make plots page ... done. Make table page ...

DHad.py table compare yields M/3.77 -t d

done.

Fit the psi(3770) width with different mass

Make plots page ... done. Make table page ...

DHad.py table psimass

02/21/2007

Fits from Peter

Evt files location:

/nfs/cor/user/ponyisi/daf9/hadD/summerconf/kspipi0/signflip/
    data_Kspipi0_multiple_candidate.evt
    data_Kspipi0_single_candidate.evt
    mc_Kspipi0_multiple_candidate.evt
    mc_Kspipi0_single_candidate.evt

DHad.py fit signflip multiple -m 203 -t d —qsub Log

DHad.py fit signflip single -m 203 -t d —qsub Log

DHad.py fit signflip multiple -m 203 -t s —qsub Log

DHad.py fit signflip single -m 203 -t s —qsub Log

Fits for the Kspipi0_UPDDALITZ

Message from Peter:

The "signflip" MC (new Dalitz structure) for Ks pi pi0 is available in
/nfs/cor/an1/ponyisi/dhad/summerconf/signal/dtuple_c_0525_p2_1104, in
all the files with "Kspipi0_UPDDALITZ" in the name.

DHad.py yield Kspipi0_UPDDALITZ -t s —test

02/22/2007

Orgnize the plots

done.

Table of the psi(3770) width with different mass

DHad.py table psiwidth

02/23/2007

Yields for the Kspipi0 UPDDALITZ

DHad.py yield Kspipi0 UPDDALITZ -t s

DHad.py yield Kspipi0 UPDDALITZ -t s —tag d

done.

Fits for the Kspipi0 UPDDALITZ

• Get the parameters

DHad.py fit floatmass fitparams diag Kspipi0 UPDDALITZ -t s —tag d —qsub

Log

DHad.py fit Kspipi0 UPDDALITZ -t s —test

02/26/2007

Fits for the Kspipi0 UPDDALITZ

plots for parameters ... done.

table for the parameters ...

DHad.py table line_shape_parameters Kspipi0 UPDDALITZ done.

DHad.py fit Kspipi0 UPDDALITZ -t s —qsub

Log

Alternative psi(3770) lineshape

Message from Anders:

HN link

 As we discussed in the DHad meeting today BES used a slightly
 different lineshape in their paper paper (hep-ex/0612056)

  I have implemented this in:

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

Xin, can you take this version and rebuild RooFitModels and try to do
the fit for the single tag yields where you fit for the width.

  To use the new lineshape you have to set the parameter, which I
think in your script is called 'mc' to 3.

  In my implementation I assumed that the psi(3770) decays only to
DDbar.
 

cd ~/work/CLEO/analysis/DHad/src/RooFitModels/
gmake -f GNUmakefile.standalone

Error log

02/27/2007

Alternative psi(3770) lineshape

log on lnx134 to compile:

ssh lnx134
cd ~/work/CLEO/analysis/DHad/src/RooFitCore/
gmake -f GNUmakefile.standalone clean
gmake -f GNUmakefile.standalone

cd ~/work/CLEO/analysis/DHad/src/RooFitModels/
gmake -f GNUmakefile.standalone clean
gmake -f GNUmakefile.standalone

Using Anders RooFitModels code, compile completed.

DHad.py fit BES2006lineshape Gamma float -t d —qsub

Log

Out put the root file for the plots

DHad.py fit regular -t d —tag d —qsub Log

DHad.py fit regular -t s —tag d —qsub Log

DHad.py fit regular -t g —tag d —qsub Log

DHad.py fit regular -t d —qsub Log

DHad.py fit regular -t s —qsub Log

DHad.py fit regular -t g —qsub Log

02/28/2007

Out put the root file for the plots

update plots link ... done.

post on hyperNews

Alternative psi(3770) lineshape

create the plots page ... done.

update the table for comparison ... done.

psi(3770) width with different mass

Noticed problem:

• Mode 201: kpipipi0 Gamma = 31.8+/- 23.9 STATUS=NOT POSDEF

Try different initial value for gamma:

DHad.py fit BES2006lineshape Gamma float -t d -m 201 —init Gamma=31.8 —qsub Log

• Mode 204: Kspipipi Gamma = 60.0 +/- 0.0 STATUS= NOT POSDEF

DHad.py fit BES2006lineshape Gamma float -t d -m 204 —init Gamma=25 —qsub Log

Gamma = 31.4 +/- 19.7

Another try: DHad.py fit BES2006lineshape Gamma float -t d -m 204 —init Gamma=31.4 —qsub Log

Also problem in the regular fit:

• Mode 203: Kspipi0 Gamma = 20.7 +/- 19.8

DHad.py fit Gamma float -t d -m 203 —init Gamma=15 —qsub Log Gamma= 20.6 +/- 9.5 OK.

Add footnote on plots page ... done.

• Mode 205: Kpipi

Gamma = 29.1 +/- 44.9

DHad.py fit Gamma float -t d -m 205 —init Gamma=20 —qsub Log Gamma = 29.2 +/- 44.9

Try again: DHad.py fit Gamma float -t d -m 205 —init Gamma=29.2 —qsub Log Gamma = 27.0 +/- 4.6 OK.

Add notes on plots page ... done.

Set R=0 for the fits and compare the width

DHad.py fit R 0 Gamma float -t d —qsub Log

Fits for the Kspipi0 UPDDALITZ

create plots page ... done. create table to compare the yields with regular ...

DHad.py table compare yields Kspipi0/UPDDALITZ -t s done.

create the process page ... done. Send message on the HyperNews.

Redo the yield to get log files

DHad.py yield -t d —qsub Log

DHad.py yield -t s —qsub Log

DHad.py yield -t g —qsub Log

DHad.py yield -t d —tag d —qsub Log

DHad.py yield -t s —tag d —qsub Log

DHad.py yield -t g —tag d —qsub Log