D Hadronic Analysis - Log - July 2006

07/03/2006

<<< 06/30/2006

Working on the single tag table

Single tag yields for generic MC

DHadTable.py -a 6

Using the New Gamma Values

Talking with Anders.

Based on the previous fits of the Mass = 3.774 GeV, we get the Gamma :

D0_to_Kpi 32.8 +/- 0.9 (+/- 5) MeV D+_to_Kpipi 28.5 +/- 0.8 (+/- 5) MeV

The +/- 5 MeV is due to the fact that the value shifted by 10 MeV when the mass changed for 2 MeV.

Now in the summary table of the vlues of psi(3770) width from existing measuremetns:

Sample	Gamma(psi-3770) (MeV)
MARK I	28 +/- 5
MARK II	24 +/- 5
DELCO	24 +/- 5
D0_to_Kpi	32.8 +/- 5 MeV
D+_to_Kpipi	28.5 +/- 5 MeV
Average	27.5 +/- 2.2 MeV

Vary Gamma 5 MeV for the fits. 22.5 - 32.5 MeV

Data Regular (Gamma = 27.5 MeV, Mass = 3.774 GeV)

Test DHadFit.py -t d -p Gamma_27.5_M_3.774 ... OK. qsub DHadFit_d_Gamma_27.5_M_3.774.sh job 428814 ... Discard (27.5->0.0275)

Data Wide Delta E

Test DHadFit.py -t d -p widede_Gamma_27.5_M_3.774 ... OK. qsub DHadFit_d_widede_Gamma_27.5_M_3.774.sh job 428815 ... Discard (27.5->0.0275)

Data Gamma = 32.5 MeV

Test DHadFit.py -t d -p Gamma_32.5_M_3.774 ... OK. qsub DHadFit_d_Gamma_32.5_M_3.774.sh job 428816 ...Discard (27.5->0.0275)

Data Gamma = 22.5 MeV

Test DHadFit.py -t d -p Gamma_22.5_M_3.774 ... OK. qsub DHadFit_d_Gamma_22.5_M_3.774.sh job 428817 ... Discard (27.5->0.0275)

Data Gamma = 27 MeV for D0 modes

Test DHadFit.py -t d -p Gamma_27_M_3.774 -m 0,1,3 ... OK. qsub DHadFit_d_Gamma_27_M_3.774_013.sh job 428818 ... Done.

Strange behavior. See plots D0D0B : Gamma 27 MeV. Ask ... Gamma value is bad. Discard (27.5->0.0275)

Data Gamma = 20 MeV for D+ modes

Test DHadFit.py -t d -p Gamma_20_M_3.774 -m 200,201,202,203,204,205 qsub DHadFit_d_Gamma_20_M_3.774_200_205.sh job 428819 ... Discard (27.5->0.0275)

Data R = 12.7

Test DHadFit.py -t d -p R_12.7_Gamma_27.5_M_3.774 ... OK qsub DHadFit_d_R_12.7_Gamma_27.5_M_3.774.sh job 428821 ... Discard (27.5->0.0275)

Signal MC Regular

Test DHadFit.py -p Gamma_27.5_M_3.774 ... OK. qsub DHadFit_Gamma_27.5_M_3.774.sh job 428822 ... Discard (27.5->0.0275)

Signal MC fix p = 0.5

Test DHadFit.py -p P_0.5_Gamma_27.5_M_3.774 ... OK. qsub DHadFit_P_0.5_Gamma_27.5_M_3.774.sh job 428823 ... Discard (27.5->0.0275)

Signal MC no FSR

Test DHadFit.py -p noFSR_Gamma_27.5_M_3.774 ... OK. qsub DHadFit_noFSR_Gamma_27.5_M_3.774.sh job 428824 ... Discard (27.5->0.0275)

Signal Wide Delta E

Test DHadFit.py -p widede_Gamma_27.5_M_3.774 ... OK.

qsub DHadFit_widede_Gamma_27.5_M_3.774.sh job 428825 ...Discard (27.5->0.0275)

Generic MC Regular

Test DHadFit.py -t g -p Gamma_27.5_M_3.774 ... OK.

qsub DHadFit_g_Gamma_27.5_M_3.774.sh job 428826 ...Discard (27.5->0.0275)

07/04/2006

Check the previous fits

Gamma Value is bad, change all the Gamma from 27.5 to 0.0275

toolbox.py -a 2

Data Regular (Gamma = 0.0275 GeV, Mass = 3.774 GeV)

Test DHadFit.py -t d -p Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_d_Gamma_0.0275_M_3.774.sh job 429409 ... OK.

• Regular Fit

