D Hadronic Analysis - Log - April 2006

04/05/2006

When p= 0.5, save the eps file as "p0.5_***.eps".

Signal MC:

c3python DHadFit.py -a 1 -t s -m 0 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 1 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 3 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 200 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 201 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 202 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 203 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 204 --tag s -p 0.5
c3python DHadFit.py -a 1 -t s -m 205 --tag s -p 0.5

Plots for MC

Data:

c3python DHadFit.py -a 1 -t d -m 0 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 1 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 3 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 200 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 201 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 202 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 203 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 204 --tag s -p 0.5
c3python DHadFit.py -a 1 -t d -m 205 --tag s -p 0.5

Plots for Data fix p=0.5

Read out the numbers, including the errors.

Better to go from the output of the fit file. Ask Peter.

Reply, from the .c file.

• Changed Gamma(psi(3770)) from the default of 30.6 MeV to 23.6 MeV and 35.0 MeV.

Look at : mbc_singletag_3s.c

Signal MC

c3python DHadFit.py -a 1 -t s -m 0   --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 1   --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 3   --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 200 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 201 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 202 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 203 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 204 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t s -m 205 --tag s -g 0.0236

Data:

c3python DHadFit.py -a 1 -t d -m 0   --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 1   --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 3   --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 200 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 201 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 202 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 203 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 204 --tag s -g 0.0236
c3python DHadFit.py -a 1 -t d -m 205 --tag s -g 0.0236

Plots Gamma 23.6 MeV

Use qsub, see Ryd's example:

/nfs/cor/user/ryd/RooFit/01-09-02

more setenv.txt
c3rel 20050417_FULL
setenv LD_LIBRARY_PATH $C3_OTHER/Root/lib/

more sub_double
...

qsub  -l arch=lx24-x86  -j y -o /cdat/dafj/ryd/RooFit/01-09-02 -N run_dt_d0_all_data.log run_dt_d0_all_data.scr

more run_doubletagfit
c3rel 20050417_FULL
hostname
uname -a
cd /cdat/dafj/ryd/RooFit/01-09-02
root -b << EOF
.x lineshapefit2d.c("${EVTFILE}","${MODE1}","${MODE2}","${PSFILE}",${LINESHAPE},"${OPT}")
.q
EOF

Ryd's post on hypernews:

 we had noted that some of the double tag fits was not looking
very good on the high side tail. I have redone some fits where
I have used different widths for the charged and neutral modes
to get a better fit. I used a psi(3770) width of 20 MeV for the
D+D- and 27 MeV for the D0D0B. You can see the plots where all
modes are summed here:

• Fit data for D0D0B mode 27 MeV
• Data D+D- 20 MeV

04/06/2006

Keep fitting by using the sigle machine at the same time:

• 35.0 MeV

c3python DHadFit.py -a 1 -t s -m 0   --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 1   --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 3   --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 200 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 201 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 202 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 203 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 204 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t s -m 205 --tag s -g 0.0350

Data:

c3python DHadFit.py -a 1 -t d -m 0   --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 1   --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 3   --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 200 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 201 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 202 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 203 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 204 --tag s -g 0.0350
c3python DHadFit.py -a 1 -t d -m 205 --tag s -g 0.0350

Plots Gamma 35.0 MeV

- Fit data for D0D0B mode 27 MeV

c3python DHadFit.py -a 1 -t d -m 0   --tag s -g 0.027
c3python DHadFit.py -a 1 -t d -m 1   --tag s -g 0.027
c3python DHadFit.py -a 1 -t d -m 3   --tag s -g 0.027

Plots Gamma 27 MeV

• Data D+D- 20 MeV

c3python DHadFit.py -a 1 -t d -m 200 --tag s -g 0.020
c3python DHadFit.py -a 1 -t d -m 201 --tag s -g 0.020
c3python DHadFit.py -a 1 -t d -m 202 --tag s -g 0.020
c3python DHadFit.py -a 1 -t d -m 203 --tag s -g 0.020
c3python DHadFit.py -a 1 -t d -m 204 --tag s -g 0.020
c3python DHadFit.py -a 1 -t d -m 205 --tag s -g 0.020

Plots Gamma 20 MeV

• Create my_mbc_singletag_3s.C

Test

c3python DHadFit.py -a 1 -t d -m 205 --tag s -g 0.01

04/07/2006

