D Hadronic Analysis - Log - May 2006

• 04/28/2006

05/01/2006

• Add more digits in the fitting results.

in my_mbc_singletag_3s.c:

table.precision(10);

Test one:

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

Works.

Resubmit all the fits.

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

qsub g236.sh
Your job 364261 ("g236.log") has been submitted.

qsub g35.sh
Your job 364262 ("g35.log") has been submitted.

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

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

qsub g20.sh
Your job 364266 ("g20.log") has been submitted.

qsub g27.sh
Your job 364267 ("g27.log") has been submitted.

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

• Inserting the plots.

  \label{fig:singletag_sig_yields_d0}
  \label{fig:singletag_sig_yields_dp1}
  \label{fig:singletag_sig_yields_dp2}
  \label{fig:singletag_gen_yields_d0}
  \label{fig:singletag_gen_yields_dp1}
  \label{fig:singletag_gen_yields_dp2}
  \label{fig:singletag_data_yields_d0}
  \label{fig:singletag_data_yields_dp1}
  \label{fig:singletag_data_yields_dp2}

CVS check in the figures. Version 1.3

05/02/2006

Data mode 202 dosen't have the correct digits. Look at the "data.log.364260"

/cdat/daf9/ponyisi/hadD/summerconf/txt/data_Single_Dp_to_Kpipi.evt is not ready.

Check, it exists!

 more /cdat/daf9/ponyisi/hadD/summerconf/txt/data_Single_Dp_to_Kpipi.evt

So, resubmit this one for testing.

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

OK. Resubmit the data.sh

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

• Construct the core table content for CBX:

Element needed:

\caption{Single tag efficiencies in our signal Monte Carlo.

Mode Generated Yield mass $\xi$ $p$ efficiency

Python code:

DHadTable.py

05/03/2006

• Check the data fit.

mode 204 , 205 still not updated. So, fit them again not use qsub.

Relink the need libs:

DHadSoftlink.py

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

• DHadTable.py

DHadTable.py -f ../../dat/fit/signal_Single_D0AndD0B_to_Kpi.txt

05/04/2006

• DHadTable.py ready for the project:

  Single tag efficiencies in signal Monte Carlo

DHadTable.py -a 1  -o ~/work/CLEO/analysis/DHad/dat/tab/signal_Single_yield_mass_xi_p.txt

; INVALID LISP CODE

05/08/2006

• Edit the pre-tex table

doc/tab/20060405.sxc

05/09/2006

Input generated mc number in the table dat/tab/signal_Single_gen.txt

Combine the signal_Single_yield_mass_xi_p.txt with signal_Single_gen.txt into signal_Single_input.txt

Now calculate the efficiency

DHadTable.py -a 2
65.72 +/- 0.27
66.78 +/- 0.27
35.34 +/- 0.12
35.90 +/- 0.12
46.90 +/- 0.16
47.02 +/- 0.16
54.88 +/- 0.18
55.36 +/- 0.18
28.30 +/- 0.18
28.36 +/- 0.18
45.67 +/- 0.18
45.21 +/- 0.18
23.22 +/- 0.19
23.21 +/- 0.19
32.01 +/- 0.24
31.65 +/- 0.24
46.99 +/- 0.36
47.04 +/- 0.36

05/10/2006

• Get the final LaTeX format table

Create an one column table for the Mode title:

dat/tab/signal_Single_mode.txt

Keep this file seperate, only use for output style, do not merge it to the input file.

05/11/2006

• Create the whole txt table for muse format

05/12/2006

DHadTable.py -a 2

See Table

05/13/2006

Update the table in the cbx

DHadTable.py -a 2

Done

05/14/2006

• Get the table from the original output:

Data N1, N2 for all modes:

DHadTable.py -a 3

Resubmit the widede jobs for data and signal MC:

Test on widede_data:

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

O

No evt file in Peter's area. But have the widede_signal instead. So, resubmit them first.

Test widede_signal:

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

OK. Now submit the qsub job:

qsub widede_signal.sh
Your job 379687 ("widede_signal.log") has been submitted.

05/15/2006

Check the widede_signal table:

DHadTable.py -a 3 -t ws

OK.

While waiting for the widede_data, go on calculating the error.

DHadTable.py -a 3 -t ws

See Wide Signal with fitting errors.

DHadTable.py -a 3 -t s
DHadTable.py -a 3 -t d

>> Tabel Signal >> Tabel Data

O

Manually construct the table of N1: >> Tabel Delta E - N1

Error Propagation

05/16/2006

• Effeciency for Data

Start from the Table and Equation

DHadTable.py -a 4

?

eff > 1 for this mode:

204  	Kspipipi  	8375 +/- 127  	8289

Use abs(1-eff) for now:

                -------------
               / eff |1-eff|
  eff_err =   / -------------
            \/       N1'

DHadTable.py -a 4

Table Delta E - N1-out

• Add effeciency for Signal MC

DHadTable.py -a 4

Table Delta E - N1-out

05/17/2006

• Combine the errors from Data and Signal MC

     a = x / y

  sa^2      sx^2         sy^2
 ------ =  ------   +   ------
   a^2       x^2          y^2

05/18/2006

• Realize it in the code

              ---------------------
             /  sx^2         sy^2
 sa   = a * /  ------   +   ------
          \/     x^2          y^2

DHadTable.py -a 4

Table Delta E - N1-out

• For N2, Join the data stream by code:

Input files: data, wide data, (This is manually got from Yields Table) signal, wide signal

Table Delta E - N2

• N1, N2 ratio table

DHadTable.py -a 4

Table Delta E - N1 & N2

05/19/2006

• Search the variable for the final photon number in the ntuple

05/20/2006

Copy out the PueoUtils.py

05/22/2006

• Identify the variable name: Means of the HadronicDNtuple variables

INTEGER variable trrichpnp[ntrack]
"Track RICH pion number of photons"

INTEGER variable trrichknp[ntrack]
"Track RICH kaon number of photons"

05/23/2006

Message from Peter:

The procedure I would use is:

In Monte Carlo events, the variables "mcdmode" and "mcdbmode" are
filled in the ntuple.  These variables encode the final state
daughters of the D and Dbar.  The allowed final state daughters are
K+, K-, K0, pi+, pi-, pi0, gamma, lepton, and other.  They are encoded
in the HadronicDNtupleProc::constructMCDmode() function.  In the code,
each power of 10 is a different particle, so "2" means 2 K-, "20"
means 2 K+, "11" means K-K+, etc.  Photons (which do not come from pi0
decays) are multiples of 1000000.  They come either from radiative
decays (eg eta' -> gamma pi pi) or from FSR.  If you restrict yourself
to the signal MC, where we know what we generated has no radiative
decays, the photons it will find come only from FSR.

• Plot mcdmode variable

  cut1 = ROOT.TCut("(mcdmode/10e6)%10==0") cut2 = ROOT.TCut("(mcdbmode/10e6)%10==0") cut3 = ROOT.TCut("(mcdmode/10e6)%10==0 && (mcdbmode/10e6)%10==0")

• Recheck the order in the message, the ones for the modes are the correct.
• Check the C++ code HadronicDNtupleProc::constructMCDmode()

In the code:

int HadronicDNtupleProc::constructMCDmode(const MCParticle& mcpart) {
// definition of decay mode:
// 1*(#K-) + 10*(#K+) + 100*(#K0) + 1000*(#pi-) + 10000*(#pi+)
// + 100000*(#pi0) + 1000000*(#gamma)
// + 10000000*(#leptons) + 100000000*(#anything else)

Next:

• get_decimal_at_position in the PueoUtils for the digits
• Place the FSR veto cuts in front of the chooseD
• Look at the yield.py for the widede cut for chooseD

05/24/2006

Test : Single_Dp_to_KKpi.root

With FSR:

DHadPlot.py -a 1
Total: 19988, Selceted: 9510, Ratio: 0.475785471283

Withour FSR:

Total: 19988, Selceted: 9336, Effeciency: 0.467080248149, FSR events: 531

Output ASCII file for this mode sign:

05/25/2006

Test in the DHad.py

DHad.py -a 3

One mode success.

All modes and sign success, except 204 and 205 Copy functions to the DHadTools, and then do 204 and 205. Done.

Do the fit for these modes:

Test,

DHadFit.py -t noFSRs -m 0 --tag s

OK. Now, set the qsub mode:

qsub noFSRs.sh
Your job 390399 ("noFSRs.log") has been submitted.

05/26/2006

Check the result of fitting.

Fatal Python error: GC object already tracked
/nfs/sge/root/default/spool/lnx187/job_scripts/390399: line 17: 21197 Aborted                 DHadFit.py -t noFSRs -m 0 --tag s

Add the .bashrc in the noFSRs.sh file. Resubmit job.

qsub noFSRs.sh
Your job 391140 ("noFSRs.log") has been submitted.

Check the log file, noFSRs.log.o391140, it works.

• Select the widede data in the qsub mode

Look at the yieldswidedeltae.py

/nfs/cor/user/ponyisi/hadD/summerconf/yieldswidedeltae.py

for mode in modes:
    modes[mode]['decutl'] *= 2
    modes[mode]['decuth'] *= 2

Add options in DHad.py ws : wide deltaE signal , test the signal first

DHad.py -t ws -m 0 --tag s --sign 1

Works now. Process all wide modes and sign with command:

DHad.py -a 4

Done.

• Edit the plots No FSR fit- signal MC

Calculate the table.

05/29/2006

• Construct the base table first

DHadTable -a 5

Master table

Final Version

Update the corresponding table 9 in the cbx.

05/31/2006

Double Tag fit

dir : /nfs/cor/user/ponyisi/hadD/summerconf/fits qsub : doubdiagjobs.sh command : c3python driver.py —doubdiag 205 data

Start from the drive.py

DHadFit.py -a 2

Procedure OK.

Float the width of the 3770 in Single tag

DHadFit.py

• 06/01/2006