D Had Log - October 2006

/home/xs32/work/CLEO/analysis/DHad/doc/fig/david/data_Double_D0_to_Kpi__D0B_to_Kpipipi_log.eps
/home/xs32/work/CLEO/analysis/DHad/doc/fig/david/data_Double_D0_to_Kpi__D0B_to_Kpipipi_linear.eps

/home/xs32/work/CLEO/analysis/DHad/doc/fig/david/data_Double_D0_to_Kpipipi__D0B_to_Kpi_log.eps
/home/xs32/work/CLEO/analysis/DHad/doc/fig/david/data_Double_D0_to_Kpipipi__D0B_to_Kpi_linear.eps

Use .root file to plot

TFile *_file0 = TFile::Open("data_Double_D0_to_Kpipipi__D0B_to_Kpi.root")
gSystem.Load("libRooFitCore.so")
gSystem.Load("libRooFitModels.so")
.ls

TFile**         data_Double_D0_to_Kpipipi__D0B_to_Kpi.root
 TFile*         data_Double_D0_to_Kpipipi__D0B_to_Kpi.root
  OBJ: RooPlot  frame(09a253c0) D^{0}#rightarrowK^{-} #pi^{+} #pi^{+} #pi^{-} : 0 at: 0x94c4d50
  OBJ: RooPlot  frame(09aba238) #bar{D}^{0}#rightarrowK^{+} #pi^{-} : 0 at: 0x95400a0
  KEY: TCanvas  c;1     DLineShape

c.Draw()

10/12/2006

Input the paper

Improve plot the generic without the signal in it by using the MC truth info.

Test the single MC

DHad.py -t s -p dir_nosig_desidebandl -m 0

OK.

DHad.py -t g -p dir_nosig_desidebandl -m 0 —qsub job 9442 done.

DHad.py -t g -p dir_nosig_desidebandh -m 0 —qsub job 9443 done.

Get the Plots macro working

DHadPlot.py -a macro -o /home/xs32/work/CLEO/analysis/DHad/doc/fig/david/test.C

10/13/2006

D0-> K- Pi+ , D0 -> K pi pi pi invariant mass plot DT,

Signal MC mode 0

DHadPlot.py -a dmraw -o s,0

Save the histgram to

 /home/xs32/work/CLEO/analysis/DHad/doc/fig/richard/dmraw_signal_Single_D0_to_Kpi.root

Signal MC mode 3 DHadPlot.py -a dmraw -o s,3

Save the histgram to

/home/xs32/work/CLEO/analysis/DHad/doc/fig/richard/dmraw_signal_Single_D0_to_Kpipipi.root

gDirectory->Get("D raw Mass").Draw()

Data mode 0 DHadPlot.py -a dmraw -o d,0

Save the histgram to

 /home/xs32/work/CLEO/analysis/DHad/doc/fig/richard/dmraw_data_Single_D0_to_Kpi.root

Data mode 3 DHadPlot.py -a dmraw -o d,3

Save the histgram to

/home/xs32/work/CLEO/analysis/DHad/doc/fig/richard/dmraw_data_Single_D0_to_Kpipipi.root

10/14/2006

Look up the D0 Mass

location of the file :

$C3_DATA/evt.pdl

*      type        name                 id  mass/GeV width/GeV max_Dm/GeV 3*charge 2*spin lifetime*c/mm Lund-KC
...
add  p Meson       D0                  421  1.8645          0       0        0       0  0.1234 105

add  p Meson       D+                  411  1.8693          0       0        3       0  0.315 104
add  p Meson       D-                 -411  1.8693          0       0       -3      0   0.315 0

10/16/2006

Plot the generic without the signal in it by using the MC truth info.

Improve the cuts in DHadTools.checkVetoes_nosig

See the memo here for the coding detail.

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)