• Make the output table straight

  mc = 0  BW shape
  mc = -1 "Broken" MC
  mc = 1  MARK III

c3python DHadFit.py -f test

Success.

• Resubmit the fits for Data

Data: Gamma 35.0 MeV

Test.

c3python DHadFit.py -a 1 -t d -m 0   --tag s -g 0.0350

Good.

qsub g35.sh

Your job 342124 ("g35.log") has been submitted.

Data: Gamma 23.6 MeV

qsub g236.sh

Your job 342125 ("g236.log") has been submitted.

Fit data for D0D0B mode 27 MeV

qsub g27.sh

Your job 342126 ("g27.log") has been submitted.

Data D+D- 20 MeV

qsub g20.sh

Your job 342127 ("g20.log") has been submitted.

04/10/2006

Get the fitting results table:

• Data: Gamma 35.0 MeV

Look at the log file "g35.log342124", it's OK. Indicating that the fit was successful.

Read out the result table from the files:

Data fit Gamma 35.0 MeV

04/11/2006

Embed the lisp in the muse file to get the table directly from the txt file.

• Data: Gamma 23.6 MeV

"g236.log342125"

Data fit Gamma 23.6 MeV

• Data: D0D0B mode Gamma 27 MeV

Data D0D0B fit Gamma 27 MeV

• Data D+D- 20 MeV

Data D+ D- fit Gamma 20 MeV

04/13/2006

• Compare the results from the tables

DHadronic_Tables

04/14/2006

• Check out the cvs

From the DHad HyperNews:

I have set up a CVS repository to keep the cbx draft for the 281/pb
analysis. To check out the cbx do:

setenv CVSROOT /home/ryd/CVSROOT
cvs co dhadcbx281

when you do this you get the tex source for the 56/pb
analysis - but without any .eps files. This version
does not yet tex. As a first step we should copy in as
many of the figures as we have from the fits into the
new CVS directory.

If you want to check out the cbx on a machine that does
not mount the lab disks, such as a laptop you can do:

ryd@mp126:~> export CVS_RSH=ssh
ryd@mp126:~> export CVSROOT=ryd@lnx134.lepp.cornell.edu:/home/ryd/CVSROOT

cdDHad
cd doc
export CVSROOT=/home/ryd/CVSROOT
cvs co dhadcbx281

Open "dhadcbx.tex"

04/19/2006

LaTeX can't process, due to the needed plots.

• Rerun the regular fits, to get the parameter table

Using Peter's evt file.

1. Data

Test:

c3python DHadFit.py -a 1 -t d -m 0  --tag s

Fix bug , go.

qsub data.sh
Your job 355691 ("data.log") has been submitted.

2. Signal

Test:

c3python DHadFit.py -a 1 -t s -m 0  --tag s

Good.

qsub sigal.sh
Your job 355692 ("signal.log") has been submitted.

3. Generic

Test:

c3python DHadFit.py -a 1 -t g -m 0  --tag s

No files in Peter's area. So use mine.

Test-2:

c3python DHadFit.py  -t g -m 0  --tag s

Good.

qsub generic.sh
Your job 355694 ("generic.log") has been submitted.

04/20/2006

• Add tables in the regular fit

  Done.

• Look at the cvs tex file, without compiling the figures.

  Works.

  dhadcbx.ps

04/25/2006

• Get the spreadsheet ( local | remote )

04/27/2006

Signal MC : Regular fit- signal MC: Update done.

For the data : Regular fit- data Update done.

G=23.6 MeV : Data fit Gamma 23.6 MeV Update done.

G=35.0 MeV : Data fit Gamma 35.0 MeV Update done.

p=0.5 : Fix p=0.5 fit Update done.

04/28/2006

• Finish the spreadsheet updating

1. Resub the p=0.5 fit

Test:

c3python DHadFit.py -a 1 -t d -m 0 --tag s -p 0.5

OK.

Save qsub job list: p05.sh

qsub p05.sh
Your job 363010 ("p05.log") has been submitted.

2. updating the spreadsheet. Done. (use the previous output)

• Input in the CBX Table 8: "Single tag efficiencies in our signal Monte Carlo."

  o D mass in the table dose not have suffecient digits. Ask in the Hyper News. For now, extract from the eps temperialy.

python readeps.py -f ../../doc/fig/data_Single_D0AndD0B_to_Kpi.eps

• 05/01/2006