Data Wide Delta E

Test DHadFit.py -t d -p widede_Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_d_widede_Gamma_0.0275_M_3.774.sh job 429410 ... OK.

• Wide delta E

Data Gamma = 0.0325 GeV

Test DHadFit.py -t d -p Gamma_0.0325_M_3.774 ... OK. qsub DHadFit_d_Gamma_0.0325_M_3.774.sh job 429411 ... OK.

• Gamma = 32.5 MeV

Data Gamma = 0.0225 GeV

Test DHadFit.py -t d -p Gamma_0.0225_M_3.774 ... OK. qsub DHadFit_d_Gamma_0.0225_M_3.774.sh job 429412 ... OK.

• Gamma = 22.5 MeV

Data Gamma = 0.027 MeV for D0 modes

Test DHadFit.py -t d -p Gamma_0.027_M_3.774 -m 0,1,3 ... OK. qsub DHadFit_d_Gamma_0.027_M_3.774_013.sh job 429413 ... OK.

• D0D0B : Gamma 27 MeV

Data Gamma = 0.020 GeV for D+ modes

Test DHadFit.py -t d -p Gamma_0.02_M_3.774 -m 200,201,202,203,204,205 qsub DHadFit_d_Gamma_0.02_M_3.774_200_205.sh job 429414 ... OK.

• D+ D- : Gamma 20 MeV

Data R = 12.7

Test DHadFit.py -t d -p R_12.7_Gamma_0.0275_M_3.774 ... OK qsub DHadFit_d_R_12.7_Gamma_0.0275_M_3.774.sh job 429415 ... OK.

• R = 12.7

Signal MC Regular

Test DHadFit.py -p Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_Gamma_0.0275_M_3.774.sh job 429416 ... OK.

• Regular Fit

Signal MC fix p = 0.5

Test DHadFit.py -p P_0.5_Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_P_0.5_Gamma_0.0275_M_3.774.sh job 429417 ... OK.

• Fit with p=0.5

Signal MC no FSR

Test DHadFit.py -p noFSR_Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_noFSR_Gamma_0.0275_M_3.774.sh job 429419 ... OK.

• No FSR Fit

Signal Wide Delta E

Test DHadFit.py -p widede_Gamma_0.0275_M_3.774 ... OK.

qsub DHadFit_widede_Gamma_0.0275_M_3.774.sh job 429420 ... OK.

• Wide delta E

Generic MC Regular

Test DHadFit.py -t g -p Gamma_0.0275_M_3.774 ... OK.

qsub DHadFit_g_Gamma_0.0275_M_3.774.sh job 429421 ... OK.

• Regular Fit

Waiting ...

07/05/2006

Table

DHadTable.py -a 6

Produce the Efficiency column

07/06/2006

Table calc

Single tag efficiency update

DHadTable.py -a 2

• Single tag efficiencies in signal MC

NoFSR Single tag efficiency update

DHadTable.py -a 5

[

Base Input

Auxiliary

Final ]

To do after the meeting

• Update all the sigle tag tables
• Create table of new widede without FSR
• Sort out the tables accordingly

07/07/2006

Update the sigle tag tables

• Base : Mode Yield Generated Background
• Aux1 : Mode Yield - Bkg Generated Efficiency
• Aux2 : Mode eff_sig eff_gen e_sig/e_gen Number 1 e_sig/e_gen -1
• Final : Mode Yield Generated Background Efficiency e_sig/e_gen -1

widede without FSR

Select the widede_noFSR evts

DHad.py -p widede_noFSR ... Done.

Do the fits for signal MC

Test DHadFit.py -p widede_noFSR_Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_widede_noFSR_Gamma_0.0275_M_3.774.sh job 431104 ... Done.

Re-select the widede evts for signal MC to have the log file

DHad.py -p widede ... Done.

Fit with p=0.5 for data

Test DHadFit.py -t d -p P_0.5_Gamma_0.0275_M_3.774 ... OK. qsub DHadFit_d_P_0.5_Gamma_0.0275_M_3.774.sh job 431102 ... Done.

07/08/2006

widede without FSR signal MC

• Wide delta E without FSR

Fit with p=0.5 for data

• Fit with p=0.5

Table of data with different fixed values of Gamma

DHadTable.py -a 3

07/09/2006

Table of data with different fixed values of Gamma

DHadTable.py -a 3

• Single tag yields for data with different fixed values of Gamma

Delta E ratios table

Ratios of single tag data yields and efficiencies with the standard Delta E cuts(Yield1 and eff1) to those with |Delta E| < 100 MeV (Yield2 and eff2)

