D Had Log - November 2006

11/01/2006

Put scatter on the hypernews

done.

11/02/2006

Update PRD draft

• Go over through Section IV. "Data Samples and Event Selection". There are several places need David Cassel's input.

• Section IX Systematics

DHad Meeting with the paper committee

• pi0 mass energy bias? How while we can mock up the shower, observation profile has changed?

11/06/2006

Post link of the comparison tables on the DHadReviews hypernews

11/07/2006

Message from Anders

HyperNews

  this is not an official summary of the discussion we had last
Thursday, but it is a list of things that we need to respond to.

  1) We discussed lineshapes at some length. There are several
     issues here.

Redo the comparison tables

• No error, compare with standard

Done. Put them on the HyperNews.

Explain table 5 and 8 the sigma in the caption

check in the CBX (version 1.44)

11/08/2006

In MC sample rescale the pi0 momentum by 0.5% and recalculate mBC.

Message from Anders:

Just take the 3-momentum of the pi0 that you reconstructed.
Multiply by 1.005 (0.5%) for each momentum component. Use
this to recalculate the beam constrained mass.

Refere page of D Ntuple Meaning. Input ecm, pi0px[npi0], pi0py[npi0], pi0pz[npi0]

• Mbc for pi0

Plot DeltaE for the Kpipi0 mode for the different data sets

DHadPlot.py -a deltaE -o data31 DHadPlot.py -a deltaE -o data32 DHadPlot.py -a deltaE -o data33 DHadPlot.py -a deltaE -o data35 DHadPlot.py -a deltaE -o data36 DHadPlot.py -a deltaE -o data37

11/09/2006

Check the plots

• delta E for different data sets of mode 1

11/10/2006

Fix The Mbc plot

DHadPlot.py -o dir_mbcpi0,s,s,1 mbc test

Refere to the code for Lorentz boost:

Code here.

11/11/2006

The Mbc plot

Refere to the defination of mBC:

DHadronic_log_200603: mbc

Plot pi0 mass in mode 1 "Kpipi0"

pi0mass[npi0]

11/13/2006

Plot mBC by calculation

DHadPlot.py -o dir_mbc_calc,s,s,1 mbc

• Mbc for pi0

DHadPlot.py -o dir_mbc_calc_pi0_1.005_test,s,s,1 mbc

11/14/2006

Plot mBC by calculation

DHadPlot.py -o dir_mbc,s,s,1 —qsub mbc job 41028

DHadPlot.py -o dir_mbc_calc,s,s,1 —qsub mbc job 40928

DHadPlot.py -o dir_mbc_diff_calc,s,s,1 —qsub mbc job 41047

DHadPlot.py -o dir_mbc_calc_pi0_1.005,s,s,1 —qsub mbc job 40935

DHadPlot.py -o dir_mbc_diff_calc_pi0_1.005,s,s,1 —qsub mbc job 41024

• Mbc for pi0

Select the mBC events

DHad.py -p dir_mbc_calc_pi0_1.005 -t s -m 1

Do the fit

DHadFit.py -p dir_mbc_calc_pi0_1.005 -t s -m 1 —qsub job 41366 done.

• Fitting result

Compare the difference

 DHadTable.py -a compare
 ~/work/CLEO/analysis/DHad/dat/fit/mbc/calc_pi0_1.005_signal_Single_D0_to_Kpipi0__D0B_to_Kpipi0.txt
 ~/work/CLEO/analysis/DHad/dat/fit/cleo200609/signal_Single_D0_to_Kpipi0__D0B_to_Kpipi0.txt ~/work/CLEO/analysis/DHad/dat/tab/compare_mbc_pi0.htm

• Fitting result

11/15/2006

Fix the table

11/19/2006

Fix the table

 DHadTable.py -a compare
 ~/work/CLEO/analysis/DHad/dat/fit/mbc/calc_pi0_1.005_signal_Single_D0_to_Kpipi0__D0B_to_Kpipi0.txt
 ~/work/CLEO/analysis/DHad/dat/fit/cleo200609/signal_Single_D0_to_Kpipi0__D0B_to_Kpipi0.txt
 ~/work/CLEO/analysis/DHad/dat/ta\b/compare_mbc_pi0.htm

• Fitting result

done.

11/20/2006

make a plot of the pi0 mass in Kpipi0

11/21/2006

make a plot of the pi0 mass in Kpipi0

Also there was a suggestion to make a plot of the pi0 mass

for the pi0's used in Kpipi0 events to see if there was a problem with the pi0 mass.

DHadPlot.py mass pi0 mode=1

11/22/2006

make a plot of the pi0 mass in Kpipi0

DHadPlot.py mass pi0 mode=1

11/25/2006

make a plot of the pi0 mass in Kpipi0

DHadPlot.py mass pi0 -m 1

11/28/2006

make a plot of the pi0 mass in Kpipi0

DHadPlot.py mass pi0 -m 1 --path=mass/pi0 --qsub

job 46033

11/29/2006

pi0 mass plot

• pi0 mass in Kpipi0

plot the kpi, kpipi mode to see the logL distribution of fa

• find the output of logL in Minuit

DHadFit.py fa  -m 1

Warning: wrong member access operator '.' FILE:/home/xs32/work/CLEO/analysis/DHad/script/Roofit/mbc_singletag_3s.C LINE:735
Error: Symbol fitres is not defined in current scope  FILE:/home/xs32/work/CLEO/analysis/DHad/script/Roofit/mbc_singletag_3s.C LINE:735
Error: Failed to evaluate fitres.minNll*** Interpreter error recovered ***

Hit fraction shrink to 60%(70,80,90) in the Ntuple

• variable in the ntuple: trhitf[ntrack]
• Select events for 60% print it out

11/30/2006

plot the kpi, kpipi mode to see the logL distribution of fa

DHadFit.py fa  -m 1

The Gas-Wire Vapor monte carlo is in

/nfs/cor/an1/ponyisi/dhad/summerconf/signal/dtuple_c_0525_gwv_p2_1104

DHad.py yield gwv --tag d -m 0,0

• Do the fit

DHad.py fit gwv --tag d -m 0,0

Nov 30, 2006