Remove the number in the 1000000*(#gamma) and then impose the cuts

Test:

DHad.py -t s -p dir_nosig/test_desidebandl -m 0

DHad.py -t g -p dir_nosig_desidebandl -m 0 —qsub job 21180

DHad.py -t g -p dir_nosig_desidebandh -m 0 —qsub job 21181

Do not use f1b and s1b for modes without pi0 in double tag fit

Test: test-DHadFit.py -t s -p dir_Nopi0Fixfbsb -m 0

10/17/2006

Check the selection job

job 21180 done. Compare the dat/fit/log/21180 with the previous one - 9442

Name	Value (21180)	(9442)
Mode	0	0
Sign	1	1
Total	17143020	17143020
Select	2776	4396
Eff	0.000161931794981	0.00025643089724

There are less events selected, which means DCSD and radiative decays have been vetoed.

Do the fits for the nosig

To save the previous plots in the nosig dir, now use nosig1 in stead.

test : DHadFit.py -t g -p dir_nosig1_desidebandh -m 0 OK.

DHadFit.py -t g -p dir_nosig1_desidebandl -m 0 —qsub job 21992 done

no signal for mode 0 (Excluding the DCSD and radiative signal)

DHadFit.py -t g -p dir_nosig1_desidebandh -m 0 —qsub job 21995 done

no signal for mode 0 (Excluding the DCSD and radiative signal)

Do not use f1b and s1b for modes without pi0

Test: test-DHadFit.py -t s -p dir_Nopi0Fixfbsb -m 0 OK.

DHadFit.py -t s -p dir_Nopi0Fixfbsb —qsub job 22019 done.

DHadFit.py -t d -p dir_Nopi0Fixfbsb —qsub job 22020 done.

DHadFit.py -t g -p dir_Nopi0Fixfbsb —qsub job 22021 done.

10/18/2006

Do not use f1b and s1b for modes without pi0 in double tag fit

Fit finished.

Mode 203 'kspipi0':


 COVARIANCE MATRIX CALCULATED SUCCESSFULLY
 FCN=-1.46105e+06 FROM HESSE     STATUS=OK             73 CALLS         453 TOTAL
                     EDM=0.000284605    STRATEGY= 1      ERROR MATRIX ACCURATE
  EXT PARAMETER                                INTERNAL      INTERNAL
  NO.   NAME      VALUE            ERROR       STEP SIZE       VALUE
   1  N1           2.29594e+04   2.47945e+02   1.38414e-04  -1.35610e+00
   2  N2           2.27281e+04   2.47999e+02   2.77411e-05  -1.35719e+00
   3  Nbkgd1       2.09719e+04   2.43912e+02   1.40743e-04  -1.36563e+00
   4  Nbkgd2       2.12125e+04   2.44914e+02   1.40478e-04  -1.36446e+00
   5  md           1.86928e+00   1.78511e-05   7.84651e-02   4.72692e-02
   6  p            1.50000e+00   8.65678e-01   8.25741e-03   1.57078e+00
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   7  sigmap1      7.98473e-03   1.23136e-04   1.88897e-01  -7.54623e-01
   8  xi          -2.44766e+01   7.07967e-01   3.78677e-04   5.37489e-01
                               ERR DEF= 0.5

Signal MC:
not use f1b and s1b for modes without pi0

Mode 203 'kspipi0':

 COVARIANCE MATRIX CALCULATED SUCCESSFULLY
 FCN=-516792 FROM HESSE     STATUS=OK             69 CALLS         641 TOTAL
                     EDM=0.00321777    STRATEGY= 1      ERROR MATRIX ACCURATE
  EXT PARAMETER                                INTERNAL      INTERNAL
  NO.   NAME      VALUE            ERROR       STEP SIZE       VALUE
   1  N1           1.39256e+04   1.24613e+02   1.43722e-05  -1.40372e+00
   2  N2           1.39313e+04   1.24312e+02   1.43571e-05  -1.40368e+00
   3  Nbkgd1       1.54082e+03   5.60689e+01   8.72199e-05  -1.51528e+00
   4  Nbkgd2       1.42277e+03   5.42703e+01   1.75736e-05  -1.51745e+00
   5  md           1.86942e+00   1.15569e-05   1.94301e-04   7.07600e-02
   6  p            3.45055e-01   6.17949e-02   9.85215e-04  -7.07481e-01
   7  sigmap1      2.00000e-03   1.41871e-06   1.98228e-01  -1.57086e+00
                                 WARNING -   - ABOVE PARAMETER IS AT LIMIT.
   8  xi          -1.17180e+01   3.50481e+00   8.45761e-04   8.74787e-01
                               ERR DEF= 0.5

Generic MC:
not use f1b and s1b for modes without pi0

   7  sigmap1      2.00000e-03   1.58035e-06  -0.00000e+00** at limit **

Reasnos? Possible solutions?

Change the fa and sigma +/- one sigma

Test: DHadFit.py -t s -p dir_faplus -m 0 OK. DHadFit.py -t s -p dir_faminus -m 0 OK.

DHadFit.py -t s -p dir_faplus —qsub job 23111

Use fa plus one sigma

DHadFit.py -t s -p dir_faminus —qsub job 23112

Use fa minus one sigma

DHadFit.py -t d -p dir_faplus —qsub job 23113

Use fa plus one sigma

DHadFit.py -t d -p dir_faminus —qsub job 23114

Use fa minus one sigma

DHadFit.py -t g -p dir_faplus —qsub job 23115

Use fa plus one sigma

DHadFit.py -t g -p dir_faminus —qsub job 23116

Use fa minus one sigma

all done.

10/19/2006

Nopi0Fixfbsb fit to compare table

Yields compare

DHadTable.py -a compare -o yields_cleo200609_Nopi0Fixfbsb

Compare yields between cleo200609 and Nopi0Fixfbsb

Change the fa +/- one sigma fit compare

Check the job - all done.
Yields compare

DHadTable.py -a compare -o yields_cleo200609_faplus

Compare yields between cleo200609 and faplus

DHadTable.py -a compare -o yields_cleo200609_faminus

Compare yields between cleo200609 and faminus

10/20/2006

Clarify the tables page

Update: Line shape parameters
Create page Tables

10/21/2006

Fit the parameters with the fixed sigmaE (2.1MeV)

DHadFit.py -t s —tag d -p dir_Nopi0Fixfbsb_diag_floatmass_fitparams —qsub job 25005

10/23/2006

Check the fitting results

job done.

Use f1b and s1b for modes with pi0

Create table

Use f1b and s1b for modes with pi0

Do the ST fit based on the table

test DHadFit.py -t d -p dir_Nopi0Fixfbsb OK.

DHadFit.py -t d -p dir_Nopi0Fixfbsb —qsub job 31306 done

DHadFit.py -t s -p dir_Nopi0Fixfbsb —qsub job 31307 done

DHadFit.py -t g -p dir_Nopi0Fixfbsb —qsub job 31308 done

10/24/2006

Check fitting results

Used the wrong table for parameters. Fix it.

Resubmit:

DHadFit.py -t d -p dir_Nopi0Fixfbsb —qsub job 31383

not use f1b and s1b for modes without pi0

DHadFit.py -t s -p dir_Nopi0Fixfbsb —qsub job 31384

not use f1b and s1b for modes without pi0

DHadFit.py -t g -p dir_Nopi0Fixfbsb —qsub job 31385

not use f1b and s1b for modes without pi0

Fit fa plus minus

test-DHadFit.py -t d -p dir_faplus OK.

DHadFit.py -t d -p dir_faplus —qsub job 31388 DHadFit.py -t d -p dir_faminus —qsub job 31389

DHadFit.py -t s -p dir_faplus —qsub job 31390 DHadFit.py -t s -p dir_faminus —qsub job 31391

DHadFit.py -t g -p dir_faplus —qsub job 31392 DHadFit.py -t g -p dir_faminus —qsub job 31403

done.

Fit sa plus minus

test-DHadFit.py -t d -p dir_saplus OK.

DHadFit.py -t d -p dir_saplus —qsub job 31404 DHadFit.py -t d -p dir_saminus —qsub job 31405

DHadFit.py -t s -p dir_saplus —qsub job 31406 DHadFit.py -t s -p dir_saminus —qsub job 31407

DHadFit.py -t g -p dir_saplus —qsub job 31408 DHadFit.py -t g -p dir_saminus —qsub job 31409

done.

Scatter Plots

deltaE vs. mbc without signal for generic MC

DHadPlot.py -a deltaE_vs_mbc -o nosig

Need to tune the selection

10/25/2006

Check the fitting results

all done.

Change the fa +/- one sigma fit compare

Check the job - all done.
Yields compare DHadTable.py -a compare -o yields_cleo200609_Nopi0Fixfbsb
Compare yields between cleo200609 and Nopi0Fixfbsb

DHadTable.py -a compare -o yields_cleo200609_faplus

Compare yields between cleo200609 and faplus

DHadTable.py -a compare -o yields_cleo200609_faminus

Compare yields between cleo200609 and faminus

DHadTable.py -a compare -o yields_cleo200609_saplus

Compare yields between cleo200609 and saplus

DHadTable.py -a compare -o yields_cleo200609_saminus

Compare yields between cleo200609 and saminus

10/26/2006

Create scatter plot

deltaE vs. mbc without signal for generic MC

DHadPlot.py -a deltaE_vs_mbc -o desidebandl

DHadPlot.py -a deltaE_vs_mbc -o desidebandl_nosig

10/30/2006

Scatter plots

DHadPlot.py -a deltaE_vs_mbc -o noDEcut

10/31/2006

Scatter plots

No Delta E cuts

DHadPlot.py -a deltaE_vs_mbc -o noDEcut

-1: 262669 +1: 257986

deltaE vs. Mbc without delta E cuts
With Delta E cuts

DHadPlot.py -a deltaE_vs_mbc

deltaE vs. Mbc with delta E cuts

Without Delta E cut, with no signal

DHadPlot.py -a deltaE_vs_mbc -o noDEcut_nosig

without delta E cuts, no signal

[ Front | MyDHad-archives | Index ]

Updated January 11, 2007

D Had Log - October 2006

Home > D-Hadronic > Archives

Archives

Links

10/04/2006

Check the previous selection

Do the fits for the nosig

Fit the two Gaussian fuction, see the yields

10/05/2006

PRD Draft

10/11/2006

Create plot for David

Use .root file to plot

10/12/2006

Input the paper

Improve plot the generic without the signal in it by using the MC truth info.

Get the Plots macro working

10/13/2006

D0-> K- Pi+ , D0 -> K pi pi pi invariant mass plot DT,

10/14/2006

Look up the D0 Mass

10/16/2006

Plot the generic without the signal in it by using the MC truth info.

Do not use f1b and s1b for modes without pi0 in double tag fit

10/17/2006

Check the selection job

Do the fits for the nosig

Do not use f1b and s1b for modes without pi0

10/18/2006

Do not use f1b and s1b for modes without pi0 in double tag fit

Change the fa and sigma +/- one sigma

10/19/2006

Nopi0Fixfbsb fit to compare table

Change the fa +/- one sigma fit compare

10/20/2006

Clarify the tables page

10/21/2006

Fit the parameters with the fixed sigmaE (2.1MeV)

10/23/2006

Check the fitting results

Create table

Do the ST fit based on the table

10/24/2006

Check fitting results

Fit fa plus minus

Fit sa plus minus

Scatter Plots

10/25/2006

Check the fitting results

Change the fa +/- one sigma fit compare

10/26/2006

Create scatter plot

10/30/2006

Scatter plots

10/31/2006

Scatter plots