• Ratios of single tag wide Delta E (Aux0, Aux1, Aux2)

Fix p table

Ratios of single tag data yields and efficiencies with a power parameter floating (Yield1 and eff1) to those with a power parameter fixed to 0.5 (Yield2 and eff2)

• Ratios of single tag P = 0.5 (Aux0, Aux1)

07/10/2006

Update the tables in the CBX

• Check out the latest version.
• Update the tables. Commit... Done.

Update the plots

toolbox.py -a 3

cvs com

Done.

Write up the report

DHadTable.py -a rpt1

DHadTable.py -a rpt2

DHadTable.py -a rpt3

DHadTable.py -a rpt4

Problem noticed on mode 201 for Gamma float data

Check the fitting log

DHadFit_d_Gamma_float_M_3.774.log.o418731

Talk with Anders, looked at here:

   9  xi          -1.03845e+01   1.10455e+00   3.15939e-04   2.19673e+00

And in the code (mbc_singletag_3s.C):

   xi= new RooRealVar("xi","xi",-20.0,-100.0,-1.0);

The lower limit for the xi is -1.0. So, we need to change this to get over it.

    xi= new RooRealVar("xi","xi",-20.0,-100.0,-0.1);

Test DHadFit.py -t d -p Gamma_float_M_3.774 -m 201 OK. qsub DHadFit_d_Gamma_float_M_3.774_201.sh job 432862 ... Done.

Fix table : Ratios of single tag yields - Delta E , noFSR

DHadTable.py -a 1

• Ratios of single tag wide Delta E (Aux0, Aux1, Aux2)

Update the CBX. Done.

07/11/2006

Check the fit for 201

log 432862 Wrong Mass value 3.77. Error happened in the sh file.

Change to DHadFit.py -t d -p Gamma_float_M_3.774 -m 201, and resubmit:

qsub DHadFit_d_Gamma_float_M_3.774_201.sh job 433784 ... Complete.

Looks OK now.

Write up the report

Combine the yields table

07/12/2006

Combine the yields table

DHadTable.py -a rpt1

Make plot for Delta E >0 and Delta E <0

Test DHad.py -p posde -m 0 OK.

DHad.py -p posde ... done.

DHad.py -t d -p posde qsub DHad_d_posde.sh job 435146 ... done.

DHad.py -p negde ... done.

DHad.py -t d -p negde qsub DHad_d_negde.sh job 435149 ... done.

Fit Kpi for delta E > 0 and delta E < 0

DHadFit.py -t d -p posde_Gamma_0.0275_M_3.774 -m 0 ...

DHadFit.py -t d -p negde_Gamma_0.0275_M_3.774 -m 0 ...

DHadFit.py -t d -p posde_Gamma_0.0275_M_3.774 -m 1 ...

DHadFit.py -t d -p negde_Gamma_0.0275_M_3.774 -m 1 ...

• positive delta E
• negtive delta E

Draft the report

• PTA July 13 (PDF)

07/13/2006

DHad plot on wide Delta E for CBX

Copy the pictures to the destination.

toolbox.py -a cp_widede

check in ... Done.

Do for MC as well ... Done.

07/17/2006

Select the Double Tag events

DHad.py —tag d ... Done.

Test DHad.py -t d —tag d OK. qsub DHad_d_d.sh job 439510 ... Done.

Test DHad.py -t g —tag d OK. qsub DHad_g_d.sh job 439511 ... Done.

Double Tag fit

Refers to the original CBX, get the fit parameters table update.

Use the option 'mfs' to fit : m - D mass f - fa, fb, sa, sb s - sigma_E

Change the mass of psi(3770)

07/18/2006

Check The Selection

Generic one still running ... Done.

Start fit for data

Test DHadFit.py -t d —tag d -p diag_floatmass_fitparams_floatsigmaE OK.

sub DHadFit_d_d_diag_floatmass_fitparams_floatsigmaE.sh job 439676 ... Done.

• Floatmass fit params float sigmaE

Notice evt difference between Peters

Sent email to Peter.

Try if pte.ecm < 3.7:

qsub DHad_d_d.sh job 439681 ... Done.

qsub DHad_g_d.sh job 439682 ...

07/19/2006

Check the Evt diff with Peter's

wc /home/xs32/work/CLEO/analysis/DHad/dat/evt/data_Double_D0_to_Kpi__D0B_to_Kpipi0.evt

    1441  4323 38907 (old values: 1445    4335   39015)

wc /nfs/cor/user/ponyisi/hadD/summerconf/txt/data_Double_D0_to_Kpi__D0B_to_Kpipi0.evt

    1441    4323   34584

Line Numbers are same, which means we got the same number of evts.

Start fit for data parameters

qsub DHadFit_d_d_diag_floatmass_fitparams_floatsigmaE.sh job 439756 ... Killed.

Construct the webpage for the previous results at the same time. Test for autumatic refresh webpage in the .sh file.

• Floatmass fit params float sigmaE

Need to add more parameter output in the .C file.

Kill the qsubjob now: qdel 439756

Then re-submit the fit :

Test output format: OK.

qsub DHadFit_d_d_diag_floatmass_fitparams_floatsigmaE.sh job job 439762 ... Done.

Investigate the fit parameters

mres = 3.77 gammad=0.025

Change? Send email to Anders. ...

Construct line shape parameters table

DHadTable.py -a line_shape_parameters

07/20/2006

Construct line shape parameters table

DHadTable.py -a line_shape_parameters

• Line shape parameters

Diag Double Tag fit for signal MC

Test DHadFit.py —tag d -p diag OK.

qsub DHadFit_d_diag.sh job 441196 ... Done.

• Diagnal Double Tag Fit For Singal MC

Diag Double Tag fit for Data

Test DHadFit.py -t d —tag d -p diag OK.

qsub DHadFit_d_d_diag.sh job 441197 ... OK.

• Diagnal Double Tag Fit For Data

Diag Double Tag fit for generic MC

Test DHadFit.py -t g —tag d -p diag OK.

qsub DHadFit_g_d_diag.sh job 441199 ... OK.

• Diagnal Double Tag Fit For Generic MC

Non-Diag Double Tag Fit

Signal MC

Test DHadFit.py —tag d -p nondiag OK. qsub DHadFit_d_nondiag.sh job 441221 ... OK.

Generic MC

Test DHadFit.py -t g —tag d -p nondiag OK. qsub DHadFit_g_d_nondiag.sh job 441222 ...OK.

Data

Test DHadFit.py -t d —tag d -p nondiag OK. qsub DHadFit_d_d_nondiag.sh job 441223 ...OK.

All the Double Tag Fit

Signal MC

Test DHadFit.py —tag d OK. qsub DHadFit_d_tag.sh job 441224 ...OK.

Generic MC

Test DHadFit.py -t g —tag d OK. qsub DHadFit_g_d.sh job 441225 ...OK.

Data

Test DHadFit.py -t d —tag d OK. qsub DHadFit_d_d.sh job 441227 ...OK.

07/21/2006

Check the fitting result

Diag Double Tag fit

• Diagnal Double Tag Fit For Data
• Diagnal Double Tag Fit For Generic MC

Non-Diag Double Tag Fit

• Non-Diagnal Double Tag Fit For Signal MC
• Non-Diagnal Double Tag Fit For Data
• Non-Diagnal Double Tag Fit For Generic MC

All the Double Tag Fit

• All Double Tag Fit For Signal MC
• All Double Tag Fit For Data
• All Double Tag Fit For Generic MC

Create Table signal_double_eff

DHadTable.py -a signal_double_eff

07/24/2006

Create Table signal_double_eff

DHadTable.py -a signal_double_eff

Can get this now:

 Mode | Generated | Signal yield | Efficiency(%)

Need the no FSR eff

Select no FSR evts for double tag signal MC

DHad.py —tag d -p noFSR

Fit the events

Test DHadFit.py —tag d -p noFSR OK.

qsub DHadFit_d_noFSR.sh job 445972 ... Done.

Create Table generic_double_eff

DHadTable.py -a generic_double_eff

07/25/2006

Check the fitting result

• no FSR all double tag fit For Signal MC

Finish the signal_double_eff table add noFSR

DHadTable.py -a signal_double_eff

• Double tag efficiencies without FSR in signal MC (Auxiliary)

Update the CBX ... OK.

Create Table generic_double_eff

DHadTable.py -a generic_double_eff

07/26/2006

Create Table generic_double_eff

DHadTable.py -a generic_double_eff

• Double tag efficiencies for generic MC (Auxiliary 1 , Auxiliary 2)

Update the CBX ... OK.

Chang Gamma for data Double fit

Gamma = 27.5 , 22.5 , 32.5 MeV

Set the default value in DHadFit to be 0.0275 GeV for double tag

Test DHadFit.py -t d —tag d ... OK. qsub DHadFit_d_d.sh job 450572 ... OK.

Test DHadFit.py -t d —tag d -p Gamma_0.0225 ... OK. qsub DHadFit_d_d_Gamma_0.0225.sh job 450573 ... OK.

Test DHadFit.py -t d —tag d -p Gamma_0.0325 ... OK. qsub DHadFit_d_d_Gamma_0.0325.sh job 450574 ... OK.

07/27/2006

Check the fitting result

• All Double Tag Gamma = 27.5 MeV
• Gamma = 22.5 MeV
• Gamma = 32.5 MeV

p = 0.5 fit

Test DHadFit.py -t d —tag d -p P_0.5 OK

qsub DHadFit_d_d_P_0.5.sh job 451980 ... OK.

• Fix p=0.5

Fit mass = 3.774 GeV For all

Data

Test DHadFit.py -t d —tag d -p M_3.774 OK. qsub DHadFit_d_d_M_3.774.sh job 451981 ... OK.

• Gamma = 27.5 MeV

Test DHadFit.py -t d —tag d -p Gamma_0.0225_M_3.774 ... OK. qsub DHadFit_d_d_Gamma_0.0225_M_3.774.sh job 451982 ... OK.

• Gamma = 22.5 MeV

Test DHadFit.py -t d —tag d -p Gamma_0.0325_M_3.774 ... OK. qsub DHadFit_d_d_Gamma_0.0325_M_3.774.sh job 451983 ... OK.

• Gamma = 32.5 MeV

Test DHadFit.py -t d —tag d -p P_0.5_M_3.774 ... OK. qsub DHadFit_d_d_P_0.5_M_3.774.sh job 451985 ... OK.

• Fix p=0.5

Signal MC

Test DHadFit.py —tag d -p M_3.774 OK. qsub DHadFit_d_tag_M_3.774.sh job 451986 ... OK.

• All Double Tag Fit For Signal MC

Test DHadFit.py —tag d -p noFSR_M_3.774 OK. qsub DHadFit_d_noFSR_M_3.774.sh job 451987 ... OK.

• no FSR all double tag fit For Signal MC

Generic MC

Test DHadFit.py -t g —tag d -p M_3.774 OK. qsub DHadFit_g_d_M_3.774.sh job 451988 ... OK.

• All Double Tag Fit For Generic MC

Create Table data_double_yields

DHadTable.py -a data_double_yields

• Double tag yields for data

Compare with the CBX ... consistent.

07/28/2006

Update the double tag MC table for Mass = 3.774 GeV

DHadTable.py -a signal_double_eff

• Double tag efficiencies with and without FSR in signal MC (Auxiliary)

Update CBX table OK.

DHadTable.py -a generic_double_eff

• Double tag efficiencies for generic MC (Auxiliary 1 , Auxiliary 2)

Update CBX table OK.

cvs commit ... Done.

07/31/2006

Background shape for single tag fit

Message from Anders:

 This has to do with the background shapes for the single
tags.

 To select the signal events we cut on delta E around 0.
But we can use the Delta E sidebands to model the backgrounds.

 So, can you select the events, separately, in the high and low
Delta E sidebands and fit the mBC distribution for these
events to the Argus background shape.
Then using the parameters you find for the shape (not yield),
fix the shape parameters for the background and fit
the delta E signal region. For most modes you can select
the region in Delta E from 50 to 100MeV as the high sideband
and -100 to -50 MeV as the low sideband. In modes with pi0's
you might want to make the sideband even smaller to make
sure you don't have signal in the sideband region.

Select the events in the high and low Delta E sidebands separately

Test DHad.py -p desidebandl OK. qsub DHad_desidebandl.sh job 464410 ... Done.

Test DHad.py -p desidebandh OK. qsub DHad_desidebandh.sh job 464415 ... Done.

Test DHad.py -t d -p desidebandl OK. qsub DHad_d_desidebandl.sh job 464416 ... Done.

Test DHad.py -t d -p desidebandh OK. qsub DHad_d_desidebandh.sh job 464427 ... Done.

Test DHad.py -t g -p desidebandl OK. qsub DHad_g_desidebandl.sh job 464430 ... Done.

Test DHad.py -t g -p desidebandh OK. qsub DHad_g_desidebandh.sh job 464434 ... Done.

Floatmass fit params float sigmaE for data with Mass = 3.774 GeV

Test DHadFit.py -t d —tag d -p diag_floatmass_fitparams_floatsigmaE_M_3.774 OK.

qsub DHadFit_d_d_diag_floatmass_fitparams_floatsigmaE_M_3.774.sh job 464473 ->

Warning message:

Warning: wrong member access operator '.' FILE:/home/xs32/work/CLEO/analysis/DHad/script/Roofit/lineshapefit2d.C LINE:319

But the output txt file is still OK.

Make individual plots for the fit

mbcFrame->SetMinimum(0.5);

>>> 08/01/2006